Mobility Modeling is the TCAD simulation of charge carrier drift mobility (μ) as a function of doping concentration, electric field, temperature, interface quality, and crystal strain — predicting the carrier transport speed that determines transistor drive current (I_on), switching speed (f_T), and energy efficiency, using Matthiessen's Rule to combine the independent contributions of phonon scattering, ionized impurity scattering, surface roughness scattering, and other mechanisms into a total effective mobility.
What Is Carrier Mobility?
Mobility quantifies how fast a carrier drifts in response to an electric field:
μ = v_drift / E (units: cm²/V·s)
Higher mobility → faster carrier response → faster transistor switching at lower supply voltage.
Matthiessen's Rule — Combining Scattering Mechanisms
Each scattering mechanism independently limits mobility. The total mobility is their harmonic sum:
1/μ_total = 1/μ_phonon + 1/μ_impurity + 1/μ_surface + 1/μ_other
The mechanism with the lowest individual mobility dominates the total (bottleneck principle).
Low-Field Mobility Models
Phonon Scattering Component (μ_phonon): Acoustic and optical phonon scattering dominate in lightly doped silicon at room temperature. Temperature dependence follows μ_phonon ∝ T^(-3/2) for acoustic phonons — mobility degrades with increasing temperature, the fundamental reason processor performance drops under thermal throttling.
Ionized Impurity Scattering Component (μ_imp): Coulomb interaction with ionized donor and acceptor atoms. Concentration dependence modeled by Masetti et al.:
μ = μ_min + (μ_max - μ_min) / (1 + (N/N_ref)^α)
Where N = total ionized impurity concentration. Mobility drops sharply above ~10¹⁷ cm⁻³ doping — the key trade-off between conductivity (needs high doping) and mobility (degraded by high doping).
Surface Roughness Scattering Component (μ_sr): Dominates in the MOSFET inversion layer under high vertical fields. The Lombardi model adds a field-dependent surface mobility component:
μ_sr ∝ 1/(E_perp)² × 1/δ_rms²
Where E_perp = perpendicular field and δ_rms = oxide interface roughness amplitude. As gate overdrive increases, E_perp increases, confining carriers tighter against the rough interface → mobility decreases. This "mobility degradation" is why measured MOSFET mobility peaks at low gate voltage and falls at high VGS.
High-Field Velocity Saturation
At high lateral electric fields, carriers emit optical phonons faster than they gain energy from the field — reaching a saturation velocity:
v_sat(Si electrons) ≈ 10⁷ cm/s
The Caughey-Thomas model transitions smoothly from ohmic to saturated velocity:
v(E) = μ_low × E / [1 + (μ_low × E / v_sat)^β]^(1/β)
Velocity saturation is the fundamental limit of drive current in nanometer-scale transistors where the entire channel is near saturation.
Quantum Confinement Corrections
In FinFETs and nanosheet FETs with body thickness < 10 nm, quantum confinement shifts the energy subbands and modifies carrier occupancy relative to bulk. Effective mass and density of states corrections to the mobility model are required to avoid overestimating drive current.
Why Mobility Modeling Matters
- Drive Current Prediction: I_on ∝ μ × Cox × (VGS - Vth) × V_drain for long channel. Mobility accuracy directly determines drive current prediction accuracy — 10% mobility error → 10% drive current error → incorrect power/performance model.
- Process Optimization: Simulation-guided mobility optimization identifies the trade-off between higher channel doping (needed to suppress short-channel effects) and lower channel mobility (consequence of higher impurity scattering). Finding the optimal pocket implant dose requires accurate mobility modeling.
- Strain Engineering Validation: The mobility enhancement from strained silicon channels must be accurately predicted to justify the process integration cost. Piezoresistance models and band structure-derived mobility enhancements are validated against measurement in simulation.
- Self-Heating Coupling: In FinFETs at high power density, junction temperature rises substantially. Since μ_phonon ∝ T^(-3/2), self-heating reduces carrier mobility, further reducing drive current — a negative feedback that simulation must capture for accurate I_on–I_off modeling under realistic operating conditions.
Tools
- Synopsys Sentaurus Device: Full mobility model library including Masetti, Lombardi surface model, high-field saturation, quantum correction, and strain-dependent piezoresistance.
- Silvaco Atlas: Device simulator with comprehensive mobility models for Si, SiGe, Ge, III-V materials.
- nextnano: k·p-based quantum transport simulation including mobility in nanostructures.
Mobility Modeling is calculating the speed limit for charge carriers — summing all the scattering forces that impede carrier drift through the transistor channel to predict the drive current and switching speed that determine whether a chip delivers its target performance, guiding process engineers to the optimal combination of doping, strain, interface quality, and geometry that maximizes carrier speed at minimum power consumption.