Electron Parallel Nonlinearity (EPN)

• Electron Parallel Nonlinearity (EPN) is the nonlinear fluctuation-induced advection of the electron distribution along the parallel direction in gyrokinetic and drift-kinetic models.
• It plays a key role in electromagnetic turbulence by triggering early saturation mechanisms and significantly reducing heat transport in kinetic ballooning regimes.
• EPN connects diverse phenomena—from reconnection heating and Langmuir wave trapping to nonlocal electron transport closures and hot‐electron beam coupling—by emphasizing electron dynamics along a privileged propagation direction.

Electron Parallel Nonlinearity (EPN) most explicitly denotes the nonlinear fluctuation-induced advection of the electron distribution in parallel phase space, represented in gyrokinetics by terms such as $$\dot v_{\parallel 1}\,\partial \delta f/\partial v_\parallel$$ or $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$ (Chen et al., 25 Sep 2025). Several adjacent literatures do not use the label explicitly, yet treat closely related electron-mediated phenomena: trapped-electron modification of Langmuir-wave response (Strozzi et al., 2012), heating and distribution reshaping by field-aligned electric fields in reconnection (Richard et al., 2024), nonlocal parallel response operators for heat flow, friction, and viscosity (Ji et al., 2019, Ji et al., 2019), and hot-electron two-beam coupling through interference-induced temperature gratings (Paul et al., 2020). This suggests that EPN is best understood as a family of nonlinear or nonlocal electron responses organized by motion, acceleration, or transport along a distinguished parallel direction.

1. Formal gyrokinetic definition

In the most direct usage, EPN is a specific nonlinear term in the electron drift-kinetic or gyrokinetic equation. In the parallel-velocity formulation, the perturbed distribution satisfies

$$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$

and the EPN is the nonlinear contribution

$$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$

In GEM, which uses the parallel canonical momentum

$$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$

the corresponding term is

$$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$

Physically, this is nonlinear coupling between fluctuation-induced parallel acceleration and fine structure in electron parallel velocity space (Chen et al., 25 Sep 2025).

This term is usually dropped under the Frieman–Chen ordering because it is formally $O(\delta^3)$, whereas the standard retained nonlinearities are $O(\delta^2)$. The 2025 GEM study is important because it shows that formal ordering smallness does not guarantee dynamical insignificance in electromagnetic turbulence. The same paper also makes clear that implementation is not limited to modifying the particle equation of motion: in a split-weight PIC formulation, EPN-consistent higher-order terms must also be retained in the electron weight equation and in the vorticity equation. A partially linearized treatment neglects precisely these higher-order contributions, including the distinction between $p_\parallel$ and $v_\parallel$ in the nonlinear terms (Chen et al., 25 Sep 2025).

A common misconception is that EPN is simply another name for parallel conductivity or parallel streaming. The gyrokinetic definition is narrower: it is a nonlinear phase-space advection term, not a constitutive coefficient. That distinction becomes important when comparing EPN to nonlocal closure models, which encode parallel electron response but do not evolve the full nonlinear electron distribution.

2. Electromagnetic turbulence and kinetic ballooning saturation

The clearest demonstration of EPN as a turbulence-saturation mechanism appears in the GEM study of the Cyclone Base Case above the kinetic ballooning mode threshold. For the global electromagnetic setup, the chosen primary mode is $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$0, corresponding to

$$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$1

and a linear $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$2 scan shows transition from ITG to KBM at about

$$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$3

At $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$4, the linear $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$5 mode has

$$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$6

The key nonlinear result is that a simulation without full EPN undergoes runaway early in the nonlinear phase, whereas the simulation including EPN saturates at a low initial amplitude and exhibits a strong reduction of heat transport relative to the runaway case (Chen et al., 25 Sep 2025).

The onset of nonlinear regulation is unusually early. With EPN, deviation from linear growth occurs at

$$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$7

whereas in the electrostatic ITG comparison the corresponding deviation occurs only near

$$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$8

The same study states that the EPN-inclusive KBM case saturates at an initial amplitude two orders of magnitude below that of the electrostatic case. By contrast, in the electrostatic regime at $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$9, including EPN does not substantially alter the usual zonal-flow-mediated saturation picture. The paper therefore argues that EPN is not generically strong in all turbulence, but becomes decisive in electromagnetic KBM regimes (Chen et al., 25 Sep 2025).

The proposed mechanism depends on magnetic flutter. For the $$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$0 KBM, the linear simulation gives

$$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$1

which might seem modest, but for electrons

$$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$2

Thus electron radial displacement is dominated by magnetic flutter rather than by $$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$3 drift. The maximum field-line width near the mode location is reported as

$$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$4

while the characteristic radial fine-structure scale of the linear eigenmode is estimated as

$$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$5

Because the flutter-induced radial displacement depends on $$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$6, radial fine structure is converted into velocity-space fine structure. The paper expresses this by differentiating the phase-space constancy relation

$$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$7

with respect to $$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$8, producing a term proportional to

$$\frac{\partial \delta f}{\partial t} + \mathbf{v}_G \cdot \nabla \delta f + \dot v_\parallel \frac{\partial \delta f}{\partial v_\parallel} = - \mathbf{v}_{G1}\cdot \nabla f_M - \dot{\mathcal E}_k \frac{\partial f_M}{\partial \mathcal E_k},$$9

The resulting estimate

$$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$0

shows how a large $$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$1 can be generated from pre-existing radial fine structure, thereby amplifying the EPN term (Chen et al., 25 Sep 2025).

The same work gives a plausibility-level nonlinear damping estimate,

$$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$2

with

$$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$3

at $$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$4, increasing to roughly

$$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$5

by $$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$6. A mass-ratio scan over

$$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$7

finds an approximate scaling

$$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$8

which the authors interpret as strong evidence for an electron-streaming-based phase-space mechanism. The paper is careful, however, to note that a quantitative confirmation would require direct diagnostics of the evolving velocity-space structure, and its model neglects $$\dot v_{\parallel 1}\,\frac{\partial \delta f}{\partial v_\parallel}.$$9 (Chen et al., 25 Sep 2025).

3. Trapping nonlinearity in Langmuir waves

In laser-plasma interaction, the phenomenon most closely corresponding to EPN is electron trapping nonlinearity in Langmuir waves. Here the wave acts on electrons with parallel velocity near the phase velocity

$$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$0

and trapped electrons execute bounce orbits inside the electrostatic potential troughs. The resulting flattening of the velocity distribution near $$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$1 reduces Landau damping and contributes to nonlinear frequency shift. The paper emphasizes a crucial distinction: it is not enough for electrons to lie instantaneously inside the separatrix; the distribution must have time to flatten, which takes of order a bounce period (Strozzi et al., 2012).

For a Langmuir wave with density perturbation amplitude

$$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$2

the bounce frequency and bounce period are

$$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$3

The paper introduces the bounce number

$$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$4

where $$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$5 is the detrapping rate. The practical onset criterion is

$$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$6

For a time-dependent wave, the dynamic bounce number is

$$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$7

and the simulations show that significant trapping response begins near

$$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$8

The authors explicitly stress that this is not a hard threshold; trapping nonlinearity increases continuously with $$p_\parallel = v_\parallel + \frac{q}{m}A_\parallel,$$9 (Strozzi et al., 2012).

The paper treats two classes of detrapping process. Convective loss arises because a finite wavepacket has finite longitudinal and transverse extent. For transverse side loss,

$$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$0

with

$$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$1

and the corresponding threshold amplitude is

$$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$2

Collisional detrapping is treated with a Fokker–Planck operator that includes electron-ion pitch-angle scattering, electron-electron pitch-angle scattering, and electron-electron parallel drag and diffusion in one framework. The collisional bounce number scales as

$$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$3

so collisions become less effective at detrapping as the wave amplitude increases (Strozzi et al., 2012).

This threshold theory is supported by 2D-2V Vlasov simulations with the code loki. In the driven-wave simulations, trapping nonlinearity increases continuously with $$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$4 and is significant for

$$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$5

For a single NIF-like $$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$6 speckle, the paper concludes that convective transverse loss is usually the most potent detrapping process, and under typical inner-cone conditions local reflectivities around $$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$7 are estimated to produce significant trapping effects. This makes trapped-electron EPN a quantitatively accessible effect in stimulated Raman backscatter, rather than a merely asymptotic kinetic correction (Strozzi et al., 2012).

4. Field-aligned electric fields, beam formation, and reconnection heating

A distinct but closely related usage arises in collisionless magnetic reconnection, where nonlinear electron dynamics are driven by magnetic-field-aligned electric fields. In a statistical study of 140 MMS magnetotail reconnection outflows, the acceleration potential

$$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$8

is inferred from the morphology of electron velocity distribution functions rather than from direct field-line integration. The paper reports that heating by $$\dot p_{\parallel 1}\,\frac{\partial \delta f_{p_\parallel}}{\partial p_\parallel}.$$9 in the reconnection outflow can reach up to ten times the inflow electron temperature, and that the potential scales approximately as

$O(\delta^3)$0

The authors interpret this as evidence that $O(\delta^3)$1 is set by the need to maintain quasi-neutrality and confine or trap electrons in the ion diffusion region (Richard et al., 2024).

The inferred nonlinear pathway is beam acceleration followed by wave-driven thermalization. The paper describes eVDFs that consist of a cold beam accelerated toward the X-line and a thermalized hot population. These non-equilibrium distributions can drive Buneman and electron two-stream instabilities, whose electrostatic waves scatter the beam and reshape the eVDF into a flat-top distribution. This links EPN, in the reconnection context, to field-aligned beam formation, trapping, pressure anisotropy, and phase-space flattening rather than to a single fluid coefficient (Richard et al., 2024).

The inference procedure is technically specific. The field-aligned part of the eVDF is fit with a generalized $O(\delta^3)$2 distribution,

$O(\delta^3)$3

using a Levenberg–Marquardt algorithm minimizing the reduced weighted $O(\delta^3)$4. Only energies satisfying

$O(\delta^3)$5

are included, and pitch-angle selection is restricted to within $O(\delta^3)$6 of $O(\delta^3)$7 and $O(\delta^3)$8 in the spacecraft frame. The “knee” of the plateau is defined using the threshold

$O(\delta^3)$9

which yields

$O(\delta^2)$0

This is therefore an inferred potential from nonlinear eVDF morphology, not a direct electrostatic potential measurement (Richard et al., 2024).

The paper also uses electron force balance along the field,

$O(\delta^2)$1

to interpret the spatial structure of the potential. In the case study, the observed drop from current-sheet center to edge is

$O(\delta^2)$2

decomposed approximately into

$O(\delta^2)$3

so the density term dominates. The observations favor an ambipolar, quasi-neutrality-enforcing interpretation, although the paper explicitly does not separate electrostatic and inductive contributions to $O(\delta^2)$4 directly. It further argues that as $O(\delta^2)$5 increases, $O(\delta^2)$6 becomes increasingly important to the ion-to-electron energy partition in reconnection (Richard et al., 2024).

5. Nonlocal closure formulations for parallel electron response

A different branch of the literature addresses EPN at the level of fluid closure rather than full nonlinear phase-space evolution. The papers on electron parallel closures formulate the parallel heat flux, friction force, and viscosity as kernel-weighted integrals of thermodynamic drives along the magnetic field line. The central claim is that when collisionality decreases, electrons free-stream along the field, so the response at one point depends on remote locations measured in units of the electron mean free path rather than on purely local gradients. The closed quantities are the parallel heat flux density $O(\delta^2)$7, the parallel friction force density $O(\delta^2)$8, and the parallel viscosity $O(\delta^2)$9, driven by the parallel temperature gradient, the relative electron-ion parallel flow, and the parallel rate-of-strain tensor (Ji et al., 2019).

In this formulation, the response is spatially nonlocal in a normalized parallel coordinate. The fitted kernels provide a compact surrogate for a very large moment system and are written in the universal form

$p_\parallel$0

For $p_\parallel$1, the fits were constructed to reproduce the 6400-moment solution in the converged regime together with the collisionless asymptotic behavior. The same work reports a concrete failure mode of local closure: Braginskii’s local closure for the temperature-gradient-driven heat flow has about $p_\parallel$2 error at

$p_\parallel$3

and about $p_\parallel$4 error at

$p_\parallel$5

For reduced models that require an EPN-aware representation of electron parallel transport, this is precisely the regime where local conductivity becomes quantitatively unreliable (Ji et al., 2019).

The extension to arbitrary effective charge is provided by the companion study for various ion charge numbers. It keeps the same analytic kernel form and tabulates parameters for

$p_\parallel$6

with smooth $p_\parallel$7 dependence suitable for interpolation to noninteger $p_\parallel$8. The recommended interpolation is linear in the parameters $p_\parallel$9, together with asymptotic constraints for selected kernels so that the correct collisionless limit is preserved. Benchmark errors of the fitted-kernel closures relative to the 6400-moment closures are reported as typically a few percent or less over the convergent regime, with worst cases under $v_\parallel$0 (Ji et al., 2019).

These closures are directly relevant to EPN-oriented fluid or extended-fluid modeling, but they should not be confused with a complete nonlinear kinetic theory. The arbitrary-collisionality closure paper explicitly states that it is a linearized, moment-based, nonlocal closure model. It captures nonlocal parallel electron transport, finite-collisionality corrections, and cross-coupling between thermal gradients, relative flows, and strain, but it does not provide a fully nonlinear electron kinetic closure, nonlinear phase-space trapping, nonlinear Landau damping, or self-consistent nonlinear evolution of the electron distribution function. A plausible implication is that these kernels are best viewed as closure ingredients for EPN-aware reduced models, rather than as a full theory of EPN in the strict nonlinear kinetic sense (Ji et al., 2019, Ji et al., 2019).

6. Hot-electron beam coupling and conceptual extensions

Outside plasma transport and reconnection, a related electron-mediated nonlinearity appears in epsilon-near-zero ITO. The two-beam coupling study does not use the term EPN explicitly, but it is highly relevant if EPN is understood as an electron-mediated nonlinear interaction between two beams in a parallel or crossing-beam geometry. The physical mechanism is hot-electron nonlinearity: optical absorption heats the conduction electrons, the electron temperature $v_\parallel$1 rises during the pulse, and the optical properties depend on $v_\parallel$2 through temperature-dependent Drude parameters. The temperature dynamics are described with a two-temperature model,

$v_\parallel$3

and the dielectric function is written as

$v_\parallel$4

The mechanism is therefore a temperature-driven Drude-parameter renormalization involving both $v_\parallel$5 and $v_\parallel$6, not an instantaneous bound-electronic Kerr response (Paul et al., 2020).

The central two-beam result is that pump and probe interference creates a spatially periodic heating term and hence a temperature grating,

$v_\parallel$7

Because $v_\parallel$8 depends on $v_\parallel$9, the medium develops a dielectric grating that diffracts pump light into the probe direction. The output field in the probe channel contains both the modified probe and the diffracted pump contribution: $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$00

$$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$01

This is a beam-coupled nonlinear response proportional to a field product rather than merely to pump intensity alone, and it produces polarization-sensitive transients that can mimic an instantaneous $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$02 unless the grating mechanism is taken into account (Paul et al., 2020).

The experiment uses a $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$03 ITO film on fused silica, near the ENZ wavelength of $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$04, with $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$05 pulses at $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$06 and a pump–probe angle of $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$07. The observed smooth hot-electron response decays with a characteristic time around $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$08–$$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$09, while the coherent two-beam-coupling signal appears only for co-polarized beams during temporal overlap. The authors fit effective nonlinear coefficients

$$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$10

and show that in the continuous-wave nondegenerate limit the delayed response renormalizes them as

$$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$11

$$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$12

At degeneracy, each nonlinear coefficient is enhanced by a factor of $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$13; away from degeneracy, refractive and absorptive nonlinearities mix. The paper also shows that chirp values of $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$14, $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$15, $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$16, and $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$17 modify the polarization-anisotropy transient because the coupling depends on $$\dot p_{\parallel 1}\,\partial \delta f_{p_\parallel}/\partial p_\parallel$$18 (Paul et al., 2020).

This optical literature defines neither a gyrokinetic EPN operator nor a Langmuir-wave trapping threshold. Nevertheless, it provides a concrete example of electron-mediated nonlinear coupling in which the internal state variable is electron temperature, the response is delayed rather than instantaneous, and beam interaction is mediated by an interference-induced grating. A plausible implication is that the broader EPN label can encompass not only parallel kinetic transport and reconnection physics, but also hot-electron beam coupling whenever the dominant nonlinearity is carried by electron dynamics along a privileged propagation geometry (Paul et al., 2020).