4. TCAD仿真基础

在SOI器件的自热效应（SHE）仿真中，TCAD工具链的选择与配置直接决定了仿真结果的物理准确性与收敛效率。本章聚焦于业界两大主流TCAD平台——Synopsys Sentaurus与Silvaco Victory/Atlas，系统介绍其工具链架构、热仿真与电学仿真耦合的物理模型配置，以及从器件结构定义到热-电耦合求解的完整工作流。

4.1 主流TCAD工具链架构对比

Sentaurus与Silvaco均采用模块化工具链设计，但各自在网格引擎、求解器策略与热模型接口上存在显著差异。下表从SOI自热仿真视角对比两大平台的核心组件：

4.2 热仿真与电学仿真耦合的物理基础

自热效应仿真必须同时求解泊松方程、载流子连续性方程与晶格热方程。在Sentaurus中，热-电耦合通过以下控制方程实现：

// 晶格热方程 (Sentaurus Thermodynamic Model)
∇·(κ ∇T) + H = ρ·c·∂T/∂t

// 其中热产生项 H 包含：
// H = J·E + (Eg + 3kT)·(R - G) + 载流子能量弛豫项

// 迁移率退化模型 (考虑晶格温度与自热)
μ_eff(T) = μ_0 · (T / T0)^(-α)
// α 典型值：电子 1.5，空穴 1.2 (SOI 薄体)

Silvaco Atlas中对应的热模型开启语法如下：

# Silvaco Atlas 热仿真配置示例
models region=1 lat.temp print
# 开启晶格热模型

material region=1 \
  taurel.el=0.3e-12 \
  taurel.ho=0.3e-12 \
  lat.temp

# 指定能量弛豫时间（影响自热升温幅度）

solve init
method newton trap maxtrap=10
solve vdrain=0.0 vgate=1.0 vstep=0.1 vfinal=1.5 \
  name=drain electtemp lat.temp
# 同时求解电学与晶格温度

4.3 工具链工作流：从结构到热-电耦合求解

下图展示了基于Sentaurus的SOI器件自热效应仿真标准流程，包含从器件结构定义到热-电耦合后处理的关键步骤：

4.4 热-电耦合的关键配置参数

在Sentaurus Device中，热-电耦合的精度高度依赖于以下参数的合理设置。下表总结了SOI自热仿真中必须校准的核心参数：

4.5 迁移率退化模型在热仿真中的集成

自热效应导致晶格温度升高，进而通过多种散射机制降低载流子迁移率。在TCAD中，迁移率退化模型必须与热方程耦合。以下为Sentaurus中常用的迁移率模型配置：

// Sentaurus Device 迁移率模型 (热耦合)
Physics {
  Mobility (
    DopingDependence
    HighFieldSaturation
    Enormal   // 考虑垂直电场 (SOI 薄体)
    TemperatureDependence  // 开启温度依赖
  )
  EffectiveIntrinsicDensity( BandGapNarrowing( OldSlotboom ) )
  Thermodynamic  // 开启热力学模型 (关键)
}

// 自定义迁移率温度退化系数 (覆盖默认值)
Material = "Silicon" {
  Mobility {
    mu0_TAU = 1.5   // 电子温度指数
    mu0_TAU_H = 1.2 // 空穴温度指数
  }
}

Silvaco Atlas中对应的迁移率-热耦合配置：

# Silvaco Atlas 迁移率热退化
models region=1 \
  conmob \
  fldmob \
  srh \
  lat.temp \
  print

# 指定温度依赖迁移率参数
material region=1 \
  mun.tau=1.5 \
  mup.tau=1.2 \
  taurel.el=0.3e-12 \
  taurel.ho=0.3e-12

# 输出晶格温度与迁移率
output region=1 \
  lat.temp \
  mobility

4.6 收敛性优化与调试策略

热-电耦合仿真常因强非线性而出现收敛困难。以下为经过验证的收敛性优化步骤：

1. 初始条件预热：先执行纯电学仿真 (不开启热方程) 获得初始电势分布，再开启热耦合。
2. 逐步加载偏置：使用较小的电压步长 (如 Vd_step=0.05V) 并配合 ramp 策略。
3. Newton 阻尼：设置阻尼因子 damp=1e-3 防止温度突变。
4. 网格优化：在 BOX 界面与沟道-栅氧界面进行网格加密 (最小网格 1nm)。
5. 时间步长控制：瞬态仿真时，初始时间步长设为 1e-12 s，最大步长不超过 1e-9 s。

4.7 本章小结

本章系统介绍了Sentaurus与Silvaco两大TCAD平台在SOI器件自热效应仿真中的工具链架构、热-电耦合物理模型、关键配置参数与收敛性优化策略。核心要点包括：

• 热-电耦合必须同时求解泊松方程、载流子连续性方程与晶格热方程，三者通过热产生项与温度依赖迁移率形成反馈回路。
• Sentaurus的Thermodynamic模型与Silvaco的Lattice Heat模型在接口语法上存在差异，但物理本质一致。
• 迁移率温度退化指数 α、BOX 热导率、能量弛豫时间 τ_e 是影响自热仿真精度的三个最敏感参数。
• 合理的网格加密、Newton 阻尼与逐步加载策略是保障热-电耦合收敛的关键工程实践。