Cosmic-ray transport in inhomogeneous media

A theory of cosmic-ray transport in multi-phase diffusive media is developed, with the specific application to cases in which the cosmic-ray diffusion coefficient has large spatial fluctuations that may be inherently multi-scale. We demonstrate that the resulting transport of cosmic rays is diffusive in the long-time limit, with an average diffusion coefficient equal to the harmonic mean of the spatially varying diffusion coefficient. Thus, cosmic-ray transport is dominated by areas of low diffusion even if these areas occupy a relatively small, but not infinitesimal, fraction of the volume. On intermediate time scales, the cosmic rays experience transient effective sub-diffusion, as a result of low-diffusion regions interrupting long flights through high-diffusion regions. In the simplified case of a two-phase medium, we show that the extent and extremity of the sub-diffusivity of cosmic-ray transport is controlled by the spectral exponent of the distribution of patch sizes of each of the phases. We finally show that, despite strongly influencing the confinement times, the multi-phase medium is only capable of altering the energy dependence of cosmic-ray transport when there is a moderate (but not excessive) level of perpendicular diffusion across magnetic-field lines.

💡 Research Summary

This paper develops a comprehensive theory of cosmic‑ray (CR) transport in a multi‑phase, highly inhomogeneous medium where the parallel diffusion coefficient varies dramatically from place to place. The authors model CR propagation along magnetic‑field lines as a one‑dimensional diffusion process with a spatially varying coefficient κ(x). They first show that on very short timescales a CR only “feels” the diffusion coefficient of its local patch, so the ensemble‑averaged diffusion coefficient is simply the arithmetic mean κ_A of κ(x). As time progresses and particles sample many patches, the transport becomes governed by the harmonic mean κ_H of the local coefficients. This result is derived analytically (Appendix B) and can be understood physically as the scattering rates (∝ κ⁻¹) adding in series, analogous to resistors in an electric circuit.

The authors then focus on a two‑phase model in which κ(x) takes only two values, κ_high and κ_low, arranged in patches whose lengths follow a power‑law distribution P(l)∝l⁻ᵅ. Using Monte‑Carlo random‑walk simulations (both in real space and in pitch‑angle space) together with a semi‑analytic scaling argument, they demonstrate a transient sub‑diffusive regime at intermediate times. During this regime the “running diffusion coefficient” \bar{κ}(t) decays from κ_A toward κ_H. The decay law depends on the spectral exponent α: for α > 2 the decay is logarithmic, while for α < 2 it follows a power law \bar{κ}(t)∝t^{-(2‑α)}. The authors provide a quantitative prediction for the decay exponent (Appendix C) that matches the simulations very well (Fig. 3). The key physical picture is that low‑diffusion patches act as bottlenecks, interrupting long flights through high‑diffusion regions and thereby slowing the overall spread.

In the limit of many patches the long‑time diffusion coefficient reduces to a simple expression involving only the filling fraction f of the low‑diffusion phase: κ_A≈(1‑f)κ_high and κ_H≈κ_low/f. Thus, detailed knowledge of the full turbulence spectrum is unnecessary for the mean‑field description; only the two phase diffusivities and their volume fractions matter.

The paper also addresses the energy dependence of CR transport. If diffusion is purely parallel (κ_⊥ = 0), the multi‑phase structure does not alter the energy scaling of κ, because every particle must cross the low‑diffusion patches regardless of energy. However, when a modest perpendicular diffusion κ_⊥ is allowed, CRs can wander around low‑diffusion islands. In this case the effective transport becomes sensitive to the energy dependence of both κ_∥(E) and κ_⊥(E). The authors argue that only a moderate level of perpendicular diffusion (neither vanishingly small nor overwhelmingly large) can imprint the patchy structure onto the overall energy dependence of the escape rate.

Mathematical details are relegated to the appendices: Appendix A proves the short‑time arithmetic‑mean result, Appendix B derives the harmonic‑mean limit, Appendix C develops the sub‑diffusive scaling theory, and Appendix D describes the Monte‑Carlo implementation. The authors conclude by emphasizing that their framework provides a clear “early‑time → sub‑diffusive → late‑time” picture of CR transport in realistic, multi‑phase astrophysical plasmas. They suggest future work to (i) infer the filling fraction and patch‑size exponent from observations of CR spectra and γ‑ray emission, (ii) extend the theory to fully three‑dimensional magnetic topologies, and (iii) explore regimes of strong perpendicular diffusion driven by plasma instabilities. Overall, the study offers a robust, analytically grounded mean‑field model that captures how low‑diffusion regions dominate CR confinement and how the statistical properties of the medium control transient transport behavior.