---
title: Field-Particle Correlation Technique
url: https://www.emergentmind.com/topics/field-particle-correlation-technique
type: topic
---

# Field-Particle Correlation Technique

Field–particle correlation is a plasma diagnostic that uses single-point measurements of electromagnetic fields and particle velocity distribution functions to isolate the secular transfer of energy between fields and particles and to reveal its structure in velocity space. In the plasma-kinetic literature, the technique is formulated by correlating electric-field fluctuations with velocity-space gradients of the distribution function, over a correlation interval long enough to suppress oscillatory “sloshing” and retain the net energy exchange. Its defining output is not only a local rate of energization, but also a velocity-space signature that can be compared with resonant conditions such as Landau, cyclotron, and related wave–particle interactions [1607.01738, 1705.06385, 2006.02563].

## 1. Vlasov basis and the definition of the correlation

The technique is derived from the collisionless Vlasov equation. For species \(s\),
\[
\frac{\partial f_s}{\partial t} + \mathbf v\cdot \nabla f_s + \frac{q_s}{m_s}\left[\mathbf E + \frac{\mathbf v\times \mathbf B}{c}\right]\cdot \frac{\partial f_s}{\partial \mathbf v}=0.
\]
Multiplication by the single-particle kinetic energy \(m_s v^2/2\) defines the phase-space energy density,
\[
\Theta_s(\mathbf r,\mathbf v,t)\equiv \frac{m_s v^2}{2}f_s(\mathbf r,\mathbf v,t),
\]
or equivalently \(w_s=\tfrac12 m_s v^2 f_s\) in the 1D–1V formulation. Its evolution contains a ballistic term, a field–particle term, and a Lorentz-force term. Under periodic or open boundary conditions the ballistic contribution vanishes after the appropriate spatial integration, and the \((\mathbf v\times \mathbf B)\) term integrates to zero over velocity, so the electric-field term is isolated as the part responsible for net energization [1607.01738, 2006.02563].

This leads to the central correlation
\[
C_E(\mathbf r,\mathbf v,t,\tau)\equiv
\left\langle
-\,q_s\,\frac{v^2}{2}\,
\left[\frac{\partial \delta f_s(\mathbf r,\mathbf v,\cdot)}{\partial \mathbf v}\right]\cdot
\mathbf E(\mathbf r,\cdot)
\right\rangle_{t-\tau/2}^{t+\tau/2},
\]
where \(\delta f_s\equiv f_s-F_{0,s}\), \(F_{0,s}\) is the time-average over the interval \(T\), and \(\langle\cdot\rangle\) denotes a uniform average over the time window \(\tau\). In the original 1D–1V electrostatic formulation this appears as
\[
C_1(v,t_i;\tau)\equiv \frac{1}{N}\sum_{j=i}^{i+N-1}
\left[-\,q_s\,\frac{v^2}{2}\,\frac{\partial \delta f_{sj}(v)}{\partial v}\,E_j\right].
\]
Provided \(\tau\) exceeds the characteristic wave period of the dominant fluctuation, the fast oscillatory exchange cancels out in the sum, leaving a nonzero mean only where secular damping or growth occurs [1607.01738].

A central conceptual point is that the method retains the full \((v,t)\) dependence of the energy-transfer term at a fixed point. The velocity integral of the correlation recovers the local field–particle work, but the unintegrated structure shows which particles gain or lose energy. This is the basis for mechanism identification from single-point data [1705.06385].

## 2. Operational forms of the technique

In magnetized plasmas the electric field is split into components parallel and perpendicular to the local magnetic field. In field-aligned coordinates, the hybrid Vlasov–Maxwell implementation defines
\[
C_{E_\parallel}(\mathbf r_0,\mathbf v,t,\tau)\equiv
\left\langle
-\,q_p\,\frac{v_\parallel^2}{2}\,
\frac{\partial \delta f_p(\mathbf r_0,\mathbf v,t')}{\partial v_\parallel}\,
E_\parallel(\mathbf r_0,t')
\right\rangle_{t'=t-\tau/2}^{t+\tau/2},
\]
and
\[
C_{E_\perp}(\mathbf r_0,\mathbf v,t,\tau)\equiv
\sum_{i=1,2}
\left\langle
-\,q_p\,\frac{v_{\perp i}^2}{2}\,
\frac{\partial \delta f_p}{\partial v_{\perp i}}\,
E_{\perp i}
\right\rangle_{t\pm\tau/2}.
\]
The discrete correlation is
\[
C(A,B;t_i,\tau)=\frac{1}{N}\sum_{j=i-N/2}^{i+N/2}A_j B_j,
\qquad \tau=N\Delta t.
\]
Velocity integration gives the local spatial energy-density transfer rates,
\[
\partial_t w_{E_\parallel}(\mathbf r_0,t,\tau)=\int d^3v\, C_{E_\parallel},
\qquad
\partial_t w_{E_\perp}(\mathbf r_0,t,\tau)=\int d^3v\, C_{E_\perp},
\]
with \(\partial_t w_E=\partial_t w_{E_\parallel}+\partial_t w_{E_\perp}\). In the Alfvén–Ion Cyclotron study, all \(w\) are normalized by the time-averaged local proton energy density
\[
w_0(\mathbf r_0)=\left\langle \int d^3v\,\frac{m_p v^2}{2}\,f_p(\mathbf r_0,\mathbf v,t)\right\rangle_T,
\]
and the accumulated transfer is
\[
\Delta w_{E_l}(\mathbf r_0,t)=\int^t dt'\,\partial_{t'} w_{E_l}(\mathbf r_0,t') .
\]
These constructions make the technique quantitative rather than merely qualitative [2006.02563].

A practically important variant avoids the velocity derivative. For 1D–1V electrostatic data one may define
\[
C_2(v,t_i;\tau)\equiv \frac{1}{N}\sum_{j=i}^{i+N-1}
\left[q_s\,v\,\delta f_{sj}(v)\,E_j\right],
\]
and in gyrokinetic form
\[
C'_{E\parallel}(v,t;\tau)=
C\!\left(q_s v_\parallel f_s(\mathbf r_0,\mathbf v,t),\,E_\parallel(\mathbf r_0,t)\right).
\]
By integration by parts in velocity, \(\int dv\,C_2=\int dv\,C_1=j_sE\), and the same holds for the gyrokinetic proxy. This derivative-free form is explicitly motivated by the difficulty of computing \(\partial_v\delta f\) from noisy, discretized spacecraft measurements [1607.01738, 1705.06385].

## 3. Velocity-space signatures and mechanism identification

The distinctive feature of the technique is that different energization mechanisms imprint different velocity-space signatures. In resonant interactions such as Landau damping, the correlation is localized near the resonant velocity \(v_{\rm res}=\omega/k\) or \(v_{\rm res}=\omega/k_\parallel\). In the electrostatic tutorial, if the wave is damped, \(C_1(v)\) is positive for \(v\gtrsim v_{\rm res}\) and negative for \(v\lesssim v_{\rm res}\), corresponding to flattening of the VDF around \(v_{\rm res}\) [1607.01738].

In gyrokinetic turbulence, where the dominant resonance is \(n=0\), the reduced parallel correlation \(C_{E\parallel}(v_\parallel)\) shows a bipolar signature in \(v_\parallel\) about \(v_{\rm res}\). For fully developed 3D turbulence at \(\beta_p=0.3,1.0,3.0\), the net \(C_{E\parallel}\) is strongly localized in \(|v_\parallel|\) around the range of \(v_{\rm res}\) corresponding to peak proton damping, and a resonant fraction
\[
R = \frac{\int_{v_1}^{v_2}|C_{E\parallel}(v_\parallel)|\,dv_\parallel}
{\int_{\rm all\,v}|C_{E\parallel}(v_\parallel)|\,dv_\parallel}
\]
is found to satisfy \(R\sim 0.7\)–\(0.9\), well above the \(\sim0.1\)–\(0.5\) expected for a uniform-in-velocity transfer [1705.06385].

The Alfvén–Ion Cyclotron turbulence study extends the method beyond purely Landau channels. There, \(C_{E_\parallel}(v_\parallel,v_\perp)\) shows the canonical bipolar signature of Landau damping around \(v_\parallel\simeq \omega/k_\parallel\), with negative correlation for \(|v_\parallel|<v_{\rm res}\) and positive for \(|v_\parallel|>v_{\rm res}\), extended out to \(v_\parallel\sim 2\)–\(3\,v_{tp}\), and with weak \(v_\perp\) dependence for \(v_\perp\lesssim 1.5\,v_{tp}\). By contrast, \(C_{E_\perp}(v_\parallel,v_\perp)\) peaks in a broad band \(1\lesssim v_\perp/v_{tp}\lesssim 3\), \(|v_\parallel|\lesssim 1.3\,v_{tp}\), forming a “quasilinear plateau” consistent with contours
\[
(v_\parallel\pm\omega/k_\parallel)^2+v_\perp^2={\rm const}.
\]
The two channels are therefore separable even when they act simultaneously [2006.02563].

The method also distinguishes unstable mechanisms. In counter-streaming electrostatic instability, each electron beam’s \(C_E(v)\) shows a sign-flip at its drift speed \(v_d\), with no sharp structure at \(\omega/k\), and the ion \(C_E(v)\) is even in \(v\) and positive overall, implying net ion heating. In bump-on-tail instability, the bump population shows \(C_E<0\) around the acoustic resonant velocity, while the core population shows \(C_E>0\) peaked at the same resonance [1701.03687].

In a perpendicular collisionless shock, the ion correlation \(C_{E_y}\) exhibits nested blue/red crescents around the reflected population, identifying shock-drift acceleration, while in a transverse-drift frame the electron \(C_{E_y''}(v_x'',v_y'')\) yields a two-lobed positive signature of adiabatic heating. These signatures were reproduced with simplified “step-shock” and “linear-ramp” models, giving a direct interpretive bridge between Eulerian correlations and Lagrangian particle dynamics [2011.13829].

## 4. Numerical, spectral, and analytic realizations

The first explicit development of the technique was in the electrostatic 1D–1V Vlasov–Poisson system, where it was introduced as a practical, single-point means to diagnose how collisionless plasma fluctuations transfer energy to particles and to discuss how the method could be implemented on spacecraft data in the solar wind [1607.01738].

It was then extended to strongly driven electromagnetic gyrokinetic turbulence. In that setting, gyrokinetics provides the rigorous low-frequency, anisotropic limit of Vlasov–Maxwell, retaining Landau damping while filtering out high-frequency cyclotron resonances. Single-point time series of the complementary gyrocenter distribution \(g_s\) and \(E_\parallel\) are recorded at fixed points, and the gyrotropic correlation \(C_{E\parallel}(v_\perp,v_\parallel,t;\tau)\) or its reduced form is constructed directly from the turbulent simulation output [1705.06385].

A spectral generalization decomposes the transfer rate mode by mode in a periodic domain. Writing
\[
E(\mathbf x,t)=\sum_{\mathbf k}E_{\mathbf k}(t)e^{i\mathbf k\cdot\mathbf x},
\qquad
f_s(\mathbf x,\mathbf v,t)=\sum_{\mathbf k}f_{s,\mathbf k}(\mathbf v,t)e^{i\mathbf k\cdot\mathbf x},
\]
one obtains
\[
C_{E_\parallel}(\mathbf k,\mathbf v,t;\tau)
=
C\!\left[
-\,q_s\left(\frac{v_\parallel^2}{2}\right)\frac{\partial g^*_{s,\mathbf k}(\mathbf v,t)}{\partial v_\parallel},
\,E_{\parallel,\mathbf k}(t)
\right],
\]
and
\[
\left(\frac{dW_s}{dt}\right)_{\mathbf k}=\int d^3v\,C_{E_\parallel}(\mathbf k,\mathbf v,t;\tau).
\]
This supplies scale-by-scale accounting of energy flow in kinetic turbulence, but the same work notes that it requires full 3D spatial information and is therefore impractical for single-spacecraft data [1904.09334].

The hybrid Vlasov–Maxwell implementation resolves both Landau and cyclotron channels. In the reported simulation, the ion Vlasov equation is integrated on a \(32^3\) grid in physical space and \(51^3\) in velocity space, with isothermal fluid electrons, box size \(L=10\pi d_p\), \(B_0\parallel \hat z\), \(\beta_p=1\), \(T_p=T_e\), and a velocity grid \(\pm 5\,v_{tp}\). Sixty-four fixed points are selected, and at each point the full 3-dimensional ion distribution and the simulation-frame fields are recorded. The analysis transforms the fields and \(\partial_v\delta f\) to the local mean-flow frame and rotates into an instantaneous field-aligned basis before computing \(C_{E_\parallel}\) and \(C_{E_\perp}\) [2006.02563].

The shock application uses a 1D–2V continuum Vlasov–Maxwell simulation in Gkeyll. There, because the self-consistent shock is quasi-stationary, the analysis uses \(\tau\to0\) and directly evaluates componentwise instantaneous correlations \(C_{E_x}\) and \(C_{E_y}\) in the shock-rest frame, after the Galilean transformation \(\mathbf v' \to \mathbf v=\mathbf v'-\mathbf U_{\rm shock}\) and \(\mathbf E'\to \mathbf E=\mathbf E'+\mathbf U_{\rm shock}\times \mathbf B\) [2011.13829].

A recent analytic realization is JET-PLUME, an extension of the PLUME hot-plasma dispersion solver. In Fourier space and for a single \((\mathbf k,\omega)\) mode in the weak-damping limit,
\[
C_s(\mathbf k,\omega;\mathbf v)=\Re\{\hat E(\mathbf k,\omega)\,\hat g_s^*(\mathbf k,\omega;\mathbf v)\},
\]
or equivalently
\[
C_s(\mathbf k,\omega;\mathbf v)=
-\,q_s\frac{v_j^2}{2}\,
\Re\left\{\left(\partial_{v_j}\hat g_s\right)\hat E_j^*\right\}.
\]
Because \(\hat g_s\) contains the resonant denominator \(\omega-k_\parallel v_\parallel-n\Omega_s\), the analytic correlation isolates individual resonances directly and can separate degenerate entropy-mode components [2606.20865].

## 5. Quantitative findings and observational reach

The Alfvén–Ion Cyclotron turbulence study supplies one of the clearest quantitative demonstrations of the technique. Averaging over 64 points and over one nonlinear turnover time yields
\[
\left\langle\,\left|\partial_t w_{E_\parallel}\right|\,/\,\left|\partial_t w_E\right|\,\right\rangle
=
0.62\pm 0.24,
\]
so \(E_\parallel\) mediates \(\sim 60\%\) of the total transfer and \(E_\perp\) \(\sim 40\%\). At every point, the velocity-space pattern is reproducible: bipolar in \(v_\parallel\) for the parallel channel and plateau-like in \(v_\perp\) for the perpendicular channel. Over the chosen correlation window \(\tau\Omega_p=22.5\), approximately \(0.7\) of the largest-scale Alfvén period, the oscillatory component is largely removed and the net secular patterns persist with small scatter across points [2006.02563].

The Fourier-space treatment shows that, in a \(\beta_i=1\), solar-wind-like, low-frequency Alfvénic turbulence simulation, all diagnosed Fourier modes covering the dissipation range display resonant energy transfer localized at the Landau resonances for each mode. In the reported scale survey, the width in \(v_\parallel\) remains narrow, \(\lesssim 0.2\,v_{te}\), the accumulated \(\Delta E_s^{(\mathrm{fp})}(k)\) peaks at the lowest \(k_\perp\) mode and drops by \(\gtrsim 10^3\) across the dissipation range, electrons and ions receive comparable transfer at \(k_\perp\rho_i\approx 1\), and electron transfer remains significant out to \(k_\perp\rho_i\sim 8\)–\(10\) whereas ions are effectively damped only around \(k_\perp\rho_i\sim 1\)–\(2\) [1904.09334].

The technique has also been tested under severe cadence limitations. In a high-resolution gyrokinetic simulation, 38 distinct Landau signatures were tracked through systematic downsampling; 26, approximately \(68\%\), remained detectable at super-Nyquist frequencies, in some cases with \(f/f_{\rm Ny}\gg 10\). A rule of thumb was extracted from the undersampling study: to recover the true time-average of the product within \(10\%\), one needs at least \(\mathcal N\approx 12\) samples spanning the correlation interval, or \(\tau/\Delta t\gtrsim 12\). Signatures appearing with fewer than twelve samples in \(\tau\) may not be statistically significant [2204.00104].

The observational motivation is explicit throughout the literature. The original proposal discussed application to spacecraft missions such as the Magnetospheric Multiscale and Solar Probe Plus missions, and the later turbulence studies state that the method is readily applicable to in situ spacecraft data because it requires only local \(f_s\) and \(E\) [1607.01738, 2006.02563]. Electron Landau damping has been found in both simulations and observations from Earth’s magnetosheath using this technique, and the undersampling analysis concludes that Parker Solar Probe should be able to identify Landau damping signatures by correlating over intervals of a few seconds [2204.00104].

## 6. Assumptions, interpretive issues, and recent extensions

The method depends critically on the choice of correlation interval \(\tau\). The interval must be long enough to cancel oscillatory exchange and short enough not to smear secular evolution. In the gyrokinetic and electrostatic treatments, \(\tau\) is typically chosen to exceed one or two wave periods of the target mode or broadband fluctuation. The Fourier-space study notes that overlapping mode spectra can complicate isolation if mode lifetimes are short [1607.01738, 1904.09334].

An important interpretive caution is that the technique isolates secular transfer, but collisionless interactions are not simply equivalent to monotonic heating. In the Alfvén–Ion Cyclotron simulation, the instantaneous sign of local transfer can reverse because collisionless interactions are reversible, even though the time-averaged velocity-space signatures persist. A plausible implication is that the correlation is best understood as a diagnostic of phase-space energization channels, with collisional thermalization remaining a separate step [2006.02563, 1904.09334].

Noise and discretization in velocity space are a recurrent practical limitation. Both the original and gyrokinetic papers therefore introduce alternative proxy correlations that avoid direct computation of \(\partial f/\partial v\), while retaining the same integrated \(j\cdot E\). This is particularly relevant for spacecraft VDFs with finite cadence and finite velocity resolution [1607.01738, 1705.06385].

Model assumptions also delimit which mechanisms can be accessed. Gyrokinetics filters out high-frequency cyclotron resonances and excludes stochastic ion heating in large-amplitude turbulence, so perpendicular energization requires extensions beyond gyrokinetics. The full-Vlasov and hybrid Vlasov–Maxwell studies were introduced precisely to diagnose cyclotron, shock, and related perpendicular channels [1705.06385, 2006.02563, 2011.13829].

Recent analytic work extends the technique from diagnostics to controlled resonance isolation. In JET-PLUME, the susceptibility tensor is used to decompose parallel energy transfer into a direct Landau term and an off-diagonal term,
\[
\int C_s\,d^3v
=
\frac{\omega_r}{8\pi}\Im[\chi_{zz,s}|E_z|^2]
+
\frac{\omega_r}{8\pi}\Im[\chi_{zy,s}E_yE_z^*].
\]
For kinetic Alfvén waves at \(k_\perp\rho_i\gtrsim 1\), the off-diagonal channel can drive a parallel ion current at large \(v_\perp\) with phase relative to \(E_\parallel\) that returns field energy, partially counteracting the direct Landau damping from \(\chi_{\parallel\parallel}\). The resulting velocity-space signature “twists” sign as \(v_\perp\) increases. Because the analytic solution is built from a sum over cyclotron harmonics \(n\), Landau (\(n=0\)) and cyclotron (\(n=\pm1\)) resonances can be isolated by retaining the corresponding term in the Bessel sums, and degenerate entropy modes can be combined or separated by keeping the mixed cross terms that survive the time average [2606.20865].

Taken together, these developments define the field–particle correlation technique as a family of closely related diagnostics that transform local measurements of \(f_s(\mathbf v,t)\) and fields into a phase-space map of collisionless energy transfer. Its distinctive contribution is that it does not stop at the scalar work rate \(j\cdot E\): it identifies which particles, at which velocities, are involved in the exchange, and thereby ties measured energization directly to the underlying physical resonance.

Source: https://www.emergentmind.com/topics/field-particle-correlation-technique