epi modeling
**Semiconductor Manufacturing Process: Epitaxy (Epi) Modeling**
**1. Introduction to Epitaxy**
Epitaxy is the controlled growth of a crystalline thin film on a crystalline substrate, where the deposited layer inherits the crystallographic orientation of the substrate.
**1.1 Types of Epitaxy**
- **Homoepitaxy**
- Same material deposited on substrate
- Example: Silicon (Si) on Silicon (Si)
- Maintains perfect lattice matching
- Used for creating high-purity device layers
- **Heteroepitaxy**
- Different material deposited on substrate
- Examples:
- Gallium Arsenide (GaAs) on Silicon (Si)
- Silicon Germanium (SiGe) on Silicon (Si)
- Gallium Nitride (GaN) on Sapphire ($\text{Al}_2\text{O}_3$)
- Introduces lattice mismatch and strain
- Enables bandgap engineering
**2. Epitaxy Methods**
**2.1 Chemical Vapor Deposition (CVD) / Vapor Phase Epitaxy (VPE)**
- **Characteristics:**
- Most common method for silicon epitaxy
- Operates at atmospheric or reduced pressure
- Temperature range: $900°\text{C} - 1200°\text{C}$
- **Common Precursors:**
- Silane: $\text{SiH}_4$
- Dichlorosilane: $\text{SiH}_2\text{Cl}_2$ (DCS)
- Trichlorosilane: $\text{SiHCl}_3$ (TCS)
- Silicon tetrachloride: $\text{SiCl}_4$
- **Key Reactions:**
$$\text{SiH}_4 \xrightarrow{\Delta} \text{Si}_{(s)} + 2\text{H}_2$$
$$\text{SiH}_2\text{Cl}_2 \xrightarrow{\Delta} \text{Si}_{(s)} + 2\text{HCl}$$
**2.2 Molecular Beam Epitaxy (MBE)**
- **Characteristics:**
- Ultra-high vacuum environment ($< 10^{-10}$ Torr)
- Extremely precise thickness control (monolayer accuracy)
- Lower growth temperatures than CVD
- Slower growth rates: $\sim 1 \, \mu\text{m/hour}$
- **Applications:**
- III-V compound semiconductors
- Quantum well structures
- Superlattices
- Research and development
**2.3 Metal-Organic CVD (MOCVD)**
- **Characteristics:**
- Standard for compound semiconductors
- Uses metal-organic precursors
- Higher throughput than MBE
- **Common Precursors:**
- Trimethylgallium: $\text{Ga(CH}_3\text{)}_3$ (TMGa)
- Trimethylaluminum: $\text{Al(CH}_3\text{)}_3$ (TMAl)
- Ammonia: $\text{NH}_3$
**2.4 Atomic Layer Epitaxy (ALE)**
- **Characteristics:**
- Self-limiting surface reactions
- Digital control of film thickness
- Excellent conformality
- Growth rate: $\sim 1$ Å per cycle
**3. Physics of Epi Modeling**
**3.1 Gas-Phase Transport**
The transport of precursor gases to the substrate surface involves multiple phenomena:
- **Governing Equations:**
- **Continuity Equation:**
$$\frac{\partial \rho}{\partial t} +
abla \cdot (\rho \mathbf{v}) = 0$$
- **Navier-Stokes Equation:**
$$\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}$$
- **Species Transport Equation:**
$$\frac{\partial C_i}{\partial t} + \mathbf{v} \cdot
abla C_i = D_i
abla^2 C_i + R_i$$
Where:
- $\rho$ = fluid density
- $\mathbf{v}$ = velocity vector
- $p$ = pressure
- $\mu$ = dynamic viscosity
- $C_i$ = concentration of species $i$
- $D_i$ = diffusion coefficient of species $i$
- $R_i$ = reaction rate term
- **Boundary Layer:**
- Stagnant gas layer above substrate
- Thickness $\delta$ depends on flow conditions:
$$\delta \propto \sqrt{\frac{
u x}{u_\infty}}$$
Where:
- $
u$ = kinematic viscosity
- $x$ = distance from leading edge
- $u_\infty$ = free stream velocity
**3.2 Surface Kinetics**
- **Adsorption Process:**
- Physisorption (weak van der Waals forces)
- Chemisorption (chemical bonding)
- **Langmuir Adsorption Isotherm:**
$$\theta = \frac{K \cdot P}{1 + K \cdot P}$$
Where:
- $\theta$ = fractional surface coverage
- $K$ = equilibrium constant
- $P$ = partial pressure
- **Surface Diffusion:**
$$D_s = D_0 \exp\left(-\frac{E_d}{k_B T}\right)$$
Where:
- $D_s$ = surface diffusion coefficient
- $D_0$ = pre-exponential factor
- $E_d$ = diffusion activation energy
- $k_B$ = Boltzmann constant ($1.38 \times 10^{-23}$ J/K)
- $T$ = absolute temperature
**3.3 Crystal Growth Mechanisms**
- **Step-Flow Growth (BCF Theory):**
- Atoms attach at step edges
- Steps advance across terraces
- Dominant at high temperatures
- **2D Nucleation:**
- New layers nucleate on terraces
- Occurs when step density is low
- Creates rougher surfaces
- **Terrace-Ledge-Kink (TLK) Model:**
- Terrace: flat regions between steps
- Ledge: step edges
- Kink: incorporation sites at step edges
**4. Mathematical Framework**
**4.1 Growth Rate Models**
**4.1.1 Reaction-Limited Regime**
At lower temperatures, surface reaction kinetics dominate:
$$G = k_s \cdot C_s$$
Where the rate constant follows Arrhenius behavior:
$$k_s = k_0 \exp\left(-\frac{E_a}{k_B T}\right)$$
**Parameters:**
- $G$ = growth rate (nm/min or μm/hr)
- $k_s$ = surface reaction rate constant
- $C_s$ = surface concentration
- $k_0$ = pre-exponential factor
- $E_a$ = activation energy
**4.1.2 Mass-Transport Limited Regime**
At higher temperatures, diffusion through the boundary layer limits growth:
$$G = \frac{h_g}{N_s} \cdot (C_g - C_s)$$
Where:
$$h_g = \frac{D}{\delta}$$
**Parameters:**
- $h_g$ = mass transfer coefficient
- $N_s$ = atomic density of solid ($\sim 5 \times 10^{22}$ atoms/cm³ for Si)
- $C_g$ = gas phase concentration
- $D$ = gas phase diffusivity
- $\delta$ = boundary layer thickness
**4.1.3 Combined Model (Grove Model)**
For the general case combining both regimes:
$$G = \frac{h_g \cdot k_s}{N_s (h_g + k_s)} \cdot C_g$$
Or equivalently:
$$\frac{1}{G} = \frac{N_s}{k_s \cdot C_g} + \frac{N_s}{h_g \cdot C_g}$$
**4.2 Strain in Heteroepitaxy**
**4.2.1 Lattice Mismatch**
$$f = \frac{a_s - a_f}{a_f}$$
Where:
- $f$ = lattice mismatch (dimensionless)
- $a_s$ = substrate lattice constant
- $a_f$ = film lattice constant (relaxed)
**Example Values:**
| System | $a_f$ (Å) | $a_s$ (Å) | Mismatch $f$ |
|--------|-----------|-----------|--------------|
| Si on Si | 5.431 | 5.431 | 0% |
| Ge on Si | 5.658 | 5.431 | -4.2% |
| GaAs on Si | 5.653 | 5.431 | -4.1% |
| InAs on GaAs | 6.058 | 5.653 | -7.2% |
**4.2.2 In-Plane Strain**
For a coherently strained film:
$$\epsilon_{\parallel} = \frac{a_s - a_f}{a_f} = f$$
The out-of-plane strain (for cubic materials):
$$\epsilon_{\perp} = -\frac{2
u}{1-
u} \epsilon_{\parallel}$$
Where $
u$ = Poisson's ratio
**4.2.3 Critical Thickness (Matthews-Blakeslee)**
The critical thickness above which misfit dislocations form:
$$h_c = \frac{b}{8\pi f (1+
u)} \left[ \ln\left(\frac{h_c}{b}\right) + 1 \right]$$
Where:
- $h_c$ = critical thickness
- $b$ = Burgers vector magnitude ($\approx \frac{a}{\sqrt{2}}$ for 60° dislocations)
- $f$ = lattice mismatch
- $
u$ = Poisson's ratio
**Approximate Solution:**
For small mismatch:
$$h_c \approx \frac{b}{8\pi |f|}$$
**4.3 Dopant Incorporation**
**4.3.1 Segregation Model**
$$C_{film} = \frac{C_{gas}}{1 + k_{seg} \cdot (G/G_0)}$$
Where:
- $C_{film}$ = dopant concentration in film
- $C_{gas}$ = dopant concentration in gas phase
- $k_{seg}$ = segregation coefficient
- $G$ = growth rate
- $G_0$ = reference growth rate
**4.3.2 Dopant Profile with Segregation**
The surface concentration evolves as:
$$C_s(t) = C_s^{eq} + (C_s(0) - C_s^{eq}) \exp\left(-\frac{G \cdot t}{\lambda}\right)$$
Where:
- $\lambda$ = segregation length
- $C_s^{eq}$ = equilibrium surface concentration
**5. Modeling Approaches**
**5.1 Continuum Models**
- **Scope:**
- Reactor-scale simulations
- Temperature and flow field prediction
- Species concentration profiles
- **Methods:**
- Computational Fluid Dynamics (CFD)
- Finite Element Method (FEM)
- Finite Volume Method (FVM)
- **Governing Physics:**
- Coupled heat, mass, and momentum transfer
- Homogeneous and heterogeneous reactions
- Radiation heat transfer
**5.2 Feature-Scale Models**
- **Applications:**
- Selective epitaxial growth (SEG)
- Trench filling
- Facet evolution
- **Key Phenomena:**
- Local loading effects:
$$G_{local} = G_0 \cdot \left(1 - \alpha \cdot \frac{A_{exposed}}{A_{total}}\right)$$
- Orientation-dependent growth rates:
$$\frac{G_{(110)}}{G_{(100)}} \approx 1.5 - 2.0$$
- **Methods:**
- Level set methods
- String methods
- Cellular automata
**5.3 Atomistic Models**
**5.3.1 Kinetic Monte Carlo (KMC)**
- **Process Events:**
- Adsorption: rate $\propto P \cdot \exp(-E_{ads}/k_BT)$
- Surface diffusion: rate $\propto \exp(-E_{diff}/k_BT)$
- Desorption: rate $\propto \exp(-E_{des}/k_BT)$
- Incorporation: rate $\propto \exp(-E_{inc}/k_BT)$
- **Master Equation:**
$$\frac{dP_i}{dt} = \sum_j \left( W_{ji} P_j - W_{ij} P_i \right)$$
Where:
- $P_i$ = probability of state $i$
- $W_{ij}$ = transition rate from state $i$ to $j$
**5.3.2 Molecular Dynamics (MD)**
- **Newton's Equations:**
$$m_i \frac{d^2 \mathbf{r}_i}{dt^2} = -
abla_i U(\mathbf{r}_1, \mathbf{r}_2, ..., \mathbf{r}_N)$$
- **Interatomic Potentials:**
- Tersoff potential (Si, C, Ge)
- Stillinger-Weber potential (Si)
- MEAM (metals and alloys)
**5.3.3 Ab Initio / DFT**
- **Kohn-Sham Equations:**
$$\left[ -\frac{\hbar^2}{2m}
abla^2 + V_{eff}(\mathbf{r}) \right] \psi_i(\mathbf{r}) = \epsilon_i \psi_i(\mathbf{r})$$
- **Applications:**
- Surface energies
- Reaction barriers
- Adsorption energies
- Electronic structure
**6. Specific Modeling Challenges**
**6.1 SiGe Epitaxy**
- **Composition Control:**
$$x_{Ge} = \frac{R_{Ge}}{R_{Si} + R_{Ge}}$$
Where $R_{Si}$ and $R_{Ge}$ are partial growth rates
- **Strain Engineering:**
- Compressive strain in SiGe on Si
- Enhances hole mobility
- Critical thickness depends on Ge content:
$$h_c(x) \approx \frac{0.5}{0.042 \cdot x} \text{ nm}$$
**6.2 Selective Epitaxy**
- **Growth Selectivity:**
- Deposition only on exposed silicon
- HCl addition for selectivity enhancement
- **Selectivity Condition:**
$$\frac{\text{Growth on Si}}{\text{Growth on SiO}_2} > 100:1$$
- **Loading Effects:**
- Pattern-dependent growth rate
- Faceting at mask edges
**6.3 III-V on Silicon**
- **Major Challenges:**
- Large lattice mismatch (4-8%)
- Thermal expansion mismatch
- Anti-phase domain boundaries (APDs)
- High threading dislocation density
- **Mitigation Strategies:**
- Aspect ratio trapping (ART)
- Graded buffer layers
- Selective area growth
- Dislocation filtering
**7. Applications and Tools**
**7.1 Industrial Applications**
| Application | Material System | Key Parameters |
|-------------|-----------------|----------------|
| FinFET/GAA Source/Drain | Embedded SiGe, SiC | Strain, selectivity |
| SiGe HBT | SiGe:C | Profile abruptness |
| Power MOSFETs | SiC epitaxy | Defect density |
| LEDs/Lasers | GaN, InGaN | Composition uniformity |
| RF Devices | GaN on SiC | Buffer quality |
**7.2 Simulation Software**
- **Reactor-Scale CFD:**
- ANSYS Fluent
- COMSOL Multiphysics
- OpenFOAM
- **TCAD Process Simulation:**
- Synopsys Sentaurus Process
- Silvaco Victory Process
- Lumerical (for optoelectronics)
- **Atomistic Simulation:**
- LAMMPS (MD)
- VASP, Quantum ESPRESSO (DFT)
- Custom KMC codes
**7.3 Key Metrics for Process Development**
- **Uniformity:**
$$\text{Uniformity} = \frac{t_{max} - t_{min}}{2 \cdot t_{avg}} \times 100\%$$
- **Defect Density:**
- Threading dislocations: target $< 10^6$ cm$^{-2}$
- Stacking faults: target $< 10^3$ cm$^{-2}$
- **Profile Abruptness:**
- Dopant transition width $< 3$ nm/decade
**8. Emerging Directions**
**8.1 Machine Learning Integration**
- **Applications:**
- Surrogate models for process optimization
- Real-time virtual metrology
- Defect classification
- Recipe optimization
- **Model Types:**
- Neural networks for growth rate prediction
- Gaussian process regression for uncertainty quantification
- Reinforcement learning for process control
**8.2 Multi-Scale Modeling**
- **Hierarchical Approach:**
```
Ab Initio (DFT)
↓ Reaction rates, energies
Kinetic Monte Carlo
↓ Surface kinetics, morphology
Feature-Scale Models
↓ Local growth behavior
Reactor-Scale CFD
↓ Process conditions
Device Simulation
```
**8.3 Digital Twins**
- **Components:**
- Real-time sensor data integration
- Physics-based + ML hybrid models
- Predictive maintenance
- Closed-loop process control
**8.4 New Material Systems**
- **2D Materials:**
- Graphene via CVD
- Transition metal dichalcogenides (TMDs)
- Van der Waals epitaxy
- **Ultra-Wide Bandgap:**
- $\beta$-Ga$_2$O$_3$ ($E_g \approx 4.8$ eV)
- Diamond ($E_g \approx 5.5$ eV)
- AlN ($E_g \approx 6.2$ eV)
**Common Constants and Conversions**
| Constant | Symbol | Value |
|----------|--------|-------|
| Boltzmann constant | $k_B$ | $1.381 \times 10^{-23}$ J/K |
| Planck constant | $h$ | $6.626 \times 10^{-34}$ J·s |
| Avogadro number | $N_A$ | $6.022 \times 10^{23}$ mol$^{-1}$ |
| Si atomic density | $N_{Si}$ | $5.0 \times 10^{22}$ atoms/cm³ |
| Si lattice constant | $a_{Si}$ | 5.431 Å |