Supernova remnants (SNRs) are believed to be the source of Galactic cosmic rays (CRs). SNR shocks accelerate CR protons and electrons which reveal key insights into the non-thermal physics by means of their synchrotron and $\gamma$-ray emission. The remnant SN 1006 is an ideal particle acceleration laboratory because it is observed across all electromagnetic wavelengths from radio to $\gamma$-rays. We perform three-dimensional (3D) magnetohydrodynamics (MHD) simulations where we include CR protons and follow the CR electron spectrum. By matching the observed morphology and non-thermal spectrum of SN 1006 in radio, X-rays and $\gamma$-rays, we gain new insight into CR electron acceleration and magnetic field amplification. 1. We show that a mixed leptonic-hadronic model is responsible for the $\gamma$-ray radiation: while leptonic inverse-Compton emission and hadronic pion-decay emission contribute equally at GeV energies observed by Fermi, TeV energies observed by imaging air Cherenkov telescopes are hadronically dominated. 2. We show that quasi-parallel acceleration (i.e., when the shock propagates at a narrow angle to the upstream magnetic field) is preferred for CR electrons and that the electron acceleration efficiency of radio-emitting GeV electrons at quasi-perpendicular shocks is suppressed at least by a factor ten. This poses a challenge for the interpretation of plasma particle-in-cell simulations of particle acceleration at shocks. 3. To match the radial emission profiles and the $\gamma$-ray spectrum, we require a volume-filling, turbulently amplified magnetic field and that the Bell-amplified magnetic field is damped in the immediate post-shock region. Our work connects micro-scale plasma physics simulations to the scale of SNRs.

Galactic cosmic rays (CRs) are believed to be accelerated at supernova remnant (SNR) shocks. In the hadronic scenario the TeV gamma-ray emission from SNRs originates from decaying pions that are produced in collisions of the interstellar gas and CRs. Using CR-magnetohydrodynamic simulations, we show that magnetic obliquity-dependent shock acceleration is able to reproduce the observed TeV gamma-ray morphology of SNRs such as Vela Jr. and SN1006 solely by varying the magnetic morphology. This implies that gamma-ray bright regions result from quasi-parallel shocks (i.e., when the shock propagates at a narrow angle to the upstream magnetic field), which are known to efficiently accelerate CR protons, and that gamma-ray dark regions point to quasi-perpendicular shock configurations. Comparison of the simulated gamma-ray morphology to observations allows us to constrain the magnetic coherence scale $\lambda_B$ around Vela Jr. and SN1006 to $\lambda_B \simeq 13_{-4.3}^{+13}$ pc and $\lambda_B >200_{-10}^{+80}$ pc, respectively, where the ambient magnetic field of SN1006 is consistent with being largely homogeneous. We find consistent pure hadronic and mixed hadronic-leptonic models that both reproduce the multi-frequency spectra from the radio to TeV gamma rays and match the observed gamma-ray morphology. Finally, to capture the propagation of a SNR shock in a clumpy interstellar medium, we study the interaction of a shock with a dense cloud with numerical simulations and analytics. We construct an analytical gamma-ray model for a core collapse SNR propagating through a structured interstellar medium, and show that the gamma-ray luminosity is only biased by 30% for realistic parameters.

Cosmic ray (CR) electrons reveal key insights into the non-thermal physics of the interstellar medium, galaxies, galaxy clusters, and active galactic nuclei by means of their inverse Compton {\gamma}-ray emission and synchrotron emission in magnetic fields. While magnetohydrodynamical (MHD) simulations with CR protons capture their dynamical impact on these systems, only few computational studies include CR electron physics because of the short cooling time-scales and complex hysteresis effects, which require a numerically expensive, high-resolution spectral treatment. Since CR electrons produce important non-thermal observational signatures, such a spectral CR electron treatment is important to link MHD simulations to observations. We present an efficient post-processing code for Cosmic Ray Electron Spectra that are evolved in Time (CREST) on Lagrangian tracer particles. The CR electron spectra are very accurately evolved on comparably large MHD time steps owing to an innovative hybrid numerical-analytical scheme. CREST is coupled to the cosmological MHD code AREPO and treats all important aspects of spectral CR electron evolution such as adiabatic expansion and compression, Coulomb losses, radiative losses in form of inverse Compton, bremsstrahlung and synchrotron processes, diffusive shock acceleration and reacceleration, and Fermi-II reacceleration. After showing various code validations of idealized one-zone simulations, we study the coupling of CREST to MHD simulations. We demonstrate that the CR electron spectra are efficiently and accurately evolved in shock-tube and Sedov-Taylor blast wave simulations. This opens up the possibility to produce self-consistent synthetic observables of non-thermal emission processes in various astrophysical environments.