Focused Transport Equation Overview

• The focused transport equation is a kinetic framework that models charged particle transport by coupling systematic magnetic focusing with stochastic pitch-angle diffusion.
• Fourier–Legendre projections and SDE methods offer efficient, norm-conserving solutions that capture ballistic, diffusive, and drift regimes in particle propagation.
• Its applications include solar energetic particle events and cosmic-ray propagation, though limitations remain in handling perpendicular diffusion and energy losses.

The focused transport equation (FTE) is a fundamental kinetic framework for describing the spatial and momentum evolution of charged particles (notably solar energetic particles and cosmic rays) as they move through spatially varying magnetic fields and undergo pitch-angle scattering by turbulence. Its central feature is the coupling of deterministic focusing effects—where particles experience systematic drift due to magnetic field inhomogeneities—with stochastic pitch-angle diffusion. The FTE underpins quantitative models of cosmic-ray propagation, solar energetic particle events, and related space- and astrophysical transport phenomena (Berg et al., 2020, Klippenstein et al., 15 Jul 2025, Merten et al., 2024).

1. Mathematical Formulation and Key Physical Content

The canonical form of the focused transport equation for the phase-space density $f(z, \mu, t)$ is

$$\frac{\partial f}{\partial t} + v\,\mu\,\frac{\partial f}{\partial z} = \frac{\partial}{\partial \mu}\left[ D_{\mu\mu}(\mu)\,\frac{\partial f}{\partial \mu} \right] - \frac{v}{2L} (1-\mu^2)\frac{\partial f}{\partial \mu},$$

where $z$ is the field-aligned coordinate, $\mu = \cos\theta$ is the pitch-angle cosine, $v$ is the particle speed, $D_{\mu\mu}$ is the pitch-angle diffusion coefficient (commonly $D_{\mu\mu} = D_0 (1-\mu^2)$), and $L = -B_0/(dB_0/dz)$ is the magnetic focusing length (Klippenstein et al., 14 Apr 2025, Berg et al., 2020).

In the conservative ("modified" or norm-conserving) form, the variable change $\tilde f = e^{z/L} f$ leads to

$$\frac{\partial \tilde f}{\partial t} + \frac{\partial}{\partial z}(v\mu\,\tilde f) + \frac{\partial}{\partial \mu}\left(\frac{1-\mu^2}{2L}\,v\,\tilde f\right) = \frac{\partial}{\partial \mu}\left[ D_{\mu\mu}(\mu)\,\frac{\partial \tilde f}{\partial \mu} \right],$$

which, unlike the standard form, guarantees exact conservation of the particle number within a magnetic flux tube (Klippenstein et al., 15 Jul 2025, Merten et al., 2024).

The two fundamental FTE forms are:

Form	Norm Conservation	Pitch-angle Isotropization in Operator	Drift Structure
Standard	No	Explicit in scattering	Additive in $\mu$
Modified	Yes	Only through scattering term	Divergence in $\mu$

Conservation in the modified form follows from integrating over $z$ and $\mu$, reducing all divergences to boundary terms: $$\frac{d}{dt}\int dz\,d\mu\,\tilde f(t, z, \mu) = 0.$$ (Klippenstein et al., 15 Jul 2025)

2. Physical Interpretation and Transport Regimes

The FTE encodes the following physics:

• Streaming: $v\mu\partial_z f$ describes deterministic propagation along magnetic field lines.
• Focusing: The $L$-dependent term induces drift in pitch angle, tending to align (or anti-align) particle velocities with the local field direction depending on field convergence/divergence. The absence of this term ($L\to\infty$) recovers pure pitch-angle scattering.
• Pitch-angle diffusion: $D_{\mu\mu}(\mu)$ drives the particle distribution toward isotropy, counteracting anisotropies generated by streaming and focusing.
• Norm Conservation and Isotropization: In the conservative (modified) form, all focusing is encoded as $\mu$-space flux, ensuring norm conservation; however, isotropization is enforced solely by diffusion, not by the focusing drift (Klippenstein et al., 15 Jul 2025).

Focusing induces systematic spatial drift ("coherent streaming") and modifies the effective parallel diffusion coefficient and propagation speed. In weak focusing, the $\mu$-averaged solution exhibits classic diffusive-convective behavior with drift velocity $v_c = \kappa/L$ and suppressed diffusion $\bar\kappa = \kappa (1-\alpha\lambda^2/L^2)$ (Klippenstein et al., 14 Apr 2025, Berg et al., 2020).

3. Analytical and Numerical Solution Strategies

Subspace Approximation (Fourier–Legendre Method)

A key modern approach is the finite Legendre subspace (Fourier–Legendre) projection:

• Fourier-transform in $z$: $f(z,\mu,t) = \int dk\,F_k(\mu, t)\,e^{ikz}$.
• Legendre expansion in $\mu$: $F_k(\mu, t) = \sum_{n=0}^N C_n(t)P_n(\mu)$.

The FTE reduces to a finite-dimensional ODE system: $\dot{\mathbf{C}} = \mathbf{M}\,\mathbf{C},$ where $\mathbf{M}$ depends on $v, D, L, k$ and couples adjacent Legendre modes. For $N=2$, closed-form solutions are possible and capture ballistic, diffusive, and drift regimes analytically (Klippenstein et al., 14 Apr 2025, Klippenstein et al., 15 Jul 2025). For $N\sim10$, numerical evaluation of $\exp(\mathbf{M}t)$ yields rapid, highly accurate, norm-conserving solutions whose accuracy and speed surpass full grid-based solvers by orders of magnitude, particularly for large parameter scans (Klippenstein et al., 15 Jul 2025).

Stochastic Differential Equation (SDE) Methods

The SDE formalism maps the FTE (in Fokker–Planck form) to coupled Itô SDEs for spatial position and pitch angle: $$\begin{aligned} ds &= v\,\mu\,dt, \ d\mu &= \left[\frac{v}{2L}(1-\mu^2) + \frac{\partial D_{\mu\mu}}{\partial \mu}\right]dt + \sqrt{2D_{\mu\mu}(\mu)}\,dW_t. \end{aligned}$$ This method, implemented in CRPropa and similar frameworks, enables flexible treatments of boundary conditions, source terms, and ensemble statistics for cosmic-ray transport in arbitrary magnetic configurations and turbulent fields (Merten et al., 2024).

Finite-Difference Grid Methods

Explicit and implicit finite-difference solvers, as exemplified by van den Berg et al., discretize $f(z,\mu,t)$ on uniform grids in $z$ and $\mu$, handling streaming, focusing, and scattering by flux-conservative differencing and robust time integration (semi-implicit for the stiff diffusive term). These approaches are generally slower but handle complex time-dependent sources and arbitrary $L(z)$ (Berg et al., 2020).

4. Approximate Reductions, Boundary Conditions, and Limiting Forms

Diffusion–Advection and Telegraph Approximations

For weak anisotropy and smooth $L$, integrating over $\mu$ and expanding in Legendre polynomials yields the diffusion–advection equation for the isotropic density $F_0$: $$\frac{\partial F_0}{\partial t} + v_c\frac{\partial F_0}{\partial z} = \bar\kappa\frac{\partial^2 F_0}{\partial z^2},$$ with effective $\bar\kappa$ and $v_c$ as above (Klippenstein et al., 14 Apr 2025, Berg et al., 2020).

To enforce finite signal speed and improve early-time accuracy, the telegraph approximation includes a second-order time derivative: $$\frac{\partial F_0}{\partial t} + \tau\frac{\partial^2 F_0}{\partial t^2} = \kappa_\parallel \frac{\partial^2 F_0}{\partial z^2} + v_c\frac{\partial F_0}{\partial z},$$ where $\tau$ is the characteristic scattering time (Litvinenko et al., 2015).

Boundary Conditions

• Reflecting: $\partial_z F_0 = 0$ at wall.
• Absorbing: Mixed condition involving both $F_0$ and its time derivative at the boundary, reflecting the telegraph equation's wave nature. These boundary prescriptions resolve the otherwise non-vanishing edge density characteristic of the telegraph model, ensuring physical consistency for cosmic-ray escape, SEPs trapped by mirrors, and similar contexts (Litvinenko et al., 2015).

5. Computational Efficiency and Physical Fidelity

Hybrid semi-analytic/numerical subspace strategies reach sub-percent level accuracy for moments and characteristic functions with moderate truncation ($N\sim 10$). The matrix-exponential approach offers order-of-magnitude speedup over (z,μ,t) explicit grid solvers, enabling large-parameter surveys and resolving the ballistic-to-diffusive regime transition accurately even for strong focusing ($v/(DL)\lesssim 2$) (Klippenstein et al., 15 Jul 2025, Klippenstein et al., 14 Apr 2025).

The modified (conservative) FTE guarantees exact norm conservation at the solver algorithmic level, with all physical observables (moments, characteristic functions) obtained by direct quadrature and inversion.

6. Applications and Limitations in Astrophysical Contexts

The FTE and its approximations are foundational for simulations of:

• Solar energetic particle events: Accurate modeling along field-aligned Parker spirals with time-dependent injection and realistic boundary conditions (Berg et al., 2020).
• Galactic cosmic-ray escape: Focusing, boundaries, and coherent streaming all play critical roles for particle densities and time profiles in the Milky Way halo (Litvinenko et al., 2015).
• Shock acceleration and magnetic mirroring: Telegraphed and focused forms capture transient intensities, anisotropies, and non-diffusive particle propagation scenarios.

Limitations include the neglect of perpendicular diffusion, momentum diffusion, energy losses, and more complex magnetic topologies (except where SDEs allow easy generalization). The diffusion–advection and telegraph reductions fail during highly anisotropic early phases and strong focusing. For non-uniform focusing ($L(z)$), the full FTE or at least high-N subspace methods are required for accurate quantitative predictions (Berg et al., 2020, Merten et al., 2024).

7. Summary and Outlook

The focused transport equation, in both standard and conservative forms, remains the rigorous basis for kinetic modeling of energetic particle transport in heliospheric and astrophysical environments. Recent advances highlight the power of subspace (Fourier–Legendre) methods and SDE-based solvers for combining computational efficiency with norm conservation and physical fidelity. The FTE captures the interplay of magnetic focusing, pitch-angle diffusion, and ballistic transport, with broad relevance to cosmic-ray modulation, space weather forecasting, and fundamental plasma astrophysics (Klippenstein et al., 15 Jul 2025, Klippenstein et al., 14 Apr 2025, Merten et al., 2024, Berg et al., 2020).