The total stellar halo mass of the Milky Way

We measure the total stellar halo luminosity using red giant branch (RGB) stars selected from Gaia data release 2. Using slices in magnitude, colour, and location on the sky, we decompose RGB stars belonging to the disc and halo by fitting two-dimensional Gaussians to the Galactic proper motion distributions. The number counts of RGB stars are converted to total stellar halo luminosity using a suite of isochrones weighted by age and metallicity, and by applying a volume correction based on the stellar halo density profile. Our method is tested and calibrated using Galaxia and N-body models. We find a total luminosity (out to 100 kpc) of Lhalo = 7.9 ± 2.0 × 108 L excluding Sgr, and Lhalo = 9.4 ± 2.4 × 108 L including Sgr. These values are appropriate for our adopted stellar halo density profile and metallicity distribution, but additional systematics related to these assumptions are quantified and discussed. Assuming a stellar mass-to-light ratio appropriate for a Kroupa initial mass function (M/L = 1.5), we estimate a stellar halo mass of M halo = 1.4 ± 0.4 × 109 M. This mass is larger than previous estimates in the literature, but is in good agreement with the emerging picture that the (inner) stellar halo is dominated by one massive dwarf progenitor. Finally, we argue that the combination of a ∼109 M mass and an average metallicity of [Fe/H]∼−1.5 for the Galactic halo points to an ancient (∼10 Gyr) merger event

The Imprint of Cosmic Reionization on the Luminosity Function of Galaxies

The (re)ionization of hydrogen in the early universe has a profound effect on the formation of the first galaxies: by raising the gas temperature and pressure, it prevents gas from cooling into small halos, thus affecting the abundance of present-day small galaxies. Using the Galform semi-analytic model of galaxy formation, we show that two key aspects of the reionization process—when reionization takes place and the characteristic scale below which it suppresses galaxy formation—are imprinted in the luminosity function of dwarf galaxies. We focus on the luminosity function of satellites of galaxies like the Milky Way and the LMC, which is easier to measure than the luminosity function of the dwarf population as a whole. Our results show that the details of these two characteristic properties of reionization determine the shape of the luminosity distribution of satellites in a unique way, and are largely independent of the other details of the galaxy formation model. Our models generically predict a bimodality in the distribution of satellites as a function of luminosity: a population of faint satellites and population of bright satellites separated by a "valley" forged by reionization. We show that this bimodal distribution is present at high statistical significance in the combined satellite luminosity function of the Milky Way and M31. We make predictions for the expected number of satellites around LMC-mass dwarfs where the bimodality may also be measurable in future observational programs. Our preferred model predicts a total of 26 ± 10 (68% confidence) satellites brighter than M V = 0 in LMC-mass systems

Can we really pick and choose? Benchmarking various selections of Gaia Enceladus/Sausage stars in observations with simulations

Large spectroscopic surveys plus Gaia astrometry have shown us that the inner stellar halo of the Galaxy is dominated by the debris of Gaia Enceladus/Sausage (GES). With the richness of data at hand, there are a myriad of ways these accreted stars have been selected. We investigate these GES selections and their effects on the inferred progenitor properties using data constructed from APOGEE and Gaia. We explore selections made in eccentricity, energy-angular momentum (E-Lz), radial action-angular momentum (Jr-Lz), action diamond, and [Mg/Mn]-[Al/Fe] in the observations, selecting between 144 and 1279 GES stars with varying contamination from in-situ and other accreted stars. We also use the Auriga cosmological hydrodynamic simulations to benchmark the different GES dynamical selections. Applying the same observational GES cuts to nine Auriga galaxies with a GES, we find that the Jr-Lz method is best for sample purity and the eccentricity method for completeness. Given the average metallicity of GES (−1.28 < [Fe/H] < −1.18), we use the z = 0 mass–metallicity relationship to find an average of ∼4 × 108 M⊙. We adopt a similar procedure and derive for the GES-like systems in Auriga and find that the eccentricity method overestimates the true by ∼2.6 × while E-Lz underestimates by ∼0.7 ×. Lastly, we estimate the total mass of GES to be using the relationship between the metallicity gradient and the GES-to-in-situ energy ratio. In the end, we cannot just ‘pick and choose’ how we select GES stars, and instead should be motivated by the science question

Can we really pick and choose? Benchmarking various selections of Gaia Enceladus/Sausage stars in observations with simulations

Large spectroscopic surveys plus Gaia astrometry have shown us that the inner stellar halo of the Galaxy is dominated by the debris of Gaia Enceladus/Sausage (GES). With the richness of data at hand, there are a myriad of ways these accreted stars have been selected. We investigate these GES selections and their effects on the inferred progenitor properties using data constructed from APOGEE and Gaia. We explore selections made in eccentricity, energy-angular momentum (E-Lz), radial action-angular momentum (Jr-Lz), action diamond, and [Mg/Mn]-[Al/Fe] in the observations, selecting between 144 and 1,279 GES stars with varying contamination from in-situ and other accreted stars. We also use the Auriga cosmological hydrodynamic simulations to benchmark the different GES dynamical selections. Applying the same observational GES cuts to nine Auriga galaxies with a GES, we find that the Jr-Lz method is best for sample purity and the eccentricity method for completeness. Given the average metallicity of GES (-1.28 < [Fe/H] < -1.18), we use the $z=0$ mass-metallicity relationship to find an average $\rm M_{\star}$ of $\sim 4 \times 10^{8}$ $\rm M_{\odot}$. We adopt a similar procedure and derive $\rm M_{\star}$ for the GES-like systems in Auriga and find that the eccentricity method overestimates the true $\rm M_{\star}$ by $\sim2.6\times$ while E-Lz underestimates by $\sim0.7\times$. Lastly, we estimate the total mass of GES to be $\rm 10^{10.5 - 11.1}~M_{\odot}$ using the relationship between the metallicity gradient and the GES-to-in-situ energy ratio. In the end, we cannot just `pick and choose' how we select GES stars, and instead should be motivated by the science question.Comment: 20 pages, 14 figures, submitted to MNRA

The effects of dynamical substructure on Milky Way mass estimates from the high velocity tail of the local stellar halo

We investigate the impact of dynamical streams and substructure on estimates of the local escape speed and total mass of Milky Way-mass galaxies from modelling the high velocity tail of local halo stars. We use a suite of high-resolution, magneto-hydrodynamical cosmological zoom-in simulations, which resolve phase space substructure in local volumes around solar-like positions. We show that phase space structure varies significantly between positions in individual galaxies and across the suite. Substructure populates the high velocity tail unevenly and leads to discrepancies in the mass estimates. We show that a combination of streams, sample noise and truncation of the high velocity tail below the escape speed leads to a distribution of mass estimates with a median that falls below the true value by $\sim 20 \%$, and a spread of a factor of 2 across the suite. Correcting for these biases, we derive a revised value for the Milky Way mass presented in Deason et al. of $1.29 ^{+0.37}_{-0.47} \times 10^{12}$ $\rm M_{\odot}$.Comment: Re-submitted to MNRAS Letters after minor revisio

The effects of dynamical substructure on Milky Way mass estimates from the high-velocity tail of the local stellar halo

We investigate the impact of dynamical streams and substructure on estimates of the local escape speed and total mass of Milky-Way-mass galaxies from modelling the high-velocity tail of local halo stars. We use a suite of high-resolution magnetohydrodynamical cosmological zoom-in simulations that resolve phase space substructure in local volumes around solar-like positions. We show that phase space structure varies significantly between positions in individual galaxies and across the suite. Substructure populates the high-velocity tail unevenly and leads to discrepancies in the mass estimates. We show that a combination of streams, sample noise, and truncation of the high-velocity tail below the escape speed leads to a distribution of mass estimates with a median that falls below the true value by ∼20 per cent ∼20 per cent ⁠, and a spread of a factor of 2 across the suite. Correcting for these biases, we derive a revised value for the Milky Way mass presented in Deason et al. of 1.29 +0.37 −0.47 × 10 12 M ⊙ 1.29−0.47+0.37×1012M⊙ ⁠

HALO7D III: Chemical Abundances of Milky Way Halo Stars from Medium Resolution Spectra

The Halo Assembly in Lambda Cold Dark Matter: Observations in 7 Dimensions (HALO7D) survey measures the kinematics and chemical properties of stars in the Milky Way (MW) stellar halo to learn about the formation of our Galaxy. HALO7D consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of MW halo main sequence turn-off (MSTO) stars in the four CANDELS fields. HALO7D consists of deep pencil beams, making it complementary to other contemporary wide-field surveys. We present the [Fe/H] and [alpha/Fe] abundances for 113 HALO7D stars in the Galactocentric radial range of $\sim10-40$ kpc. Using the full 7D chemodynamical data (3D positions, 3D velocities, and abundances) of HALO7D, we measure the velocity anisotropy, $\beta$, of the halo velocity ellipsoid for each field and for different metallicity-binned subsamples. We find that two of the four fields have stars on very radial orbits while the remaining two have stars on more isotropic orbits. Separating the stars into high, mid, and low [Fe/H] bins at $-2.2$ dex and $-1.1$ dex for each field separately, we find differences in the anisotropies between the fields and between the bins; some fields appear dominated by radial orbits in all bins while other fields show variation between the [Fe/H] bins. These chemodynamical differences are evidence that the HALO7D fields have different fractional contributions from the progenitors that built up the MW stellar halo. Our results highlight the additional information that is available on smaller spatial scales when compared to results from a spherical average of the stellar halo.Comment: 32 pages, 15 figure

Constraining the shape of dark matter haloes with globular clusters and diffuse stellar light in the E-MOSAICS simulations

We explore how diffuse stellar light and globular clusters (GCs) can be used to trace the matter distribution of their host halo using an observational methodology. For this, we use 117 simulated dark matter (DM) haloes from the periodic volume of the E-MOSAICS project. For each halo, we compare the stellar surface brightness and GC projected number density maps to the surface density of DM. We find that the dominant structures identified in the stellar light and GCs correspond closely with those from the DM. Our method is unaffected by the presence of satellites and its precision improves with fainter GC samples. We recover tight relations between the dimensionless profiles of stellar-to-DM surface density and GC-to-DM surface density, suggesting that the profile of DM can be accurately recovered from the stars and GCs (σ ≤ 0.5 dex). We quantify the projected morphology of DM, stars, and GCs and find that the stars and GCs are more flattened than the DM. Additionally, the semimajor axes of the distribution of stars and GCs are typically misaligned by ∼10 degrees from that of DM. We demonstrate that deep imaging of diffuse stellar light and GCs can place constraints on the shape, profile, and orientation of their host halo. These results extend down to haloes with central galaxies M⋆ ≥ 1010 M⊙, and the analysis will be applicable to future data from the Euclid, Roman, and the Rubin observatorie