@article{WerhahnPfrommerGirichidisetal.2021,
  author    = {Werhahn, Maria and Pfrommer, Christoph and Girichidis, Philipp and Puchwein, Ewald and Pakmor, R{\"u}diger},
  title     = {Cosmic rays and non-thermal emission in simulated galaxies},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {505},
  journal   = {Monthly notices of the Royal Astronomical Society},
  number    = {3},
  publisher = {Oxford University Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/stab1324},
  pages     = {3273 -- 3294},
  year      = {2021},
  abstract  = {Current-day cosmic ray (CR) propagation studies use static Milky Way models and fit parametrized source distributions to data. Instead, we use three-dimensional magnetohydrodynamic (MHD) simulations of isolated galaxies with the moving-mesh code arepo that self-consistently accounts for hydrodynamic effects of CR protons. In post-processing, we calculate their steady-state spectra, taking into account all relevant loss processes. We show that this steady-state assumption is well justified in the disc and generally for regions that emit non-thermal radio and gamma rays. Additionally, we model the spectra of primary electrons, accelerated by supernova remnants, and secondary electrons and positrons produced in hadronic CR proton interactions with the gas. We find that proton spectra above 10 GeV only weakly depend on galactic radius, while they acquire a radial dependence at lower energies due to Coulomb interactions. Radiative losses steepen the spectra of primary CR electrons in the central galactic regions, while diffusive losses dominate in the outskirts. Secondary electrons exhibit a steeper spectrum than primaries because they originate from the transported steeper CR proton spectra. Consistent with Voyager-1 and AMS-02 data, our models (i) show a turnover of proton spectra below GeV energies due to Coulomb interactions so that electrons start to dominate the total particle spectra and (ii) match the shape of the positron fraction up to 10 GeV. We conclude that our steady-state CR modelling in MHD CR galaxy simulations is sufficiently realistic to capture the dominant transport effects shaping their spectra, arguing for a full MHD treatment to accurately model CR transport in the future.},
  language  = {en}
}
@article{WerhahnPfrommerGirichidis2021,
  author    = {Werhahn, Maria and Pfrommer, Christoph and Girichidis, Philipp},
  title     = {Cosmic rays and non-thermal emission in simulated galaxies - III. Probing cosmic-ray calorimetry with radio spectra and the FIR-radio correlation},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {508},
  journal   = {Monthly notices of the Royal Astronomical Society},
  number    = {3},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/stab2535},
  pages     = {4072 -- 4095},
  year      = {2021},
  abstract  = {An extinction-free estimator of the star formation rate (SFR) of galaxies is critical for understanding the high-redshift universe. To this end, the nearly linear, tight correlation of far-infrared (FIR), and radio luminosity of star-forming galaxies is widely used. While the FIR is linked to massive star formation, which also generates shock-accelerated cosmic-ray (CR) electrons and radio synchrotron emission, a detailed understanding of the underlying physics is still lacking. Hence, we perform three-dimensional magnetohydrodynamical (MHD) simulations of isolated galaxies over a broad range of halo masses and SFRs using the moving-mesh code AREPO, and evolve the CR proton energy density self-consistently. In post-processing, we calculate the steady-state spectra of primary, shock-accelerated and secondary CR electrons, which result from hadronic CR proton interactions with the interstellar medium. The resulting total radio luminosities correlate with the FIR luminosities as observed and are dominated by primary CR electrons if we account for anisotropic CR diffusion. The increasing contribution of secondary emission up to 30 per cent in starbursts is compensated by the larger bremsstrahlung and Coulomb losses. CR electrons are in the calorimetric limit and lose most of their energy through inverse Compton interactions with star light and cosmic microwave background (CMB) photons while less energy is converted into synchrotron emission. This implies steep steady-state synchrotron spectra in starbursts. Interestingly, we find that thermal free-free emission flattens the total radio spectra at high radio frequencies and reconciles calorimetric theory with observations while free-free absorption explains the observed low-frequency flattening towards the central regions of starbursts.},
  language  = {en}
}