Feedback In Realistic Environments Research Articles

Unravelling jet quenching criteria across L* galaxies and massive cluster ellipticals

ABSTRACT In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. Active galactic nucleus jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work concludes from simulations of a halo with $10^{14} \,\mathrm{ M}_\odot$ that a successful jet model should have an energy flux comparable to the free-fall energy flux at the cooling radius and should inflate a sufficiently wide cocoon with a long enough cooling time. In this paper, we investigate three jet modes with constant fluxes satisfying the criteria, including high-temperature thermal jets, cosmic ray (CR)-dominant jets, and widely precessing kinetic jets in $10^{12}-10^{15}\, {\rm M}_{\odot }$ haloes using high-resolution, non-cosmological magnetohydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We find that scaling the jet energy according to the free-fall energy at the cooling radius can successfully suppress the cooling flows and quench galaxies without violating observational constraints. On the contrary, if we scale the energy flux based on the total cooling rate within the cooling radius, strong interstellar medium cooling dominates this scaling, resulting in a jet flux exceeding what is needed. Among the three jet types, the CR-dominant jet is most effective in suppressing cooling flows across all surveyed halo masses due to enhanced CR pressure support. We confirm that the criteria for a successful jet model work across a wider range, encompassing halo masses of $10^{12}-10^{15} {\rm M_\odot }$.

Dense stellar clump formation driven by strong quasar winds in the FIRE cosmological hydrodynamic simulations

ABSTRACT We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star-forming galaxy at z ∼ 2 under the presence of strong quasar winds. Our simulations include multiphase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disc winds. We show that powerful quasar winds can have a global negative impact on galaxy growth while in the strongest cases triggering the formation of an off-centre clump with stellar mass ${\rm M}_{\star }\sim 10^{7}\, {\rm M}_{\odot }$, effective radius ${\rm R}_{\rm 1/2\, \rm Clump}\sim 20\, {\rm pc}$, and surface density $\Sigma _{\star } \sim 10^{4}\, {\rm M}_{\odot }\, {\rm pc}^{-2}$. The clump progenitor gas cloud is originally not star-forming, but strong ram pressure gradients driven by the quasar winds (orders of magnitude stronger than experienced in the absence of winds) lead to rapid compression and subsequent conversion of gas into stars at densities much higher than the average density of star-forming gas. The AGN-triggered star-forming clump reaches ${\rm SFR} \sim 50\, {\rm M}_{\odot }\, {\rm yr}^{-1}$ and $\Sigma _{\rm SFR} \sim 10^{4}\, {\rm M}_{\odot }\, {\rm yr}^{-1}\, {\rm kpc}^{-2}$, converting most of the progenitor gas cloud into stars in ∼2 Myr, significantly faster than its initial free-fall time and with stellar feedback unable to stop star formation. In contrast, the same gas cloud in the absence of quasar winds forms stars over a much longer period of time (∼35 Myr), at lower densities, and losing spatial coherency. The presence of young, ultra-dense, gravitationally bound stellar clumps in recently quenched galaxies could thus indicate local positive feedback acting alongside the strong negative impact of powerful quasar winds, providing a plausible formation scenario for globular clusters.

Synchrotron emission on FIRE: equipartition estimators of magnetic fields in simulated galaxies with spectrally resolved cosmic rays

ABSTRACT Synchrotron emission is one of few observable tracers of galactic magnetic fields (B) and cosmic rays (CRs). Much of our understanding of B in galaxies comes from utilizing synchrotron observations in conjunction with several simplifying assumptions of equipartition models, however, it remains unclear how well these assumptions hold, and what B these estimates physically represent. Using Feedback in Realistic Environments project simulations which self-consistently evolve CR proton, electron, and positron spectra from MeV to TeV energies, we present the first synthetic synchrotron emission predictions from simulated L* galaxies with ‘live’ spectrally resolved CR-magnetohydrodynamic. We find that synchrotron emission can be dominated by relatively cool and dense gas, resulting in equipartition estimates of B with fiducial assumptions underestimating the ‘true’ B in the gas that contributes the most emission by factors of 2–3 due to small volume-filling factors. Motivated by our results, we present an analytical framework that expands upon equipartition models for estimating B in a multiphase medium. Comparing our spectrally resolved synchrotron predictions to simpler spectral assumptions used in galaxy simulations with CRs, we find that spectral evolution can be crucial for accurate synchrotron calculations towards galactic centres, where loss terms are large.

C ii] 158 μm emission as an indicator of galaxy star formation rate

ABSTRACT Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C\, {\small II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C\, {\small II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C\, {\small II}]}{}/{}\rm SFR$ than local SFGs, including the infrared (IR)-luminous, starburst galaxies at low and high redshifts as well as some moderately SFGs at the epoch of re-ionization (EoR). The origins of this ‘$\rm [C\, {\small II}]$ deficit’ is unclear. In this work, we study the $L_{\rm [C\, {\small II}]}$–SFR relation of galaxies using a sample of z = 0–8 galaxies with $M_*\approx 10^7-5\times 10^{11}\, \mathrm{M}_\odot$ extracted from cosmological volume and zoom-in simulations from the Feedback in Realistic Environments (fire) project. We find a simple analytic expression for $L_{\rm [C\, {\small II}]}$/SFR of galaxies in terms of the following parameters: mass fraction of $\rm [C\, {\small II}]$-emitting gas ($f_{\rm [C\, {\small II}]}$), gas metallicity (Zgas), gas density (ngas), and gas depletion time ($t_{\rm dep}{}={}M_{\rm gas}{}/{}\rm SFR$). We find two distinct physical regimes: $\rm H_2$-rich galaxies, where tdep is the main driver of the $\rm [C\, {\small II}]$ deficit and $\rm H_2$-poor galaxies where Zgas is the main driver. The observed $\rm [C\, {\small II}]$ 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 the $\rm [C\, {\small II}]$ deficit is a common phenomenon of galaxies, and caution needs to be taken when applying a constant $L_{\rm [C\, {\small II}]}$-to-SFR conversion factor derived from local SFGs to estimate cosmic SFR density at high redshifts and interpret data from upcoming $\rm [C\, {\small II}]$ line intensity mapping experiments.

Bursty Star Formation Naturally Explains the Abundance of Bright Galaxies at Cosmic Dawn

Recent discoveries of a significant population of bright galaxies at cosmic dawn have enabled critical tests of cosmological galaxy formation models. In particular, the bright end of the galaxys’ UV luminosity functions (UVLFs) appear higher than predicted by many models. Using approximately 25,000 galaxy snapshots at 8 ≤ z ≤ 12 in a suite of FIRE-2 cosmological “zoom-in” simulations from the Feedback in Realistic Environments (FIRE) project, we show that the observed abundance of UV-bright galaxies at cosmic dawn is reproduced in these simulations with a multichannel implementation of standard stellar feedback processes, without any fine-tuning. Notably, we find no need to invoke previously suggested modifications, such as a nonstandard cosmology, a top-heavy stellar initial mass function, or a strongly enhanced star formation efficiency. We contrast the UVLFs predicted by bursty star formation in these original simulations to those derived from star formation histories (SFHs) smoothed over prescribed timescales (e.g., 100 Myr). The comparison demonstrates that the strongly time-variable SFHs predicted by the FIRE simulations play a key role in correctly reproducing the observed, bright-end UVLFs at cosmic dawn: the bursty SFHs induce order-or-magnitude changes in the abundance of UV-bright (M UV ≲ −20) galaxies at z ≳ 10. The predicted bright-end UVLFs are consistent with both the spectroscopically confirmed population and the photometrically selected candidates. We also find good agreement between the predicted and observationally inferred integrated UV luminosity densities, which evolve more weakly with redshift in FIRE than suggested by some other models.

The inefficiency of stellar feedback in driving galactic outflows in massive galaxies at high redshift

ABSTRACT Recent observations indicate that galactic outflows are ubiquitous in high-redshift (high-z) 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 and galactic outflow properties of galaxies in massive haloes ($10^{12}\, {\rm M}_{\odot }\ \lt\ M_{\rm vir}\ \lt\ 5\times 10^{12}\, {\rm M}_{\odot }$) at z ≳ 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. The simulated galaxies resemble z > 4 DSFGs, with star-formation rates of $\sim\!{1000}\ {\rm M}_{\odot }\, \rm yr^{-1}$ and molecular gas masses of Mmol ∼ 1010 M⊙. However, the simulated galaxies are characterized 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. 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 per cent. Finally, compared with FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high-z sample. We argue that both effects originate from the higher total and gas surface densities that characterize high-z massive systems.

Efficient simulations of ionized ISM emission lines: a detailed comparison between the FIRE high-redshift suite and observations

ABSTRACT The Atacama Large Millimeter/Submillimeter Array (ALMA) in the submillimetre and the JWST in the infrared have achieved robust spectroscopic detections of emission lines from the interstellar medium (ISM) in some of the first galaxies. These unprecedented measurements provide valuable information regarding the ISM properties, stellar populations, galaxy morphologies, and kinematics in these high-redshift galaxies and, in principle, offer powerful tests of state of the art galaxy formation models, as implemented in hydrodynamical simulations. To facilitate direct comparisons between simulations and observations, we develop a fast post-processing pipeline to predict line emission from the H ii regions around simulated star particles, accounting for spatial variations in the surrounding gas density, metallicity, and incident radiation spectrum. Our ISM line emission model currently captures H α, H β, and all of the [O iii] and [O ii] lines targeted by ALMA and JWST at z > 6. We illustrate the power of this approach by applying our line emission model to the publicly available Feedback In Realistic Environments (FIRE) high-z simulation suite and perform a detailed comparison with current observations. We show that the FIRE mass–metallicity relation is in 1σ agreement with ALMA/JWST measurements after accounting for the inhomogeneities in the ISM properties. We also quantitatively validate the description of the one-zone model, which is widely used for interpreting [O iii] and H β line luminosity measurements. This model is publicly available and can be implemented on top of a broad range of galaxy formation simulations for comparison with JWST and ALMA measurements.

Local positive feedback in the overall negative: the impact of quasar winds on star formation in the FIRE cosmological simulations

ABSTRACT Negative feedback from accreting supermassive black holes is considered crucial in suppressing star formation and quenching massive galaxies. However, several models and observations suggest that black hole feedback may have a positive effect, triggering star formation by compressing interstellar medium gas to higher densities. We investigate the dual role of black hole feedback using cosmological hydrodynamic simulations from the Feedback In Realistic Environment (FIRE) project, incorporating a novel implementation of hyper-refined accretion-disc winds. Focusing on a massive, star-forming galaxy at z ∼ 2 ($M_{\rm halo} \sim 10^{12.5}\, {\rm M}_{\odot }$), we demonstrate that strong quasar winds with a kinetic power of ∼1046 erg s−1, persisting for over 20 Myr, drive the formation of a central gas cavity and significantly reduce the surface density of star formation across the galaxy’s disc. The suppression of star formation primarily occurs by limiting the availability of gas for star formation rather than by evacuating the pre-existing star-forming gas reservoir (preventive feedback dominates over ejective feedback). Despite the overall negative impact of quasar winds, we identify several potential indicators of local positive feedback, including (1) the spatial anticorrelation between wind-dominated regions and star-forming clumps, (2) higher local star formation efficiency in compressed gas at the edge of the cavity, and (3) increased contribution of outflowing material to local star formation. Moreover, stars formed under the influence of quasar winds tend to be located at larger radial distances. Our findings suggest that both positive and negative AGN feedback can coexist within galaxies, although the local positive triggering of star formation has a minor influence on global galaxy growth.

Dwarf Galaxy Formation with and without Dark Matter–Baryon Streaming Velocities

We study how supersonic streaming velocities of baryons relative to dark matter—a large-scale effect imprinted at recombination and coherent over ∼3 Mpc scales—affect the formation of dwarf galaxies at z ≳ 5. We perform cosmological hydrodynamic simulations, including and excluding streaming velocities, in regions centered on halos with M vir(z = 0) ≈ 1010 M ⊙; the simulations are part of the Feedback In Realistic Environments (FIRE) project and run with FIRE-3 physics. Our simulations comprise many thousands of systems with halo masses between M vir = 2 × 105 M ⊙ and 2 × 109 M ⊙ in the redshift range z = 20–5. A few hundred of these galaxies form stars and have stellar masses ranging from 100 to 107 M ⊙. While star formation is globally delayed by approximately 50 Myr in the streaming relative to nonstreaming simulations and the number of luminous galaxies is correspondingly suppressed at high redshift in the streaming runs, these effects decay with time. By z = 5, the properties of the simulated galaxies are nearly identical in the streaming versus nonstreaming runs, indicating that any effects of streaming velocities on the properties of galaxies at the mass scale of classical dwarfs and larger do not persist to z = 0.

Molecular Gas and Star Formation in Nearby Starburst Galaxy Mergers

We employ the Feedback In Realistic Environments (FIRE-2) physics model to study how the properties of giant molecular clouds (GMCs) evolve during galaxy mergers. We conduct a pixel-by-pixel analysis of molecular gas properties in both the simulated control galaxies and galaxy major mergers. The simulated GMC pixels in the control galaxies follow a similar trend in a diagram of velocity dispersion (σ v ) versus gas surface density (Σmol) to the one observed in local spiral galaxies in the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) survey. For GMC pixels in simulated mergers, we see a significant increase of a factor of 5–10 in both Σmol and σ v , which puts these pixels above the trend of PHANGS galaxies in the σ v versus Σmol diagram. This deviation may indicate that GMCs in the simulated mergers are much less gravitationally bound compared with simulated control galaxies with virial parameters (α vir) reaching 10–100. Furthermore, we find that the increase in α vir happens at the same time as the increase in global star formation rate, which suggests that stellar feedback is responsible for dispersing the gas. We also find that the gas depletion time is significantly lower for high-α vir GMCs during a starburst event. This is in contrast to the simple physical picture that low-α vir GMCs are easier to collapse and form stars on shorter depletion times. This might suggest that some other physical mechanisms besides self-gravity are helping the GMCs in starbursting mergers collapse and form stars.

The impact of AGN-driven winds on physical and observable galaxy sizes

ABSTRACT Without active galactic nucleus (AGN) feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z ≲ 2. In this paper, we perform high-resolution re-simulations of a massive ($M_{\star }\sim 10^{11}\, \rm {{\rm M}_{\odot }}$) galaxy at z ∼ 2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half-mass radius. We experiment with driving mechanical AGN winds, using a state-of-the-art hyper-Lagrangian refinement technique to increase particle resolution. These winds reduce the gas surface density in the inner regions of the galaxy, suppressing the compact starburst and maintaining an approximately constant half-mass radius. Using radiative transfer, we study the impact of AGN feedback on the magnitude and extent of the multiwavelength continuum emission. When AGN winds are included, the suppression of the compact, dusty starburst results in lowered flux at FIR wavelengths (due to decreased star formation) but increased flux at optical-to-near-IR wavelengths (due to decreased dust attenuation, in spite of the lowered star formation rate), relative to the case without AGN winds. The FIR half-light radius decreases from ∼1 to $\sim 0.1\, \rm {kpc}$ in $\lesssim 40\, \rm {Myr}$ when AGN winds are not included, but increases to $\sim 2\, \rm {kpc}$ when they are. Interestingly, the half-light radius at optical-NIR wavelengths remains approximately constant over $35\, \rm {Myr}$, for simulations with and without AGN winds. In the case without winds, this occurs despite the rapid compaction, and is due to heavy dust obscuration in the inner regions of the galaxy. This work highlights the importance of forward-modelling when comparing simulated and observed galaxy populations.

Public Data Release of the FIRE-2 Cosmological Zoom-in Simulations of Galaxy Formation

We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation (available at http://flathub.flatironinstitute.org/fire) from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multiphase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Type Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from three suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way (MW)–mass galaxies, five SMC/LMC-mass galaxies, and four lower-mass galaxies including one ultrafaint; we release 39 snapshots across z = 0–10. The second comprises four massive galaxies, with 19 snapshots across z = 1–10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5–10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all MW-mass galaxies, we release the formation coordinates and an “ex situ” flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations.

Spiral arms are metal freeways: azimuthal gas-phase metallicity variations in flocculent discs in the FIRE-2 cosmological zoom-in simulations

ABSTRACT We examine the azimuthal variations in gas-phase metallicity profiles in simulated Milky Way-mass disc galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite, which includes a sub-grid turbulent metal mixing model. We produce spatially resolved maps of the discs at z ≈ 0 with pixel sizes ranging from 250 to 750 pc, analogous to modern integral field unit galaxy surveys, mapping the gas-phase metallicities in both the cold and dense gas and the ionized gas correlated with H ii regions. We report that the spiral arms alternate in a pattern of metal rich and metal poor relative to the median metallicity of the order of ≲0.1 dex, appearing generally in this sample of flocculent spirals. The pattern persists even in a simulation with different strengths of metal mixing, indicating that the pattern emerges from physics above the sub-grid scale. Local enrichment does not appear to be the dominant source of the azimuthal metallicity variations at z ≈ 0: there is no correlation with local star formation on these spatial scales. Rather, the arms are moving radially inwards and outwards relative to each other, carrying their local metallicity gradients with them radially before mixing into the larger-scale interstellar medium. We propose that the arms act as freeways channeling relatively metal poor gas radially inwards, and relatively enriched gas radially outwards.

Constraining cosmic ray transport with observations of the circumgalactic medium

ABSTRACT Recent theoretical studies predict that the circumgalactic medium (CGM) around low-redshift, ∼L* galaxies could have substantial non-thermal pressure support in the form of cosmic rays. However, these predictions are sensitive to the specific model of cosmic ray transport employed, which is theoretically and observationally underconstrained. In this work, we propose a novel observational constraint for calculating the lower limit of the radially averaged, effective cosmic ray transport rate, ${\kappa _{\rm eff}^{\rm min}}$. Under a wide range of assumptions (so long as cosmic rays do not lose a significant fraction of their energy in the galactic disc, regardless of whether the cosmic ray pressure is important or not in the CGM), we demonstrate a well-defined relationship between ${\kappa _{\rm eff}^{\rm min}}$ and three observable galaxy properties: the total hydrogen column density, the average star formation rate, and the gas circular velocity. We use a suite of Feedback in Realistic Environments 2 galaxy simulations with a variety of cosmic ray transport physics to demonstrate that our analytical model of ${\kappa _{\rm eff}^{\rm min}}$ is a robust lower limit of the true cosmic ray transport rate. We then apply our new model to calculate ${\kappa _{\rm eff}^{\rm min}}$ for galaxies in the COS-Halos sample, and confirm this already reveals strong evidence for an effective transport rate that rises rapidly away from the interstellar medium to values ${\kappa _{\rm eff}^{\rm min}}\gtrsim 10^{30\!-\!31}\, {\rm cm}^2\, {\rm s}^{-1}$ (corresponding to anisotropic streaming velocities of $v^{\rm stream}_{\rm eff} \gtrsim 1000\, {\rm km}\, {\rm s}^{-1}$) in the diffuse CGM, at impact parameters larger than 50–100 kpc. We discuss how future observations can provide qualitatively new constraints in our understanding of cosmic rays in the CGM and intergalactic medium.

Exploring supermassive black hole physics and galaxy quenching across halo mass in FIRE cosmological zoom simulations

ABSTRACT Feedback from accreting supermassive black holes (SMBHs) is thought to be a primary driver of quenching in massive galaxies, but how to best implement SMBH physics into galaxy formation simulations remains ambiguous. As part of the Feedback in Realistic Environments (FIRE) project, we explore the effects of different modelling choices for SMBH accretion and feedback in a suite of ∼500 cosmological zoom-in simulations across a wide range of halo mass (1010–1013 M⊙). Within the suite, we vary the numerical schemes for BH accretion and feedback, accretion efficiency, and the strength of mechanical, radiative, and cosmic ray feedback independently. We then compare the outcomes to observed galaxy scaling relations. We find several models satisfying observational constraints for which the energetics in different feedback channels are physically plausible. Interestingly, cosmic rays accelerated by SMBHs play an important role in many plausible models. However, it is non-trivial to reproduce scaling relations across halo mass, and many model variations produce qualitatively incorrect results regardless of parameter choices. The growth of stellar and BH mass are closely related: for example, overmassive BHs tend to overquench galaxies. BH mass is most strongly affected by the choice of accretion efficiency in high-mass haloes, but by feedback efficiency in low-mass haloes. The amount of star formation suppression by SMBH feedback in low-mass haloes is determined primarily by the time-integrated feedback energy. For massive galaxies, the ‘responsiveness’ of a model (how quickly and powerfully the BH responds to gas available for accretion) is an additional important factor for quenching.

Shapes of Milky-Way-mass galaxies with self-interacting dark matter

ABSTRACT Self-interacting dark matter (SIDM) models offer one way to reconcile inconsistencies between observations and predictions from collisionless cold dark matter (CDM) models on dwarf-galaxy scales. In order to incorporate the effects of both baryonic and SIDM interactions, we study a suite of cosmological-baryonic simulations of Milky-Way (MW)-mass galaxies from the Feedback in Realistic Environments (FIRE-2) project where we vary the SIDM self-interaction cross-section σ/m. We compare the shape of the main dark matter (DM) halo at redshift z = 0 predicted by SIDM simulations (at σ/m = 0.1, 1, and 10 cm2 g−1) with CDM simulations using the same initial conditions. In the presence of baryonic feedback effects, we find that SIDM models do not produce the large differences in the inner structure of MW-mass galaxies predicted by SIDM-only models. However, we do find that the radius where the shape of the total mass distribution begins to differ from that of the stellar mass distribution is dependent on σ/m. This transition could potentially be used to set limits on the SIDM cross-section in the MW.

Extinguishing the FIRE: environmental quenching of satellite galaxies around Milky Way-mass hosts in simulations

ABSTRACT The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in z = 0 galaxies at $M_{*}=10^{5\!-\!10}\, \rm {M}_{\odot }$ around 14 MW-mass hosts from the Feedback In Realistic Environments 2 (FIRE-2) simulations. Lower mass satellites ($M_{*}\lesssim 10^7\, \rm {M}_{\odot }$) are mostly quiescent and higher mass satellites ($M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$) are mostly star forming, with intermediate-mass satellites ($M_{*}\approx 10^{7\!-\!8}\, \rm {M}_{\odot }$) split roughly equally between quiescent and star forming. Hosts with more gas in their circumgalactic medium have a higher quiescent fraction of massive satellites ($M_{*}=10^{8\!-\!9}\, \rm {M}_{\odot }$). We find no significant dependence on isolated versus paired (LG-like) host environments, and the quiescent fractions of satellites around MW-mass and Large Magellanic Cloud (LMC)-mass hosts from the FIRE-2 simulations are remarkably similar. Environmental effects that lead to quenching can also occur as pre-processing in low-mass groups prior to MW infall. Lower mass satellites typically quenched before MW infall as central galaxies or rapidly during infall into a low-mass group or a MW-mass galaxy. Most intermediate- to high-mass quiescent satellites have experienced ≥1–2 pericentre passages (≈2.5–5 Gyr) within a MW-mass halo. Most galaxies with $M_{*}\gtrsim 10^{6.5}\, \rm {M}_{\odot }$ did not quench before falling into a host, indicating a possible upper mass limit for isolated quenching. The simulations reproduce the average trend in the LG quiescent fraction across the full range of satellite stellar masses. Though the simulations are consistent with the Satellites Around Galactic Analogs (SAGA) survey’s quiescent fraction at $M_{*}\gtrsim 10^8\, \rm {M}_{\odot }$, they do not generally reproduce SAGA’s turnover at lower masses.

The In Situ Origins of Dwarf Stellar Outskirts in FIRE-2

Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies (108.5 < M ⋆/M ⊙ < 109.6) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structures, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in situ, largely via the outward migration of disk stars. However, there also exists a large population of “nondisky” dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These nondisky dwarfs tend to be either more gas-poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy’s intrinsic shape—which is a direct quantification of the distribution of the galaxy’s stellar content—to interrogate the feedback implementation in simulated galaxies.

Relating the Diverse Merger Histories and Satellite Populations of Nearby Galaxies

We investigate whether the considerable diversity in the satellite populations of nearby Milky Way (MW)-mass galaxies is connected with the diversity in their host’s merger histories. Analyzing eight nearby galaxies with extensive observations of their satellite populations and stellar halos, we characterize each galaxy’s merger history using the metric of its most dominant merger, M ⋆,Dom, defined as the greater of either its total accreted stellar mass or most massive current satellite. We find an unexpectedly tight relationship between these galaxies’ number of M V < − 9 satellites within 150 kpc (N Sat) and M ⋆,Dom. This relationship remains even after accounting for differences in galaxy mass. Using the star formation and orbital histories of satellites around the MW and M81, we demonstrate that both likely evolved along the M ⋆,Dom–N Sat relation during their current dominant mergers with the Large Magellanic Cloud and M82, respectively. We investigate the presence of this relation in galaxy formation models, including using the Feedback In Realistic Environments (FIRE) simulations to directly compare to the observations. We find no relation between M ⋆,Dom and N Sat in FIRE, and a universally large scatter in N Sat with M ⋆,Dom across simulations—in direct contrast with the tightness of the empirical relation. This acute difference in the observed and predicted scaling relation between two fundamental galaxy properties signals that current simulations do not sufficiently reproduce diverse merger histories and their effects on satellite populations. Explaining the emergence of this relation is therefore essential for obtaining a complete understanding of galaxy formation.

The effects of LMC-mass environments on their dwarf satellite galaxies in the FIRE simulations

ABSTRACTCharacterizing the predicted environments of dwarf galaxies like the Large Magellanic Cloud (LMC) is becoming increasingly important as next-generation surveys push sensitivity limits into this low-mass regime at cosmological distances. We study the environmental effects of LMC-mass haloes (M200m ∼ 1011 M⊙) on their populations of satellites (M⋆ ≥ 104 M⊙) using a suite of zoom-in simulations from the Feedback In Realistic Environments (FIRE) project. Our simulations predict significant hot coronas with T ∼ 105 K and Mgas ∼ 109.5 M⊙. We identify signatures of environmental quenching in dwarf satellite galaxies, particularly for satellites with intermediate mass (M⋆ = 106–107 M⊙). The gas content of such objects indicates ram pressure as the likely quenching mechanism, sometimes aided by star formation feedback. Satellites of LMC-mass hosts replicate the stellar mass dependence of the quiescent fraction found in satellites of Milky Way-mass hosts (i.e. that the quiescent fraction increases as stellar mass decreases). Satellites of LMC-mass hosts have a wider variety of quenching times when compared to the strongly bimodal distribution of quenching times of nearby centrals. Finally, we identify significant tidal stellar structures around four of our six LMC analogues, suggesting that stellar streams may be common. These tidal features originated from satellites on close orbits, extend to ∼80 kpc from the central galaxy, and contain ∼106–107 M⊙ of stars.