Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCosmological structure formation shocks and cosmic rays in hydrodynamical simulations
116
0
0.0
(
0
)
Authors
Christoph Pfrommer
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
We present new results characterizing cosmological shocks within adaptive mesh refinement N-Body/hydrodynamic simulations that are used to predict non-thermal components of large-scale structure. This represents the first study of shocks using adaptive mesh refinement. We propose a modified algorithm for finding shocks from those used on unigrid simulations that reduces the shock frequency of low Mach number shocks by a factor of ~3. We then apply our new technique to a large, (512 Mpc/h)^3, cosmological volume and study the shock Mach number (M) distribution as a function of pre-shock temperature, density, and redshift. Because of the large volume of the simulation, we have superb statistics that results from having thousands of galaxy clusters. We find that the Mach number evolution can be interpreted as a method to visualize large-scale structure formation. Shocks with Mach<5 typically trace mergers and complex flows, while 5<Mach<20 and Mach>20 generally follow accretion onto filaments and galaxy clusters, respectively. By applying results from nonlinear diffusive shock acceleration models using the first-order Fermi process, we calculate the amount of kinetic energy that is converted into cosmic ray protons. The acceleration of cosmic ray protons is large enough that in order to use galaxy clusters as cosmological probes, the dynamic response of the gas to the cosmic rays must be included in future numerical simulations.
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)
We study the properties of two bars formed in fully cosmological hydrodynamical simulations of the formation of Milky Way-mass galaxies. In one case, the bar formed in a system with disc, bulge and halo components and is relatively strong and long, as could be expected for a system where the spheroid strongly influences the evolution. The second bar is less strong, shorter, and formed in a galaxy with no significant bulge component. We study the strength and length of the bars, the stellar density profiles along and across the bars and the velocity fields in the bar region. We compare them with the results of dynamical (idealised) simulations and with observations, and find, in general, a good agreement, although we detect some important differences as well. Our results show that more or less realistic bars can form naturally in a $Lambda$CDM cosmology, and open up the possibility to study the bar formation process in a more consistent way than previously done, since the host galaxies grow, accrete matter and significantly evolve during the formation and evolution of the bar.
We present a study of hydrodynamic drag forces in smoothed particle simulations. In particular, the deceleration of a resolution-limited cold clump of gas moving through a hot medium is examined. It is found that the drag for subsonic velocities exceeds that predicted by simple physical approximations. This is shown to be a result of the hydrodynamical method which encourages the accretion of particles from the hot medium onto a shell around the cold clump, effectively increasing the radius of the clump. For sonic and supersonic velocities, the drag is shown to be dependent on the effective cross section of the clump. The consequences for cosmological simulations are discussed.
State-of-the-art cosmological hydrodynamical simulations of galaxy formation have reached the point at which their outcomes result in galaxies with ever more realism. Still, the employed sub-grid models include several free parameters such as the density threshold, $n$, to localize the star-forming gas. In this work, we investigate the possibilities to utilize the observed clustered nature of star formation (SF) in order to refine SF prescriptions and constrain the density threshold parameter. To this end, we measure the clustering strength, correlation length and power-law index of the two-point correlation function of young ($tau<50$ Myr) stellar particles and compare our results to observations from the HST Legacy Extragalactic UV Survey (LEGUS). Our simulations reveal a clear trend of larger clustering signal and power-law index and lower correlation length as the SF threshold increases with only mild dependence on galaxy properties such as stellar mass or specific star formation rate. In conclusion, we find that the observed clustering of SF is inconsistent with a low threshold for SF ($n<1$ cm$^{-3}$) and strongly favours a high value for the density threshold of SF ($n>10$ cm$^{-3}$), as for example employed in the NIHAO project.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
