No Arabic abstract
We use the Bullock & Johnston suite of simulations to study the density profiles of L*-type galaxy stellar haloes. Observations of the Milky Way and M31 stellar haloes show contrasting results: the Milky Way has a `broken profile, where the density falls off more rapidly beyond ~ 25 kpc, while M31 has a smooth profile out to 100 kpc with no obvious break. Simulated stellar haloes, built solely by the accretion of dwarf galaxies, also exhibit this behavior: some haloes have breaks, while others dont. The presence or absence of a break in the stellar halo profile can be related to the accretion history of the galaxy. We find that a break radius is strongly related to the build up of stars at apocentres. We relate these findings to observations, and find that the `break in the Milky Way density profile is likely associated with a relatively early (~ 6-9 Gyr ago) and massive accretion event. In contrast, the absence of a break in the M31 stellar halo profile suggests that its accreted satellites have a wide range of apocentres. Hence, it is likely that M31 has had a much more prolonged accretion history than the Milky Way.
We introduce the ARTEMIS simulations, a new set of 42 zoomed-in, high-resolution (baryon particle mass of ~ 2x10^4 Msun/h), hydrodynamical simulations of galaxies residing in haloes of Milky Way mass, simulated with the EAGLE galaxy formation code with re-calibrated stellar feedback. In this study, we analyse the structure of stellar haloes, specifically the mass density, surface brightness, metallicity, colour and age radial profiles, finding generally very good agreement with recent observations of local galaxies. The stellar density profiles are well fitted by broken power laws, with inner slopes of ~ -3, outer slopes of ~ -4 and break radii that are typically ~ 20-40 kpc. The break radii generally mark the transition between in situ formation and accretion-driven formation of the halo. The metallicity, colour and age profiles show mild large-scale gradients, particularly when spherically-averaged or viewed along the major axes. Along the minor axes, however, the profiles are nearly flat, in agreement with observations. Overall, the structural properties can be understood by two factors: that in situ stars dominate the inner regions and that they reside in a spatially-flattened distribution that is aligned with the disc. Observations targeting both the major and minor axes of galaxies are thus required to obtain a complete picture of stellar haloes.
We investigate the chemical and kinematic properties of the diffuse stellar haloes of six simulated Milky Way-like galaxies from the Aquarius Project. Binding energy criteria are adopted to defined two dynamically distinct stellar populations: the diffuse inner and outer haloes, which comprise different stellar sub-populations with particular chemical and kinematic characteristics. Our simulated inner- and outer-halo stellar populations have received contributions from debris stars (formed in sub-galactic systems while they were outside the virial radius of the main progenitor galaxies) and endo-debris stars (those formed in gas-rich sub-galactic systems inside the dark matter haloes). The inner haloes possess an additional contribution from disc-heated stars in the range $sim 3 - 30 %$, with a mean of $sim 20% $. Disc-heated stars might exhibit signatures of kinematical support, in particular among the youngest ones. Endo-debris plus disc-heated stars define the so-called insitu stellar populations. In both the inner- and outer-halo stellar populations, we detect contributions from stars with moderate to low [$alpha$/Fe] ratios, mainly associated with the endo-debris or disc-heated sub-populations. The observed abundance gradients in the inner-halo regions are influenced by both the level of chemical enrichment and the relative contributions from each stellar sub-population. Steeper abundance gradients in the inner-halo regions are related to contributions from the disc-heated and endo-debris stars, which tend to be found at lower binding energies than debris stars. (Abridged).
We have examined the resolved stellar populations at large galactocentric distances along the minor axis (from 10 kpc up to between 40 and 75 kpc), with limited major axis coverage, of six nearby highly-inclined Milky Way-mass disc galaxies using HST data from the GHOSTS survey. We select red giant branch stars to derive stellar halo density profiles. The projected minor axis density profiles can be approximated by power laws with projected slopes of between $-2$ and $-3.7$ and a diversity of stellar halo masses of $1-6times 10^{9}M_{odot}$, or $2-14%$ of the total galaxy stellar masses. The typical intrinsic scatter around a smooth power law fit is $0.05-0.1$ dex owing to substructure. By comparing the minor and major axis profiles, we infer projected axis ratios $c/a$ at $sim 25$ kpc between $0.4-0.75$. The GHOSTS stellar haloes are diverse, lying between the extremes charted out by the (rather atypical) haloes of the Milky Way and M31. We find a strong correlation between the stellar halo metallicities and the stellar halo masses. We compare our results with cosmological models, finding good agreement between our observations and accretion-only models where the stellar haloes are formed by the disruption of dwarf satellites. In particular, the strong observed correlation between stellar halo metallicity and mass is naturally reproduced. Low-resolution hydrodynamical models have unrealistically high stellar halo masses. Current high-resolution hydrodynamical models appear to predict stellar halo masses somewhat higher than observed but with reasonable metallicities, metallicity gradients and density profiles.
We study the formation of stellar haloes in three Milky Way-mass galaxies using cosmological smoothed particle hydrodynamics simulations, focusing on the subset of halo stars that form in situ, as opposed to those accreted from satellites. In situ stars in our simulations dominate the stellar halo out to 20 kpc and account for 30-40 per cent of its total mass. We separate in situ halo stars into three straightforward, physically distinct categories according to their origin: stars scattered from the disc of the main galaxy (heated disc), stars formed from gas smoothly accreted on to the halo (smooth gas) and stars formed in streams of gas stripped from infalling satellites (stripped gas). We find that most belong to the stripped gas category. Those originating in smooth gas outside the disc tend to form at the same time and place as the stripped-gas population, suggesting that their formation is associated with the same gas-rich accretion events. The scattered disc star contribution is negligible overall but significant in the Solar neighbourhood, where ~90 per cent of stars on eccentric orbits once belonged to the disc. However, the distinction between halo and thick disc in this region is highly ambiguous. The chemical and kinematic properties of the different components are very similar at the present day, but the global properties of the in situ halo differ substantially between the three galaxies in our study. In our simulations, the hierarchical buildup of structure is the driving force behind not only the accreted stellar halo, but also those halo stars formed in situ.
We present a comprehensive study of the chemical properties of the stellar haloes of Milky-Way mass galaxies, analysing the transition between the inner to the outer haloes. We find the transition radius between the relative dominance of the inner-halo and outer-halo stellar populations to be ~15-20 kpc for most of our haloes, similar to that inferred for the Milky Way from recent observations. While the number density of stars in the simulated inner-halo populations decreases rapidly with distance, the outer-halo populations contribute about 20-40 per cent in the fiducial solar neighborhood, in particular at the lowest metallicities. We have determined [Fe/H] profiles for our simulated haloes; they exhibit flat or mild gradients, in the range [-0.002, -0.01 ] dex/kpc. The metallicity distribution functions exhibit different features, reflecting the different assembly history of the individual stellar haloes. We find that stellar haloes formed with larger contributions from massive subgalactic systems have steeper metallicity gradients. Very metal-poor stars are mainly contributed to the halo systems by lower-mass satellites. There is a clear trend among the predicted metallicity distribution functions that a higher fraction of low-metallicity stars are found with increasing radius. These properties are consistent with the range of behaviours observed for stellar haloes of nearby galaxies.