Category: Uncategorized

Keating et al., available on arXiv

Abstract: 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. 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.

Santistevan et al., available on arXiv

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 a 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 FIRE-2 cosmological simulations, including 6 isolated MW/M31-mass galaxies and 6 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 a MW/M31-mass galaxy emerged typically at z ~ 3-4 (11.6-12.2 Gyr ago), while stars in the bulge region (inner 2 kpc) at z = 0 formed primarily in a single main progenitor at z < 5 (< 12.6 Gyr ago). Compared with isolated hosts, the main progenitors of LG-like paired hosts emerged significantly earlier (Delta z ~ 2, Delta t ~ 1.6 Gyr), with ~ 4x higher stellar mass at all z > 4 (> 12.2 Gyr ago). This highlights the importance of environment in MW/M31-mass galaxy formation, especially at early times. Overall, about 100 galaxies with M_star > 10^5 M_sun formed a typical MW/M31-mass system. Thus, surviving satellites represent a highly incomplete census (by ~ 5x) of the progenitor population.

Yu et al., available on arXiv

Abstract: We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that 5-40% of the outer (50-300 kpc) stellar halo in each system consists of in-situ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at z=0.In the outer halo, outflow stars are more spatially homogeneous, metal rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to ~10% of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as ~40% for the most metal-rich local halo stars ([Fe/H] > -0.5). We conclude that the Milky Way stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell & Sandange, who suggested that halo stars formed in rapidly infalling gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme velocity stars in the Milky Way.

Benincasa et al., available on arXiv

Abstract: 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-mass galaxies. We track GMCs with total gas mass >~10^5 Msun at high spatial (~1 pc), mass (7100 Msun), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the Milky Way and nearby galaxies. We find GMC lifetimes of 5-7 Myr, or 1-2 freefall times, on average, with less than 1% of clouds living longer than 20 Myr. We find 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.

Orr et al., available on arXiv

Abstract: We study the spatially resolved (sub-kpc) gas velocity dispersion (sigma)–star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way mass disk galaxies at late times. In agreement with observations, we find a relatively flat relationship, with sigma ~15-30 km/s in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant sigma. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher sigma at constant SFR. The limits of the sigma-Sigma_SFR relation correspond to the onset of strong outflows. We see evidence of “on-off” cycles of star formation in the simulations, corresponding to feedback injection timescales of 10-100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disk (Toomre’s Q=1) model predictions, and the data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e., 1% per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of sigma, it appears to be strongly subdominant to stellar feedback in the simulated galaxy disks.

Hani et al., available on arXiv

Abstract: The correlation between galaxies’ integrated stellar masses and star formation rates (the `star formation main sequence’; SFMS) is a well-established scaling relation. Recently, surveys have found a relationship between the star formation rate and stellar mass surface densities on kpc and sub-kpc scales (the `resolved SFMS’; rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in FIRE-2 zoom-in simulations of Milky Way-mass galaxies. We make SFR and stellar mass maps of the simulated galaxies at a variety of spatial resolutions and star formation averaging time-scales and fit the rSFMS using multiple methods from the literature. While the absolute value of the SFMS slope depends on the fitting method, the slope is steeper for longer star formation time-scales and lower spatial resolutions regardless of the fitting method employed. We present a toy model that quantitatively captures the dependence of the simulated galaxies’ rSFMS slope on spatial resolution and use it to illustrate how this dependence can be used to constrain the characteristic mass of star-forming clumps.

Hafen et al., available on arXiv

Abstract: We analyze the different fates of the circumgalactic medium (CGM) in FIRE-2 cosmological simulations, focusing on the redshifts z=0.25 and z=2 representative of recent surveys. Our analysis includes 21 zoom-in simulations covering the halo mass range Mh(z=0) ~ 10^10 – 10^12 Msun. We analyze both where the gas ends up after first leaving the CGM (its “proximate” fate), as well as its location at z=0 (its “ultimate” fate). Of the CGM at z=2, about half is found in the ISM or stars of the central galaxy by z=0 in Mh(z=2) ~ 5e11 Msun halos, but most of the CGM in lower-mass halos ends up in the IGM. This is so even though most of the CGM in M_h(z=2) ~ 5e10 Msun halos first accretes onto the central galaxy before being ejected into the IGM. On the other hand, most of the CGM mass at z=0.25 remains in the CGM by z=0 at all halo masses analyzed. Of the CGM gas that subsequently accretes onto the central galaxy in the progenitors of Mh(z=0) ~10^12 Msun halos, most of it is cool (T~10^4 K) at z=2 but hot (~Tvir) at z=0.25, consistent with the expected transition from cold mode to hot mode accretion. Despite the transition in accretion mode, at both z=0.25 and z=2 >~80% of the cool gas in Mh >~ 10^11 Msun halos will accrete onto a galaxy. We find that the metallicity of CGM gas is typically a poor predictor of both its proximate and ultimate fates. This is because there is in general little correlation between the origin of CGM gas and its fate owing to substantial mixing while in the CGM.

Wellons et al., available on arXiv

Abstract: Advances in instrumentation have recently extended detailed measurements of gas kinematics to large samples of high-redshift galaxies. Relative to most nearby, thin disk 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 modeled, these effects can strongly bias dynamical mass measurements. We use high-resolution FIRE-2 cosmological zoom-in simulations to analyze the physical effects that must be considered to correctly infer dynamical masses from gas kinematics. Our analysis covers a wide range of galaxy properties, from low-redshift Milky-Way-mass galaxies to massive high-redshift galaxies (M_* > 10^11 M_sun at z=1). Selecting only snapshots where a well-ordered disk is present, we calculate the rotational profile (r) of the cool (10^3.5 K < T < 10^4.5 K) gas and compare it to the circular velocity v_c=sqrt(GM_enc/r) assuming spherical symmetry. In the simulated massive high-redshift galaxies, the gas rotation traces the circular velocity reasonably well at intermediate radii r~1-3 kpc, but the two quantities diverge significantly outside that range. At larger radii, gradients in the turbulent pressure can bias dynamical mass measurements low by ~10-40%. In the interior, the assumption of a spherically-symmetric gravitational potential becomes increasingly poor owing to a massive disk component, reducing the gas rotational velocities by >~10%. Finally, in 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.

Ji et al., available on arXiv

Abstract: We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity kappa~3×10^29 cm^2 s^-1 chosen to match gamma-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at z<~1-2. The gas in these ``CR-dominated'' halos differs significantly from runs without CRs: the gas is primarily cool (a few ~10^4 K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the ``low'' and ``mid'' ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns >~10^14.5 cm^-2 at distances >~150 kpc. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot (T>~10^5 K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.

Ma et al., available on arXiv

Abstract: We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM~M^-2. The cluster formation efficiency is similar across galaxies with stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.