skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


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
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 Surveys of the Milky Way (MW) and M31 enable detailed studies of stellar populations across ages and metallicities, with the goal of reconstructing formation histories across cosmic time. These surveys motivate key questions for galactic archaeology in a cosmological context: When did the main progenitor of an MW/M31-mass galaxy form, and what were the galactic building blocks that formed it? We investigate the formation times and progenitor galaxies of MW/M31-mass galaxies using the Feedback In Realistic Environments-2 cosmological simulations, including six isolated MW/M31-mass galaxies and six galaxies in Local Group (LG)-like pairs at z = 0. We examine main progenitor ‘formation’ based on two metrics: (1) transition from primarily ex-situ to in-situ stellar mass growth and (2) mass dominance compared to other progenitors. We find that the main progenitor of an MW/M31-mass galaxy emerged typically at z ∼ 3–4 ($$11.6\!\!-\!\!12.2\, \rm {Gyr}$$ ago), while stars in the bulge region (inner 2 kpc) at z = 0 formed primarily in a single main progenitor at z ≲ 5 ($${\lesssim} \!12.6\, \rm {Gyr}$$ ago). Compared with isolated hosts, the main progenitors of LG-like paired hosts emerged significantly earlier (Δz ∼ 2, $$\Delta t\!\sim \!1.6\, \rm {Gyr}$$), with ∼4× higher stellar mass at all z ≳ 4 ($${\gtrsim} \!12.2\, \rm {Gyr}$$ ago). This highlights the importance of environment in MW/M31-mass galaxy formation, especially at early times. On average, about 100 galaxies with $$\rm {\it{ M}}_\rm {star}\!\gtrsim \!10^5\, \rm {M}_\odot$$ went into building a typical MW/M31-mass system. Thus, surviving satellites represent a highly incomplete census (by ∼5×) of the progenitor population. 
    more » « less
  3. 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
  4. 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
  5. ABSTRACT Recent observational studies have uncovered a small number of very metal-poor (VMP) stars with cold kinematics in the Galactic disc and bulge. However, their origins remain enigmatic. We select a total of 138 Milky Way (MW) analogues from the TNG50 cosmological simulation based on their z = 0 properties: discy morphology, stellar mass, and local environment. In order to make more predictive statements for the MW, we further limit the spatial volume coverage of stellar populations in galaxies to that targeted by the upcoming 4MOST high-resolution survey of the Galactic disc and bulge. We find that across all galaxies, ∼20 per cent of VMP ([Fe/H] < −2) stars belong to the disc, with some analogues reaching 30 per cent. About 50 ± 10 per cent of the VMP disc stars are, on average, older than 12.5 Gyr and ∼70 ± 10 per cent come from accreted satellites. A large fraction of the VMP stars belong to the halo (∼70) and have a median age of 12 Gyr. Our results with the TNG50 cosmological simulation confirm earlier findings with simulations of fewer individual galaxies, and suggest that the stellar disc of the MW is very likely to host significant amounts of very- and extremely-metal-poor stars that, although mostly of ex situ origin, can also form in situ, reinforcing the idea of the existence of a primordial Galactic disc. 
    more » « less