Part 8, Chapter 4

Laser-Plasma Interaction

High-intensity laser physics and parametric processes

4.1 The Ponderomotive Force

When a charged particle oscillates in a spatially inhomogeneous electromagnetic wave, it experiences a time-averaged force that pushes it away from regions of high field intensity. This is the ponderomotive force, the fundamental nonlinear mechanism in laser-plasma interactions.

Consider an electron oscillating in a laser field with slowly varying envelope. The quiver velocity is $\mathbf{v}_{osc} = e\mathbf{E}/(m_e\omega)$. To second order, the $\mathbf{v} \times \mathbf{B}$ force and the spatial variation of$\mathbf{E}$ produce a time-averaged force:

$$\mathbf{F}_p = -\frac{e^2}{4m_e\omega^2}\nabla E^2 = -\frac{m_e c^2}{4}\nabla a_0^2$$

where $a_0 = eE_0/(m_e\omega c)$ is the normalized vector potential. The ponderomotive potential is:

$$\Phi_p = \frac{e^2 E^2}{4m_e\omega^2} = \frac{m_e c^2}{4} a_0^2$$

For ions, the ponderomotive force is smaller by a factor of $m_e/m_i$ and is usually negligible for optical-frequency lasers. However, for low-frequency or very long-pulse lasers, ion ponderomotive effects can become important (hole boring, radiation pressure acceleration).

4.2 Critical Density and Light Propagation

An electromagnetic wave can only propagate in a plasma if its frequency exceeds the local plasma frequency. The dispersion relation is:

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

The wave becomes evanescent when $\omega = \omega_{pe}$, which defines the critical density:

$$n_c = \frac{m_e \omega^2 \epsilon_0}{e^2} = \frac{1.1 \times 10^{21}}{\lambda_{\mu m}^2} \;\text{cm}^{-3}$$

For a 1 micrometer Nd:glass laser, $n_c \approx 1.1 \times 10^{21}$ cm^-3. For a CO₂ laser at 10.6 micrometers, $n_c \approx 10^{19}$ cm^-3. The region where $n_e = n_c$ is the critical surface, where the laser is reflected or absorbed.

4.3 Stimulated Raman and Brillouin Scattering

A high-intensity laser in a plasma can decay into daughter waves through parametric instabilities. These are three-wave processes satisfying frequency and wavevector matching:

$$\omega_0 = \omega_1 + \omega_2, \qquad \mathbf{k}_0 = \mathbf{k}_1 + \mathbf{k}_2$$

Stimulated Raman Scattering (SRS): The laser (EM wave) decays into a scattered EM wave plus an electron plasma wave (EPW). This requires$n_e < n_c/4$ (so the EPW frequency$\omega_{pe}$ is less than half the laser frequency). The growth rate is:

$$\gamma_{SRS} = \frac{k_{EPW} v_{osc}}{4}\sqrt{\frac{\omega_{pe}}{\omega_0}}$$

Stimulated Brillouin Scattering (SBS): The laser decays into a scattered EM wave plus an ion acoustic wave (IAW). Since the IAW frequency is much lower than the laser frequency, the scattered light is nearly the same frequency as the laser. SBS can occur up to the critical surface:

$$\gamma_{SBS} = \frac{k_{IAW} v_{osc}}{4}\sqrt{\frac{n_e}{n_c}\frac{\omega_0}{\omega_{IAW}}}$$

Both SRS and SBS are critical concerns for inertial confinement fusion, as they scatter laser energy away from the target and can produce energetic electrons (from SRS) that preheat the fuel.

4.4 Self-Focusing and Filamentation

The ponderomotive force expels electrons from the high-intensity center of a laser beam, creating a density channel. Since the refractive index of a plasma is$\eta = \sqrt{1 - n_e/n_c}$, the reduced density on axis increases the refractive index, focusing the beam further. This positive feedback leads to relativistic and ponderomotive self-focusing.

Self-focusing occurs when the laser power exceeds the critical power:

$$P_c = 17\left(\frac{n_c}{n_e}\right) \;\text{GW}$$

For relativistic self-focusing, the electron mass increase$\gamma = \sqrt{1 + a_0^2/2}$ reduces the effective plasma frequency on axis, also increasing the refractive index. Filamentation is the transverse breakup of a wide beam into multiple self-focused filaments, driven by small transverse perturbations in intensity.

Interactive Simulation

Ponderomotive Potential and Critical Density

Python

script.py83 lines

import numpy as np

# Constants
m_e = 9.109e-31
e = 1.602e-19
c = 3e8
eps0 = 8.854e-12
h = 6.626e-34

print("=== Ponderomotive Potential ===")
print()

# Gaussian beam profile: E(r) = E0 * exp(-r^2 / w0^2)
# Ponderomotive potential: Phi_p = e^2 E^2 / (4 m_e omega^2)

wavelength = 1.064e-6  # Nd:glass laser (m)
omega = 2 * np.pi * c / wavelength
w0 = 50e-6   # beam waist (m)

# Intensities to explore
I_values_Wcm2 = [1e14, 1e15, 1e16, 1e17, 1e18]

print(f"Laser wavelength: {wavelength*1e6:.3f} um")
print(f"Beam waist: {w0*1e6:.0f} um")
print(f"omega = {omega:.3e} rad/s")
print()

print(f"{'I (W/cm2)':>12} {'a0':>8} {'E0 (V/m)':>12} {'Phi_p (eV)':>12} {'v_osc/c':>10}")
print("-" * 58)

for I_Wcm2 in I_values_Wcm2:
    I_Wm2 = I_Wcm2 * 1e4
    E0 = np.sqrt(2 * I_Wm2 / (eps0 * c))
    a0 = e * E0 / (m_e * omega * c)
    Phi_p_J = e**2 * E0**2 / (4 * m_e * omega**2)
    Phi_p_eV = Phi_p_J / e
    v_osc = e * E0 / (m_e * omega)
    print(f"{I_Wcm2:12.0e} {a0:8.3f} {E0:12.3e} {Phi_p_eV:12.1f} {v_osc/c:10.4f}")

# Ponderomotive potential profile for I = 1e16 W/cm^2
I0 = 1e16 * 1e4  # W/m^2
E0 = np.sqrt(2 * I0 / (eps0 * c))
r = np.linspace(-3*w0, 3*w0, 80)
E_r = E0 * np.exp(-r**2 / w0**2)
Phi_r = e**2 * E_r**2 / (4 * m_e * omega**2) / e  # in eV

print()
print(f"\nPonderomotive potential profile (I0 = 1e16 W/cm2):")
print(f"{'r/w0':>8} {'Phi_p (eV)':>12}")
print("-" * 24)
for i in range(0, len(r), 8):
    print(f"{r[i]/w0:8.2f} {Phi_r[i]:12.1f}")

# Critical density vs wavelength
print()
print("=== Critical Density vs Laser Wavelength ===")
print()
wavelengths_um = np.array([0.25, 0.35, 0.53, 0.80, 1.064, 3.0, 10.6])
print(f"{'lambda (um)':>12} {'n_c (cm^-3)':>14} {'Application':>25}")
print("-" * 55)
labels = ["4w Nd:glass", "3w Nd:glass", "2w Nd:glass", "Ti:sapphire",
          "Nd:glass", "Mid-IR", "CO2"]
for lam_um, label in zip(wavelengths_um, labels):
    lam_m = lam_um * 1e-6
    om = 2 * np.pi * c / lam_m
    nc = m_e * om**2 * eps0 / e**2  # in m^-3
    nc_cm3 = nc * 1e-6
    print(f"{lam_um:12.3f} {nc_cm3:14.2e} {label:>25}")

# Critical power for self-focusing
print()
print("=== Critical Power for Self-Focusing ===")
n_e_values = [1e19, 1e20, 1e21]  # cm^-3
nc_1um = m_e * (2*np.pi*c/1.064e-6)**2 * eps0 / e**2 * 1e-6  # cm^-3
print(f"n_c(1.064 um) = {nc_1um:.2e} cm^-3")
print()
print(f"{'n_e (cm^-3)':>14} {'n_e/n_c':>8} {'P_c (GW)':>10}")
print("-" * 36)
for ne in n_e_values:
    ratio = ne / nc_1um
    Pc = 17 * (nc_1um / ne)
    print(f"{ne:14.1e} {ratio:8.3f} {Pc:10.1f}")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

← Previous Next →