Author: Claude-Andre Faucher-Giguere

Hopkins et al., available on arXiv

Abstract: As they grow, galaxies can transition from irregular/spheroidal with ‘bursty’ star formation histories (SFHs), to disky with smooth SFHs. But even in simulations, the direct physical cause of such transitions remains unclear. We therefore explore this in a large suite of numerical experiments re-running portions of cosmological simulations with widely varied physics, further validated with existing FIRE simulations. We show that gas supply, cooling/thermodynamics, star formation model, Toomre scale, galaxy dynamical times, and feedback properties do not have a direct causal effect on these transitions. Rather, both the formation of disks and cessation of bursty star formation are driven by the gravitational potential, but in different ways. Disk formation is promoted when the mass profile becomes sufficiently centrally-concentrated in shape (relative to circularization radii): we show that this provides a well-defined dynamical center, ceases to support the global ‘breathing modes’ which can persist indefinitely in less-concentrated profiles and efficiently destroy disks, promotes orbit mixing to form a coherent angular momentum, and stabilizes the disk. Smooth SF is promoted by the potential or escape velocity (not circular velocity) becoming sufficiently large at the radii of star formation that cool, mass-loaded (momentum-conserving) outflows are trapped/confined near the galaxy, as opposed to escaping after bursts. We discuss the detailed physics, how these conditions arise in cosmological contexts, their relation to other correlated phenomena (e.g. inner halo virialization, vertical disk ‘settling’), and observations.

Parul et al., available on arXiv

Abstract: Milky Way-mass galaxies in the FIRE-2 simulations demonstrate two main modes of star formation. At high redshifts star formation occurs in a series of short and intense bursts, while at low redshifts star formation proceeds at a steady rate with a transition from one mode to another at times ranging from 3 to 7 Gyr ago for different galaxies. We analyse how the mode of star formation affects iron and alpha-element abundance. We find that the early bursty regime imprints a measurable pattern in stellar elemental abundances in the form of a “sideways chevron” shape on the [Fe/H] – [O/Fe] plane and the scatter in [O/Fe] at a given stellar age is higher than when a galaxy is in the steady regime. That suggests that the evolution of [O/Fe] scatter with age provides an estimate of the end of the bursty phase. We investigate the feasibility of observing of this effect by adding mock observational errors to a simulated stellar survey and find that the transition between the bursty and steady phase should be detectable in the Milky Way, although larger observational uncertainties make the transition shallower. We apply our method to observations of the Milky Way from the Second APOKASC Catalog and estimate that the transition to steady star formation in the Milky Way happened 7-8 Gyrs ago, earlier than transition times measured in the simulations.

McCluskey et al., available on arXiv

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 disk settling and post-formation dynamical heating. We identify three eras: a Pre-Disk Era (typically >8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disk Era (~ 8 – 4 Gyr ago), when stars started to form on rotation-dominated orbits but with high velocity dispersion, sigma_form; and a Late-Disk Era (< 4 Gyr ago), when stars formed with low sigma_form. sigma_form increased with time during the Pre-Disk Era, peaking ~ 8 Gyr ago, then decreased throughout the Early-Disk Era as the disk settled and remained low throughout the Late-Disk Era. By contrast, the velocity dispersion measured today, sigma_now, increases monotonically with age because of stronger post-formation heating for Pre-Disk stars. Importantly, most of sigma_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, sigma_R,form > sigma_phi,form > sigma_Z,form. Post-formation heating primarily increased sigma_R at ages < 4 Gyr but acted nearly isotropically for older stars. The lookback time that the disk began to settle correlates with its dynamical state today: earlier-settling galaxies currently form colder disks. Young stars in FIRE-2 are kinematically hotter than the MW but broadly agree with M31 and M33. Including stellar cosmic-ray feedback does not significantly change the amount of disk rotational support at fixed stellar mass.

Liang et al., available on arXiv

Abstract: Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity (L[CII]) and star formation rate (SFR), suggesting that L[CII] may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower L[CII]/SFR than the local SFGs, including the infrared-luminous, starburst galaxies at low and high redshifts, as well as some moderately star-forming galaxies at the epoch of re-ionization (EoR). The origin of this ‘CII] deficit’ is unclear. In this work, we study the L[CII]-SFR relation of galaxies using a sample of z=0−8 galaxies with Mstar~10^7−5×10^11 Msun extracted from cosmological volume and zoom-in simulations from the Feedback in Realistic Environments (FIRE) project. We find a simple analytic expression for L[CII]/SFR of galaxies in terms of the following parameters: mass fraction of [CII]-emitting gas (f[CII]), gas metallicity (Zgas), gas density (ngas) and gas depletion time (tdep=Mgas/SFR). We find two distinct physical regimes, where tdep (Zgas) is the main driver of the [CII] deficit in H2-rich (H2-poor) galaxies. The observed [CII] deficit of IR-luminous galaxies and early EoR galaxies, corresponding to the two different regimes, is due to short gas depletion time and low gas metallicity, respectively. Our result indicates that [CII] deficit is a common phenomenon of galaxies, and caution needs to be taken when applying a constant L[CII]-to-SFR conversion factor derived from local SFGs to estimate cosmic SFR density at high redshifts and interpret data from upcoming [CII] line intensity mapping experiments.

Barry et al., available on arXiv

Abstract: A variety of observational campaigns seek to test dark-matter models by measuring dark-matter subhaloes at low masses. Despite their predicted lack of stars, these subhaloes may be detectable through gravitational lensing or via their gravitational perturbations on stellar streams. To set measurable expectations for subhalo populations within LambdaCDM, we examine 11 Milky Way (MW)-mass haloes from the FIRE-2 baryonic simulations, quantifying the counts and orbital fluxes for subhaloes with properties relevant to stellar stream interactions: masses down to 10^6 Msun, distances < 50 kpc of the galactic center, across z = 0 - 1 (lookback time 0 - 8 Gyr). We provide fits to our results and their dependence on subhalo mass, distance, and lookback time, for use in (semi)analytic models. A typical MW-mass halo contains ~16 subhaloes >10^7 Msun (~1 subhalo >10^8 Msun) within 50 kpc at z = 0. We compare our results with dark-matter-only versions of the same simulations: because they lack a central galaxy potential, they overpredict subhalo counts by 2-10x, more so at smaller distances. Subhalo counts around a given MW-mass galaxy declined over time, being ~10x higher at z = 1 than at z = 0. Subhaloes have nearly isotropic orbital velocity distributions at z = 0. Across our simulations, we also identified 4 analogs of Large Magellanic Cloud satellite passages; these analogs enhance subhalo counts by 1.4-2.7 times, significantly increasing the expected subhalo population around the MW today. Our results imply an interaction rate of ~5 per Gyr for a stream like GD-1, sufficient to make subhalo-stream interactions a promising method of measuring dark subhaloes.

Nguyen et al., available on arXiv

Abstract: The third data release (DR3) of Gaia has provided a five-fold increase in the number of radial velocity measurements of stars, as well as a stark improvement in parallax and proper motion measurements. To help with studies that seek to test models and interpret Gaia DR3, we present nine Gaia synthetic surveys, based on three solar positions in three Milky-Way-mass galaxies of the Latte suite of the FIRE-2 cosmological simulations. These synthetic surveys match the selection function, radial velocity measurements, and photometry of Gaia DR3, adapting the code base Ananke, previously used to match the Gaia DR2 release in Sanderson et al. 2020. The synthetic surveys are publicly available and can be found at this http URL. Similarly to the previous release of Ananke, these surveys are based on cosmological simulations and thus able to model non-equilibrium dynamical effects, making them a useful tool in testing and interpreting Gaia DR3.

Samuel et al., available on arXiv

Abstract: Low-mass galaxies are highly susceptible to environmental effects that can efficiently quench star formation. We explore the role of ram pressure in quenching low-mass galaxies (M∗∼10^(5-9) Msun) within 2 Mpc of Milky Way (MW) hosts using the FIRE-2 simulations. Ram pressure is highly variable across different environments, within individual MW haloes, and for individual low-mass galaxies over time. The impulsiveness of ram pressure — the maximum ram pressure scaled to the integrated ram pressure prior to quenching — correlates with whether a galaxy is quiescent or star-forming. The time-scale between maximum ram pressure and quenching is anticorrelated with impulsiveness, such that high impulsiveness corresponds to quenching time-scales <1 Gyr. Galaxies in low-mass groups (M∗,host∼10^(7-9) Msun) outside of MW haloes experience typical ram pressure only slightly lower than ram pressure on MW satellites, helping to explain effective quenching via group pre-processing. Ram pressure on MW satellites rises sharply with decreasing distance to the host, and, at a fixed physical distance, more recent pericentre passages are typically associated with higher ram pressure because of greater gas density in the inner host halo at late times. Furthermore, the inner gas density of Local Group-like paired host haloes is larger at small angles (<~30 deg) off the host galaxy disc, compared to isolated hosts. The ram pressure and quiescent fraction of satellites within these low latitude regions are correspondingly elevated, signaling anisotropic quenching via ram pressure around paired hosts.

Baptista et al., available on arXiv

Abstract: The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way’s DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FIRE-2 Milky Way-mass galaxies to explore the evolution of the DM halo’s orientation with radius and time, with or without a major merger with a Large Magellanic Cloud (LMC) analog, and when varying the DM model. In three of the four CDM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20 degrees beyond about 30 galactocentric kpc, reaching a maximum of 30–90 degrees depending on the individual halo’s formation history. In identical simulations using a model of self-interacting DM with sigma=1 cm^2 g^−1, the halo remains aligned with the stellar disk out to ~200–400 kpc. Interactions with massive satellites (M>~4×10^10 Msun at pericenter; M>~3.3×10^10 Msun at infall) affect the orientation of the halo significantly, aligning the halo’s major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of SIDM if the effects of massive satellites can be accounted for.

Pandya et al., available on arXiv

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 co-evolve 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 is remarkably successful in reproducing the baryon cycles of the simulated galaxies and their CGM. In particular, we find that, for these parameters, the phase transition occurs at high-redshift in ultrafaint progenitors and at low redshift in classical Mvir∼10^11 Msun dwarfs, while Milky Way-mass halos undergo the transition at z~0.5, in agreement with the simulations. We discuss the ways in which our model is complementary to existing approaches for modeling the CGM–galaxy connection and possible future directions.

Bassini et al., available on arXiv

Abstract: Recent observations indicate that galactic outflows are ubiquitous in high redshift galaxies, including normal star forming galaxies, quasar hosts, and dusty star forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star formation histories (SFH) and galactic outflow properties of galaxies in massive haloes (10^12 Msun~5.5 in three zoom-in cosmological simulations from the MassiveFIRE suite, as part of the Feedback In Realistic Environments (FIRE) project. The simulations were run with the FIRE-2 model, which does not include feedback from active galactic nuclei (AGN). The simulated galaxies resemble z>4 DSFGs, with SFRs of ~1000 Msun yr^−1 and molecular gas masses of Mmol~10^10 Msun. However, the simulated galaxies are characterised by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of ~0.1, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium (ISM). We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of 30%. Finally, compared to FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high redshift sample. We argue that both effects originate from the higher total and gas surface densities that characterise high−z massive systems.