Semiconductor Manufacturing Process Metrology: Science, Mathematics, and Modeling

A comprehensive exploration of the physics, mathematics, and computational methods underlying nanoscale measurement in semiconductor fabrication.

1. The Fundamental Challenge

Modern semiconductor manufacturing produces structures with critical dimensions of just a few nanometers. At leading-edge nodes (3nm, 2nm), we are measuring features only 10–20 atoms wide.

Key Requirements

• Sub-angstrom precision in measurement
• Complex 3D architectures: FinFETs, Gate-All-Around (GAA) transistors, 3D NAND (200+ layers)
• High throughput: seconds per measurement in production
• Multi-parameter extraction: distinguish dozens of correlated parameters

Metrology Techniques Overview

2. Physics Foundation

2.1 Maxwell's Equations

At the heart of optical metrology lies the solution to Maxwell's equations:

$$

abla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t} $$

$$

abla \times \mathbf{H} = \mathbf{J} + \frac{\partial \mathbf{D}}{\partial t} $$

$$

abla \cdot \mathbf{D} = \rho $$

$$

abla \cdot \mathbf{B} = 0 $$

Where:

• $\mathbf{E}$ = Electric field vector
• $\mathbf{H}$ = Magnetic field vector
• $\mathbf{D}$ = Electric displacement field
• $\mathbf{B}$ = Magnetic flux density
• $\mathbf{J}$ = Current density
• $\rho$ = Charge density

2.2 Constitutive Relations

For linear, isotropic media:

$$ \mathbf{D} = \varepsilon_0 \varepsilon_r \mathbf{E} = \varepsilon_0 (1 + \chi_e) \mathbf{E} $$

$$ \mathbf{B} = \mu_0 \mu_r \mathbf{H} $$

The complex dielectric function:

$$ \tilde{\varepsilon}(\omega) = \varepsilon_1(\omega) + i\varepsilon_2(\omega) = \tilde{n}^2 = (n + ik)^2 $$

Where:

• $n$ = Refractive index
• $k$ = Extinction coefficient

2.3 Fresnel Equations

At an interface between media with refractive indices $\tilde{n}_1$ and $\tilde{n}_2$:

s-polarization (TE):

$$ r_s = \frac{n_1 \cos\theta_i - n_2 \cos\theta_t}{n_1 \cos\theta_i + n_2 \cos\theta_t} $$

$$ t_s = \frac{2 n_1 \cos\theta_i}{n_1 \cos\theta_i + n_2 \cos\theta_t} $$

p-polarization (TM):

$$ r_p = \frac{n_2 \cos\theta_i - n_1 \cos\theta_t}{n_2 \cos\theta_i + n_1 \cos\theta_t} $$

$$ t_p = \frac{2 n_1 \cos\theta_i}{n_2 \cos\theta_i + n_1 \cos\theta_t} $$

With Snell's law:

$$ n_1 \sin\theta_i = n_2 \sin\theta_t $$

3. Mathematics of Inverse Problems

3.1 Problem Formulation

Metrology is fundamentally an inverse problem:

We seek parameters $\mathbf{p}$ that minimize the difference between model $M(\mathbf{p})$ and data $\mathbf{D}$:

$$ \min_{\mathbf{p}} \left\| M(\mathbf{p}) - \mathbf{D} \right\|^2 $$

Or with weighted least squares:

$$ \chi^2 = \sum_{k=1}^{N} \frac{\left( M_k(\mathbf{p}) - D_k \right)^2}{\sigma_k^2} $$

3.2 Levenberg-Marquardt Algorithm

The workhorse optimization algorithm interpolates between gradient descent and Gauss-Newton:

$$ \left( \mathbf{J}^T \mathbf{J} + \lambda \mathbf{I} \right) \delta\mathbf{p} = \mathbf{J}^T \left( \mathbf{D} - M(\mathbf{p}) \right) $$

Where:

• $\mathbf{J}$ = Jacobian matrix (sensitivity matrix)
• $\lambda$ = Damping parameter
• $\delta\mathbf{p}$ = Parameter update step

The Jacobian elements:

$$ J_{ij} = \frac{\partial M_i}{\partial p_j} $$

Algorithm behavior:

• Large $\lambda$ → Gradient descent (robust, slow)
• Small $\lambda$ → Gauss-Newton (fast near minimum)

3.3 Regularization Techniques

For ill-posed problems, regularization is essential:

Tikhonov Regularization (L2):

$$ \min_{\mathbf{p}} \left\| M(\mathbf{p}) - \mathbf{D} \right\|^2 + \alpha \left\| \mathbf{p} - \mathbf{p}_0 \right\|^2 $$

LASSO Regularization (L1):

$$ \min_{\mathbf{p}} \left\| M(\mathbf{p}) - \mathbf{D} \right\|^2 + \alpha \left\| \mathbf{p} \right\|_1 $$

Bayesian Inference:

$$ P(\mathbf{p} | \mathbf{D}) = \frac{P(\mathbf{D} | \mathbf{p}) \cdot P(\mathbf{p})}{P(\mathbf{D})} $$

Where:

• $P(\mathbf{p} | \mathbf{D})$ = Posterior probability
• $P(\mathbf{D} | \mathbf{p})$ = Likelihood
• $P(\mathbf{p})$ = Prior probability

4. Thin Film Optics

4.1 Ellipsometry Fundamentals

Ellipsometry measures the change in polarization state upon reflection:

$$ \rho = \tan(\Psi) \cdot e^{i\Delta} = \frac{r_p}{r_s} $$

Where:

• $\Psi$ = Amplitude ratio angle
• $\Delta$ = Phase difference
• $r_p, r_s$ = Complex reflection coefficients

4.2 Transfer Matrix Method

For multilayer stacks, the characteristic matrix for layer $j$:

$$ \mathbf{M}_j = \begin{pmatrix} \cos\delta_j & \frac{i \sin\delta_j}{\eta_j} \\ i\eta_j \sin\delta_j & \cos\delta_j \end{pmatrix} $$

Where the phase thickness:

$$ \delta_j = \frac{2\pi}{\lambda} \tilde{n}_j d_j \cos\theta_j $$

And the optical admittance:

$$ \eta_j = \begin{cases} \tilde{n}_j \cos\theta_j & \text{(s-pol)} \\ \frac{\tilde{n}_j}{\cos\theta_j} & \text{(p-pol)} \end{cases} $$

Total system matrix:

$$ \mathbf{M}_{total} = \mathbf{M}_1 \cdot \mathbf{M}_2 \cdot \ldots \cdot \mathbf{M}_N = \begin{pmatrix} m_{11} & m_{12} \\ m_{21} & m_{22} \end{pmatrix} $$

Reflection coefficient:

$$ r = \frac{\eta_0 m_{11} + \eta_0 \eta_s m_{12} - m_{21} - \eta_s m_{22}}{\eta_0 m_{11} + \eta_0 \eta_s m_{12} + m_{21} + \eta_s m_{22}} $$

4.3 Dispersion Models

Lorentz Oscillator Model:

$$ \varepsilon(\omega) = \varepsilon_\infty + \sum_j \frac{A_j}{\omega_j^2 - \omega^2 - i\gamma_j \omega} $$

Tauc-Lorentz Model (for amorphous semiconductors):

$$ \varepsilon_2(E) = \begin{cases} \frac{A E_0 C (E - E_g)^2}{(E^2 - E_0^2)^2 + C^2 E^2} \cdot \frac{1}{E} & E > E_g \\ 0 & E \leq E_g \end{cases} $$

With $\varepsilon_1$ obtained via Kramers-Kronig relations:

$$ \varepsilon_1(E) = \varepsilon_{1,\infty} + \frac{2}{\pi} \mathcal{P} \int_{E_g}^{\infty} \frac{\xi \varepsilon_2(\xi)}{\xi^2 - E^2} d\xi $$

5. Scatterometry and RCWA

5.1 Rigorous Coupled-Wave Analysis

For a grating with period $\Lambda$, electromagnetic fields are expanded in Fourier orders:

$$ E(x,z) = \sum_{m=-M}^{M} E_m(z) \exp(i k_{xm} x) $$

Where the diffracted wave vectors:

$$ k_{xm} = k_{x0} + \frac{2\pi m}{\Lambda} = k_0 \left( n_1 \sin\theta_i + \frac{m\lambda}{\Lambda} \right) $$

5.2 Eigenvalue Problem

In each layer, the field satisfies:

$$ \frac{d^2 \mathbf{E}}{dz^2} = \mathbf{\Omega}^2 \mathbf{E} $$

Where $\mathbf{\Omega}^2$ is a matrix determined by the Fourier components of the permittivity:

$$ \varepsilon(x) = \sum_n \varepsilon_n \exp\left( i \frac{2\pi n}{\Lambda} x \right) $$

The eigenvalue decomposition:

$$ \mathbf{\Omega}^2 = \mathbf{W} \mathbf{\Lambda} \mathbf{W}^{-1} $$

Provides propagation constants (eigenvalues $\lambda_m$) and field profiles (eigenvectors in $\mathbf{W}$).

5.3 S-Matrix Formulation

For numerical stability, use the scattering matrix formulation:

$$ \begin{pmatrix} \mathbf{a}_1^- \\ \mathbf{a}_N^+ \end{pmatrix} = \mathbf{S} \begin{pmatrix} \mathbf{a}_1^+ \\ \mathbf{a}_N^- \end{pmatrix} $$

Where $\mathbf{a}^+$ and $\mathbf{a}^-$ represent forward and backward propagating waves.

The S-matrix is built recursively:

$$ \mathbf{S}_{1 \to j+1} = \mathbf{S}_{1 \to j} \star \mathbf{S}_{j,j+1} $$

Using the Redheffer star product $\star$.

6. Statistical Process Control

6.1 Control Charts

$\bar{X}$ Chart (Mean):

$$ UCL = \bar{\bar{X}} + A_2 \bar{R} $$

$$ LCL = \bar{\bar{X}} - A_2 \bar{R} $$

R Chart (Range):

$$ UCL_R = D_4 \bar{R} $$

$$ LCL_R = D_3 \bar{R} $$

EWMA (Exponentially Weighted Moving Average):

$$ Z_t = \lambda X_t + (1 - \lambda) Z_{t-1} $$

With control limits:

$$ UCL = \mu_0 + L \sigma \sqrt{\frac{\lambda}{2 - \lambda} \left[ 1 - (1-\lambda)^{2t} \right]} $$

6.2 Process Capability Indices

$C_p$ (Process Capability):

$$ C_p = \frac{USL - LSL}{6\sigma} $$

$C_{pk}$ (Centered Process Capability):

$$ C_{pk} = \min \left( \frac{USL - \mu}{3\sigma}, \frac{\mu - LSL}{3\sigma} \right) $$

$C_{pm}$ (Taguchi Capability):

$$ C_{pm} = \frac{USL - LSL}{6\sqrt{\sigma^2 + (\mu - T)^2}} $$

Where:

• $USL$ = Upper Specification Limit
• $LSL$ = Lower Specification Limit
• $T$ = Target value
• $\mu$ = Process mean
• $\sigma$ = Process standard deviation

6.3 Gauge R&R Analysis

Total measurement variance decomposition:

$$ \sigma^2_{total} = \sigma^2_{part} + \sigma^2_{gauge} $$

$$ \sigma^2_{gauge} = \sigma^2_{repeatability} + \sigma^2_{reproducibility} $$

Precision-to-Tolerance Ratio:

$$ P/T = \frac{6 \sigma_{gauge}}{USL - LSL} \times 100\% $$

7. Uncertainty Quantification

7.1 Fisher Information Matrix

The Fisher Information Matrix for parameter estimation:

$$ F_{ij} = \sum_{k=1}^{N} \frac{1}{\sigma_k^2} \frac{\partial M_k}{\partial p_i} \frac{\partial M_k}{\partial p_j} $$

Or equivalently:

$$ F_{ij} = -E \left[ \frac{\partial^2 \ln L}{\partial p_i \partial p_j} \right] $$

Where $L$ is the likelihood function.

7.2 Cramér-Rao Lower Bound

The covariance matrix of any unbiased estimator is bounded:

$$ \text{Cov}(\hat{\mathbf{p}}) \geq \mathbf{F}^{-1} $$

For a single parameter:

$$ \text{Var}(\hat{\theta}) \geq \frac{1}{I(\theta)} $$

Interpretation:

• Diagonal elements of $\mathbf{F}^{-1}$ give minimum variance for each parameter
• Off-diagonal elements indicate parameter correlations
• Large condition number of $\mathbf{F}$ indicates ill-conditioning

7.3 Correlation Coefficient

$$ \rho_{ij} = \frac{F^{-1}_{ij}}{\sqrt{F^{-1}_{ii} F^{-1}_{jj}}} $$

7.4 GUM Framework

According to the Guide to the Expression of Uncertainty in Measurement:

Combined standard uncertainty:

$$ u_c^2(y) = \sum_{i=1}^{N} \left( \frac{\partial f}{\partial x_i} \right)^2 u^2(x_i) + 2 \sum_{i=1}^{N-1} \sum_{j=i+1}^{N} \frac{\partial f}{\partial x_i} \frac{\partial f}{\partial x_j} u(x_i, x_j) $$

Expanded uncertainty:

$$ U = k \cdot u_c(y) $$

Where $k$ is the coverage factor (typically $k=2$ for 95% confidence).

8. Machine Learning in Metrology

8.1 Neural Network Surrogate Models

Replace expensive physics simulations with trained neural networks:

$$ M_{NN}(\mathbf{p}; \mathbf{W}) \approx M_{physics}(\mathbf{p}) $$

Training objective:

$$ \mathcal{L} = \frac{1}{N} \sum_{i=1}^{N} \left\| M_{NN}(\mathbf{p}_i) - M_{physics}(\mathbf{p}_i) \right\|^2 + \lambda \left\| \mathbf{W} \right\|^2 $$

Speedup: Typically $10^4$ – $10^6 \times$ faster than RCWA/FEM.

8.2 Physics-Informed Neural Networks (PINNs)

Incorporate physical laws into the loss function:

$$ \mathcal{L}_{total} = \mathcal{L}_{data} + \lambda_{physics} \mathcal{L}_{physics} $$

Where:

$$ \mathcal{L}_{physics} = \left\| abla \times \mathbf{E} + \frac{\partial \mathbf{B}}{\partial t} \right\|^2 + \ldots $$

8.3 Gaussian Process Regression

A non-parametric Bayesian approach:

$$ f(\mathbf{x}) \sim \mathcal{GP}\left( m(\mathbf{x}), k(\mathbf{x}, \mathbf{x}') \right) $$

Common kernel (RBF/Squared Exponential):

$$ k(\mathbf{x}, \mathbf{x}') = \sigma_f^2 \exp\left( -\frac{\left\| \mathbf{x} - \mathbf{x}' \right\|^2}{2\ell^2} \right) $$

Posterior prediction:

$$ \mu_ = \mathbf{k}_^T (\mathbf{K} + \sigma_n^2 \mathbf{I})^{-1} \mathbf{y} $$

$$ \sigma_^2 = k_{**} - \mathbf{k}_^T (\mathbf{K} + \sigma_n^2 \mathbf{I})^{-1} \mathbf{k}_* $$

Advantages:

• Provides uncertainty estimates naturally
• Works well with limited training data
• Interpretable hyperparameters

8.4 Virtual Metrology

Predict wafer properties from equipment sensor data:

$$ \hat{y} = f(FDC_1, FDC_2, \ldots, FDC_n) $$

Where $FDC_i$ are Fault Detection and Classification sensor readings.

Common approaches:

• Partial Least Squares (PLS) regression
• Random Forests
• Gradient Boosting (XGBoost, LightGBM)
• Deep neural networks

9. Advanced Topics and Frontiers

9.1 3D Metrology Challenges

Modern structures require 3D measurement:

9.2 Hybrid Metrology

Combining multiple techniques to break parameter correlations:

$$ \chi^2_{total} = \sum_{techniques} w_t \chi^2_t $$

Example combination:

• OCD for periodic structure parameters
• Ellipsometry for film optical constants
• XRR for density and interface roughness

Mathematical framework:

$$ \mathbf{F}_{hybrid} = \sum_t \mathbf{F}_t $$

Reduces off-diagonal elements, improving condition number.

9.3 Atomic-Scale Considerations

At the 2nm node and beyond:

Line Edge Roughness (LER):

$$ \sigma_{LER} = \sqrt{\frac{1}{L} \int_0^L \left[ x(z) - \bar{x} \right]^2 dz} $$

Power Spectral Density:

$$ PSD(f) = \frac{\sigma^2 \xi}{1 + (2\pi f \xi)^{2(1+H)}} $$

Where:

• $\xi$ = Correlation length
• $H$ = Hurst exponent (roughness character)

Quantum Effects:

• Tunneling through thin barriers
• Discrete dopant effects
• Wave function penetration

9.4 Model-Measurement Circularity

A fundamental epistemological challenge:

-
┌──────────────┐      ┌──────────────┐
│   Physical   │ ───► │   Measured   │
│   Structure  │      │    Signal    │
└──────────────┘      └──────────────┘
       ▲                     │
       │                     ▼
       │              ┌──────────────┐
       │              │    Model     │
       └────────────◄─┤   Inversion  │
                      └──────────────┘

Key questions:

• How do we validate models when "truth" requires modeling?
• Reference metrology (TEM) also requires interpretation
• What does it mean to "know" a dimension at atomic scale?

Key Symbols and Notation

Computational Complexity

Where:

• $N$ = Number of layers / mesh elements
• $M$ = Number of Fourier orders
• $L$ = Number of layers
• $T$ = Number of time steps