Search

Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox-HR re-simulates the cosmic volume (L = 22.1 cMpc) of the original FIREbox run with eight times higher mass resolution (m_b ~ 7800 M_sun), but with identical physics, down to z ~ 6. FIREbox-HR predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at z ~ 6 - 14, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox-HR, the SFE - halo mass relation for intermediate mass halos (M_halo ~ 10^9 - 10^11 M_sun) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE - halo mass relation lead to a larger contribution from lower mass halos at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE - halo mass relation inferred from FIREbox-HR allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at z > 12 will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn., Comment: 28 pages, 14 figures, 5 tables, revised to match version accepted by MNRAS, Appendix D added, UV LF shown in Fig. 1 provided as txt file

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 analyze 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., Comment: 6 pages, 4 figures. Submitted for publication on MNRAS

Line-intensity mapping (LIM) experiments coming online now will survey fluctuations in aggregate emission in the [C II] ionized carbon line from galaxies at the end of reionization. Experimental progress must be matched by theoretical reassessments of approaches to modelling and the information content of the signal. We present a new model for the halo-[C II] connection, building upon results from the FIRE simulations suggesting that gas mass and metallicity most directly determine [C II] luminosity. Applying our new model to an ensemble of peak-patch halo lightcones, we generate new predictions for the [C II] LIM signal at $z\gtrsim6$. We expect a baseline 4000-hour LIM survey from the CCAT facility to have the fundamental sensitivity to detect the [C II] power spectrum at a significance of $5\sigma$ at $z\sim6$, with an extended or successor Stage 2 experiment improving significance to $48\sigma$ at $z\sim6$ and achieving $11\sigma$ at $z\sim7.5$. Cross-correlation through stacking, simulated against a mock narrow-band Lyman-break galaxy survey, would yield a strong detection of the radial profile of cosmological [C II] emission surrounding star-forming galaxies. We also analyse the role of a few of our model's parameters through the pointwise relative entropy (PRE) of the distribution of [C II] intensities. While the PRE signature of different model parameters can become degenerate or diminished after factoring in observational distortions, various parameters do imprint themselves differently on the one-point statistics of the intrinsic signal. Further work can pave the way to access this information and distinguish different sources of non-Gaussianity in the [C II] LIM observation., Comment: 16 pages + acknowledgements/bibliography/appendix (18 pages total); 14 figures, one table; v3 incorporates minor changes reflecting version accepted for publication in MNRAS

Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-luminous, starburst galaxies at low and high redshifts, as well as some moderately star-forming galaxies at the epoch of re-ionization (EoR). The origin of this `$\rm [C_{II}]$ deficit' is unclear. In this work, we study the $L_{\rm [C_{II}]}$-SFR relation of galaxies using a sample of $z=0-8$ galaxies with $M_*\approx10^7-5\times10^{11}\,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_{II}]}$/SFR of galaxies in terms of the following parameters: mass fraction of $\rm [C_{II}]$-emitting gas ($f_{\rm [C_{II}]}$), gas metallicity ($Z_{\rm gas}$), gas density ($n_{\rm gas}$) and gas depletion time ($t_{\rm dep}{}={}M_{\rm gas}{}/{}\rm SFR$). We find two distinct physical regimes, where $t_{\rm dep}$ ($Z_{\rm gas}$) is the main driver of the $\rm [C_{II}]$ deficit in $\rm H_2$-rich ($\rm H_2$-poor) galaxies. The observed $\rm [C_{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 $\rm [C_{II}]$ deficit is a common phenomenon of galaxies, and caution needs to be taken when applying a constant $L_{\rm [C_{II}]}$-to-SFR conversion factor derived from local SFGs to estimate cosmic SFR density at high redshifts and interpret data from upcoming $\rm [C_{II}]$ line intensity mapping experiments., Comment: 44 page, 24 figures, 8 tables and 8 appendices, accepted for publication in MNRAS. Key figures: fig. 16 & 17. Comments are welcome

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 Mstar > 10^8 Msun at z=0) at a resolution (~20 pc, m_b ~ 6x10^4 Msun) 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 (~10^6) and the inclusion of multi-channel stellar feedback. Here, we focus on validating the simulation predictions by comparing to observational data. We find that simulated galaxies with Mstar < 10^{10.5-11} Msun have star formation rates, gas masses, and metallicities in broad agreement with observations. These galaxy scaling relations extend to low masses (Mstar ~ 10^7 Msun) 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 $\Lambda$CDM Universe while also highlighting modeling challenges to be addressed in next-generation galaxy simulations., Comment: 31 pages, 17 figures, 3 tables, accepted for publication in MNRAS

Galaxy sizes correlate closely with the sizes of their parent dark matter haloes, suggesting a link between halo formation and galaxy growth. However, the precise nature of this relation and its scatter remains to be understood fully, especially for low-mass galaxies. We analyse the galaxy-halo size relation for low-mass ($M_\star \sim 10^{7-9} {\rm M_\odot}$) central galaxies over the past 12.5 billion years with the help of cosmological volume simulations (FIREbox) from the Feedback in Realistic Environments (FIRE) project. We find a nearly linear relationship between the half-stellar mass galaxy size $R_{1/2}$ and the parent dark matter halo virial radius $R_{\rm vir}$. This relation evolves only weakly since redshift $z = 5$: $R_{1/2} {\rm kpc} = (0.053\pm0.002)(R_{\rm vir}/35 {\rm kpc})^{0.934\pm0.054}$, with a nearly constant scatter $\langle \sigma \rangle = 0.084 [{\rm dex}]$. Whilst this ratio is similar to what is expected from models where galaxy disc sizes are set by halo angular momentum, the low-mass galaxies in our sample are not angular momentum supported, with stellar rotational to circular velocity ratios $v_{\rm rot} / v_{\rm circ} \sim 0.15$. Introducing redshift as another parameter to the GHSR does not decrease the scatter. Furthermore, this scatter does not correlate with any of the halo properties we investigate -- including spin and concentration -- suggesting that baryonic processes and feedback physics are instead critical in setting the scatter in the galaxy-halo size relation. Given the relatively small scatter and the weak dependence of the galaxy-halo size relation on redshift and halo properties for these low-mass central galaxies, we propose using galaxy sizes as an independent method from stellar masses to infer halo masses., Comment: 14 pages, 9 figures, 3 tables + 3 appendices (5 pages, 5 figures, 2 tables). Accepted for publication in MNRAS. Reuploaded after proofing stage

We study the interstellar medium in a sample of 27 high-redshift quasar host galaxies at z>6, using the [CII] 158um emission line and the underlying dust continuum observed at ~1kpc resolution with ALMA. By performing uv-plane spectral stacking of both the high and low spatial resolution data, we investigate the spatial and velocity extent of gas, and the size of the dust-emitting regions. We find that the average surface brightness profile of both the [CII] and the dust continuum emission can be described by a steep component within a radius of 2kpc, and a shallower component with a scale length of 2kpc, detected up to ~10kpc. The surface brightness of the extended emission drops below ~1% of the peak at radius of ~5kpc, beyond which it constitutes 10-20% of the total measured flux density. Although the central component of the dust continuum emission is more compact than that of the [CII] emission, the extended components have equivalent profiles. The observed extended components are consistent with those predicted by hydrodynamical simulations of galaxies with similar infrared luminosities, where the dust emission is powered by star formation. The [CII] spectrum measured in the mean uv-plane stacked data can be described by a single Gaussian, with no observable [CII] broad-line emission (velocities in excess of >500km/s), that would be indicative of outflows. Our findings suggest that we are probing the interstellar medium and associated star formation in the quasar host galaxies up to radii of 10kpc, whereas we find no evidence for halos or outflows., Comment: 32 pages, 20 figures, 1 table, accepted for publication in ApJ

The relation between infrared excess (IRX) and UV spectral slope ($\beta_{\rm 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-$\beta_{\rm UV}$ relation with a sample of $z=2-6$ galaxies ($M_*\approx 10^9-10^{12}\,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-$\beta_{\rm 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 $\beta_{\rm UV}$ is set both by stellar properties, UV optical depth, and the dust extinction law. Overall, much of the scatter in the IRX-$\beta_{\rm UV}$ relation of our sample is found to be driven by variations of the intrinsic UV spectral slope. We further assess how the IRX-$\beta_{\rm 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., Comment: 32 pages, 28 figures, 3 tables, submitted to MNRAS (comments are welcomed)

Wind-driven outflows are observed around a broad range of accreting objects throughout the Universe, ranging from forming low-mass stars to super-massive black holes. We study the interaction between a central isotropic wind and an infalling, rotating, envelope, determining the steady-state cavity shape formed at their interface under the assumption of weak mixing. The shape of the resulting wind-blown cavity is elongated and self-similar, with a physical size determined by the ratio between wind ram pressure and envelope thermal pressure. We compute the growth of a warm turbulent mixing-layer between the shocked wind and the deflected envelope, and calculate the resultant broad line profile, under the assumption of a linear (Couette-type) velocity profile across the layer. We then test our model against the warm broad velocity component observed in CO $J$=16--15 by Herschel/HIFI in the protostar Serpens-Main SMM1. Given independent observational constraints on the temperature and density of the dust envelope around SMM1, we find an excellent match to all its observed properties (line profile, momentum, temperature) and to the SMM1 outflow cavity width for a physically reasonable set of parameters: a ratio of wind to infall mass-flux $\simeq 4\%$, a wind speed $v_{\rm w} \simeq 30$ km/s, an interstellar abundance of CO and H$_2$, and a turbulent entrainment efficiency consistent with laboratory experiments. The inferred ratio of ejection to disk accretion rate, $\simeq 6-20\%$, is in agreement with current disk wind theories. Thus, the model provides a new framework to reconcile the modest outflow cavity widths in protostars with the large observed flow velocities. Being self-similar, it is applicable over a broader range of astrophysical contexts as well., Comment: 31 pages, 21 figures, accepted to ApJ for publication (comments are welcome)

Dust temperature is an important property of the interstellar medium (ISM) of galaxies. It is required when converting (sub)millimeter broadband flux to total infrared luminosity (L_IR), and hence star formation rate, in high-z 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* > 10^10 Msun) 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 does not strongly evolve with redshift over z=2-6 at fixed IR luminosity but is tightly correlated with L_IR at fixed z. The mass-weighted temperature is important for accurately estimating the total dust mass. We also analyse an 'equivalent' dust temperature for converting (sub)millimeter 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., Comment: 26 pages, 15 figures, accepted for publication in MNRAS

We present a suite of 34 high-resolution cosmological zoom-in simulations consisting of thousands of halos up to M_halo~10^12 M_sun (M_star~10^10.5 M_sun) 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 extinction, 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 (beta_UV), despite the patchy, clumpy dust geometry shown in our simulations. We find that the IRX-beta_UV relation is mainly determined by the shape of the extinction curve 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 (LFs) 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 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., Comment: 21 pages, 17 figures, key results shown in Figs. 5, 8, 10, 11, MNRAS accepted. Data products including mock galaxy SEDs and images are available upon request

Recent long wavelength observations on the thermal dust continuum suggest that the Rayleigh-Jeans (RJ) tail can be used as a time-efficient quantitative probe of the dust and ISM mass in high-$z$ galaxies. We use high-resolution cosmological simulations from the Feedback in Realistic Environment (FIRE) project to analyze the dust emission of $M_*>10^{10}\;M_{\odot}$ galaxies at $z=2-4$. Our simulations (MassiveFIRE) explicitly include various forms of stellar feedback, and they produce the stellar masses and star formation rates of high-$z$ galaxies in agreement with observations. Using radiative transfer modelling, we show that sub-millimeter (sub-mm) luminosity and molecular ISM mass are tightly correlated and that the overall normalization is in quantitative agreement with observations. Notably, sub-mm luminosity traces molecular ISM mass even during starburst episodes as dust mass and mass-weighted temperature evolve only moderately between $z=4$ and $z=2$, including during starbursts. Our finding supports the empirical approach of using broadband sub-mm flux as a proxy for molecular gas content in high-$z$ galaxies. We thus expect single-band sub-mm observations with ALMA to dramatically increase the sample size of high-$z$ galaxies with reliable ISM masses in the near future., Comment: 5 pages, 3 figures, accepted by MNRAS

The thermal and chemical properties of the hot diffuse intragroup medium (IGrM) provide important constraints on the feedback processes associated with massive galaxy formation and evolution. Here we explore these constraints via a detailed analysis of the global properties of simulated z<3 galaxy groups from a cosmological simulation that includes a well-constrained prescription for stellar/supernovae-powered galactic outflows but no AGN feedback. Our aims are to (a) establish a baseline against which we will compare future models; (b) identify model successes due to stellar/supernovae-powered outflows; and (c) pinpoint features that signal the need for, and constrain the nature of, AGN feedback. Our simulation successfully reproduces key observed z=0 group IGrM properties, including the various X-ray Lx - Tx - entropy scaling relations, for all but the most massive groups. The z<1 redshift evolution of these also agree with observations. Contrary to expectations, the simulated groups' IGrM does not suffer catastrophic cooling. Yet, the z=0 group stellar mass is ~ 2X too large. This is due to the build-up of cold gas in the massive galaxies before they are incorporated inside groups. This in turn indicates that other feedback mechanisms must activate in real galaxies once their stellar masses grow to a few X 10^{10} M_sun. We show that these must be powerful enough to expel a significant fraction of the gas from the galactic halos. Gentle maintenance-mode (quenching) AGN feedback, as seen in galaxy clusters, will not do. Just as importantly, we find that the stellar/supernovae-powered winds are essential for understanding the IGrM metal abundances. Our simulation is able to reproduce the observed relationship between the global IGrM iron and silicon abundance and the group X-ray temperature, and these results ought to be relatively insensitive to the addition of AGN feedback., Comment: 26 pages, 16 figures, accepted by MNRAS