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: Born this way: thin disc, thick disc, and isotropic spheroid formation in FIRE-2 Milky Way–mass galaxy simulations
ABSTRACT We investigate the formation of Milky Way–mass galaxies using FIRE-2 ΛCDM cosmological zoom-in simulations by studying the orbital evolution of stars formed in the main progenitor of the galaxy, from birth to the present day. We classify in situ stars as isotropic spheroid, thick-disc, and thin-disc according to their orbital circularities and show that these components are assembled in a time-ordered sequence from early to late times, respectively. All simulated galaxies experience an early phase of bursty star formation that transitions to a late-time steady phase. This transition coincides with the time that the inner CGM virializes. During the early bursty phase, galaxies have irregular morphologies and new stars are born on radial orbits; these stars evolve into an isotropic spheroidal population today. The bulk of thick-disc stars form at intermediate times, during a clumpy-disc ‘spin-up’ phase, slightly later than the peak of spheroid formation. At late times, once the CGM virializes and star formation ‘cools down,’ stars are born on circular orbits within a narrow plane. Those stars mostly inhabit thin discs today. Broadly speaking, stars with disc-like or spheroid-like orbits today were born that way. Mergers on to discs and secular processes do affect kinematics in our simulations, but play only secondary roles in populating thick-disc and in situ spheroid populations at z = 0. The age distributions of spheroid, thick disc, and thin disc populations scale self-similarly with the steady-phase transition time, which suggests that morphological age dating can be linked to the CGM virialization time in galaxies.  more » « less
Award ID(s):
2108962 1910346 1752913 2108318 1713353 2045928 1652522 2108230 2307327
PAR ID:
10428709
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
523
Issue:
4
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 6220-6238
Size(s):
p. 6220-6238
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. ABSTRACT Recent observations and simulations indicate substantial evolution in the properties of galaxies with time, wherein rotationally supported and steady thin discs (like those frequently observed in the local Universe) emerge from galaxies that are clumpy, irregular, and have bursty star formation rates (SFRs). To better understand the progenitors of local disc galaxies, we carry out an analysis of three FIRE-2 simulated galaxies with a mass similar to the Milky Way at redshift z = 0. We show that all three galaxies transition from bursty to steady SFRs at a redshift between z = 0.5 and z = 0.8, and that this transition coincides with the rapid (≲1 Gyr) emergence of a rotationally supported interstellar medium (ISM). In the late phase with steady SFR, the rotational energy comprises $${\gtrsim }90{{\ \rm per\ cent}}$$ of the total kinetic + thermal energy in the ISM, and is roughly half the gravitational energy. By contrast, during the early bursty phase, the ISM initially has a quasi-spheroidal morphology and its energetics are dominated by quasi-isotropic in- and outflows out of virial equilibrium. The subdominance of rotational support and out-of-equilibrium conditions at early times challenge the application of standard equilibrium disc models to high-redshift progenitors of Milky Way-like galaxies. We further find that the formation of a rotationally-supported ISM coincides with the onset of a thermal pressure supported inner circumgalactic medium (CGM). Before this transition, there is no clear boundary between the ISM and the inner CGM. 
    more » « less
  3. ABSTRACT We study the kinematics of stars both at their formation and today within 14 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. We quantify the relative importance of cosmological disc settling and post-formation dynamical heating. We identify three eras: a Pre-Disc Era (typically ≳ 8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disc Era (≈8–4 Gyr ago), when stars started to form on rotation-dominated orbits but with high velocity dispersion, σform; and a Late-Disc Era (≲ 4 Gyr ago), when stars formed with low σform. σform increased with time during the Pre-Disc Era, peaking ≈8 Gyr ago, then decreased throughout the Early-Disc Era as the disc settled and remained low throughout the Late-Disc Era. By contrast, the dispersion measured today, σnow, increases monotonically with age because of stronger post-formation heating for Pre-Disc stars. Importantly, most of σnow was in place at formation, not added post-formation, for stars younger than ≈10 Gyr. We compare the evolution of the three velocity components: at all times, σR, form > σϕ, form > σZ, form. Post-formation heating primarily increased σR at ages ≲ 4 Gyr but acted nearly isotropically for older stars. The kinematics of young stars in FIRE-2 broadly agree with the range observed across the MW, M31, M33, and PHANGS-MUSE galaxies. The lookback time that the disc began to settle correlates with its dynamical state today: earlier-settling galaxies currently form colder discs. Including stellar cosmic-ray feedback does not significantly change disc rotational support at fixed stellar mass. 
    more » « less
  4. ABSTRACT We study the growth of stellar discs of Milky Way-sized galaxies using a suite of cosmological simulations. We calculate the half-mass axis lengths and axis ratios of stellar populations split by age in galaxies with stellar mass $$M_{*}=10^7\!-\!10^{10}\, \mathrm{M}_{\odot }$$ at redshifts z > 1.5. We find that in our simulations stars always form in relatively thin discs, and at ages below 100 Myr are contained within half-mass height z1/2 ∼ 0.1 kpc and short-to-long axial ratio z1/2/x1/2 ∼ 0.15. Disc thickness increases with the age of stellar population, reaching median z1/2 ∼ 0.8 kpc and z1/2/x1/2 ∼ 0.6 for stars older than 500 Myr. We trace the same group of stars over the simulation snapshots and show explicitly that their intrinsic shape grows more spheroidal over time. We identify a new mechanism that contributes to the observed disc thickness: rapid changes in the orientation of the galactic plane mix the configuration of young stars. The frequently mentioned ‘upside-down’ formation scenario of galactic discs, which posits that young stars form in already thick discs at high redshift, may be missing this additional mechanism of quick disc inflation. The actual formation of stars within a fairly thin plane is consistent with the correspondingly flat configuration of dense molecular gas that fuels star formation. 
    more » « less
  5. null (Ed.)
    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