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 仿真的常用选择。
蒙特卡洛方法直接求解玻尔兹曼输运方程,不依赖迁移率模型,物理最精确,但计算代价极高,通常只用于小尺寸器件或关键物理验证。
二、各模型的适用场景与选择依据
选择模型时,先问三个问题:器件特征尺寸是否接近或小于载流子能量弛豫长度(约 50~100 nm)?是否存在显著自加热或高电流密度?是否需要描述碰撞电离、雪崩击穿等非平衡效应?
对于常规 180 nm 及以上 CMOS 工艺的 I-V 特性、阈值电压、亚阈值摆幅分析,DD 模型配合 Lombardi 迁移率模型和 SRH 复合模型通常足够,且收敛快、参数少。对于功率 LDMOS 或 IGBT 的导通电阻与热分布分析,TD 模型能给出晶格温度分布,帮助判断热失效点。对于 40 nm 以下器件的驱动电流、跨导、衬底电流,HD 模型能捕捉速度过冲和能量输运,但需要仔细校准能量弛豫时间 τw 等参数。
三、计算代价的量化对比与实操建议
以一个 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}。
四、模型层级与校准策略的对应关系
不同模型层级需要校准的参数集不同。DD 模型主要校准低场迁移率、高场饱和速度、SRH 寿命和接触电阻;TD 模型额外需要校准热导率、热边界热阻;HD 模型还需要校准能量弛豫时间、能量相关迁移率参数。校准顺序应遵循“从低到高”:先用 DD 校准基础电学特性,再逐步引入 TD 或 HD 校准热效应和非平衡效应。
在实际项目中,一个高效的流程是:用 DD 完成器件结构、掺杂和基础 I-V 的快速迭代;确认电学行为合理后,切换到 TD 或 HD 进行关键物理点的验证;最后用 MC 对最关心的物理机制(如碰撞电离系数)做抽样验证。这样既控制了计算成本,又保证了关键结论的物理可信度。