2. 器件物理模型体系总览

在 Sentaurus 中做器件仿真,第一步不是“把网格画细”,而是先决定用哪一层物理模型。模型层级决定了方程组的未知量个数、耦合强度和求解难度,也直接决定了你能否在合理时间内得到可信结果。Sentaurus 的器件物理模型大致可以按“载流子输运描述的精细程度”排成一条阶梯:漂移扩散(Drift-Diffusion, DD)→ 热力学(Thermodynamic, TD)→ 流体动力学(Hydrodynamic, HD)→ 蒙特卡洛/玻尔兹曼求解(MC/BTE)。越往上,物理越精细,但未知量、非线性耦合和计算代价也急剧上升。

一、模型层级的物理本质:未知量与方程

漂移扩散模型求解三个基本未知量:静电势 ψ、电子浓度 n、空穴浓度 p,对应泊松方程与两个连续性方程。载流子温度默认等于晶格温度,迁移率用场依赖模型(如 Canali、Lombardi)修正。它假设载流子与晶格处于热平衡,因此无法描述速度过冲、能量输运等非平衡效应。

热力学模型在 DD 基础上增加晶格热流方程,未知量增加晶格温度 TL。它考虑自加热效应,适合功率器件、大电流密度下的热效应分析,但载流子仍假设与晶格热平衡。

流体动力学模型进一步把电子、空穴视为独立流体,引入载流子温度 Tn、Tp,未知量增至 5~6 个。它能描述载流子能量输运、速度过冲、碰撞电离增强等效应,是亚 0.1 μm 器件、雪崩击穿、ESD 仿真的常用选择。

蒙特卡洛方法直接求解玻尔兹曼输运方程,不依赖迁移率模型,物理最精确,但计算代价极高,通常只用于小尺寸器件或关键物理验证。

Sentaurus 器件物理模型层级与计算代价 物理精细度 ↑ 蒙特卡洛 / 玻尔兹曼求解 未知量:分布函数 f(r,k,t)|无迁移率假设 适用:小尺寸器件物理验证、热载流子 流体动力学 (HD) 未知量:ψ, n, p, Tn, Tp|5~6 个耦合方程 适用:亚 0.1 μm、雪崩击穿、ESD、速度过冲 热力学 (TD) 未知量:ψ, n, p, TL|增加晶格热流方程 适用:功率器件自加热、大电流热效应 漂移扩散 (DD) 未知量:ψ, n, p|3 个基本方程 适用:常规 CMOS、长沟道、初步探索 计算代价与耦合强度 → 计算代价:极高(小时~天) 内存:极大|收敛:困难|统计噪声大 通常仅用于关键物理点验证 计算代价:高(分钟~小时) 内存:大|收敛:较难|需良好初值 适合中小规模 2D/3D 结构 计算代价:中(秒~分钟) 内存:中|收敛:中等|热耦合增加非线性 适合功率器件与热分析 计算代价:低(秒级) 内存:小|收敛:容易|鲁棒性好 适合快速迭代与参数扫描 选择原则:先用 DD 快速探索 → 需要热效应加 TD → 需要非平衡输运换 HD → 关键物理点用 MC 验证 模型越精细 ≠ 结果越正确:参数校准质量往往比模型层级更影响最终精度

二、各模型的适用场景与选择依据

选择模型时,先问三个问题:器件特征尺寸是否接近或小于载流子能量弛豫长度(约 50~100 nm)?是否存在显著自加热或高电流密度?是否需要描述碰撞电离、雪崩击穿等非平衡效应?

对于常规 180 nm 及以上 CMOS 工艺的 I-V 特性、阈值电压、亚阈值摆幅分析,DD 模型配合 Lombardi 迁移率模型和 SRH 复合模型通常足够,且收敛快、参数少。对于功率 LDMOS 或 IGBT 的导通电阻与热分布分析,TD 模型能给出晶格温度分布,帮助判断热失效点。对于 40 nm 以下器件的驱动电流、跨导、衬底电流,HD 模型能捕捉速度过冲和能量输运,但需要仔细校准能量弛豫时间 τw 等参数。

经验法则:如果 DD 与 HD 的仿真结果差异小于 10%,且你关心的是电路级指标(如 Ion、Vth),优先用 DD 以节省时间;如果差异超过 20%,或你关心的是击穿电压、热载流子退化,则必须用 HD 或 TD 重新评估。

三、计算代价的量化对比与实操建议

以一个 2D MOSFET 结构(栅长 100 nm,网格约 2 万节点)为例,在相同网格和收敛精度下,各模型的典型计算代价如下表所示。表中数据为教学示例的合理估算,实际值随结构、网格和求解器设置变化。

模型 未知量数 典型迭代次数 相对计算时间 内存占用 收敛难度
DD 3 10~20 1× 低 容易
TD 4 20~40 3~5× 中 中等
HD 5~6 40~80 8~15× 大 较难
MC/BTE 分布函数 — 100×以上 极大 困难

在 Sentaurus 中切换模型主要通过 Physics 段实现。例如,启用 DD 模型时使用 Mobility(DopingDep HighFieldSaturation) 和 Recombination(SRH Auger);启用 HD 模型时增加 Hydrodynamic(eTemperature) 并指定 eTemperature 的初值策略;启用 TD 模型时增加 Thermodynamic 并设置热边界条件。以下是一个 HD 模型的 Physics 段示例:

Physics {
  Hydrodynamic(eTemperature)
  Mobility(
    DopingDep
    HighFieldSaturation
    Enormal
  )
  Recombination(
    SRH(DopingDep)
    Auger
    Avalanche(vanOverstraeten)
  )
  EffectiveIntrinsicDensity(OldSlotboom)
}

注意 HD 模型对初值非常敏感。建议先用 DD 模型求解到收敛,再以 DD 结果为初值切换到 HD,可显著提高收敛率。在 Solve 段中,可以先执行 Poisson 和 Coupled{Poisson Electron Hole},再执行 Coupled{Poisson Electron Hole eTemperature}。

常见误区:不要为了“更精确”而盲目启用 HD 或 MC。如果迁移率、复合、带隙等基础参数没有校准,HD 模型只会让错误的结果更复杂。模型层级的选择应与校准数据的可获得性匹配。

四、模型层级与校准策略的对应关系

不同模型层级需要校准的参数集不同。DD 模型主要校准低场迁移率、高场饱和速度、SRH 寿命和接触电阻;TD 模型额外需要校准热导率、热边界热阻;HD 模型还需要校准能量弛豫时间、能量相关迁移率参数。校准顺序应遵循“从低到高”:先用 DD 校准基础电学特性,再逐步引入 TD 或 HD 校准热效应和非平衡效应。

在实际项目中,一个高效的流程是:用 DD 完成器件结构、掺杂和基础 I-V 的快速迭代;确认电学行为合理后,切换到 TD 或 HD 进行关键物理点的验证;最后用 MC 对最关心的物理机制(如碰撞电离系数)做抽样验证。这样既控制了计算成本,又保证了关键结论的物理可信度。