Kinetic or thermal AGN feedback in simulations of isolated and merging disc galaxies calibrated by the M-  relation (original) (raw)

Large-eddy simulations of isolated disc galaxies with thermal and turbulent feedback

We present a subgrid-scale model for the Multi-phase Interstellar medium, Star formation, and Turbulence (MIST) and explore its behaviour in high-resolution large-eddy simulations of isolated disc galaxies. MIST follows the evolution of a clumpy cold and a diffuse warm component of the gas within a volume element which exchange mass and energy via various cooling, heating and mixing processes. The star formation rate is dynamically computed from the state of the gas in the cold phase. An important feature of MIST is the treatment of unresolved turbulence in the two phases and its interaction with star formation and feedback by supernovae. This makes MIST a particularly suitable model for the interstellar medium in galaxy simulations. We carried out a suite of simulations varying fundamental parameters of our feedback implementation. Several observational properties of galactic star formation are reproduced in our simulations, such as an average star formation efficiency ∼1 per cent, a typical velocity dispersion around ∼10 km s −1 in star-forming regions, and an almost linear relationship between the column densities of star formation and dense molecular gas.

Cold Gas in Massive Galaxies as A Critical Test of Black Hole Feedback Models

2021

Black hole feedback has been widely implemented as the key recipe to quench star formation in massive galaxies in modern semi-analytic models and hydrodynamical simulations. As the theoretical details surrounding the accretion and feedback of black holes continue to be refined, various feedback models have been implemented across simulations, with notable differences in their outcomes. Yet, most of these simulations have successfully reproduced some observations, such as stellar mass function and star formation rate density in the local Universe. We use the recent observation on the change of neutral hydrogen gas mass (including both H2 and H1) with star formation rate of massive central disc galaxies as a critical test of black hole feedback models across several simulations. We find that IllustrisTNG reproduces the observations significantly better than the other models that we have tested. This favors IllustrisTNG's treatment of active galactic nuclei---where kinetic winds ar...

Wind-driven gas networks and star formation in galaxies: reaction-advection hydrodynamic simulations

Monthly Notices of the Royal Astronomical Society, 2001

The effects of wind-driven star formation feedback on the spatio-temporal organization of stars and gas in galaxies is studied using two-dimensional intermediaterepresentational quasi-hydrodynamical simulations. The model retains only a reduced subset of the physics, including mass and momentum conservation, fully nonlinear fluid advection, inelastic macroscopic interactions, threshold star formation, and momentum forcing by winds from young star clusters on the surrounding gas. Expanding shells of swept-up gas evolve through the action of fluid advection to form a "turbulent" network of interacting shell fragments whose overall appearance is a web of filaments (in two dimensions). A new star cluster is formed whenever the column density through a filament exceeds a critical threshold based on the gravitational instability criterion for an expanding shell, which then generates a new expanding shell after some time delay. A filament-finding algorithm is developed to locate the potential sites of new star formation.

Feedback and recycled wind accretion: assembling the z= 0 galaxy mass function

Monthly Notices of the Royal Astronomical Society, 2010

We analyse cosmological hydrodynamic simulations that include theoretically and observationally motivated prescriptions for galactic outflows. If these simulated winds accurately represent winds in the real Universe, then material previously ejected in winds provides the dominant source of gas infall for new star formation at redshifts z < 1. This recycled wind accretion, or wind mode, provides a third physically distinct accretion channel in addition to the 'hot' and 'cold' modes emphasized in recent theoretical studies. The recycling time of wind material (t rec) is shorter in higher mass systems owing to the interaction between outflows and the increasingly higher gas densities in and around higher mass haloes. This differential recycling plays a central role in shaping the present-day galaxy stellar mass function (GSMF), because declining t rec leads to increasing wind mode galaxy growth in more massive haloes. For the three feedback models explored, the wind mode dominates above a threshold mass that primarily depends on wind velocity; the shape of the GSMF therefore can be directly traced back to the feedback prescription used. If we remove all particles that were ever ejected in a wind, then the predicted GSMFs are much steeper than observed. In this case, galaxy masses are suppressed both by the ejection of gas from galaxies and by the hydrodynamic heating of their surroundings, which reduces subsequent infall. With wind recycling included, the simulation that incorporates our favoured momentum-driven wind scalings reproduces the observed GSMF for stellar masses 10 9 M ≤ M ≤ 5 × 10 10 M. At higher masses, wind recycling leads to excessive galaxy masses and star formation rates relative to observations. In these massive systems, some quenching mechanism must suppress not only the direct accretion from the primordial intergalactic medium but the re-accretion of gas ejected from star-forming galaxies. In short, as has long been anticipated, the form of the GSMF is governed by outflows; the unexpected twist here for our simulated winds is that it is not primarily the ejection of material but how the ejected material is re-accreted that governs the GSMF.

Large-eddy simulations of isolated disk galaxies with thermal and turbulent feedback

Monthly Notices of the Royal Astronomical Society

We present a subgrid-scale model for the Multi-phase Interstellar medium, Star formation, and Turbulence (MIST) and explore its behavior in high-resolution large-eddy simulations of isolated disk galaxies. MIST follows the evolution of a clumpy cold and a diffuse warm component of the gas within a volume element which exchange mass and energy via various cooling, heating and mixing processes. The star formation rate is dynamically computed from the state of the gas in the cold phase. An important feature of MIST is the treatment of unresolved turbulence in the two phases and its interaction with star formation and feedback by supernovae. This makes MIST a particularly suitable model for the interstellar medium in galaxy simulations. We carried out a suite of simulations varying fundamental parameters of our feedback implementation. Several observational properties of galactic star formation are reproduced in our simulations, such as an average star formation efficiency ~1%, a typica...

A thermal-kinetic subgrid model for supernova feedback in simulations of galaxy formation

Cornell University - arXiv, 2022

We present a subgrid model for supernova feedback designed for simulations of galaxy formation. The model uses thermal and kinetic channels of energy injection, which are built upon the stochastic kinetic and thermal models for stellar feedback used in the and simulations, respectively. In the thermal channel, the energy is distributed statistically isotropically and injected stochastically in large amounts per event, which minimizes spurious radiative energy losses. In the kinetic channel, we inject the energy in small portions by kicking gas particles in pairs in opposite directions. The implementation of kinetic feedback is designed to conserve energy, linear momentum and angular momentum, and is statistically isotropic. To test and validate the model, we run simulations of isolated Milky Way-mass and dwarf galaxies, in which the gas is allowed to cool down to 10 K. Using the thermal and kinetic channels together, we obtain smooth star formation histories and powerful galactic winds with realistic mass loading factors. Furthermore, the model produces spatially resolved star formation rates and velocity dispersions that are in agreement with observations. We vary the numerical resolution by several orders of magnitude and find excellent convergence of the global star formation rates and the mass loading of galactic winds. We show that large thermal-energy injections generate a hot phase of the interstellar medium (ISM) and modulate the star formation by ejecting gas from the disc, while the low-energy kicks increase the turbulent velocity dispersion in the neutral ISM, which in turn helps suppress star formation.

Cosmological Simulations of the Intergalactic Medium Evolution. II. Galaxy Model and Feedback

The Astrophysical Journal, 2015

We investigate models of self-consistent chemical enrichment of the intergalactic medium (IGM) from z = 6.0 → 1.5, based on hydrodynamic simulations of structure formation that explicitly incorporate outflows from star-forming galaxies. Our main result is that outflow parametrizations derived from observations of local starburst galaxies, in particular momentum-driven wind scenarios, provide the best agreement with observations of C IV absorption at z ∼ 2-5. Such models sufficiently enrich the high-z IGM to produce a global mass density of C IV absorbers that is relatively invariant from z = 5.5 → 1.5, in agreement with observations. This occurs despite continual IGM enrichment causing an increase in volumeaveraged metallicity by ∼ × 5-10 over this redshift range, because energy input accompanying the enriching outflows causes a drop in the global ionization fraction of C IV. Comparisons to observed C IV column density and linewidth distributions and C IV-based pixel optical depth ratios provide significant constraints on wind models. Our best-fitting outflow models show mean IGM temperatures only slightly above our no-outflow case, metal filling factors of just a few per cent with volume-weighted metallicities around 10 −3 at z ∼ 3, significant amounts of collisionally ionized C IV absorption and a metallicity-density relationship that rises rapidly at low overdensities and flattens at higher ones. In general, we find that outflow speeds must be high enough to enrich the low-density IGM at early times but low enough not to overheat it, and concurrently must significantly suppress early star formation while still producing enough early metals. It is therefore non-trivial that locally calibrated momentum-driven wind scenarios naturally yield the desired strength and evolution of outflows, and suggest that such models represent a significant step towards understanding the impact of galactic outflows on galaxies and the IGM across cosmic time.

A warm mode of gas accretion on forming galaxies

2012

We present results from high--resolution cosmological hydrodynamical simulations of a Milky--Way-sized halo, aimed at studying the effect of feedback on the nature of gas accretion. Simulations include a model of inter-stellar medium and star formation, in which SN explosions provide effective thermal feedback. We distinguish between gas accretion onto the halo, which occurs when gas particles cross the halo virial radius, and gas accretion onto the central galaxy, which takes place when gas particles cross the inner one-tenth of the virial radius. Gas particles can be accreted through three different channels, depending on the maximum temperature value, T_ max, reached during the particles' past evolution: a cold channel for T_ max<2.5 × 10^5 K, a hot one for T>10^6K, and a warm one for intermediate values of T_ max. We find that the warm channel is at least as important as the cold one for gas accretion onto the central galaxy. This result is at variance with previous fi...

Generating Hot Gas in Simulations of Disk-Galaxy Major Mergers

Astrophysical Journal, 2004

We report on the merger-induced generation of a shock-heated gas wind and formation of a remnant gas halo in simulations of colliding disk galaxies. The simulations use cosmologically motivated initial conditions and include the effects of radiative cooling, star formation, stellar feedback and the non-adiabatic heating of gas. The non-adiabatic heating, i.e. shocks, generated in the final merger forces gas out of the central region of the merger remnant and into the dark-matter halo. We demonstrate that the amount of heating depends on the size of the progenitor disk galaxy as well as the initial orbit the galaxies are placed on. Based upon these dependencies, we motivate a possible recipe for including this effect in semi-analytic models of galaxy formation.

The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations

Monthly Notices of the Royal Astronomical Society, 2015

We present results from 13 cosmological simulations that explore the parameter space of the 'Evolution and Assembly of GaLaxies and their Environments' (EAGLE) simulation project. Four of the simulations follow the evolution of a periodic cube L = 50 cMpc on a side, and each employs a different subgrid model of the energetic feedback associated with star formation. The relevant parameters were adjusted so that the simulations each reproduce the observed galaxy stellar mass function at z = 0.1. Three of the simulations fail to form disc galaxies as extended as observed, and we show analytically that this is a consequence of numerical radiative losses that reduce the efficiency of stellar feedback in high-density gas. Such losses are greatly reduced in the fourth simulation-the EAGLE reference model-by injecting more energy in higher density gas. This model produces galaxies with the observed size distribution, and also reproduces many galaxy scaling relations. In the remaining nine simulations, a single parameter or process of the reference model was varied at a time. We find that the properties of galaxies with stellar mass M (the 'knee' of the galaxy stellar mass function) are largely governed by feedback associated with star formation, while those of more massive galaxies are also controlled by feedback from accretion on to their central black holes. Both processes must be efficient in order to reproduce the observed galaxy population. In general, simulations that have been calibrated to reproduce the low-redshift galaxy stellar mass function will still not form realistic galaxies, but the additional requirement that galaxy sizes be acceptable leads to agreement with a large range of observables.