Validating weak particle-field interaction theory for cosmic-ray parallel diffusion: higher-order 90-degree scattering corrections, empirical validity conditions, and the diffusive-to-coherent transition

This record audits the open questions left by the weak particle-field interaction (quasilinear) theory of cosmic-ray diffusion in random magnetic fields developed by Alpers, Hasselmann & Kunstmann (1978) and its companion paper Kunstmann & Alpers (1977). The theory expands the particle orbit in powers of the small random-field amplitude and yields a resonant pitch-angle diffusion coefficient $D$; under a two-timing approximation the parallel spatial diffusion coefficient $\kappa_{\parallel}$ follows from an integral over $D$ that becomes singular at 90 degrees pitch angle, where the lowest-order $D$ vanishes or develops a delta-function term. Three linked issues remain open in this program: whether real interplanetary magnetic irregularities satisfy the weak-interaction smallness and time-scale-separation conditions at resonant wavenumbers; whether the sketched higher-order correction $D^{(3)}$ produces a relaxation time near 90 degrees short enough to validate the two-timing approximation and recover a finite $\kappa_{\parallel}$; and the quantitative spectral criteria separating diffusive from coherent propagation when the lowest-order $\kappa_{\parallel}$ integral diverges - an open determination posed explicitly by Kunstmann & Alpers (1977) within the same joint program. Later literature resolved the idealized slab-model version of the higher-order question (second-order quasilinear theory plus test-particle simulations now describe 90-degree scattering) and established that realistic turbulence levels violate weak-interaction smallness. The surviving open core is the quantitative completion and observational confrontation of the program for realistic three-dimensional, dynamical, strong-amplitude turbulence.

Background

Energetic charged particles (galactic cosmic rays and solar energetic particles) propagate through the heliosphere and interstellar space in trajectories dominated by a large-scale, nearly uniform magnetic field $B_0$, with a turbulent random component $\delta B$ superposed. The pitch angle (cosine $\mu=v_{\parallel}/v$, with $v$ the particle speed) is the angle between the particle velocity and the mean field; the particle gyrates about the field with gyrofrequency $\Omega=e|B_0|/(\gamma mc)$ and is slowly scattered in pitch angle by the random component. The scattering arises from gyroresonance, the condition $k_{\parallel}v_{\parallel}+n\Omega=0$ (integer harmonic $n$) between the particle motion along the field and plane-wave components of the turbulence with field-parallel wavenumber $k_{\parallel}$. Quasilinear theory (QLT), introduced for cosmic-ray transport by Jokipii (1966) and developed by the Hamburg school as a weak particle-field interaction theory, expands the orbit in powers of the small parameter $\delta B/B_0$ and keeps the first non-vanishing term of the pitch-angle diffusion coefficient (Fokker-Planck coefficient) $D$, obtained from resonant wave contributions. Under a two-timing approximation (pitch-angle scattering fast compared with spatial transport), the parallel spatial diffusion coefficient follows from the Hasselmann-Wibberenz integral, which in the source notation reads $\kappa_{\parallel}=\frac{1}{8U}\int_{-U}^{+U}\frac{v_{\perp}^{4}}{D(v_{\parallel})}\,\mathrm{d}v_{\parallel}$ with $U$ the total particle speed (the same quantity as $v$ above) and $v_{\parallel}$, $v_{\perp}$ its components along and across the mean field; the source paper denotes the field-parallel component $U_{\parallel}$ when stating its resonance and validity conditions (Hasselmann & Wibberenz 1970).

At pitch angle 90 degrees ($v_{\parallel}=0$) this construction becomes singular: the resonance condition degenerates (the $n=0$ term resonates for all $k_{\parallel}$), Fisk et al. (1974) found a delta-function contribution to $D$ concentrated at 90 degrees, and for several turbulence spectra the lowest-order $D$ vanishes (or behaves pathologically) near $v_{\parallel}=0$, making the $\kappa_{\parallel}$ integral diverge. Whether this signals a breakdown of the weak-interaction framework was a live 1970s debate (Völk 1973; Jones, Kaiser & Birmingham 1973; Owens 1974; Goldstein, Klimas & Sandri 1975).

Alpers, Hasselmann & Kunstmann (1978, AHK) re-analyzed the question. They derived the complete lowest-order $D$ including the non-resonant $n=0$ contribution proportional to $\delta(U_{\parallel})$, argued this singular term is purely local and does not affect pitch-angle scattering or $\kappa_{\parallel}$ for averaged particle properties, and derived explicit conditions for the validity of the resonance asymptotics: the observation time must greatly exceed the gyroperiod ($t\gg\Omega^{-1}$) and the spectral-smoothness time $2\pi/(|U_{\parallel}|\,\delta k_{\parallel})$, where $\delta k_{\parallel}$ is the wavenumber scale over which the magnetic spectral tensor varies across the resonance surface $k_{\parallel}U_{\parallel}=-n\Omega$, while remaining below the perturbative bound $t\ll U^{2}/(2D)$. Near 90 degrees, where the lowest-order $D\to0$, AHK noted that an extension to higher order is necessary: under the Gaussian hypothesis for the field fluctuations the second-order term vanishes, and the leading correction is the third-order quantity $D^{(3)}=\frac{1}{2t}\langle 2v_{\parallel}^{(1)}v_{\parallel}^{(3)}+v_{\parallel}^{(2)}v_{\parallel}^{(2)}\rangle$, which they sketched but did not evaluate. The quantitative program that follows from this critical validity discussion - confronting the validity inequalities with measured interplanetary spectra, completing and evaluating the higher-order corrections, and mapping when transport is diffusive rather than coherent, that is, nearly scatter-free, ballistic-like propagation - remains to be carried out. The companion paper Kunstmann & Alpers (1977) developed the coherent-propagation side for power-law spectra and explicitly left open the determination of which mode of transport is present for a specific spectral index actually measured in interplanetary space.

Subsequent decades reframed the problem. Spacecraft measurements established that solar-wind turbulence at resonant scales is strong, $\delta B/B_0\sim O(1)$, and anisotropic, dominated by a composite of slab (field-aligned wavevectors) and two-dimensional (field-perpendicular wavevectors, magnetostatic) components. Nonlinear corrections to QLT were systematized - resonance broadening, in which the finite duration of the resonant wave-particle interaction replaces the sharp resonance condition by a smoothed kernel, and the nonlinear guiding-center and weakly nonlinear theories, which build such higher-order effects into closed transport equations - culminating in second-order quasilinear theory (SOQLT, the systematic second-order extension of QLT; Shalchi 2005) and the monograph of Shalchi (2009). Test-particle simulations became the arbiter of the 90-degree behavior (Qin & Shalchi 2009, 2014; Tautz et al. 2013; van den Berg et al. 2024). The 90-degree scattering problem in slab turbulence is now generally regarded as resolved by SOQLT, while realistic three-dimensional, dynamical, imbalanced (with unequal fluxes of oppositely propagating wave modes) and strong-amplitude turbulence remains an active frontier.

Problem Statement

The problem left open by the joint program of Alpers, Hasselmann & Kunstmann (1978) and Kunstmann & Alpers (1977) is to complete the quantitative validation program of the weak particle-field interaction theory of cosmic-ray parallel diffusion in random magnetic fields. Its three components share one scientific objective - a quantitatively validated theory of the parallel diffusion of energetic charged particles, with validity conditions expressed in terms of observable turbulence properties. The task therefore comprises three linked components (1)-(3), which jointly define the research target: the empirical validity boundary (1), the higher-order 90-degree completion (2), and the diffusive-to-non-diffusive regime criteria (3), including the specific spectral-index determination nested in component 3.

1. Empirical validity: measured interplanetary magnetic-field irregularities are already known to violate the weak-interaction smallness condition (small fluctuation amplitude $\delta B/B_0$) at the resonant wavenumbers; given that violation, quantify whether and where their spectral shapes at those wavenumbers satisfy the time-scale-separation conditions $t \gg \Omega^{-1}$, $t \gg 2\pi/(|U_{\parallel}|\,\delta k_{\parallel})$ and $t \ll U^{2}/(2D)$, equivalently mapping the validity boundary of the quasilinear description in terms of measured spectral parameters and turbulence levels and quantifying by how much observed amplitudes and spectral shapes at the resonant wavenumbers lie outside that boundary.

2. Higher-order completion: carry out the higher-order perturbation calculation of the pitch-angle diffusion coefficient near 90 degrees pitch angle (the third-order term $D^{(3)}=\frac{1}{2t}\langle 2v_{\parallel}^{(1)}v_{\parallel}^{(3)}+v_{\parallel}^{(2)}v_{\parallel}^{(2)}\rangle$ in the source expansion, or a higher-order/nonlinear correction demonstrably equivalent to it in the sense of carrying the same third-order contribution to $D$ near $U_{\parallel}=0$) in full, and determine quantitatively whether the finite relaxation time it produces near $U_{\parallel}=0$ is short enough, under realistic solar-wind spectra and fluctuation amplitudes, to validate the two-timing approximation and recover a finite parallel diffusion coefficient $\kappa_{\parallel}$.

3. Regime criteria: when the lowest-order integral for $\kappa_{\parallel}$ diverges, establish the quantitative transition criteria - spectral thresholds at the resonant wavenumbers, the role of finite resonance widths, turbulence geometry and dynamics - separating the diffusive regime from the non-diffusive regimes - coherent (nearly scatter-free, ballistic) propagation, subdiffusive transport, and superdiffusive transport - for observed magnetic spectra (the determination left open by Kunstmann & Alpers 1977: which mode of transport is present for a specific spectral index actually measured in interplanetary space).

References

1. Alpers, W., Hasselmann, K. & Kunstmann, J. (1978). On the validity of weak particle-field interaction theory for the description of cosmic-ray particle diffusion in random magnetic fields. Astrophysics and Space Science 58, 259-271. DOI: 10.1007/BF00644516. https://doi.org/10.1007/BF00644516
2. Jokipii, J. R. (1966). Cosmic-Ray Propagation. I. Charged Particles in a Random Magnetic Field. The Astrophysical Journal 146, 480. DOI: 10.1086/148912. https://doi.org/10.1086/148912
3. Hasselmann, K. & Wibberenz, G. (1970). A Note on the Parallel Diffusion Coefficient. The Astrophysical Journal 162, 1049. DOI: 10.1086/150736. https://doi.org/10.1086/150736
4. Jones, F. C., Kaiser, T. B. & Birmingham, T. J. (1973). New Approach to Cosmic-Ray Diffusion Theory. Physical Review Letters 31, 485-487. DOI: 10.1103/PhysRevLett.31.485. https://doi.org/10.1103/PhysRevLett.31.485
5. Völk, H. J. (1973). Nonlinear perturbation theory for cosmic ray propagation in random magnetic fields. Astrophysics and Space Science 25, 471-490. DOI: 10.1007/BF00649186. https://doi.org/10.1007/BF00649186
6. Fisk, L. A., Goldstein, M. L., Klimas, A. J. & Sandri, G. (1974). The Fokker-Planck Coefficient for Pitch-Angle Scattering of Cosmic Rays. The Astrophysical Journal 190, 417. DOI: 10.1086/152893. https://doi.org/10.1086/152893
7. Owens, A. J. (1974). The Effects of Nonlinear Terms in Cosmic-Ray Diffusion Theory. The Astrophysical Journal 191, 235. DOI: 10.1086/152960. https://doi.org/10.1086/152960
8. Goldstein, M. L., Klimas, A. J. & Sandri, G. (1975). Mirroring in the Fokker-Planck coefficient for cosmic-ray pitch-angle scattering in homogeneous magnetic turbulence. The Astrophysical Journal 195, 787. DOI: 10.1086/153383. https://doi.org/10.1086/153383
9. Kunstmann, J. E. & Alpers, W. (1977). A perturbation approach to coherent propagation of energetic charged particles in random magnetic fields. The Astrophysical Journal 211, 587. DOI: 10.1086/154967. https://doi.org/10.1086/154967
10. Bieber, J. W., Earl, J. A., Green, G., Kunow, H., Müller-Mellin, R. & Wibberenz, G. (1980). Interplanetary pitch angle scattering and coronal transport of solar energetic particles: New information from Helios. Journal of Geophysical Research 85, 2313-2323. DOI: 10.1029/JA085IA05P02313. https://doi.org/10.1029/JA085IA05P02313
11. Bieber, J. W., Smith, C. W. & Matthaeus, W. H. (1988). Cosmic-ray pitch-angle scattering in isotropic turbulence. The Astrophysical Journal 334, 470. DOI: 10.1086/166851. https://doi.org/10.1086/166851
12. Felice, G. M. & Kulsrud, R. M. (2001). Cosmic-Ray Pitch-Angle Scattering through 90°. The Astrophysical Journal 553, 198-210. DOI: 10.1086/320651. https://doi.org/10.1086/320651
13. Shalchi, A. (2005). Second-order quasilinear theory of cosmic ray transport. Physics of Plasmas 12, 052905. DOI: 10.1063/1.1895805. https://doi.org/10.1063/1.1895805
14. Shalchi, A., Bieber, J. W. & Matthaeus, W. H. (2008). Pitch-angle scattering in pure two-dimensional and two-component turbulence. Astronomy & Astrophysics 483, 371-381. DOI: 10.1051/0004-6361:20079099. https://doi.org/10.1051/0004-6361:20079099
15. Shalchi, A., Škoda, T., Tautz, R. C. & Schlickeiser, R. (2009). Analytical description of nonlinear cosmic ray scattering: isotropic and quasilinear regimes of pitch-angle diffusion. Astronomy & Astrophysics 507, 589-597. DOI: 10.1051/0004-6361:200912755. https://doi.org/10.1051/0004-6361:200912755
16. Qin, G. & Shalchi, A. (2009). Pitch-angle diffusion coefficients of charged particles from computer simulations. The Astrophysical Journal 707, 61-66. DOI: 10.1088/0004-637X/707/1/61. https://doi.org/10.1088/0004-637X/707/1/61
17. Shalchi, A. (2009). Nonlinear Cosmic Ray Diffusion Theories. Astrophysics and Space Science Library, Springer. DOI: 10.1007/978-3-642-00309-7. https://doi.org/10.1007/978-3-642-00309-7
18. Shalchi, A., Webb, G. M. & le Roux, J. A. (2012). Parallel transport of cosmic rays for non-diffusive pitch-angle scattering: I. Using the standard Fokker-Planck equation. Physica Scripta 85, 065901. DOI: 10.1088/0031-8949/85/06/065901. https://doi.org/10.1088/0031-8949/85/06/065901
19. Mace, R. L., Dalena, S. & Matthaeus, W. H. (2012). Velocity space diffusion of charged particles in weak magnetostatic fields: Nonlinear effects, model constraints, and implications for simulations. Physics of Plasmas 19, 032309. DOI: 10.1063/1.3693379. https://doi.org/10.1063/1.3693379
20. Tautz, R. C., Dosch, A., Effenberger, F., Fichtner, H. & Kopp, A. (2013). Pitch-angle scattering in magnetostatic turbulence. Astronomy & Astrophysics 558, A147. DOI: 10.1051/0004-6361/201322142. https://doi.org/10.1051/0004-6361/201322142
21. Qin, G. & Shalchi, A. (2014). Detailed numerical investigation of 90° scattering of energetic particles interacting with magnetic turbulence. Physics of Plasmas 21, 042906. DOI: 10.1063/1.4873895. https://doi.org/10.1063/1.4873895
22. Weidl, M. S., Jenko, F., Teaca, B. & Schlickeiser, R. (2015). Cosmic-ray pitch-angle scattering in imbalanced MHD turbulence simulations. The Astrophysical Journal 811, 8. DOI: 10.1088/0004-637X/811/1/8. https://doi.org/10.1088/0004-637X/811/1/8
23. Ostryakov, V. M. (2015). How to measure pitch-angle diffusion coefficient at θ ~ 90°. Geomagnetism and Aeronomy 55, 976-978. DOI: 10.1134/S001679321507018X. https://doi.org/10.1134/S001679321507018X
24. Reichherzer, P., Becker Tjus, J., Zweibel, E. G., Merten, L. & Pueschel, M. J. (2020). Turbulence-level dependence of cosmic ray parallel diffusion. Monthly Notices of the Royal Astronomical Society 498, 5051-5064. DOI: 10.1093/mnras/staa2533. https://doi.org/10.1093/mnras/staa2533
25. van den Berg, J. P., Els, P. L. & Engelbrecht, N. E. (2024). An Evaluation of Different Numerical Methods to Calculate the Pitch-angle Diffusion Coefficient from Full-orbit Simulations: Disentangling a Rope of Sand. The Astrophysical Journal 977, 174. DOI: 10.3847/1538-4357/ad8b4a. https://doi.org/10.3847/1538-4357/ad8b4a
26. Subashchandar, N. S. M., Zhao, L., Shalchi, A., Zank, G. P., le Roux, J., Li, H., Zhu, X., Silwal, A. & Alonso Guzman, J. G. (2025). Parallel and perpendicular diffusion of energetic particles in the near-Sun solar wind observed by Parker Solar Probe. Accepted for publication in The Astrophysical Journal Letters; arXiv:2509.10648. https://arxiv.org/abs/2509.10648
27. Yan, H., Zhao, S. & Zhang, M. (2026). Pitch-Angle Scattering of Cosmic Rays: Confronting Theory with Observations. Accepted in The Astrophysical Journal Letters; arXiv:2605.24143. https://arxiv.org/abs/2605.24143