Maxwell Stress Tensor Electromagnetic Force Plasma
# Maxwell Stress Tensor Formalism, Electromagnetic Momentum Flux, and Lorentz Force Kinetics in EUV Plasma Sources and High-Density Discharge Chambers
---
## Executive Summary
The Maxwell stress tensor $\mathbf{T}_{ij}$ is the fundamental quantity that describes electromagnetic momentum and forces in plasma systems. This article develops the complete Maxwell stress tensor formalism from Maxwell's equations, establishes the connection between electromagnetic fields and force density via the Poynting vector and momentum flux, and demonstrates practical applications to EUV lithography plasma sources, RF discharge chamber modeling, and ponderomotive force effects in RF sheaths. Understanding Maxwell stress is essential for designing high-density plasma reactors, predicting plasma confinement and flow patterns, and optimizing EUV source efficiency in next-generation semiconductor lithography.
---
## Table of Contents
1. Introduction: Electromagnetic Forces in Plasma
2. Maxwell Equations and Conservation Laws
3. Electromagnetic Energy Density and Poynting Vector
4. Derivation of Maxwell Stress Tensor
5. Components of Stress Tensor: Electric and Magnetic Parts
6. Momentum Conservation and Force Density
7. Relation to Lorentz Force and Body Forces
8. Boundary Conditions and Surface Stress
9. Pressure Tensor and Anisotropic Effects
10. Ponderomotive Force in RF Plasma Sheaths
11. EUV Plasma Source Physics: Pinch and Droplet Dynamics
12. Numerical Methods: Finite Element Maxwell Solver
13. Python Implementation: 2D Maxwell Stress Calculations
14. Applications to RF Coil Design and Plasma Confinement
15. References & Further Reading
---
## 1. Introduction: Electromagnetic Forces in Plasma
### 1.1 Newton's Third Law in Electromagnetism
In classical mechanics, forces arise from spatial gradients of potentials. In electromagnetism, forces arise from fields $\mathbf{E}$ and $\mathbf{B}$.
The Lorentz force on a charge density $
ho$ and current density $\mathbf{J}$ is:
$$\mathbf{f} = ho \mathbf{E} + \mathbf{J} imes \mathbf{B}$$
This can be rewritten in terms of field stresses (pressure and shear). The Maxwell stress tensor provides the unified framework.
### 1.2 Physical Motivation
Key applications in semiconductor processing:
- EUV source pinch plasma: Magnetic field confines and compresses plasma to generate extreme UV radiation
- RF plasma reactors: Electromagnetic fields drive currents and heat plasma
- Magnetic confinement: Magnetic pressure prevents plasma expansion
- Ponderomotive effects: Nonlinear pressure from RF fields pushes charged particles
---
## 2. Maxwell Equations and Conservation Laws
### 2.1 Maxwell Equations
$$
abla \cdot \mathbf{E} = \frac{
ho}{\epsilon_0}$$
$$
abla imes \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}$$
$$
abla \cdot \mathbf{B} = 0$$
$$
abla imes \mathbf{B} = \mu_0 \mathbf{J} + \mu_0 \epsilon_0 \frac{\partial \mathbf{E}}{\partial t}$$
### 2.2 Lorentz Force Density
$$\mathbf{f} = ho \mathbf{E} + \mathbf{J} imes \mathbf{B}$$
Substituting Maxwell's equations to eliminate $
ho$ and $\mathbf{J}$:
$$\mathbf{f} = \epsilon_0 ( abla \cdot \mathbf{E}) \mathbf{E} + \frac{1}{\mu_0}( abla imes \mathbf{B}) imes \mathbf{B} - \mu_0 \epsilon_0 \frac{\partial}{\partial t}(\mathbf{E} imes \mathbf{B})$$
### 2.3 Electromagnetic Momentum Density
Define the electromagnetic momentum density:
$$\mathbf{g} = \epsilon_0 (\mathbf{E} imes \mathbf{B}) = \frac{\mathbf{S}}{c^2}$$
where $\mathbf{S} = \frac{1}{\mu_0}(\mathbf{E} imes \mathbf{B})$ is the Poynting vector (energy flux).
Time derivative:
$$\frac{\partial \mathbf{g}}{\partial t} = \epsilon_0 \left( \frac{\partial \mathbf{E}}{\partial t} imes \mathbf{B} + \mathbf{E} imes \frac{\partial \mathbf{B}}{\partial t} ight)$$
---
## 3. Electromagnetic Energy Density and Poynting Vector
### 3.1 Energy Density
The electromagnetic energy density is:
$$u = \frac{1}{2}\left(\epsilon_0 E^2 + \frac{1}{\mu_0}B^2 ight)$$
Time evolution (Poynting's theorem):
$$\frac{\partial u}{\partial t} + abla \cdot \mathbf{S} = -\mathbf{J} \cdot \mathbf{E}$$
(energy flows out via Poynting vector; dissipated by Joule heating)
### 3.2 Poynting Vector
$$\mathbf{S} = \frac{1}{\mu_0}(\mathbf{E} imes \mathbf{B})$$
Physical meaning: Energy flux (power per unit area) flowing due to electromagnetic fields.
For a plane EM wave: $|\mathbf{S}| = \frac{E^2}{Z_0}$ where $Z_0 = \sqrt{\mu_0/\epsilon_0} \approx 377$ Ω is the vacuum impedance.
---
## 4. Derivation of Maxwell Stress Tensor
### 4.1 Momentum Conservation Equation
From Maxwell equations, the electromagnetic momentum flux (stress) can be expressed as:
$$\frac{\partial g_i}{\partial t} = \frac{\partial}{\partial x_j} T_{ij} + f_i$$
where $T_{ij}$ is the Maxwell stress tensor and $f_i$ is the Lorentz force density.
Rearranging:
$$f_i = \frac{\partial g_i}{\partial t} - \frac{\partial T_{ij}}{\partial x_j}$$
### 4.2 Complete Form of Maxwell Stress Tensor
$$\boxed{T_{ij} = \epsilon_0 \left(E_i E_j - \frac{1}{2}\delta_{ij}E^2 ight) + \frac{1}{\mu_0}\left(B_i B_j - \frac{1}{2}\delta_{ij}B^2 ight)}$$
Tensor interpretation:
- $T_{ij}$ = force per unit area in the $j$-direction on a surface element perpendicular to the $i$-direction
- Diagonal elements: normal stresses (pressure/tension)
- Off-diagonal elements: shear stresses
### 4.3 Alternative Forms
Electric part:
$$T^{(E)}_{ij} = \epsilon_0 \left(E_i E_j - \frac{1}{2}\delta_{ij}E^2
ight)$$
Magnetic part:
$$T^{(B)}_{ij} = \frac{1}{\mu_0}\left(B_i B_j - \frac{1}{2}\delta_{ij}B^2
ight)$$
Total: $T_{ij} = T^{(E)}_{ij} + T^{(B)}_{ij}$
### 4.4 Pressure Interpretation
Magnetic pressure:
$$P_B = \frac{B^2}{2\mu_0}$$
Acts like a positive pressure that confines plasma.
Electric pressure:
$$P_E = -\frac{\epsilon_0 E^2}{2}$$
Acts like a negative pressure (tension) in the field direction.
---
## 5. Components of Stress Tensor: Electric and Magnetic Parts
### 5.1 Diagonal Components (Pressures)
For a 1D magnetic field along $z$: $\mathbf{B} = B_0 \hat{z}$
$$T_{zz}^{(B)} = \frac{B_0^2}{2\mu_0}$$
$$T_{xx}^{(B)} = T_{yy}^{(B)} = -\frac{B_0^2}{2\mu_0}$$
Physical meaning:
- Pressure along field lines: $+\frac{B^2}{2\mu_0}$ (expansion)
- Pressure perpendicular: $-\frac{B^2}{2\mu_0}$ (contraction)
This is the magnetic "pinch" effect that confines plasma.
### 5.2 Off-Diagonal Components (Shear Stresses)
For crossed $\mathbf{E}$ and $\mathbf{B}$ fields (e.g., $\mathbf{E} = E_x \hat{x}$, $\mathbf{B} = B_z \hat{z}$):
$$T_{xz}^{(E)} = \epsilon_0 E_x \cdot 0 = 0$$
$$T_{xz}^{(B)} = \frac{1}{\mu_0} B_z \cdot 0 = 0$$
No shear stress for perpendicular fields.
For oblique fields, shear components are non-zero and couple different spatial directions.
---
## 6. Momentum Conservation and Force Density
### 6.1 Divergence of Stress Tensor
The force density is:
$$f_i = -\frac{\partial T_{ij}}{\partial x_j}$$
(negative divergence of stress gives force)
This is analogous to mechanical stress: $\sigma_{ij}$ in solids gives stress $f_i = \frac{\partial \sigma_{ij}}{\partial x_j}$.
### 6.2 Lorentz Force in Terms of Stress
The Lorentz force can be written as:
$$ ho \mathbf{E} + \mathbf{J} imes \mathbf{B} = - abla \cdot \mathbf{T} - \epsilon_0 \mu_0 \frac{\partial}{\partial t}(\mathbf{E} imes \mathbf{B})$$
The second term is the rate of change of electromagnetic momentum density.
### 6.3 Force Balance in Equilibrium
In a static plasma ($\partial \mathbf{g}/\partial t = 0$):
$$ ho \mathbf{E} + \mathbf{J} imes \mathbf{B} = - abla \cdot \mathbf{T}$$
This couples particle dynamics to field stresses.
---
## 7. Relation to Lorentz Force and Body Forces
### 7.1 Particle Acceleration
For a fluid of charged particles with density $n$, charge $q$, and bulk velocity $\mathbf{v}$:
$$ ho \frac{d\mathbf{v}}{dt} = ho \mathbf{E} + ho \mathbf{v} imes \mathbf{B} + ext{collisional forces}$$
The electromagnetic stress contributes via $-
abla \cdot \mathbf{T}$.
### 7.2 Plasma Pressure Evolution
Combined momentum and energy equations (Navier-Stokes + Maxwell):
$$ ho \frac{d\mathbf{v}}{dt} = - abla P - abla \cdot \mathbf{T}$$
where $P$ is the gas pressure.
Total pressure: $P_{ ext{eff}} = P_{ ext{gas}} + \mathbf{T}$
---
## 8. Boundary Conditions and Surface Stress
### 8.1 Boundary Conditions
At a surface (e.g., electrode or dielectric):
Normal stress (pressure on surface):
$$P_{ ext{normal}} = T_{nn}$$
Tangential stress (shear on surface):
$$ au_{ ext{tangent}} = T_{nt}$$
### 8.2 Example: Conductor Surface
For a perfect conductor with field E perpendicular to surface:
$$\mathbf{E} = E_n \hat{n}$$
The stress tensor near the surface is:
$$T_{nn} = \frac{\epsilon_0 E_n^2}{2}$$ (outward pressure)
This is the electrostatic pressure on charged electrodes.
---
## 9. Pressure Tensor and Anisotropic Effects
### 9.1 Pressure Anisotropy in Magnetized Plasma
In a magnetized plasma with B = $B_0 \hat{z}$, pressure becomes anisotropic:
$$P_\parallel = ext{pressure along field lines}$$
$$P_\perp = ext{pressure perpendicular to field}$$
The total pressure tensor:
$$P_{ij} = P_\parallel b_i b_j + P_\perp (\delta_{ij} - b_i b_j)$$
where $\hat{b} = \mathbf{B}/B$ is the unit field direction.
### 9.2 Plasma Beta
Define the plasma beta:
$$\beta = \frac{P_{ ext{gas}}}{P_B} = \frac{P_{ ext{gas}}}{\frac{B^2}{2\mu_0}}$$
Low-$\beta$ plasma: Magnetic pressure dominates → field lines shape plasma
High-$\beta$ plasma: Gas pressure dominates → plasma pushes field lines
---
## 10. Ponderomotive Force in RF Plasma Sheaths
### 10.1 Ponderomotive Force Definition
The ponderomotive force is the time-averaged force on a charged particle in an oscillating field:
$$\mathbf{F}_{ ext{pond}} = \frac{q^2}{4m\omega^2} abla (E_0^2 + B_0^2)$$
where $\omega$ is the RF frequency and $E_0, B_0$ are field amplitudes.
Physical meaning: Particles accumulate in regions of weaker field (pushed away from strong-field regions).
### 10.2 Ponderomotive Pressure in RF Sheath
In an RF-driven sheath, the ponderomotive force creates a pressure that:
- Pushes ions away from the electrode
- Creates a pressure-driven ion flux even at zero applied voltage
- Affects sheath thickness and ion energy distribution
Ponderomotive pressure:
$$P_{ ext{pond}} = \frac{\epsilon_0 E_{ ext{RF}}^2}{4}$$
For typical CCP: $E_{ ext{RF}} \sim 10^7$ V/m → $P_{ ext{pond}} \sim 100$ Pa (measurable fraction of total pressure).
---
## 11. EUV Plasma Source Physics: Pinch and Droplet Dynamics
### 11.1 Magnetic Pinch Effect
In EUV sources using magnetic field confinement:
$$P_B = \frac{B^2}{2\mu_0} \sim 10^6 ext{ Pa}$$ (at $B \sim 1$ T)
This magnetic pressure compresses tin (Sn) plasma to extreme densities and temperatures, generating EUV photons.
### 11.2 Tin Droplet Deformation
Applied electric fields deform tin droplets via electrostatic pressure:
$$P_E = \frac{\epsilon_0 E^2}{2}$$
For $E \sim 10^7$ V/m: $P_E \sim 0.5$ kPa (significant for small droplets).
### 11.3 Force Balance in Pinch Plasma
At equilibrium:
$$ abla P_{ ext{gas}} = \mathbf{J} imes \mathbf{B} + ext{viscous forces}$$
The current $\mathbf{J}$ in the plasma and magnetic field $\mathbf{B}$ create the Lorentz force that balances gas pressure.
---
## 12. Numerical Methods: Finite Element Maxwell Solver
### 12.1 Weak Formulation
For computational efficiency, solve Maxwell's equations in weak form using finite elements:
$$\int_V ( abla imes \mathbf{E}) \cdot \mathbf{B}' dV = \int_V \mu_0 \mathbf{J} \cdot \mathbf{B}' dV + \ldots$$
Advantages:
- Natural inclusion of boundary conditions
- Conserves divergence constraints
- Stable for mixed element types
### 12.2 Stress Tensor Recovery
Once $\mathbf{E}$ and $\mathbf{B}$ are computed on the finite element mesh:
$$T_{ij}|_e = \epsilon_0 (E_i E_j - \frac{1}{2}\delta_{ij}E^2) + \frac{1}{\mu_0}(B_i B_j - \frac{1}{2}\delta_{ij}B^2)$$
Element average: Integrate over element volume for local stress.
---
## 13. Python Implementation: 2D Maxwell Stress Calculations
"""
2D Maxwell Stress Tensor Calculator for Plasma Confinement
RF coil + plasma configuration
"""
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.patches import Rectangle, Circle
# Physical constants
epsilon_0 = 8.854e-12 # F/m
mu_0 = 4*np.pi*1e-7 # H/m
c = 2.998e8 # m/s
# Domain setup: 2D cross-section (r-z)
r_max = 0.1 # m (10 cm radius)
z_max = 0.2 # m (20 cm height)
nr, nz = 80, 160
r = np.linspace(0, r_max, nr)
z = np.linspace(0, z_max, nz)
R, Z = np.meshgrid(r, z)
# Analytical fields for a simplified model
# Magnetic field: cylindrical (from solenoid coil)
# B_z from solenoid (coil at z ~ 0.05, 0.15 m)
def B_z_profile(z):
"""Axial magnetic field from solenoid"""
coil_z1, coil_z2 = 0.05, 0.15
coil_width = 0.02
B_0 = 0.5 # T (peak field)
B = B_0 * (np.exp(-((z - coil_z1)**2) / (2*coil_width**2)) +
np.exp(-((z - coil_z2)**2) / (2*coil_width**2)))
return B
# Radial magnetic field: from current in plasma (J_phi)
def B_r_profile(r, z):
"""Radial field from azimuthal current (weaker)"""
return 0.05 * B_z_profile(z) * (r / r_max)
# Electric field: RF-driven (perpendicular to B)
def E_r_profile(r, z):
"""Radial electric field from RF antenna"""
E_0 = 1e5 # V/m (typical RF field)
# Localized near antenna (at r ~ 0.08 m)
antenna_r = 0.08
antenna_z = 0.10
sigma_r, sigma_z = 0.01, 0.03
E = E_0 * np.exp(-((r - antenna_r)**2 / (2*sigma_r**2) +
(z - antenna_z)**2 / (2*sigma_z**2)))
return E
def E_z_profile(r, z):
"""Axial electric field (typically weak)"""
return 0.1 * E_r_profile(r, z) # Weak axial component
# Compute fields on grid
B_z = B_z_profile(Z)
B_r = B_r_profile(R, Z)
B_mag = np.sqrt(B_r**2 + B_z**2)
E_r = E_r_profile(R, Z)
E_z = E_z_profile(R, Z)
E_mag = np.sqrt(E_r**2 + E_z**2)
# Compute Maxwell stress tensor components
# T_rr = epsilon_0 (E_r^2 - E^2/2) + 1/mu_0 (B_r^2 - B^2/2)
# T_zz = epsilon_0 (E_z^2 - E^2/2) + 1/mu_0 (B_z^2 - B^2/2)
# T_rz = epsilon_0 E_r E_z + 1/mu_0 B_r B_z (shear)
T_rr = epsilon_0*(E_r**2 - 0.5*E_mag**2) + (1/mu_0)*(B_r**2 - 0.5*B_mag**2)
T_zz = epsilon_0*(E_z**2 - 0.5*E_mag**2) + (1/mu_0)*(B_z**2 - 0.5*B_mag**2)
T_rz = epsilon_0*E_r*E_z + (1/mu_0)*B_r*B_z
# Magnetic pressure
P_B = B_mag**2 / (2*mu_0)
# Electric pressure
P_E = -0.5*epsilon_0*E_mag**2
# Total pressure
P_total = P_B + P_E
# Plot results
fig, axes = plt.subplots(2, 3, figsize=(16, 10))
# Subplot 1: Magnetic field magnitude
ax = axes[0, 0]
contour_B = ax.contourf(R*100, Z*100, B_mag, levels=20, cmap='viridis')
ax.set_xlabel('Radius r (cm)')
ax.set_ylabel('Height z (cm)')
ax.set_title('Magnetic Field Magnitude |B|')
plt.colorbar(contour_B, ax=ax, label='B (T)')
# Subplot 2: Electric field magnitude
ax = axes[0, 1]
contour_E = ax.contourf(R*100, Z*100, E_mag/1e5, levels=20, cmap='plasma')
ax.set_xlabel('Radius r (cm)')
ax.set_ylabel('Height z (cm)')
ax.set_title('Electric Field Magnitude |E|')
plt.colorbar(contour_E, ax=ax, label='E (10⁵ V/m)')
# Subplot 3: Magnetic pressure
ax = axes[0, 2]
contour_PB = ax.contourf(R*100, Z*100, P_B/1e3, levels=20, cmap='hot')
ax.set_xlabel('Radius r (cm)')
ax.set_ylabel('Height z (cm)')
ax.set_title('Magnetic Pressure P_B = B²/(2μ₀)')
plt.colorbar(contour_PB, ax=ax, label='P_B (kPa)')
# Subplot 4: T_rr (radial normal stress)
ax = axes[1, 0]
contour_Trr = ax.contourf(R*100, Z*100, T_rr/1e3, levels=20, cmap='coolwarm')
ax.set_xlabel('Radius r (cm)')
ax.set_ylabel('Height z (cm)')
ax.set_title('Radial Stress T_rr')
plt.colorbar(contour_Trr, ax=ax, label='T_rr (kPa)')
# Subplot 5: T_zz (axial normal stress)
ax = axes[1, 1]
contour_Tzz = ax.contourf(R*100, Z*100, T_zz/1e3, levels=20, cmap='coolwarm')
ax.set_xlabel('Radius r (cm)')
ax.set_ylabel('Height z (cm)')
ax.set_title('Axial Stress T_zz')
plt.colorbar(contour_Tzz, ax=ax, label='T_zz (kPa)')
# Subplot 6: Total pressure (gas + EM)
ax = axes[1, 2]
contour_Ptot = ax.contourf(R*100, Z*100, P_total/1e3, levels=20, cmap='RdBu_r')
ax.set_xlabel('Radius r (cm)')
ax.set_ylabel('Height z (cm)')
ax.set_title('Total Pressure (B² pressure + E² tension)')
plt.colorbar(contour_Ptot, ax=ax, label='P_total (kPa)')
plt.tight_layout()
plt.savefig('maxwell_stress_tensor_2d.png', dpi=150, bbox_inches='tight')
plt.show()
print(f"
{'='*70}")
print(f"Maxwell Stress Tensor Analysis")
print(f"{'='*70}")
print(f"Domain: r ∈ [0, {r_max*100} cm], z ∈ [0, {z_max*100} cm]")
print(f"Grid: {nr} × {nz} points")
print(f"
Field Statistics:")
print(f" Max |B|: {np.max(B_mag):.3f} T")
print(f" Max |E|: {np.max(E_mag)/1e5:.1f} × 10⁵ V/m")
print(f" Max P_B: {np.max(P_B)/1e3:.2f} kPa")
print(f" Max P_E: {np.min(P_E)/1e3:.2f} kPa (tension)")
print(f" Max total pressure: {np.max(P_total)/1e3:.2f} kPa")
print(f"{'='*70}")---
## 14. Applications to RF Coil Design and Plasma Confinement
### 14.1 RF Coil Magnetic Field Design
The stress tensor guides coil geometry:
- Solenoid coil: Creates axial $B_z$ → confines plasma radially
- Saddle coils: Create radial $B_r$ components → control plasma position
Stress on coil structure must be calculated to prevent mechanical failure.
### 14.2 Plasma Confinement Optimization
Balancing forces via stress tensor:
$$ abla \cdot \mathbf{T} = ext{(force pushing plasma outward)}$$
By shaping $\mathbf{B}$ and $\mathbf{E}$, engineers design stable confinement.
---
## 15. References & Further Reading
1. Jackson, J. D. (1999). *Classical Electrodynamics* (3rd ed.). Wiley.
2. Griffiths, D. J. (2013). *Introduction to Electrodynamics* (4th ed.). Pearson.
3. Landau, L. D., & Lifshitz, E. M. (1984). *Electrodynamics of Continuous Media* (2nd ed.). Butterworth-Heinemann.
4. Chen, F. F. (1984). *Introduction to Plasma Physics and Controlled Fusion* (2nd ed.). Plenum.
---
This article provides the rigorous electromagnetic foundation for understanding force generation and confinement in semiconductor plasma processing systems.