Abstract

We review numerical methods for simulations of cosmic ray (CR) propagation on galactic and larger scales. We present the development of algorithms designed for phenomenological and self-consistent models of CR propagation in kinetic description based on numerical solutions of the Fokker–Planck equation. The phenomenological models assume a stationary structure of the galactic interstellar medium and incorporate diffusion of particles in physical and momentum space together with advection, spallation, production of secondaries and various radiation mechanisms. The self-consistent propagation models of CRs include the dynamical coupling of the CR population to the thermal plasma. The CR transport equation is discretized and solved numerically together with the set of MHD equations in various approaches treating the CR population as a separate relativistic fluid within the two-fluid approach or as a spectrally resolved population of particles evolving in physical and momentum space. The relevant processes incorporated in self-consistent models include advection, diffusion and streaming propagation as well as adiabatic compression and several radiative loss mechanisms. We discuss, applications of the numerical models for the interpretation of CR data collected by various instruments. We present example models of astrophysical processes influencing galactic evolution such as galactic winds, the amplification of large-scale magnetic fields and instabilities of the interstellar medium.

Highlights

  • 1.1 IntroductionCosmic rays (CRs) are charged particles with non-thermal energy distributions (Strong et al 2007; Grenier et al 2015; Gabici et al 2019)

  • If the coupling of CRs is reduced, due to wave damping mechanisms, the effective propagation speed of the CR population might be higher than the Alfven speed, it is argued that the heating rate is still given by the expression (8), because all momentum and energy transfer between the CRs and gas is mediated by hydromagnetic waves which propagate at the speed vA (e.g. Zweibel 2013)

  • We have reviewed the numerical treatment of CRs assuming a fluid description of this high-energy component

Read more

Summary

Introduction

Cosmic rays (CRs) are charged particles with non-thermal energy distributions (Strong et al 2007; Grenier et al 2015; Gabici et al 2019). There are both hadronic and leptonic cosmic rays among which protons and electrons are the most abundant particles. Naab and Ostriker 2017; Zweibel 2017) In addition they might have an impact on the distribution of gas in the galactic disc and alter the star formation process on molecular cloud scales, even though their dynamical impact on molecular clouds is expected to be much weaker compared to other drivers in the ISM like radiation or supernovae. If CRs are efficiently coupled to the gas or penetrate into regions dense enough such that direct particle-particle collisions become relevent, they provide a temperature floor, which directly impacts the fragmentation scale and the seeds of star formation

Origin of CRs
Numerical approaches of CR transport and scope of this review
Basics and approaches
Particle versus kinetic versus fluid approach
Grey versus spectrally resolved CR fluids
The canonical CR propagation equation
Page 10 of 92
Streaming and diffusion
Page 12 of 92
Reacceleration
Page 14 of 92
Weighted slabs and leaky boxes
Page 16 of 92
Explicit models
GALPROP
Page 18 of 92
Page 20 of 92
Numerical solution of the propagation equation
Page 22 of 92
Page 24 of 92
Other codes
The system of equations
CR number density
CR energy density
Page 26 of 92
Two-fluid diffusion-advection models
Page 28 of 92
Riemann problem
Page 30 of 92
Extensions of the Riemann problem including CRs
CR injection in SN shocks
Page 34 of 92
Unified Alfven wave regulated CR transport
Page 36 of 92
Page 38 of 92
Momentum-dependent diffusion-advection models
One-moment approaches
Page 42 of 92
Two-moment approaches
Page 44 of 92
Page 46 of 92
Advances in theory of CR propagation
Page 48 of 92
Spatial diffusion of spectrally resolved CRs
D31 D32 D33
Page 50 of 92
Numerical treatment of CR diffusion
Anisotropic diffusion on a regular mesh
Page 52 of 92
Page 54 of 92
A semi-implicit extension
Page 56 of 92
Diffusion on a Voronoi mesh
Evolution of the isotropic CR spectrum on the momentum grid
Page 58 of 92
Page 60 of 92
Conversion between quantities
Examples of GALPROP model predictions
Diffusive reacceleration and alternatives
CRs in star formation
Page 66 of 92
CRs in the interstellar medium
CR driven large-scale instabilities of the ISM
Page 68 of 92
Cosmic ray driven galactic dynamo
Page 70 of 92
Page 72 of 92
Cosmic ray driven galactic winds
Page 74 of 92
Cosmic rays in galaxy clusters
Page 78 of 92
Conclusions and outlook
Page 80 of 92
Page 82 of 92
Page 84 of 92
Page 86 of 92
Page 88 of 92
Page 90 of 92
Findings
Page 92 of 92
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call