Search

Fuelling star formation in large, discy galaxies requires a continuous supply of gas accreting into star-forming regions. Previously, we characterized this accretion in four Milky Way mass galaxies (⁠|$M_{\rm halo}\sim 10^{12}{\rm M}_{\odot }$|⁠) in the FIRE-2 cosmological zoom-in simulations. At |$z\sim 0$|⁠ , we found that gas within the inner circumgalactic medium (iCGM) approaches the disc with comparable angular momentum (AM) to the disc edge, joining in the outer half of the gaseous disc. Within the disc, gas moves inwards at velocities of |$\sim$| 1–5 km s |$^{-1}$| while fully rotationally supported. In this study, we analyse the torques that drive these flows. In all cases studied, we find that the torques in discs enable gas accreted near the disc edge to transport inwards and fuel star formation in the central few kpc. The primary sources of torque come from gravity, hydrodynamical forces, and the sub-grid |$P \mathrm{ d}V$| work done by supernova (SN) remnants interacting with gas on |$\lesssim$| 10 pc scales. These SNe remnant interactions induce negative torques within the inner disc and positive torques in the outer disc. The gas–gas gravitational, hydro, and 'feedback' torques transfer AM outwards to where accreting gas joins the disc, playing an important role in driving inflows and regulating disc structure. Gravitational torques from stars and dark matter provide an AM sink within the innermost regions of the disc and iCGM, respectively. Feedback torques are dominant within the disc, while gravitational and hydrodynamical torques have similar significance depending on the system/region. Torques from viscous shearing, magnetic forces, stellar winds, and radiative transfer are less significant. [ABSTRACT FROM AUTHOR]

Atomic hydrogen (H  i) serves a crucial role in connecting galactic-scale properties such as star formation with the large-scale structure of the Universe. While recent numerical simulations have successfully matched the observed covering fraction of H  i near Lyman Break Galaxies (LBGs) and in the foreground of luminous quasars at redshifts |$z \lesssim 3$|⁠ , the low-mass end remains as-of-yet unexplored in observational and computational surveys. We employ a cosmological, hydrodynamical simulation (FIREbox) supplemented with zoom-in simulations (MassiveFIRE) from the Feedback In Realistic Environments (FIRE) project to investigate the H  i covering fraction of Lyman Limit Systems (⁠|$N_{{\text{H}}\, \rm{{\small I}}} \gtrsim 10^{17.2}$|  cm |$^{-2}$|⁠) across a wide range of redshifts (⁠|$z=0-6$|⁠) and halo masses (⁠|$10^8-10^{13} \, \,\mathrm{ M}_{\odot }$| at |$z=0$|⁠ , |$10^8-10^{11}\, \,\mathrm{ M}_{\odot }$| at |$z=6$|⁠) in the absence of feedback from active galactic nuclei. We find that the covering fraction inside haloes exhibits a strong increase with redshift, with only a weak dependence on halo mass for higher mass haloes. For massive haloes (⁠|$M_{\mathrm{vir}} \sim 10^{11}-10^{12} \,\mathrm{ M}_{\odot }$|⁠), the radial profiles showcase scale-invariance and remain independent of mass. The radial dependence is well captured by a fitting function. The covering fractions in our simulations are in good agreement with measurements of the covering fraction in LBGs. Our comprehensive analysis unveils a complex dependence with redshift and halo mass for haloes with |$M_{\mathrm{vir}} \lesssim 10^{10} \,\mathrm{ M}_{\odot }$| that future observations aim to constrain, providing key insights into the physics of structure formation and gas assembly. [ABSTRACT FROM AUTHOR]

Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H |$_2$| gas mass of the region using the CO-to-H |$_2$| conversion factor (⁠|$X_{\rm{{\small CO}}}$|⁠). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulations of isolated disc galaxies, ranging from dwarfs to Milky Way-mass galaxies, that were run using the fire-2 subgrid models coupled to the chimes non-equilibrium chemistry solver. We use the non-equilibrium abundances from the simulations, and we also compare to results using abundances assuming equilibrium, which we calculate from the simulation in post-processing. Our non-equilibrium simulations are able to reproduce the relation between CO and H |$_2$| column densities, and the relation between |$X_{\rm{{\small CO}}}$| and metallicity, seen within observations of the Milky Way. We also compare to the xCOLD GASS survey, and find agreement with their data to our predicted CO luminosities at fixed star formation rate. We also find the multivariate function used by xCOLD GASS overpredicts the H |$_2$| mass for our simulations, motivating us to suggest an alternative multivariate function of our fitting, though we caution that this fitting is uncertain due to the limited range of galaxy conditions covered by our simulations. We also find that the non-equilibrium chemistry has little effect on the conversion factor (<5 per cent) for our high-mass galaxies, though still affects the H |$_2$| mass and |$L_{\rm{{\small CO}}}$| by |$\approx$| 25  per cent. [ABSTRACT FROM AUTHOR]

Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, feedback, and cosmological inflows. While theoretical models ascribe the shape of the MZR to the lower efficiency of galactic outflows in more massive galaxies, what drives its evolution remains an open question. In this letter, we analyse how the MZR evolves over z  = 0–3, combining results from the FIREbox cosmological volume simulation with analytical models. Contrary to a frequent assertion in the literature, we find that the evolution of the gas fraction does not contribute significantly to the redshift evolution of the MZR. Instead, we show that the latter is driven by the redshift dependence of the inflow metallicity, outflow metallicity, and mass loading factor, whose relative importance depends on stellar mass. These findings also suggest that the evolution of the MZR is not explained by galaxies moving along a fixed surface in the space spanned by stellar mass, gas-phase metallicity, and star formation rate. [ABSTRACT FROM AUTHOR]

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. [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 study the kinematics of stars both at their formation and today within 14 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. We quantify the relative importance of cosmological disc settling and post-formation dynamical heating. We identify three eras: a Pre-Disc Era (typically ≳ 8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disc Era (≈8–4 Gyr ago), when stars started to form on rotation-dominated orbits but with high velocity dispersion, σform; and a Late-Disc Era (≲ 4 Gyr ago), when stars formed with low σform. σform increased with time during the Pre-Disc Era, peaking ≈8 Gyr ago, then decreased throughout the Early-Disc Era as the disc settled and remained low throughout the Late-Disc Era. By contrast, the dispersion measured today, σnow, increases monotonically with age because of stronger post-formation heating for Pre-Disc stars. Importantly, most of σnow was in place at formation, not added post-formation, for stars younger than ≈10 Gyr. We compare the evolution of the three velocity components: at all times, σR, form > σϕ, form > σZ, form. Post-formation heating primarily increased σR at ages ≲ 4 Gyr but acted nearly isotropically for older stars. The kinematics of young stars in FIRE-2 broadly agree with the range observed across the MW, M31, M33, and PHANGS-MUSE galaxies. The lookback time that the disc began to settle correlates with its dynamical state today: earlier-settling galaxies currently form colder discs. Including stellar cosmic-ray feedback does not significantly change disc rotational support at fixed stellar mass. [ABSTRACT FROM AUTHOR]

Both observations and simulations have shown strong evidence for highly time-variable star formation in low-mass and/or high-redshift galaxies, which has important observational implications because high-redshift galaxy samples are rest-ultraviolet (rest-UV) selected and therefore particularly sensitive to the recent star formation. Using a suite of cosmological 'zoom-in' simulations at z > 5 from the Feedback in Realistic Environments project, we examine the implications of bursty star formation histories for observations of high-redshift galaxies with JWST. We characterize how the galaxy observability depends on the star formation history. We also investigate selection effects due to bursty star formation on the physical properties measured, such as the gas fraction, specific star formation rate, and metallicity. We find the observability to be highly time-dependent for galaxies near the survey's limiting flux due to the star formation rate variability: as the star formation rate fluctuates, the same galaxy oscillates in and out of the observable sample. The observable fraction |$f_\mathrm{obs} = 50~{{\ \rm per\ cent}}$| at z ∼ 7 and M ⋆ ∼ 108.5– |$10^{9}\, {\rm M}_{\odot }$| for a JWST /NIRCam survey reaching a limiting magnitude of |$m^\mathrm{lim}_\mathrm{AB} \sim 29{\!-\!}30$|⁠ , representative of surveys such as JADES and CEERS. JWST -detectable galaxies near the survey limit tend to have properties characteristic of galaxies in the bursty phase: on average, they show approximately 2.5 times higher cold, dense gas fractions and 20 times higher specific star formation rates at a given stellar mass than galaxies below the rest-UV detection threshold. Our study represents a first step in quantifying selection effects and the associated biases due to bursty star formation in studying high-redshift galaxy properties. [ABSTRACT FROM AUTHOR]

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

Low-mass galaxies are highly susceptible to environmental effects that can efficiently quench star formation. We explore the role of ram pressure in quenching low-mass galaxies (⁠|$M_{*}\sim 10^{5}{-}10^{9}\, \rm {M}_{\odot }$|⁠) within 2 Mpc of Milky Way (MW) hosts using the FIRE-2 simulations. Ram pressure is highly variable across different environments, within individual MW haloes, and for individual low-mass galaxies over time. The impulsiveness of ram pressure – the maximum ram pressure scaled to the integrated ram pressure prior to quenching – correlates with whether a galaxy is quiescent or star forming. The time-scale between maximum ram pressure and quenching is anticorrelated with impulsiveness, such that high impulsiveness corresponds to quenching time-scales <1 Gyr. Galaxies in low-mass groups (⁠|$M_\mathrm{*,host}10^{7}{-}10^{9}\, \rm {M}_{\odot }$|⁠) outside of MW haloes experience typical ram pressure only slightly lower than ram pressure on MW satellites, helping to explain effective quenching via group preprocessing. Ram pressure on MW satellites rises sharply with decreasing distance to the host, and, at a fixed physical distance, more recent pericentre passages are typically associated with higher ram pressure because of greater gas density in the inner host halo at late times. Furthermore, the ram pressure and gas density in the inner regions of Local Group-like paired host haloes are higher at small angles off the host galaxy disc compared to isolated hosts. The quiescent fraction of satellites within these low-latitude regions is also elevated in the simulations and observations, signaling possible anisotropic quenching via ram pressure around MW-mass hosts. [ABSTRACT FROM AUTHOR]

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

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

Understanding the star formation rate (SFR) variability and how it depends on physical properties of galaxies is important for developing and testing the theory of galaxy formation. We investigate how statistical measurements of the extragalactic background light (EBL) can shed light on this topic and complement traditional methods based on observations of individual galaxies. Using semi-empirical models of galaxy evolution and SFR indicators sensitive to different star formation time-scales (e.g. H α and ultraviolet continuum luminosities), we show that the SFR variability, quantified by the joint probability distribution of the SFR indicators (i.e. the bivariate conditional luminosity function), can be characterized as a function of galaxy mass and redshift through the cross-correlation between deep, near-infrared maps of the EBL and galaxy distributions. As an example, we consider combining upcoming SPHEREx maps of the EBL with galaxy samples from Rubin Observatory Legacy Survey of Space and Time. We demonstrate that their cross-correlation over a sky fraction of f sky ∼ 0.5 can constrain the joint SFR indicator distribution at high significance up to z ∼ 2.5 for mass-complete samples of galaxies down to |$M_{*}\sim 10^9\, {\rm M}_{\odot }$|⁠. These constraints not only allow models of different SFR variability to be distinguished, but also provide unique opportunities to investigate physical mechanisms that require large number statistics such as environmental effects. The cross-correlations investigated illustrate the power of combining cosmological surveys to extract information inaccessible from each data set alone, while the large galaxy populations probed capture ensemble-averaged properties beyond the reach of targeted observations towards individual galaxies. [ABSTRACT FROM AUTHOR]

Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated haloes. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution function only require resolutions sufficient to keep the jet-inflated cavities stable. Comparing different model variations, as well as an independent jet model using a different hydrodynamics code, we show that the dominant uncertainties are the choices of jet properties within a given model. Independent of implementation, we find that light, thermal jets with low momentum flux tend to delay the onset of a cooling flow more efficiently on a 50 Myr time-scale than heavy, kinetic jets. The delay of the cooling flow originates from a displacement and boost in entropy of the central gas. If the jet kinetic luminosity depends on accretion rate, collimated, light, hydrodynamic jets are able to reduce cooling flows in haloes, without a need for jet precession or wide opening angles. Comparing the jet feedback with a 'kinetic wind' implementation shows that equal amounts of star formation rate reduction can be achieved by different interactions with the halo gas: the jet has a larger effect on the hot halo gas while leaving the denser, star-forming phase in place, while the wind acts more locally on the star-forming phase, which manifests itself in different time-variability properties. [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]

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

Recent observations and simulations indicate substantial evolution in the properties of galaxies with time, wherein rotationally supported and steady thin discs (like those frequently observed in the local Universe) emerge from galaxies that are clumpy, irregular, and have bursty star formation rates (SFRs). To better understand the progenitors of local disc galaxies, we carry out an analysis of three FIRE-2 simulated galaxies with a mass similar to the Milky Way at redshift z  = 0. We show that all three galaxies transition from bursty to steady SFRs at a redshift between z  = 0.5 and z  = 0.8, and that this transition coincides with the rapid (≲1 Gyr) emergence of a rotationally supported interstellar medium (ISM). In the late phase with steady SFR, the rotational energy comprises |${\gtrsim }90{{\ \rm per\ cent}}$| of the total kinetic + thermal energy in the ISM, and is roughly half the gravitational energy. By contrast, during the early bursty phase, the ISM initially has a quasi-spheroidal morphology and its energetics are dominated by quasi-isotropic in- and outflows out of virial equilibrium. The subdominance of rotational support and out-of-equilibrium conditions at early times challenge the application of standard equilibrium disc models to high-redshift progenitors of Milky Way-like galaxies. We further find that the formation of a rotationally-supported ISM coincides with the onset of a thermal pressure supported inner circumgalactic medium (CGM). Before this transition, there is no clear boundary between the ISM and the inner CGM. [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]

The properties of young star clusters formed within a galaxy are thought to vary in different interstellar medium conditions, but the details of this mapping from galactic to cluster scales are poorly understood due to the large dynamic range involved in galaxy and star cluster formation. We introduce a new method for modelling cluster formation in galaxy simulations: mapping giant molecular clouds (GMCs) formed self-consistently in a FIRE-2 magnetohydrodynamic galaxy simulation on to a cluster population according to a GMC-scale cluster formation model calibrated to higher resolution simulations, obtaining detailed properties of the galaxy's star clusters in mass, metallicity, space, and time. We find |$\sim 10{{\ \rm per\ cent}}$| of all stars formed in the galaxy originate in gravitationally bound clusters overall, and this fraction increases in regions with elevated Σgas and ΣSFR, because such regions host denser GMCs with higher star formation efficiency. These quantities vary systematically over the history of the galaxy, driving variations in cluster formation. The mass function of bound clusters varies – no single Schechter-like or power-law distribution applies at all times. In the most extreme episodes, clusters as massive as 7 × 106 M⊙ form in massive, dense clouds with high star formation efficiency. The initial mass–radius relation of young star clusters is consistent with an environmentally dependent 3D density that increases with Σgas and ΣSFR. The model does not reproduce the age and metallicity statistics of old (⁠|$\gt 11\rm Gyr$|⁠) globular clusters found in the Milky Way, possibly because it forms stars more slowly at z > 3. [ABSTRACT FROM AUTHOR]

Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbulence, metallicity, interstellar radiation field, simulation geometry and turbulent driving. In our fiducial cloud, the raw simulation data reproduces the observed multiplicity fractions for Solar-type and higher mass stars, similar to previous works. However, after correcting for observational incompleteness the simulation underpredicts these values. The discrepancy is likely due to the lack of disc fragmentation, as the simulation only resolves multiples that form either through capture or core fragmentation. The raw mass distribution of companions is consistent with randomly drawing from the initial mass function for the companions of |$\gt 1\, \mathrm{M}_{\rm \odot }$| stars. However, accounting for observational incompleteness produces a flatter distribution similar to observations. We show that stellar multiplicity changes as the cloud evolves and anticorrelates with stellar density. This relationship also explains most multiplicity variations between runs, i.e. variations in the initial conditions that increase stellar density (increased surface density, reduced turbulence) also act to decrease multiplicity. While other parameters, such as metallicity, interstellar radiation, and geometry significantly affect the star formation history or the IMF, varying them produces no clear trend in stellar multiplicity properties. [ABSTRACT FROM AUTHOR]

Type Ia supernovae are critical for feedback and elemental enrichment in galaxies. Recent surveys like the All-Sky Automated Survey for Supernova (ASAS-SN) and the Dark Energy Survey (DES) find that the specific supernova Ia rate at z ∼ 0 may be ≲20–50× higher in lower mass galaxies than at Milky Way-mass. Independently, observations show that the close-binary fraction of solar-type Milky Way stars is higher at lower metallicity. Motivated by these observations, we use the FIRE-2 cosmological zoom-in simulations to explore the impact of metallicity-dependent rate models on galaxies of |$M_* \sim 10^7\!-\!10^{11}\, \rm {M}_{\odot }$|⁠. First, we benchmark our simulated star formation histories against observations, and show that the assumed stellar mass functions play a major role in determining the degree of tension between observations and metallicity-independent rate models, potentially causing ASAS-SN and DES observations to agree more than might appear. Models in which the supernova Ia rate increases with decreasing metallicity (⁠|$\propto Z^{-0.5 \,\, \rm {to} \,\, -1}$|⁠) provide significantly better agreement with observations. Encouragingly, these rate increases (≳10× in low-mass galaxies) do not significantly impact galaxy masses and morphologies, which remain largely unaffected except for our most extreme models. We explore implications for both [Fe/H] and [ |$\alpha /\rm {Fe}$| ] enrichment; metallicity-dependent rate models can improve agreement with the observed stellar mass–metallicity relations in low-mass galaxies. Our results demonstrate that a range of metallicity-dependent rate models are viable for galaxy formation and motivate future work. [ABSTRACT FROM AUTHOR]

We investigate stellar elemental abundance patterns at |$z$|  = 0 in eight low-mass (⁠|$M_{*}=10^{6}{-}10^{9}\ \text{M}_{\odot }$|⁠) galaxies in the Feedback in Realistic Environments cosmological simulations. Using magnesium (Mg) as a representative α-element, we explore stellar abundance patterns in magnesium-to-iron ([Mg/Fe]) versus iron-to-hydrogen ([Fe/H]), which follow an overall monotonic trend that evolved slowly over time. Additionally, we explore three notable secondary features in enrichment (in three different case-study galaxies) that arise from a galaxy merger or bursty star formation. First, we observe a secondary track with a lower [Mg/Fe] than the main trend. At |$z$|  = 0, stars from this track are predominantly found within 2–6 kpc of the centre; they were accreted in a 1:3 total-mass-ratio merger ∼0.4 Gyr ago. Second, we find a distinct elemental bimodality that forms following a strong burst in star formation in a galaxy at |$t_{\text{lookback}}\, \sim 10$| Gyr. This burst quenched star formation for ∼0.66 Gyr, allowing Type Ia supernovae to enrich the system with iron (Fe) before star formation resumed. Third, we examine stripes in enrichment that run roughly orthogonal to the dominant [Mg/Fe] versus [Fe/H] trend; these stripes correspond to short bursts of star formation during which core-collapse supernovae enrich the surrounding medium with Mg (and Fe) on short time-scales. If observed, these features would substantiate the utility of elemental abundances in revealing the assembly and star-formation histories of dwarf galaxies. We explore the observability of these features for upcoming spectroscopic studies. Our results show that precise measurements of elemental abundance patterns can reveal critical events in the formation histories of low-mass galaxies. [ABSTRACT FROM AUTHOR]

We analyse the first giant molecular cloud (GMC) simulation to follow the formation of individual stars and their feedback from jets, radiation, winds, and supernovae, using the STARFORGE framework in the GIZMO code. We evolve the GMC for |$\sim 9 \rm Myr$|⁠ , from initial turbulent collapse to dispersal by feedback. Protostellar jets dominate feedback momentum initially, but radiation and winds cause cloud disruption at |$\sim 8{{\ \rm per\ cent}}$| star formation efficiency (SFE), and the first supernova at |$8.3\, \rm Myr$| comes too late to influence star formation significantly. The per-free-fall SFE is dynamic, accelerating from 0 per cent to |$\sim 18{{\ \rm per\ cent}}$| before dropping quickly to <1 per cent, but the estimate from YSO counts compresses it to a narrower range. The primary cluster forms hierarchically and condenses to a brief (⁠|$\sim 1\, \mathrm{Myr}$|⁠) compact (⁠|$\sim 1\, \rm pc$|⁠) phase, but does not virialize before the cloud disperses, and the stars end as an unbound expanding association. The initial mass function resembles the Chabrier (2005) form with a high-mass slope α = −2 and a maximum mass of 55 M⊙. Stellar accretion takes |$\sim 400\, \rm kyr$| on average, but |$\gtrsim 1\,\rm Myr$| for >10 M⊙ stars, so massive stars finish growing latest. The fraction of stars in multiples increase as a function of primary mass, as observed. Overall, the simulation much more closely resembles reality, compared to previous versions that neglected different feedback physics entirely. But more detailed comparison with synthetic observations will be needed to constrain the theoretical uncertainties. [ABSTRACT FROM AUTHOR]

Simulations of isolated giant molecular clouds (GMCs) are an important tool for studying the dynamics of star formation, but their turbulent initial conditions (ICs) are uncertain. Most simulations have either initialized a velocity field with a prescribed power spectrum on a smooth density field (failing to model the full structure of turbulence) or 'stirred' turbulence with periodic boundary conditions (which may not model real GMC boundary conditions). We develop and test a new GMC simulation setup (called turbsphere) that combines advantages of both approaches: we continuously stir an isolated cloud to model the energy cascade from larger scales, and use a static potential to confine the gas. The resulting cloud and surrounding envelope achieve a quasi-equilibrium state with the desired hallmarks of supersonic ISM turbulence (e.g. density PDF and a ∼ k −2 velocity power spectrum), whose bulk properties can be tuned as desired. We use the final stirred state as initial conditions for star formation simulations with self-gravity, both with and without continued driving and protostellar jet feedback, respectively. We then disentangle the respective effects of the turbulent cascade, simulation geometry, external driving, and gravity/MHD boundary conditions on the resulting star formation. Without external driving, the new setup obtains results similar to previous simple spherical cloud setups, but external driving can suppress star formation considerably in the new setup. Periodic box simulations with the same dimensions and turbulence parameters form stars significantly slower, highlighting the importance of boundary conditions and the presence or absence of a global collapse mode in the results of star formation calculations. [ABSTRACT FROM AUTHOR]

Previous studies of fueling black holes in galactic nuclei have argued (on scales |${\sim}0.01{-}1000\,$| pc) accretion is dynamical with inflow rates |$\dot{M}\sim \eta \, M_{\rm gas}/t_{\rm dyn}$| in terms of gas mass M gas, dynamical time t dyn, and some η. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency f wind = M ejected/ M total scales with the gravitational acceleration as |$(1-f_{\rm wind})/f_{\rm wind}\sim \bar{a}_{\rm grav}/\langle \dot{p}/m_{\ast }\rangle \sim \Sigma _{\rm eff}/\Sigma _{\rm crit}$| where |$\bar{a}_{\rm grav}\equiv G\, M_{\rm tot}(\lt r)/r^{2}$| and |$\langle \dot{p}/m_{\ast }\rangle$| is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, |$\eta \rightarrow \eta \, (1-f_{\rm wind})$|⁠ , we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting the slope and normalization of the M BH − σ and M BH − M bulge relation, L AGN −SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient (f wind ∼ 1) below total-mass surface density |$M_{\rm tot}/\pi \, r^{2}\lt \Sigma _{\rm crit}\sim 3\times 10^{9}\, \mathrm{M}_{\odot }\, {\rm kpc^{-2}}$| (where |$\Sigma _{\rm crit}=\langle \dot{p}/m_{\ast }\rangle /(\pi \, G)$|⁠), BH mass is predicted to specifically trace host galaxy properties above a critical surface brightness Σcrit (B-band |$\mu _{\rm B}^{\rm crit}\sim 19\, {\rm mag\, arcsec^{-2}}$|⁠). This naturally explains why BH masses preferentially reflect bulge properties or central surface densities (e.g. |$\Sigma _{1\, {\rm kpc}}$|⁠), not 'total' galaxy properties. [ABSTRACT FROM AUTHOR]

We characterize mass, momentum, energy, and metal outflow rates of multiphase galactic winds in a suite of FIRE-2 cosmological 'zoom-in' simulations from the Feedback in Realistic Environments (FIRE) project. We analyse simulations of low-mass dwarfs, intermediate-mass dwarfs, Milky Way-mass haloes, and high-redshift massive haloes. Consistent with previous work, we find that dwarfs eject about 100 times more gas from their interstellar medium (ISM) than they form in stars, while this mass 'loading factor' drops below one in massive galaxies. Most of the mass is carried by the hot phase (>105 K) in massive haloes and the warm phase (103−105 K) in dwarfs; cold outflows (<103 K) are negligible except in high-redshift dwarfs. Energy, momentum, and metal loading factors from the ISM are of order unity in dwarfs and significantly lower in more massive haloes. Hot outflows have 2−5 × higher specific energy than needed to escape from the gravitational potential of dwarf haloes; indeed, in dwarfs, the mass, momentum, and metal outflow rates increase with radius whereas energy is roughly conserved, indicating swept up halo gas. Burst-averaged mass loading factors tend to be larger during more powerful star formation episodes and when the inner halo is not virialized, but we see effectively no trend with the dense ISM gas fraction. We discuss how our results can guide future controlled numerical experiments that aim to elucidate the key parameters governing galactic winds and the resulting associated preventative feedback. [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 a large suite of magnetohydrodynamic simulations of turbulent, star-forming giant molecular clouds (GMCs) with stellar feedback, extending previous work by simulating 10 different random realizations for each point in the parameter space of cloud mass and size. It is found that once the clouds disperse due to stellar feedback, both self-gravitating star clusters and unbound stars generally remain, which arise from the same underlying continuum of substructured stellar density, i.e. the hierarchical cluster formation scenario. The fraction of stars that are born within gravitationally bound star clusters is related to the overall cloud star formation efficiency set by stellar feedback, but has significant scatter due to stochastic variations in the small-scale details of the star-forming gas flow. We use our numerical results to calibrate a model for mapping the bulk properties (mass, size, and metallicity) of self-gravitating GMCs on to the star cluster populations they form, expressed statistically in terms of cloud-level distributions. Synthesizing cluster catalogues from an observed GMC catalogue in M83, we find that this model predicts initial star cluster masses and sizes that are in good agreement with observations, using only standard IMF and stellar evolution models as inputs for feedback. Within our model, the ratio of the strength of gravity to stellar feedback is the key parameter setting the masses of star clusters, and of the various feedback channels direct stellar radiation (photon momentum and photoionization) is the most important on GMC scales. [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 use FIRE-2 simulations to examine 3D variations of gas-phase elemental abundances of [O/H], [Fe/H], and [N/H] in 11 MW and M31-mass galaxies across their formation histories at z ≤ 1.5 (⁠|$t_{\rm lookback} \le 9.4 \, \rm {Gyr}$|⁠), motivated by characterizing the initial conditions of stars for chemical tagging. Gas within |$1 \, \rm {kpc}$| of the disc mid-plane is vertically homogeneous to |$\lesssim 0.008 \, \rm {dex}$| at all z ≤ 1.5. We find negative radial gradients (metallicity decreases with galactocentric radius) at all times, which steepen over time from |$\approx \! -0.01 \, \rm {dex}\, \rm {kpc}^{-1}$| at z = 1 (⁠|$t_{\rm lookback} = 7.8 \, \rm {Gyr}$|⁠) to |$\approx \! -0.03 \, \rm {dex}\, \rm {kpc}^{-1}$| at z = 0, and which broadly agree with observations of the MW, M31, and nearby MW/M31-mass galaxies. Azimuthal variations at fixed radius are typically |$0.14 \, \rm {dex}$| at z = 1, reducing to |$0.05 \, \rm {dex}$| at z = 0. Thus, over time radial gradients become steeper while azimuthal variations become weaker (more homogeneous). As a result, azimuthal variations were larger than radial variations at z ≳ 0.8 (⁠|$t_{\rm lookback} \gtrsim 6.9 \, \rm {Gyr}$|⁠). Furthermore, elemental abundances are measurably homogeneous (to ≲0.05 dex) across a radial range of |$\Delta R \approx 3.5 \, \rm {kpc}$| at z ≳ 1 and |$\Delta R \approx 1.7 \, \rm {kpc}$| at z = 0. We also measure full distributions of elemental abundances, finding typically negatively skewed normal distributions at z ≳ 1 that evolve to typically Gaussian distributions by z = 0. Our results on gas abundances inform the initial conditions for stars, including the spatial and temporal scales for applying chemical tagging to understand stellar birth in the MW. [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 study the impact of cosmic rays (CRs) on the structure of virial shocks, using a large suite of high-resolution cosmological FIRE-2 simulations accounting for CR injection by supernovae. In Milky Way-mass, low-redshift (z ≲ 1−2) haloes, which are expected to form 'hot haloes' with slowly cooling gas in quasi-hydrostatic equilibrium (with a stable virial shock), our simulations without CRs do exhibit clear virial shocks. The cooler phase condensing out from inflows becomes pressure confined to overdense clumps, embedded in low-density, volume-filling hot gas with volume-weighted cooling time longer than inflow time. The gas thus transitions sharply from cool free-falling inflow, to hot and thermal-pressure supported at approximately the virial radius (≈ R vir), and the shock is quasi-spherical. With CRs, we previously argued that haloes in this particular mass and redshift range build up CR-pressure-dominated gaseous haloes. Here, we show that when CR pressure dominates over thermal pressure, there is no significant virial shock. Instead, inflowing gas is gradually decelerated by the CR pressure gradient and the gas is relatively subsonic out to and even beyond R vir. Rapid cooling also maintains subvirial temperatures in the inflowing gas within ∼ R vir. [ABSTRACT FROM AUTHOR]

This document is a correction to a previous article titled "STARFORGE: Towards a comprehensive numerical model of star cluster formation and feedback." The correction addresses an error in Equation 44, which provides an incorrect expression for the wind mass loss rate of a main-sequence star in the STARFORGE numerical model. The correct formula, as used in the simulations and subsequent works, is provided in Equation (1). The correction is authored by Michael Y Grudić, Dávid Guszejnov, Philip F Hopkins, Stella S R Offner, and Claude-André Faucher-Giguère. [Extracted from the article]

The initial mass function (IMF) of stars is a key quantity affecting almost every field of astrophysics, yet it remains unclear what physical mechanisms determine it. We present the first runs of the STAR FORmation in Gaseous Environments project, using a new numerical framework to follow the formation of individual stars in giant molecular clouds (GMCs) using the gizmo code. Our suite includes runs with increasingly complex physics, starting with isothermal ideal magnetohydrodynamics (MHD) and then adding non-isothermal thermodynamics and protostellar outflows. We show that without protostellar outflows the resulting stellar masses are an order of magnitude too high, similar to the result in the base isothermal MHD run. Outflows disrupt the accretion flow around the protostar, allowing gas to fragment and additional stars to form, thereby lowering the mean stellar mass to a value similar to that observed. The effect of jets upon global cloud evolution is most pronounced for lower mass GMCs and dense clumps, so while jets can disrupt low-mass clouds, they are unable to regulate star formation in massive GMCs, as they would turn an order unity fraction of the mass into stars before unbinding the cloud. Jets are also unable to stop the runaway accretion of massive stars, which could ultimately lead to the formation of stars with masses |${\gt}500\, \mathrm{M}_{\rm \odot }$|⁠. Although we find that the mass scale set by jets is insensitive to most cloud parameters (i.e. surface density, virial parameter), it is strongly dependent on the momentum loading of the jets (which is poorly constrained by observations) as well as the temperature of the parent cloud, which predicts slightly larger IMF variations than observed. We conclude that protostellar jets play a vital role in setting the mass scale of stars, but additional physics are necessary to reproduce the observed IMF. [ABSTRACT FROM AUTHOR]

The relation between infrared excess (IRX) and UV spectral slope (βUV) is an empirical probe of dust properties of galaxies. The shape, scatter, and redshift evolution of this relation are not well understood, however, leading to uncertainties in estimating the dust content and star formation rates (SFRs) of galaxies at high redshift. In this study, we explore the nature and properties of the IRX–βUV relation with a sample of z = 2–6 galaxies (⁠|$M_*\approx 10^9\!-\!10^{12}\, \mathrm{M}_\odot$|⁠) extracted from high-resolution cosmological simulations (MassiveFIRE) of the Feedback in Realistic Environments (FIRE) project. The galaxies in our sample show an IRX–βUV relation that is in good agreement with the observed relation in nearby galaxies. IRX is tightly coupled to the UV optical depth, and is mainly determined by the dust-to-star geometry instead of total dust mass, while βUV is set both by stellar properties, UV optical depth, and the dust extinction law. Overall, much of the scatter in the IRX–βUV relation of our sample is found to be driven by variations of the intrinsic UV spectral slope. We further assess how the IRX–βUV relation depends on viewing direction, dust-to-metal ratio, birth-cloud structures, and the dust extinction law and we present a simple model that encapsulates most of the found dependencies. Consequently, we argue that the reported 'deficit' of the infrared/sub-millimetre bright objects at z ≳ 5 does not necessarily imply a non-standard dust extinction law at those epochs. [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]

Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question by using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disc galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyse how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the mid-plane. Bulk flows (e.g. inflows and fountains) are important at a few disc scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total mid-plane pressure is well-predicted by the weight of the disc gas and we show that it also scales linearly with the star formation rate surface density (ΣSFR). These results support the notion that the Kennicutt–Schmidt relation arises because ΣSFR and the gas surface density (Σg) are connected via the ISM mid-plane pressure. [ABSTRACT FROM AUTHOR]

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

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 revisit the question of 'hot mode' versus 'cold mode' accretion on to galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally supported, R circ, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of |$\sim 10^{12}\, {\rm M_{\odot }}$| derived for baryon-complete haloes. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of |${\dot{M}}_{\rm crit}\approx 0.7(v_{\rm c}/100\, \rm km\ s^{-1})^{5.4}(R_{\rm circ}/10\, {\rm kpc})(Z/\, {\rm Z_{\odot }})^{-0.9}\, {\rm M_{\odot }}\, {\rm yr}^{-1}$|⁠ , where Z and v c are the metallicity and circular velocity measured at R circ. For accretion rates |$\gtrsim {\dot{M}}_{\rm crit}$| the volume-filling gas phase can in principle be 'transonic' – virialized in the outer halo but cool and free-falling near the galaxy. We compare |${\dot{M}}_{\rm crit}$| to the average star formation rate (SFR) in haloes at 0 < z < 10 implied by the stellar-mass–halo-mass relation. For a plausible metallicity evolution with redshift, we find that |${\rm SFR}\lesssim {\dot{M}}_{\rm crit}$| at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of |$\sim 10^{12}\, {\rm M_{\odot }}$|⁠. [ABSTRACT FROM AUTHOR]

The correlation between galaxies' integrated stellar masses and star formation rates (the 'star formation main sequence', SFMS) is a well-established scaling relation. Recently, surveys have found a relationship between the star formation rate (SFR) and stellar mass surface densities on kpc and sub-kpc scales (the 'resolved SFMS', rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in Feedback In Realistic Environments 2 (FIRE-2) zoom-in simulations of Milky Way-mass galaxies. We make SFR and stellar mass maps of the simulated galaxies at a variety of spatial resolutions and star formation averaging time-scales and fit the rSFMS using multiple methods from the literature. While the absolute value of the SFMS slope (αMS) depends on the fitting method, the slope is steeper for longer star formation time-scales and lower spatial resolutions regardless of the fitting method employed. We present a toy model that quantitatively captures the dependence of the simulated galaxies' αMS on spatial resolution and use it to illustrate how this dependence can be used to constrain the characteristic mass of star-forming clumps. [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 explore the radial variation of star formation histories (SFHs) in dwarf galaxies simulated with Feedback In Realistic Environments (FIRE) physics. The sample contains 26 field dwarf galaxies with M star = 105–109 M⊙. We find age gradients are common in our dwarfs, with older stars dominant at large radii. The strength of the gradient correlates with overall galaxy age such that earlier star formation produces a more pronounced gradient. The relation between formation time and strength of the gradient is driven by both mergers and star formation feedback. Mergers can both steepen and flatten the age gradient depending on the timing of the merger and SFHs of the merging galaxy. In galaxies without significant mergers, feedback pushes stars to the outskirts. The strength of the age gradient is determined by the subsequent evolution of the galaxy. Galaxies with weak age gradients constantly grow to z  = 0, meaning that young star formation occurs at a similar radius to which older stars are heated to. In contrast, galaxies with strong age gradients tend to maintain a constant half-mass radius over time. If real galaxies have age gradients as we predict, stellar population studies that rely on sampling a limited fraction of a galaxy can give a biased view of its global SFH. Central fields can be biased young by Gyrs while outer fields are biased old. Fields positioned near the 2D half-light radius will provide the least biased measure of a dwarf galaxy's global SFH. [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]

Dust temperature is an important property of the interstellar medium (ISM) of galaxies. It is required when converting (sub)millimetre broad-band flux to total infrared luminosity (L IR), and hence star formation rate, in high-redshift galaxies. However, different definitions of dust temperatures have been used in the literature, leading to different physical interpretations of how ISM conditions change with, e.g. redshift and star formation rate. In this paper, we analyse the dust temperatures of massive (⁠|$M_{\rm star} \gt 10^{10}\, \mathrm{M}_{\odot }$|⁠) |$z$| = 2–6 galaxies with the help of high-resolution cosmological simulations from the Feedback in Realistic Environments (fire) project. At |$z$| ∼ 2, our simulations successfully predict dust temperatures in good agreement with observations. We find that dust temperatures based on the peak emission wavelength increase with redshift, in line with the higher star formation activity at higher redshift, and are strongly correlated with the specific star formation rate. In contrast, the mass-weighted dust temperature, which is required to accurately estimate the total dust mass, does not strongly evolve with redshift over |$z$| = 2–6 at fixed IR luminosity but is tightly correlated with L IR at fixed |$z$|⁠. We also analyse an 'equivalent' dust temperature for converting (sub)millimetre flux density to total IR luminosity, and provide a fitting formula as a function of redshift and dust-to-metal ratio. We find that galaxies of higher equivalent (or higher peak) dust temperature ('warmer dust') do not necessarily have higher mass-weighted temperatures. A 'two-phase' picture for interstellar dust can explain the different scaling relations of the various dust temperatures. [ABSTRACT FROM AUTHOR]

Measurements of the star formation efficiency (SFE) of giant molecular clouds (GMCs) in the Milky Way generally show a large scatter, which could be intrinsic or observational. We use magnetohydrodynamic simulations of GMCs (including feedback) to forward-model the relationship between the true GMC SFE and observational proxies. We show that individual GMCs trace broad ranges of observed SFE throughout collapse, star formation, and disruption. Low measured SFEs (⁠|${\ll} 1\hbox{ per cent}$|⁠) are 'real' but correspond to early stages; the true 'per-freefall' SFE where most stars actually form can be much larger. Very high (⁠|${\gg} 10\hbox{ per cent}$|⁠) values are often artificially enhanced by rapid gas dispersal. Simulations including stellar feedback reproduce observed GMC-scale SFEs, but simulations without feedback produce 20× larger SFEs. Radiative feedback dominates among mechanisms simulated. An anticorrelation of SFE with cloud mass is shown to be an observational artefact. We also explore individual dense 'clumps' within GMCs and show that (with feedback) their bulk properties agree well with observations. Predicted SFEs within the dense clumps are ∼2× larger than observed, possibly indicating physics other than feedback from massive (main-sequence) stars is needed to regulate their collapse. [ABSTRACT FROM AUTHOR]

We present a suite of 34 high-resolution cosmological zoom-in simulations consisting of thousands of haloes up to |$M_{\rm halo}\sim 10^{12}\, \mathrm{M}_{\odot }$| (⁠|$M_{\ast }\sim 10^{10.5}\, \mathrm{M}_{\odot }$|⁠) at z ≥ 5 from the Feedback in Realistic Environments project. We post-process our simulations with a three-dimensional Monte Carlo dust radiative transfer code to study dust attenuation, dust emission, and dust temperature within these simulated z ≥ 5 galaxies. Our sample forms a tight correlation between infrared excess (IRX ≡ F IR/ F UV) and ultraviolet (UV)-continuum slope (βUV), despite the patchy, clumpy dust geometry shown in our simulations. We find that the IRX–βUV relation is mainly determined by the shape of the attenuation law and is independent of its normalization (set by the dust-to-gas ratio). The bolometric IR luminosity (L IR) correlates with the intrinsic UV luminosity and the star formation rate (SFR) averaged over the past 10 Myr. We predict that at a given L IR, the peak wavelength of the dust spectral energy distributions for z ≥ 5 galaxies is smaller by a factor of 2 (due to higher dust temperatures on average) than at z  = 0. The higher dust temperatures are driven by higher specific SFRs and SFR surface densities with increasing redshift. We derive the galaxy UV luminosity functions (UVLFs) at z  = 5–10 from our simulations and confirm that a heavy attenuation is required to reproduce the observed bright-end UVLFs. We also predict the IR luminosity functions (IRLFs) and UV luminosity densities at z  = 5–10. We discuss the implications of our results on current and future observations probing dust attenuation and emission in z ≥ 5 galaxies. [ABSTRACT FROM AUTHOR]