No Arabic abstract
The propagation of charged particles, including cosmic rays, in a partially ordered magnetic field is characterized by a diffusion tensor whose components depend on the particles Larmor radius $R_L$ and the degree of order in the magnetic field. Most studies of the particle diffusion presuppose a scale separation between the mean and random magnetic fields (e.g., there being a pronounced minimum in the magnetic power spectrum at intermediate scales). Scale separation is often a good approximation in laboratory plasmas, but not in most astrophysical environments such as the interstellar medium (ISM). Modern simulations of the ISM have numerical resolution of order 1 pc, so the Larmor radius of the cosmic rays that dominate in energy density is at least $10^{6}$ times smaller than the resolved scales. Large-scale simulations of cosmic ray propagation in the ISM thus rely on oversimplified forms of the diffusion tensor. We take the first steps towards a more realistic description of cosmic ray diffusion for such simulations, obtaining direct estimates of the diffusion tensor from test particle simulations in random magnetic fields (with the Larmor radius scale being fully resolved), for a range of particle energies corresponding to $10^{-2}lesssim R_L/l_c lesssim 10^{3}$, where $l_c$ is the magnetic correlation length. We obtain explicit expressions for the cosmic ray diffusion tensor for $R_L/l_c ll 1$, that might be used in a sub-grid model of cosmic ray diffusion. The diffusion coefficients obtained are closely connected with existing transport theories that include the random walk of magnetic lines.
As the fundamental physical process with many astrophysical implications, the diffusion of cosmic rays (CRs) is determined by their interaction with magnetohydrodynamic (MHD) turbulence. We consider the magnetic mirroring effect arising from MHD turbulence on the diffusion of CRs. Due to the intrinsic superdiffusion of turbulent magnetic fields, CRs with large pitch angles that undergo mirror reflection, i.e., bouncing CRs, are not trapped between magnetic mirrors, but move diffusively along the magnetic field, leading to a new type of parallel diffusion. This diffusion is in general slower than the diffusion of non-bouncing CRs with small pitch angles that undergo gyroresonant scattering. The critical pitch angle at the balance between magnetic mirroring and pitch-angle scattering is important for determining the diffusion coefficients of both bouncing and non-bouncing CRs and their scalings with the CR energy. We find non-universal energy scalings of diffusion coefficients, depending on the properties of MHD turbulence.
Cosmic ray propagation is diffusive because of pitch angle scattering by waves. We demonstrate that if the high-amplitude magnetohydrodynamic turbulence with $tilde B/langle Brangle sim 1$ is present on top of the mean field gradient, the diffusion becomes asymmetric. As an example, we consider the vertical transport of cosmic rays in our Galaxy propagating away from a point-like source. We solve this diffusion problem analytically using a one-dimensional Markov chain analysis. We obtained that the cosmic ray density markedly differs from the standard diffusion prediction and has a sizable effect on their distribution throughout the galaxy. The equation for the continuous limit is also derived, which shows limitations of the convection-diffusion equation.
Fluid approximations to cosmic ray (CR) transport are often preferred to kinetic descriptions in studies of the dynamics of the interstellar medium (ISM) of galaxies, because they allow simpler analytical and numerical treatments. Magnetohydrodynamic (MHD) simulations of the ISM usually incorporate CR dynamics as an advection-diffusion equation for CR energy density, with anisotropic, magnetic field-aligned diffusion with the diffusive flux assumed to obey Ficks law. We compare test-particle and fluid simulations of CRs in a random magnetic field. We demonstrate that a non-Fickian prescription of CR diffusion, which corresponds to the telegraph equation for the CR energy density, can be easily calibrated to match the test particle simulations with great accuracy. In particular, we consider a random magnetic field in the fluid simulation that has a lower spatial resolution than that used in the particle simulation to demonstrate that an appropriate choice of the diffusion tensor can account effectively for the unresolved (subgrid) scales of the magnetic field. We show that the characteristic time which appears in the telegraph equation can be physically interpreted as the time required for the particles to reach a diffusive regime and we stress that the Fickian description of the CR fluid is unable to describe complex boundary or initial conditions for the CR energy flux.
Modelling of cosmic ray transport and interpretation of cosmic ray data ultimately rely on a solid understanding of the interactions of charged particles with turbulent magnetic fields. The paradigm over the last 50 years has been the so-called quasi-linear theory, despite some well-known issues. In the absence of a widely accepted extension of quasi-linear theory, wave-particle interactions must also be studied in numerical simulations where the equations of motion are directly solved in a realisation of the turbulent magnetic field. The applications of such test particle simulations of cosmic rays are manifold: testing transport theories, computing parameters like diffusion coefficients or making predictions for phenomena beyond standard diffusion theories, e.g. for cosmic ray small-scale anisotropies. In this review, we seek to give a low-level introduction to test particle simulations of cosmic rays, enabling readers to perform their own test particle simulations. We start with a review of quasi-linear theory, highlighting some of its issues and suggested extensions. Next, we summarise the state-of-the-art in test particle simulations and give concrete recipes for generating synthetic turbulence. We present a couple of examples for applications of such simulations and comment on an important conceptual detail in the backtracking of particles.
We investigate the shock acceleration of particles in massive galaxy mergers or collisions, and show that cosmic rays (CRs) can be accelerated up to the second knee energy ~0.1-1 EeV and possibly beyond, with a hard spectral index Gamma ~ 2. Such CRs lose their energy via hadronuclear interactions within a dynamical timescale of the merger shock, producing gamma rays and neutrinos as a by-product. If ~ 10 % of the shock dissipated energy goes into CR acceleration, some local merging galaxies will produce gamma-ray counterparts detectable by CTA. Also, based on the concordance cosmology, where a good fraction of the massive galaxies experience a major merger in a cosmological timescale, the neutrino counterparts can constitute ~ 20-60 % of the isotropic background detected by IceCube.