Electromagnetic Waves in Magnetized Plasmas & Dielectric Tensor

1. Maxwell's Wave Equation in Plasma

Transverse electromagnetic waves possess oscillating electric and magnetic fields perpendicular to the wavevector ($\vec{k}\cdot\vec{E}_1 = 0, \vec{k}\cdot\vec{B}_1 = 0$). Taking the curl of Faraday's law $\nabla \times \vec{E} = -\frac{\partial \vec{B}}{\partial t}$ and substituting the Ampère-Maxwell law $\nabla \times \vec{B} = \mu_0 \vec{J} + \frac{1}{c^2}\frac{\partial \vec{E}}{\partial t}$:

$$\nabla \times (\nabla \times \vec{E}_1) = -\mu_0 \frac{\partial \vec{J}_1}{\partial t} - \frac{1}{c^2}\frac{\partial^2 \vec{E}_1}{\partial t^2}$$

Using the vector identity $\nabla \times (\nabla \times \vec{E}_1) = \nabla(\nabla\cdot\vec{E}_1) - \nabla^2 \vec{E}_1 = -\nabla^2 \vec{E}_1$ for transverse waves:

$$\nabla^2 \vec{E}_1 - \frac{1}{c^2}\frac{\partial^2 \vec{E}_1}{\partial t^2} = \mu_0 \frac{\partial \vec{J}_1}{\partial t}$$

Assuming harmonic plane waves $\propto e^{i(\vec{k}\cdot\vec{r} - \omega t)}$:

$$-k^2 \vec{E}_1 + \frac{\omega^2}{c^2}\vec{E}_1 = -i\omega \mu_0 \vec{J}_1$$

2. Dispersion Relation & Refractive Index

In a cold unmagnetized plasma, the high-frequency electron current is $\vec{J}_1 = -e n_0 \vec{u}_{e1}$. From the linearized electron momentum equation $m_e(-i\omega)\vec{u}_{e1} = -e\vec{E}_1$:

$$\vec{J}_1 = -e n_0 \left( \frac{e}{i\omega m_e}\vec{E}_1 \right) = i \frac{n_0 e^2}{\omega m_e}\vec{E}_1$$

Substituting $\vec{J}_1$ into the wave equation and using $\mu_0 = 1/(\varepsilon_0 c^2)$:

$$\left( \frac{\omega^2}{c^2} - k^2 \right)\vec{E}_1 = -i\omega \left(\frac{1}{\varepsilon_0 c^2}\right) \left( i\frac{n_0 e^2}{\omega m_e}\vec{E}_1 \right) = \frac{n_0 e^2}{\varepsilon_0 m_e c^2}\vec{E}_1 = \frac{\omega_{pe}^2}{c^2}\vec{E}_1$$

Canceling $\vec{E}_1$ yields the foundational dispersion relation:

$$\omega^2 = \omega_{pe}^2 + c^2 k^2$$

The index of refraction $N \equiv \frac{c k}{\omega} = \frac{c}{v_{ph}}$ is given by:

$$N^2 = 1 - \frac{\omega_{pe}^2}{\omega^2}$$

Phase and group velocities satisfy:

$$v_{ph} = \frac{c}{\sqrt{1 - \omega_{pe}^2/\omega^2}} > c, \quad v_g = \frac{d\omega}{dk} = c \sqrt{1 - \frac{\omega_{pe}^2}{\omega^2}} < c$$ $$v_{ph} v_g = c^2$$

1. Propagation, Cutoff & Evanescence

The index of refraction $N^2 = 1 - \omega_{pe}^2 / \omega^2$ dictates three regimes:

• Overdense / Propagating Regime ($\omega > \omega_{pe}$): $N^2 > 0$, $k$ is real. Transverse electromagnetic waves propagate freely through the plasma.
• Cutoff Condition ($\omega = \omega_{pe}$): $N = 0$, $k = 0$, $v_{ph} \to \infty$. The plasma frequency acts as a strict low-frequency cutoff.
• Underdense / Evanescent Regime ($\omega < \omega_{pe}$): $N^2 < 0$, $k$ is purely imaginary ($k = i\kappa$). The wave cannot propagate: $$\vec{E}_1(x) = \vec{E}_0 e^{-\kappa x} = \vec{E}_0 \exp\left( -\frac{x}{\delta} \right)$$ The wave penetrates only an exponential plasma skin depth $\delta$: $$\delta = \frac{1}{\kappa} = \frac{c}{\sqrt{\omega_{pe}^2 - \omega^2}}$$ Incoming radiation is 100% reflected back at the cutoff boundary.

2. Ionospheric Radio Reflection

The Earth's ionosphere has peak electron densities $n_e \sim 10^{12}\text{ m}^{-3}$, corresponding to a plasma cutoff frequency:

$$f_{pe} \approx 8.98\sqrt{10^{12}}\text{ Hz} \approx 9.0\text{ MHz}$$

High-frequency (HF) shortwave radio signals ($3\text{ to }10\text{ MHz}$) incident on the ionosphere encounter $\omega < \omega_{pe}$ and reflect back to Earth, enabling global over-the-horizon telecommunications. Spacecraft satellite signals (e.g., GPS at $1.5\text{ GHz}$) operate far above the cutoff ($\omega \gg \omega_{pe}$) and transmit cleanly through the ionospheric layer.

1. The Anisotropic Dielectric Tensor

In the presence of an ambient magnetic field $\vec{B}_0 = B_0 \hat{z}$, the plasma response is anisotropic. The induced RF current $\vec{J}_1$ is related to $\vec{E}_1$ via the conductivity tensor $\mathbf{\sigma}$, giving the equivalent dielectric tensor $\mathbf{K}$:

$$\mathbf{K} = \mathbf{I} + \frac{i}{\varepsilon_0 \omega}\mathbf{\sigma} = \begin{pmatrix} S & -i D & 0 \\ i D & S & 0 \\ 0 & 0 & P \end{pmatrix}$$

2. The Stix Notation

Following Thomas Stix's standard notation, the components are defined in terms of right-hand ($R$) and left-hand ($L$) circular response functions:

$$R \equiv 1 - \sum_\alpha \frac{\omega_{p\alpha}^2}{\omega(\omega + \Omega_\alpha)}, \quad L \equiv 1 - \sum_\alpha \frac{\omega_{p\alpha}^2}{\omega(\omega - \Omega_\alpha)}, \quad P \equiv 1 - \sum_\alpha \frac{\omega_{p\alpha}^2}{\omega^2}$$ $$S \equiv \frac{1}{2}(R + L), \quad D \equiv \frac{1}{2}(R - L)$$

where $\Omega_\alpha = q_\alpha B_0 / m_\alpha$ is the signed cyclotron frequency ($\Omega_e = -\omega_{ce} < 0$, $\Omega_i = +\omega_{ci} > 0$).

1. Parallel Wave Equation (k parallel B0)

For wave propagation parallel to the magnetic field ($\vec{k} = k \hat{z} \parallel \vec{B}_0$), Maxwell's equations decouple into two independent circularly polarized eigenmodes:

$$\nabla \times (\nabla \times \vec{E}_1) = \frac{\omega^2}{c^2}\mathbf{K}\cdot\vec{E}_1 \implies \begin{pmatrix} k^2 - \frac{\omega^2}{c^2}S & i\frac{\omega^2}{c^2}D \\ -i\frac{\omega^2}{c^2}D & k^2 - \frac{\omega^2}{c^2}S \end{pmatrix} \begin{pmatrix} E_x \\ E_y \end{pmatrix} = 0$$

Setting the determinant to zero yields the two refractive indices:

$$N^2 = S \pm D$$

2. Right-Hand (R) and Left-Hand (L) Circular Modes

1. Right-Hand Circular Wave (R-wave, $N^2 = R$): $$N_R^2 = 1 - \frac{\omega_{pe}^2}{\omega(\omega - \omega_{ce})}$$ The electric field rotates clockwise, in the same direction as gyrating electrons. As $\omega \to \omega_{ce}$, the denominator vanishes: $N_R^2 \to \infty$, exhibiting electron cyclotron resonance!
2. Left-Hand Circular Wave (L-wave, $N^2 = L$): $$N_L^2 = 1 - \frac{\omega_{pe}^2}{\omega(\omega + \omega_{ce})}$$ The electric field rotates counter-clockwise, in the same direction as gyrating positive ions. It exhibits ion cyclotron resonance at $\omega = \omega_{ci}$.

1. Low-Frequency Right-Hand Branch: Whistler Waves

Consider the R-wave in the intermediate frequency window between the ion and electron cyclotron frequencies:

$$\omega_{ci} \ll \omega \ll \omega_{ce} < \omega_{pe}$$

In this regime, the R-wave index of refraction simplifies dramatically:

$$N_R^2 = 1 - \frac{\omega_{pe}^2}{\omega(\omega - \omega_{ce})} \approx 1 + \frac{\omega_{pe}^2}{\omega(\omega_{ce} - \omega)} \approx \frac{\omega_{pe}^2}{\omega \omega_{ce}}$$

Using $N = c k / \omega$:

$$\frac{c^2 k^2}{\omega^2} \approx \frac{\omega_{pe}^2}{\omega \omega_{ce}} \implies \omega = \frac{c^2 \omega_{ce}}{\omega_{pe}^2} k^2$$

This quadratic dispersion relation ($\omega \propto k^2$) characterizes whistler waves (also called helicons in laboratory solid-state and plasma processing devices).

2. Atmospheric Whistler Chirp Phenomenon

The group velocity of whistler waves is:

$$v_g = \frac{d\omega}{dk} = 2 \frac{c^2 \omega_{ce}}{\omega_{pe}^2} k = 2 c \frac{\sqrt{\omega \omega_{ce}}}{\omega_{pe}} \propto \sqrt{\omega}$$

The group velocity is directly proportional to the square root of frequency: higher frequencies travel faster than lower frequencies!

When lightning strikes the Earth's surface, it emits an impulsive, broadband electromagnetic burst. The pulse couples into the magnetosphere as a whistler wave guided along geomagnetic dipole field lines into the conjugate hemisphere. Because high-frequency components arrive earlier than low-frequency components, an audio receiver detects a distinctive descending tone: a "whistle" descending in pitch over a duration of 1 to 3 seconds.

1. Birefringence of Magnetized Plasma

A linearly polarized electromagnetic wave can be mathematically decomposed into equal-amplitude right-hand and left-hand circularly polarized components:

$$\vec{E}_1 = \frac{E_0}{2}(\hat{x} + i\hat{y})e^{i(k_R z - \omega t)} + \frac{E_0}{2}(\hat{x} - i\hat{y})e^{i(k_L z - \omega t)}$$

Because $N_R \ne N_L$ ($k_R \ne k_L$), the two circular modes propagate with unequal phase velocities. As the wave traverses a distance $d$ through the magnetized plasma, a net phase difference $\Delta\Phi = (k_L - k_R) d$ accumulates between the two components.

2. The Faraday Rotation Angle

Recombining the two modes, the polarization plane of the wave rotates through the Faraday rotation angle $\Delta\theta_F$:

$$\Delta\theta_F = \frac{1}{2}(k_L - k_R) d = \frac{\omega}{2c}\int_0^d (N_L - N_R) dz$$

In the high-frequency limit ($\omega \gg \omega_{ce}, \omega_{pe}$):

$$N_{R,L} \approx 1 - \frac{\omega_{pe}^2}{2\omega^2}\left( 1 \pm \frac{\omega_{ce}}{\omega} \right) \implies N_L - N_R \approx \frac{\omega_{pe}^2 \omega_{ce}}{\omega^3}$$

Substituting $\omega_{pe}^2 = \frac{n_e e^2}{\varepsilon_0 m_e}$ and $\omega_{ce} = \frac{e B_\parallel}{m_e}$:

$$\Delta\theta_F = \frac{e^3}{2\varepsilon_0 m_e^2 c \omega^2} \int_0^d n_e(z) B_\parallel(z) dz = \lambda^2 \left[ \frac{e^3}{8\pi^2 \varepsilon_0 m_e^2 c^3} \int_0^d n_e(z) B_\parallel(z) dz \right]$$

Defining the Rotation Measure (RM):

$$\Delta\theta_F = \text{RM} \cdot \lambda^2, \quad \text{RM} = \frac{e^3}{8\pi^2 \varepsilon_0 m_e^2 c^3} \int_0^d n_e B_\parallel dz$$

Measuring $\Delta\theta_F$ at multiple radio wavelengths $\lambda$ allows astrophysicists to determine the line-of-sight magnetic field of galaxies, pulsars, and the interstellar medium.

1. Perpendicular Propagation (k perp B0)

When wavevector $\vec{k} = k\hat{x}$ is perpendicular to $\vec{B}_0 = B_0 \hat{z}$, two distinct polarization modes emerge:

1. Ordinary Wave (O-mode, $\vec{E}_1 \parallel \vec{B}_0$): The electric field is parallel to the background magnetic field. Because electrons accelerate purely along $\vec{B}_0$, the magnetic Lorentz force vanishes: $\vec{v}_1 \times \vec{B}_0 = 0$. The wave behaves identically to an EM wave in an unmagnetized plasma: $$N_O^2 = P = 1 - \frac{\omega_{pe}^2}{\omega^2}$$ Cutoff occurs at $\omega = \omega_{pe}$; no resonance exists.
2. Extraordinary Wave (X-mode, $\vec{E}_1 \perp \vec{B}_0$): The electric field lies in the $xy$-plane perpendicular to $\vec{B}_0$. The wave is elliptically polarized and drives cyclotron gyration, yielding: $$N_X^2 = \frac{R L}{S} = \frac{(S - D)(S + D)}{S} = \frac{S^2 - D^2}{S}$$ The X-mode exhibits:
   • Resonance ($N_X^2 \to \infty$): When $S = 0$, giving the Upper Hybrid resonance: $$\omega^2 = \omega_{UH}^2 = \omega_{pe}^2 + \omega_{ce}^2$$
   • Cutoffs ($N_X^2 = 0$): When $R = 0$ (Right-hand cutoff $\omega_R$) or $L = 0$ (Left-hand cutoff $\omega_L$): $$\omega_{R,L} = \sqrt{\omega_{pe}^2 + \frac{\omega_{ce}^2}{4}} \pm \frac{\omega_{ce}}{2}$$

2. The Clemmow-Mullaly-Allis (CMA) Diagram

The CMA diagram provides a master topological classification of all electromagnetic wave modes in a cold, magnetized two-fluid plasma by plotting the parameter space:

$$X \equiv \frac{\omega_{pe}^2}{\omega^2} \quad \text{versus} \quad Y \equiv \frac{\omega_{ce}}{\omega}$$

Boundaries in the CMA plane correspond to cutoffs ($R=0, L=0, P=0$) and resonances ($S=0, R=\infty, L=\infty$), partitioning the parameter space into distinct topological regions where wave normal surfaces (phase velocity polar plots) exhibit characteristic propagating or evanescent topologies.