Tracing the stellar halo with BHB stars Guillaume THOMAS - - PowerPoint PPT Presentation

tracing the stellar halo with bhb stars
SMART_READER_LITE
LIVE PREVIEW

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


slide-1
SLIDE 1

Tracing the stellar halo with BHB stars

Guillaume THOMAS

@Thomas_gft

Stellar Halos Across the Cosmos 3rd July 2018

slide-2
SLIDE 2
  • Mon. Not. R. Astron. Soc. 000, 1–?? (2018)

Printed 30 May 2018 (MN L

A

T 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¨

  • nigstuhl 17, 69117 Heidelberg, Germany
4Departement of Astronomy and Astrophysics, University of Toronto, Toronto, ON M5S 3H4, Canada 5 Departement of Physics and Astronomy, University of Victoria, Victoria, BC, V8P 1A1, Canada 6Institut UTINAM, CNRS UMR6213, Univ. Bourgogne Franche-Comt´

e, OSU THETA Franche-Compt´ e-Bourgogne, Observatoire de Besan¸ con, BP 1615, 25010 Besan¸ con Cedex, France

  • 7Dept. of Physics and Astronomy, University of British Columbia, Vancouver, B.C., V6T 1Z1, Canada
8Leibniz Institute for Astrophysics Potsdam (AIP), An der Sternwarte 16, D-14482 Potsdam, Germany

30 May 2018

ABSTRACT

We present the stellar density profile of the outer halo of the Galaxy traced over a range

  • f Galactocentric radii from 15 < RGC < 250 kpc by blue horizontal branch (BHB)
  • stars. These stars are identified photometrically using u−band imaging from the new

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

slide-3
SLIDE 3

CFIS & UNIONS

slide-4
SLIDE 4

CFIS

Ibata et al., 2017

  • u-band : 10,000 deg2
  • r-band : 5,000 deg2
  • u-band : 3 mag deeper than SDSS
slide-5
SLIDE 5

CFIS

Ibata et al., 2017

  • u-band : 10,000 deg2
  • r-band : 5,000 deg2
  • u-band : 3 mag deeper than SDSS

u-SDSS u-CFIS

slide-6
SLIDE 6

UNIONS

The Ultraviolet Near-Infrared Optical Northern Survey

  • MoU between CFIS and Pan-STARRS
  • u W r i z photometric bands
slide-7
SLIDE 7

What is the stellar halo? How does it form?

slide-8
SLIDE 8

How was the Stellar Halo formed?

  • Two mechanisms:
  • In-situ stars: formed initially in the

stellar halo or kicked out of the disc * metal rich ([Fe/H]>-1.0) * dominant < 20 kpc

  • Accreted stars: coming from the

accreted galaxies/globular clusters * metal poor ([Fe/H]<-1.0) * dominate the outer stellar halo

Bullock & Johnston 2005

slide-9
SLIDE 9

How was the Stellar Halo formed?

  • Two mechanisms:
  • In-situ stars: formed initially in the

stellar halo or kicked out of the disc * metal rich ([Fe/H]>-1.0) * dominant < 20 kpc

  • Accreted stars: coming from the

accreted galaxies/globular clusters * metal poor ([Fe/H]<-1.0) * dominate the outer stellar halo

Bullock & Johnston 2005

slide-10
SLIDE 10

How was the Stellar Halo formed?

  • Two mechanisms:
  • In-situ stars: formed initially in the

stellar halo or kicked out of the disc * metal rich ([Fe/H]>-1.0) * dominant < 20 kpc

  • Accreted stars: coming from the

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

slide-11
SLIDE 11

How was the Stellar Halo formed?

Pillepich et al., 2014, 2018

  • Correlation between the slope of the

stellar halo and the total mass of a galaxy

  • Correlation number of

principle progenitors

slide-12
SLIDE 12

How was the Stellar Halo formed?

Pillepich et al., 2014, 2018 Ibata et al., 2014

  • Correlation between the slope of the

stellar halo and the total mass of a galaxy

  • Presence of substructures bias the slope
  • Correlation number of

principle progenitors

  • Find substructures
slide-13
SLIDE 13

How was the Stellar Halo formed?

Pillepich et al., 2014, 2018 Ibata et al., 2014

  • Correlation between the slope of the

stellar halo and the total mass of a galaxy

  • Presence of substructures bias the slope

Need accurate distance

  • Correlation number of

principle progenitors

  • Find substructures
slide-14
SLIDE 14

The BHB stars How to use them as distance tracers?

slide-15
SLIDE 15

The BHB stars

  • The Blue Horizontal Branch stars:
  • Hot stars 7400 < Teff < 9300 K
  • Member of the Horizontal Branch

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

slide-16
SLIDE 16

The BHB stars

  • The Blue Horizontal Branch stars:
  • Hot stars 7400 < Teff < 9300 K
  • Member of the Horizontal Branch

Accurate photometric distance (5% of precision)

  • Contaminated by the Blue Stragglers (BS)

−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

slide-17
SLIDE 17

The BHB selection with CFIS

  • Disentangle the BHB and the BS

with hydrogen lines sensitive to the surface gravity:

  • Balmer lines : u-band
  • Paschen lines : z-band
slide-18
SLIDE 18

The BHB selection with CFIS

  • Disentangle the BHB and the BS

with hydrogen lines sensitive to the surface gravity:

  • Balmer lines : u-band
  • Paschen lines : z-band
  • Bell et al., 2010 :

74% pure and 57% complete (ugr)

  • Vickers et al., 2012 :

77% pure and 51% complete (griz)

slide-19
SLIDE 19

The BHB selection with CFIS

  • Disentangle the BHB and the BS

with hydrogen lines sensitive to the surface gravity:

  • Balmer lines : u-band
  • Paschen lines : z-band
  • Bell et al., 2010 :

74% pure and 57% complete (ugr)

  • Vickers et al., 2012 :

77% pure and 51% complete (griz)

  • Principal Component Analysis (PCA)
  • 75% pure and 71% complete
  • Deeper thanks to CFIS (up to 240 kpc)

u + griz CFIS PS1

−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

slide-20
SLIDE 20

The BHB selection with CFIS

  • Disentangle the BHB and the BS

with hydrogen lines sensitive to the surface gravity:

  • Balmer lines : u-band
  • Paschen lines : z-band
  • Bell et al., 2010 :

74% pure and 57% complete (ugr)

  • Vickers et al., 2012 :

77% pure and 51% complete (griz)

  • Principal Component Analysis (PCA)
  • 75% pure and 71% complete
  • Deeper thanks to CFIS (up to 240 kpc)

u + griz CFIS PS1

−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)

slide-21
SLIDE 21

The BHB stars Radial profile of the stellar halo

slide-22
SLIDE 22

The radial profile of the BHBs

  • Find the radial profile that

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)

slide-23
SLIDE 23

The radial profile of the BHBs

  • Find the radial profile that

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)

slide-24
SLIDE 24

The radial profile of the BHBs

  • Find the radial profile that

fit the BHB distribution :

Sarea(l, b) = ⇢ 1 if (l, b) in CFIS 0 otherwise

  • Selection function :
  • Footprint

Souter halo(RGC) = ⇢ 1 if 15 < RGC < 240 kpc 0 otherwise

  • Distance (15< RGC < 240 kpc)

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

  • Contamination of

compact objects

8 9 10 11 12 13 14 5.A. (deg) 38 39 40 41 42 43 44
  • Dec. (deg)
22.5 23.0 23.5 24.0 24.5 5.A. (deg) 29.5 30.0 30.5 31.0 31.5 32.0

M 31 M33

SSgr(l, b) = ⇢ 0 if | ˜ B| < 10.0 deg 1 otherwise

  • Presence of the Sgr stream

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) (

  • Completeness
slide-25
SLIDE 25

The radial profile of the BHBs

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:

  • Constant oblateness
  • q(RGC) :

γ = 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

slide-26
SLIDE 26

The radial profile of the BHBs

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:

  • Constant oblateness
  • q(RGC) :

γ = 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)

slide-27
SLIDE 27

Comparison with previous results

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)

slide-28
SLIDE 28

Comparison with previous results

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

slide-29
SLIDE 29

Comparison with previous results

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

slide-30
SLIDE 30

Conclusions

  • Identified the BHBs with an

unprecedented precision by their photometry

  • Radial profile of the stellar halo

traced by the BHB is similar to the halo traced by the giants …

  • Traced the radial profile of the

stellar halo up to ~240 kpc

  • … but shallower than traced by the

RR Lyrae

slide-31
SLIDE 31

Conclusions

  • Identified the BHBs with an

unprecedented precision by their photometry

  • Radial profile of the stellar halo

traced by the BHB is similar to the halo traced by the giants …

  • Traced the radial profile of the

stellar halo up to ~240 kpc

  • … but shallower than traced by the

RR Lyrae

  • Compare the distribution of these

different populations with cosmological simulations (Auriga?) and/or the Galactic Besançon Model

  • Study the dynamics of the BHBs

=> Mass of the MW and the 3D distribution of the DM halo at large distances

Future works

slide-32
SLIDE 32

Completeness of the sample

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

  • Assume that the limiting

magnitude is not the u-band

  • Studied the completeness of the

griz bands via HSC-SSP

slide-33
SLIDE 33

Completeness of the sample

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

  • Assume that the limiting

magnitude is not the u-band

  • Studied the completeness of the

griz bands via HSC-SSP Limiting magnitude = z-band

slide-34
SLIDE 34

Completeness of the sample

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

  • Assume that the limiting

magnitude is not the u-band

  • Studied the completeness of the

griz bands via HSC-SSP Limiting magnitude = z-band

  • Used the Luminosity

Function to see the spatial variation of the completeness Field of reference