No Arabic abstract
Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic-ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficients dependencies on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely-used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in the present work. Consequently, a proper description of cosmic-ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic-ray gradient problem.
Cosmic-ray transport in astrophysical environments is often dominated by the diffusion of particles in a magnetic field composed of both a turbulent and a mean component. This process needs to be understood in order to properly model cosmic-ray signatures. One of the most important aspects in the modeling of cosmic-ray diffusion is that fully resonant scattering, the most effective such process, is only possible if the wave spectrum covers the entire range of propagation angles. By taking the wave spectrum boundaries into account, we quantify cosmic-ray diffusion parallel and perpendicular to the guide field direction at turbulence levels above 5% of the total magnetic field. We apply our results of the parallel and perpendicular diffusion coefficient to the Milky Way. We show that simple purely diffusive transport is in conflict with observations of the inner Galaxy, but that just by taking a Galactic wind into account, data can be matched in the central 5 kpc zone. Further comparison shows that the outer Galaxy at $>5$ kpc, on the other hand, should be dominated by perpendicular diffusion, likely changing to parallel diffusion at the outermost radii of the Milky Way.
We have performed magnetohydrodynamical simulations to study the amplification of magnetic fields in the precursors of shock waves. Strong magnetic fields are required in the precursors of the strong shocks that occur in supernova remnants. Observations also suggest that magnetic field amplification takes place in the weak shocks that occur in galaxy clusters and that produce so-called radio relics. Here, we extend the study of magnetic field amplification by cosmic-ray driven turbulence to weak shocks. The amplification is driven by turbulence that is produced by the cosmic-ray pressure acting on the density inhomogeneities in the upstream fluid. The clumping that has been inferred from X-ray data for the outskirts of galaxy clusters could provide some of the seed inhomogeneities. Magnetic field power spectra and Faraday maps are produced. Furthermore, we investigate how the synchrotron emission in the shock precursor can be used to verify the existence of this instability and constrain essential plasma parameters.
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.
This work has the main objective to provide a detailed investigation of cosmic ray propagation in magnetohydrodynamic turbulent fields generated by forcing the fluid velocity field at large scales. It provides a derivation of the particle mean free path dependences in terms of the turbulence level described by the Alfvenic Mach number and in terms of the particle rigidity. We use an upgrade version of the magnetohydrodynamic code {tt RAMSES} which includes a forcing module and a kinetic module and solve the Lorentz equation for each particle. The simulations are performed using a 3 dimension periodical box in the test-particle and magnetostatic limits. The forcing module is implemented using an Ornstein-Uhlenbeck process. An ensemble average over a large number of particle trajectories is applied to reconstruct the particle mean free paths. We derive the cosmic ray mean free paths in terms of the Alfvenic Mach numbers and particle reduced rigidities in different turbulence forcing geometries. The reduced particle rigidity is $rho=r_L/L$ where $r_L$ is the particle Larmor radius and $L$ is the simulation box length related to the turbulence coherence or injection scale $L_{inj}$ by $L sim 5 L_{inj}$. We have investigated with a special attention compressible and solenoidal forcing geometries. We find that compressible forcing solutions are compatible with the quasi-linear theory or more advanced non-linear theories which predict a rigidity dependence as $rho^{1/2}$ or $rho^{1/3}$. Solenoidal forcing solutions at least at low or moderate Alfvenic numbers are not compatible with the above theoretical expectations and require more refined arguments to be interpreted. It appears especially for Alfvenic Mach numbers close to one that the wandering of field lines controls perpendicular mean free path solutions whatever the forcing geometry.
Low-energy cosmic rays, in particular protons with energies below 1 GeV, are significant drivers of the thermochemistry of molecular clouds. However, these cosmic rays are also greatly impacted by energy losses and magnetic field transport effects in molecular gas. Explaining cosmic ray ionization rates of $10^{-16}$ s$^{-1}$ or greater in dense gas requires either a high external cosmic ray flux, or local sources of MeV-GeV cosmic ray protons. We present a new local source of low-energy cosmic rays in molecular clouds: first order Fermi-acceleration of protons in regions undergoing turbulent reconnection in molecular clouds. We show from energetic-based arguments there is sufficient energy within the magneto-hydrodynamic turbulent cascade to produce ionization rates compatible with inferred ionization rates in molecular clouds. As turbulent reconnection is a volume-filling process, the proposed mechanism can produce a near-homogeneous distribution of low-energy cosmic rays within molecular clouds.