skip to main content


Title: The origin of metal-poor stars on prograde disc orbits in FIRE simulations of Milky Way-mass galaxies
ABSTRACT In hierarchical structure formation, metal-poor stars in and around the Milky Way (MW) originate primarily from mergers of lower mass galaxies. A common expectation is therefore that metal-poor stars should have isotropic, dispersion-dominated orbits that do not correlate strongly with the MW disc. However, recent observations of stars in the MW show that metal-poor ($\rm {[Fe/H]}\lesssim -2$) stars are preferentially on prograde orbits with respect to the disc. Using the Feedback In Realistic Environments 2 (FIRE-2) suite of cosmological zoom-in simulations of MW/M31-mass galaxies, we investigate the prevalence and origin of prograde metal-poor stars. Almost all (11 of 12) of our simulations have metal-poor stars on preferentially prograde orbits today and throughout most of their history: we thus predict that this is a generic feature of MW/M31-mass galaxies. The typical prograde-to-retrograde ratio is ∼2:1, which depends weakly on stellar metallicity at $\rm {[Fe/H]}\lesssim -1$. These trends predicted by our simulations agree well with MW observations. Prograde metal-poor stars originate largely from a single Large/Small Magellanic Cloud (LMC/SMC)-mass gas-rich merger $7\!-\!12.5\, \rm {Gyr}$ ago, which deposited existing metal-poor stars and significant gas on an orbital vector that sparked the formation of and/or shaped the orientation of a long-lived stellar disc, giving rise to a prograde bias for all low-metallicity stars. We find subdominant contributions from in situ stars formed in the host galaxy before this merger, and in some cases, additional massive mergers. We find few clear correlations between any properties of our MW/M31-mass galaxies at z = 0 and the degree of this prograde bias as a result of diverse merger scenarios.  more » « less
Award ID(s):
1715216
NSF-PAR ID:
10278874
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
505
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
921 to 938
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We use FIRE-2 simulations to examine 3D variations of gas-phase elemental abundances of [O/H], [Fe/H], and [N/H] in 11 MW and M31-mass galaxies across their formation histories at z ≤ 1.5 ($t_{\rm lookback} \le 9.4 \, \rm {Gyr}$), motivated by characterizing the initial conditions of stars for chemical tagging. Gas within $1 \, \rm {kpc}$ of the disc mid-plane is vertically homogeneous to $\lesssim 0.008 \, \rm {dex}$ at all z ≤ 1.5. We find negative radial gradients (metallicity decreases with galactocentric radius) at all times, which steepen over time from $\approx \! -0.01 \, \rm {dex}\, \rm {kpc}^{-1}$ at z = 1 ($t_{\rm lookback} = 7.8 \, \rm {Gyr}$) to $\approx \! -0.03 \, \rm {dex}\, \rm {kpc}^{-1}$ at z = 0, and which broadly agree with observations of the MW, M31, and nearby MW/M31-mass galaxies. Azimuthal variations at fixed radius are typically $0.14 \, \rm {dex}$ at z = 1, reducing to $0.05 \, \rm {dex}$ at z = 0. Thus, over time radial gradients become steeper while azimuthal variations become weaker (more homogeneous). As a result, azimuthal variations were larger than radial variations at z ≳ 0.8 ($t_{\rm lookback} \gtrsim 6.9 \, \rm {Gyr}$). Furthermore, elemental abundances are measurably homogeneous (to ≲0.05 dex) across a radial range of $\Delta R \approx 3.5 \, \rm {kpc}$ at z ≳ 1 and $\Delta R \approx 1.7 \, \rm {kpc}$ at z = 0. We also measure full distributions of elemental abundances, finding typically negatively skewed normal distributions at z ≳ 1 that evolve to typically Gaussian distributions by z = 0. Our results on gas abundances inform the initial conditions for stars, including the spatial and temporal scales for applying chemical tagging to understand stellar birth in the MW. 
    more » « less
  2. ABSTRACT

    We characterize the 3D spatial variations of [Fe/H], [Mg/H], and [Mg/Fe] in stars at the time of their formation, across 11 simulated Milky Way (MW)- and M31-mass galaxies in the FIRE-2 simulations, to inform initial conditions for chemical tagging. The overall scatter in [Fe/H] within a galaxy decreased with time until $\approx 7 \, \rm {Gyr}$ ago, after which it increased to today: this arises from a competition between a reduction of azimuthal scatter and a steepening of the radial gradient in abundance over time. The radial gradient is generally negative, and it steepened over time from an initially flat gradient $\gtrsim 12 \, \rm {Gyr}$ ago. The strength of the present-day abundance gradient does not correlate with when the disc ‘settled’; instead, it best correlates with the radial velocity dispersion within the galaxy. The strength of azimuthal variation is nearly independent of radius, and the 360 deg scatter decreased over time, from $\lesssim 0.17 \, \rm {dex}$ at $t_{\rm lb} = 11.6 \, \rm {Gyr}$ to $\sim 0.04 \, \rm {dex}$ at present-day. Consequently, stars at $t_{\rm lb} \gtrsim 8 \, \rm {Gyr}$ formed in a disc with primarily azimuthal scatter in abundances. All stars formed in a vertically homogeneous disc, Δ[Fe/H]$\le 0.02 \, \rm {dex}$ within $1 \, \rm {kpc}$ of the galactic mid-plane, with the exception of the young stars in the inner $\approx 4 \, \rm {kpc}$ at z ∼ 0. These results generally agree with our previous analysis of gas-phase elemental abundances, which reinforces the importance of cosmological disc evolution and azimuthal scatter in the context of stellar chemical tagging. We provide analytic fits to our results for use in chemical-tagging analyses.

     
    more » « less
  3. null (Ed.)
    ABSTRACT Carbon enhanced metal poor (CEMP)-no stars, a subset of CEMP stars ($\rm [C/Fe]\ge 0.7$ and $\rm [Fe/H]\lesssim -1$) have been discovered in ultra-faint dwarf (UFD) galaxies, with $M_{\rm vir}\approx 10^8{\, \mathrm{ M}_\odot }$ and $M_{\ast }\approx 10^3-10^4{\, \mathrm{ M}_\odot }$ at z = 0, as well as in the halo of the Milky Way (MW). These CEMP-no stars are local fossils that may reflect the properties of the first (Pop III) and second (Pop II) generation of stars. However, cosmological simulations have struggled to reproduce the observed level of carbon enhancement of the known CEMP-no stars. Here, we present new cosmological hydrodynamic zoom-in simulations of isolated UFDs that achieve a gas mass resolution of $m_{\rm gas}\approx 60{\, \mathrm{ M}_\odot }$. We include enrichment from Pop III faint supernovae (SNe), with ESN = 0.6 × 1051 erg, to understand the origin of CEMP-no stars. We confirm that Pop III and Pop II stars are mainly responsible for the formation of CEMP and C-normal stars, respectively. New to this study, we find that a majority of CEMP-no stars in the observed UFDs and the MW halo can be explained by Pop III SNe with normal explosion energy (ESN = 1.2 × 1051 erg) and Pop II enrichment, but faint SNe might also be needed to produce CEMP-no stars with $\rm [C/Fe]\gtrsim 2$, corresponding to the absolute carbon abundance of $\rm A(C)\gtrsim 6.0$. Furthermore, we find that while we create CEMP-no stars with high carbon ratio $\rm [C/Fe]\approx 3-4$, by adopting faint SNe, it is still challenging to reproduce CEMP-no stars with extreme level of carbon abundance of $\rm A(C)\approx 7.0-7.5$, observed both in the MW halo and UFDs. 
    more » « less
  4. ABSTRACT

    We investigate the formation (spin-up) of galactic discs in the artemis simulations of Milky Way (MW)-mass galaxies. In almost all galaxies, discs spin up at higher [Fe/H] than the MW. Those galaxies that contain an analogue of the Gaia Sausage-Enceladus (GSE) spin up at a lower average metallicity than those without. We identify six galaxies with spin-up metallicity similar to that of the MW, which formed their discs ∼8–11 Gyr ago. Five of these experience a merger similar to the GSE. The spin-up times correlate with the halo masses at early times: galaxies with early spin-up have larger virial masses at a lookback time tL = 12 Gyr. The fraction of stars accreted from outside the host galaxy is smaller in galaxies with earlier spin-ups. Accreted fractions small enough to be comparable to the MW are only found in galaxies with the earliest disc formation and large initial virial masses (M200c ≈ 2 × 1011 M⊙ at tL = 12 Gyr). We find that discs form when the halo’s virial mass reaches a threshold of M200c ≈ (6 ± 3) × 1011 M⊙, independent of the spin-up time. However, the failure to form a disc in other galaxies appears to be instead related to mergers at early times. We also find that discs form when the central potential is not particularly steep. Our results indicate that the MW assembled its mass and formed its disc earlier than the average galaxy of a similar mass.

     
    more » « less
  5. null (Ed.)
    ABSTRACT In this work, we combine spectroscopic information from the SkyMapper survey for Extremely Metal-Poor stars and astrometry from Gaia DR2 to investigate the kinematics of a sample of 475 stars with a metallicity range of $-6.5 \le \rm [Fe/H] \le -2.05$ dex. Exploiting the action map, we identify 16 and 40 stars dynamically consistent with the Gaia Sausage and Gaia Sequoia accretion events, respectively. The most metal poor of these candidates have metallicities of $\rm [Fe/H]=-3.31\, \mathrm{ and }\, -3.74$, respectively, helping to define the low-metallicity tail of the progenitors involved in the accretion events. We also find, consistent with other studies, that ∼21 per cent of the sample have orbits that remain confined to within 3 kpc of the Galactic plane, that is, |Zmax| ≤ 3 kpc. Of particular interest is a subsample (∼11 per cent of the total) of low |Zmax| stars with low eccentricities and prograde motions. The lowest metallicity of these stars has [Fe/H] = –4.30 and the subsample is best interpreted as the very low-metallicity tail of the metal-weak thick disc population. The low |Zmax|, low eccentricity stars with retrograde orbits are likely accreted, while the low |Zmax|, high eccentricity pro- and retrograde stars are plausibly associated with the Gaia Sausage system. We find that a small fraction of our sample (∼4 per cent of the total) is likely escaping from the Galaxy, and postulate that these stars have gained energy from gravitational interactions that occur when infalling dwarf galaxies are tidally disrupted. 
    more » « less