No Arabic abstract
Recent efforts in cosmic ray (CR) confinement and transport theory are discussed. Three problems are addressed as being crucial for understanding the present day observations and their possible telltale signs of the CR origin. The first problem concerns CR behavior right after their release from a source, such as a supernova remnant (SNR). At this phase the CRs are confined near the source by self-emitted Alfven waves. The second is the problem of diffusive propagation of CRs through the turbulent ISM. This is a seemingly straightforward and long-resolved problem, but it remains controversial and reveals paradoxes. A resolution based on the Chapman-Enskog asymptotic CR transport analysis, that also includes magnetic focusing, is suggested. The third problem is about a puzzling sharp ($sim10^{circ}$) anisotropies in the CR arrival directions that might bear on important clues of their transport between the source and observer. The overarching goal is to improve our understanding of all aspects of the CRs source escape and ensuing propagation through the galaxy to the level at which their sources can be identified observationally.
The microphysics of ~GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant-diffusivity). We explore different physically-motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological FIRE simulations of dwarf and ~$L_{ast}$ galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfven-wave damping, etc. We self-consistently predict observables including $gamma$-rays ($L_{gamma}$), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multi-phase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or $L_{gamma}$), as transport is rate-limited by the ionized WIM and inner CGM gaseous halo ($10^{4}-10^{6}$ K gas within 10-30 kpc), and $L_{gamma}$ can be dominated by trapping in small patches. Most physical ET models contribute negligible scattering of ~1-10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely-assumed scalings for SC models produce excessive confinement by factors >100 in the WIM and inner CGM, where turbulent and Landau damping dominate. This suggests either a breakdown of quasi-linear theory used to derive the CR transport parameters in SC, or that other novel damping mechanisms dominate in intermediate-density ionized gas.
This Commentary discusses the widespread impact of the milestone 1980 paper by Sunyaev and Titarchuk on Comptonization and points out a new possible application to the temporal broadening of the most distant Gamma Ray Bursts.
Supernova remnants (SNRs) are the prime candidates for the acceleration of the Galactic Cosmic Rays. Tracers for interactions of Cosmic Rays with ambient material are gamma rays at TeV energies, which can be observed with ground based Cherenkov telescopes like H.E.S.S. In the recent years H.E.S.S. has detected several SNRs and interactions of SNRs with molecular clouds. Here the current results of these observations are presented and possible leptonic and hadronic scenarios are discussed. It is shown that it is likely that SNRs are the sources of Galactic Cosmic Rays.
This article aims at establishing new benchmark scenarios for Galactic cosmic-ray propagation in the GV-TV rigidity range, based on fits to the AMS-02 B/C data with the USINE v3.5 propagation code. We employ a new fitting procedure, cautiously taking into account data systematic error correlations in different rigidity bins and considering Solar modulation potential and leading nuclear cross-section as nuisance parameters. We delineate specific low, intermediate, and high-rigidity ranges that can be related to both features in the data and peculiar microphysics mechanisms resulting in spectral breaks. We single out a scenario which yields excellent fits to the data and includes all the presumably relevant complexity, the BIG model. This model has two limiting regimes: (i) the SLIM model, a minimal diffusion-only setup, and (ii) the QUAINT model, a convection-reacceleration model where transport is tuned by non-relativistic effects. All models lead to robust predictions in the high-energy regime ($gtrsim10$GV), i.e. independent of the propagation scenario: at $1sigma$, the diffusion slope $delta$ is $[0.43-0.53]$, whereas $K_{10}$, the diffusion coefficient at 10GV, is $[0.26-0.36]$kpc$^2$Myr$^{-1}$; we confirm the robustness of the high-energy break, with a typical value $Delta_hsim 0.2$. We also find a hint for a similar (reversed) feature at low rigidity around the B/C peak ($sim 4$GV) which might be related to some effective damping scale in the magnetic turbulence.
Numerical solutions of the cosmic-ray (CR) magneto-hydrodynamic equations are dogged by a powerful numerical instability, which arises from the constraint that CRs can only stream down their gradient. The standard cure is to regularize by adding artificial diffusion. Besides introducing ad-hoc smoothing, this has a significant negative impact on either computational cost or complexity and parallel scalings. We describe a new numerical algorithm for CR transport, with close parallels to two moment methods for radiative transfer under the reduced speed of light approximation. It stably and robustly handles CR streaming without any artificial diffusion. It allows for both isotropic and field-aligned CR streaming and diffusion, with arbitrary streaming and diffusion coefficients. CR transport is handled explicitly, while source terms are handled implicitly. The overall time-step scales linearly with resolution (even when computing CR diffusion), and has a perfect parallel scaling. It is given by the standard Courant condition with respect to a constant maximum velocity over the entire simulation domain. The computational cost is comparable to that of solving the ideal MHD equation. We demonstrate the accuracy and stability of this new scheme with a wide variety of tests, including anisotropic streaming and diffusion tests, CR modified shocks, CR driven blast waves, and CR transport in multi-phase media. The new algorithm opens doors to much more ambitious and hitherto intractable calculations of CR physics in galaxies and galaxy clusters. It can also be applied to other physical processes with similar mathematical structure, such as saturated, anisotropic heat conduction.