Semiconductor Manufacturing Process: Materials Science & Mathematical Modeling
A comprehensive guide to the physics, chemistry, and mathematics underlying modern semiconductor fabrication.
1. Overview
Modern semiconductor manufacturing is one of the most complex and precise engineering endeavors ever undertaken. Key characteristics include:
- Feature sizes: Leading-edge nodes at 3nm, 2nm, and research into sub-nm
- Precision requirements: Atomic-level control (angstrom tolerances)
- Process steps: Hundreds of sequential operations per chip
- Yield sensitivity: Parts-per-billion defect control
1.1 Core Process Steps
- Crystal Growth
- Czochralski (CZ) process
- Float-zone (FZ) refining
- Epitaxial growth
- Pattern Definition
- Photolithography (DUV, EUV)
- Electron-beam lithography
- Nanoimprint lithography
- Material Addition
- Chemical Vapor Deposition (CVD)
- Physical Vapor Deposition (PVD)
- Atomic Layer Deposition (ALD)
- Epitaxy (MBE, MOCVD)
- Material Removal
- Wet etching (isotropic)
- Dry/plasma etching (anisotropic)
- Chemical Mechanical Polishing (CMP)
- Doping
- Ion implantation
- Thermal diffusion
- Plasma doping
- Thermal Processing
- Oxidation
- Annealing (RTA, spike, laser)
- Silicidation
2. Materials Science Foundations
2.1 Silicon Properties
- Crystal structure: Diamond cubic (Fd3m space group)
- Lattice constant: $a = 5.431 \text{ Å}$
- Bandgap: $E_g = 1.12 \text{ eV}$ (indirect, at 300K)
- Intrinsic carrier concentration:
$$n_i = \sqrt{N_c N_v} \exp\left(-\frac{E_g}{2k_B T}\right)$$
At 300K: $n_i \approx 1.0 \times 10^{10} \text{ cm}^{-3}$
2.2 Crystal Defects
- Point Defects
- Vacancies (V): Missing lattice atoms
- Self-interstitials (I): Extra Si atoms in interstitial sites
- Substitutional impurities: Dopants (B, P, As, Sb)
- Interstitial impurities: Fast diffusers (Fe, Cu, Au)
- Line Defects
- Edge dislocations: Extra half-plane of atoms
- Screw dislocations: Helical atomic arrangement
- Dislocation density target: $< 100 \text{ cm}^{-2}$ for device wafers
- Planar Defects
- Stacking faults: ABCABC → ABCBCABC
- Twin boundaries: Mirror symmetry planes
- Grain boundaries: (avoided in single-crystal wafers)
2.3 Dielectric Materials
| Material | Dielectric Constant ($\kappa$) | Bandgap (eV) | Application |
|----------|-------------------------------|--------------|-------------|
| SiO₂ | 3.9 | 9.0 | Traditional gate oxide |
| Si₃N₄ | 7.5 | 5.3 | Spacers, hard masks |
| HfO₂ | ~25 | 5.8 | High-κ gate dielectric |
| Al₂O₃ | 9 | 8.8 | ALD dielectric |
| ZrO₂ | ~25 | 5.8 | High-κ gate dielectric |
Equivalent Oxide Thickness (EOT):
$$\text{EOT} = t_{\text{high-}\kappa} \cdot \frac{\kappa_{\text{SiO}_2}}{\kappa_{\text{high-}\kappa}} = t_{\text{high-}\kappa} \cdot \frac{3.9}{\kappa_{\text{high-}\kappa}}$$
2.4 Interconnect Materials
- Evolution: Al/SiO₂ → Cu/low-κ → Cu/air-gap → (future: Ru, Co)
- Electromigration - Black's equation for mean time to failure:
$$\text{MTTF} = A \cdot j^{-n} \exp\left(\frac{E_a}{k_B T}\right)$$
Where:
- $j$ = current density
- $n$ ≈ 1-2 (current exponent)
- $E_a$ ≈ 0.7-0.9 eV for Cu
3. Crystal Growth Modeling
3.1 Czochralski Process Physics
The Czochralski process involves pulling a single crystal from a melt. Key phenomena:
- Heat transfer (conduction, convection, radiation)
- Fluid dynamics (buoyancy-driven and forced convection)
- Mass transport (dopant distribution)
- Phase change (solidification at the interface)
3.2 Heat Transfer Equation
$$\rho c_p \frac{\partial T}{\partial t} =
abla \cdot (k
abla T) + Q$$
Where:
- $\rho$ = density [kg/m³]
- $c_p$ = specific heat capacity [J/(kg·K)]
- $k$ = thermal conductivity [W/(m·K)]
- $Q$ = volumetric heat source [W/m³]
3.3 Stefan Problem (Phase Change)
At the solid-liquid interface, the Stefan condition applies:
$$k_s \frac{\partial T_s}{\partial n} - k_\ell \frac{\partial T_\ell}{\partial n} = \rho L v_n$$
Where:
- $k_s$, $k_\ell$ = thermal conductivity of solid and liquid
- $L$ = latent heat of fusion [J/kg]
- $v_n$ = interface velocity normal to the surface [m/s]
3.4 Melt Convection (Navier-Stokes with Boussinesq Approximation)
$$\rho \left( \frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot
abla \mathbf{v} \right) = -
abla p + \mu
abla^2 \mathbf{v} + \rho \mathbf{g} \beta (T - T_0)$$
Dimensionless parameters:
- Grashof number: $Gr = \frac{g \beta \Delta T L^3}{
u^2}$
- Prandtl number: $Pr = \frac{
u}{\alpha}$
- Rayleigh number: $Ra = Gr \cdot Pr$
3.5 Dopant Segregation
Equilibrium segregation coefficient:
$$k_0 = \frac{C_s}{C_\ell}$$
Effective segregation coefficient (Burton-Prim-Slichter model):
$$k_{\text{eff}} = \frac{k_0}{k_0 + (1 - k_0) \exp\left(-\frac{v \delta}{D}\right)}$$
Where:
- $v$ = crystal pull rate [m/s]
- $\delta$ = boundary layer thickness [m]
- $D$ = diffusion coefficient in melt [m²/s]
Dopant concentration along crystal (normal freezing):
$$C_s(f) = k_{\text{eff}} C_0 (1 - f)^{k_{\text{eff}} - 1}$$
Where $f$ = fraction solidified.
4. Diffusion Modeling
4.1 Fick's Laws
First Law (flux proportional to concentration gradient):
$$\mathbf{J} = -D
abla C$$
Second Law (conservation equation):
$$\frac{\partial C}{\partial t} =
abla \cdot (D
abla C)$$
For constant $D$ in 1D:
$$\frac{\partial C}{\partial t} = D \frac{\partial^2 C}{\partial x^2}$$
4.2 Analytical Solutions
Constant surface concentration (predeposition):
$$C(x,t) = C_s \cdot \text{erfc}\left(\frac{x}{2\sqrt{Dt}}\right)$$
Fixed total dose (drive-in):
$$C(x,t) = \frac{Q}{\sqrt{\pi D t}} \exp\left(-\frac{x^2}{4Dt}\right)$$
Where:
- $C_s$ = surface concentration
- $Q$ = total dose [atoms/cm²]
- $\text{erfc}(z) = 1 - \text{erf}(z)$ = complementary error function
4.3 Temperature Dependence
Diffusion coefficient follows Arrhenius behavior:
$$D = D_0 \exp\left(-\frac{E_a}{k_B T}\right)$$
| Dopant | $D_0$ (cm²/s) | $E_a$ (eV) |
|--------|---------------|------------|
| B | 0.76 | 3.46 |
| P | 3.85 | 3.66 |
| As | 0.32 | 3.56 |
| Sb | 0.214 | 3.65 |
4.4 Point-Defect Mediated Diffusion
Dopants diffuse via interactions with point defects. The total diffusivity:
$$D_{\text{eff}} = D_I \frac{C_I}{C_I^} + D_V \frac{C_V}{C_V^}$$
Where:
- $D_I$, $D_V$ = interstitial and vacancy components
- $C_I^$, $C_V^$ = equilibrium concentrations
Coupled defect-dopant equations:
$$\frac{\partial C_I}{\partial t} = D_I
abla^2 C_I + G_I - k_{IV} C_I C_V$$
$$\frac{\partial C_V}{\partial t} = D_V
abla^2 C_V + G_V - k_{IV} C_I C_V$$
Where:
- $G_I$, $G_V$ = generation rates
- $k_{IV}$ = I-V recombination rate constant
4.5 Transient Enhanced Diffusion (TED)
After ion implantation, excess interstitials cause enhanced diffusion:
- "+1" model: Each implanted ion creates ~1 net interstitial
- TED factor: Can enhance diffusion by 10-1000×
- Decay time: τ ~ seconds at high T, hours at low T
5. Ion Implantation
5.1 Range Statistics
Gaussian approximation (light ions, amorphous target):
$$n(x) = \frac{\phi}{\sqrt{2\pi} \Delta R_p} \exp\left(-\frac{(x - R_p)^2}{2 \Delta R_p^2}\right)$$
Where:
- $\phi$ = implant dose [ions/cm²]
- $R_p$ = projected range [nm]
- $\Delta R_p$ = range straggle (standard deviation) [nm]
Pearson IV distribution (heavier ions, includes skewness and kurtosis):
$$n(x) = \frac{\phi}{\Delta R_p} \cdot f\left(\frac{x - R_p}{\Delta R_p}; \gamma, \beta\right)$$
5.2 Stopping Power
Total stopping power (LSS theory):
$$S(E) = -\frac{1}{N}\frac{dE}{dx} = S_n(E) + S_e(E)$$
Where:
- $S_n(E)$ = nuclear stopping (elastic collisions with nuclei)
- $S_e(E)$ = electronic stopping (inelastic interactions with electrons)
- $N$ = atomic density of target
Nuclear stopping (screened Coulomb potential):
$$S_n(E) = \frac{\pi a^2 \gamma E}{1 + M_2/M_1}$$
Where:
- $a$ = screening length
- $\gamma = 4 M_1 M_2 / (M_1 + M_2)^2$
Electronic stopping (velocity-proportional regime):
$$S_e(E) = k_e \sqrt{E}$$
5.3 Monte Carlo Simulation (BCA)
The Binary Collision Approximation treats each collision as isolated:
1. Free flight: Ion travels until next collision
2. Collision: Classical two-body scattering
3. Energy loss: Nuclear + electronic contributions
4. Repeat: Until ion stops ($E < E_{\text{threshold}}$)
Scattering angle (center of mass frame):
$$\theta_{cm} = \pi - 2 \int_{r_{min}}^{\infty} \frac{b \, dr}{r^2 \sqrt{1 - V(r)/E_{cm} - b^2/r^2}}$$
5.4 Damage Accumulation
Kinchin-Pease model for displacement damage:
$$N_d = \frac{0.8 E_d}{2 E_{th}}$$
Where:
- $N_d$ = number of displaced atoms
- $E_d$ = damage energy deposited
- $E_{th}$ = displacement threshold (~15 eV for Si)
Amorphization: Occurs when damage density exceeds ~10% of atomic density
6. Thermal Oxidation
6.1 Deal-Grove Model
The oxide thickness $x$ as a function of time $t$:
$$x^2 + A x = B(t + \tau)$$
Or solved for thickness:
$$x = \frac{A}{2} \left( \sqrt{1 + \frac{4B(t + \tau)}{A^2}} - 1 \right)$$
6.2 Rate Constants
Parabolic rate constant (diffusion-limited):
$$B = \frac{2 D C^*}{N_1}$$
Where:
- $D$ = diffusion coefficient of O₂ in SiO₂
- $C^*$ = equilibrium concentration at surface
- $N_1$ = number of oxidant molecules per unit volume of oxide
Linear rate constant (reaction-limited):
$$\frac{B}{A} = \frac{k_s C^*}{N_1}$$
Where $k_s$ = surface reaction rate constant
6.3 Limiting Cases
Thin oxide ($x \ll A$): Linear regime
$$x \approx \frac{B}{A}(t + \tau)$$
Thick oxide ($x \gg A$): Parabolic regime
$$x \approx \sqrt{B(t + \tau)}$$
6.4 Temperature and Pressure Dependence
$$B = B_0 \exp\left(-\frac{E_B}{k_B T}\right) \cdot \frac{p}{p_0}$$
$$\frac{B}{A} = \left(\frac{B}{A}\right)_0 \exp\left(-\frac{E_{B/A}}{k_B T}\right) \cdot \frac{p}{p_0}$$
| Condition | $E_B$ (eV) | $E_{B/A}$ (eV) |
|-----------|------------|----------------|
| Dry O₂ | 1.23 | 2.0 |
| Wet O₂ (H₂O) | 0.78 | 2.05 |
7. Chemical Vapor Deposition (CVD)
7.1 Reactor Transport Equations
Continuity equation:
$$
abla \cdot (\rho \mathbf{v}) = 0$$
Momentum equation (Navier-Stokes):
$$\rho \left( \frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot
abla \mathbf{v} \right) = -
abla p + \mu
abla^2 \mathbf{v} + \rho \mathbf{g}$$
Energy equation:
$$\rho c_p \left( \frac{\partial T}{\partial t} + \mathbf{v} \cdot
abla T \right) =
abla \cdot (k
abla T) + \sum_i H_i R_i$$
Species transport:
$$\frac{\partial (\rho Y_i)}{\partial t} +
abla \cdot (\rho \mathbf{v} Y_i) =
abla \cdot (\rho D_i
abla Y_i) + M_i \sum_j
u_{ij} r_j$$
Where:
- $Y_i$ = mass fraction of species $i$
- $D_i$ = diffusion coefficient
- $
u_{ij}$ = stoichiometric coefficient
- $r_j$ = reaction rate of reaction $j$
7.2 Surface Reaction Kinetics
Langmuir-Hinshelwood mechanism:
$$R_s = \frac{k_s K_1 K_2 p_1 p_2}{(1 + K_1 p_1 + K_2 p_2)^2}$$
First-order surface reaction:
$$R_s = k_s C_s = k_s \cdot h_m (C_g - C_s)$$
At steady state:
$$C_s = \frac{h_m C_g}{h_m + k_s}$$
7.3 Step Coverage
Thiele modulus for feature filling:
$$\Phi = L \sqrt{\frac{k_s}{D_{\text{Kn}}}}$$
Where:
- $L$ = feature depth
- $D_{\text{Kn}}$ = Knudsen diffusion coefficient
Step coverage behavior:
- $\Phi \ll 1$: Reaction-limited → conformal deposition
- $\Phi \gg 1$: Transport-limited → poor step coverage
7.4 Growth Rate
$$G = \frac{M_f}{\rho_f} \cdot R_s = \frac{M_f}{\rho_f} \cdot \frac{h_m k_s C_g}{h_m + k_s}$$
Where:
- $M_f$ = molecular weight of film
- $\rho_f$ = film density
8. Atomic Layer Deposition (ALD)
8.1 Self-Limiting Surface Reactions
ALD relies on sequential, self-saturating surface reactions.
Surface site model:
$$\frac{d\theta}{dt} = k_{\text{ads}} p (1 - \theta) - k_{\text{des}} \theta$$
At steady state:
$$\theta_{eq} = \frac{K p}{1 + K p}$$
Where $K = k_{\text{ads}} / k_{\text{des}}$ = equilibrium constant
8.2 Growth Per Cycle (GPC)
$$\text{GPC} = \Gamma_{\text{max}} \cdot \theta \cdot \frac{M_f}{\rho_f N_A}$$
Where:
- $\Gamma_{\text{max}}$ = maximum surface site density [sites/cm²]
- $\theta$ = surface coverage (0 to 1)
- $N_A$ = Avogadro's number
Typical GPC values:
- Al₂O₃ (TMA/H₂O): ~1.1 Å/cycle
- HfO₂ (HfCl₄/H₂O): ~1.0 Å/cycle
- TiN (TiCl₄/NH₃): ~0.4 Å/cycle
8.3 Conformality in High Aspect Ratio Features
Penetration depth:
$$\Lambda = \sqrt{\frac{D_{\text{Kn}}}{k_s \Gamma_{\text{max}}}}$$
Conformality factor:
$$\text{CF} = \frac{1}{\sqrt{1 + (L/\Lambda)^2}}$$
For 100% conformality: Require $L \ll \Lambda$
9. Plasma Etching
9.1 Plasma Fundamentals
Electron energy balance:
$$n_e \frac{\partial}{\partial t}\left(\frac{3}{2} k_B T_e\right) =
abla \cdot (\kappa_e
abla T_e) + P_{\text{abs}} - P_{\text{loss}}$$
Debye length (shielding distance):
$$\lambda_D = \sqrt{\frac{\epsilon_0 k_B T_e}{n_e e^2}}$$
Plasma frequency:
$$\omega_{pe} = \sqrt{\frac{n_e e^2}{\epsilon_0 m_e}}$$
9.2 Sheath Physics
Child-Langmuir law (collisionless sheath):
$$J_i = \frac{4 \epsilon_0}{9} \sqrt{\frac{2e}{M_i}} \frac{V_s^{3/2}}{d^2}$$
Where:
- $J_i$ = ion current density
- $V_s$ = sheath voltage
- $d$ = sheath thickness
- $M_i$ = ion mass
Bohm criterion (ion velocity at sheath edge):
$$v_B = \sqrt{\frac{k_B T_e}{M_i}}$$
9.3 Etch Rate Modeling
Ion-enhanced etching:
$$R = R_{\text{chem}} + R_{\text{ion}} = k_n n_{\text{neutral}} + Y \cdot \Gamma_{\text{ion}}$$
Where:
- $R_{\text{chem}}$ = chemical (isotropic) component
- $R_{\text{ion}}$ = ion-enhanced (directional) component
- $Y$ = sputter yield
- $\Gamma_{\text{ion}}$ = ion flux
Anisotropy:
$$A = 1 - \frac{R_{\text{lateral}}}{R_{\text{vertical}}}$$
- $A = 0$: Isotropic
- $A = 1$: Perfectly anisotropic
9.4 Feature-Scale Modeling
Level set equation for surface evolution:
$$\frac{\partial \phi}{\partial t} + F |
abla \phi| = 0$$
Where:
- $\phi(\mathbf{x}, t)$ = level set function
- $F$ = local velocity (etch or deposition rate)
- Surface defined by $\phi = 0$
10. Lithography
10.1 Resolution Limits
Rayleigh criterion:
$$R = k_1 \frac{\lambda}{NA}$$
Depth of focus:
$$DOF = k_2 \frac{\lambda}{NA^2}$$
Where:
- $\lambda$ = wavelength (193 nm DUV, 13.5 nm EUV)
- $NA$ = numerical aperture
- $k_1$, $k_2$ = process-dependent factors
| Technology | λ (nm) | NA | Minimum k₁ | Resolution (nm) |
|------------|--------|-----|------------|-----------------|
| DUV (ArF) | 193 | 1.35 | 0.25 | ~36 |
| EUV | 13.5 | 0.33 | 0.25 | ~10 |
| High-NA EUV | 13.5 | 0.55 | 0.25 | ~6 |
10.2 Aerial Image Formation
Coherent illumination:
$$I(x,y) = \left| \mathcal{F}^{-1} \left\{ \tilde{M}(f_x, f_y) \cdot H(f_x, f_y) \right\} \right|^2$$
Where:
- $\tilde{M}$ = Fourier transform of mask transmission
- $H$ = optical transfer function (pupil function)
Partially coherent illumination (Hopkins formulation):
$$I(x,y) = \iint \iint TCC(f_1, g_1, f_2, g_2) \cdot \tilde{M}(f_1, g_1) \cdot \tilde{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
10.3 Photoresist Chemistry
Chemically Amplified Resists (CARs):
Photoacid generation:
$$\frac{\partial [\text{PAG}]}{\partial t} = -C \cdot I \cdot [\text{PAG}]$$
Acid diffusion and reaction:
$$\frac{\partial [H^+]}{\partial t} = D_H
abla^2 [H^+] + k_{\text{gen}} - k_{\text{neut}}[H^+][Q]$$
Deprotection kinetics:
$$\frac{\partial [M]}{\partial t} = -k_{\text{amp}} [H^+] [M]$$
Where:
- $[\text{PAG}]$ = photoacid generator concentration
- $[H^+]$ = acid concentration
- $[Q]$ = quencher concentration
- $[M]$ = protected site concentration
10.4 Stochastic Effects in EUV
Photon shot noise:
$$\sigma_N = \sqrt{N}$$
Line Edge Roughness (LER):
$$\sigma_{\text{LER}} \propto \frac{1}{\sqrt{\text{dose}}} \propto \frac{1}{\sqrt{N_{\text{photons}}}}$$
Stochastic defect probability:
$$P_{\text{defect}} = 1 - \exp(-\lambda A)$$
Where $\lambda$ = defect density, $A$ = feature area
11. Chemical Mechanical Polishing (CMP)
11.1 Preston Equation
$$\frac{dh}{dt} = K_p \cdot P \cdot v$$
Where:
- $dh/dt$ = material removal rate [nm/s]
- $K_p$ = Preston coefficient [nm/(Pa·m)]
- $P$ = applied pressure [Pa]
- $v$ = relative velocity [m/s]
11.2 Contact Mechanics
Greenwood-Williamson model for asperity contact:
$$A_{\text{real}} = \pi n \beta \sigma \int_{d}^{\infty} (z - d) \phi(z) \, dz$$
$$F = \frac{4}{3} n E^* \sqrt{\beta} \int_{d}^{\infty} (z - d)^{3/2} \phi(z) \, dz$$
Where:
- $n$ = asperity density
- $\beta$ = asperity radius
- $\sigma$ = RMS roughness
- $\phi(z)$ = height distribution
- $E^*$ = effective elastic modulus
11.3 Pattern-Dependent Effects
Dishing (in metal features):
$$\Delta h_{\text{dish}} \propto w^2$$
Where $w$ = line width
Erosion (in dielectric):
$$\Delta h_{\text{erosion}} \propto \rho_{\text{metal}}$$
Where $\rho_{\text{metal}}$ = local metal pattern density
12. Device Simulation (TCAD)
12.1 Poisson Equation
$$
abla \cdot (\epsilon
abla \psi) = -q(p - n + N_D^+ - N_A^-)$$
Where:
- $\psi$ = electrostatic potential [V]
- $\epsilon$ = permittivity
- $n$, $p$ = electron and hole concentrations
- $N_D^+$, $N_A^-$ = ionized donor and acceptor concentrations
12.2 Drift-Diffusion Equations
Current densities:
$$\mathbf{J}_n = q \mu_n n \mathbf{E} + q D_n
abla n$$
$$\mathbf{J}_p = q \mu_p p \mathbf{E} - q D_p
abla p$$
Einstein relation:
$$D_n = \frac{k_B T}{q} \mu_n, \quad D_p = \frac{k_B T}{q} \mu_p$$
Continuity equations:
$$\frac{\partial n}{\partial t} = \frac{1}{q}
abla \cdot \mathbf{J}_n + G - R$$
$$\frac{\partial p}{\partial t} = -\frac{1}{q}
abla \cdot \mathbf{J}_p + G - R$$
12.3 Carrier Statistics
Boltzmann approximation:
$$n = N_c \exp\left(\frac{E_F - E_c}{k_B T}\right)$$
$$p = N_v \exp\left(\frac{E_v - E_F}{k_B T}\right)$$
Fermi-Dirac (degenerate regime):
$$n = N_c \mathcal{F}_{1/2}\left(\frac{E_F - E_c}{k_B T}\right)$$
Where $\mathcal{F}_{1/2}$ = Fermi-Dirac integral of order 1/2
12.4 Recombination Models
Shockley-Read-Hall (SRH):
$$R_{\text{SRH}} = \frac{pn - n_i^2}{\tau_p(n + n_1) + \tau_n(p + p_1)}$$
Auger recombination:
$$R_{\text{Auger}} = (C_n n + C_p p)(pn - n_i^2)$$
Radiative recombination:
$$R_{\text{rad}} = B(pn - n_i^2)$$
13. Advanced Mathematical Methods
13.1 Level Set Methods
Evolution equation:
$$\frac{\partial \phi}{\partial t} + F |
abla \phi| = 0$$
Reinitialization (maintain signed distance function):
$$\frac{\partial \phi}{\partial \tau} = \text{sign}(\phi_0)(1 - |
abla \phi|)$$
Curvature:
$$\kappa =
abla \cdot \left( \frac{
abla \phi}{|
abla \phi|} \right)$$
13.2 Kinetic Monte Carlo (KMC)
Rate catalog:
$$r_i =
u_0 \exp\left(-\frac{E_i}{k_B T}\right)$$
Event selection (Bortz-Kalos-Lebowitz algorithm):
1. Calculate total rate: $R_{\text{tot}} = \sum_i r_i$
2. Generate random $u \in (0,1)$
3. Select event $j$ where $\sum_{i=1}^{j-1} r_i < u \cdot R_{\text{tot}} \leq \sum_{i=1}^{j} r_i$
Time advancement:
$$\Delta t = -\frac{\ln(u')}{R_{\text{tot}}}$$
13.3 Phase Field Methods
Free energy functional:
$$F[\phi] = \int \left[ f(\phi) + \frac{\epsilon^2}{2} |
abla \phi|^2 \right] dV$$
Allen-Cahn equation (non-conserved order parameter):
$$\frac{\partial \phi}{\partial t} = -M \frac{\delta F}{\delta \phi} = M \left[ \epsilon^2
abla^2 \phi - f'(\phi) \right]$$
Cahn-Hilliard equation (conserved order parameter):
$$\frac{\partial \phi}{\partial t} =
abla \cdot \left( M
abla \frac{\delta F}{\delta \phi} \right)$$
13.4 Density Functional Theory (DFT)
Kohn-Sham equations:
$$\left[ -\frac{\hbar^2}{2m}
abla^2 + V_{\text{eff}}(\mathbf{r}) \right] \psi_i(\mathbf{r}) = \epsilon_i \psi_i(\mathbf{r})$$
Effective potential:
$$V_{\text{eff}}(\mathbf{r}) = V_{\text{ext}}(\mathbf{r}) + V_H(\mathbf{r}) + V_{xc}(\mathbf{r})$$
Where:
- $V_{\text{ext}}$ = external (ionic) potential
- $V_H = e^2 \int \frac{n(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|} d\mathbf{r}'$ = Hartree potential
- $V_{xc} = \frac{\delta E_{xc}[n]}{\delta n}$ = exchange-correlation potential
Electron density:
$$n(\mathbf{r}) = \sum_i f_i |\psi_i(\mathbf{r})|^2$$
14. Current Frontiers
14.1 Extreme Ultraviolet (EUV) Lithography
- Challenges:
- Stochastic effects at low photon counts
- Mask defectivity and pellicle development
- Resist trade-offs (sensitivity vs. resolution vs. LER)
- Source power and productivity
- High-NA EUV:
- NA = 0.55 (vs. 0.33 current)
- Anamorphic optics (4× magnification in one direction)
- Sub-8nm half-pitch capability
14.2 3D Integration
- Through-Silicon Vias (TSVs):
- Via-first, via-middle, via-last approaches
- Cu filling and barrier requirements
- Thermal-mechanical stress modeling
- Hybrid Bonding:
- Cu-Cu direct bonding
- Sub-micron alignment requirements
- Surface preparation and activation
14.3 New Materials
- 2D Materials:
- Graphene (zero bandgap)
- Transition metal dichalcogenides (MoS₂, WS₂, WSe₂)
- Hexagonal boron nitride (hBN)
- Wide Bandgap Semiconductors:
- GaN: $E_g = 3.4$ eV
- SiC: $E_g = 3.3$ eV (4H-SiC)
- Ga₂O₃: $E_g = 4.8$ eV
14.4 Novel Device Architectures
- Gate-All-Around (GAA) FETs:
- Nanosheet and nanowire channels
- Superior electrostatic control
- Samsung 3nm, Intel 20A/18A
- Complementary FET (CFET):
- Vertically stacked NMOS/PMOS
- Reduced footprint
- Complex fabrication
- Backside Power Delivery (BSPD):
- Power rails on wafer backside
- Reduced IR drop
- Intel PowerVia
14.5 Machine Learning in Semiconductor Manufacturing
- Virtual Metrology: Predict wafer properties from tool sensor data
- Defect Detection: CNN-based wafer map classification
- Process Optimization: Bayesian optimization, reinforcement learning
- Surrogate Models: Neural networks replacing expensive simulations
- OPC (Optical Proximity Correction): ML-accelerated mask design
Physical Constants
| Constant | Symbol | Value |
|----------|--------|-------|
| Boltzmann constant | $k_B$ | $1.381 \times 10^{-23}$ J/K |
| Elementary charge | $e$ | $1.602 \times 10^{-19}$ C |
| Planck constant | $h$ | $6.626 \times 10^{-34}$ J·s |
| Electron mass | $m_e$ | $9.109 \times 10^{-31}$ kg |
| Permittivity of free space | $\epsilon_0$ | $8.854 \times 10^{-12}$ F/m |
| Avogadro's number | $N_A$ | $6.022 \times 10^{23}$ mol⁻¹ |
| Thermal voltage (300K) | $k_B T/q$ | 25.85 mV |
Multiscale Modeling Hierarchy
| Level | Method | Length Scale | Time Scale | Application |
|-------|--------|--------------|------------|-------------|
| 1 | Ab initio (DFT) | Å | fs | Reaction mechanisms, band structure |
| 2 | Molecular Dynamics | nm | ps-ns | Defect dynamics, interfaces |
| 3 | Kinetic Monte Carlo | nm-μm | ns-s | Growth, etching, diffusion |
| 4 | Continuum (PDE) | μm-mm | s-hr | Process simulation (TCAD) |
| 5 | Compact Models | Device | — | Circuit simulation |
| 6 | Statistical | Die/Wafer | — | Yield prediction |