Part 4, Chapter 5

Whistler Waves

Right-hand circularly polarized electromagnetic waves in magnetized plasma

5.1 Discovery and Physical Picture

Whistler waves were first detected during World War I as descending-tone audio signals in military telephone lines. The characteristic "whistle" -- a tone that drops in pitch over a fraction of a second -- is caused by lightning-generated electromagnetic waves propagating through the Earth's magnetosphere, where higher frequencies travel faster than lower frequencies.

Whistlers are the right-hand circularly polarized (R-wave) branch of electromagnetic wave propagation parallel to the magnetic field, in the frequency range:

$$\Omega_i \ll \omega \ll \Omega_e$$

where $\Omega_e = eB/m_e$ is the electron cyclotron frequency and $\Omega_i = eB/m_i$ is the ion cyclotron frequency. The wave electric field rotates in the same sense as electron gyration, which is why electrons interact strongly with this mode near the cyclotron resonance.

5.2 Dispersion Relation

For right-hand circularly polarized waves propagating along$\mathbf{B}_0$, the cold-plasma dispersion relation is:

$$\boxed{\frac{k^2 c^2}{\omega^2} = 1 - \frac{\omega_{pe}^2}{\omega(\omega - \Omega_e)}}$$

In the whistler regime ($\omega \ll \Omega_e$ and$\omega \ll \omega_{pe}$), this simplifies to:

$$k^2 c^2 \approx \frac{\omega_{pe}^2 \omega}{\Omega_e}$$

Solving for the frequency:

$$\omega \approx \frac{k^2 c^2 \Omega_e}{\omega_{pe}^2} = k^2 \frac{B}{\mu_0 n_e e}$$

The quadratic dependence $\omega \propto k^2$ leads to a frequency-dependent phase velocity and group velocity:

$$v_\phi = \frac{\omega}{k} = \frac{kc^2\Omega_e}{\omega_{pe}^2} \propto \sqrt{\omega}$$$$v_g = \frac{\partial\omega}{\partial k} = \frac{2kc^2\Omega_e}{\omega_{pe}^2} = 2v_\phi \propto \sqrt{\omega}$$

Both velocities increase with frequency. This is why a broadband lightning pulse is received as a descending tone: the high-frequency components arrive first.

5.3 Magnetospheric Propagation

In the Earth's magnetosphere, whistler waves propagate along magnetic field lines from one hemisphere to the other. The time delay between two frequency components$f_1$ and $f_2$ after traveling a path length $L$ is characterized by the dispersion measure:

$$D = \frac{1}{2c}\int_0^L \frac{\omega_{pe}^2}{\Omega_e^{1/2}}\,ds$$

The arrival time at frequency $f$ follows:

$$t(f) = t_0 + \frac{D}{\sqrt{f}}$$

This $1/\sqrt{f}$ dependence produces the characteristic whistling tone on a spectrogram. By measuring the dispersion, one can infer the integrated electron density along the magnetic field line, making whistlers a natural diagnostic for the plasmasphere.

5.4 Electron Cyclotron Resonance

As the whistler frequency approaches $\Omega_e$, the wave enters the electron cyclotron resonance regime. The refractive index diverges (in the cold-plasma limit), indicating strong wave-particle interaction. Near the resonance:

$$n^2 = \frac{k^2c^2}{\omega^2} \approx \frac{\omega_{pe}^2}{\omega(\Omega_e - \omega)} \to \infty \text{ as } \omega \to \Omega_e$$

In reality, thermal effects prevent the true divergence. Cyclotron damping absorbs the wave energy, heating the electrons. This is exploited in:

Electron cyclotron resonance heating (ECRH): A primary heating method in modern tokamaks and stellarators, using gyrotrons at ~100-170 GHz.
Chorus emissions: Naturally occurring whistler-mode waves near $0.5\Omega_e$ that scatter radiation belt electrons, contributing to satellite damage.

5.5 Applications in Space and Laboratory

Whistler waves are central to several areas of plasma physics research:

Radiation belt dynamics: Whistler-mode chorus and hiss waves are the dominant mechanism for scattering energetic electrons in the Van Allen radiation belts, controlling the lifetime of trapped particles and the radiation environment for spacecraft.
Magnetic reconnection: Whistler waves appear in the electron diffusion region during magnetic reconnection, where they mediate the decoupling of electron and ion motion at scales between the ion and electron inertial lengths.
Helicon plasma sources: In laboratory plasmas, helicon waves (bounded whistler modes) are used to create high-density plasma sources for materials processing and electric propulsion, achieving densities exceeding$10^{19}\;\text{m}^{-3}$.

5.5 Computational Exploration

The code below plots the full whistler dispersion relation and the group velocity as a function of frequency, illustrating the strong frequency dispersion responsible for the whistling tone.

Whistler Waves Simulation

Python

script.py58 lines

import numpy as np
import matplotlib.pyplot as plt

# Parameters (normalized: Omega_e = 1)
Omega_e = 1.0
omega_pe = 5.0  # typical magnetospheric value omega_pe >> Omega_e
c = 1.0

omega = np.linspace(0.001, 0.99 * Omega_e, 1000)

# Full R-wave dispersion: k^2 c^2 / omega^2 = 1 - omega_pe^2 / (omega*(omega - Omega_e))
n_sq = 1.0 - omega_pe**2 / (omega * (omega - Omega_e))
# Only keep propagating solutions (n^2 > 0)
n_sq_prop = np.where(n_sq > 0, n_sq, np.nan)
k_vals = omega * np.sqrt(n_sq_prop) / c

# Group velocity: dw/dk via implicit differentiation
# d(k^2c^2)/domega = d[omega^2 * n^2]/domega
# vg = domega/dk = (2kc^2) / d(omega^2 n^2)/domega
n2_omega2 = omega**2 * n_sq_prop
dn2w2 = np.gradient(n2_omega2, omega)
vg = np.where(np.isfinite(k_vals) & (k_vals > 0), 2 * k_vals * c**2 / dn2w2, np.nan)

# Simplified whistler approximation
omega_approx = np.linspace(0.001, 0.5 * Omega_e, 500)
k_approx = np.sqrt(omega_pe**2 * omega_approx / (c**2 * Omega_e))
vg_approx = 2 * c * np.sqrt(Omega_e * omega_approx) / omega_pe

fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(12, 5))

# Dispersion relation
ax1.plot(k_vals, omega / Omega_e, 'orange', linewidth=2, label='Full R-wave')
ax1.plot(k_approx, omega_approx / Omega_e, '--', color='cyan', linewidth=1.5, label='Whistler approx')
ax1.axhline(y=1.0, color='red', linestyle=':', alpha=0.5, label=r'$\Omega_e$')
ax1.set_xlabel(r'$k$ (normalized)', fontsize=13)
ax1.set_ylabel(r'$\omega / \Omega_e$', fontsize=13)
ax1.set_title('Whistler Dispersion Relation', fontsize=14)
ax1.set_xlim(0, 30)
ax1.set_ylim(0, 1.2)
ax1.legend(fontsize=10)
ax1.grid(True, alpha=0.3)

# Group velocity
ax2.plot(omega / Omega_e, np.abs(vg) / c, 'orange', linewidth=2, label='Full R-wave $v_g$')
ax2.plot(omega_approx / Omega_e, vg_approx / c, '--', color='cyan', linewidth=1.5, label='Whistler approx')
ax2.set_xlabel(r'$\omega / \Omega_e$', fontsize=13)
ax2.set_ylabel(r'$v_g / c$', fontsize=13)
ax2.set_title('Group Velocity vs Frequency', fontsize=14)
ax2.set_xlim(0, 1.0)
ax2.set_ylim(0, 0.25)
ax2.legend(fontsize=10)
ax2.grid(True, alpha=0.3)

plt.tight_layout()
plt.savefig('output.png', dpi=100, bbox_inches='tight', facecolor='#0a0a1a')
plt.close()
print("Whistler wave dispersion and group velocity plotted.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

← Previous Next →