Search

We present an analysis of the R ≲ 1.5 kpc core regions of seven simulated Milky Way-mass galaxies, from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, for a finely sampled period (Δt = 2.2 Myr) of 22 Myr at z ≈ 0, and compare them with star formation rate (SFR) and gas surface density observations of the Milky Way's Central Molecular Zone (CMZ). Despite not being tuned to reproduce the detailed structure of the CMZ, we find that four of these galaxies are consistent with CMZ observations at some point during this 22 Myr period. The galaxies presented here are not homogeneous in their central structures, roughly dividing into two morphological classes; (a) several of the galaxies have very asymmetric gas and SFR distributions, with intense (compact) starbursts occurring over a period of roughly 10 Myr, and structures on highly eccentric orbits through the CMZ, whereas (b) others have smoother gas and SFR distributions, with only slowly varying SFRs over the period analyzed. In class (a) centers, the orbital motion of gas and star-forming complexes across small apertures (R ≲ 150 pc, analogously |l| < 1° in the CMZ observations) contributes as much to tracers of star formation/dense gas appearing in those apertures, as the internal evolution of those structures does. These asymmetric/bursty galactic centers can simultaneously match CMZ gas and SFR observations, demonstrating that time-varying star formation can explain the CMZ's low star formation efficiency.

We present models of CO(1-0) emission from Milky-Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which, in our models, sets the attenuation of the incident UV radiation field. At the resolution of these simulations, commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H2 abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of ~3 pc in cold gas (T < 300 K) for both CO and H2, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H2 conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model, which controls the amount of radiation that penetrates giant molecular clouds.

We present the first measurement of the lifetimes of giant molecular clouds (GMCs) in cosmological simulations at z = 0, using the Latte suite of FIRE-2 simulations of Milky Way (MW) mass galaxies. We track GMCs with total gas mass ≳105 M⊙ at high spatial (∼1 pc), mass (7100 M⊙), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the MW and nearby galaxies. We find GMC lifetimes of 5-7 Myr, or 1-2 freefall times, on average, with less than 2 per cent of clouds living longer than 20 Myr. We find decreasing GMC lifetimes with increasing virial parameter, and weakly increasing GMC lifetimes with galactocentric radius, implying that environment affects the evolutionary cycle of GMCs. However, our GMC lifetimes show no systematic dependence on GMC mass or amount of star formation. These results are broadly consistent with inferences from the literature and provide an initial investigation into ultimately understanding the physical processes that govern GMC lifetimes in a cosmological setting.

Giant molecular clouds (GMCs) are well studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of GMCs in a Milky Way-like galaxy and an Large Magellanic Cloud (LMC)-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic centre occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures.With the exception of metallicity, we find no significant differences in cloud properties between the MilkyWay-like anddwarf galaxies.These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties.

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.

While many tensions between Local Group (LG) satellite galaxies and ⋀ cold dark matter cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with M∗ > 105 M☉ around eight isolated Milky Way (MW) mass host galaxies and four hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with M∗ ≲ 105 M☉. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances ≾100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic discs of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhaloes in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at M∗ ≳ 105 M☉: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogues in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.

Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of ∼3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.