Kinetic Monte Carlo (KMC) is a stochastic simulation method that models the time evolution of a system by statistically sampling transitions between discrete states based on their transition rates — enabling simulation of diffusion, crystal growth, defect annealing, and surface phenomena over timescales of microseconds to hours that Molecular Dynamics (limited to nanoseconds) cannot reach, while preserving the atomic-scale resolution that continuum models sacrifice.
What Is Kinetic Monte Carlo?
KMC treats the system as a collection of possible events (atomic transitions), each with a rate derived from physics:
The KMC Algorithm (Bortz-Kalos-Lebowitz, BKL method):
1. Catalog Events: List all possible transitions from the current system state. Each event i has a rate Rᵢ (units: s⁻¹), computed from an Arrhenius expression: Rᵢ = ν₀ exp(−Eₐ/kT), where ν₀ is an attempt frequency (~10¹³ s⁻¹ for lattice vibrations), Eₐ is the activation energy, k is Boltzmann's constant, T is temperature.
2. Select Event: Choose event j with probability Rⱼ / ΣRᵢ (normalized by total rate) using a uniform random number.
3. Advance Time: The time increment is drawn from an exponential distribution: Δt = −ln(u) / ΣRᵢ, where u is a uniform random number. This ensures Poisson-distributed event times consistent with the physical process.
4. Update State: Execute the selected event (move atom, form cluster, annihilate defect pair).
5. Repeat: Accumulate statistics over millions of KMC steps.
Why KMC Bridges the Timescale Gap
The fundamental challenge in semiconductor simulation is the gap between:
- MD timescales (~10 ns maximum): Too short to observe diffusion at processing temperatures.
- Continuum TCAD timescales (~seconds to hours): Accurate for gradual processes but loses atomic-scale mechanism.
KMC fills this gap by advancing time event-by-event rather than step-by-step at fixed time increments. When the system state is static (no events occur for long periods), KMC idle time is skipped automatically — allowing rapid simulation of arbitrarily long time periods while maintaining atomic resolution during the active events.
Applications in Semiconductor Processing
Transient Enhanced Diffusion (TED):
The primary application in TCAD. Implant damage creates excess silicon interstitials that form clusters ({311} defects, Frank loops). KMC tracks the emission of single interstitials from these clusters, their diffusion to the surface, and their enhancement of dopant diffusion. KMC TED models provide the physical basis for the empirical parameters in commercial TCAD diffusion simulators.
Thin Film Deposition (CVD/ALD/MBE):
Adsorption, surface diffusion, nucleation island formation, and layer-by-layer vs. 3D growth transitions are naturally simulated by KMC on a surface lattice — capturing roughness evolution and step flow dynamics that continuum models of film growth cannot resolve.
Dopant-Defect Cluster Evolution:
Formation and dissolution of boron-interstitial clusters (BnIm), phosphorus-vacancy clusters, and arsenic clusters during annealing determine the fraction of electrically active dopant. KMC directly simulates cluster growth/shrinkage kinetics.
Electromigration in Interconnects:
Void nucleation and growth in copper interconnects under electromigration stress is a discrete event process accurately modeled by KMC with activation energies derived from DFT.
Coupling in the Multiscale Hierarchy
KMC occupies the critical middle layer in semiconductor multiscale simulation:
DFT/MD → (activation energies, attempt frequencies) → KMC → (effective diffusivities, cluster size distributions) → Continuum TCAD
Tools
- DADOS / DADOS3D: University of Murcia KMC simulator for dopant-defect interaction in silicon — widely used in academic TED research.
- LKMC (Lattice KMC): Generic framework for surface growth and diffusion simulations.
- Synopsys Sentaurus Process (KMC mode): Commercial TCAD with KMC-based diffusion for advanced node TED and cluster simulation.
Kinetic Monte Carlo is simulating time by jumping between events — the stochastic method that bridges the nanosecond limit of molecular dynamics and the second-scale reach of continuum models, preserving atomic-scale physics while enabling simulation of the microsecond-to-millisecond thermal processes that govern dopant activation and diffusion in modern semiconductor manufacturing.