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. null (Ed.)
    ABSTRACT We investigate thin and thick stellar disc formation in Milky Way-mass galaxies using 12 FIRE-2 cosmological zoom-in simulations. All simulated galaxies experience an early period of bursty star formation that transitions to a late-time steady phase of near-constant star formation. Stars formed during the late-time steady phase have more circular orbits and thin-disc-like morphology at z = 0, while stars born during the bursty phase have more radial orbits and thick-disc structure. The median age of thick-disc stars at z = 0 correlates strongly with this transition time. We also find that galaxies with an earlier transition from bursty to steady star formation have a higher thin-disc fractions at z = 0. Three of our systems have minor mergers with Large Magellanic Cloud-size satellites during the thin-disc phase. These mergers trigger short starbursts but do not destroy the thin disc nor alter broad trends between the star formation transition time and thin/thick-disc properties. If our simulations are representative of the Universe, then stellar archaeological studies of the Milky Way (or M31) provide a window into past star formation modes in the Galaxy. Current age estimates of the Galactic thick disc would suggest that the Milky Way transitioned from bursty to steady phase ∼6.5 Gyr ago; prior to that time the Milky Way likely lacked a recognizable thin disc. 
    more » « less