📚 课程目录第3章器件仿真基础

3. 器件仿真基础：漂移-扩散模型、载流子传输方程、PN结与MOSFET的SDevice仿真设置

本章是TCAD器件仿真的核心理论基础与实操起点。我们将从半导体器件物理的基本方程出发，深入理解漂移-扩散模型，并掌握Sentaurus SDevice中PN结与MOSFET的仿真设置流程。准确理解这些基础是后续进行参数校准与提取的前提。

3.1 漂移-扩散模型与载流子传输方程

在TCAD器件仿真中，漂移-扩散模型是最经典且应用最广泛的载流子传输模型。它描述了电子和空穴在半导体中如何在外加电场和浓度梯度作用下运动。该模型由三个核心方程组成：泊松方程、电子连续性方程和空穴连续性方程。

3.1.1 核心方程体系

方程名称	物理意义	数学表达式
泊松方程	描述静电势与空间电荷密度的关系	`∇²ψ = -q/ε (p - n + N_D⁺ - N_A⁻)`
电子连续性方程	描述电子浓度随时间的变化	`∂n/∂t = (1/q) ∇·J_n - R_net`
空穴连续性方程	描述空穴浓度随时间的变化	`∂p/∂t = -(1/q) ∇·J_p - R_net`

其中，ψ 为静电势，n 和 p 分别为电子和空穴浓度，J_n 和 J_p 为电子和空穴电流密度，R_net 为净复合率。

3.1.2 电流密度方程

在漂移-扩散模型中，电流密度由漂移分量和扩散分量线性叠加而成：

电子电流密度： J_n = q μ_n n E + q D_n ∇n
空穴电流密度： J_p = q μ_p p E - q D_p ∇p

式中，μ 为迁移率，D 为扩散系数，二者通过爱因斯坦关系关联：D/μ = k_BT/q。

    关键理解： 漂移项（q μ n E）由电场驱动，扩散项（q D ∇n）由浓度梯度驱动。在短沟道器件中，扩散分量往往占据主导地位。

仿真提示： 在SDevice中，默认即采用漂移-扩散模型。对于先进小尺寸器件（如亚10nm节点），可能需要开启量子修正模型（如密度梯度模型），但漂移-扩散模型仍是所有高级模型的基础。

3.1.3 漂移-扩散模型求解流程

3.2 PN结的SDevice仿真设置

PN结是最基本的半导体器件结构，其I-V特性（整流特性）是验证仿真模型正确性的第一步。下面以一个均匀掺杂的硅PN结为例，展示SDevice的完整设置流程。

3.2.1 器件结构与网格定义

在SDevice中，器件结构通常由Sentaurus Structure Editor (SDE) 生成，或通过命令直接定义。以下是一个1D PN结的SDevice输入文件片段：

Device PN_Junction {
    Electrode {
        { Name="anode"    Voltage=0.0 }
        { Name="cathode"  Voltage=0.0 }
    }
    File {
        Grid = "pn_1d_msh.tdr"
        Doping = "pn_1d_msh.tdr"
    }
    Physics {
        AreaFactor = 1e-4
        Mobility( DopingDependence )
        Recombination( SRH(DopingDependence) Auger )
    }
    Plot {
        eDensity hDensity Potential ElectricField
        eMobility hMobility SRH Auger
    }
    CurrentPlot {
        Electrode( [anode cathode] )
    }
}

常见错误： 忘记定义 AreaFactor 会导致电流值偏离预期。对于2D仿真，AreaFactor 用于将单位深度的电流转换为实际电流值。

3.2.2 偏置扫描与I-V特性提取

PN结的I-V特性通过扫描阳极电压获得。以下为SDevice中的Math和Solve部分设置：

Math {
    NumberOfThreads = 4
    Extrapolate
    RelErrControl
    Digits = 5
    Notdamped = 50
    Iterations = 20
}
Solve {
    Coupled( Iterations=100 ) { Poisson }
    Coupled( Iterations=100 ) { Poisson Electron Hole }
    Quasistationary (
        InitialStep=0.01 Increment=1.5 MinStep=1e-8 MaxStep=0.05
        Goal { Parameter=anode Voltage=1.0 }
    ) { Coupled { Poisson Electron Hole } }
}

3.2.3 PN结能带图与内建电势

_C E_V E_F V_bi 空间电荷区

3.3 MOSFET的SDevice仿真设置

MOSFET是现代集成电路的核心器件。其仿真设置比PN结复杂，需要定义栅极、源极、漏极和衬底四个电极，并考虑栅氧化层和沟道掺杂分布。

3.3.1 器件结构与关键参数

参数名称	典型值	说明
沟道长度 (L_g)	45 nm	栅极长度，决定器件开关速度
栅氧化层厚度 (T_ox)	2 nm	等效氧化层厚度，影响栅控能力
沟道掺杂浓度 (N_ch)	1×10¹⁸ cm^-3	P型衬底掺杂，调节阈值电压
源/漏掺杂浓度 (N_S/D)	1×10²⁰ cm^-3	重掺杂N型，形成欧姆接触
栅极功函数 (Φ_M)	4.5 eV	金属栅极功函数，影响V_th

3.3.2 SDevice输入文件核心设置

以下是一个45nm NMOSFET的SDevice仿真设置片段，重点展示物理模型和接触定义：

Device NMOS {
    Electrode {
        { Name="source"   Voltage=0.0 }
        { Name="drain"    Voltage=0.0 }
        { Name="gate"     Voltage=0.0 }
        { Name="substrate" Voltage=0.0 }
    }
    File {
        Grid = "nmos_45nm_msh.tdr"
        Doping = "nmos_45nm_msh.tdr"
    }
    Physics {
        AreaFactor = 1e-4
        Mobility( DopingDependence HighFieldSaturation Enormal )
        EffectiveIntrinsicDensity( BandGapNarrowing( OldSlotboom ) )
        Recombination( SRH( DopingDependence ) Auger )
        Fermi
    }
    Physics( Material="Silicon" ) {
        Piezo( kT = 0.0259 )
    }
    Contact( gate ) {
        Workfunction = 4.5
    }
    Plot {
        eDensity hDensity Potential ElectricField
        eMobility hMobility eVelocity hVelocity
        SRH Auger BandGap
    }
    CurrentPlot {
        Electrode( [source drain gate substrate] )
    }
}

物理模型选择： 对于短沟道MOSFET，必须开启 HighFieldSaturation（高场饱和）和 Enormal（垂直电场迁移率退化）模型，否则会严重高估饱和电流。

3.3.3 转移特性与输出特性扫描

MOSFET的特性仿真通常包含两组扫描：

转移特性 (I_D-V_G)： 固定V_D（如0.05V和1.0V），扫描V_G从0到1.2V，提取阈值电压V_th、亚阈值摆幅SS。
输出特性 (I_D-V_D)： 固定V_G（如0.4V, 0.6V, 0.8V, 1.0V），扫描V_D从0到1.2V，提取饱和电流I_Dsat。

Solve {
    * 初始解
    Coupled( Iterations=100 ) { Poisson }
    Coupled( Iterations=100 ) { Poisson Electron Hole }

    * 转移特性扫描 (Vd=0.05V)
    Quasistationary (
        InitialStep=0.01 Increment=1.5 MinStep=1e-8 MaxStep=0.05
        Goal { Parameter=gate Voltage=1.2 }
    ) {
        Coupled { Poisson Electron Hole }
    }

    * 输出特性扫描 (Vg=1.0V)
    Set( gate Voltage=1.0 )
    Quasistationary (
        InitialStep=0.01 Increment=1.5 MinStep=1e-8 MaxStep=0.05
        Goal { Parameter=drain Voltage=1.2 }
    ) {
        Coupled { Poisson Electron Hole }
    }
}

3.3.4 MOSFET沟道电势与载流子分布

_S V_D V_G V_B

3.4 仿真收敛性调试与常见问题

在SDevice仿真中，收敛失败是最常见的工程问题。以下是一些实用的调试策略：

问题现象	可能原因	解决方案
Newton迭代不收敛	初始步长过大；网格质量差	减小 `InitialStep`；加密网格
电流不连续或跳变	网格过渡不均匀；接触定义错误	检查网格渐变比；确认电极名称与SDE一致
负微分电阻出现	自热效应未开启；迁移率模型不匹配	添加 `Thermodynamic` 模型；检查 `HighFieldSaturation`
阈值电压偏移过大	功函数设置错误；界面陷阱未考虑	校准功函数；添加 `FixedCharge` 或 `InterfaceTraps`

    工程经验： 在开始复杂器件仿真前，始终先运行一个简单的PN结或长沟道MOSFET（Lg > 1μm）来验证物理模型和数值设置的正确性。这可以节省大量调试时间。

3.5 本章小结

漂移-扩散模型是TCAD器件仿真的基石，由泊松方程和载流子连续性方程耦合构成。
PN结仿真验证了基本的整流特性和内建电势，是模型校准的第一步。
MOSFET仿真需要精确设置栅极功函数、迁移率模型和接触定义，才能准确预测I-V特性。
收敛性调试是仿真工程师的核心技能，需要结合物理理解和数值技巧。