Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCosmic ray feedback in hydrodynamical simulations of galaxy formation
178
0
0.0
(
0
)
Authors
Martin Jubelgas
Ask ChatGPT about the research
No Arabic abstract
It is well known that cosmic rays (CRs) contribute significantly to the pressure of the interstellar medium in our own Galaxy, suggesting that they may play an important role in regulating star formation during the formation and evolution of galaxies. We here discuss a novel numerical treatment of the physics of CRs and its implementation in the parallel smoothed particle hydrodynamics code GADGET-2. In our methodology, the non-thermal CR population of each gaseous fluid element is approximated by a simple power law spectrum in particle momentum, characterized by an amplitude, a cut-off, and a fixed slope. Adiabatic compression, and a number of physical source and sink terms are modelled which modify the CR pressure of each particle. The most important sources considered are injection by supernovae and diffusive shock acceleration, while the primary sinks are thermalization by Coulomb interactions, and catastrophic losses by hadronic interactions. We also include diffusion of CRs. Our scheme allows us to carry out the first cosmological structure formation simulations that self-consistently account for CR physics. In simulations of isolated galaxies, we find that CRs can significantly reduce the star formation efficiencies of small galaxies, with virial velocities below ~80 km/s, an effect that becomes progressively stronger towards low mass scales. In cosmological simulations at high redshift, the total mass-to-light ratio of small halos and the faint-end of the luminosity function are strongly affected. When CR acceleration in shocks is followed as well, up to ~40% of the energy dissipated at structure formation shocks can appear as CR pressure at z~3-6, but this fraction drops to ~10% at low redshifts when the shock distribution becomes increasingly dominated by lower Mach numbers. (abridged)
rate research
Read More
It is well known that cosmic rays (CRs) contribute significantly to the pressure of the interstellar medium in our own Galaxy, suggesting that they may play an important role in regulating star formation during the formation and evolution of galaxies. We will present a novel numerical treatment of the physics of CRs and its implementation in the parallel smoothed particle hydrodynamics (SPH) code GADGET-2. In our methodology, the non-thermal CR population is treated self-consistently in order to assess its dynamical impact on the thermal gas as well as other implications on cosmological observables. In simulations of galaxy formation, we find that CRs can significantly reduce the star formation efficiencies of small galaxies. This effect becomes progressively stronger towards low mass scales. In cosmological simulations of the formation of dwarf galaxies at high redshift, we find that the total mass-to-light ratio of small halos and the faint-end of the luminosity function are affected. In high resolution simulations of galaxy clusters, we find lower contributions of CR pressure, due to the smaller CR injection efficiencies at low Mach number flow shocks inside halos, and the softer adiabatic index of CRs, which disfavours them when a composite of thermal gas and CRs is adiabatically compressed. Within cool core regions, the CR pressure reaches equipartition with the thermal pressure leading to an enhanced compressibility of the central intra-cluster medium, an effect that increases the central density and pressure of the gas. While the X-ray luminosity in low mass cool core clusters is boosted, the integrated Sunyaev-Zeldovich effect is only slightly changed. The resolved Sunyaev-Zeldovich maps, however, show a larger variation with an increased central flux decrement.
We investigate a numerical model for AGN feedback where for the first time a relativistic particle population in AGN-inflated bubbles is followed within a full cosmological context. In our high-resolution simulations of galaxy cluster formation, we assume that BH accretion is accompanied by energy feedback that occurs in two different modes, depending on the accretion rate itself. Unlike in previous work, we inject a non-thermal particle population of relativistic protons into the AGN bubbles, instead of adopting a purely thermal heating. We then follow the subsequent evolution of the cosmic ray (CR) plasma inside the bubbles, considering both its hydrodynamical interactions and dissipation processes relevant for the CR population. Due to the different buoyancy of relativistic plasma and the comparatively long CR dissipation timescale we find substantial changes in the evolution of clusters as a result of CR feedback. In particular, the non-thermal population can provide significant pressure support in central cluster regions at low thermal temperatures, providing a natural explanation for the decreasing temperature profiles found in cool core clusters. At the same time, the morphologies of the bubbles and of the induced X-ray cavities show a striking similarity to observational findings. AGN feedback with CRs also proves efficient in regulating cluster cooling flows so that the total baryon fraction in stars becomes limited to realistic values of the order of 10%. We find that the partial CR support of the intracluster gas also affects the expected signal of the thermal Sunyaev-Zeldovich effect, with typical modifications of the integrated Compton-y parameter within the virial radius of the order of 10%. [Abridged]
Cosmological shock waves during structure formation not only play a decisive role for the thermalization of gas in virializing structures but also for the acceleration of relativistic cosmic rays (CRs) through diffusive shock acceleration. We discuss a novel numerical treatment of the physics of cosmic rays in combination with a formalism for identifying and measuring the shock strength on-the-fly during a smoothed particle hydrodynamics simulation. In our methodology, the non-thermal CR population is treated self-consistently in order to assess its dynamical impact on the thermal gas as well as other implications on cosmological observables. Using this formalism, we study the history of the thermalization process in high-resolution hydrodynamic simulations of the Lambda cold dark matter model. Collapsed cosmological structures are surrounded by shocks with high Mach numbers up to 1000, but they play only a minor role in the energy balance of thermalization. However, this finding has important consequences for our understanding of the spatial distribution of CRs in the large-scale structure. In high resolution simulations of galaxy clusters, we find a low contribution of the averaged CR pressure, due to the small acceleration efficiency of lower Mach numbers of flow shocks inside halos and the softer adiabatic index of CRs. However, within cool core regions, the CR pressure reaches equipartition with the thermal pressure leading there to a lower effective adiabatic index and thus to an enhanced compressibility of the central intracluster medium. This effect increases the central density and pressure of the cluster and thus the resulting X-ray emission and the central Sunyaev-Zeldovich flux decrement. The integrated Sunyaev-Zeldovich effect, however, is only slightly changed.
It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in Galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In Eagle, stars form preferentially out of low specific angular momentum gas in the interstellar medium (ISM) due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between Galform and Eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in Eagle is slower than is assumed in Galform. Finally, by matching individual galaxies between Eagle and Galform, we find that an artificial dependence of AGN feedback and gas infall rates on halo mass doubling events in Galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the Galform semi-analytic model can be significantly improved in light of recent advances.
Various heuristic approaches to model unresolved supernova (SN) feedback in galaxy formation simulations exist to reproduce the formation of spiral galaxies and the overall inefficient conversion of gas into stars. Some models, however, require resolution dependent scalings. We present a sub-resolution model representing the three major phases of supernova blast wave evolution $-$free expansion, energy conserving Sedov-Taylor, and momentum conserving snowplow$-$ with energy scalings adopted from high-resolution interstellar-medium simulations in both uniform and multiphase media. We allow for the effects of significantly enhanced SN remnant propagation in a multiphase medium with the cooling radius scaling with the hot volume fraction, $f_{mathrm{hot}}$, as $(1 - f_{mathrm{hot}})^{-4/5}$. We also include winds from young massive stars and AGB stars, Stromgren sphere gas heating by massive stars, and a gas cooling limiting mechanism driven by radiative recombination of dense HII regions. We present initial tests for isolated Milky-Way like systems simulated with the GADGET based code SPHgal with improved SPH prescription. Compared to pure thermal SN input, the model significantly suppresses star formation at early epochs, with star formation extended both in time and space in better accord with observations. Compared to models with pure thermal SN feedback, the age at which half the stellar mass is assembled increases by a factor of 2.4, and the mass loading parameter and gas outflow rate from the galactic disk increase by a factor of 2. Simulation results are converged for a two order of magnitude variation in particle mass in the range (1.3$-$130)$times 10^4$ solar masses.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
