Tracing the stellar halo with BHB stars
Guillaume THOMAS
@Thomas_gft
Stellar Halos Across the Cosmos 3rd July 2018
Tracing the stellar halo with BHB stars Guillaume THOMAS - - PowerPoint PPT Presentation
Tracing the stellar halo with BHB stars Guillaume THOMAS @Thomas_gft Stellar Halos Across the Cosmos 3 rd July 2018 Mon. Not. R. Astron. Soc. 000 , 1 ?? (2018) Printed 30 May 2018 (MN L A T EX style file v2.2) Blue horizontal branch
@Thomas_gft
Stellar Halos Across the Cosmos 3rd July 2018
Printed 30 May 2018 (MN L
AT EX style file v2.2)
Blue horizontal branch stars in the Canada-France Imaging Survey I. The stellar halo of the Milky Way traced to large radius
Guillaume F. Thomas1, Alan McConnachie1, Rodrigo A. Ibata2, Patrick Cˆ
e1, Nicolas Martin2,3, Nicholas Fantin1, S´ ebastien Fabbro1, Benoit Famaey2, Raymond Carlberg4, Stephen Gwyn1, Vincent Henault-Brunet1, Khyati Malhan2, Julio Navarro5, Annie C. Robin6, Douglas Scott7, Else Starkenburg8
1NRC Herzberg Astronomy and Astrophysics, 5071 West Saanich Road, Victoria, BC, V9E 2E7, Canada 2Observatoire astronomique de Strasbourg, Universit´e de Strasbourg, CNRS, UMR 7550, 11 rue de l’Universit´ e, F-67000 Strasbourg, France
3Max-Planck-Institut f¨ur Astronomie, K¨
e, OSU THETA Franche-Compt´ e-Bourgogne, Observatoire de Besan¸ con, BP 1615, 25010 Besan¸ con Cedex, France
30 May 2018
ABSTRACT
We present the stellar density profile of the outer halo of the Galaxy traced over a range
Canada-France-Imaging-Survey (CFIS), which reaches 24.5 mag, combined with griz bands from Pan-STARRS 1, covering a total of ∼ 4200 deg2 of the northern sky. We present a new method to select BHB stars that has low contamination from blue stragglers and high completeness. We use this sample to measure and parameterize the three dimensional density profile of the outer stellar halo, using both a simple power-law with a constant flattening, and a flattening that varies as a function of Galactocentric radius. In the case of a constant flattening, we find that the density profile is well described by a slope of γ = 3.42±0.02 and an oblateness of q = 1.06±0.2, consistent with the recent result of Fukushima et al. (2017). In the case of the radius- dependent flattening, we find that the inner halo is more oblate (q0 = 0.96 ± 0.03) than at large distance (q∞ = 1.25+0.07
−0.06), and has a power-law slope of γ = 3.60 ± 0.04.
With these two models, the profile of the stellar halo trace by BHB stars is shallower than when traced by RR Lyrae, a surprising result given the similarity of these stellar populations. Key words: stars: horizontal branch – stars: distances – stars: statistics – Galaxy: structure – Galaxy: halo
1 INTRODUCTION It is now generally accepted that large galaxies, like the Milky Way, have been formed by a succession of mergers and via the accretion of smaller galaxies, in a process called hierarchical formation. In the case of accretions, the smaller galaxy is disrupted due to the tidal effects generated by the larger (host) galaxy. This leads to the formation of stellar streams clearly visible around many massive galaxies of the Local Group (e.g. Mart´ ınez-Delgado et al. 2010; Martin et al. 2013; Grillmair & Carlin 2016; Bernard et al. 2016; Malhan et al. 2018). Although these structures stay spatially coher- ent for many Gyr (Johnston et al. 2008), they tend to be eventually destroyed by mixing effects and are in turn as- similated to form part of the “smooth” stellar halo.
c 2018 RAS
Ibata et al., 2017
Ibata et al., 2017
u-SDSS u-CFIS
The Ultraviolet Near-Infrared Optical Northern Survey
stellar halo or kicked out of the disc * metal rich ([Fe/H]>-1.0) * dominant < 20 kpc
accreted galaxies/globular clusters * metal poor ([Fe/H]<-1.0) * dominate the outer stellar halo
Bullock & Johnston 2005
stellar halo or kicked out of the disc * metal rich ([Fe/H]>-1.0) * dominant < 20 kpc
accreted galaxies/globular clusters * metal poor ([Fe/H]<-1.0) * dominate the outer stellar halo
Bullock & Johnston 2005
stellar halo or kicked out of the disc * metal rich ([Fe/H]>-1.0) * dominant < 20 kpc
accreted galaxies/globular clusters * metal poor ([Fe/H]<-1.0) * dominate the outer stellar halo
Bullock & Johnston 2005
“…a study of these subsystems allows us partially to reconstruct the Galactic past…” Olin Eggen
Pillepich et al., 2014, 2018
stellar halo and the total mass of a galaxy
principle progenitors
Pillepich et al., 2014, 2018 Ibata et al., 2014
stellar halo and the total mass of a galaxy
principle progenitors
Pillepich et al., 2014, 2018 Ibata et al., 2014
stellar halo and the total mass of a galaxy
Need accurate distance
principle progenitors
Accurate photometric distance (5% of precision)
−0.4 −0.2 0.0 0.2 0.4 g0 − r0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 u0 − g0 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 1.25 1.50 u0 − z0 −0.6 −0.5 −0.4 −0.3 −0.2 −0.1 0.0 0.1 r0 − i0
Accurate photometric distance (5% of precision)
−0.4 −0.2 0.0 0.2 0.4 g0 − r0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 u0 − g0 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 1.25 1.50 u0 − z0 −0.6 −0.5 −0.4 −0.3 −0.2 −0.1 0.0 0.1 r0 − i0 −0.4 −0.2 0.0 0.2 0.4 g0 − r0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 u0 − g0 −0.50 −0.25 0.00 0.25 0.50 0.75 1.00 1.25 1.50 u0 − z0 −0.6 −0.5 −0.4 −0.3 −0.2 −0.1 0.0 0.1 r0 − i0
with hydrogen lines sensitive to the surface gravity:
with hydrogen lines sensitive to the surface gravity:
74% pure and 57% complete (ugr)
77% pure and 51% complete (griz)
with hydrogen lines sensitive to the surface gravity:
74% pure and 57% complete (ugr)
77% pure and 51% complete (griz)
−0.3 −0.2 −0.1 0.0 0.1 0.2 0.3 P1 −0.20 −0.15 −0.10 −0.05 0.00 0.05 0.10 0.15 0.20 P2 − P3 BHB arHa BS arHa
with hydrogen lines sensitive to the surface gravity:
74% pure and 57% complete (ugr)
77% pure and 51% complete (griz)
−0.3 −0.2 −0.1 0.0 0.1 0.2 0.3 P1 −0.20 −0.15 −0.10 −0.05 0.00 0.05 0.10 0.15 0.20 P2 − P3 BHB arHa BS arHa
50 100 150 200 250 300 350 l (deg) −75 −50 −25 25 50 75 b (deg)
fit the BHB distribution :
p(Di|θ) = ρBHB(Di|θ) |J| S(li, bi, Di) R R R ρBHB(l, b, D|θ) |J| S(l, b, D) dl db dD . (19)
fit the BHB distribution :
p(Di|θ) = ρBHB(Di|θ) |J| S(li, bi, Di) R R R ρBHB(l, b, D|θ) |J| S(l, b, D) dl db dD . (19)
fit the BHB distribution :
Sarea(l, b) = ⇢ 1 if (l, b) in CFIS 0 otherwise
Souter halo(RGC) = ⇢ 1 if 15 < RGC < 240 kpc 0 otherwise
Sconta(l, b) = 8 > > > > < > > > > : 0 if dM31 < 4.0 deg 0 if dM33 < 2.0 deg 0 if dNGC5466 < 0.4 deg 0 if dDraco < 0.5 deg 1 otherwise
compact objects
8 9 10 11 12 13 14 5.A. (deg) 38 39 40 41 42 43 44M 31 M33
SSgr(l, b) = ⇢ 0 if | ˜ B| < 10.0 deg 1 otherwise
p(Di|θ) = ρBHB(Di|θ) |J| S(li, bi, Di) R R R ρBHB(l, b, D|θ) |J| S(l, b, D) dl db dD . (19)
Scomp(l, b, rhelio) = Cz (zBHB zlim(l, b) + zlim,ref) (
100 20 30 40 50 60 70 80 90 200 5GC (NSc) 1 10 100 1000 1b Rf BHBs SRwHr law SRwHr law q(RGC) BHB with 6gr strHam BHB withRut 6gr strHam
Two models with a single power law:
γ = 3.42 ± 0.02
γ = 3.60 ± 0.04 q = 1.06 ± 0.02
q0 = 0.96 ± 0.03
q∞ = 1.25+0.07
−0.06
rq = 46.4+12.0
−8.3
100 20 30 40 50 60 70 80 90 200 5GC (NSc) 1 10 100 1000 1b Rf BHBs SRwHr law SRwHr law q(RGC) BHB with 6gr strHam BHB withRut 6gr strHam
Two models with a single power law:
γ = 3.42 ± 0.02
γ = 3.60 ± 0.04 q = 1.06 ± 0.02
q0 = 0.96 ± 0.03
q∞ = 1.25+0.07
−0.06
rq = 46.4+12.0
−8.3
50 100 150 200 250 5GC (kpc) 0.95 1.00 1.05 1.10 1.15 1.20 1.25 1.30 q(5GC)
100 20 30 40 50 60 70 80 90 200 300 5GC (kSF) 10−4 10−3 10−2 10−1 100 DHnsLWy (ArbLWrDry XnLWs) WDWkLns HW Dl., 2009 (55 LyrDH) DH 3rRSrLs HW Dl., 2010 (BHB) DHDsRn HW Dl., 2014 (BHB) XXH HW Dl., 2015 (.-gLDnWs) 6lDWHr HW Dl., 2016 (GLDnWs) CRKHn HW Dl., 2017 (55 LyrDH) )XkXsKLPD HW Dl., 2017 (BHB) HHrnLWsFKHk HW Dl., 2018 (55 LyrDH) SRwHr lDw (BHB)
100 20 30 40 50 60 70 80 90 200 300 5GC (kSF) 10−4 10−3 10−2 10−1 100 DHnsLWy (ArbLWrDry XnLWs) WDWkLns HW Dl., 2009 (55 LyrDH) DH 3rRSrLs HW Dl., 2010 (BHB) DHDsRn HW Dl., 2014 (BHB) XXH HW Dl., 2015 (.-gLDnWs) 6lDWHr HW Dl., 2016 (GLDnWs) CRKHn HW Dl., 2017 (55 LyrDH) )XkXsKLPD HW Dl., 2017 (BHB) HHrnLWsFKHk HW Dl., 2018 (55 LyrDH) SRwHr lDw (BHB)
Profile traced by the BHB is closer to the profile traced by the giants than by the RR Lyraes
100 20 30 40 50 60 70 80 90 200 300 5GC (kSF) 10−4 10−3 10−2 10−1 100 DHnsLWy (ArbLWrDry XnLWs) WDWkLns HW Dl., 2009 (55 LyrDH) DH 3rRSrLs HW Dl., 2010 (BHB) DHDsRn HW Dl., 2014 (BHB) XXH HW Dl., 2015 (.-gLDnWs) 6lDWHr HW Dl., 2016 (GLDnWs) CRKHn HW Dl., 2017 (55 LyrDH) )XkXsKLPD HW Dl., 2017 (BHB) HHrnLWsFKHk HW Dl., 2018 (55 LyrDH) SRwHr lDw (BHB)
Profile traced by the BHB is closer to the profile traced by the giants than by the RR Lyraes
unprecedented precision by their photometry
traced by the BHB is similar to the halo traced by the giants …
stellar halo up to ~240 kpc
RR Lyrae
unprecedented precision by their photometry
traced by the BHB is similar to the halo traced by the giants …
stellar halo up to ~240 kpc
RR Lyrae
different populations with cosmological simulations (Auriga?) and/or the Galactic Besançon Model
=> Mass of the MW and the 3D distribution of the DM halo at large distances
17 18 19 20 21 22 23 24 25 26 0ag 0.0 0.2 0.4 0.6 0.8 1.0 Completeness gPS rPS iPS zPS
magnitude is not the u-band
griz bands via HSC-SSP
17 18 19 20 21 22 23 24 25 26 0ag 0.0 0.2 0.4 0.6 0.8 1.0 Completeness gPS rPS iPS zPS
magnitude is not the u-band
griz bands via HSC-SSP Limiting magnitude = z-band
17 18 19 20 21 22 23 24 25 26 0ag 0.0 0.2 0.4 0.6 0.8 1.0 Completeness gPS rPS iPS zPS
50 100 150 200 250 300 350 5.A. (deg) 20 40 60 Dec (deg) 21.6 21.8 22.0 22.2 22.4 0Dg limit
magnitude is not the u-band
griz bands via HSC-SSP Limiting magnitude = z-band
Function to see the spatial variation of the completeness Field of reference