Category: Uncategorized

Su et al, available on arXiv

Abstract: Using high-resolution simulations from the FIRE-2 (Feedback In Realistic Environments) project, we study the effects of discreteness in stellar feedback processes on the evolution of galaxies and the properties of the interstellar medium (ISM). We specifically consider the discretization of supernovae (SNe), including hypernovae (HNe), and sampling the initial mass function (IMF). We study these processes in cosmological simulations of dwarf galaxies with z=0 stellar masses Mstar~10^4-3×10^6 Msun (halo masses ~ 10^9-10^10 Msun). We show that the discrete nature of individual SNe (as opposed to a model in which their energy/momentum deposition is continuous over time, similar to stellar winds) is crucial in generating a reasonable ISM structure and galactic winds and in regulating dwarf stellar masses. However, once SNe are discretized, accounting for the effects of IMF sampling on continuous mechanisms such as radiative feedback and stellar mass-loss (as opposed to adopting IMF-averaged rates) has weak effects on galaxy-scale properties. We also consider the effects of rare HNe events with energies ~10^53 erg. The effects of HNe are similar to the effects of clustered explosions of SNe — which are already captured in our default simulation setup — and do not quench star formation (provided that the HNe do not dominate the total SNe energy budget), which suggests that HNe yield products should be observable in ultra-faint dwarfs today.

Chan et al., available on arXiv

Abstract: We test if the cosmological zoom-in simulations of isolated galaxies from the FIRE project reproduce the properties of ultra diffuse galaxies. We show that stellar feedback-generated outflows that dynamically heat galactic stars, together with a passively aging stellar population after imposed quenching (from e.g. infall into a galaxy cluster), naturally reproduce the observed population of red UDGs, without the need for high spin halos or dynamical influence from their host cluster. We reproduce the range of surface brightness, radius and absolute magnitude of the observed z=0 red UDGs by quenching simulated galaxies at a range of different times. They represent a mostly uniform population of dark matter-dominated galaxies with M_star ~1e8 Msun, low metallicity and a broad range of ages. The most massive simulated UDGs require earliest quenching and are therefore the oldest. Our simulations provide a good match to the central enclosed masses and the velocity dispersions of the observed UDGs (20-50 km/s). The enclosed masses of the simulated UDGs remain largely fixed across a broad range of quenching times because the central regions of their dark matter halos complete their growth early. A typical UDG forms in a dwarf halo mass range of Mh~4e10-1e11 Msun. The most massive red UDG in our sample requires quenching at z~3 when its halo reached Mh ~ 1e11 Msun. If it, instead, continues growing in the field, by z=0 its halo mass reaches > 5e11 Msun, comparable to the halo of an L* galaxy. If our simulated dwarfs are not quenched, they evolve into bluer low-surface brightness galaxies with mass-to-light ratios similar to observed field dwarfs. While our simulation sample covers a limited range of formation histories and halo masses, we predict that UDG is a common, and perhaps even dominant, galaxy type around Ms~1e8 Msun, both in the field and in clusters.

Escala et al., available at

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

Su et al., available on arXiv

Abstract: Using high-resolution magnetohydrodynamic simulations of idealized, non-cosmological galaxies, we investigate how cooling, star formation, and stellar feedback affect galactic magnetic fields. We find that the amplification histories, saturation values, and morphologies of the magnetic fields vary considerably depending on the baryonic physics employed, primarily because of differences in the gas density distribution. In particular, adiabatic runs and runs with a sub-grid (effective equation of state) stellar feedback model yield lower saturation values and morphologies that exhibit greater large-scale order compared with runs that adopt explicit stellar feedback and runs with cooling and star formation but no feedback. The discrepancies mostly lie in gas denser than the galactic average, which requires cooling and explicit fragmentation to capture. Independent of the baryonic physics included, the magnetic field strength scales with gas density as B~n^(2/3), suggesting isotropic flux freezing or equipartition between the magnetic and gravitational energies during the field amplification. We conclude that accurate treatments of cooling, star formation, and stellar feedback are crucial for obtaining the correct magnetic field strength and morphology in dense gas, which, in turn, is essential for properly modeling other physical processes that depend on the magnetic field, such as cosmic ray feedback.

Ma et al., available on arXiv

Abstract: We study the morphologies and sizes of galaxies at z>5 using high-resolution cosmological zoom-in simulations from the Feedback In Realistic Environments project. The galaxies show a variety of morphologies, from compact to clumpy to irregular. The simulated galaxies have more extended morphologies and larger sizes when measured using rest-frame optical B-band light than rest-frame UV light; sizes measured from stellar mass surface density are even larger. The UV morphologies are usually dominated by several small, bright young stellar clumps that are not always associated with significant stellar mass. The B-band light traces stellar mass better than the UV, but it can also be biased by the bright clumps. At all redshifts, galaxy size correlates with stellar mass/luminosity with large scatter. The half-light radii range from 0.01 to 0.2 arcsec (0.05-1 kpc physical) at fixed magnitude. At z>5, the size of galaxies at fixed stellar mass/luminosity evolves as (1+z)^{-m}, with m~1-2. For galaxies less massive than M_star~10^8 M_sun, the ratio of the half-mass radius to the halo virial radius is ~10% and does not evolve significantly at z=5-10; this ratio is typically 1-5% for more massive galaxies. A galaxy’s “observed” size decreases dramatically at shallower surface brightness limits. This effect may account for the extremely small sizes of z>5 galaxies measured in the Hubble Frontier Fields. We provide predictions for the cumulative light distribution as a function of surface brightness for typical galaxies at z=6.

Orr et al., available on arXiv

Abstract: In a recent work based on 3200 stacked Hα maps of galaxies at z~1, Nelson et al. find evidence for “coherent star formation”: the stacked star formation rate (SFR) profiles of galaxies above (below) the “star formation main sequence” (MS) are above (below) that of galaxies on the MS at all radii. One might interpret this result as inconsistent with highly bursty star formation and evidence that galaxies evolve smoothly along the MS rather than crossing it many times. We analyze six simulated galaxies at z~1 from the Feedback in Realistic Environments (FIRE) project in a manner analogous to the observations to test whether the above interpretations are correct. The trends in stacked SFR profiles are qualitatively consistent with those observed. However, SFR profiles of individual galaxies are much more complex than the stacked profiles: the former can be flat or even peak at large radii because of the highly clustered nature of star formation in the simulations. Moreover, the SFR profiles of individual galaxies above (below) the MS are not systematically above (below) those of MS galaxies at all radii. We conclude that the time-averaged coherent star formation evident stacks of observed galaxies is consistent with highly bursty, clumpy star formation of individual galaxies and is not evidence that galaxies evolve smoothly along the MS.

Hopkins et al., available on arXiv.

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

Angles-Alcazar et al., available on arXiv.

Abstract: We introduce massive black holes (BHs) in the Feedback In Realistic Environments project and perform high-resolution cosmological hydrodynamic simulations of quasar-mass halos (M_halo(z=2)~10^12.5 Msun) down to z=1. These simulations model stellar feedback by supernovae, stellar winds, and radiation, and BH growth using a gravitational torque-based prescription tied to resolved properties of galactic nuclei. We do not include BH feedback. We show that early BH growth occurs through short (<~1 Myr) accretion episodes that can reach or even exceed the Eddington rate. In this regime, BH growth is limited by bursty stellar feedback continuously evacuating gas from galactic nuclei, and BHs remain under-massive relative to the local M_BH-M_bulge relation. BH growth is more efficient at later times, when the nuclear stellar potential retains a significant gas reservoir, star formation becomes less bursty, and galaxies settle into a more ordered state, with BHs rapidly converging onto the scaling relation when the host reaches M_bulge~10^10 Msun. Our results are not sensitive to the details of the accretion model so long as BH growth is tied to the gas content within ~100 pc of the BH. Our simulations imply that bursty stellar feedback has strong implications for BH and AGN demographics, especially in the early Universe and for low-mass galaxies.

Hopkins et al., available on arXiv.

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

Price et al., available on arXiv.

Abstract: Accurate measurements of galaxy masses and sizes are key to tracing galaxy evolution over time. Cosmological zoom-in simulations provide an ideal test bed for assessing the recovery of galaxy properties from observations. Here, we utilize galaxies with Mstar ~ 10^10-10^11.5 Msun at z ~ 1.7-2 from the MassiveFIRE cosmological simulation suite, part of the Feedback in Realistic Environments (FIRE) project. Using mock multi-band images, we compare intrinsic galaxy masses and sizes to observational estimates. We find that observations accurately recover stellar masses, with a slight average underestimate of ~0.06 dex and a ~0.15 dex scatter. Recovered half-light radii agree well with intrinsic half-mass radii when averaged over all viewing angles, with a systematic offset of ~0.1 dex (with the half-light radii being larger) and a scatter of ~0.2 dex. When using color gradients to account for mass-to-light variations, recovered half-mass radii also exceed the intrinsic half-mass radii by ~0.1 dex. However, if not properly accounted for, aperture effects can bias size estimates by ~0.1 dex. No differences are found between the mass and size offsets for star-forming and quiescent galaxies. Variations in viewing angle are responsible for ~25% of the scatter in the recovered masses and sizes. Our results thus suggest that the intrinsic scatter in the mass–size relation may have previously been overestimated by ~25%. Moreover, orientation-driven scatter causes the number density of very massive galaxies to be overestimated by ~0.5 dex at Mstar ~ 10^11.5 Msun.