The nature of cosmic ray (CR) transport in the Milky Way remains elusive. The predictions of current microphysical CR transport models in magnetohydrodynamic (MHD) turbulence are drastically different from what is observed. These models usually focus on MHD turbulence with a strong guide field and ignore the impact of turbulent intermittency on particle propagation. This motivates our studying the alternative regime of large-amplitude turbulence with δB/B0 ≫ 1, in which intermittent small-scale magnetic field reversals are ubiquitous. We study particle transport in such turbulence by integrating trajectories in stationary snapshots. To quantify spatial diffusion, we use a set-up with continuous particle injection and escape, which we term the turbulent leaky box. We find that particle transport is very different from the strong guide-field case. Low-energy particles are better confined than high-energy particles, despite less efficient pitch-angle isotropization at small energies. In the limit of weak guide field, energy-dependent confinement is driven by the energy-dependent (in)ability to follow reversing magnetic field lines exactly and by the scattering in regions of ‘resonant curvature’, where the field line bends on a scale that is of the order of the local particle gyro-radius. We derive a heuristic model of particle transport in magnetic folds that approximately reproduces the energy dependence of transport found numerically. We speculate that CR propagation in the Galaxy is regulated by the intermittent field reversals highlighted here and discuss the implications of our findings for CR transport in the Milky Way.more » « less
- NSF-PAR ID:
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Page Range / eLocation ID:
- p. 4985-4998
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
Abstract 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, which is two-fold turbulent mixing in that the particle motion is stochastic with respect to the field lines, 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$$ > 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.more » « less
Phenomenological models of cosmic ray (CR) transport in the Milky Way can reproduce a wide range of observations assuming that CRs scatter off of magnetic-field fluctuations with spectrum ∝ k−δ and δ ∼ [1.4, 1.67]. We study the extent to which such models can be reconciled with current microphysical theories of CR transport, specifically self-confinement due to the streaming instability and/or extrinsic turbulence due to a cascade of magnetohydrodynamic (MHD) fast modes. We first review why it is that on their own neither theory is compatible with observations. We then highlight that CR transport is a strong function of local plasma conditions in the multiphase interstellar medium, and may be diffusive due to turbulence in some regions and streaming due to self-confinement in others. A multiphase combination of scattering mechanisms can in principle reproduce the main trends in the proton spectrum and the boron-to-carbon ratio. However, models with a combination of scattering by self-excited waves and fast-mode turbulence require significant fine-tuning due to fast-mode damping, unlike phenomenological models that assume undamped Kolmogorov turbulence. The assumption that fast modes follow a weak cascade is also not well justified theoretically, as the weak cascade is suppressed by wave steepening and weak-shock dissipation even in subsonic turbulence. These issues suggest that there may be a significant theoretical gap in our understanding of MHD turbulence. We discuss a few topics at the frontier of MHD turbulence theory that bear on this (possible) gap and that may be relevant for CR scattering.
Subsonic, compressive turbulence transfers energy to cosmic rays (CRs), a process known as nonresonant reacceleration. It is often invoked to explain the observed ratios of primary to secondary CRs at ∼GeV energies, assuming wholly diffusive CR transport. However, such estimates ignore the impact of CR self-confinement and streaming. We study these issues in stirring box magnetohydrodynamic (MHD) simulations using Athena++, with field-aligned diffusive and streaming CR transport. For diffusion only, we find CR reacceleration rates in good agreement with analytic predictions. When streaming is included, reacceleration rates depend on plasma
β. Due to streaming-modified phase shifts between CR and gas variables, they are slower than canonical reacceleration rates in low- βenvironments like the interstellar medium but remain unchanged in high- βenvironments like the intracluster medium. We also quantify the streaming energy-loss rate in our simulations. For sub-Alfvénic turbulence, it is resolution dependent (hence unconverged in large-scale simulations) and heavily suppressed compared to the isotropic loss rate vA· ∇ PCR/ PCR∼ vA/ L0, due to misalignment between the mean field and isotropic CR gradients. Unlike acceleration efficiencies, CR losses are almost independent of magnetic field strength over β∼ 1–100 and are, therefore, not the primary factor behind lower acceleration rates when streaming is included. While this paper is primarily concerned with how turbulence affects CRs, in a follow-up paper we consider how CRs affect turbulence by diverting energy from the MHD cascade, altering the pathway to gas heating and steepening the turbulent spectrum.
The coarse-grained propagation of galactic cosmic rays (CRs) is traditionally constrained by phenomenological models of Milky Way CR propagation fit to a variety of direct and indirect observables; however, constraining the fine-grained transport of CRs along individual magnetic field lines—for instance, diffusive vs streaming transport models—is an unsolved challenge. Leveraging a recent training set of magnetohydrodynamic turbulent box simulations, with CRs spanning a range of transport parameters, we use convolutional neural networks (CNNs) trained solely on gas density maps to classify CR transport regimes. We find that even relatively simple CNNs can quite effectively classify density slices to corresponding CR transport parameters, distinguishing between streaming and diffusive transport, as well as magnitude of diffusivity, with class accuracies between 92% and 99%. As we show, the transport-dependent imprints that CRs leave on the gas are not all tied to the resulting density power spectra: classification accuracies are still high even when image spectra are flattened (85%–98% accuracy), highlighting CR transport-dependent changes to turbulent phase information. We interpret our results with saliency maps and image modifications, and we discuss physical insights and future applications.
Direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers ( 𝑆𝑡=0.01,0.1,0.5,1,2 ) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.more » « less