Plasma is often described as the "fourth state of matter." Physically, it is an ionized gas consisting of free electrons and ions. Unlike a neutral gas, a plasma is dominated by electromagnetic forces. To understand how electromagnetic waves (like light or radio waves) travel through this medium, we must return to the fundamental microscopic form of Maxwell's Equations.
1. Maxwell’s Equations (Microscopic Form)
In a dielectric solid, we often use macroscopic fields (\(\mathbf{D}\) and \(\mathbf{H}\)) to hide the complex atomic behavior. In plasma physics, because the charges are free and mobile, it is clearer to use the vacuum equations and treat the charge density \(\rho\) and current density \(\mathbf{J}\) explicitly.
The governing equations are:
- Gauss's Law: \[ \nabla \cdot \mathbf{E} = \frac{\rho}{\varepsilon_0} \]
- Gauss's Law for Magnetism: \[ \nabla \cdot \mathbf{B} = 0 \]
- Faraday's Law: \[ \nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t} \]
- Ampère–Maxwell Law: \[ \nabla \times \mathbf{B} = \mu_0 \mathbf{J} + \mu_0 \varepsilon_0 \frac{\partial \mathbf{E}}{\partial t} \]
Note: The physics of the plasma enters entirely through how \(\rho\) and \(\mathbf{J}\) respond to the fields \(\mathbf{E}\) and \(\mathbf{B}\).
2. Charge and Current Models
A plasma consists of species \(s\) (electrons \(e\) and ions \(i\)). The total source terms are sums over these species:
\[ \rho = \sum_s q_s n_s \quad \text{and} \quad \mathbf{J} = \sum_s q_s n_s \mathbf{v}_s \]The Quasi-Neutral Approximation
In equilibrium, a plasma is quasi-neutral, meaning the number of electrons roughly equals the number of ions (\(n_e \approx n_i\)).
\[ \rho_{total} \approx 0 \]However, this does not mean electric fields vanish. It implies that charge separation occurs only over very small scales (the Debye length). For electromagnetic waves, the current density \(\mathbf{J}\) is the critical factor.
3. Dynamical Equation: The Equation of Motion
To "close" Maxwell's equations, we need to know how the charges move. The equation of motion is given by the Lorentz Force Law:
\[ m_s \frac{d\mathbf{v}_s}{dt} = q_s (\mathbf{E} + \mathbf{v}_s \times \mathbf{B}) \]Simplifying Assumptions (Cold Plasma Model)
- Electrons dominate: Ions are massive \(m_i \gg m_e\) and respond too slowly to high-frequency waves. We treat ions as a fixed background and focus only on electron motion.
- Linearization: For small oscillations, the velocity \(\mathbf{v}\) is small. The magnetic force term \(\mathbf{v} \times \mathbf{B}\) is a second-order small quantity (product of two small terms) and is neglected.
- Cold Plasma: We ignore thermal motion (pressure gradient terms).
The equation of motion for an electron simplifies to:
\[ m_e \frac{\partial \mathbf{v}}{\partial t} = -e \mathbf{E} \]Deriving the Current Response
Current density for electrons is \(\mathbf{J} = -e n_0 \mathbf{v}\). Differentiating this with respect to time:
\[ \frac{\partial \mathbf{J}}{\partial t} = -e n_0 \frac{\partial \mathbf{v}}{\partial t} \]Substitute the acceleration \(\partial \mathbf{v}/\partial t = -e\mathbf{E}/m_e\):
\[ \frac{\partial \mathbf{J}}{\partial t} = -e n_0 \left( -\frac{e\mathbf{E}}{m_e} \right) = \frac{n_0 e^2}{m_e} \mathbf{E} \]Using the definition of the Plasma Frequency (\(\omega_p\)), where \(\omega_p^2 = \frac{n_0 e^2}{m_e \varepsilon_0}\):
4. Deriving the Wave Equation
Now we combine the electromagnetic laws with the plasma response.
Step A: Take the curl of Faraday’s Law:
\[ \nabla \times (\nabla \times \mathbf{E}) = -\frac{\partial}{\partial t} (\nabla \times \mathbf{B}) \]Step B: Substitute Ampère’s Law and use the vector identity \(\nabla \times (\nabla \times \mathbf{E}) = \nabla(\nabla \cdot \mathbf{E}) - \nabla^2 \mathbf{E}\). Assuming \(\nabla \cdot \mathbf{E} = 0\) (transverse waves):
\[ -\nabla^2 \mathbf{E} = -\mu_0 \frac{\partial \mathbf{J}}{\partial t} - \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2} \]Step C: Rearrange into the standard inhomogeneous wave equation:
\[ \nabla^2 \mathbf{E} - \frac{1}{c^2} \frac{\partial^2 \mathbf{E}}{\partial t^2} = \mu_0 \frac{\partial \mathbf{J}}{\partial t} \]Step D: Substitute the Plasma Current Response derived in Section 3:
\[ \nabla^2 \mathbf{E} - \frac{1}{c^2} \frac{\partial^2 \mathbf{E}}{\partial t^2} = \mu_0 (\varepsilon_0 \omega_p^2 \mathbf{E}) \]Since \(\mu_0 \varepsilon_0 = 1/c^2\):
5. Dispersion Relation and Permittivity
We assume a plane wave solution of the form \(\mathbf{E} \sim e^{i(\mathbf{k}\cdot\mathbf{r} - \omega t)}\).
- Spatial derivative \(\nabla^2 \rightarrow -k^2\)
- Time derivative \(\partial^2/\partial t^2 \rightarrow -\omega^2\)
Substituting these into the wave equation:
\[ -k^2 \mathbf{E} - \frac{1}{c^2}(-\omega^2 \mathbf{E}) = \frac{\omega_p^2}{c^2} \mathbf{E} \] \[ \frac{\omega^2}{c^2} - k^2 = \frac{\omega_p^2}{c^2} \]Multiply by \(c^2\) to get the fundamental Dispersion Relation:
Effective Permittivity \(\varepsilon(\omega)\)
We can rewrite the wave equation to look like a standard dielectric medium by defining a frequency-dependent permittivity. From the dispersion relation \(k^2 = \frac{\omega^2}{c^2} \left( 1 - \frac{\omega_p^2}{\omega^2} \right)\), we identify:
\[ \varepsilon(\omega) = \varepsilon_0 \left( 1 - \frac{\omega_p^2}{\omega^2} \right) \]This single equation explains why plasma behaves like a frequency-selective medium rather than a standard dielectric.
6. Physical Interpretation
The behavior of the wave depends entirely on the frequency \(\omega\) relative to the plasma frequency \(\omega_p\).
Case 1: \(\omega < \omega_p\) (Cutoff / Reflection)
- The term \((1 - \omega_p^2/\omega^2)\) is negative.
- The permittivity \(\varepsilon(\omega)\) is negative.
- The wave number \(k\) becomes imaginary.
- Result: The wave cannot propagate. It decays exponentially (evanescent wave) and is reflected by the plasma. This is why the ionosphere reflects AM radio waves.
Case 2: \(\omega > \omega_p\) (Transmission)
- The permittivity is positive.
- \(k\) is real.
- Result: The wave propagates through the plasma. However, the phase velocity \(v_p = \omega/k\) is greater than \(c\). This is the regime for FM radio, TV, and GPS signals passing through the ionosphere.
7. A Note on Magnetized Plasma
If an external static magnetic field \(\mathbf{B}_0\) is present, the physics becomes significantly more complex. The term \(\mathbf{v} \times \mathbf{B}_0\) can no longer be neglected.
- Electron motion becomes anisotropic (different along vs. across field lines).
- The scalar permittivity \(\varepsilon(\omega)\) becomes a tensor \(\bar{\bar{\varepsilon}}\).
- This splits the wave into two modes: the Ordinary (O) wave and the Extraordinary (X) wave, each with different refractive indices.