Search

Cosmic rays (CRs) are dynamically important in the evolution of galaxies by regulating star formation and powering galactic outflows. However, to what extent CRs regulate galaxy formation depends on the coupling strength of CRs with the ambient plasma and the effective CR transport speed. Moreover, both properties sensitively depend on the CR momentum, which is largely unexplored in three-dimensional hydrodynamical simulations. We perform magnetohydrodynamical simulations of entire galaxies with masses ranging from 1010 to |$10^{12}\, \mathrm{M}_\odot$| and compare dynamically coupled CRs in the grey approximation with a spectrally resolved model that includes CR momenta from |$0.1\, \mathrm{GeV} c^{-1}$| to |$100\, \mathrm{TeV} c^{-1}$|⁠. We find that hadronic cooling of CRs dominates over Alfvén cooling, with the latter emulating CR losses as a result of streaming of CRs down their pressure gradient. While star formation rates and galaxy morphologies are only mildly affected by the spectral CR modelling, mass loading factors of galactic outflows can differ by up to a factor of 4 in dwarf galaxies. All simulated low-mass haloes (M  = 1010, 1011, and |$3\times 10^{11}\, \mathrm{M}_\odot$|⁠) drive strong outflows, where CR transport is temporally dominated by advection. In contrast, the Milky Way-mass galaxy with |$M=10^{12}\, \mathrm{M}_\odot$| does not drive sustained outflows, so that CR transport is entirely dominated by diffusion. The effective energy weighted diffusion coefficients vary by two orders of magnitude from the canonical energy-weighted values of |$\langle {D}\rangle _{e_\mathrm{cr}}\sim 10^{28}\, \mathrm{cm^2\, s^{-1}}$| in the disc up to |$3\times 10^{29}\, \mathrm{cm^2\, s^{-1}}$| in the circumgalactic medium, where we observe substantial temperature and CR pressure differences between our grey and spectral CR models. [ABSTRACT FROM AUTHOR]

Cosmic ray (CR) protons are an important component in many astrophysical systems. Processes like CR injection, cooling, adiabatic changes as well as active CR transport through the medium strongly modify the CR momentum distribution and have to be taken into account in hydrodynamical simulations. We present an efficient novel numerical scheme to accurately compute the evolution of the particle distribution function by solving the Fokker–Planck equation with a low number of spectral bins (10–20), which is required to include a full spectrum for every computational fluid element. The distribution function is represented by piecewise power laws and is not forced to be continuous, which enables an optimal representation of the spectrum. The Fokker–Planck equation is solved with a two-moment approach evolving the CR number and energy density. The low numerical diffusion of the scheme reduces the numerical errors by orders of magnitude in comparison to classical schemes with piecewise constant spectral representations. With this method not only the spectral evolution of CRs can be computed accurately in magnetohydrodynamic simulations but also their dynamical impact as well as CR ionization. This allows for more accurate models for astrophysical plasmas, like the interstellar medium, and direct comparisons with observations. [ABSTRACT FROM AUTHOR]

We present three-dimensional radiation-hydrodynamical simulations of the impact of stellar winds, photoelectric heating, photodissociating and photoionizing radiation, and supernovae on the chemical composition and star formation in a stratified disc model. This is followed by a sink-based model for star clusters with populations of individual massive stars. Stellar winds and ionizing radiation regulate the star formation rate at a factor of ~10 below the simulation with only supernova feedback due to their immediate impact on the ambient interstellar medium after star formation. Ionizing radiation (with winds and supernovae) significantly reduces the ambient densities for most supernova explosions to ρ < 10-25 g cm-3, compared to 10-23g cm-3 for the model with only winds and supernovae. Radiation from massive stars reduces the amount of molecular hydrogen and increases the neutral hydrogen mass and volume filling fraction. Only this model results in a molecular gas depletion time-scale of 2 Gyr and shows the best agreement with observations. In the radiative models, the Hα emission is dominated by radiative recombination as opposed to collisional excitation (the dominant emission in non-radiative models), which only contributes ~1-10 per cent to the total Hα emission. Individual massive stars (M ≥ 30M☉) with short lifetimes are responsible for significant fluctuations in the Hα luminosities. The corresponding inferred star formation rates can underestimate the true instantaneous star formation rate by a factor of ~10. [ABSTRACT FROM AUTHOR]

ABSTRACT We investigate the impact of different initial conditions for the initial density profile and the initial turbulence on the formation process of protostellar clusters. We study the collapse of dense molecular cloud cores with three-dimensional adaptive mesh refinement simulations. We focus our discussion on the distribution of the gas among the protostellar objects in the turbulent dynamical cores. Despite the large variations in the initial configurations and the resulting gas and cluster morphology we find that all stellar clusters follow a very similar gas accretion behaviour. Once secondary protostars begin to form, the central region of a cluster is efficiently shielded from further accretion. Hence, objects located close to the centre are starved of material, as indicated by a strong decrease of the central accretion rate. This fragmentation-induced starvation occurs not only in rotationally supported discs and filaments, but also in more spherically symmetric clusters with complex chaotic motions. [ABSTRACT FROM AUTHOR]