No Arabic abstract
We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at $z=0$.In the outer halo, outflow stars are more spatially homogeneous, metal rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to $sim 10 %$ of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as $sim 40%$ for the most metal-rich local halo stars ([Fe/H] $> -0.5$). We conclude that the Milky Way stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell $&$ Sandange, who suggested that halo stars formed in rapidly $textit{infalling}$ gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme velocity stars in 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 2-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 L_halo = 7.9 +/- 2.0 x 10^8 L_Sun excluding Sgr, and L_halo = 9.4 +/- 2.4 x 10^8 L_Sun 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 x 10^9 M_Sun. 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 ~10^9 M_Sun mass and an average metallicity of <[Fe/H]> ~ -1.5 for the Galactic halo points to an ancient (~10 Gyr) merger event.
In the $Gaia$ era stellar kinematics are extensively used to study Galactic halo stellar populations, to search for halo structures, and to characterize the interface between the halo and hot disc populations. We use distribution function-based models of modern datasets with 6D phase space data to qualitatively describe a variety of kinematic spaces commonly used in the study of the Galactic halo. Furthermore, we quantitatively assess how well each kinematic space can separate radially anisotropic from isotropic halo populations. We find that scaled action space (the ``action diamond) is superior to other commonly used kinematic spaces at this task. We present a new, easy to implement selection criterion for members of the radially-anisotropic $Gaia$-Enceladus merger remnant, which we find achieves a sample purity of 82 per cent in our models with respect to contamination from the more isotropic halo. We compare this criterion to literature criteria, finding that it produces the highest purity in the resulting samples, at the expense of a modest reduction in completeness. We also show that selection biases that underlie nearly all contemporary spectroscopic datasets can noticeably impact the $E-L_{z}$ distribution of samples in a manner that may be confused for real substructure. We conclude by providing recommendations for how authors should use stellar kinematics in the future to study the Galactic stellar halo.
We present chemical abundances of 57 metal-poor stars that are likely constituents of the outer stellar halo in the Milky Way. Almost all of the sample stars have an orbit reaching a maximum vertical distance (Z_max) of >5 kpc above and below the Galactic plane. High-resolution, high signal-to-noise spectra for the sample stars obtained with Subaru/HDS are used to derive chemical abundances of Na, Mg, Ca, Ti, Cr, Mn, Fe, Ni, Zn, Y and Ba with an LTE abundance analysis code. The resulting abundance data are combined with those presented in literature that mostly targeted at smaller Z_max stars, and both data are used to investigate any systematic trends in detailed abundance patterns depending on their kinematics. It is shown that, in the metallicity range of -2<[Fe/H]<-1, the [Mg/Fe] ratios for the stars with Z_max>5 kpc are systematically lower (~0.1 dex) than those with smaller Z_max. This result of the lower [alpha/Fe] for the assumed outer halo stars is consistent with previous studies that found a signature of lower [alpha/Fe] ratios for stars with extreme kinematics. A distribution of the [Mg/Fe] ratios for the outer halo stars partly overlaps with that for stars belonging to the Milky Way dwarf satellites in the metallicity interval of -2<[Fe/H]<-1 and spans a range intermediate between the distributions for the inner halo stars and the stars belonging to the satellites. Our results confirm inhomogeneous nature of chemical abundances within the Milky Way stellar halo depending on kinematic properties of constituent stars as suggested by earlier studies. Possible implications for the formation of the Milky Way halo and its relevance to the suggested dual nature of the halo are discussed.
We analyse the structure of the local stellar halo of the Milky Way using $sim$ 60000 stars with full phase space coordinates extracted from the SDSS--{it Gaia} catalogue. We display stars in action space as a function of metallicity in a realistic axisymmetric potential for the Milky Way Galaxy. The metal-rich population is more distended towards high radial action $J_R$ as compared to azimuthal or vertical action, $J_phi$ or $J_z$. It has a mild prograde rotation $(langle v_phi rangle approx 25$ km s$^{-1}$), is radially anisotropic and highly flattened with axis ratio $q approx 0.6 - 0.7$. The metal-poor population is more evenly distributed in all three actions. It has larger prograde rotation $(langle v_phi rangle approx 50$ km s$^{-1}$), a mild radial anisotropy and a roundish morphology ($qapprox 0.9$). We identify two further components of the halo in action space. There is a high energy, retrograde component that is only present in the metal-rich stars. This is suggestive of an origin in a retrograde encounter, possibly the one that created the stripped dwarf galaxy nucleus, $omega$Centauri. Also visible as a distinct entity in action space is a resonant component, which is flattened and prograde. It extends over a range of metallicities down to [Fe/H] $approx -3$. It has a net outward radial velocity $langle v_R rangle approx 12$ km s$^{-1}$ within the Solar circle at $|z| <3.5$ kpc. The existence of resonant stars at such extremely low metallicities has not been seen before.
We report on the global structure of the Milky Way (MW) stellar halo up to its outer boundary based on the analysis of blue-horizontal branch stars (BHBs). These halo tracers are extracted from the $(g,r,i,z)$ band multi-photometry in the internal data release of the on-going Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) surveyed over $sim550$~deg$^2$ area. In order to select most likely BHBs by removing blue straggler stars (BSs) and other contamination in a statistically significant manner, we have developed and applied an extensive Bayesian method, instead of the simple color cuts adopted in our previous work, where each of the template BHBs and non-BHBs obtained from the available catalogs is represented as a mixture of multiple Gaussian distributions in the color-color diagrams. We found from the candidate BHBs in the range of 18.5<g<23.5 mag that the radial density distribution over a Galactocentric radius of r=36-360 kpc can be approximated as a single power-law profile with an index of $alpha=3.74^{+0.21}_{-0.22}$ or a broken power-law profile with an index of $alpha_{rm in}=2.92^{+0.33}_{-0.33}$ at $r$ below a broken radius of $r_{rm b}=160^{+18}_{-19}$~kpc and a very steep slope of $alpha_{rm out}=15.0^{+3.7}_{-4.5}$ at $r>r_{rm b}$. The latter profile with a prolate shape having an axial ratio of $q=1.72^{+0.44}_{-0.28}$ is most likely and this halo may hold a rather sharp boundary at r=160kpc. The slopes of the halo density profiles are compared with those from the suite of hydrodynamical simulations for the formation of stellar halos. This comparison suggests that the MW stellar halo may consist of the two overlapping components: the in situ. inner halo as probed by RR Lyrae stars showing a relatively steep radial density profile and the ex situ. outer halo with a shallow profile probed by BHBs here, which is made by accretion of small stellar systems.