Derivation of Dispersion Relations and Wave Velocity

In the study of plasma physics, waves are fundamental to understanding how energy is transported and how instabilities evolve. The mathematical cornerstone of wave theory is the dispersion relation, a functional relationship that links the angular frequency $\omega$ of a wave to its wave vector $\mathbf{k}$. This relationship, typically expressed in the implicit form $D(\omega, \mathbf{k}) = 0$, serves as the bridge between the microscopic particle dynamics and the macroscopic collective behavior of the plasma.

To derive a dispersion relation, one typically begins with a set of governing equations—such as the linearized fluid equations (continuity and momentum equations) or the kinetic Vlasov equation. By assuming a small perturbation from equilibrium and employing the plane-wave ansatz:
[
\delta \psi \propto \exp!\big[i(\mathbf{k}\cdot\mathbf{r} - \omega t)\big]
]
the differential equations are transformed into algebraic equations. The condition for a non-trivial solution to these algebraic equations yields the dispersion relation.

Key Wave Characteristics

Once the dispersion relation is established, two critical velocities define the wave's behavior:

• Phase Velocity ($v_{\mathrm{ph}}$): This represents the speed at which a point of constant phase (such as a crest or a trough) moves through space. It is defined as:
  [
  v_{\mathrm{ph}} = \frac{\omega}{|\mathbf{k}|}
  ]
• Group Velocity ($v_{\mathrm{g}}$): This represents the speed at which the wave envelope (and thus the energy or information) propagates. It is defined as the gradient of the frequency in $\mathbf{k}$-space:
  [
  v_{\mathrm{g}} = \nabla_{\mathbf{k}}\omega(\mathbf{k})
  ]
• Dispersion: A medium is said to be dispersive if the phase velocity depends on the wavenumber $k$. In such media, different spectral components of a wave packet travel at different speeds, causing the packet to spread out or "disperse" over time.

Taxonomy of Common Plasma Waves

Plasma supports a vast spectrum of waves, ranging from high-frequency electromagnetic oscillations to low-frequency magnetohydrodynamic (MHD) modes. The following table summarizes several prototypical dispersion relations:

Methodological Framework for Derivation

The systematic derivation of wave velocities from a dispersion relation follows a rigorous mathematical procedure:

1. Solving the Dispersion Equation: After substituting the plane-wave solution into the linearized system, we obtain $D(\omega, \mathbf{k}) = 0$. For complex systems, this may yield multiple branches of solutions, each corresponding to a different physical mode.
2. Determining Phase Velocity: The phase velocity is calculated directly via $v_{\mathrm{ph}} = \omega/k$. In inhomogeneous plasmas, this is often calculated using a "local approximation" where the medium is assumed to be locally uniform.
3. Determining Group Velocity:
   • If $\omega$ can be isolated as an explicit function of $\mathbf{k}$, we use: $v_{\mathrm{g}} = \frac{\partial \omega}{\partial k}$.
   • If the relation is implicit, we apply the implicit function theorem:
     [
     \frac{\partial D}{\partial \omega} \frac{\partial \omega}{\partial k} + \frac{\partial D}{\partial k} = 0 \implies v_{\mathrm{g}} = -\frac{\partial D / \partial k}{\partial D / \partial \omega}
     ]
4. Analyzing Dispersion Properties: By evaluating $\partial v_{\mathrm{ph}} / \partial k$, we can determine if the wave is non-dispersive ($\partial v_{\mathrm{ph}} / \partial k = 0$) or dispersive.

Detailed Case Studies

Case I: Langmuir Waves in a Warm Plasma

Langmuir waves are high-frequency longitudinal oscillations of electrons. Starting from the linearized electron fluid equations and accounting for thermal pressure, the dispersion relation is:
[
\omega^{2} = \omega_{pe}^{2} + 3k^{2}v_{th}^{2}
]

Phase Velocity Analysis:
[
v_{\mathrm{ph}} = \frac{\omega}{k} = \sqrt{\frac{\omega_{pe}^{2}}{k^{2}} + 3v_{th}^{2}}
]
As $k \to 0$, the phase velocity $v_{\mathrm{ph}} \to \infty$. This indicates that at very long wavelengths, the oscillation is nearly spatially uniform.

Group Velocity Analysis:
[
v_{\mathrm{g}} = \frac{\partial \omega}{\partial k} = \frac{3k v_{th}^{2}}{\sqrt{\omega_{pe}^{2} + 3k^{2}v_{th}^{2}}}
]
In the limit of small wavenumbers ($k \ll \omega_{pe}/v_{th}$), the group velocity simplifies to $v_{\mathrm{g}} \approx 3k v_{th}^{2} / \omega_{pe}$. Here, $v_{\mathrm{g}}$ is proportional to $k$, demonstrating significant dispersion. This explains why a pulse of Langmuir waves will broaden as it propagates through the plasma.

Case II: Electromagnetic Waves in Cold Plasma

In a cold plasma (where thermal effects are neglected), the interaction between Maxwell's equations and electron motion leads to the following relation:
[
\frac{k^{2}c^{2}}{\omega^{2}} = 1 - \frac{\omega_{pe}^{2}}{\omega^{2}}
]

Phase and Group Velocity Relationship:
The phase velocity is given by:
[
v_{\mathrm{ph}} = \frac{\omega}{k} = \frac{c}{\sqrt{1 - \omega_{pe}^{2}/\omega^{2}}}
]
Note that as $\omega$ approaches the cutoff frequency $\omega_{pe}$ from above, $v_{\mathrm{ph}} \to \infty$.

The group velocity is:
[
v_{\mathrm{g}} = \frac{\partial \omega}{\partial k} = c\sqrt{1 - \frac{\omega_{pe}^{2}}{\omega^{2}}}
]
An elegant result emerges here: $v_{\mathrm{ph}} \cdot v_{\mathrm{g}} = c^2$. This relationship is characteristic of electromagnetic waves in dispersive media. Furthermore, if $\omega < \omega_{pe}$, the wavenumber $k$ becomes imaginary, meaning the wave cannot propagate and instead becomes an evanescent wave.

Summary and Practical Implications

The derivation and analysis of dispersion relations are indispensable tools in plasma science. While the phase velocity is critical for understanding phase-matching conditions in nonlinear processes and antenna design, the group velocity is the primary determinant for energy transport, wave-packet evolution, and the study of plasma instabilities.

By mastering the transition from linearized fluid/kinetic equations to the explicit calculation of $v_{\mathrm{ph}}$ and $v_{\mathrm{g}}$, researchers can accurately predict how waves will behave in complex environments, ranging from laboratory fusion devices to the vast reaches of the solar corona.