TY - JOUR A1 - Werhahn, Maria A1 - Pfrommer, Christoph A1 - Girichidis, Philipp A1 - Puchwein, Ewald A1 - Pakmor, Rüdiger T1 - Cosmic rays and non-thermal emission in simulated galaxies BT - I. Electron and proton spectra compared to Voyager-1 data JF - Monthly notices of the Royal Astronomical Society N2 - 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. KW - astroparticle physics KW - MHD KW - methods: numerical KW - cosmic rays KW - local KW - interstellar matter Y1 - 2021 U6 - https://doi.org/10.1093/mnras/stab1324 SN - 0035-8711 SN - 1365-2966 VL - 505 IS - 3 SP - 3273 EP - 3294 PB - Oxford University Press CY - Oxford ER - TY - JOUR A1 - Werhahn, Maria A1 - Pfrommer, Christoph A1 - Girichidis, Philipp T1 - Cosmic rays and non-thermal emission in simulated galaxies - III. Probing cosmic-ray calorimetry with radio spectra and the FIR-radio correlation JF - Monthly notices of the Royal Astronomical Society N2 - 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. KW - MHD KW - methods: numerical KW - cosmic rays KW - galaxies: magnetic fields KW - galaxies: starburst KW - radio continuum: galaxies Y1 - 2021 U6 - https://doi.org/10.1093/mnras/stab2535 SN - 0035-8711 SN - 1365-2966 VL - 508 IS - 3 SP - 4072 EP - 4095 PB - Oxford Univ. Press CY - Oxford ER - TY - THES A1 - Werhahn, Maria T1 - Simulating galaxy evolution with cosmic rays: the multi-frequency view T1 - Galaxienentwicklung in Simulationen mit komischer Strahlung und deren Strahlungsprozessen N2 - Cosmic rays (CRs) constitute an important component of the interstellar medium (ISM) of galaxies and are thought to play an essential role in governing their evolution. In particular, they are able to impact the dynamics of a galaxy by driving galactic outflows or heating the ISM and thereby affecting the efficiency of star-formation. Hence, in order to understand galaxy formation and evolution, we need to accurately model this non-thermal constituent of the ISM. But except in our local environment within the Milky Way, we do not have the ability to measure CRs directly in other galaxies. However, there are many ways to indirectly observe CRs via the radiation they emit due to their interaction with magnetic and interstellar radiation fields as well as with the ISM. In this work, I develop a numerical framework to calculate the spectral distribution of CRs in simulations of isolated galaxies where a steady-state between injection and cooling is assumed. Furthermore, I calculate the non-thermal emission processes arising from the modelled CR proton and electron spectra ranging from radio wavelengths up to the very high-energy gamma-ray regime. I apply this code to a number of high-resolution magneto-hydrodynamical (MHD) simulations of isolated galaxies, where CRs are included. This allows me to study their CR spectra and compare them to observations of the CR proton and electron spectra by the Voyager-1 satellite and the AMS-02 instrument in order to reveal the origin of the measured spectral features. Furthermore, I provide detailed emission maps, luminosities and spectra of the non-thermal emission from our simulated galaxies that range from dwarfs to Milk-Way analogues to starburst galaxies at different evolutionary stages. I successfully reproduce the observed relations between the radio and gamma-ray luminosities with the far-infrared (FIR) emission of star-forming (SF) galaxies, respectively, where the latter is a good tracer of the star-formation rate. I find that highly SF galaxies are close to the limit where their CR population would lose all of their energy due to the emission of radiation, whereas CRs tend to escape low SF galaxies more quickly. On top of that, I investigate the properties of CR transport that are needed in order to match the observed gamma-ray spectra. Furthermore, I uncover the underlying processes that enable the FIR-radio correlation (FRC) to be maintained even in starburst galaxies and find that thermal free-free-emission naturally explains the observed radio spectra in SF galaxies like M82 and NGC 253 thus solving the riddle of flat radio spectra that have been proposed to contradict the observed tight FRC. Lastly, I scrutinise the steady-state modelling of the CR proton component by investigating for the first time the influence of spectrally resolved CR transport in MHD simulations on the hadronic gamma-ray emission of SF galaxies revealing new insights into the observational signatures of CR transport both spectrally and spatially. N2 - Kosmische Strahlung (CR) ist ein essentieller Bestandteil des interstellaren Mediums (ISM) von Galaxien und spielt eine wichtige Rolle in deren Entwicklung. Insbesondere ist sie in der Lage, die Dynamik einer Galaxie zu beeinflussen, indem sie galaktische Ausflüsse treibt oder das ISM aufheizt und sich dadurch auf die Effizienz der Sternentstehung auswirkt. Um Galaxienentstehung zu verstehen ist es daher notwendig, diesen nicht-thermischen Bestandteil des ISM genau zu modellieren. Aber außerhalb unserer lokalen Umgebung innerhalb der Milchstraße haben wir keine Möglichkeit, um CRs in anderen Galaxien direkt zu messen. Allerdings gibt es viele Möglichkeiten, CRs indirekt über die Strahlung zu beobachten, die sie auf Grund ihrer Interaktion mit Magnetfeldern und interstellarer Strahlung sowie mit dem ISM emittieren. In dieser Arbeit habe ich einen numerischen Code entwickelt, der die spektrale Verteilung der CRs in Simulationen von isolierten Galaxien berechnet, wobei ein stationäres Verhältnis zwischen Injektion und Kühlen angenommen wird. Des Weiteren berechnet er die nicht-thermischen Strahlungsprozesse, die aus den modellierten CR Protonen- und Elektronenspektren hervorgehen. Diese reichen von Radiowellenlängen bis hin zu hochenergetischer Gammastrahlung. Ich wende diesen Code auf eine Vielzahl von hoch aufgelösten, magneto-hydrodynamischen Simulationen von isolierten Galaxien an, die CRs beinhalten. Das ermöglicht es mir, ihre CR Spektren zu untersuchen und mit Beobachtungen des Voyager-1 Satelliten sowie des AMS-02 Instruments von CR Protonen- und Elektronenspektren zu vergleichen, um dem Ursprung von den gemessenen spektralen Besonderheiten nachzugehen. Außerdem lege ich detaillierte Emissionskarten, Leuchtkräfte und Spektren der nicht-thermischen Strahlung unserer simulierten Galaxien vor, die von Zwerggalaxien über Milchstraßen-ähnliche Galaxien bis hin zu Starburst-Galaxien bei verschiedensten Entwicklungsstadien reichen. Damit kann ich erfolgreich die beobachteten Zusammenhänge zwischen jeweils der Radio- und Gammastrahlungsleuchtkraft mit der Ferninfrarot (FIR) Strahlung der sternbildenden Galaxien reproduzieren, wobei die FIR Strahlung ein guter Indikator für die Rate der Sternentstehung ist. Dabei finde ich heraus, dass Galaxien mit einer hohen Rate an Sternentstehung sehr nah an dem Limit sind, in dem ihre CR Population all ihre Energie an die Produktion von Strahlung verlieren würde, während CRs dazu tendieren, Galaxien mit einer niedrigen Sternentstehungsrate schneller zu verlassen. Zusätzlich untersuche ich die Eigenschaften des Transports von CRs, die benötigt werden, um die beobachteten Spektren der Gammastrahlung zu reproduzieren. Außerdem decke ich die zugrundeliegenden physikalischen Prozesse auf, durch die die Korrelation zwischen der FIR- und Radioleuchtkraft auch in Starburst-Galaxien aufrecht erhalten werden kann und finde heraus, dass die thermische Emission natürlicherweise die beobachteten Radiospektren in Galaxien wie M82 und NGC 253 erklärt, wodurch sich das Rätsel der flachen Radiospektren löst, die scheinbar im Widerspruch zum beobachteten engen Zusammenhang zwischen der FIR- und Radioleuchtkraft standen. Zuletzt hinterfrage ich die Annahme eines stationären Zustandes bei der Modellierung der CR Protonenspektren, indem ich zum ersten Mal den Einfluss von spektral aufgelöstem Transport von CR Protonen in magneto-hydrodynamischen Simulationen auf die hadronische Gammastrahlung von sternbildenden Galaxien untersuche, was neue Einblicke in beobachtbare Signaturen, sowohl spektral als auch räumlich, von CR-Transport ermöglicht. KW - galaxies KW - cosmic rays KW - numerical astrophysics KW - kosmische Strahlung KW - Galaxien KW - numerische Astrophysik Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-572851 ER -