advanced mathematics
**Advanced Mathematics in Semiconductor Manufacturing**
**1. Lithography & Optical Physics**
This is arguably the most mathematically demanding area of semiconductor manufacturing.
**1.1 Fourier Optics & Partial Coherence Theory**
The foundation of photolithography treats optical imaging as a spatial frequency filtering problem.
- **Key Concept**: The mask pattern is decomposed into spatial frequency components
- **Optical System**: Acts as a low-pass filter on spatial frequencies
- **Hopkins Formulation**: Describes partially coherent imaging
The aerial image intensity $I(x,y)$ is given by:
$$
I(x,y) = \iint\iint TCC(f_1, g_1, f_2, g_2) \cdot M(f_1, g_1) \cdot M^*(f_2, g_2) \cdot e^{2\pi i[(f_1-f_2)x + (g_1-g_2)y]} \, df_1 \, dg_1 \, df_2 \, dg_2
$$
Where:
- $TCC$ = Transmission Cross-Coefficient
- $M(f,g)$ = Mask spectrum (Fourier transform of mask pattern)
- $M^*$ = Complex conjugate of mask spectrum
**SOCS Decomposition** (Sum of Coherent Systems):
$$
TCC(f_1, g_1, f_2, g_2) = \sum_{k=1}^{N} \lambda_k \phi_k(f_1, g_1) \phi_k^*(f_2, g_2)
$$
- Eigenvalue decomposition makes computation tractable
- $\lambda_k$ are eigenvalues (typically only 10-20 terms needed)
- $\phi_k$ are eigenfunctions
**1.2 Inverse Lithography Technology (ILT)**
Given a desired wafer pattern $T(x,y)$, find the optimal mask $M(x,y)$.
**Mathematical Framework**:
- **Objective Function**:
$$
\min_{M} \left\| I[M](x,y) - T(x,y) \right\|^2 + \alpha R[M]
$$
- **Key Methods**:
- Variational calculus and gradient descent in function spaces
- Level-set methods for topology optimization:
$$
\frac{\partial \phi}{\partial t} + v|
abla\phi| = 0
$$
- Tikhonov regularization: $R[M] = \|
abla M\|^2$
- Total-variation regularization: $R[M] = \int |
abla M| \, dx \, dy$
- Adjoint methods for efficient gradient computation
**1.3 EUV & Rigorous Electromagnetics**
At $\lambda = 13.5$ nm, scalar diffraction theory fails. Full vector Maxwell's equations are required.
**Maxwell's Equations** (time-harmonic form):
$$
abla \times \mathbf{E} = -i\omega\mu\mathbf{H}
$$
$$
abla \times \mathbf{H} = i\omega\varepsilon\mathbf{E}
$$
**Numerical Methods**:
- **RCWA** (Rigorous Coupled-Wave Analysis):
- Eigenvalue problem for each diffraction order
- Transfer matrix for multilayer stacks:
$$
\begin{pmatrix} E^+ \\ E^- \end{pmatrix}_{out} = \mathbf{T} \begin{pmatrix} E^+ \\ E^- \end{pmatrix}_{in}
$$
- **FDTD** (Finite-Difference Time-Domain):
- Yee grid discretization
- Leapfrog time integration:
$$
E^{n+1} = E^n + \frac{\Delta t}{\varepsilon}
abla \times H^{n+1/2}
$$
- **Multilayer Thin-Film Optics**:
- Fresnel coefficients at each interface
- Transfer matrix method for $N$ layers
**1.4 Aberration Theory**
Optical aberrations characterized using **Zernike Polynomials**:
$$
W(\rho, \theta) = \sum_{n,m} Z_n^m R_n^m(\rho) \cdot
\begin{cases}
\cos(m\theta) & \text{(even)} \\
\sin(m\theta) & \text{(odd)}
\end{cases}
$$
Where $R_n^m(\rho)$ are radial polynomials:
$$
R_n^m(\rho) = \sum_{k=0}^{(n-m)/2} \frac{(-1)^k (n-k)!}{k! \left(\frac{n+m}{2}-k\right)! \left(\frac{n-m}{2}-k\right)!} \rho^{n-2k}
$$
**Common Aberrations**:
| Zernike Term | Name | Effect |
|--------------|------|--------|
| $Z_4^0$ | Defocus | Uniform blur |
| $Z_3^1$ | Coma | Asymmetric distortion |
| $Z_4^0$ | Spherical | Halo effect |
| $Z_2^2$ | Astigmatism | Directional blur |
**2. Quantum Mechanics & Device Physics**
As transistors reach sub-5nm dimensions, classical models break down.
**2.1 Schrödinger Equation & Quantum Transport**
**Time-Independent Schrödinger Equation**:
$$
\hat{H}\psi = E\psi
$$
$$
\left[-\frac{\hbar^2}{2m}
abla^2 + V(\mathbf{r})\right]\psi(\mathbf{r}) = E\psi(\mathbf{r})
$$
**Non-Equilibrium Green's Function (NEGF) Formalism**:
- Retarded Green's function:
$$
G^R(E) = \left[(E + i\eta)I - H - \Sigma_L - \Sigma_R\right]^{-1}
$$
- Self-energy $\Sigma$ incorporates:
- Contact coupling
- Scattering mechanisms
- Electron-phonon interaction
- Current calculation:
$$
I = \frac{2e}{h} \int T(E) [f_L(E) - f_R(E)] \, dE
$$
- Transmission function:
$$
T(E) = \text{Tr}\left[\Gamma_L G^R \Gamma_R G^A\right]
$$
**Wigner Function** (bridging quantum and semiclassical):
$$
W(x,p) = \frac{1}{2\pi\hbar} \int \psi^*\left(x + \frac{y}{2}\right) \psi\left(x - \frac{y}{2}\right) e^{ipy/\hbar} \, dy
$$
**2.2 Band Structure Theory**
**$k \cdot p$ Perturbation Theory**:
$$
H_{k \cdot p} = \frac{p^2}{2m_0} + V(\mathbf{r}) + \frac{\hbar}{m_0}\mathbf{k} \cdot \mathbf{p} + \frac{\hbar^2 k^2}{2m_0}
$$
**Effective Mass Tensor**:
$$
\frac{1}{m^*_{ij}} = \frac{1}{\hbar^2} \frac{\partial^2 E}{\partial k_i \partial k_j}
$$
**Tight-Binding Hamiltonian**:
$$
H = \sum_i \varepsilon_i |i\rangle\langle i| + \sum_{\langle i,j \rangle} t_{ij} |i\rangle\langle j|
$$
- $\varepsilon_i$ = on-site energy
- $t_{ij}$ = hopping integral (Slater-Koster parameters)
**2.3 Semiclassical Transport**
**Boltzmann Transport Equation**:
$$
\frac{\partial f}{\partial t} + \mathbf{v} \cdot
abla_r f + \frac{\mathbf{F}}{\hbar} \cdot
abla_k f = \left(\frac{\partial f}{\partial t}\right)_{coll}
$$
- 6D phase space $(x, y, z, k_x, k_y, k_z)$
- Collision integral (scattering):
$$
\left(\frac{\partial f}{\partial t}\right)_{coll} = \sum_{k'} [S(k',k)f(k')(1-f(k)) - S(k,k')f(k)(1-f(k'))]
$$
**Drift-Diffusion Equations** (moment expansion):
$$
\mathbf{J}_n = q\mu_n n\mathbf{E} + qD_n
abla n
$$
$$
\mathbf{J}_p = q\mu_p p\mathbf{E} - qD_p
abla p
$$
**3. Process Simulation PDEs**
**3.1 Dopant Diffusion**
**Fick's Second Law** (concentration-dependent):
$$
\frac{\partial C}{\partial t} =
abla \cdot (D(C,T)
abla C) + G - R
$$
**Coupled Point-Defect System**:
$$
\begin{aligned}
\frac{\partial C_A}{\partial t} &=
abla \cdot (D_A
abla C_A) + k_{AI}C_AC_I - k_{AV}C_AC_V \\
\frac{\partial C_I}{\partial t} &=
abla \cdot (D_I
abla C_I) + G_I - k_{IV}C_IC_V \\
\frac{\partial C_V}{\partial t} &=
abla \cdot (D_V
abla C_V) + G_V - k_{IV}C_IC_V
\end{aligned}
$$
Where:
- $C_A$ = dopant concentration
- $C_I$ = interstitial concentration
- $C_V$ = vacancy concentration
- $k_{ij}$ = reaction rate constants
**3.2 Oxidation & Film Growth**
**Deal-Grove Model**:
$$
x_{ox}^2 + Ax_{ox} = B(t + \tau)
$$
- $A$ = linear rate constant (surface reaction limited)
- $B$ = parabolic rate constant (diffusion limited)
- $\tau$ = time offset for initial oxide
**Moving Boundary (Stefan) Problem**:
$$
D\frac{\partial C}{\partial x}\bigg|_{x=s(t)} = C^* \frac{ds}{dt}
$$
**3.3 Ion Implantation**
**Binary Collision Approximation** (Monte Carlo):
- Screened Coulomb potential:
$$
V(r) = \frac{Z_1 Z_2 e^2}{r} \phi\left(\frac{r}{a}\right)
$$
- Scattering angle from two-body collision integral
**As-Implanted Profile** (Pearson IV distribution):
$$
f(x) = f_0 \left[1 + \left(\frac{x-R_p}{b}\right)^2\right]^{-m} \exp\left[-r \tan^{-1}\left(\frac{x-R_p}{b}\right)\right]
$$
Parameters: $R_p$ (projected range), $\Delta R_p$ (straggle), skewness, kurtosis
**3.4 Plasma Etching**
**Electron Energy Distribution** (Boltzmann equation):
$$
\frac{\partial f}{\partial t} + \mathbf{v} \cdot
abla f - \frac{e\mathbf{E}}{m} \cdot
abla_v f = C[f]
$$
**Child-Langmuir Law** (sheath ion flux):
$$
J = \frac{4\varepsilon_0}{9} \sqrt{\frac{2e}{M}} \frac{V^{3/2}}{d^2}
$$
**3.5 Chemical-Mechanical Polishing (CMP)**
**Preston Equation**:
$$
\frac{dh}{dt} = K_p \cdot P \cdot V
$$
- $K_p$ = Preston coefficient
- $P$ = local pressure
- $V$ = relative velocity
**Pattern-Density Dependent Model**:
$$
P_{local} = P_{avg} \cdot \frac{A_{total}}{A_{contact}(\rho)}
$$
**4. Electromagnetic Simulation**
**4.1 Interconnect Modeling**
**Capacitance Extraction** (Laplace equation):
$$
abla^2 \phi = 0 \quad \text{(dielectric regions)}
$$
$$
abla \cdot (\varepsilon
abla \phi) = -\rho \quad \text{(with charges)}
$$
**Boundary Element Method**:
$$
c(\mathbf{r})\phi(\mathbf{r}) = \int_S \left[\phi(\mathbf{r}') \frac{\partial G}{\partial n'} - G(\mathbf{r}, \mathbf{r}') \frac{\partial \phi}{\partial n'}\right] dS'
$$
Where $G(\mathbf{r}, \mathbf{r}') = \frac{1}{4\pi|\mathbf{r} - \mathbf{r}'|}$ (free-space Green's function)
**4.2 Partial Inductance**
**PEEC Method** (Partial Element Equivalent Circuit):
$$
L_{p,ij} = \frac{\mu_0}{4\pi} \frac{1}{a_i a_j} \int_{V_i} \int_{V_j} \frac{d\mathbf{l}_i \cdot d\mathbf{l}_j}{|\mathbf{r}_i - \mathbf{r}_j|}
$$
**5. Statistical & Stochastic Methods**
**5.1 Process Variability**
**Multivariate Gaussian Model**:
$$
p(\mathbf{x}) = \frac{1}{(2\pi)^{n/2}|\Sigma|^{1/2}} \exp\left(-\frac{1}{2}(\mathbf{x}-\boldsymbol{\mu})^T \Sigma^{-1} (\mathbf{x}-\boldsymbol{\mu})\right)
$$
**Principal Component Analysis**:
$$
\mathbf{X} = \mathbf{U}\mathbf{S}\mathbf{V}^T
$$
- Transform to uncorrelated variables
- Dimensionality reduction: retain components with largest singular values
**Polynomial Chaos Expansion**:
$$
Y(\boldsymbol{\xi}) = \sum_{k=0}^{P} y_k \Psi_k(\boldsymbol{\xi})
$$
- $\Psi_k$ = orthogonal polynomial basis (Hermite for Gaussian inputs)
- Enables uncertainty quantification without Monte Carlo
**5.2 Yield Modeling**
**Poisson Defect Model**:
$$
Y = e^{-D \cdot A}
$$
- $D$ = defect density (defects/cm²)
- $A$ = critical area
**Negative Binomial** (clustered defects):
$$
Y = \left(1 + \frac{DA}{\alpha}\right)^{-\alpha}
$$
**5.3 Reliability Physics**
**Weibull Distribution** (lifetime):
$$
F(t) = 1 - \exp\left[-\left(\frac{t}{\eta}\right)^\beta\right]
$$
- $\eta$ = scale parameter (characteristic life)
- $\beta$ = shape parameter (failure mode indicator)
**Black's Equation** (electromigration):
$$
MTTF = A \cdot J^{-n} \cdot \exp\left(\frac{E_a}{k_B T}\right)
$$
**6. Optimization & Inverse Problems**
**6.1 Design of Experiments**
**Response Surface Methodology**:
$$
y = \beta_0 + \sum_i \beta_i x_i + \sum_i \beta_{ii} x_i^2 + \sum_{i E_g \\
0 & E \leq E_g
\end{cases}
$$
**7. Computational Geometry & Graph Theory**
**7.1 VLSI Physical Design**
**Graph Partitioning** (min-cut):
$$
\min_{P} \sum_{(u,v) \in E : u \in P, v
otin P} w(u,v)
$$
- Kernighan-Lin algorithm
- Spectral methods using Fiedler vector
**Placement** (quadratic programming):
$$
\min_{\mathbf{x}, \mathbf{y}} \sum_{(i,j) \in E} w_{ij} \left[(x_i - x_j)^2 + (y_i - y_j)^2\right]
$$
**Steiner Tree Problem** (routing):
- Given pins to connect, find minimum-length tree
- NP-hard; use approximation algorithms (RSMT, rectilinear Steiner)
**7.2 Mask Data Preparation**
- **Boolean Operations**: Union, intersection, difference of polygons
- **Polygon Clipping**: Sutherland-Hodgman, Vatti algorithms
- **Fracturing**: Decompose complex shapes into trapezoids for e-beam writing
**8. Thermal & Mechanical Analysis**
**8.1 Heat Transport**
**Fourier Heat Equation**:
$$
\rho c_p \frac{\partial T}{\partial t} =
abla \cdot (k
abla T) + Q
$$
**Phonon Boltzmann Transport** (nanoscale):
$$
\frac{\partial f}{\partial t} + \mathbf{v}_g \cdot
abla f = \frac{f_0 - f}{\tau}
$$
- Required when feature size $<$ phonon mean free path
- Non-Fourier effects: ballistic transport, thermal rectification
**8.2 Thermo-Mechanical Stress**
**Linear Elasticity**:
$$
\sigma_{ij} = C_{ijkl} \varepsilon_{kl}
$$
**Equilibrium**:
$$
abla \cdot \boldsymbol{\sigma} + \mathbf{f} = 0
$$
**Thin Film Stress** (Stoney Equation):
$$
\sigma_f = \frac{E_s h_s^2}{6(1-
u_s) h_f} \cdot \frac{1}{R}
$$
- $R$ = wafer curvature radius
- $h_s$, $h_f$ = substrate and film thickness
**Thermal Stress**:
$$
\varepsilon_{thermal} = \alpha \Delta T
$$
$$
\sigma_{thermal} = E(\alpha_{film} - \alpha_{substrate})\Delta T
$$
**9. Multiscale & Atomistic Methods**
**9.1 Molecular Dynamics**
**Equation of Motion**:
$$
m_i \frac{d^2 \mathbf{r}_i}{dt^2} = -
abla_i U(\{\mathbf{r}\})
$$
**Interatomic Potentials**:
- **Tersoff** (covalent, e.g., Si):
$$
V_{ij} = f_c(r_{ij})[f_R(r_{ij}) + b_{ij} f_A(r_{ij})]
$$
- **Embedded Atom Method** (metals):
$$
E_i = F_i(\rho_i) + \frac{1}{2}\sum_{j
eq i} \phi_{ij}(r_{ij})
$$
**Velocity Verlet Integration**:
$$
\mathbf{r}(t+\Delta t) = \mathbf{r}(t) + \mathbf{v}(t)\Delta t + \frac{\mathbf{a}(t)}{2}\Delta t^2
$$
$$
\mathbf{v}(t+\Delta t) = \mathbf{v}(t) + \frac{\mathbf{a}(t) + \mathbf{a}(t+\Delta t)}{2}\Delta t
$$
**9.2 Kinetic Monte Carlo**
**Master Equation**:
$$
\frac{dP_i}{dt} = \sum_j (W_{ji} P_j - W_{ij} P_i)
$$
**Transition Rates** (Arrhenius):
$$
W_{ij} =
u_0 \exp\left(-\frac{E_a}{k_B T}\right)
$$
**BKL Algorithm**:
1. Compute all rates $\{r_i\}$
2. Total rate: $R = \sum_i r_i$
3. Select event $j$ with probability $r_j / R$
4. Advance time: $\Delta t = -\ln(u) / R$ where $u \in (0,1)$
**9.3 Ab Initio Methods**
**Kohn-Sham Equations** (DFT):
$$
\left[-\frac{\hbar^2}{2m}
abla^2 + V_{eff}(\mathbf{r})\right]\psi_i(\mathbf{r}) = \varepsilon_i \psi_i(\mathbf{r})
$$
$$
V_{eff} = V_{ext} + V_H[n] + V_{xc}[n]
$$
Where:
- $V_H[n] = \int \frac{n(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|} d\mathbf{r}'$ (Hartree potential)
- $V_{xc}[n] = \frac{\delta E_{xc}[n]}{\delta n}$ (exchange-correlation)
**10. Machine Learning & Data Science**
**10.1 Virtual Metrology**
**Regression Models**:
- Linear: $y = \mathbf{w}^T \mathbf{x} + b$
- Kernel Ridge Regression:
$$
\mathbf{w} = (\mathbf{K} + \lambda \mathbf{I})^{-1} \mathbf{y}
$$
- Neural Networks: $y = f_L \circ f_{L-1} \circ \cdots \circ f_1(\mathbf{x})$
**10.2 Defect Detection**
**Convolutional Neural Networks**:
$$
(f * g)[n] = \sum_m f[m] \cdot g[n-m]
$$
- Feature extraction through learned filters
- Pooling for translation invariance
**Anomaly Detection**:
- Autoencoders: $\text{loss} = \|x - D(E(x))\|^2$
- Isolation Forest: anomaly score based on path length
**10.3 Process Optimization**
**Bayesian Optimization**:
$$
x_{next} = \arg\max_x \alpha(x | \mathcal{D})
$$
**Acquisition Functions**:
- Expected Improvement: $\alpha_{EI}(x) = \mathbb{E}[\max(f(x) - f^*, 0)]$
- Upper Confidence Bound: $\alpha_{UCB}(x) = \mu(x) + \kappa \sigma(x)$
**Summary**
| Domain | Key Mathematical Topics |
|--------|-------------------------|
| **Lithography** | Fourier analysis, inverse problems, PDEs, optimization |
| **Device Physics** | Quantum mechanics, functional analysis, group theory |
| **Process Simulation** | Nonlinear PDEs, Monte Carlo, stochastic processes |
| **Metrology** | Inverse problems, electromagnetics, statistical inference |
| **Yield/Reliability** | Probability theory, extreme value statistics |
| **Physical Design** | Graph theory, combinatorial optimization, ILP |
| **Thermal/Mechanical** | Continuum mechanics, FEM, tensor analysis |
| **Atomistic Modeling** | Statistical mechanics, DFT, stochastic simulation |
| **Machine Learning** | Neural networks, Bayesian inference, optimization |