Mathematical Modeling of Epitaxy in…

Home› Knowledge Base› Mathematical Modeling of Epitaxy in Semiconductor Front-End Processing (FEP)

Mathematical Modeling of Epitaxy in Semiconductor Front-End Processing (FEP)

<svg viewBox="0 0 760 470" xmlns="http://www.w3.org/2000/svg" font-family="-apple-system,BlinkMacSystemFont,Segoe UI,Roboto,sans-serif">
  <rect x="0" y="0" width="760" height="470" fill="#0d1117"/>
  <text x="380" y="32" fill="#e6edf3" font-size="20" font-weight="700" text-anchor="middle">Epitaxial Growth Mechanics — Lattice Matching &amp; Strain Engineering</text>
  <text x="380" y="52" fill="#8b98a5" font-size="12" text-anchor="middle">Pseudomorphic Lattice Distortion, Critical Thickness Thresholds &amp; Selective Growth (SEG)</text>

  <!-- Left Panel: Atomic Lattice Matching & Pseudomorphic Strain -->
  <g transform="translate(25, 70)">
    <rect width="345" height="330" fill="#161b22" stroke="#30363d" stroke-width="1.5" rx="8"/>
    <text x="172.5" y="24" fill="#60a5fa" font-size="13" font-weight="700" text-anchor="middle">1. Lattice Matching &amp; Strain States</text>
    
    <!-- Sub-panel A: Homoepitaxy -->
    <g transform="translate(15, 38)">
      <rect width="315" height="125" fill="#0d1117" stroke="#3b82f6" stroke-width="1" rx="5"/>
      <text x="157.5" y="18" fill="#60a5fa" font-size="10.5" font-weight="700" text-anchor="middle">A. Homoepitaxy (Si on Si — Matched Lattice)</text>
      
      <!-- Atomic Grid -->
      <g fill="#60a5fa">
        <circle cx="50" cy="40" r="4"/><circle cx="90" cy="40" r="4"/><circle cx="130" cy="40" r="4"/><circle cx="170" cy="40" r="4"/><circle cx="210" cy="40" r="4"/><circle cx="250" cy="40" r="4"/>
        <circle cx="50" cy="65" r="4"/><circle cx="90" cy="65" r="4"/><circle cx="130" cy="65" r="4"/><circle cx="170" cy="65" r="4"/><circle cx="210" cy="65" r="4"/><circle cx="250" cy="65" r="4"/>
      </g>
      <g stroke="#3b82f6" stroke-width="1" stroke-opacity="0.5">
        <line x1="50" y1="40" x2="250" y2="40"/><line x1="50" y1="65" x2="250" y2="65"/>
        <line x1="50" y1="40" x2="50" y2="105"/><line x1="90" y1="40" x2="90" y2="105"/><line x1="130" y1="40" x2="130" y2="105"/><line x1="170" y1="40" x2="170" y2="105"/><line x1="210" y1="40" x2="210" y2="105"/><line x1="250" y1="40" x2="250" y2="105"/>
      </g>
      <!-- Substrate Grid -->
      <g fill="#475569">
        <circle cx="50" cy="90" r="4"/><circle cx="90" cy="90" r="4"/><circle cx="130" cy="90" r="4"/><circle cx="170" cy="90" r="4"/><circle cx="210" cy="90" r="4"/><circle cx="250" cy="90" r="4"/>
      </g>
      <text x="157.5" y="118" fill="#8b98a5" font-size="8.5" text-anchor="middle">a_film = a_sub | Zero Strain (ε = 0) | Perfect Crystal Continuity</text>
    </g>

    <!-- Sub-panel B: Heteroepitaxy Pseudomorphic Strain -->
    <g transform="translate(15, 175)">
      <rect width="315" height="140" fill="#0d1117" stroke="#f59e0b" stroke-width="1" rx="5"/>
      <text x="157.5" y="18" fill="#fbbf24" font-size="10.5" font-weight="700" text-anchor="middle">B. Heteroepitaxy (Compressive SiGe on Si)</text>
      
      <!-- Strained Film Atoms (Vertical elongation) -->
      <g fill="#fbbf24">
        <circle cx="50" cy="35" r="4.5"/><circle cx="90" cy="35" r="4.5"/><circle cx="130" cy="35" r="4.5"/><circle cx="170" cy="35" r="4.5"/><circle cx="210" cy="35" r="4.5"/><circle cx="250" cy="35" r="4.5"/>
        <circle cx="50" cy="65" r="4.5"/><circle cx="90" cy="65" r="4.5"/><circle cx="130" cy="65" r="4.5"/><circle cx="170" cy="65" r="4.5"/><circle cx="210" cy="65" r="4.5"/><circle cx="250" cy="65" r="4.5"/>
      </g>
      <!-- Substrate Atoms -->
      <g fill="#475569">
        <circle cx="50" cy="95" r="4"/><circle cx="90" cy="95" r="4"/><circle cx="130" cy="95" r="4"/><circle cx="170" cy="95" r="4"/><circle cx="210" cy="95" r="4"/><circle cx="250" cy="95" r="4"/>
      </g>
      <!-- Compressive arrows -->
      <line x1="30" y1="50" x2="42" y2="50" stroke="#ef4444" stroke-width="1.5"/>
      <polygon points="45,50 40,47 40,53" fill="#ef4444"/>
      <line x1="270" y1="50" x2="258" y2="50" stroke="#ef4444" stroke-width="1.5"/>
      <polygon points="255,50 260,47 260,53" fill="#ef4444"/>

      <text x="157.5" y="114" fill="#fca5a5" font-size="8.5" font-weight="700" text-anchor="middle">Compressive Strain: a_parallel = a_sub, a_perp &gt; a_sub</text>
      <text x="157.5" y="128" fill="#8b98a5" font-size="8" text-anchor="middle">Coherent Pseudomorphic Growth below Critical Thickness h_c</text>
    </g>
  </g>

  <!-- Right Panel: Selective Epitaxial Growth (SEG) in Transistors -->
  <g transform="translate(390, 70)">
    <rect width="345" height="330" fill="#161b22" stroke="#30363d" stroke-width="1.5" rx="8"/>
    <text x="172.5" y="24" fill="#34d399" font-size="13" font-weight="700" text-anchor="middle">2. Selective Epitaxy (SEG) in FinFET/GAA</text>

    <!-- Transistor Structure with e-SiGe S/D -->
    <g transform="translate(15, 38)">
      <rect width="315" height="140" fill="#0d1117" stroke="#10b981" stroke-width="1.2" rx="6"/>
      <text x="157.5" y="18" fill="#3fb950" font-size="11" font-weight="700" text-anchor="middle">Embedded SiGe Source/Drain (Uniaxial Strain)</text>

      <!-- Gate Oxide & Gate Electrodes -->
      <rect x="120" y="32" width="75" height="25" fill="#1e1b4b" stroke="#818cf8" stroke-width="1" rx="2"/>
      <text x="157.5" y="48" fill="#c7d2fe" font-size="9" font-weight="700" text-anchor="middle">Gate (HKMG)</text>
      
      <!-- Spacer Walls -->
      <rect x="105" y="32" width="15" height="25" fill="#334155"/>
      <rect x="195" y="32" width="15" height="25" fill="#334155"/>

      <!-- Silicon Substrate & Channel -->
      <rect x="25" y="57" width="265" height="70" fill="#1e293b" stroke="#64748b" rx="2"/>
      <text x="157.5" y="78" fill="#6ee7b7" font-size="9" font-weight="700" text-anchor="middle">Si Channel (Compressively Strained)</text>
      
      <!-- Embedded SiGe S/D Facets -->
      <polygon points="35,57 95,57 105,95 25,95" fill="#d97706" stroke="#fbbf24"/>
      <text x="65" y="80" fill="#ffffff" font-size="8.5" font-weight="700" text-anchor="middle">e-SiGe</text>
      
      <polygon points="210,57 280,57 290,95 200,95" fill="#d97706" stroke="#fbbf24"/>
      <text x="245" y="80" fill="#ffffff" font-size="8.5" font-weight="700" text-anchor="middle">e-SiGe</text>

      <!-- Strain Force Arrows towards Channel -->
      <line x1="95" y1="72" x2="120" y2="72" stroke="#ef4444" stroke-width="2"/>
      <polygon points="123,72 117,69 117,75" fill="#ef4444"/>

      <line x1="210" y1="72" x2="185" y2="72" stroke="#ef4444" stroke-width="2"/>
      <polygon points="182,72 188,69 188,75" fill="#ef4444"/>

      <text x="157.5" y="118" fill="#34d399" font-size="8" font-weight="700" text-anchor="middle">Boosts pMOS Hole Mobility &gt; 50% via Longitudinal Strain</text>
    </g>

    <!-- Sub-panel B: Selectivity & Faceting Mechanics -->
    <g transform="translate(15, 190)">
      <rect width="315" height="125" fill="#0d1117" stroke="#30363d" rx="6"/>
      <text x="157.5" y="18" fill="#e6edf3" font-size="10.5" font-weight="700" text-anchor="middle">Growth Kinetics &amp; HCl Chemical Selectivity</text>
      
      <text x="25" y="42" fill="#8b98a5" font-size="8.5">• <tspan fill="#38bdf8" font-weight="700">Gas Phase Precursors</tspan>: DCS (SiH₂Cl₂) + GeH₄ + HCl + B₂H₆</text>
      <text x="25" y="60" fill="#8b98a5" font-size="8.5">• <tspan fill="#fbbf24" font-weight="700">HCl Role</tspan>: Etches amorphous nuclei on SiO₂/SiN mask</text>
      <text x="25" y="78" fill="#8b98a5" font-size="8.5">• <tspan fill="#34d399" font-weight="700">Faceting</tspan>: {111} planes have slowest growth rate</text>
      
      <rect x="25" y="90" width="265" height="22" fill="#0f172a" stroke="#60a5fa" stroke-width="0.8" rx="3"/>
      <text x="157.5" y="105" fill="#60a5fa" font-size="8.5" font-weight="700" text-anchor="middle">Atomic-Level Doping: In-situ Boron / Phosphorus Insertion</text>
    </g>
  </g>

  <!-- Bottom Key Insight Bar -->
  <rect x="25" y="415" width="710" height="22" rx="3" fill="#0b1220" stroke="#233043" stroke-width="0.8"/>
  <text x="380" y="430" fill="#fbbf24" font-size="9.5" font-weight="700" text-anchor="middle">Misfit Strain ε = (a_film - a_sub)/a_sub | Pseudomorphic growth below Critical Thickness h_c drives Channel Mobility</text>

  <!-- Footer -->
  <text x="380" y="458" fill="#6b7684" font-size="11" text-anchor="middle">Foundational FEOL process for GAA nanosheets, SiGe pMOS channels &amp; InP/GaAs RF photonics</text>
</svg>

1. Overview

Epitaxy is a critical Front-End Process (FEP) step where crystalline films are grown on crystalline substrates with precise control of:

Thickness
Composition
Doping concentration
Defect density

Mathematical modeling enables:

Process optimization
Defect prediction
Virtual fabrication
Equipment design

1.1 Types of Epitaxy

Homoepitaxy: Same material as substrate (e.g., Si on Si)
Heteroepitaxy: Different material from substrate (e.g., GaAs on Si, SiGe on Si)

1.2 Epitaxy Methods

Vapor Phase Epitaxy (VPE) / Chemical Vapor Deposition (CVD)

Atmospheric Pressure CVD (APCVD)

Low Pressure CVD (LPCVD)

Metal-Organic CVD (MOCVD)

Molecular Beam Epitaxy (MBE)
Liquid Phase Epitaxy (LPE)
Solid Phase Epitaxy (SPE)

2. Fundamental Thermodynamic Framework

2.1 Driving Force for Growth

The supersaturation provides the thermodynamic driving force:

\Delta \mu = k_B T \ln\left(\frac{P}{P_{eq}}\right)

Where:

$\Delta \mu$ = chemical potential difference (driving force)
$k_B$ = Boltzmann's constant ($1.38 \times 10^{-23}$ J/K)
$T$ = absolute temperature (K)
$P$ = actual partial pressure of precursor
$P_{eq}$ = equilibrium vapor pressure

2.2 Free Energy of Mixing (Multi-component Systems)

For systems like SiGe alloys:

\Delta G_{mix} = RT\left(x \ln x + (1-x) \ln(1-x)\right) + \Omega x(1-x)

Where:

$R$ = universal gas constant (8.314 J/mol$\cdot$K)
$x$ = mole fraction of component
$\Omega$ = interaction parameter (regular solution model)

2.3 Gibbs Free Energy of Formation

\Delta G = \Delta H - T\Delta S

For spontaneous growth: $\Delta G < 0$

3. Growth Rate Kinetics

3.1 The Two-Regime Model

Epitaxial growth rate is governed by two competing mechanisms:

Overall growth rate equation:

G = \frac{k_s \cdot h_g \cdot C_g}{k_s + h_g}

Where:

$G$ = growth rate (nm/min or $\mu$m/min)
$k_s$ = surface reaction rate constant
$h_g$ = gas-phase mass transfer coefficient
$C_g$ = gas-phase reactant concentration

3.2 Temperature Dependence

The surface reaction rate follows Arrhenius behavior:

k_s = A \exp\left(-\frac{E_a}{k_B T}\right)

Where:

$A$ = pre-exponential factor (frequency factor)
$E_a$ = activation energy (eV or J/mol)

3.3 Growth Rate Regimes

Temperature Regime	Limiting Factor	Growth Rate Expression	Temperature Dependence
Low T	Surface reaction	$G \approx k_s \cdot C_g$	Strong (exponential)
High T	Mass transport	$G \approx h_g \cdot C_g$	Weak (~$T^{1.5-2}$)

3.4 Boundary Layer Analysis

For horizontal CVD reactors, the boundary layer thickness evolves as:

\delta(x) = \sqrt{\frac{ u \cdot x}{v_{\infty}}}

Where:

$\delta(x)$ = boundary layer thickness at position $x$
$

u$ = kinematic viscosity (m²/s)

$x$ = distance from gas inlet (m)
$v_{\infty}$ = free stream gas velocity (m/s)

The mass transfer coefficient:

h_g = \frac{D_{gas}}{\delta}

Where $D_{gas}$ is the gas-phase diffusion coefficient.

4. Surface Kinetics: BCF Theory

The Burton-Cabrera-Frank (BCF) model describes atomic-scale growth mechanisms.

4.1 Surface Diffusion Equation

D_s \nabla^2 n_s - \frac{n_s - n_{eq}}{\tau_s} + J_{ads} = 0

Where:

$n_s$ = adatom surface density (atoms/cm²)
$D_s$ = surface diffusion coefficient (cm²/s)
$n_{eq}$ = equilibrium adatom density
$\tau_s$ = mean adatom lifetime before desorption (s)
$J_{ads}$ = adsorption flux (atoms/cm²$\cdot$s)

4.2 Characteristic Diffusion Length

\lambda_s = \sqrt{D_s \tau_s}

This parameter determines the growth mode:

Step-flow growth: $\lambda_s > L$ (terrace width)
2D nucleation growth: $\lambda_s < L$

4.3 Surface Diffusion Coefficient

D_s = D_0 \exp\left(-\frac{E_m}{k_B T}\right)

Where:

$D_0$ = pre-exponential factor (~$10^{-3}$ cm²/s)
$E_m$ = migration energy barrier (eV)

4.4 Step Velocity

v_{step} = \frac{2 D_s (n_s - n_{eq})}{\lambda_s} \tanh\left(\frac{L}{2\lambda_s}\right)

Where $L$ is the inter-step spacing (terrace width).

4.5 Growth Rate from Step Flow

G = \frac{v_{step} \cdot h_{step}}{L}

Where $h_{step}$ is the step height (monolayer thickness).

5. Heteroepitaxy and Strain Modeling

5.1 Lattice Mismatch

f = \frac{a_{film} - a_{substrate}}{a_{substrate}}

Where:

$f$ = lattice mismatch (dimensionless, often expressed as %)
$a_{film}$ = lattice constant of film material
$a_{substrate}$ = lattice constant of substrate

Example values:

System	Lattice Mismatch
Si₀.₇Ge₀.₃ on Si	~1.2%
Ge on Si	~4.2%
GaAs on Si	~4.0%
InAs on GaAs	~7.2%
GaN on Sapphire	~16%

5.2 Strain Components

For biaxial strain in (001) films:

\varepsilon_{xx} = \varepsilon_{yy} = \varepsilon_{\parallel} = \frac{a_s - a_f}{a_f} \approx -f

\varepsilon_{zz} = \varepsilon_{\perp} = -\frac{2C_{12}}{C_{11}} \varepsilon_{\parallel}

Where $C_{11}$ and $C_{12}$ are elastic constants.

5.3 Elastic Energy

For a coherently strained film:

E_{elastic} = \frac{2G(1+ u)}{1- u} f^2 h = M f^2 h

Where:

$G$ = shear modulus (Pa)
$

u$ = Poisson's ratio

$h$ = film thickness
$M$ = biaxial modulus = $\frac{2G(1+

u)}{1- u}$

5.4 Critical Thickness (Matthews-Blakeslee)

h_c = \frac{b}{8\pi f(1+ u)} \left[\ln\left(\frac{h_c}{b}\right) + 1\right]

Where:

$h_c$ = critical thickness for dislocation formation
$b$ = Burgers vector magnitude
$f$ = lattice mismatch
$

u$ = Poisson's ratio

5.5 People-Bean Approximation (for SiGe)

Empirical formula:

h_c \approx \frac{0.55}{f^2} \text{ (nm, with } f \text{ as a decimal)}

Or equivalently:

h_c \approx \frac{5500}{x^2} \text{ (nm, for Si}_{1-x}\text{Ge}_x\text{)}

5.6 Threading Dislocation Density

Above critical thickness, dislocation density evolves:

\rho_{TD}(h) = \rho_0 \exp\left(-\frac{h}{h_0}\right) + \rho_{\infty}

Where:

$\rho_{TD}$ = threading dislocation density (cm⁻²)
$\rho_0$ = initial density
$h_0$ = characteristic decay length
$\rho_{\infty}$ = residual density

6. Reactor-Scale Modeling

6.1 Coupled Transport Equations

6.1.1 Momentum Conservation (Navier-Stokes)

\rho\left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v}\right) = -\nabla p + \mu \nabla^2 \mathbf{v} + \rho \mathbf{g}

Where:

$\rho$ = gas density (kg/m³)
$\mathbf{v}$ = velocity vector (m/s)
$p$ = pressure (Pa)
$\mu$ = dynamic viscosity (Pa$\cdot$s)
$\mathbf{g}$ = gravitational acceleration

6.1.2 Continuity Equation

\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0

6.1.3 Species Transport

\frac{\partial C_i}{\partial t} + \mathbf{v} \cdot \nabla C_i = D_i \nabla^2 C_i + R_i

Where:

$C_i$ = concentration of species $i$ (mol/m³)
$D_i$ = diffusion coefficient of species $i$ (m²/s)
$R_i$ = net reaction rate (mol/m³$\cdot$s)

6.1.4 Energy Conservation

\rho c_p \left(\frac{\partial T}{\partial t} + \mathbf{v} \cdot \nabla T\right) = k \nabla^2 T + \sum_j \Delta H_j r_j

Where:

$c_p$ = specific heat capacity (J/kg$\cdot$K)
$k$ = thermal conductivity (W/m$\cdot$K)
$\Delta H_j$ = enthalpy of reaction $j$ (J/mol)
$r_j$ = rate of reaction $j$ (mol/m³$\cdot$s)

6.2 Silicon CVD Chemistry

6.2.1 From Silane (SiH₄)

Gas phase decomposition:

\text{SiH}_4 \xrightarrow{k_1} \text{SiH}_2 + \text{H}_2

Surface reaction:

\text{SiH}_2(g) + * \xrightarrow{k_2} \text{Si}(s) + \text{H}_2(g)

Where $*$ denotes a surface site.

6.2.2 From Dichlorosilane (DCS)

\text{SiH}_2\text{Cl}_2 \rightarrow \text{SiCl}_2 + \text{H}_2

\text{SiCl}_2 + \text{H}_2 \rightarrow \text{Si}(s) + 2\text{HCl}

6.2.3 Rate Law

r_{dep} = k_2 P_{SiH_2} (1 - \theta)

Where:

$P_{SiH_2}$ = partial pressure of SiH₂
$\theta$ = surface site coverage

6.3 Dimensionless Numbers

Number	Definition	Physical Meaning
Reynolds	$Re = \frac{\rho v L}{\mu}$	Inertia vs. viscous forces
Prandtl	$Pr = \frac{\mu c_p}{k}$	Momentum vs. thermal diffusivity
Schmidt	$Sc = \frac{\mu}{\rho D}$	Momentum vs. mass diffusivity
Damköhler	$Da = \frac{k_s L}{D}$	Reaction rate vs. diffusion rate
Grashof

u^2}$ | Buoyancy vs. viscous forces |

7. Selective Epitaxial Growth (SEG) Modeling

7.1 Overview

In SEG, growth occurs on exposed Si but not on dielectric (SiO₂/Si₃N₄).

7.2 Loading Effect Model

G_{local} = G_0 \left(1 + \alpha \cdot \frac{A_{mask}}{A_{Si}}\right)

Where:

$G_{local}$ = local growth rate
$G_0$ = baseline growth rate
$\alpha$ = pattern sensitivity factor
$A_{mask}$ = dielectric (mask) area
$A_{Si}$ = exposed silicon area

7.3 Pattern-Dependent Growth

Sources of non-uniformity:

Local depletion of reactants over Si regions
Species reflected/desorbed from mask contribute to nearby Si
Gas-phase diffusion length effects

7.4 Selectivity Condition

For selective growth on Si vs. oxide:

r_{deposition,Si} > 0 \quad \text{and} \quad r_{deposition,oxide} < r_{etching,oxide}

Achieved by adding HCl:

\text{Si}(nuclei) + 2\text{HCl} \rightarrow \text{SiCl}_2 + \text{H}_2

Nuclei on oxide are etched before they can grow, maintaining selectivity.

7.5 Faceting Model

Growth rate depends on crystallographic orientation:

G_{(hkl)} = G_0 \cdot f(hkl) \cdot \exp\left(-\frac{E_{a,(hkl)}}{k_B T}\right)

Typical growth rate hierarchy:

G_{(100)} > G_{(110)} > G_{(111)}

8. Dopant Incorporation

8.1 Segregation Coefficient

Equilibrium segregation coefficient:

k_0 = \frac{C_{solid}}{C_{liquid/gas}}

Effective segregation coefficient:

k_{eff} = \frac{k_0}{k_0 + (1-k_0)\exp\left(-\frac{G\delta}{D_l}\right)}

Where:

$k_0$ = equilibrium segregation coefficient
$G$ = growth rate
$\delta$ = boundary layer thickness
$D_l$ = diffusivity in liquid/gas phase

8.2 Dopant Concentration in Film

C_{film} = k_{eff} \cdot C_{gas}

8.3 Dopant Profile Abruptness

The transition width is limited by:

Surface segregation length: $\lambda_{seg}$
Diffusion during growth: $L_D = \sqrt{D \cdot t}$
Autodoping from substrate

\Delta z_{transition} \approx \sqrt{\lambda_{seg}^2 + L_D^2}

8.4 Common Dopants for Si Epitaxy

Dopant	Type	Precursor	Segregation Behavior
B	p-type	B₂H₆, BCl₃	Low segregation
P	n-type	PH₃, PCl₃	Moderate segregation
As	n-type	AsH₃	Strong segregation
Sb	n-type	SbH₃	Very strong segregation

9. Atomistic Simulation Methods

9.1 Kinetic Monte Carlo (KMC)

9.1.1 Event Rates

Each atomic event has a rate following Arrhenius:

\Gamma_i = u_0 \exp\left(-\frac{E_i}{k_B T}\right)

Where:

$\Gamma_i$ = rate of event $i$ (s⁻¹)
$

u_0$ = attempt frequency (~10¹²-10¹³ s⁻¹)

$E_i$ = activation energy for event $i$

9.1.2 Events Modeled

Adsorption: $\Gamma_{ads} = \frac{P}{\sqrt{2\pi m k_B T}} \cdot s$
Desorption: $\Gamma_{des} =

u_0 \exp(-E_{des}/k_B T)$

Surface diffusion: $\Gamma_{diff} =

u_0 \exp(-E_m/k_B T)$

Step attachment: $\Gamma_{attach}$
Step detachment: $\Gamma_{detach}$

9.1.3 Time Advancement

\Delta t = -\frac{\ln(r)}{\Gamma_{total}} = -\frac{\ln(r)}{\sum_i \Gamma_i}

Where $r$ is a uniform random number in $(0,1]$.

9.2 Density Functional Theory (DFT)

Provides input parameters for KMC:

Adsorption energies
Migration barriers
Surface reconstruction energetics
Reaction pathways

Kohn-Sham equation:

\left[-\frac{\hbar^2}{2m}\nabla^2 + V_{eff}(\mathbf{r})\right]\psi_i(\mathbf{r}) = \varepsilon_i \psi_i(\mathbf{r})

9.3 Molecular Dynamics (MD)

Newton's equations:

m_i \frac{d^2 \mathbf{r}_i}{dt^2} = -\nabla_i U(\mathbf{r}_1, \mathbf{r}_2, ..., \mathbf{r}_N)

Where $U$ is the interatomic potential (e.g., Stillinger-Weber, Tersoff for Si).

10. Nucleation Theory

10.1 Classical Nucleation Theory (CNT)

10.1.1 Gibbs Free Energy Change

\Delta G(r) = -\frac{4}{3}\pi r^3 \cdot \frac{\Delta \mu}{\Omega} + 4\pi r^2 \gamma

Where:

$r$ = nucleus radius
$\Delta \mu$ = supersaturation (driving force)
$\Omega$ = atomic volume
$\gamma$ = surface energy

10.1.2 Critical Nucleus Radius

Setting $\frac{d(\Delta G)}{dr} = 0$:

r^* = \frac{2\gamma \Omega}{\Delta \mu}

10.1.3 Free Energy Barrier

\Delta G^* = \frac{16 \pi \gamma^3 \Omega^2}{3 (\Delta \mu)^2}

10.1.4 Nucleation Rate

J = Z \beta^* N_s \exp\left(-\frac{\Delta G^*}{k_B T}\right)

Where:

$J$ = nucleation rate (nuclei/cm²$\cdot$s)
$Z$ = Zeldovich factor (~0.01-0.1)
$\beta^*$ = attachment rate to critical nucleus
$N_s$ = surface site density

10.2 Growth Modes

Mode	Surface Energy Condition	Growth Behavior	Example
Frank-van der Merwe	$\gamma_s \geq \gamma_f + \gamma_{int}$	Layer-by-layer (2D)	Si on Si
Volmer-Weber	$\gamma_s < \gamma_f + \gamma_{int}$	Island (3D)	Metals on oxides
Stranski-Krastanov	Intermediate	2D then 3D islands	InAs/GaAs QDs

10.3 2D Nucleation

Critical island size (atoms):

i^* = \frac{\pi \gamma_{step}^2 \Omega}{(\Delta \mu)^2 k_B T}

11. TCAD Process Simulation

11.1 Overview

Tools: Synopsys Sentaurus Process, Silvaco Victory Process

11.2 Diffusion-Reaction System

\frac{\partial C_i}{\partial t} = \nabla \cdot (D_i \nabla C_i - \mu_i C_i \nabla \phi) + G_i - R_i

Where:

First term: Fickian diffusion
Second term: Drift in electric field (for charged species)
$G_i$ = generation rate
$R_i$ = recombination rate

11.3 Point Defect Dynamics

Vacancy concentration:

\frac{\partial C_V}{\partial t} = D_V \nabla^2 C_V + G_V - k_{IV} C_I C_V

Interstitial concentration:

\frac{\partial C_I}{\partial t} = D_I \nabla^2 C_I + G_I - k_{IV} C_I C_V

Where $k_{IV}$ is the recombination rate constant.

11.4 Stress Evolution

Equilibrium equation:

\nabla \cdot \boldsymbol{\sigma} = 0

Constitutive relation:

\boldsymbol{\sigma} = \mathbf{C} : (\boldsymbol{\varepsilon} - \boldsymbol{\varepsilon}^{thermal} - \boldsymbol{\varepsilon}^{intrinsic})

Where:

$\boldsymbol{\sigma}$ = stress tensor
$\mathbf{C}$ = elastic stiffness tensor
$\boldsymbol{\varepsilon}$ = total strain
$\boldsymbol{\varepsilon}^{thermal}$ = thermal strain = $\alpha \Delta T$
$\boldsymbol{\varepsilon}^{intrinsic}$ = intrinsic strain (lattice mismatch)

11.5 Level Set Method for Interface Tracking

\frac{\partial \phi}{\partial t} + v_n |\nabla \phi| = 0

Where:

$\phi$ = level set function (interface at $\phi = 0$)
$v_n$ = interface normal velocity

12. Advanced Topics

12.1 Atomic Layer Epitaxy (ALE) / Atomic Layer Deposition (ALD)

Self-limiting surface reactions modeled as Langmuir kinetics:

\theta = \frac{K \cdot P \cdot t}{1 + K \cdot P \cdot t} \rightarrow 1 \quad \text{as } t \rightarrow \infty

Growth per cycle (GPC):

GPC = \theta_{sat} \cdot d_{monolayer}

Typical GPC values: 0.5-1.5 Å/cycle

12.2 III-V on Silicon Integration

Challenges and models:

Anti-phase boundaries (APBs): Form at single-step terraces
Threading dislocations: $\rho_{TD} \propto f^2$ initially
Thermal mismatch stress: $\sigma_{thermal} = \frac{E \Delta \alpha \Delta T}{1-

u}$

12.3 Quantum Dot Formation (Stranski-Krastanov)

Critical thickness for islanding:

h_{SK} \approx \frac{\gamma}{M f^2}

Island density:

n_{island} \propto \exp\left(-\frac{E_{island}}{k_B T}\right) \cdot F^{1/3}

Where $F$ is the deposition flux.

12.4 Machine Learning in Epitaxy Modeling

Physics-Informed Neural Networks (PINNs):

\mathcal{L}_{total} = \mathcal{L}_{data} + \lambda_{PDE}\mathcal{L}_{physics} + \lambda_{BC}\mathcal{L}_{boundary}

Where:

$\mathcal{L}_{data}$ = data fitting loss
$\mathcal{L}_{physics}$ = PDE residual loss
$\mathcal{L}_{boundary}$ = boundary condition loss
$\lambda$ = weighting parameters

Applications:

Surrogate models for reactor optimization
Inverse problems (parameter extraction)
Process window optimization
Defect prediction

13. Key Equations

Phenomenon	Key Equation	Primary Parameters
Growth rate (dual regime)	$G = \frac{k_s h_g C_g}{k_s + h_g}$	Temperature, pressure, flow
Surface diffusion length	$\lambda_s = \sqrt{D_s \tau_s}$	Temperature
Lattice mismatch	$f = \frac{a_f - a_s}{a_s}$	Material system
Critical thickness

u)}\left[\ln\frac{h_c}{b}+1\right]$ | Mismatch, Burgers vector |

Elastic strain energy	$E = M f^2 h$	Mismatch, thickness, modulus
Nucleation rate	$J \propto \exp(-\Delta G^*/k_BT)$	Supersaturation, surface energy
Species transport	$\frac{\partial C}{\partial t} + \mathbf{v}\cdot\nabla C = D\nabla^2 C + R$	Diffusivity, velocity, reactions
KMC event rate

u_0 \exp(-E_a/k_BT)$ | Activation energy, temperature |

Physical Constants

Constant	Symbol	Value
Boltzmann constant	$k_B$	$1.38 \times 10^{-23}$ J/K
Gas constant	$R$	8.314 J/mol$\cdot$K
Planck constant	$h$	$6.63 \times 10^{-34}$ J$\cdot$s
Electron charge	$e$	$1.60 \times 10^{-19}$ C
Si lattice constant	$a_{Si}$	5.431 Å
Ge lattice constant	$a_{Ge}$	5.658 Å
GaAs lattice constant	$a_{GaAs}$	5.653 Å

Single-Wafer Epi Reactor Cross-Section. The dominant reactor architecture for advanced logic epi is the cold-wall, single-wafer, lamp-heated chamber — a design that can ramp from 400$^\circ$C to 1150$^\circ$C in under 30 seconds and cool back in 60 seconds, enabling the tight thermal budgets that GAA nanosheet superlattices demand. Gas enters through a horizontal injector, flows across the rotating wafer in a laminar sheet, and exits through an exhaust port on the opposite side. The wafer sits on a SiC-coated graphite susceptor heated by banks of tungsten-halogen lamps above and below the quartz chamber walls.

CVD Chemistry and Thermal Budget. The epi process is a chemical vapor deposition where the substrate temperature determines whether growth is limited by surface kinetics (low T) or by gas-phase mass transport (high T). At 1050–1150$^\circ$C with trichlorosilane (SiHCl$_3$/H$_2$), the growth rate reaches 0.5–4 $\mu$m/min in the mass-transport regime — fast enough for substrate epi layers of 2–10 $\mu$m in under 5 minutes. At 500–700$^\circ$C with dichlorosilane (SiH$_2$Cl$_2$/GeH$_4$/HCl), growth drops to 5–30 nm/min but enables selective epitaxy: HCl etches polycrystalline nuclei on oxide and nitride surfaces while leaving single-crystal growth on exposed silicon intact, achieving selectivity exceeding 100:1. The thermal budget constraint is absolute — at 900$^\circ$C, Ge interdiffusion reaches 1.9 nm/hour ($D = 10^{-17}$ cm$^2$/s), which smears a 5 nm SiGe/Si interface into a graded transition. This is why GAA nanosheet superlattice epi operates at 500–650$^\circ$C despite the 10$\times$ lower growth rate: interface abruptness ($<$1 nm transition width) matters more than throughput for the 2 nm node.

GAA Nanosheet Superlattice — The Defining Epi Challenge of the 2 nm Node. Gate-all-around transistors require a perfectly periodic Si/SiGe superlattice grown by epitaxy: alternating layers of Si (5–7 nm, future channel) and Si$_{0.7}$Ge$_{0.3}$ (8–12 nm, sacrificial, later removed by selective etch). A typical 2 nm node stack has 4–8 pairs, totaling 60–120 nm, with each layer controlled to $\pm$0.3 nm thickness and Ge composition held at $\pm$1 atomic percent uniformity across 300 mm. The growth sequence alternates SiH$_4$ (Si layers) and SiH$_2$Cl$_2$/GeH$_4$ (SiGe layers) at 500–650$^\circ$C with H$_2$ carrier gas at 10–80 Torr. Interface abruptness demands gas switching in $<$2 seconds (purge between layers) and minimal thermal exposure after growth. Any interdiffusion wider than 1 nm at the Si/SiGe boundary creates a graded composition that shifts the selective etch endpoint by 0.5–2 nm — enough to fail the channel thickness specification.

Selective Epitaxial Growth (SEG) for Raised Source/Drain. Since the 45 nm node, PMOS transistors use compressive-strained SiGe grown selectively in etched recesses adjacent to the gate. The epi fills only the crystalline Si surface while HCl in the gas stream continuously etches any polycrystalline nuclei that form on surrounding SiO$_2$ or Si$_3$N$_4$ — achieving greater than 100:1 selectivity without a mask. At 650$^\circ$C with SiH$_2$Cl$_2$/GeH$_4$/HCl/H$_2$ at 10 Torr, growth proceeds at 10–20 nm/min with Ge content of 25–50 atomic percent. The resulting SiGe exerts uniaxial compressive strain on the Si channel between the source/drain regions, increasing hole mobility by 50–100% — the performance boost that kept planar CMOS scaling alive from 90 nm through 22 nm before FinFET took over.

EPI Chamber Cross-Section — Hardware Subsystems. The single-wafer epi chamber integrates thermal, chemical, mechanical, and optical subsystems into a compact cold-wall reactor optimized for rapid thermal cycling. Unlike plasma etch chambers that require RF generators and vacuum in the millitorr regime, the epi reactor operates at 10–760 Torr with purely thermal activation — the lamp bank replaces the plasma as the energy source for breaking precursor bonds on the wafer surface.

EPI Chamber Control Schematic — Temperature, Gas, and Thickness Feedback. The epi reactor is a multi-input, multi-output control system where lamp power (11+ independent zones), gas flow (4–8 MFC channels), susceptor rotation, and chamber pressure must all coordinate to deliver ±0.3 nm thickness and ±1 at% composition uniformity. Unlike etch where RF power is the primary control variable, epi control is dominated by temperature — because growth rate has an Arrhenius dependence with activation energy 1.5–3.0 eV, meaning a 1°C error at 600°C changes the rate by 0.3–0.5%.

EPI Chamber Process Environment — No Plasma, Pure Thermal Activation. Unlike etch and PECVD chambers that use plasma to dissociate precursors, the epi reactor relies entirely on thermal energy at the wafer surface to decompose gas molecules. At 1050°C, SiHCl$_3$ pyrolyzes on the Si surface with an activation energy of 1.8 eV — the surface temperature provides sufficient energy to break the Si–H and Si–Cl bonds, releasing HCl and incorporating Si into the crystal lattice. At 600°C for selective SiGe, the lower activation energy of GeH$_4$ decomposition (0.9 eV) enables Ge incorporation while SiH$_2$Cl$_2$ decomposition (1.5 eV) proceeds more slowly — this differential sets the Ge fraction. The absence of plasma means no ion bombardment, no radiation damage, no charging — enabling perfect crystalline growth with defect densities below $10^2$ cm$^{-2}$.

EPI Process Metrics — What the Fab Measures. The epi process is qualified by six metrics that collectively determine whether the grown layer meets transistor specifications: (1) thickness uniformity ($\pm$0.5% across 300 mm, $\pm$0.3 nm for nanosheets), (2) composition uniformity (Ge $\pm$1 at% for SiGe), (3) defect density ($<$0.1 defects/cm$^2$ for substrate epi, $<$100/cm$^2$ for selective), (4) resistivity uniformity ($\pm$3% for doped layers), (5) surface roughness ($<$0.1 nm RMS by AFM), and (6) interface abruptness ($<$1 nm Ge transition for superlattice). Metrology uses spectroscopic ellipsometry (thickness/composition), four-point probe (resistivity), haze inspection (particles), X-ray diffraction (strain/composition), and cross-section TEM (interface verification). Every wafer gets inline ellipsometry; TEM sampling runs at 1 per lot (25 wafers) for process monitoring.

EPI Typical Process — Step-by-Step Sequence. A representative selective SiGe source/drain epi process on a 300 mm wafer runs the following sequence in a single chamber recipe lasting 4–8 minutes total: (1) Load wafer onto susceptor at 400°C, purge chamber with H$_2$ at 100 Torr (30 s). (2) Ramp to 800°C in H$_2$ for pre-bake surface clean — removes native oxide via H$_2$ reduction (60 s). (3) Cool to 650°C stabilization temperature (30 s). (4) Introduce SiH$_2$Cl$_2$ + GeH$_4$ + HCl + B$_2$H$_6$ (dopant) at 10 Torr — selective growth at 15 nm/min, Ge = 35%, boron $2 \times 10^{20}$ cm$^{-3}$ (120–240 s for 30–60 nm). (5) Purge all precursors, ramp to 700°C for 10 s anneal (optional, for dopant activation). (6) Cool to 400°C, unload (60 s). Total thermal budget: 650°C peak for 4 minutes — compatible with HKMG gate-last integration. Chamber conditions between wafers: 30 s H$_2$ purge + lamp idle at 400°C. Throughput: 5–8 WPH per chamber, 20–32 WPH on a 4-chamber cluster.

EPI Typical Productivity Improvements (2015–2024). The epi equipment industry has delivered consistent productivity gains through hardware and process innovation, reducing cost-per-wafer-pass by approximately 8% per year: (1) Multi-zone lamp PID (11→16 zones) reduced center-to-edge temperature non-uniformity from $\pm$3°C to $\pm$1°C, eliminating the need for rework and increasing first-pass yield from 92% to 99%. (2) Fast gas switching valves ($<$1 s actuation vs $<$5 s legacy) enabled superlattice epi without dedicated purge steps, cutting cycle time by 15%. (3) Higher lamp power density (100 kW peak vs 60 kW) enabled 30 s ramps instead of 60 s — directly adding 30 s throughput per wafer. (4) In-situ reflectometry closed the thickness loop wafer-by-wafer, reducing the metrology burden and enabling APC (advanced process control) that compensates for susceptor aging across 5,000-wafer campaigns. (5) Selective epi without the pre-clean step (replacing ex-situ HF dip with in-situ H$_2$ bake at 800°C) eliminated a wet-bench tool from the flow, saving 2 minutes of queue time and one cross-contamination risk. (6) Cluster tool integration (epi + pre-clean + cool-down in vacuum) removed atmosphere exposure between steps, reducing interface oxygen from $10^{12}$ to $<10^{10}$ atoms/cm$^2$. Net result: cost per epi wafer-pass dropped from $\sim$45 USD (2015) to $\sim$28 USD (2024) while the process specification tightened 3$\times$ — the definition of a mature equipment learning curve.

epiepitaxyepitaxialepitaxial layerepi layerepi process

Explore 500+ Semiconductor & AI Topics

From EUV lithography to CUDA optimization — search the full knowledge base or chat with our AI assistant.

🔍 Search Topics 💬 Ask CFSGPT 📚 Browse All

Related Topics

Explore 500+ Semiconductor & AI Topics