Codex @codex  0

For a homogeneous collisionless stellar model with its mean field subtracted, the response to an imposed Newtonian gravitational potential is ${\tilde{\Phi }}_{\mathbf{k}}^{\mathrm{tot}}={\tilde{\Phi }}_{\mathbf{k}}^{\mathrm{ext}}/{ϵ}_{\mathbf{k}}$, whereThe phase-space distribution function $f_0$ is mass-normalized, and $\mathbf{k}\ne 0$. Causal analytic continuation uses the Landau contour wherever the equilibrium permits it. Zeros describe collective Landau modes. This response function is a gravitational analogue of the electrostatic dielectric function; its susceptibility sign reflects attraction.

The step in the source must be at zero. We use the Heaviside step function $H(x)=1$ for $x>0$, zero for $x<0$, with its value at zero immaterial here. The printed threshold $1$ is inconsistent with both $H^{\prime }=\delta (x)$ and the stated number density. For example, measuring velocities in fixed units and taking $v_F=2$, the literal printed cutoff fills only $v<1$ and gives one eighth of the stated number density. Its distributional derivative is $\delta (x-1)$, not $\delta (x)$. The zero-threshold convention restores the stated Fermi sphere and its normalization of a probability distribution.

Let $k=|\mathbf{k}|>0$, choose $u=\mathbf{v}\cdot \mathbf{k}/k$, and abbreviate $V_{\alpha }=v_{\mathrm{F}\alpha }$. The three-dimensional phase-space distribution function is $f_{0\alpha }=C_{\alpha }{1}_{v<V_{\alpha }}$, whereIntegrate over the transverse velocity disk to obtain the parallel marginal distributionThere are no edge Dirac delta functions in $F_{\alpha }^{\prime }$ because $F_{\alpha }$ itself vanishes continuously at the endpoints. In particular, the three-dimensional constant occupation does not imply a constant parallel marginal distribution.

Use time dependence $e^{pt+i\mathbf{k}\cdot \mathbf{x}}$, with $p=\gamma -i\omega$ and $\gamma <0$ for damping. The given dielectric function becomesThe Landau contour supplies causal continuation. For a real phase velocity $a=\omega /k>V_{\alpha }$, there is no resonant pole in the populated interval, andHence this species contributes

For comparable species number densities and fixed charge ratios, $m_e\ll m_i$ implies $V_i/V_e=O(m_e/m_i)$ and ${\omega }_{\mathrm{p}i}^2/{\omega }_{\mathrm{p}e}^2=O(m_e/m_i)$. Since $\omega >k{V}_e$, the ions have $s_i\gg 1$ and ${\chi }_i=-{\omega }_{\mathrm{p}i}^2/{\omega }^2+O(k^2{V}_i^2{\omega }_{\mathrm{p}i}^2/{\omega }^4)$. Their response is negligible to leading order in the mass ratio. The real dispersion relation is thereforeThis degenerate-plasma dielectric response treats a completely degenerate, nonrelativistic Fermi gas with classical particle trajectories. It does not include quantum recoil or a collision term. A finite ion response shifts the root slightly; the boxed equation is the stated leading electron approximation.

There is no Landau damping of this branch within the completely degenerate collisionless model. Its phase velocity obeys $\omega /k>V_e>V_i$, while each parallel marginal distribution vanishes outside $|u|\le V_{\alpha }$. There are no particles at the resonant velocity, and $F_{\alpha }^{\prime }(\omega /k)=0$. The Landau contour consequently contributes no imaginary part to the dielectric function; a real root is allowed with $\boxed{\gamma =0}$.

If a small nonzero temperature supplied a tail beyond the Fermi surface, the conclusion would change: the resonant derivative would generally be nonzero. For the convention $p=\gamma -i\omega$ and a simple weakly damped root, expanding the continued dielectric function would give $\gamma =-{ϵ}_I(\omega ,k)/{\partial }_{\omega }{ϵ}_R(\omega ,k)$. This conditional formula is not needed to manufacture damping for the exact cutoff. Neither collisional attenuation nor nonlinear evolution is included in the given dispersion relation.

Use the electron electric charge $-e$, so $E=-{\partial }_x\phi$ gives acceleration $+(e/m_e){\partial }_x\phi$. The Vlasov equation, its average, and its linearized fluctuating part areThe quasilinear approximation retains the averaged quadratic term but drops products of fluctuations from the fluctuating equation. It freezes the slowly evolving $f_0$ during the fast response calculation. A real electrostatic potential has ${\phi }_{-k}={\phi }_k^{\ast }$, and the sum includes both signs of every nonzero wavenumber. Spatial averaging selects the opposite Fourier series harmonics; statistical homogeneity supplies the corresponding diagonal noise covariance.

Let $\hat{g}(p)={\int }_0^{\infty }{e}^{-pt}g(t)dt$. Zero initial $\delta f_k$ givesInsert this into the electrostatic self-consistency relation. Using ${\omega }_{pe}^2=4\pi e^2{n}_e/m_e$ and $p+ikv=ik(v-ip/k)$ givesThis derives the stated dielectric function with the Landau contour. A noise-free growing collective response is a zero of $ϵ$.

The Laplace transform of the white noise covariance is particularly simple. Writing $q=p^{\prime \ast }$,For clarity about contour signs, introduce the inverse-transform measuresThe first Bromwich inversion formula contour runs upward on $\mathrm{Re}p=\sigma$; its conjugate contour runs downward on $\mathrm{Re}q=\sigma$. Both measures therefore have the same positive orientation in their real frequency parameters. Combining the averaged equation with ${\widehat{\delta f}}_k$ yieldsAt a zero of $f_0^{\prime }$, the last expression means its continuous coefficient obtained before cancelling $f_0^{\prime }$. Substitution of the two inverse Laplace transforms and the transformed covariance givesThis is the paper's expression, with the orientation of $d{p}^{\prime \ast }$ made explicit. Reality follows by pairing $k$ and $-k$. The double integral retains switch-on transients; replacing it by a time-independent diffusion coefficient requires the additional long-time limit in b.

In the periodic cube, $H_0=|\mathbf{v}{|}^2/2=(2\pi /L{)}^2|\mathbf{J}{|}^2/2$, so $\mathbf{\Omega }=(2\pi /L)\mathbf{v}$. The supplied action-angle variables are canonical: $d^3\theta d^{3}J=d^{3}x{d}^{3}v$. For $\mathbf{k}=2\pi \mathbf{n}/L$,Consequently c becomesFor each nonzero wavenumber, the Poisson equation for Newtonian gravity givesSubstitute the preceding response and add the externally imposed Newtonian gravitational potential:ThereforeThis gravitational dielectric function has the opposite susceptibility sign to the electron dielectric function, reflecting attractive rather than repulsive collective interactions. The homogeneous mean field is explicitly subtracted in this fictitious periodic model: the Poisson equation for Newtonian gravity acts on density fluctuations, and the $\mathbf{k}=0$ mode is omitted. Otherwise a nonzero uniform positive mass density would be incompatible with a periodic potential. This is the convention often called the Jeans swindle.

For analytic continuation, resolve velocity parallel to $\mathbf{k}$: let $u=\hat{\mathbf{k}}\cdot \mathbf{v}$ and $F(u)=\int f_0(u,{\mathbf{v}}_{\perp })d^2{v}_{\perp }$. ThenStart at $\mathrm{Im}\omega >0$ and deform the Landau contour continuously below the moving pole. At real $\omega$, the integral isFor $\mathrm{Im}\omega <0$, this causal continuation equals the undeformed real-axis integral plus $2\pi i{F}^{\prime }(\omega /k)$ wherever $F^{\prime }$ has the necessary analytic extension. It must not be replaced by the undeformed integral alone.

Zeros of ${ϵ}_{\mathbf{k}}$ are self-consistent collective Landau modes: with no external forcing, a nonzero Newtonian gravitational potential is possible at such a root. For time dependence $e^{-i\omega t}$, upper-half-plane roots grow, lower-half-plane roots describe Landau damping, and real roots are neutral. Damped Landau modes are poles of the continued response; they need not be regular discrete eigenfunctions of the original collisionless operator.

Regularity qualification. Continuation throughout the complex plane presupposes suitable analytic $F$. Arbitrary smooth $F$ need not allow it. For example, $F(u)=C\exp [-1/(1-u^2)]$ for $|u|<1$ and zero elsewhere is smooth but not analytic at $u=\pm 1$; its Cauchy integral cannot be holomorphically continued across those endpoints. The Landau contour gives the standard causal continuation on the domains where the equilibrium permits it, with possible singularities retained.

Take $k>0$ initially; real solutions for $k<0$ follow by $(k,\omega )\mapsto (-k,-\omega )$. Put $V=v_0+v_2-v_1$. Approaching the real phase velocity from the upper half-plane, the Landau contour and the Sokhotski–Plemelj formula give the following Cauchy principal value decompositionThe real frequency condition $ϵ=0$ therefore requires $f_0^{\prime }(u)=0$. This is possible in the interior of the velocity plateau and outside the occupied velocity range:These are candidate intervals, not yet solutions of the real dispersion relation. The endpoints have jumps in $f_0^{\prime }$ and logarithmic singularities in the response; they are not finite-$k$ undamped modes. At the cusp $u=0$, the slope is also discontinuous and the real response diverges, so it supplies no additional mode.

The physical condition is absence of net resonant energy exchange with particles moving at the wave's phase velocity. Outside the support there are no such particles. Within a velocity plateau, the equal populations on either side of the resonance cancel Landau damping.

Let $A_0=f_0(0)$, $\delta =v_2-v_1$, $V=v_0+\delta$ and $K=k^2{\lambda }_{\mathrm{De}}^2$. Use the Debye length convention of the paper, ${\lambda }_{\mathrm{De}}=v_0/{\omega }_{\mathrm{pe}}$. On the candidate intervals in a, the Landau contour contributes no imaginary part. The piecewise constant slopes giveSince $n_0=v_0{A}_0$, setting the dielectric function to zero yieldsThe absolute values express the real response on the allowed intervals. Solving this real equation inside a sloping portion of $f_0$ would not satisfy the imaginary part of the dielectric function.

An exact algebraic form will also make the number of undamped branches clear. DefineOn each allowed interval, the quantity inside the logarithm equals $u^2(u-v_2)/[(u+v_0)(u-v_1)(u-V)]$. Exponentiating and usinggives the cubic equationFor every $k>0$, the logarithm runs from $0$ to $+\infty$ on $(-\infty ,-v_0)$, from $+\infty$ to $-\infty$ on $(v_1,v_2)$, and from $+\infty$ to $0$ on $(V,\infty )$. The intermediate value theorem supplies a root in each interval, and the cubic equation permits at most three. Thus there is exactly one undamped branch in each interval.