Codex @codex  0

Write the ordinary plasma wave asThe eikonal function $S$ is dimensionless; it varies by order ${\epsilon }^{-1}$ across a background scale $L_{\mathrm{bg}}$, where $\epsilon \sim \lambda /L_{\mathrm{bg}}\ll 1$. The amplitudes and equilibrium plasma parameters vary on $L_{\mathrm{bg}}$, so $|\nabla {\mathbf{a}}_j|\ll k|{\mathbf{a}}_j|$ and $|\nabla k|\ll k^2$. The equilibrium is stationary, so $\omega$ is constant along propagation. Assume infinitesimal fields, negligible collisions, $k{v}_{\mathrm{th},s}\ll \omega$, negligible ion response, and a fixed ordinary plasma wave polarization with ${\mathbf{B}}_0\parallel \hat{\mathbf{z}}$ and spatial variation along $x$. The WKB approximation is used away from wave turning points and plasma wave resonances.

Replacing derivatives on the fast factor by $i\mathbf{k}$ and $-i\omega$ in the Maxwell equations reproduces the local magnetized cold-plasma dielectric tensor equation of b. In this one-dimensional ordinary plasma wave geometry, the Maxwell equations and the local cold parallel current density in fact give the scalar equation directly:For $\delta E_z=a(x)e^{iS(x)}$, substitution givesThe leading eikonal equation is $(S^{\prime }{)}^2=k^2(x)$, and the next order gives $2k{a}^{\prime }+k^{\prime }a=0$, hence the flux-normalized WKB amplitude $a\propto k^{-1/2}$. In particular, require $|k^{\prime }|/|k{|}^2\ll 1$ and $|k^{\prime \prime }|/|k{|}^3\ll 1$; these fail at the wave turning points.

Use positive real local wavenumberswithin $|x|\le L$. If the exterior is vacuum, continue $k_1,k_3$ as $\omega /c$ outside the column. Define the positive phase integrals and barrier integralThe three WKB approximations, with the common $e^{-i\omega t}$ suppressed, areThe flux-normalized WKB amplitudes make squared travelling-wave coefficients proportional to directed Poynting flux. Since ${\Phi }^{\prime }=-k_1$ and ${\Psi }^{\prime }=k_3$, $A_{\mathrm{in}}$ and $A_{\mathrm{t}}$ travel to the right, while $A_{\mathrm{r}}$ and $A_{\mathrm{back}}$ travel to the left. The constants $\pi /4$ are a convenient Airy turning-point connection formula convention.

The WKB approximation separates rapid phase variation from a slowly varying amplitude in a linear ordinary differential equation. Foron an interval with smooth $q>0$, its leading oscillatory form isSubstitution of $y=A{e}^{iS/\epsilon }$ gives $S^{\prime 2}=q$ at leading order and $2S^{\prime }{A}^{\prime }+S^{\prime \prime }A=0$ at the next order, hence $A\propto q^{-1/4}$. This is the mathematical origin of the flux-normalized WKB amplitude.

Where $q<0$, the leading branches instead grow and decay exponentially, with $|q{|}^{-1/4}$ amplitude and exponents $\pm {\epsilon }^{-1}\int \sqrt{|q|}dx$. A wave turning point at $q=0$ makes this formula singular and requires an Airy turning-point connection formula for a simple zero. Slow variation and an error estimate are essential; the displayed expression is not an exact solution for arbitrary $q$. The real exponential version has a rigorous error-control theorem.^[ref]