An official website of the United States government
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
Rapid disc settling and the transition from bursty to steady star formation in Milky Way-mass galaxies
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
- PAR ID:
- 10389545
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 519
- Issue:
- 2
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 2598-2614
- Size(s):
- p. 2598-2614
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
ABSTRACT Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-α (Lyα), and Balmer line (Hα and Hβ) radiative transfer study of a high-resolution isolated Milky Way simulation using the state-of-the-art Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry accounting for molecular hydrogen, and dust grain evolution in the interstellar medium (ISM). We extend our publicly available Cosmic Lyα Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer and various methodology improvements for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust ($$f_\text{abs} \approx 27.5\,\rm{per\,\,cent}$$), helium ionization ($$f_\text{He} \approx 8.7\,\rm{per\,\,cent}$$), and anisotropic escape fractions ($$f_\text{esc} \approx 7.9\,\rm{per\,\,cent}$$), as these reduce the available budget for hydrogen line emission ($$f_\text{H} \approx 55.9\,\rm{per\,\,cent}$$). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas-phase structure, and spatial spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local Hα surveys, but would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Lyα emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve ≲10 pc multiphase ISM and ≲1 kpc CGM structures will provide crucial insights and predictions for current and upcoming Lyα observations.more » « less
-
Abstract The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM coevolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and turbulent kinetic energy of the CGM. Our model accounts for heating and turbulence driven by both supernova winds and cosmic accretion as well as radiative cooling, turbulence dissipation, and halo outflows due to CGM overpressurization. We demonstrate that, depending on parameters, the CGM can undergo a phase transition (“thermalization”) from a cool, turbulence-supported phase to a virial-temperature, thermally supported phase. This CGM phase transition is largely determined by the ability of radiative cooling to balance heating from supernova winds and turbulence dissipation. We perform an initial calibration of our model to the FIRE-2 cosmological hydrodynamical simulations and show that it can approximately reproduce the baryon cycles of the simulated halos. In particular, we find that, for these parameters, the phase transition occurs at high redshift in ultrafaint progenitors and at low redshift in classicalMvir∼ 1011M⊙dwarfs, while Milky Way–mass halos undergo the transition atz≈ 0.5. We see a similar transition in the simulations though it is more gradual, likely reflecting radial dependence and multiphase gas not captured by our model. We discuss these and other limitations of the model and possible future extensions.more » « less
-
null (Ed.)ABSTRACT Advances in instrumentation have recently extended detailed measurements of gas kinematics to large samples of high-redshift galaxies. Relative to most nearby, thin disc galaxies, in which gas rotation accurately traces the gravitational potential, the interstellar medium (ISM) of $$z$$ ≳ 1 galaxies is typically more dynamic and exhibits elevated turbulence. If not properly modelled, these effects can strongly bias dynamical mass measurements. We use high-resolution FIRE-2 cosmological zoom-in simulations to analyse the physical effects that must be considered to correctly infer dynamical masses from gas kinematics. Our analysis covers a range of galaxy properties from low-redshift Milky-Way-mass galaxies to massive high-redshift galaxies (M⋆ > 1011 M⊙ at $$z$$ = 1). Selecting only snapshots where a disc is present, we calculate the rotational profile $$\bar{v}_\phi (r)$$ of the cool ($$10^{3.5}\,\lt {\it T}\lt 10^{4.5}~\rm {K}$$) gas and compare it to the circular velocity $$v_{\rm c}=\sqrt{GM_{\rm enc}/r}$$. In the simulated galaxies, the gas rotation traces the circular velocity at intermediate radii, but the two quantities diverge significantly in the centre and in the outer disc. Our simulations appear to over-predict observed rotational velocities in the centres of massive galaxies (likely from a lack of black hole feedback), so we focus on larger radii. Gradients in the turbulent pressure at these radii can provide additional radial support and bias dynamical mass measurements low by up to 40 per cent. In both the interior and exterior, the gas’ motion can be significantly non-circular due to e.g. bars, satellites, and inflows/outflows. We discuss the accuracy of commonly used analytic models for pressure gradients (or ‘asymmetric drift’) in the ISM of high-redshift galaxies.more » « less
