4. Sentaurus Device物理模型体系

在Sentaurus Device中，载流子输运模型的选择直接决定了仿真能捕捉到哪些物理效应，也决定了计算代价与收敛难度。功率半导体器件的核心特征——高电压、大电流、强电场、显著自热——使得模型选择不再是"精度越高越好"的简单问题，而是要在物理保真度与数值可行性之间做出工程判断。本节围绕漂移扩散（Drift-Diffusion, DD）、热力学（Thermodynamic, TD）、流体动力学（Hydrodynamic, HD）三类模型，说明其物理内涵、适用边界与典型功率器件场景。

三类模型的物理内涵差异

漂移扩散模型是最基础的输运框架。它假设载流子与晶格处于热平衡，载流子温度等于晶格温度，电流由漂移项与扩散项构成：

Jn = q·μn·n·E + q·Dn·∇n
Jp = q·μp·p·E − q·Dp·∇p

其中迁移率μ和扩散系数D通过爱因斯坦关系关联。DD模型求解泊松方程与载流子连续性方程，未知量为电势ψ、电子浓度n、空穴浓度p。它的优势是变量少、雅可比矩阵规模小、收敛稳健，适合绝大多数中低压器件的静态与瞬态分析。

热力学模型在DD基础上引入晶格热流方程，把晶格温度T_L作为第四个未知量。它仍然假设载流子与晶格瞬时热平衡，但允许晶格温度在空间上分布不均。对于功率器件，这意味着可以计算自热引起的温升、温度梯度以及温度对迁移率和本征载流子浓度的影响。TD模型是功率器件电热耦合仿真的默认起点。

流体动力学模型进一步把载流子温度T_n、T_p从晶格温度中分离出来，为电子和空穴分别增加能量平衡方程。它能够描述载流子被强电场加热后远离热平衡的状态，即"热载流子"效应。在功率器件中，雪崩击穿、栅氧注入、热载流子退化等场景下，载流子温度可显著高于晶格温度，此时DD和TD都会低估碰撞电离率与注入电流。HD模型的代价是未知量增至六个以上，方程强耦合、非线性强，收敛对初值和网格质量极为敏感。

选择依据：从物理效应反推模型

模型选择的第一步不是看器件类型，而是看目标物理效应是否被模型覆盖。如果仿真目标只涉及导通电阻、阈值电压、体二极管正向压降，且功率密度不高、自热可忽略，DD模型足够。若器件工作电流密度超过约1×10³ A/cm²，或需要评估稳态结温、热阻网络，则必须启用TD模型。若目标涉及雪崩击穿电压的精确值、栅氧注入电流或开关瞬态中的过冲电流，则需要HD模型。

一个实用的判断流程是：先问"晶格温度是否显著变化"，若是则排除纯DD；再问"载流子温度是否显著高于晶格温度"，若是则必须用HD。对于硅器件，当峰值电场超过约3×10⁵ V/cm时，载流子加热效应开始显著；对于碳化硅器件，由于临界击穿场更高，这一阈值可放宽到约1×10⁶ V/cm以上，但高场下的热载流子效应同样不可忽略。

功率器件典型场景的模型匹配

以1200V SiC MOSFET的雪崩耐量仿真为例。若使用DD模型，碰撞电离率由局部电场通过Chynoweth公式直接计算，忽略了载流子温度对电离系数的影响，得到的击穿电压通常偏高，且无法正确描述雪崩过程中的电流丝形成。切换到HD模型后，电子温度在耗尽区边缘显著升高，电离系数增大，击穿电压下降，更接近实验值。但HD模型对网格要求极高，需要在结附近加密到纳米级，且时间步长需控制在皮秒量级，否则能量方程不收敛。

对于IGBT的短路耐受仿真，TD模型通常是首选。短路时电流密度可达数百A/cm²，晶格温度在微秒内上升数百开尔文，自热效应主导器件行为。此时载流子温度虽也有升高，但晶格温度的变化对迁移率和本征载流子浓度的影响更为关键。TD模型可以在可接受的收敛性下给出结温随时间的变化曲线，用于评估短路耐受时间。

参数与数值设置的差异

三类模型在参数配置上的差异主要体现在迁移率模型、复合模型和能量弛豫时间上。DD模型通常使用DopingDependence和HighFieldSaturation迁移率模型；TD模型需额外启用晶格热流方程，并设置热导率、热容和热边界条件；HD模型则必须指定能量弛豫时间τ_w，该参数对结果影响显著，硅中电子典型值为0.1~0.5 ps，空穴为0.2~0.8 ps，需根据文献或校准确定。

模型选择的决策要点

综合来看，模型选择应遵循"效应驱动、分级验证"的原则。首先明确仿真要回答的工程问题，列出必须捕捉的物理效应；然后对照三类模型的覆盖范围，选择能够覆盖这些效应的最低阶模型；最后通过网格收敛性检查和实验数据对比验证模型适用性。对于功率器件，TD模型是电热耦合分析的基线，HD模型则用于击穿、注入和瞬态过冲等局部强场场景。盲目使用HD模型不仅浪费计算资源，还可能因收敛失败而得不到任何结果。