Search for: All records

The early growth of black holes (BHs) in high-redshift galaxies is likely feedback regulated. While radiative feedback has been extensively studied, the role of mechanical feedback has received less scrutiny to date. Here, we use high-resolution parsec-scale hydrodynamical simulations to study jet propagation and its effect on 100 M⊙ BH accretion in the dense, low-metallicity gas expected in early protogalaxies. As the jet propagates, it shocks the surrounding gas forming a jet cocoon. The cocoon consists of a rapidly cooling cold phase at the interface with the background gas and an overpressured subsonic phase of reverse shock-heated gas filling the interior. We vary the background gas density and temperature, BH feedback efficiency, and the jet model. We found that the width of the jet cocoon roughly follows a scaling derived by assuming momentum conservation in the jet-propagation direction and energy conservation in the lateral directions. Depending on the assumed gas and jet properties, the cocoon either stays elongated to large radii or isotropizes before reaching the Bondi radius, forming a nearly spherical bubble. Lower jet velocities and higher background gas densities result in self-regulation to higher momentum fluxes and elongated cocoons. In all cases, the outward cocoon momentum flux balancesmore »the inward inflowing gas momentum flux near the Bondi radius, which ultimately regulates BH accretion. The time-averaged accretion rate always remains below the Bondi rate, and exceeds the Eddington rate only if the ambient medium is dense and cold, and/or the jet is weak (low velocity and mass loading).

« less

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

« less

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 furthermore »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.

« less

Interstellar chemistry is important for galaxy formation, as it determines the rate at which gas can cool, and enables us to make predictions for observable spectroscopic lines from ions and molecules. We explore two central aspects of modelling the chemistry of the interstellar medium (ISM): (1) the effects of local stellar radiation, which ionizes and heats the gas, and (2) the depletion of metals on to dust grains, which reduces the abundance of metals in the gas phase. We run high-resolution (400 M⊙ per baryonic particle) simulations of isolated disc galaxies, from dwarfs to Milky Way-mass, using the fire galaxy formation models together with the chimes non-equilibrium chemistry and cooling module. In our fiducial model, we couple the chemistry to the stellar fluxes calculated from star particles using an approximate radiative transfer scheme; and we implement an empirical density-dependent prescription for metal depletion. For comparison, we also run simulations with a spatially uniform radiation field, and without metal depletion. Our fiducial model broadly reproduces observed trends in H i and H2 mass with stellar mass, and in line luminosity versus star formation rate for [C ii]$_{158 \rm {\mu m}}$, [O i]$_{63 \rm {\mu m}}$, [O iii]$_{88 \rm {\mu m}}$, [N ii]$_{122 \rm {\mu m}}$, andmore »H α6563Å. Our simulations with a uniform radiation field predict fainter luminosities, by up to an order of magnitude for [O iii]$_{88 \rm {\mu m}}$ and H α6563Å, while ignoring metal depletion increases the luminosity of carbon and oxygen lines by a factor ≈ 2. However, the overall evolution of the galaxy is not strongly affected by local stellar fluxes or metal depletion, except in dwarf galaxies where the inclusion of local fluxes leads to weaker outflows and hence higher gas fractions.

« less

ABSTRACT 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’more »physics, including spectrally resolved cosmic rays and supermassive black hole growth and feedback.« less

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

« less

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

« less

Galaxy mergers are crucial to understanding galaxy evolution, therefore we must determine their observational signatures to select them from large IFU galaxy samples such as MUSE and SAMI. We employ 24 high-resolution idealized hydrodynamical galaxy merger simulations based on the ‘Feedback In Realistic Environment’ (FIRE-2) model to determine the observability of mergers to various configurations and stages using synthetic images and velocity maps. Our mergers cover a range of orbital configurations at fixed 1:2.5 stellar mass ratio for two gas rich spirals at low redshift. Morphological and kinematic asymmetries are computed for synthetic images and velocity maps spanning each interaction. We divide the interaction sequence into three: (1) the pair phase; (2) the merging phase; and (3) the post-coalescence phase. We correctly identify mergers between first pericentre passage and 500 Myr after coalescence using kinematic asymmetry with 66 per cent completeness, depending upon merger phase and the field of view of the observation. We detect fewer mergers in the pair phase (40 per cent) and many more in the merging and post-coalescence phases (97 per cent). We find that merger detectability decreases with field of view, except in retrograde mergers, where centrally concentrated asymmetric kinematic features enhances their detectability. Using a cut-off derived from a combinationmore »of photometric and kinematic asymmetry, we increase these detections to 89 per cent overall, 79 per cent in pairs, and close to 100 per cent in the merging and post-coalescent phases. By using this combined asymmetry cut-off we mitigate some of the effects caused by smaller fields of view subtended by massively multiplexed integral field spectroscopy programmes.

« less

We compare an analytic model for the evolution of supernova-driven superbubbles with observations of local and high-redshift galaxies, and the properties of intact Hishells in local star-forming galaxies. Our model correctly predicts the presence of superwinds in local star-forming galaxies (e.g., NGC 253) and the ubiquity of outflows nearz∼ 2. We find that high-redshift galaxies may “capture” 20%–50% of their feedback momentum in the dense ISM (with the remainder escaping into the nearby CGM), whereas local galaxies may contain ≲10% of their feedback momentum from the central starburst. Using azimuthally averaged galaxy properties, we predict that most superbubbles stall and fragmentwithinthe ISM, and that this occurs at, or near, the gas scale height. We find a consistent interpretation in the observed Hibubble radii and velocities, and predict that most will fragment within the ISM, and that those able to break out originate from short dynamical time regions (where the dynamical time is shorter than feedback timescales). Additionally, we demonstrate that models with constant star cluster formation efficiency per Toomre mass are inconsistent with the occurrence of outflows from high-zstarbursts and local circumnuclear regions.

With Σ_SFR∼ 4200M_⊙yr⁻¹kpc⁻², SPT 0346–52 (z= 5.7) is the most intensely star-forming galaxy discovered by the South Pole Telescope. In this paper, we expand on previous spatially resolved studies, using ALMA observations of dust continuum, [Nii] 205μm, [Cii] 158μm, [Oi] 146μm, and undetected [Nii] 122μm and [Oi] 63μm emission to study the multiphase interstellar medium (ISM) in SPT 0346–52. We use pixelated, visibility-based lens modeling to reconstruct the source-plane emission. We also model the source-plane emission using the photoionization codecloudyand find a supersolar metallicity system. We calculateT_dust= 48.3 K andλ_peak= 80μm and see line deficits in all five lines. The ionized gas is less dense than comparable galaxies, withn_e< 32 cm⁻³, while ∼20% of the [Cii] 158μm emission originates from the ionized phase of the ISM. We also calculate the masses of several phases of the ISM. We find that molecular gas dominates the mass of the ISM in SPT 0346–52, with the molecular gas mass ∼4× higher than the neutral atomic gas mass and ∼100× higher than the ionized gas mass.