Semiconductor Manufacturing Process Thermal Dynamics
1. Introduction and Fundamental Importance
Thermal dynamics govern nearly every step in semiconductor fabrication. Temperature control determines chemical reaction rates, diffusion velocities, film properties, stress states, and ultimately device performance.
1.1 The Arrhenius Relationship
The fundamental equation governing thermally-activated processes:
$$
k = A \cdot e^{-\frac{E_a}{k_B T}}
$$
Where:
- $k$ = reaction rate constant
- $A$ = pre-exponential factor (frequency factor)
- $E_a$ = activation energy (eV or J/mol)
- $k_B$ = Boltzmann constant ($8.617 \times 10^{-5}$ eV/K)
- $T$ = absolute temperature (K)
Key Implication: A temperature variation of just 10°C can change reaction rates by 20-30%.
1.2 Diffusion Fundamentals
Dopant diffusion follows Fick's Laws with temperature-dependent diffusivity:
$$
D = D_0 \cdot e^{-\frac{E_a}{k_B T}}
$$
Fick's First Law (steady-state diffusion):
$$
J = -D \frac{\partial C}{\partial x}
$$
Fick's Second Law (time-dependent diffusion):
$$
\frac{\partial C}{\partial t} = D \frac{\partial^2 C}{\partial x^2}
$$
Where:
- $J$ = diffusion flux (atoms/cm²·s)
- $D$ = diffusivity (cm²/s)
- $C$ = concentration (atoms/cm³)
- $D_0$ = pre-exponential diffusion coefficient
2. Key Thermal Processes in Semiconductor Manufacturing
2.1 Thermal Oxidation
Silicon dioxide growth follows the Deal-Grove Model:
$$
x_{ox}^2 + A \cdot x_{ox} = B(t + \tau)
$$
Where:
- $x_{ox}$ = oxide thickness
- $A$, $B$ = rate constants (temperature-dependent)
- $t$ = oxidation time
- $\tau$ = time offset for initial oxide
Oxidation Reactions:
- Dry oxidation: $\text{Si} + \text{O}_2 \rightarrow \text{SiO}_2$ (800–1200°C)
- Wet oxidation: $\text{Si} + 2\text{H}_2\text{O} \rightarrow \text{SiO}_2 + 2\text{H}_2$
Critical Parameters:
- Temperature uniformity requirement: $\pm 0.5°C$
- Typical temperature range: 800–1200°C
- Ramp rate affects interface quality and stress
2.2 Chemical Vapor Deposition (CVD)
Deposition Rate Temperature Dependence:
$$
R_{dep} = R_0 \cdot e^{-\frac{E_a}{k_B T}} \cdot P_{reactant}^n
$$
| CVD Type | Temperature Range | Pressure |
|----------|-------------------|----------|
| LPCVD | 400–900°C | 0.1–10 Torr |
| PECVD | 200–400°C | 0.1–10 Torr |
| APCVD | 300–500°C | 760 Torr |
| ALD | 150–400°C | 0.1–10 Torr |
Temperature affects:
- Deposition rate
- Film composition and stoichiometry
- Step coverage conformality
- Intrinsic film stress
- Grain structure and crystallinity
2.3 Rapid Thermal Processing (RTP)
Heat Balance Equation:
$$
\rho c_p V \frac{dT}{dt} = \alpha_{abs} P_{lamp} A - \varepsilon \sigma A (T^4 - T_{amb}^4) - h A (T - T_{amb})
$$
Where:
- $\rho$ = density (kg/m³)
- $c_p$ = specific heat capacity (J/kg·K)
- $V$ = wafer volume
- $\alpha_{abs}$ = optical absorptivity
- $P_{lamp}$ = lamp power density (W/m²)
- $\varepsilon$ = emissivity
- $\sigma$ = Stefan-Boltzmann constant ($5.67 \times 10^{-8}$ W/m²·K⁴)
- $h$ = convective heat transfer coefficient
RTP Specifications:
- Ramp rates: 50–400°C/s
- Peak temperatures: up to 1100°C
- Soak times: 0–60 seconds
- Spike anneal: ~1050°C, 0 second soak
2.4 Ion Implantation and Annealing
Implant Damage Annealing:
$$
f_{activated} = 1 - e^{-\left(\frac{t}{\tau}\right)^n}
$$
Where $\tau$ is the characteristic annealing time (temperature-dependent).
Annealing Methods:
| Method | Temperature | Time | Application |
|--------|-------------|------|-------------|
| Furnace Anneal | 800–1000°C | 30–60 min | Bulk damage repair |
| RTP Spike | 1000–1100°C | ~1 s | USJ activation |
| Flash Anneal | 1200–1350°C | 1–20 ms | Minimal diffusion |
| Laser Anneal | 1300–1414°C | 0.1–10 μs | Maximum activation |
3. Heat Transfer Mechanisms
3.1 Conduction
Fourier's Law:
$$
\vec{q} = -k
abla T
$$
3D Heat Equation:
$$
\rho c_p \frac{\partial T}{\partial t} = k
abla^2 T + \dot{Q}
$$
Or in Cartesian coordinates:
$$
\rho c_p \frac{\partial T}{\partial t} = k \left( \frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} + \frac{\partial^2 T}{\partial z^2} \right) + \dot{Q}
$$
Silicon Thermal Properties:
| Property | Value | Temperature Dependence |
|----------|-------|------------------------|
| Thermal conductivity | ~150 W/m·K @ 300K | $k \propto T^{-1.3}$ |
| Thermal diffusivity | ~0.9 cm²/s @ 300K | Decreases with T |
| Specific heat | ~700 J/kg·K @ 300K | Increases with T |
3.2 Radiation
Stefan-Boltzmann Law:
$$
q_{rad} = \varepsilon \sigma (T_s^4 - T_{surr}^4)
$$
Planck's Distribution:
$$
E_b(\lambda, T) = \frac{2\pi h c^2}{\lambda^5} \cdot \frac{1}{e^{\frac{hc}{\lambda k_B T}} - 1}
$$
Wien's Displacement Law:
$$
\lambda_{max} \cdot T = 2897.8 \text{ } \mu\text{m} \cdot \text{K}
$$
Or equivalently: $\lambda_{max} = \frac{2897.8}{T} \text{ } \mu\text{m}$ (where $T$ is in Kelvin)
Silicon Emissivity Considerations:
- Heavily doped Si: $\varepsilon \approx 0.7$
- Lightly doped Si: $\varepsilon \approx 0.3$ (semi-transparent in IR)
- With oxide film: interference effects modify $\varepsilon$
- Temperature dependent: $\varepsilon$ changes with $T$
3.3 Convection
Newton's Law of Cooling:
$$
q_{conv} = h(T_s - T_\infty)
$$
Nusselt Number Correlations:
For forced convection over a wafer:
$$
Nu = \frac{hL}{k_f} = C \cdot Re^m \cdot Pr^n
$$
Where:
- $Re = \frac{\rho v L}{\mu}$ (Reynolds number)
- $Pr = \frac{c_p \mu}{k_f}$ (Prandtl number)
4. Temperature Measurement
4.1 Pyrometry Fundamentals
Monochromatic Pyrometry:
$$
T = \frac{c_2}{\lambda \ln\left( \frac{\varepsilon c_1}{\lambda^5 L} + 1 \right)}
$$
Where:
- $c_1 = 3.742 \times 10^{-16}$ W·m²
- $c_2 = 1.439 \times 10^{-2}$ m·K
- $L$ = measured spectral radiance
- $\varepsilon$ = spectral emissivity
Two-Color (Ratio) Pyrometry:
$$
T = \frac{c_2 \left( \frac{1}{\lambda_1} - \frac{1}{\lambda_2} \right)}{\ln\left( \frac{L_1 \lambda_1^5}{L_2 \lambda_2^5} \cdot \frac{\varepsilon_2}{\varepsilon_1} \right)}
$$
Measurement Challenges:
- Unknown emissivity (varies with films, doping, temperature)
- Reflected radiation from chamber walls
- Transmission through silicon at certain wavelengths ($\lambda > 1.1$ μm)
- Pattern effects causing local emissivity variation
4.2 Contact Methods
- Thermocouples: $V = S_{AB} \cdot \Delta T$ (Seebeck coefficient)
- RTDs: $R(T) = R_0[1 + \alpha(T - T_0)]$
5. Thermal Stress Analysis
5.1 Thermal Stress Equations
Biaxial Thermal Stress in Thin Film:
$$
\sigma_{th} = \frac{E_f}{1 -
u_f} (\alpha_s - \alpha_f)(T - T_{dep})
$$
Where:
- $E_f$ = film Young's modulus
- $
u_f$ = film Poisson's ratio
- $\alpha_s$ = substrate CTE
- $\alpha_f$ = film CTE
- $T_{dep}$ = deposition temperature
Wafer Bow (Stoney's Equation):
$$
\sigma_f = \frac{E_s t_s^2}{6(1-
u_s) t_f} \cdot \frac{1}{R}
$$
Where:
- $t_s$ = substrate thickness
- $t_f$ = film thickness
- $R$ = radius of curvature
5.2 Slip Dislocation Criterion
Slip occurs when resolved shear stress exceeds critical value:
$$
\tau_{resolved} = \sigma \cdot \cos\phi \cdot \cos\lambda > \tau_{CRSS}(T)
$$
Critical Temperature: Slip typically begins above ~1050°C in silicon.
Temperature Gradient Stress:
$$
\sigma_{gradient} \approx \frac{E \alpha \Delta T}{1 -
u}
$$
6. Nanoscale Thermal Transport
6.1 Phonon Transport
When feature sizes approach phonon mean free path ($\Lambda_{mfp} \approx 100-300$ nm in Si at 300K):
Ballistic Transport Regime:
$$
q = \frac{1}{4} C v_{ph} \Delta T \quad \text{(when } L < \Lambda_{mfp}\text{)}
$$
Modified Thermal Conductivity:
$$
k_{eff} = k_{bulk} \cdot \frac{1}{1 + \frac{\Lambda_{mfp}}{L}}
$$
6.2 Interface Thermal Resistance (Kapitza Resistance)
$$
R_{th,interface} = \frac{\Delta T}{q} = R_{Kapitza}
$$
Acoustic Mismatch Model:
$$
R_{Kapitza} \propto \frac{(\rho_1 v_1 - \rho_2 v_2)^2}{(\rho_1 v_1 + \rho_2 v_2)^2}
$$
Where $\rho v$ is the acoustic impedance.
7. Equipment and Process Parameters
7.1 Batch Furnace Specifications
- Temperature uniformity: $\pm 0.5°C$ across wafer zone
- Ramp rates: 1–10°C/min
- Maximum temperature: 1200°C
- Batch size: 50–150 wafers
7.2 RTP System Parameters
- Lamp types:
- Tungsten-halogen: $\lambda_{peak} \approx 1$ μm
- Arc lamps: broadband emission
- Ramp rates: 50–400°C/s
- Temperature uniformity target: $\pm 2°C$
7.3 Laser Annealing Parameters
| Parameter | Excimer Laser | CW Laser |
|-----------|---------------|----------|
| Wavelength | 308 nm (XeCl) | 532 nm, 808 nm |
| Pulse duration | 10–100 ns | Continuous |
| Melt depth | 10–100 nm | 1–10 μm |
| Peak temperature | >1414°C (melt) | 1200–1414°C |
8. Process Integration Considerations
8.1 Thermal Budget
Cumulative Thermal Budget:
$$
D_t = \sum_i D_0 \cdot e^{-\frac{E_a}{k_B T_i}} \cdot t_i
$$
Where $D_t$ is the total diffusion length squared.
Effective $D \cdot t$:
$$
(Dt)_{eff} = \int_0^{t_{process}} D(T(t')) dt'
$$
8.2 Junction Depth Estimation
For constant-source diffusion:
$$
x_j = 2\sqrt{Dt} \cdot \text{erfc}^{-1}\left(\frac{C_B}{C_s}\right)
$$
Where:
- $x_j$ = junction depth
- $C_B$ = background concentration
- $C_s$ = surface concentration
9. Key Equations
| Process | Key Equation | Critical Parameters |
|---------|--------------|---------------------|
| Reaction Rate | $k = A e^{-E_a/k_B T}$ | $E_a$, $T$ |
| Diffusion | $D = D_0 e^{-E_a/k_B T}$ | $D_0$, $E_a$ |
| Oxidation | $x^2 + Ax = B(t+\tau)$ | $A$, $B$ (T-dependent) |
| Radiation | $q = \varepsilon \sigma T^4$ | $\varepsilon$, $T$ |
| Thermal Stress | $\sigma = \frac{E}{1-
u}\Delta\alpha\Delta T$ | CTE mismatch |
| Heat Conduction | $q = -k
abla T$ | $k(T)$ |