Plasma Etch Uniformity Analysis

Within-Wafer Non-Uniformity (WIWNU) Characterization

Mean Etch Rate
0.00
nm/min
Std Deviation
0.00
nm/min
WIWNU
0.00
%
Edge Exclusion
5.0
mm
3D Wafer Surface with Etch Depth Variation
2D Uniformity Heatmap
Radial Etch Rate Profile
EM Field Intensity Distribution
Gas Flow Velocity Vectors
Temperature Distribution
Uniformity Statistics Histogram

Physics of Etch Uniformity

1. Electromagnetic Standing Wave Patterns

In capacitively coupled plasma (CCP) reactors, the RF electric field creates standing wave patterns that vary with radial position. The electric field distribution follows Bessel function behavior:

E(r,θ) = E₀ × J₀(2πr/λ) × cos(mθ)

where J₀ is the zeroth-order Bessel function, λ is the wavelength at 13.56 MHz (~22m in vacuum, but modified by plasma permittivity), and m represents azimuthal mode number. The radial variation in field intensity directly affects local ionization rates and ion flux to the wafer surface.

2. Sheath Voltage Variation

The voltage drop across the plasma sheath varies radially due to edge effects and non-uniform plasma density. This affects ion bombardment energy:

V_s(r) = V₀ × [1 + α(r/R)²]

where V₀ is the center sheath voltage, R is wafer radius (150mm for 300mm wafers), and α is the edge enhancement factor (typically 0.1-0.3). Higher sheath voltage at the edge increases ion energy and etch rate, contributing to center-fast or edge-fast profiles depending on chemistry.

3. Ion Flux Distribution

Ion flux to the wafer surface exhibits radial non-uniformity due to diffusion and drift processes:

Φ_ion(r) = Φ₀ × exp(-β(r-R)²)

where Φ₀ is the peak flux and β controls the width of the distribution. Edge effects cause flux depletion near the wafer periphery due to recombination at chamber walls and reduced confinement.

4. Gas Residence Time Effects

The residence time of reactive species varies with radial position, affecting local etch chemistry:

τ(r) = V_chamber / Q_flow × f(r)

where V_chamber is the chamber volume, Q_flow is the gas flow rate, and f(r) is a radial distribution function. Longer residence times at the center lead to greater depletion of reactive species, while the edge receives fresher gas supply. This creates center-slow etch profiles in reaction-rate-limited regimes.

5. Temperature Gradient Effects

Wafer temperature varies radially due to non-uniform heating from plasma and cooling from the backside:

T(r) = T_center + ΔT × (r/R)^n

where ΔT is the center-to-edge temperature difference (typically 5-20°C) and n is typically 2-4. Temperature affects etch rate through Arrhenius-type surface reaction kinetics. A 10°C increase can change etch rate by 5-15% depending on activation energy.

6. WIWNU Calculation and Metrics

Within-wafer non-uniformity is quantified as the relative standard deviation:

WIWNU = (σ/μ) × 100%

where σ is the standard deviation of etch rates across all measurement zones and μ is the mean. Advanced nodes require WIWNU < 2% for critical layers. Alternative metrics include range (max-min)/2×mean and site-to-site uniformity for specific measurement patterns (49-point, 121-point).

7. Edge Exclusion Zone Modeling

The outer 3-5mm of the wafer (edge exclusion zone) is typically excluded from uniformity calculations due to extreme edge effects. This region experiences:

8. Electromagnetic Skin Depth

The penetration depth of RF fields into the plasma is limited by skin effect:

δ = c / (ω √(ε_r μ_r))

where c is speed of light, ω is angular frequency (2π × 13.56 MHz), and ε_r is the relative permittivity of the plasma (negative for frequencies below plasma frequency). Skin depth of 1-10 cm means that large-area reactors can have significant field attenuation across the wafer, contributing to non-uniformity.

9. Gas Depletion in High-Density Regions

In regions with high plasma density and etch rate, reactive species can be depleted faster than they are replenished by gas flow. This creates a feedback loop where high-rate regions become reaction-limited. The depletion factor D is given by:

D(r) = 1 - (k_etch × n_wafer × A(r)) / (k_flow × Q_flow)

where k_etch is the etch rate constant, n_wafer is wafer site density, A(r) is the local area element, and k_flow accounts for gas mixing efficiency.

10. Wafer Bow and Thermal Expansion

Wafer bow (curvature) affects gap spacing to the upper electrode, creating non-uniform sheath width:

gap(r) = gap₀ + h_bow × [1 - (r/R)²]

where h_bow is the bow height at center (can be 50-200 μm). Thermal expansion during processing can increase bow by 20-50%, significantly impacting uniformity. Sheath width scales with gap^0.5 to gap^0.75, so even small gap variations translate to measurable etch rate differences.