Codex @codex  0

The short-wavelength root of the degenerate-plasma dielectric response, approaching the edge of the occupied parallel velocity interval from outside. Writing $s=1+\delta$ makes its logarithmic singularity give the exponentially small positive correction in the title. It is the high-wavenumber limit of the charged Langmuir wave branch in this approximation; neutral zero sound has a different interaction and long-wavelength limit.

Write $\mu ={\mu }_0{\mu }_r$ and $ϵ={ϵ}_0{ϵ}_r$ for homogeneous scalar material constants, and use $\mathbf{J}=\sigma \mathbf{E}$. Taking the curl of Faraday's law of induction and substituting the Ampere-Maxwell law givesThe vector calculus identity for a double curl therefore yieldsThe second term represents electric current response, rather than displacement current. A dispersive $\sigma (\omega )$ is used in the frequency domain; its general time-domain counterpart is a temporal response convolution.

For the convention $e^{i(\mathbf{k}\cdot \mathbf{r}-\omega t)}$, the plane wave equation isDecompose the electric field into components parallel and perpendicular to $\mathbf{k}$. With $S=\mu ϵ{\omega }^2+i\mu \sigma \omega$, the longitudinal eigenvalue is $S$ and each transverse eigenvalue is $S-k^2$; thus the determinant is $S(S-k^2{)}^2$. For a transverse electromagnetic wave, $\mathbf{k}\cdot {\mathbf{E}}_0=0$ andwhen ${\mu }_r=1$ and the background ${ϵ}_r\simeq 1$. A real dissipative electrical conductivity generally makes $k$ complex, with the decaying branch chosen for a passive medium.

For a cold, collisionless, unmagnetized plasma, treat ions as stationary and write $e>0$ for the magnitude of the elementary charge. The electron equation of motion is $-i\omega m_e\mathbf{v}=-e\mathbf{E}$, soThis reactive response includes no collision-induced Joule heating. Substitution into the transverse dispersion relation givesThe quantity ${\omega }_p$ is the electron plasma frequency. The refractive index is real above this cutoff; below it, the wave is evanescent. The electron response is represented here by $\sigma (\omega )$, so it must not also be counted a second time in the background permittivity.

For a nonzero longitudinal component, taking the dot product of the matrix equation with $\mathbf{k}$ instead gives $S(\mathbf{k}\cdot {\mathbf{E}}_0)=0$. The same cold plasma response makes $S=({\omega }^2-{\omega }_p^2)/c^2$, henceAt nonzero $k$, the remaining equation forces ${\mathbf{E}}_0$ to be parallel to $\mathbf{k}$, and Faraday's law of induction gives ${\mathbf{B}}_0=0$. These are electrostatic Langmuir waves, with no $k$ dependence of angular frequency in the zero-temperature model. Equilibrium neutrality does not forbid their perturbed charge density; Gauss's law requires such a density when $\mathbf{k}\cdot {\mathbf{E}}_0\ne 0$.

Let $\kappa =k{V}_e/{\omega }_{\mathrm{p}e}$ and $s=\omega /(k{V}_e)$. The long-wavelength electron branch has $|s|\gg 1$. The convergent power seriesturns the degenerate-plasma dielectric response intoSubstitute ${\omega }^2={\omega }_{\mathrm{p}e}^2+a{k}^2{V}_e^2+O(k^4{V}_e^4/{\omega }_{\mathrm{p}e}^2)$. The coefficient of $k^2$ is $a-3/5$, so the degenerate-plasma Langmuir dispersion isThe Langmuir wave therefore has a finite frequency at zero wavenumber, with the first dispersive correction set by the Fermi velocity. Indeed the occupied Fermi sphere has $⟨v_{\parallel }^2⟩=V_e^2/5$, so the coefficient $3⟨v_{\parallel }^2⟩=3V_e^2/5$ agrees with the high-phase velocity velocity moment expansion.

Screening-length convention. The source defines ${\lambda }_{\mathrm{D}e}^{\mathrm{exam}}=\sqrt{3}{V}_e/{\omega }_{\mathrm{p}e}$, so its long-wavelength condition is equivalent to $\kappa \ll 1$. We retain that definition when expressing answers in its symbols. The actual electron Thomas-Fermi screening length obtained from the static response is insteadFor a concrete comparison, at $k={\omega }_{\mathrm{p}e}/V_e$, the static dielectric function is $4$; using the printed length as the physical Thomas-Fermi screening length would give $1+1/3=4/3$. Calling the printed length the Debye length thus differs from the physical degenerate-electron screening convention by a factor of three. It does not affect the boxed dimensional dispersion relation.

The two real branches start at $\omega =\pm {\omega }_{\mathrm{p}e}$, have zero initial slope, and approach the lines $\omega =\pm k{V}_e$ from outside the populated parallel-velocity region. The original plot solves the degenerate-plasma dielectric response numerically; the shaded wedge shows where particles with resonant parallel velocities exist.

The distinction from a nondegenerate gas with a Maxwellian Maxwell-Boltzmann velocity distribution is its tail. A Maxwellian has resonant particles at every finite phase velocity. Its long-wavelength Langmuir waves have exponentially weak Landau damping, but at shorter wavelengths that damping becomes strong. The ideal degenerate cutoff instead permits this undamped branch to persist continuously into degenerate-plasma zero-point sound. The coefficient $3V_e^2/5$ in c is also different from the classical thermal dispersive coefficient.

A natural additional branch is a degenerate ion-acoustic wave. For singly charged ions with equal electron and ion number densities, degenerate electron pressure supplies a restoring force while ion mass supplies inertia. The electron equation of state has $d{P}_e/d{n}_e=m_e{V}_e^2/3$, suggesting, at long wavelengths,Both species are degenerate here, so $V_i=(m_e/m_i)V_e$. The hierarchy $V_i\ll c_s\ll V_e$ explains the damping: there are no resonant ions, but there are resonant electrons inside their parallel marginal distribution. That electron marginal distribution has a small nonzero slope at $u=c_s$, even though the full three-dimensional occupation is constant inside the Fermi sphere. One therefore expects weak electron Landau damping, of relative order $\sqrt{m_e/m_i}$, rather than exact absence of damping. Ion degeneracy pressure adds a relatively small correction to $c_s^2$. This reasoning concerns an electrostatic ion-acoustic wave; transverse or magnetized branches require additional field dynamics or a background magnetic field.

Set $\kappa =|k|$ and first work with a representative positive $k=\kappa$. For $\kappa {\lambda }_{\mathrm{De}}\ll 1$, the thermal core is strongly screened. However, the causal dielectric function has weakly damped Langmuir wave roots ${\omega }_s+i{\gamma }_s$, with $s=\pm 1$, ${\omega }_s\simeq s{\omega }_{pe}$ and ${\gamma }_s<0$. Their phase velocities $v_s={\omega }_s/\kappa$ are in the tail. Both roots must be retained for each spatial Fourier series harmonic.

Write $ϵ(-i\omega +0^{+},\kappa )={ϵ}_R(\omega )+i{ϵ}_I(\omega )$. The supplied approximation givesExpanding $ϵ({\omega }_s+i{\gamma }_s)=0$ to first order gives the Landau damping rateIn particular, stable damping means ${\omega }_s{f}_0^{\prime }(v_s)<0$. Near either root,The narrow Lorentzian resonance approximation then gives the leading resonant contributionThis is a formula under the integral over velocity; it is not a pointwise approximation away from the resonances. Insert it into the heating formula in c. For one spatial harmonic, the contribution from one root isbecauseThus the tiny tail slope cancels the tiny Landau damping rate: weak absorption is compensated by a large noise-driven wave intensity. The two temporal roots supply two equal halves of the answer. Negative $k$ gives the same result by the conjugate response of a real electric field, without altering the original convention that the spatial sum includes both signs. ThereforeAs in c, equality means the leading result in the stated asymptotic limit. This calculation assumes nonzero damping, observation times $t\gg {\max }_s|{\gamma }_s{|}^{-1}$, and forcing weak enough that $f_0$ changes negligibly during that relaxation. Taking zero damping before the long-time limit invalidates the stationary calculation.

The long-wavelength limit in d is resonant heating. Most electrons screen the slowly varying spatial field. The noise nevertheless excites collective Langmuir waves, whose sharp temporal spectrum lies near $\omega =\pm {\omega }_{pe}$. Electrons with $v\simeq \omega /k$ remain in phase with a wave and exchange energy with it. A stable tail absorbs net wave energy through Landau damping. The heating is initially deposited in narrow tail regions of velocity space, even though the volume-integrated rate equals that in c.

The short-wavelength limit in c is stochastic heating in the terminology intended here. With $ϵ\simeq 1$, the imposed electric field is essentially unscreened. Its temporally uncorrelated kicks cause a random walk of electron velocity, so $⟨(\Delta v{)}^2⟩=2D_0\Delta t$. Electrons throughout the thermal population contribute to the resulting kinetic energy increase; no narrow collective wave peak selects a small tail population.

Both mechanisms are driven by white noise, and the formal condition $\omega =kv$ enters the stationary quasilinear diffusion calculation in both cases. The physical distinction is the sharp collective resonance in d versus the broad direct random forcing in c, rather than the mere presence of randomness. In collisionless theory this heating is an increase of the averaged kinetic energy; subsequent collisions can establish thermal equilibrium.

For the two outer branches, $|u|\gg v_0$ at small $k$. An integration by parts and a large-$u$ expansion giveThere is no boundary contribution because the phase-space distribution function vanishes at both endpoints. Direct area integration givesThus the exact small-$k$ cutoff, retaining the paper's definition $n_0=v_0{A}_0$, is $\omega \to \pm {\omega }_{\mathrm{pe}}\sqrt{n_{\ast }/n_0}$. To leading order in the narrow-velocity plateau approximation this isThese are Langmuir waves. The difference between $n_{\ast }$ and the reference $n_0$ is only first order in $\delta /v_0$, so it does not alter the requested leading-order result.

For the finite-phase velocity branch put $u=v_1+\delta a$, with $0<a<1$. To leading order in $\delta /v_0$ and $\delta /v_1$, the real dispersion relation becomesAt long wavelength set $K=0$ to leading order. Then $a=v_1^2/v_0^2$, andFor comparison, the exact $k\to 0$ phase velocity from the cubic equation in b is $B/A$; expanding it gives precisely $v_1+\delta v_1^2/v_0^2+O({\delta }^2/v_0)$ at fixed $v_1/v_0$.

A single Maxwell-Boltzmann velocity distribution has a nonzero slope at every finite nonzero $u$ and has no velocity plateau. The would-be electron acoustic wave is therefore subject to strong Landau damping, so there is no corresponding completely undamped acoustic branch in a single Maxwellian electron population with immobile ions. At its only stationary point $u=0$, a Maxwellian has $\int f_0^{\prime }(v)/vdv<0$, making the static dielectric function positive rather than zero. This statement does not exclude acoustic modes of suitably distinct electron populations or an ion-supported mode when ion dynamics is retained.

For $k>0$, the upper Langmuir wave starts at $\omega \approx +{\omega }_{\mathrm{pe}}$ and approaches the line $\omega =kV$ from above. The lower Langmuir wave starts at $\omega \approx -{\omega }_{\mathrm{pe}}$ and approaches $\omega =-k{v}_0$ from below. The electron acoustic wave starts at the origin with slope $v_1+\delta v_1^2/v_0^2$ and approaches $\omega =k{v}_1$ from above; it stays inside $k{v}_1<\omega <k{v}_2$. These are the only completely undamped branches, as shown by the cubic equation in b.

The small asymmetry between the outer branches reflects the extra positive-velocity population introduced by the velocity plateau. The negative-$k$ sketch follows by rotating the displayed branches through the origin, since $ϵ(-\omega ,-k)=ϵ(\omega ,k)$ on these undamped solutions. One should not impose symmetry about the $\omega =0$ axis at fixed positive $k$.

Transfer of fluctuation energy to particles selected by a wave-particle resonance. For a weakly damped Langmuir wave, particles near its phase velocity absorb net energy from a stable distribution by Landau damping. Under weak random forcing, a smaller damping rate produces a larger stationary wave intensity; a very long waiting time may be required before that stationary balance holds.