Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userH2 chemistry in galaxy simulations: an improved supernova feedback model
78
0
0.0
(
0
)
Authors
Alessandro Lupi
Ask ChatGPT about the research
No Arabic abstract
In this study, we present and validate a variation of recently-developed physically motivated sub-grid prescriptions for supernova feedback that account for the unresolved energy-conserving phase of the bubble expansion. Our model builds upon the implementation publicly available in the mesh-less hydrodynamic code GIZMO, and is coupled with the chemistry library KROME. Here, we test it against different setups, to address how it affects the formation/dissociation of molecular hydrogen (H$_2$). First, we explore very idealised conditions, to show that it can accurately reproduce the terminal momentum of the blast-wave independent of resolution. Then, we apply it to a suite of numerical simulations of an isolated Milky Way-like galaxy and compare it with a similar run employing the delayed-cooling sub-grid prescription. We find that the delayed-cooling model, by pressurising ad-hoc the gas, is more effective in suppressing star formation. However, to get this effect, it must maintain the gas warm/hot at densities where it is expected to cool efficiently, artificially changing the thermo-chemical state of the gas, and reducing the H$_2$ abundance even in dense gas. Mechanical feedback, on the other hand, is able to reproduce the H$_2$ column densities without altering the gas thermodynamics, and, at the same time, drives more powerful outflows. However, being less effective in suppressing star formation, it over-predicts the Kennicutt-Schmidt relation by a factor of about 2.5. Finally, we show that the model is consistent at different resolution levels, with only mild differences.
rate research
Read More
We study various implementations of supernova feedback model and present the results of our `Osaka feedback model using isolated galaxy simulations performed by the smoothed particle hydrodynamics (SPH) code {small GADGET-3}. Our model is a modified version of Stinson et al.s work, and we newly add the momentum kick for SN feedback rather than only thermal feedback. We incorporate the physical properties of SN remnants from the results of Chevalier and McKee & Ostriker, such as the effective radius of SN bubble and the remnant life-time, in the form of Sedov-Taylor (ST)-like solutions with the effect of radiative cooling. Our model utilizes the local, physical parameters such as density and temperature of the ISM rather than galactic or halo properties to determine the galactic wind velocity or mass-loading factor. The Osaka model succeeds in self-regulating star formation, and naturally produces galactic outflow with variable velocities depending on the local environment and available SN energy as a function of time.An important addition to our previous work by Aoyama et al. is the implementation of the {small CELib} chemistry library which allows us to deal with the time-dependent input of energy and metal yields for type Ia & II supernovae (SNe) and asymptotic giant branch (AGB) stars. As initial tests of our model, we apply it to isolated galaxy simulations, and examine various galactic properties and compare with observational data including metal abundances.
Cosmological simulations of galaxies have typically produced too many stars at early times. We study the global and morphological effects of radiation pressure (RP) in eight pairs of high-resolution cosmological galaxy formation simulations. We find that the additional feedback suppresses star formation globally by a factor of ~2. Despite this reduction, the simulations still overproduce stars by a factor of ~2 with respect to the predictions provided by abundance matching methods for halos more massive than 5E11 Msun/h (Behroozi, Wechsler & Conroy 2013). We also study the morphological impact of radiation pressure on our simulations. In simulations with RP the average number of low mass clumps falls dramatically. Only clumps with stellar masses Mclump/Mdisk <= 5% are impacted by the inclusion of RP, and RP and no-RP clump counts above this range are comparable. The inclusion of RP depresses the contrast ratios of clumps by factors of a few for clump masses less than 5% of the disk masses. For more massive clumps, the differences between and RP and no-RP simulations diminish. We note however, that the simulations analyzed have disk stellar masses below about 2E10 Msun/h. By creating mock Hubble Space Telescope observations we find that the number of clumps is slightly reduced in simulations with RP. However, since massive clumps survive the inclusion of RP and are found in our mock observations, we do not find a disagreement between simulations of our clumpy galaxies and observations of clumpy galaxies. We demonstrate that clumps found in any single gas, stellar, or mock observation image are not necessarily clumps found in another map, and that there are few clumps common to multiple maps.
In this study, we present a suite of high-resolution numerical simulations of an isolated galaxy to test a sub-grid framework to consistently follow the formation and dissociation of H$_2$ with non-equilibrium chemistry. The latter is solved via the package KROME, coupled to the mesh-less hydrodynamic code GIZMO. We include the effect of star formation (SF), modelled with a physically motivated prescription independent of H$_2$, supernova feedback and mass losses from low-mass stars, extragalactic and local stellar radiation, and dust and H$_2$ shielding, to investigate the emergence of the observed correlation between H$_2$ and SF rate surface densities. We present two different sub-grid models and compare them with on-the-fly radiative transfer (RT) calculations, to assess the main differences and limits of the different approaches. We also discuss a sub-grid clumping factor model to enhance the H$_2$ formation, consistent with our SF prescription, which is crucial, at the achieved resolution, to reproduce the correlation with H$_2$. We find that both sub-grid models perform very well relative to the RT simulation, giving comparable results, with moderate differences, but at much lower computational cost. We also find that, while the Kennicutt-Schmidt relation for the total gas is not strongly affected by the different ingredients included in the simulations, the H$_2$-based counterpart is much more sensitive, because of the crucial role played by the dissociating radiative flux and the gas shielding.
We introduce a simple analytic model of galaxy formation that links the growth of dark matter haloes in a cosmological background to the build-up of stellar mass within them. The model aims to identify the physical processes that drive the galaxy-halo co-evolution through cosmic time. The model restricts the role of baryonic astrophysics to setting the relation between galaxies and their haloes. Using this approach, galaxy properties can be directly predicted from the growth of their host dark matter haloes. We explore models in which the effective star formation efficiency within haloes is a function of mass (or virial temperature) and independent of time. Despite its simplicity, the model reproduces self-consistently the shape and evolution of the cosmic star formation rate density, the specific star formation rate of galaxies, and the galaxy stellar mass function, both at the present time and at high redshifts. By systematically varying the effective star formation efficiency in the model, we explore the emergence of the characteristic shape of the galaxy stellar mass function. The origin of the observed double Schechter function at low redshifts is naturally explained by two efficiency regimes in the stellar to halo mass relation, namely, a stellar feedback regulated stage, and a supermassive black hole regulated stage. By providing a set of analytic differential equations, the model can be easily extended and inverted, allowing the roles and impact of astrophysics and cosmology to be explored and understood.
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 SN 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 solar masses, 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. However, in both idealized single-SN simulations and cosmological galaxy-formation simulations, the FIRE-2 algorithm converges much faster than other sub-grid models without re-tuning parameters.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
