Reading the 1967 plasma wave echo paper

Page 219: the argument in one paragraph

The right column of the first page gives the whole mechanism before any equation. The first pulse leaves a term $f_1(v)\exp(-ik_1x+ik_1vt)$; the second, at time $\tau$, leaves a second-order term whose velocity coefficient vanishes at $t=\tau k_2/(k_2-k_1)$. Note the sign convention: the paper writes the first pulse as $e^{-ik_1x}$ and the second as $e^{ik_2x}$, so its echo has wave number $k_3=k_2-k_1$ with both $k$ positive. Chapter 6 used the same convention.

The paper then fixes its model: one dimension, a uniform ion background, an initially uniform electron distribution $f_0(v)$, and two impulsive external potentials. Footnote 4 says $\Phi_{k_1}$ and $\Phi_{k_2}$ carry units of potential because the delta functions take $\omega_pt$ as their argument. In plasma units, with $\omega_p=1$, the impulse is the same as the potential kick used in the echo Lab.

Equations (1) to (4): from the transformed equations to the echo integral

Equations (1) and (2) are the Vlasov and Poisson equations after a Fourier transform in space and a Laplace transform in time; the external pulses appear as sources in Eq. (2). Expanding in powers of the two potentials, the first order is ordinary Landau damping after each pulse. Equation (3) is the second-order potential at $k_3$, a double contour integral with three dielectric functions in the denominator.

Equation (4) follows by picking up two poles, under stated assumptions: $\tau$ long compared with a damping period, $|\gamma(k_1)\tau|$ and the other damping exponents large, and $k_1/k_3\approx1$. The integral stops phase mixing when $k_3(t-\tau)\approx k_1\tau$, the echo time again. The paper adds that setting the dielectric functions to one gives the free-streaming result, “an easy limit in which to check Eq. (4)”. The project benchmark does that check: with the self-consistent field switched off, the simulated echo-mode density matches an independent quadrature of the free-streaming solution.

Equation (5) and Fig. 1: an echo that rises and falls at different rates

Equation (5) evaluates Eq. (4) for a Maxwellian. Before the echo the envelope grows as $\exp[\gamma(k_1)(k_3/k_1)(\tau^{\prime}-t)]$; after it, it decays as $\exp[\gamma(k_3)(t-\tau^{\prime})]$, where $\gamma(k)$ is the negative Landau damping rate. Fig. 1 sketches three panels: the response to the first pulse, the response to the second, and the echo between them, drawn for $k_3\approx k_1\approx k_2/2$.

Read the figure as a schematic. Its caption says “approximate variation”, and the text says the curves are drawn for $k_1\simeq k_3$. Its time axis has no scale and the vertical axes have no units, so there is nothing to digitise. What can be checked are the two envelope rates.

Collisions and the suggestion at the end

The last column says the calculation fails if collisions destroy the phase information before the echo. Small-angle Coulomb collisions act through a Fokker–Planck operator, which multiplies the collision rate by about $(kv\tau)^2\approx(\omega_p\tau)^2$ for a perturbation like $e^{ikv\tau}$. The authors suggest using this sensitivity to measure the collision rate. Chapter 7 turns the remark into an estimate and compares it with the simulation; the estimate overstates the suppression, and the cause of that gap remains an open check.

A checklist for the next primary source

Write down the version and date, the sign convention for every wave, and the units of every applied quantity, including footnotes such as Footnote 4. List the stated assumptions of each derived equation before comparing it with a simulation, and note whether your test case satisfies them. For a schematic figure, compare rates or ratios, not drawn amplitudes. Finally, give each claim one of four labels: published theory, published experiment, our numerical result, or still to be reproduced.