Search

The observationally inferred size versus stellar–mass relationship (SMR) for low-mass galaxies provides an important test for galaxy formation models. However, the relationship relies on assumptions that relate observed luminosity profiles to underlying stellar mass profiles. Here we use the Feedback in Realistic Environments simulations of low-mass galaxies to explore how the predicted SMR changes depending on whether one uses star-particle counts directly or mock observations. We reproduce the SMR found in The Exploration of Local Volume Satellites survey remarkably well only when we infer stellar masses and sizes using mock observations. However, when we use star particles to directly infer stellar masses and half-mass radii, we find that our galaxies are too large and obey an SMR with too little scatter compared to observations. This discrepancy between the 'true' galaxy size and mass and those derived in the mock observation approach is twofold. First, our simulated galaxies have higher and more varied mass-to-light ratios (MLR) at a fixed colour than those commonly adopted, which tends to underestimate their stellar masses compared to their true, simulated values. Second, our galaxies have radially increasing MLR gradients therefore using a single MLR tends to underpredict the mass in the outer regions. Similarly, the true half-mass radius is larger than the half-light radius because the light is more concentrated than the mass. If our simulations are accurate representations of the real Universe, then the relationship between galaxy size and stellar mass is even tighter for low-mass galaxies than is commonly inferred from observed relations. [ABSTRACT FROM AUTHOR]

Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the H  i properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP- Pathfinder and Fire box). This set of simulations comprises galaxies evolved self-consistently across cosmic time with different baryonic sub-grid physics: three different star formation models [constant star formation efficiency (SFE) with different star formation eligibility criteria, and an environmentally dependent, turbulence-based SFE] and two different feedback prescriptions, where only one sub-sample includes early stellar feedback. We use these simulations to assess the impact of different baryonic physics on the H  i content of galaxies. We find that the galaxy-wide H  i properties agree with each other and with observations. However, differences appear for small-scale properties. The thin H  i discs observed in the local universe are only reproduced with a turbulence-dependent SFE and/or early stellar feedback. Furthermore, we find that the morphology of H  i discs is particularly sensitive to the different physics models: galaxies simulated with a turbulence-based SFE have discs that are smoother and more rotationally symmetric, compared to those simulated with a constant SFE; galaxies simulated with early stellar feedback have more regular discs than supernova-feedback-only galaxies. We find that the rotational asymmetry of the H  i discs depends most strongly on the underlying physics model, making this a promising observable for understanding the physics responsible for shaping the interstellar medium of galaxies. [ABSTRACT FROM AUTHOR]

Using cosmological hydrodynamical zoom-in simulations, we explore the properties of subhalos in Milky Way analogs that contain a subcomponent of atomic dark matter (ADM). ADM differs from cold dark matter (CDM) due to the presence of self-interactions that lead to energy dissipation, analogous to standard model baryons. This model can arise in dark sectors that are natural and theoretically motivated extensions to the standard model. The simulations used in this work were carried out using GIZMO and utilize the FIRE-2 galaxy formation physics in the standard model baryonic sector. For the parameter points we consider, the ADM gas cools efficiently, allowing it to collapse to the center of subhalos. This increases a subhalo's central density and affects its orbit, with more subhalos surviving small pericentric passages. The subset of subhalos that host satellite galaxies have cuspier density profiles and smaller stellar half-mass radii relative to CDM. The entire population of dwarf galaxies produced in the ADM simulations is more compact than those seen in CDM simulations, unable to reproduce the entire diversity of observed dwarf galaxy structures. Additionally, we also identify a population of highly compact subhalos that consist nearly entirely of ADM and form in the central region of the host, where they can leave distinctive imprints in the baryonic disk. This work presents the first detailed exploration of subhalo properties in a strongly dissipative dark matter scenario, providing intuition for how other regions of ADM parameter space, as well as other dark sector models, would impact galactic-scale observables. [ABSTRACT FROM AUTHOR]

Cosmic rays (CRs) may drive outflows and alter the phase structure of the circumgalactic medium, with potentially important implications on galaxy formation. However, these effects ultimately depend on the dominant mode of transport of CRs within and around galaxies, which remains highly uncertain. To explore potential observable constraints on CR transport, we investigate a set of cosmological fire -2 CR-magnetohydrodynamic simulations of L* galaxies which evolve CRs with transport models motivated by self-confinement (SC) and extrinsic turbulence (ET) paradigms. To first order, the synchrotron properties diverge between SC and ET models due to a CR physics-driven hysteresis. SC models show a higher tendency to undergo 'ejective' feedback events due to a runaway buildup of CR pressure in dense gas due to the behaviour of SC transport scalings at extremal CR energy densities. The corresponding CR wind-driven hysteresis results in brighter, smoother, and more extended synchrotron emission in SC runs relative to ET and constant diffusion runs. The differences in synchrotron arise from different morphology, interstellar medium gas, and B properties, potentially ruling out SC as the dominant mode of CR transport in typical star-forming L* galaxies, and indicating the prospect for non-thermal radio continuum observations to constrain CR transport physics. [ABSTRACT FROM AUTHOR]

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 (Z gas), gas density (n gas), and gas depletion time (⁠|$t_{\rm dep}{}={}M_{\rm gas}{}/{}\rm SFR$|⁠). We find two distinct physical regimes: |$\rm H_2$| -rich galaxies, where t dep is the main driver of the |$\rm [C\, {\small II}]$| deficit and |$\rm H_2$| -poor galaxies where Z gas 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. [ABSTRACT FROM AUTHOR]

We use high-resolution, hydrodynamic, galaxy simulations from the Latte suite of FIRE-2 simulations to investigate the inherent variation of dark matter in sub-sampled regions around the Solar Circle of a Milky Way-type analogue galaxy and its impact on direct dark matter detection. These simulations show that the baryonic back reaction, as well as the assembly history of substructures, has lasting impacts on the dark matter's spatial and velocity distributions. These are experienced as 'gusts' of dark matter wind around the Solar Circle, potentially complicating interpretations of direct detection experiments on Earth. We find that the velocity distribution function in the galactocentric frame shows strong deviations from the Maxwell Boltzmann form typically assumed in the fiducial Standard Halo Model, indicating the presence of high-velocity substructures. By introducing a new numerical integration technique that removes any dependencies on the Standard Halo Model, we generate event-rate predictions for both single-element Germanium and compound Sodium Iodide detectors, and explore how the variability of dark matter around the Solar Circle influences annual modulation signal predictions. We find that these velocity substructures contribute additional astrophysical uncertainty to the interpretation of event rates, although their impact on summary statistics, such as the peak day of annual modulation, is generally low. [ABSTRACT FROM AUTHOR]

We introduce a suite of cosmological volume simulations to study the evolution of galaxies as part of the Feedback in Realistic Environments project. FIREbox, the principal simulation of the present suite, provides a representative sample of galaxies (∼1000 galaxies with |$M_{\rm star}\gt 10^8\, M_\odot$| at z  = 0) at a resolution (⁠|$\Delta {}x\sim {}20\, {\rm pc}$| , |$m_{\rm b}\sim {}6\times {}10^4\, M_\odot$|) comparable to state-of-the-art galaxy zoom-in simulations. FIREbox captures the multiphase nature of the interstellar medium in a fully cosmological setting (L  = 22.1 Mpc) thanks to its exceptionally high dynamic range (≳106) and the inclusion of multichannel stellar feedback. Here, we focus on validating the simulation predictions by comparing to observational data. We find that star formation rates, gas masses, and metallicities of simulated galaxies with |$M_{\rm star}\lt 10^{10.5-11}\, M_\odot$| broadly agree with observations. These galaxy scaling relations extend to low masses (⁠|$M_{\rm star}\sim {}10^7\, M_\odot$|) and follow a (broken) power-law relationship. Also reproduced are the evolution of the cosmic HI density and the HI column density distribution at z ∼ 0–5. At low z , FIREbox predicts a peak in the stellar-mass–halo-mass relation but also a higher abundance of massive galaxies and a higher cosmic star formation rate density than observed, showing that stellar feedback alone is insufficient to reproduce the properties of massive galaxies at late times. Given its high resolution and sample size, FIREbox offers a baseline prediction of galaxy formation theory in a ΛCDM Universe while also highlighting modelling challenges to be addressed in next-generation galaxy simulations. [ABSTRACT FROM AUTHOR]

Many recent numerical studies have argued that cosmic rays (CRs) from supernovae (SNe) or active galactic nuclei (AGNs) could play a crucial role in galaxy formation, in particular by establishing a CR-pressure-dominated circumgalactic medium (CGM). But explicit CR-magnetohydrodynamics (CR-MHD) remains computationally expensive, and it is not clear whether those results can be applied to simulations that do not explicitly treat magnetic fields or resolved interstellar medium phase structure. We therefore present an intentionally extremely simplified 'sub-grid' model for CRs, which attempts to capture the key qualitative behaviors of greatest interest for those interested in simulations or semi-analytical models including some approximate CR effects on galactic (≳ kpc) scales, while imposing negligible computational overhead. The model is numerically akin to some recently developed sub-grid models for radiative feedback, and allows for a simple constant parametrization of the CR diffusivity and/or streaming speed; it allows for an arbitrary distribution of sources (proportional to black hole accretion rates or star–particle SNe rates or gas/galaxy star formation rates), and interpolates between the limits where CRs escape the galaxies with negligible losses and those where CRs lose most of their energy catastrophically before escape (relevant in e.g. starburst galaxies). The numerical equations are solved trivially alongside gravity in most codes. We compare this to explicit CR-MHD simulations and discuss where the (many) sub-grid approximations break down, and what drives the major sources of uncertainty. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Several recent simulations of galaxy formation predict two main phases of supermassive black hole (BH) accretion: an early, highly intermittent phase (during which BHs are undermassive relative to local scaling relations), followed by a phase of accelerated growth. We investigate physical factors that drive the transition in BH accretion in cosmological zoom-in simulations from the FIRE project, ranging from dwarf galaxies to galaxies sufficiently massive to host luminous quasars. The simulations model multichannel stellar feedback, but neglect AGN feedback. We show that multiple physical properties, including halo mass, galaxy stellar mass, and depth of the central gravitational potential correlate with accelerated BH fuelling: constant thresholds in these properties are typically crossed within ∼0.1 Hubble time of accelerated BH fuelling. Black hole masses increase sharply when the stellar surface density in the inner 1 kpc crosses a threshold |$\Sigma^\star _{1\,\rm kpc}\approx 10^{9.5} \, {\rm M_{\odot }}\,{\rm kpc}^{-2}$|⁠ , a characteristic value above which gravity prevents stellar feedback from ejecting gas, and similar to the value above which galaxies are observed to quench. We further show that accelerated BH growth correlates with the emergence of long-lived thin gas discs, as well as with virialization of the inner circumgalactic medium. The halo mass M halo ∼ 1012 M⊙ and stellar mass M * ∼ 1010.5 M⊙ at which BH growth accelerates correspond to ∼ L ⋆ galaxies. The fact that stellar feedback becomes inefficient at ejecting gas from the nucleus above this mass scale may play an important role in explaining why AGN feedback appears to be most important in galaxies above L ⋆. [ABSTRACT FROM AUTHOR]

Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the Feedback In Realistic Environment (FIRE)-2 version of the FIRE project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numerical methods, here we update the stellar evolution tracks used to determine stellar feedback inputs, e.g. stellar mass-loss (O/B and AGB), spectra (luminosities and ionization rates), and supernova rates (core-collapse and Ia), as well as detailed mass-dependent yields. We also update the low-temperature cooling and chemistry, to enable improved accuracy at |$T \lesssim 10^{4}\,$| K and densities |$n\gg 1\, {\rm cm^{-3}}$|⁠ , and the meta-galactic ionizing background. All of these synthesize newer empirical constraints on these quantities and updated stellar evolution and yield models from a number of groups, addressing different aspects of stellar evolution. To make the updated models as accessible as possible, we provide fitting functions for all of the relevant updated tracks, yields, etc, in a form specifically designed so they can be directly 'plugged in' to existing galaxy formation simulations. We also summarize the default FIRE-3 implementations of 'optional' physics, including spectrally resolved cosmic rays and supermassive black hole growth and feedback. [ABSTRACT FROM AUTHOR]

We present the first simulations evolving resolved spectra of cosmic rays (CRs) from MeV–TeV energies (including electrons, positrons, (anti)protons, and heavier nuclei), in live kinetic-magnetohydrodynamics galaxy simulations with star formation and feedback. We utilize new numerical methods including terms often neglected in historical models, comparing Milky Way analogues with phenomenological scattering coefficients ν to Solar-neighbourhood [Local interstellar medium (LISM)] observations (spectra, B/C, e+/e−, |$\mathrm{\bar{p}}/\mathrm{p}$|⁠ , 10Be/9Be, ionization, and γ-rays). We show it is possible to reproduce observations with simple single-power-law injection and scattering coefficients (scaling with rigidity R), similar to previous (non-dynamical) calculations. We also find: (1) The circumgalactic medium in realistic galaxies necessarily imposes an |$\sim 10\,$| kpc CR scattering halo, influencing the required ν(R). (2) Increasing the normalization of ν(R) re-normalizes CR secondary spectra but also changes primary spectral slopes, owing to source distribution and loss effects. (3) Diffusive/turbulent reacceleration is unimportant and generally sub-dominant to gyroresonant/streaming losses, which are sub-dominant to adiabatic/convective terms dominated by |$\sim 0.1-1\,$|  kpc turbulent/fountain motions. (4) CR spectra vary considerably across galaxies; certain features can arise from local structure rather than transport physics. (5) Systematic variation in CR ionization rates between LISM and molecular clouds (or Galactic position) arises naturally without invoking alternative sources. (6) Abundances of CNO nuclei require most CR acceleration occurs around when reverse shocks form in SNe, not in OB wind bubbles or later Sedov–Taylor stages of SNe remnants. [ABSTRACT FROM AUTHOR]

The physics of magnetic fields (B) and cosmic rays (CRs) have recently been included in simulations of galaxy formation. However, significant uncertainties remain in how these components affect galaxy evolution. To understand their common observational tracers, we analyse the magnetic fields in a set of high-resolution, magnetohydrodynamic, cosmological simulations of Milky-Way-like galaxies from the FIRE-2 project. We compare mock observables of magnetic field tracers for simulations with and without CRs to observations of Zeeman splitting and rotation/dispersion measures. We find reasonable agreement between simulations and observations in both the neutral and the ionized interstellar medium (ISM). We find that the simulated galaxies with CRs show weaker ISM | B | fields on average compared to their magnetic-field-only counterparts. This is a manifestation of the effects of CRs in the diffuse, low density inner circumgalactic medium (CGM). We find that equipartition between magnetic and cosmic ray energy densities may be valid at large (> 1 kpc) scales for typical ISM densities of Milky-Way-like galaxies, but not in their haloes. Within the ISM, the magnetic fields in our simulated galaxies follow a power-law scaling with gas density. The scaling extends down to neutral hydrogen number densities < 300 cm−3, in contrast to observationally derived models, but consistent with the observational measurements. Finally, we generate synthetic rotation measure (RM) profiles for projections of the simulated galaxies and compare to observational constraints in the CGM. While consistent with upper limits, improved data are needed to detect the predicted CGM RMs at 10–200 kpc and better constrain theoretical predictions. [ABSTRACT FROM AUTHOR]

We use FIRE simulations to study disc formation in z ∼ 0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar discs is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation prior to joining the galaxy. Among galaxies with a high fraction (⁠|$\gt 70{{\ \rm per\ cent}}$|⁠) of their young stars in a thin disc (h / R ∼ 0.1), we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (≳20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to |$\lesssim 10^4\, {\rm K}$|⁠ ; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disc, and while cooling transitions from a quasi-spherical spatial configuration to a more-flattened, disc-like configuration. We show that the existence of this 'rotating cooling flow' accretion mode is strongly correlated with the fraction of stars forming in a thin disc, using a sample of 17 z ∼ 0 galaxies spanning a halo mass range of 1010.5 M⊙ ≲ M h ≲ 1012 M⊙ and stellar mass range of 108 M⊙ ≲ M ⋆ ≲ 1011 M⊙. Notably, galaxies with a thick disc or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion on to the galaxy are likely a necessary condition for the formation of thin, star-forming disc galaxies in a ΛCDM universe. [ABSTRACT FROM AUTHOR]

Recent strides have been made developing dust evolution models for galaxy formation simulations but these approaches vary in their assumptions and degree of complexity. Here, we introduce and compare two separate dust evolution models (labelled 'Elemental' and 'Species'), based on recent approaches, incorporated into the gizmo code and coupled with fire -2 stellar feedback and interstellar medium physics. Both models account for turbulent dust diffusion, stellar production of dust, dust growth via gas-dust accretion, and dust destruction from time-resolved supernovae, thermal sputtering in hot gas, and astration. The 'Elemental' model tracks the evolution of generalized dust species and utilizes a simple, 'tunable' dust growth routine, while the 'Species' model tracks the evolution of specific dust species with set chemical compositions and incorporates a physically motivated, two-phase dust growth routine. We test and compare these models in an idealized Milky Way-mass galaxy and find that while both produce reasonable galaxy-integrated dust-to-metals (D/Z) ratios and predict gas-dust accretion as the main dust growth mechanism, a chemically motivated model is needed to reproduce the observed scaling relation between individual element depletions and D/Z with column density and local gas density. We also find the inclusion of theoretical metallic iron and O-bearing dust species are needed in the case of specific dust species in order to match observations of O and Fe depletions, and the integration of a sub-resolution dense molecular gas/CO scheme is needed to both match observed C depletions and ensure carbonaceous dust is not overproduced in dense environments. [ABSTRACT FROM AUTHOR]

The thermal Sunyaev–Zel'dovich (tSZ) effect is a powerful tool with the potential for constraining directly the properties of the hot gas that dominates dark matter halos because it measures pressure and thus thermal energy density. Studying this hot component of the circumgalactic medium (CGM) is important because it is strongly impacted by star formation and active galactic nucleus (AGN) activity in galaxies, participating in the feedback loop that regulates star and black hole mass growth in galaxies. We study the tSZ effect across a wide halo-mass range using three cosmological hydrodynamical simulations: Illustris-TNG, EAGLE, and FIRE-2. Specifically, we present the scaling relation between the tSZ signal and halo mass and the (mass-weighted) radial profiles of gas density, temperature, and pressure for all three simulations. The analysis includes comparisons to Planck tSZ observations and to the thermal pressure profile inferred from the Atacama Cosmology Telescope (ACT) measurements. We compare these tSZ data to simulations to interpret the measurements in terms of feedback and accretion processes in the CGM. We also identify as-yet unobserved potential signatures of these processes that may be visible in future measurements, which will have the capability of measuring tSZ signals to even lower masses. We also perform internal comparisons between runs with different physical assumptions. We conclude (1) there is strong evidence for the impact of feedback at R500, but that this impact decreases by 5R500, and (2) the thermodynamic profiles of the CGM are highly dependent on the implemented model, such as cosmic-ray or AGN feedback prescriptions. [ABSTRACT FROM AUTHOR]

Without additional heating, radiative cooling of the halo gas of massive galaxies (Milky Way-mass and above) produces cold gas or stars exceeding that observed. Heating from active galactic nucleus (AGN) jets is likely required, but the jet properties remain unclear. This is particularly challenging for galaxy simulations, where the resolution is orders-of-magnitude insufficient to resolve jet formation and evolution. On such scales, the uncertain parameters include the jet energy form [kinetic, thermal, cosmic ray (CR)]; energy, momentum, and mass flux; magnetic fields; opening angle; precession; and duty cycle. We investigate these parameters in a |$10^{14}\, {\rm M}_{\odot }$| halo using high-resolution non-cosmological magnetohydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We explore which scenarios qualitatively meet observational constraints on the halo gas and show that CR-dominated jets most efficiently quench the galaxy by providing CR pressure support and modifying the thermal instability. Mildly relativistic (∼MeV or ∼1010K) thermal plasma jets work but require ∼10 times larger energy input. For fixed energy flux, jets with higher specific energy (longer cooling times) quench more effectively. For this halo mass, kinetic jets are inefficient at quenching unless they have wide opening or precession angles. Magnetic fields also matter less except when the magnetic energy flux reaches ≳ 1044 erg s−1 in a kinetic jet model, which significantly widens the jet cocoon. The criteria for a successful jet model are an optimal energy flux and a sufficiently wide jet cocoon with a long enough cooling time at the cooling radius. [ABSTRACT FROM AUTHOR]

We present the first set of cosmological baryonic zoom-in simulations of galaxies including dissipative self-interacting dark matter (dSIDM). These simulations utilize the Feedback In Realistic Environments galaxy formation physics, but allow the dark matter to have dissipative self-interactions analogous to standard model forces, parametrized by the self-interaction cross-section per unit mass, (σ/ m), and the dimensionless degree of dissipation, 0 < f diss < 1. We survey this parameter space, including constant and velocity-dependent cross-sections, and focus on structural and kinematic properties of dwarf galaxies with |$M_{\rm halo} \sim 10^{10-11}{\, \rm M_\odot }$| and |$M_{\ast } \sim 10^{5-8}{\, \rm M_\odot }$|⁠. Central density profiles (parametrized as ρ ∝  r α) of simulated dwarfs become cuspy when |$(\sigma /m)_{\rm eff} \gtrsim 0.1\, {\rm cm^{2}\, g^{-1}}$| (and f diss = 0.5 as fiducial). The power-law slopes asymptote to α ≈ −1.5 in low-mass dwarfs independent of cross-section, which arises from a dark matter 'cooling flow'. Through comparisons with dark matter only simulations, we find the profile in this regime is insensitive to the inclusion of baryons. However, when |$(\sigma /m)_{\rm eff} \ll 0.1\, {\rm cm^{2}\, g^{-1}}$|⁠ , baryonic effects can produce cored density profiles comparable to non-dissipative cold dark matter (CDM) runs but at smaller radii. Simulated galaxies with |$(\sigma /m) \gtrsim 10\, {\rm cm^{2}\, g^{-1}}$| and the fiducial f diss develop significant coherent rotation of dark matter, accompanied by halo deformation, but this is unlike the well-defined thin 'dark discs' often attributed to baryon-like dSIDM. The density profiles in this high cross-section model exhibit lower normalizations given the onset of halo deformation. For our surveyed dSIDM parameters, halo masses and galaxy stellar masses do not show appreciable difference from CDM, but dark matter kinematics and halo concentrations/shapes can differ. [ABSTRACT FROM AUTHOR]

We present STARFORGE (STAR FORmation in Gaseous Environments): a new numerical framework for 3D radiation magnetohydrodynamic (MHD) simulations of star formation that simultaneously follow the formation, accretion, evolution, and dynamics of individual stars in massive giant molecular clouds (GMCs), while accounting for stellar feedback, including jets, radiative heating and momentum, stellar winds, and supernovae. We use the gizmo code with the MFM mesh-free Lagrangian MHD method, augmented with new algorithms for gravity, time-stepping, sink particle formation and accretion, stellar dynamics, and feedback coupling. We survey a wide range of numerical parameters/prescriptions for sink formation and accretion and find very small variations in star formation history and the IMF (except for intentionally unphysical variations). Modules for mass-injecting feedback (winds, SNe, and jets) inject new gas elements on the fly, eliminating the lack of resolution in diffuse feedback cavities otherwise inherent in Lagrangian methods. The treatment of radiation uses GIZMO's radiative transfer solver to track five frequency bands (IR, optical, NUV, FUV, ionizing), coupling direct stellar emission and dust emission with gas heating and radiation pressure terms. We demonstrate accurate solutions for SNe, winds, and radiation in problems with known similarity solutions, and show that our jet module is robust to resolution and numerical details, and agrees well with previous AMR simulations. STARFORGE can scale up to massive (>105 M⊙) GMCs on current supercomputers while predicting the stellar (≳0.1 M⊙) range of the IMF, permitting simulations of both high- and low-mass cluster formation in a wide range of conditions. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We investigate the spatial structure and evolution of star formation and the interstellar medium (ISM) in interacting galaxies. We use an extensive suite of parsec-scale galaxy-merger simulations (stellar mass ratio = 2.5:1), which employs the 'Feedback In Realistic Environments-2' model (fire-2). This framework resolves star formation, feedback processes, and the multiphase structure of the ISM. We focus on the galaxy-pair stages of interaction. We find that close encounters substantially augment cool (H  i) and cold-dense (H2) gas budgets, elevating the formation of new stars as a result. This enhancement is centrally concentrated for the secondary galaxy, and more radially extended for the primary. This behaviour is weakly dependent on orbital geometry. We also find that galaxies with elevated global star formation rate (SFR) experience intense nuclear SFR enhancement, driven by high levels of either star formation efficiency (SFE) or available cold-dense gas fuel. Galaxies with suppressed global SFR also contain a nuclear cold-dense gas reservoir, but low SFE levels diminish SFR in the central region. Concretely, in the majority of cases, SFR enhancement in the central kiloparsec is fuel-driven (55 per cent for the secondary, 71 per cent for the primary) – while central SFR suppression is efficiency-driven (91 per cent for the secondary, 97 per cent for the primary). Our numerical predictions underscore the need of substantially larger, and/or merger-dedicated, spatially resolved galaxy surveys – capable of examining vast and diverse samples of interacting systems – coupled with multiwavelength campaigns aimed to capture their internal ISM structure. [ABSTRACT FROM AUTHOR]

Observations of ultraviolet (UV) metal absorption lines have provided insight into the structure and composition of the circumgalactic medium (CGM) around galaxies. We compare these observations with the low-redshift (z ≤ 0.3) CGM around dwarf galaxies in high-resolution cosmological zoom-in runs in the FIRE-2 (Feedback In Realistic Environments) simulation suite. We select simulated galaxies that match the halo mass, stellar mass, and redshift of the observed samples. We produce absorption measurements using trident for UV transitions of C  iv , O  vi , Mg  ii , and Si  iii. The FIRE equivalent width (EW) distributions and covering fractions for the C  iv ion are broadly consistent with observations inside 0.5 R vir, but are underpredicted for O  vi , Mg  ii , and Si  iii. The absorption strengths of the ions in the CGM are moderately correlated with the masses and star formation activity of the galaxies. The correlation strengths increase with the ionization potential of the ions. The structure and composition of the gas from the simulations exhibit three zones around dwarf galaxies characterized by distinct ion column densities: the discy interstellar medium, the inner CGM (the wind-dominated regime), and the outer CGM (the IGM accretion-dominated regime). We find that the outer CGM in the simulations is nearly but not quite supported by thermal pressure, so it is not in hydrostatic equilibrium, resulting in halo-scale bulk inflow and outflow motions. The net gas inflow rates are comparable to the star formation rate of the galaxy, but the bulk inflow and outflow rates are greater by an order of magnitude, with velocities comparable to the virial velocity of the halo. These roughly virial velocities (⁠|${\sim } 100 \, \rm km\, s^{-1}$|⁠) produce large EWs in the simulations. This supports a picture for dwarf galaxies in which the dynamics of the CGM at large scales are coupled to the small-scale star formation activity near the centre of their haloes. [ABSTRACT FROM AUTHOR]

We study the spatially resolved (sub-kpc) gas velocity dispersion (σ)–star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way-mass disc galaxies at late times (z ≈ 0). In agreement with observations, we find a relatively flat relationship, with σ ≈ 15–30 km s−1 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 σ. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher σ at constant SFR. The limits of the σ–Σ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 time-scales 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 disc (Toomre's Q  = 1) model predictions, and the simulation data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e. 1 per cent per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of σ, it appears to be subdominant to stellar feedback in the simulated galaxy discs at z ≈ 0. [ABSTRACT FROM AUTHOR]

It has been established for decades that rotation curves deviate from the Newtonian gravity expectation given baryons alone below a characteristic acceleration scale |$g_{\dagger }\sim 10^{-8}\, \rm {cm\, s^{-2}}$|⁠ , a scale promoted to a new fundamental constant in MOND. In recent years, theoretical and observational studies have shown that the star formation efficiency (SFE) of dense gas scales with surface density, SFE ∼ Σ/Σcrit with |$\Sigma _{\rm crit} \sim \langle \dot{p}/m_{\ast }\rangle /(\pi \, G)\sim 1000\, \rm {M_{\odot }\, pc^{-2}}$| (where |$\langle \dot{p}/m_{\ast }\rangle$| is the momentum flux output by stellar feedback per unit stellar mass in a young stellar population). We argue that the SFE, more generally, should scale with the local gravitational acceleration, i.e. that SFE |${\sim}g_{\rm tot}/g_{\rm crit}\equiv (G\, M_{\rm tot}/R^{2}) / \langle \dot{p}/m_{\ast }\rangle$|⁠ , where M tot is the total gravitating mass and |$g_{\rm crit}=\langle \dot{p}/m_{\ast }\rangle = \pi \, G\, \Sigma _{\rm crit} \approx 10^{-8}\, \rm {cm\, s^{-2}} \approx \mathit{ g}_{\dagger }$|⁠. Hence, the observed g † may correspond to the characteristic acceleration scale above which stellar feedback cannot prevent efficient star formation, and baryons will eventually come to dominate. We further show how this may give rise to the observed acceleration scaling |$g_{\rm obs}\sim (g_{\rm baryon}\, g_{\dagger })^{1/2}$| (where g baryon is the acceleration due to baryons alone) and flat rotation curves. The derived characteristic acceleration g † can be expressed in terms of fundamental constants (gravitational constant, proton mass, and Thomson cross-section): |$g_{\dagger }\sim 0.1\, G\, m_{\mathrm{ p}}/\sigma _{\rm T}$|⁠. [ABSTRACT FROM AUTHOR]

Protostars and young stars are strongly spatially 'clustered' or 'correlated' within their natal giant molecular clouds. We demonstrate that such clustering leads to the conclusion that the incident bolometric radiative flux upon a random young star/disc is enhanced (relative to volume-averaged fluxes) by a factor that increases with the total stellar mass of the complex. Because the Galactic cloud mass function is top-heavy, the typical star in our Galaxy experienced a much stronger radiative environment than those forming in well-observed nearby (but relatively small) clouds, exceeding fluxes in the Orion Nebular Cluster by factors of ≳30. Heating of the circumstellar disc around a median young star is dominated by this external radiation beyond |$\sim \! 50$| au. And if discs are not well shielded by ambient dust, external ultraviolet irradiation can dominate over the host star down to sub-au scales. Another consequence of stellar clustering is an extremely broad Galaxy-wide distribution of incident flux (spanning >10 decades), with half the Galactic star formation in a substantial 'tail' towards even more intense background radiation. We also show that the strength of external irradiation is amplified superlinearly in high-density environments such as the Galactic Centre, starbursts, or high-redshift galaxies. [ABSTRACT FROM AUTHOR]

We study stellar-halo formation using six Milky-Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that |$5{-}40{{\ \rm per\ cent}}$| 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 superbubble 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 |$\sim \!10 {{\ \rm per\ cent}}$| of our kinematically identified halo stars were born in outflows; the fraction rises to as high as |$\sim \!40{{\ \rm per\ cent}}$| for the most metal-rich local halo stars ([Fe/H] >−0.5). Such stars can be retrograde and create features similar to the recently discovered Milky Way 'Splash' in phase space. 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. [ABSTRACT FROM AUTHOR]

We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z ≥ 5 galaxies from the Feedback In Realistic Environments project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over |$10^4\, \mathrm{ M}_{\odot }\, {\rm 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 d N /d M ∼ M −2. The cluster formation efficiency is similar across galaxies with stellar masses of ∼107– |$10^{10}\, \mathrm{ M}_{\odot }$| at z ≥ 5. The age spread of cluster stars is typically a few Myr and increases with cluster mass. The metallicity dispersion of cluster members is ∼0.08 dex in |$\rm [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. [ABSTRACT FROM AUTHOR]

We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf (⁠|$M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$|⁠ , |$M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$|⁠) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved (⁠|$\gtrsim 1\,$|  pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs (⁠|$M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$|⁠ , |$M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$|⁠), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses (⁠|$M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$|⁠) and z ≲ 1–2, CRs can suppress SF and stellar masses by factors ∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients |$\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$|⁠. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around |$\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$|⁠ , CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR 'heating' (from collisional and streaming losses) is never dominant. [ABSTRACT FROM AUTHOR]

Radiative feedback (RFB) from stars plays a key role in galaxies, but remains poorly understood. We explore this using high-resolution, multifrequency radiation-hydrodynamics (RHD) simulations from the Feedback In Realistic Environments (FIRE) project. We study ultrafaint dwarf through Milky Way mass scales, including H+He photoionization; photoelectric, Lyman Werner, Compton, and dust heating; and single+multiple scattering radiation pressure (RP). We compare distinct numerical algorithms: ray-based LEBRON (exact when optically thin) and moments-based M1 (exact when optically thick). The most important RFB channels on galaxy scales are photoionization heating and single-scattering RP: in all galaxies, most ionizing/far-UV luminosity (∼1/2 of lifetime-integrated bolometric) is absorbed. In dwarfs, the most important effect is photoionization heating from the UV background suppressing accretion. In MW-mass galaxies, metagalactic backgrounds have negligible effects; but local photoionization and single-scattering RP contribute to regulating the galactic star formation efficiency and lowering central densities. Without some RFB (or other 'rapid' FB), resolved GMCs convert too-efficiently into stars, making galaxies dominated by hyperdense, bound star clusters. This makes star formation more violent and 'bursty' when SNe explode in these hyperclustered objects: thus, including RFB 'smoothes' SFHs. These conclusions are robust to RHD methods, but M1 produces somewhat stronger effects. Like in previous FIRE simulations, IR multiple-scattering is rare (negligible in dwarfs, |$\sim 10{{\ \rm per\ cent}}$| of RP in massive galaxies): absorption occurs primarily in 'normal' GMCs with AV  ∼ 1. [ABSTRACT FROM AUTHOR]

We study a suite of extremely high-resolution cosmological Feedback in Realistic Environments simulations of dwarf galaxies (⁠|$M_{\rm halo} \lesssim 10^{10}\rm \, M_{\odot }$|⁠), run to z  = 0 with |$30\, \mathrm{M}_{\odot }$| resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with |$M_{\rm halo} \gtrsim 10^{8.6}\, \mathrm{M}_{\odot }$| is populated by a resolved stellar galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; |$M_{\ast }\lt 10^{5}\, \mathrm{M}_{\odot }$|⁠) have their star formation (SF) truncated early (z ≳ 2), likely by reionization, while classical dwarfs (⁠|$M_{\ast }\gt 10^{5}\, \mathrm{M}_{\odot }$|⁠) continue forming stars to z < 0.5. The systems have bursty star formation histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with M */ M halo > 10−4 to form a dark matter core |${\gt}200\rm \, pc$|⁠ , while lower mass UFDs exhibit cusps down to |${\lesssim}100\rm \, pc$|⁠ , as expected from energetic arguments. Our dwarfs with |$M_{\ast }\gt 10^{4}\, \mathrm{M}_{\odot }$| have half-mass radii (R 1/2) in agreement with Local Group (LG) dwarfs (dynamical mass versus R 1/2 and stellar rotation also resemble observations). The lowest mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model). [ABSTRACT FROM AUTHOR]

We present a suite of FIRE -2 cosmological zoom-in simulations of isolated field dwarf galaxies, all with masses of |$M_{\rm halo} \approx 10^{10}\, {\rm M}_{\odot }$| at z = 0, across a range of dark matter models. For the first time, we compare how both self-interacting dark matter (SIDM) and/or warm dark matter (WDM) models affect the assembly histories as well as the central density structure in fully hydrodynamical simulations of dwarfs. Dwarfs with smaller stellar half-mass radii (r 1/2 < 500 pc) have lower σ⋆/ V max ratios, reinforcing the idea that smaller dwarfs may reside in haloes that are more massive than is naively expected. The majority of dwarfs simulated with self-interactions actually experience contraction of their inner density profiles with the addition of baryons relative to the cores produced in dark-matter-only runs, though the simulated dwarfs are always less centrally dense than in ΛCDM. The V 1/2– r 1/2 relation across all simulations is generally consistent with observations of Local Field dwarfs, though compact objects such as Tucana provide a unique challenge. Overall, the inclusion of baryons substantially reduces any distinct signatures of dark matter physics in the observable properties of dwarf galaxies. Spatially resolved rotation curves in the central regions (<400 pc) of small dwarfs could provide a way to distinguish between CDM, WDM, and SIDM, however: at the masses probed in this simulation suite, cored density profiles in dwarfs with small r 1/2 values can only originate from dark matter self-interactions. [ABSTRACT FROM AUTHOR]

We study star formation histories (SFHs) of 500 dwarf galaxies (stellar mass |$M_\ast =10^5\!-\!10^9\, \rm {M}_\odot$|⁠) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher mass dwarfs quench later (if at all), regardless of environment. We also identify differences between the environments, both in terms of 'satellite versus central' and 'LG versus individual MW versus isolated dwarf central.' Around the individual MW-mass hosts, we recover the result expected from environmental quenching: central galaxies in the 'near field' have more extended SFHs than their satellite counterparts, with the former more closely resemble isolated (true field) dwarfs (though near-field centrals are still somewhat earlier forming). However, this difference is muted in the LG-like environments, where both near-field centrals and satellites have similar SFHs, which resemble satellites of single MW-mass hosts. This distinction is strongest for M * = 106– |$10^7\, \rm {M}_\odot$| but exists at other masses. Our results suggest that the paired halo nature of the LG may regulate star formation in dwarf galaxies even beyond the virial radii of the MW and Andromeda. Caution is needed when comparing zoom-in simulations targeting isolated dwarf galaxies against observed dwarf galaxies in the LG. [ABSTRACT FROM AUTHOR]

We present the late-time evolution of m12m, a cosmological simulation of a Milky Way-like galaxy from the FIRE project. The simulation forms a bar after redshift z  = 0.2. We show that the evolution of the model exhibits behaviours typical of kinematic fractionation, with a bar weaker in older populations, an X-shape traced by the younger, metal-rich populations, and a prominent X-shape in the edge-on mean metallicity map. Because of the late formation of the bar in m12m, stars forming after |$10\mbox{$\:{\rm Gyr}$}$| (z  = 0.34) significantly contaminate the bulge, at a level higher than is observed at high latitudes in the Milky Way, implying that its bar cannot have formed as late as in m12m. We also study the model's vertex deviation of the velocity ellipsoid as a function of stellar metallicity and age in the equivalent of Baade's Window. The formation of the bar leads to a non-zero vertex deviation. We find that metal-rich stars have a large vertex deviation (∼40°), which becomes negligible for metal-poor stars, a trend also found in the Milky Way, despite not matching in detail. We demonstrate that the vertex deviation also varies with stellar age and is large for stars as old as |$9 \mbox{$\:{\rm Gyr}$}$|⁠, while |$13\mbox{$\:{\rm Gyr}$}$| old stars have negligible vertex deviation. When we exclude stars that have been accreted, the vertex deviation is not significantly changed, demonstrating that the observed variation of vertex deviation with metallicity is not necessarily due to an accreted population. [ABSTRACT FROM AUTHOR]

We present a comprehensive suite of high-resolution (parsec-scale), idealized (non-cosmological) galaxy merger simulations (24 runs, stellar mass ratio ∼2.5:1) to investigate the connection between interaction-induced star formation and the evolution of the interstellar medium (ISM) in various temperature–density regimes. We use the gizmo code and the second version of the 'Feedback in Realistic Environments' model (FIRE-2), which captures the multiphase structure of the ISM. Our simulations are designed to represent galaxy mergers in the local Universe. In this work, we focus on the 'galaxy-pair period' between first and second pericentric passage. We split the ISM into four regimes: hot, warm, cool, and cold-dense, motivated by the hot, ionized, atomic and molecular gas phases observed in real galaxies. We find that, on average, interactions enhance the star formation rate of the pair (⁠|${\sim }30{{\ \rm per\ cent}}$|⁠, merger-suite sample average) and elevate their cold-dense gas content (⁠|${\sim }18{{\ \rm per\ cent}}$|⁠). This is accompanied by a decrease in warm gas (⁠|${\sim }11{{\ \rm per\ cent}}$|⁠), a negligible change in cool gas (⁠|${\sim }4{{\ \rm per\ cent}}$| increase), and a substantial increase in hot gas (⁠|${\sim }400{{\ \rm per\ cent}}$|⁠). The amount of cold-dense gas with densities above 1000 cm−3 (the cold ultra-dense regime) is elevated significantly (⁠|${\sim }240{{\ \rm per\ cent}}$|⁠), but only accounts for |${\sim }0.15{{\ \rm per\ cent}}$| (on average) of the cold-dense gas budget. [ABSTRACT FROM AUTHOR]

We study the impact of a warm dark matter (WDM) cosmology on dwarf galaxy formation through a suite of cosmological hydrodynamical zoom-in simulations of |$M_{\rm halo} \approx 10^{10}\, \mathrm{ M}_{\odot }$| dark matter haloes as part of the Feedback in Realistic Environments (FIRE) project. A main focus of this paper is to evaluate the combined effects of dark matter physics and stellar feedback on the well-known small-scale issues found in cold dark matter (CDM) models. We find that the |$z$|  = 0 stellar mass of a galaxy is strongly correlated with the central density of its host dark matter halo at the time of formation, |$z$| f, in both CDM and WDM models. WDM haloes follow the same M ⋆(⁠|$z$|  = 0)– V max(⁠|$z$| f) relation as in CDM, but they form later, are less centrally dense, and therefore contain galaxies that are less massive than their CDM counterparts. As a result, the impact of baryonic effects on the central gravitational potential is typically diminished relative to CDM. However, the combination of delayed formation in WDM and energy input from stellar feedback results in dark matter profiles with lower overall densities. The WDM galaxies studied here have a wider diversity of star formation histories (SFHs) than the same systems simulated in CDM, and the two lowest M ⋆ WDM galaxies form all of their stars at late times. The discovery of young ultrafaint dwarf galaxies with no ancient star formation – which do not exist in our CDM simulations – would therefore provide evidence in support of WDM. [ABSTRACT FROM AUTHOR]

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code ('FIRE-1') for consistency. Motivated by the development of more accurate numerics – including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms – and exploration of new physics (e.g. magnetic fields), we introduce 'FIRE-2', an updated numerical implementation of FIRE physics for the gizmo code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star formation algorithm, cooling physics, and chemistry have weak effects provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media. Central (∼kpc) mass concentrations in massive (> L *) galaxies are sensitive to numerics (via trapping/recycling of winds in hot haloes). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion on to dwarfs and instantaneous star formation in discs. We provide all initial conditions and numerical algorithms used. [ABSTRACT FROM AUTHOR]

We study the implementation of mechanical feedback from supernovae (SNe) and stellar mass loss in galaxy simulations, within the Feedback In Realistic Environments (FIRE) project. We present the FIRE-2 algorithm for coupling mechanical feedback, which can be applied to any hydrodynamics method (e.g. fixed-grid, moving-mesh, and mesh-less methods), and black hole as well as stellar feedback. This algorithm ensures manifest conservation of mass, energy, and momentum, and avoids imprinting 'preferred directions' on the ejecta. We show that it is critical to incorporate both momentum and thermal energy of mechanical ejecta in a self-consistent manner, accounting for SNe cooling radii when they are not resolved. Using idealized simulations of single SN explosions, we show that the FIRE-2 algorithm, independent of resolution, reproduces converged solutions in both energy and momentum. In contrast, common 'fully thermal' (energy-dump) or 'fully kinetic' (particle-kicking) schemes in the literature depend strongly on resolution: when applied at mass resolution ≳100 M⊙, they diverge by orders of magnitude from the converged solution. In galaxy-formation simulations, this divergence leads to orders-of-magnitude differences in galaxy properties, unless those models are adjusted in a resolution-dependent way. We show that all models that individually time-resolve SNe converge to the FIRE-2 solution at sufficiently high resolution (<100 M⊙). However, in both idealized single-SN simulations and cosmological galaxy-formation simulations, the FIRE-2 algorithm converges much faster than other sub-grid models without re-tuning parameters. [ABSTRACT FROM AUTHOR]

Using a state-of-the-art cosmological simulation of merging proto-galaxies at high redshift from the FIRE project, with explicit treatments of star formation and stellar feedback in the interstellar medium, we investigate the formation of star clusters and examine one of the formation hypotheses of present-day metal-poor globular clusters.We find that frequent mergers in high-redshift proto-galaxies could provide a fertile environment to produce long-lasting bound star clusters. The violent merger event disturbs the gravitational potential and pushes a large gas mass of ≳ 105-6M☉ collectively to high density, at which point it rapidly turns into stars before stellar feedback can stop star formation. The high dynamic range of the reported simulation is critical in realizing such dense star-forming clouds with a small dynamical time-scale, tff ≲3 Myr, shorter than most stellar feedback time-scales. Our simulation then allows us to trace how clusters could become virialized and tightly bound to survive for up to ~420 Myr till the end of the simulation. Because the cluster's tightly bound core was formed in one short burst, and the nearby older stars originally grouped with the cluster tend to be preferentially removed, at the end of the simulation the cluster has a small age spread. [ABSTRACT FROM AUTHOR]

We investigate stellar metallicity distribution functions (MDFs), including Fe and α-element abundances, in dwarf galaxies from the Feedback in Realistic Environment (FIRE) project.We examine both isolated dwarf galaxies and those that are satellites of a MilkyWay-mass galaxy. In particular, we study the effects of including a sub-grid turbulent model for the diffusion of metals in gas. Simulations that include diffusion have narrower MDFs and abundance ratio distributions, because diffusion drives individual gas and star particles towards the average metallicity. This effect provides significantly better agreement with observed abundance distributions in dwarf galaxies in the Local Group, including small intrinsic scatter in [α/Fe] versus [Fe/H] of ≲0.1 dex. This small intrinsic scatter arises in our simulations because the interstellar medium in dwarf galaxies is well mixed at nearly all cosmic times, such that stars that form at a given time have similar abundances to ≲0.1 dex. Thus, most of the scatter in abundances at ɀ = 0 arises from redshift evolution and not from instantaneous scatter in the ISM.We find similar MDF widths and intrinsic scatter for satellite and isolated dwarf galaxies, which suggests that environmental effects play a minor role compared with internal chemical evolution in our simulations. Overall, with the inclusion of metal diffusion, our simulations reproduce abundance distribution widths of observed low-mass galaxies, enabling detailed studies of chemical evolution in galaxy formation. [ABSTRACT FROM AUTHOR]

We compare a suite of four simulated dwarf galaxies formed in 1010 M☉ haloes of collisionless cold dark matter (CDM) with galaxies simulated in the same haloes with an identical galaxy formation model but a non-zero cross-section for DM self-interactions. These cosmological zoom-in simulations are part of the Feedback In Realistic Environments (FIRE) project and utilize the FIRE-2 model for hydrodynamics and galaxy formation physics. We find the stellar masses of the galaxies formed in self-interacting dark matter (SIDM) with σ/m = 1 cm2 g-1 are very similar to those in CDM (spanning M* ≈ 105.7-7.0 M☉) and all runs lie on a similar stellar mass-size relation. The logarithmic DM density slope (α = d log ρ/d log r) in the central 250-500 pc remains steeper than α = -0.8 for the CDM-Hydro simulations with stellar mass M* ~ 106.6M☉ and core-like in the most massive galaxy. In contrast, every SIDM hydrodynamic simulation yields a flatter profile, withα >-0.4. Moreover, the central density profiles predicted in SIDM runs without baryons are similar to the SIDM runs that include FIRE- 2 baryonic physics. Thus, SIDM appears to bemuch more robust to the inclusion of (potentially uncertain) baryonic physics than CDM on this mass scale, suggesting that SIDM will be easier to falsify than CDM using low-mass galaxies. Our FIRE simulations predict that galaxies less massive than M* ≲ 3 × 106 M☉ provide potentially ideal targets for discriminating models, with SIDM producing substantial cores in such tiny galaxies and CDM producing cusps. [ABSTRACT FROM AUTHOR]

We present a suite of 15 cosmological zoom-in simulations of isolated dark matter haloes, all with masses of Mhalo ≈ 1010M⊚ at z = 0, in order to understand the relationship among halo assembly, galaxy formation and feedback’s effects on the central density structure in dwarf galaxies. These simulations are part of the Feedback in Realistic Environments (FIRE) project and are performed at extremely high resolution (mbaryon = 500M⊚, mdm = 2500M⊚). The resultant galaxies have stellar masses that are consistent with rough abundance matching estimates, coinciding with the faintest galaxies that can be seen beyond the virial radius of the Milky Way (M*/M⊚ ≈ 105 - 107). This non-negligible spread in stellar mass at z = 0 in haloes within a narrow range of virial masses is strongly correlated with central halo density or maximum circular velocity Vmax, both of which are tightly linked to halo formation time. Much of this dependence of M* on a second parameter (beyond Mhalo) is a direct consequence of the Mhalo ∼ 1010M⊚ mass scale coinciding with the threshold for strong reionization suppression: the densest, earliest-forming haloes remain above the UV-suppression scale throughout their histories while late-forming systems fall below the UV-suppression scale over longer periods and form fewer stars as a result. In fact, the latest-forming, lowest-concentration halo in our suite fails to form any stars. Haloes that form galaxies with M* ≳ 2 × 106 M⊚ have reduced central densities relative to dark-matter-only simulations, and the radial extent of the density modifications is well-approximated by the galaxy half-mass radius r1/2. Lower-mass galaxies do not modify their host dark matter haloes at the mass scale studied here. This apparent stellar mass threshold of M* ≈ 2 × 106-2 × 10-4 Mhalo is broadly consistent with previous work and provides a testable prediction of FIRE feedback models in Λcold dark matter. [ABSTRACT FROM AUTHOR]

Using high-resolution simulations with explicit treatment of stellar feedback physics based on the FIRE (Feedback In Realistic Environments) project, we study how galaxy formation and the interstellar medium (ISM) are affected by magnetic fields, anisotropic Spitzer-Braginskii conduction and viscosity, and sub-grid metal diffusion from unresolved turbulence. We consider controlled simulations of isolated (non-cosmological) galaxies but also a limited set of cosmological 'zoom-in' simulations. Although simulations have shown significant effects from these physics with weak or absent stellar feedback, the effects are much weaker than those of stellar feedback when the latter is modelled explicitly. The additional physics have no systematic effect on galactic star formation rates (SFRs). In contrast, removing stellar feedback leads to SFRs being overpredicted by factors of ~10-100. Without feedback, neither galactic winds nor volume-filling hot-phase gas exist, and discs tend to runaway collapse to ultra-thin scaleheights with unphysically dense clumps congregating at the galactic centre. With stellar feedback, a multi-phase, turbulent medium with galactic fountains and winds is established. At currently achievable resolutions and for the investigated halo mass range 1010-1013Mʘ, the additional physics investigated here (magnetohydrodynamic, conduction, viscosity, metal diffusion) have only weak (~10 per cent-level) effects on regulating SFR and altering the balance of phases, outflows or the energy in ISM turbulence, consistent with simple equipartition arguments. We conclude that galactic star formation and the ISM are primarily governed by a combination of turbulence, gravitational instabilities and feedback. We add the caveat that active galactic nucleus feedback is not included in the present work. [ABSTRACT FROM AUTHOR]

We use cosmological simulations from the FIRE (Feedback In Realistic Environments) project to study the baryon cycle and galaxy mass assembly for central galaxies in the halo mass range Mhalo ~ 1010-1013 M☉. By tracing cosmic inflows, galactic outflows, gas recycling and merger histories, we quantify the contribution of physically distinct sources of material to galaxy growth. We show that in situ star formation fuelled by fresh accretion dominates the early growth of galaxies of all masses, while the re-accretion of gas previously ejected in galactic winds often dominates the gas supply for a large portion of every galaxy's evolution. Externally processedmaterial contributes increasingly to the growth of central galaxies at lower redshifts. This includes stars formed ex situ and gas delivered by mergers, as well as smooth intergalactic transfer of gas from other galaxies, an important but previously underappreciated growth mode. By z=0, wind transfer, i.e. the exchange of gas between galaxies via winds, can dominate gas accretion on to ~L* galaxies over fresh accretion and standard wind recycling. Galaxies of all masses re-accrete ≳50 per cent of the gas ejected in winds and recurrent recycling is common. The total mass deposited in the intergalactic medium per unit stellar mass formed increases in lower mass galaxies. Re-accretion of wind ejecta occurs over a broad range of time-scales, with median recycling times (~100-350 Myr) shorter than previously found. Wind recycling typically occurs at the scale radius of the halo, independent of halo mass and redshift, suggesting a characteristic recycling zone around galaxies that scales with the size of the inner halo and the galaxy's stellar component. [ABSTRACT FROM AUTHOR]

We present an analysis of the flow of metals through the circumgalactic medium (CGM) in the Feedback in Realistic Environments (FIRE) simulations of galaxy formation, ranging from isolated dwarfs to L* galaxies. We find that nearly all metals produced in high-redshift galaxies are carried out in winds that reach 0.25Rvir. When measured at 0.25Rvir the metallicity of outflows is slightly higher than the interstellar medium (ISM) metallicity. Many metals thus reside in the CGM. Cooling and recycling from this reservoir determine the metal budget in the ISM. The outflowing metal flux decreases by a factor of ∼2–5 between 0.25Rvir and Rvir. Furthermore, outflow metallicity is typically lower at Rvir owing to dilution of the remaining outflow by metal-poor material swept up from the CGM. The inflow metallicity at Rvir is generally low, but outflow and inflow metallicities are similar in the inner halo. At low redshift, massive galaxies no longer generate outflows that reach the CGM, causing a divergence in CGM and ISM metallicity. Dwarf galaxies continue to generate outflows, although they preferentially retain metal ejecta. In all but the least massive galaxy considered, amajority of themetals are within the halo at z=0. We measure the fraction ofmetals in CGM, ISM and stars, and quantify the thermal state of CGM metals in each halo. The total amount of metals in the low-redshift CGM of two simulated L* galaxies is consistent with estimates from the Cosmic Origin Spectrograph haloes survey, while for the other two it appears to be lower. [ABSTRACT FROM AUTHOR]

We perform a systematic Bayesian analysis of rotation versus dispersion support (vrot/σ) in 40 dwarf galaxies throughout the local volume (LV) over a stellar mass range of 103.5M⊙ < M* < 108M⊙. We find that the stars in ~80 per cent of the LV dwarf galaxies studied - both satellites and isolated systems - are dispersion-supported. In particular, we show that 6/10 isolated dwarfs in our sample have vrot/σ ≤ 1.0, while all have vrot/σ ≤ 2.0. These results challenge the traditional view that the stars in gas-rich dwarf irregulars (dIrrs) are distributed in cold, rotationally supported stellar discs, while gas-poor dwarf spheroidals (dSphs) are kinematically distinct in having dispersion-supported stars. We see no clear trend between vrot/σ and distance to the closest L* galaxy, nor between vrot/σ and M* within our mass range. We apply the same Bayesian analysis to four FIRE hydrodynamic zoom-in simulations of isolated dwarf galaxies (109M⊙ < Mvir < 1010M⊙) and show that the simulated isolated dIrr galaxies have stellar ellipticities and stellar vrot/σ ratios that are consistent with the observed population of dIrrs and dSphs without the need to subject these dwarfs to any external perturbations or tidal forces. We posit that most dwarf galaxies form as puffy, dispersiondominated systems, rather than cold, angular-momentum-supported discs. If this is the case, then transforming a dIrr into a dSph may require little more than removing its gas. [ABSTRACT FROM AUTHOR]

We present the mufasa suite of cosmological hydrodynamic simulations, which employs the gizmo meshless finite mass (MFM) code including H2-based star formation, nine-element chemical evolution, two-phase kinetic outflows following scalings from the Feedback in Realistic Environments zoom simulations, and evolving halo mass-based quenching. Our fiducial (50h-1 Mpc)³ volume is evolved to z = 0 with a quarter billion elements. The predicted galaxy stellar mass functions (GSMFs) reproduces observations from z = 4 → 0 to ... 1.2σ in cosmic variance, providing an unprecedented match to this key diagnostic. The cosmic star formation history and stellar mass growth show general agreement with data, with a strong archaeological downsizing trend such that dwarf galaxies form the majority of their stars after z ~ 1. We run 25 and 12.5 h-1 Mpc volumes to z = 2 with identical feedback prescriptions, the latter resolving all hydrogen-cooling haloes, and the three runs display fair resolution convergence. The specific star formation rates broadly agree with data at z = 0, but are underpredicted at z ~ 2 by a factor of 3, re-emphasizing a longstanding puzzle in galaxy evolution models. We compare runs using MFM and two flavours of smoothed particle hydrodynamics, and show that the GSMF is sensitive to hydrodynamics methodology at the ~ x 2 level, which is sub-dominant to choices for parametrizing feedback. [ABSTRACT FROM AUTHOR]

A popular theory of star formation is gravito-turbulent fragmentation, in which self-gravitating structures are created by turbulence-driven density fluctuations. Simple theories of isothermal fragmentation successfully reproduce the core mass function (CMF) which has a very similar shape to the initial mass function (IMF) of stars. However, numerical simulations of isothermal turbulent fragmentation thus far have not succeeded in identifying a fragment mass scale that is independent of the simulation resolution. Moreover, the fluid equations for magnetized, selfgravitating, isothermal turbulence are scale-free, and do not predict any characteristic mass. In this paper we show that, although an isothermal self-gravitating flow does produce a CMF with a mass scale imposed by the initial conditions, this scale changes as the parent cloud evolves. In addition, the cores that form undergo further fragmentation and after sufficient time forget about their initial conditions, yielding a scale-free pure power-law distribution dN/dM ∝ M-2 for the stellar IMF. We show that this problem can be alleviated by introducing additional physics that provides a termination scale for the cascade. Our candidate for such physics is a simple model for stellar radiation feedback. Radiative heating, powered by accretion on to forming stars, arrests the fragmentation cascade and imposes a characteristic mass scale that is nearly independent of the time-evolution or initial conditions in the star-forming cloud, and that agrees well with the peak of the observed IMF. In contrast, models that introduce a stiff equation of state for denser clouds but that do not explicitly include the effects of feedback do not yield an invariant IMF. [ABSTRACT FROM AUTHOR]

We use high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environment (FIRE) project to study the galaxy mass-metallicity relations (MZR) from z = 0-6. These simulations include explicit models of the multiphase ISM, star formation, and stellar feedback. The simulations cover halo masses Mhalo = 109-1013M⊙ and stellar masses M﹡ = 104-1011M⊙ at z=0 and have been shown to produce many observed galaxy properties from z=0-6. For the first time, our simulations agree reasonably well with the observed mass- metallicity relations at z = 0-3 for a broad range of galaxy masses. We predict the evolution of the MZR from z = 0-6, as log(Zgas/Z⊙) = 12 + log(O/H) - 9.0 = 0.35[log(M﹡/M⊙) - 10] + 0.93 exp(-0.43z) - 1.05 and log(Z﹡/Z⊙) = [Fe/H] + 0.2 = 0.40[log(M﹡/M⊙) - 10] + 0.67 exp(-0.50z) - 1.04, for gas-phase and stellar metallicity, respectively. Our simulations suggest that the evolution of MZR is associated with the evolution of stellar/gas mass fractions at different redshifts, indicating the existence of a universal metallicity relation between stellar mass, gas mass, and metallicities. In our simulations, galaxies above M﹡ = 106M⊙ are able to retain a large fraction of their metals inside the halo, because metal-rich winds fail to escape completely and are recycled into the galaxy. This resolves a longstanding discrepancy between 'subgrid' wind models (and semi-analytic models) and observations, where common subgrid models cannot simultaneously reproduce the MZR and the stellar mass functions. [ABSTRACT FROM AUTHOR]

We present a series of high-resolution (20-2000 M⊙, 0.1-4 pc) cosmological zoom-in simulations at z ≳ 6 from the Feedback In Realistic Environment (FIRE) project. These simulations cover halo masses 109-1011 M⊙ and rest-frame ultraviolet magnitude MUV = -9 to -19. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback, which produce reasonable galaxy properties at z = 0-6. We post-process the snapshots with a radiative transfer code to evaluate the escape fraction (fesc) of hydrogen ionizing photons. We find that the instantaneous fesc has large time variability (0.01-20 percent), while the time-averaged fesc over long time-scales generally remains ≲ 5 percent, considerably lower than the estimate in many reionization models. We find no strong dependence of fesc on galaxy mass or redshift. In our simulations, the intrinsic ionizing photon budgets are dominated by stellar populations younger than 3 Myr, which tend to be buried in dense birth clouds. The escaping photons mostly come from populations between 3 and 10 Myr, whose birth clouds have been largely cleared by stellar feedback. However, these populations only contribute a small fraction of intrinsic ionizing photon budgets according to standard stellar population models. We show that fesc can be boosted to high values, if stellar populations older than 3 Myr produce more ionizing photons than standard stellar population models (as motivated by, e.g. models including binaries). By contrast, runaway stars with velocities suggested by observations can enhance fesc by only a small fraction. We show that 'sub-grid' star formation models, which do not explicitly resolve star formation in dense clouds with n » 1 cm-3, will dramatically overpredict fesc. [ABSTRACT FROM AUTHOR]