Part 4, Chapter 4

Lower Hybrid Waves

Electrostatic waves coupling ion and electron dynamics at perpendicular propagation

4.1 The Lower Hybrid Frequency

The lower hybrid frequency is a characteristic frequency of magnetized plasmas that lies between the ion and electron cyclotron frequencies. It arises from the coupling of ion and electron dynamics in the direction perpendicular to the magnetic field.

At this frequency, ions undergo unmagnetized oscillations (since $\omega_{LH} \gg \Omega_i$) while electrons remain magnetized (since $\omega_{LH} \ll \Omega_e$). The ions respond to the electric field directly, while electrons drift via$\mathbf{E}\times\mathbf{B}$ and polarization drifts.

The exact lower hybrid frequency is:

$$\boxed{\frac{1}{\omega_{LH}^2} = \frac{1}{\Omega_i\Omega_e + \omega_{pi}^2} + \frac{1}{\omega_{pi}^2 + \Omega_i^2}}$$

which can also be written as:

$$\omega_{LH}^2 = \frac{\omega_{pi}^2}{1 + \omega_{pe}^2/\Omega_e^2}$$

In the common limit where $\omega_{pe} \gg \Omega_e$ (high density), this simplifies to the well-known geometric mean approximation:

$$\omega_{LH} \approx \sqrt{\Omega_e \Omega_i} = \Omega_i\sqrt{\frac{m_i}{m_e}}$$

For a hydrogen plasma in a 1 T field, $\Omega_e \approx 1.76 \times 10^{11}\;\text{rad/s}$ and$\Omega_i \approx 9.58 \times 10^7\;\text{rad/s}$, giving$f_{LH} \approx 650\;\text{MHz}$, which is in the GHz radio frequency range.

4.2 Derivation from Cold Plasma Tensor

For electrostatic waves propagating perpendicular to $\mathbf{B}_0$($\mathbf{k} \perp \mathbf{B}_0$), the dispersion relation comes from setting the perpendicular dielectric function to zero:

$$\epsilon_\perp(\omega) = 1 - \sum_s \frac{\omega_{ps}^2}{\omega^2 - \Omega_s^2} = 0$$

Expanding the sum over electrons and ions:

$$1 - \frac{\omega_{pe}^2}{\omega^2 - \Omega_e^2} - \frac{\omega_{pi}^2}{\omega^2 - \Omega_i^2} = 0$$

This equation has two solutions. The lower hybrid frequency corresponds to the solution satisfying$\Omega_i \ll \omega \ll \Omega_e$. In this regime,$\omega^2 \ll \Omega_e^2$ so the electron term simplifies to$\omega_{pe}^2/\Omega_e^2$, and$\omega^2 \gg \Omega_i^2$ so the ion term simplifies to$\omega_{pi}^2/\omega^2$:

$$1 + \frac{\omega_{pe}^2}{\Omega_e^2} - \frac{\omega_{pi}^2}{\omega^2} = 0 \implies \omega_{LH}^2 = \frac{\omega_{pi}^2}{1 + \omega_{pe}^2/\Omega_e^2}$$

4.3 Applications

Lower hybrid waves have important applications in both laboratory and space plasmas:

Lower hybrid current drive (LHCD): Waves launched at frequencies near $\omega_{LH}$ can resonantly interact with electrons via Landau damping, driving non-inductive current in tokamaks. This is a key tool for steady-state operation in fusion reactors.
Lower hybrid heating: The waves can also heat ions through mode conversion and stochastic ion heating at the lower hybrid resonance layer.
Space physics: Lower hybrid drift instabilities occur at thin current sheets (e.g., the magnetopause) and play a role in magnetic reconnection by providing anomalous resistivity.
Ionospheric heating: High-power radio waves at the lower hybrid frequency can excite parametric instabilities in the ionosphere, used in ionospheric modification experiments.

4.4 Computational Exploration

The code below computes the lower hybrid frequency as a function of magnetic field strength and plasma density, comparing the exact result with the geometric mean approximation.

Lower Hybrid Waves Simulation

Python

script.py71 lines

import numpy as np
import matplotlib.pyplot as plt

# Physical constants
e = 1.602e-19
m_e = 9.109e-31
m_i = 1.673e-27  # proton mass
eps0 = 8.854e-12

# Scan over magnetic field for fixed density
n_e = 1e19  # m^-3 (typical tokamak edge)
B_vals = np.linspace(0.1, 6.0, 200)  # Tesla

omega_pe = np.sqrt(n_e * e**2 / (m_e * eps0))
omega_pi = np.sqrt(n_e * e**2 / (m_i * eps0))

f_LH_exact = []
f_LH_approx = []

for B in B_vals:
    Omega_e = e * B / m_e
    Omega_i = e * B / m_i
    wLH2 = omega_pi**2 / (1 + omega_pe**2 / Omega_e**2)
    f_LH_exact.append(np.sqrt(wLH2) / (2 * np.pi))
    f_LH_approx.append(np.sqrt(Omega_e * Omega_i) / (2 * np.pi))

f_LH_exact = np.array(f_LH_exact)
f_LH_approx = np.array(f_LH_approx)

# Scan over density for fixed B
B_fixed = 3.0  # Tesla
n_vals = np.logspace(17, 21, 200)

Omega_e_f = e * B_fixed / m_e
Omega_i_f = e * B_fixed / m_i

f_LH_dens = []
for n in n_vals:
    wpe = np.sqrt(n * e**2 / (m_e * eps0))
    wpi = np.sqrt(n * e**2 / (m_i * eps0))
    wLH2 = wpi**2 / (1 + wpe**2 / Omega_e_f**2)
    f_LH_dens.append(np.sqrt(wLH2) / (2 * np.pi))

f_LH_dens = np.array(f_LH_dens)

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

# Left: f_LH vs B
ax1.plot(B_vals, f_LH_exact / 1e9, 'orange', linewidth=2, label='Exact')
ax1.plot(B_vals, f_LH_approx / 1e9, '--', color='cyan', linewidth=1.5, label=r'$\sqrt{\Omega_e\Omega_i}/(2\pi)$')
ax1.set_xlabel('Magnetic Field B (T)', fontsize=13)
ax1.set_ylabel('$f_{LH}$ (GHz)', fontsize=13)
ax1.set_title(f'Lower Hybrid Frequency vs B\n($n_e = 10^{{19}}$ m$^{{-3}}$)', fontsize=13)
ax1.legend(fontsize=10)
ax1.grid(True, alpha=0.3)

# Right: f_LH vs density
ax2.loglog(n_vals, f_LH_dens / 1e9, 'orange', linewidth=2, label='Exact')
ax2.axhline(y=np.sqrt(Omega_e_f * Omega_i_f) / (2 * np.pi * 1e9), color='cyan',
            linestyle='--', linewidth=1.5, label=r'$\sqrt{\Omega_e\Omega_i}/(2\pi)$')
ax2.set_xlabel('Electron Density $n_e$ (m$^{-3}$)', fontsize=13)
ax2.set_ylabel('$f_{LH}$ (GHz)', fontsize=13)
ax2.set_title(f'Lower Hybrid Frequency vs Density\n(B = {B_fixed} T)', fontsize=13)
ax2.legend(fontsize=10)
ax2.grid(True, alpha=0.3, which='both')

plt.tight_layout()
plt.savefig('output.png', dpi=100, bbox_inches='tight', facecolor='#0a0a1a')
plt.close()
print("Lower hybrid frequency vs B and density plotted.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

4.5 Finite Temperature and Kinetic Effects

The cold-plasma treatment captures the essential physics, but thermal effects modify the lower hybrid wave in important ways:

Finite $k_\parallel$: When the wave has a component along $\mathbf{B}$, electrons can interact via parallel Landau damping, which is the mechanism behind lower hybrid current drive.
Mode conversion: At the lower hybrid resonance layer in an inhomogeneous plasma, the fast magnetosonic wave converts to a slow electrostatic wave, dramatically enhancing the local wave field and energy deposition.
Perpendicular ion heating: At high wave amplitudes, ion orbits become stochastic when the wave amplitude exceeds a threshold, leading to rapid perpendicular ion heating -- observed in both tokamaks and space plasmas.
Parametric decay: High-power lower hybrid waves can undergo parametric decay into ion acoustic waves and other lower hybrid waves, redistributing energy across the spectrum.

← Previous Next →