从TCAD到SPICE:器件参数提取仿真与模型卡生成 - 第2节

2.1 Sentaurus TCAD 工具套件概览

Sentaurus TCAD 是 Synopsys 公司推出的业界标准半导体工艺与器件仿真平台。它覆盖从工艺步骤模拟、器件结构构建到电学特性提取的全链条,是连接物理制造与电路仿真的核心桥梁。本章聚焦其核心组件、网格生成策略及物理模型选择原则。

学习目标:掌握 Sentaurus 套件中 SProcess、SDevice、Sentaurus Workbench 等关键工具的功能边界;理解网格密度对仿真精度与收敛性的影响;学会根据器件类型(MOSFET、HEMT、IGBT)选择恰当的物理模型组合。

2.1.1 工具套件核心组件

Sentaurus 套件采用模块化架构,各工具通过统一的工程管理平台 Sentaurus Workbench (SWB) 进行集成与数据传递。下表列出主要组件及其职责:

工具名称 功能定位 输入/输出 典型应用场景
Sentaurus Process (SProcess) 工艺仿真:离子注入、扩散、氧化、刻蚀、沉积 输入:工艺步骤参数;输出:掺杂分布、应力场 CMOS 前段工艺开发、新型存储器工艺优化
Sentaurus Device (SDevice) 器件电学仿真:DC/AC/瞬态特性、击穿、自热 输入:网格+掺杂+物理模型;输出:I-V/C-V 曲线 MOSFET 阈值电压提取、功率器件雪崩击穿分析
Sentaurus Structure Editor (SDE) 几何建模与网格生成:2D/3D 结构构建、掺杂定义 输入:几何参数;输出:边界文件 (.bnd) + 网格文件 (.msh) FinFET 鳍片结构构建、IGBT 元胞设计
Sentaurus Visual (SVISUAL) 后处理可视化:电场/载流子分布、能带图、电流密度 输入:仿真结果 (.plt, .tdr);输出:图像/动画 热斑定位、沟道电场峰值分析
Sentaurus Workbench (SWB) 项目管理与参数扫描:多任务调度、DOE 设计 输入:工程文件 (.prj);输出:批量仿真结果 工艺参数敏感性分析、模型校准自动化
工程实践提示:在 TCAD 到 SPICE 的流程中,SDevice 是生成器件 I-V 数据的核心引擎,而 SProcess 负责提供准确的掺杂分布。若工艺参数未知,可直接在 SDE 中定义解析掺杂分布,跳过工艺仿真以加速模型卡开发。
Sentaurus TCAD 套件架构与数据流 Sentaurus Workbench (SWB) 统一工程管理 · 参数扫描 · DOE SProcess 工艺仿真 SDE 结构编辑 & 网格生成 SDevice 器件电学仿真 掺杂分布 网格 + 掺杂 输出仿真数据 I-V / C-V 曲线 .plt / .tdr 仿真结果 读取仿真结果 SVISUAL 后处理 电场 / 载流子分布可视化 参数提取 SPICE 模型卡生成

2.2 网格生成:精度与效率的平衡艺术

网格(Mesh)是 TCAD 仿真的空间离散化基础。网格密度直接影响求解精度和计算时间。Sentaurus 使用非结构化网格(三角形/四面体),支持局部加密以捕捉高梯度区域(如 PN 结界面、沟道反型层)。

2.2.1 网格类型与选择策略

网格类型 特点 适用场景 Sentaurus 实现
Box 网格 矩形单元,结构简单,生成快 MOSFET 沟道区、均匀掺杂区域 Grid: Box 命令
三角形网格 灵活适应复杂边界,支持局部加密 FinFET 鳍片、SOI 器件、功率器件终端 Grid: Triangle 命令
四面体网格 3D 结构标准网格,自由度最高 GAA 纳米片、3D NAND、IGBT 元胞 Grid: Tetrahedron 命令
混合网格 不同区域使用不同网格类型 大功率器件(粗网格体区+细网格结区) 区域定义 + Refine 命令
网格陷阱:过度加密会导致仿真时间呈指数增长(尤其 3D 瞬态仿真),而网格过粗则可能漏掉击穿路径或错误预测阈值电压。建议在关键结区至少保证 3~5 层网格,并在首次仿真后通过 SVISUAL 检查电场梯度,对高梯度区域进行二次加密。

2.2.2 网格生成流程(SDE 示例)

以下代码展示在 SDE 中为平面 MOSFET 定义网格的典型流程:

# 定义全局网格参数

Grid: Set::MaxEdgeLength = 0.5   # 最大边长 (μm)

Grid: Set::MinEdgeLength = 0.01  # 最小边长 (μm)



# 定义区域并指定网格密度

Region: "Channel" {

    Grid: Refine::MaxEdgeLength = 0.02

    Grid: Refine::MinEdgeLength = 0.005

}



Region: "SourceDrain" {

    Grid: Refine::MaxEdgeLength = 0.1

    Grid: Refine::MinEdgeLength = 0.02

}



# 在 PN 结界面自动加密

Grid: Refine::Interface {

    Regions = ["PBody", "Ndrift"]

    MaxEdgeLength = 0.01

    MinEdgeLength = 0.002

}



# 生成网格并输出

Grid: Generate

Grid: Write::Mesh = "mosfet_2d.msh"
网格质量指标:检查网格时关注 长宽比(Aspect Ratio) 和 内角(Internal Angle)。Sentaurus 推荐三角形网格内角 > 20°,长宽比 < 10:1,否则可能导致矩阵病态,收敛困难。
网格生成策略:全局粗网格 + 局部加密 全局粗网格 (MaxEdge=0.5μm) 体区 体区 局部加密 (MinEdge=0.01μm) 沟道区 PN 结区 Refine

2.3 物理模型选择:从漂移扩散到量子修正

物理模型是 TCAD 仿真的灵魂。Sentaurus 提供从经典到量子的多层次模型库,选择不当会导致仿真结果偏离物理实际或计算发散。

2.3.1 核心物理模型分类

模型类别 典型模型 适用器件 Sentaurus 关键字
载流子输运 漂移-扩散 (DD)、流体动力学 (HD)、蒙特卡洛 (MC) DD:长沟道 MOSFET;HD:短沟道;MC:亚 10nm Physics: Hydrodynamic
迁移率退化 Lombardi 表面散射、Caughey-Thomas 高场饱和、掺杂依赖 所有 MOSFET、HEMT Mobility: Lombardi
产生-复合 Shockley-Read-Hall (SRH)、Auger、碰撞电离 (Impact Ionization) SRH:漏电分析;Auger:功率器件;碰撞电离:击穿 Recombination: SRH
量子效应 密度梯度 (DG)、薛定谔泊松耦合 (Schrodinger-Poisson) DG:FinFET 阈值电压修正;SP:量子阱器件 Quantum: DensityGradient
热效应 热力学模型 (Thermodynamic)、晶格加热 SOI 自热、功率 HEMT 热管理 Physics: Thermodynamic

2.3.2 模型选择决策流程

以下决策树帮助工程师根据器件特征快速锁定模型组合:

  1. 判断沟道长度:Lg > 1μm → 漂移-扩散模型;Lg < 100nm → 必须启用流体动力学或密度梯度。
  2. 判断工作电压:Vds > 击穿电压的 80% → 必须开启碰撞电离模型。
  3. 判断衬底类型:SOI 或超薄体 → 开启热力学模型以捕捉自热效应。
  4. 判断材料体系:SiGe/SiC/GaN → 使用专用迁移率模型(如 Caughey-Thomas 参数化)。
  5. 验证收敛性:若仿真发散,先关闭非必要模型(如量子修正),从纯 DD 模型开始逐步添加。
模型校准黄金法则:在开始 SPICE 模型卡提取前,务必使用至少一组实测 I-V 数据对 TCAD 物理模型参数(如迁移率系数、SRH 寿命)进行校准。未校准的 TCAD 仿真输出不能直接用于模型卡生成,否则误差会传递到电路级。
物理模型选择决策流程 开始:器件结构定义 Lg > 1μm ? 是 否 漂移-扩散 (DD) 流体动力学 (HD) Vds > 80% BV ? Lg < 100nm ? 是 是 否 否 开启碰撞电离 开启密度梯度 模型组合确定 补充判据:SOI / 超薄体 → 热力学模型;SiGe / SiC / GaN → 专用迁移率模型 收敛性验证:若发散,先关闭量子修正等非必要模型,从纯 DD 起步逐步添加

2.4 本章小结

  • Sentaurus 套件以 SWB 为集成平台,SProcess/SDE/SDevice 分别负责工艺、结构、电学仿真,SVISUAL 用于结果分析。
  • 网格生成遵循“全局粗网格 + 局部加密”原则,关键结区和沟道区最小边长建议 ≤ 0.01μm,并检查长宽比和内角。
  • 物理模型选择需根据沟道长度、工作电压、材料体系分层决策,从漂移-扩散基础模型逐步添加流体动力学、碰撞电离、量子修正等高级模型。
  • 模型校准是 TCAD 到 SPICE 流程的关键质量门,未经实测数据校准的仿真结果不能用于模型卡生成。
下一步实践:在 Sentaurus Workbench 中创建一个简单的 2D MOSFET 工程,分别使用 DD 和 HD 模型仿真转移特性曲线,观察迁移率退化模型对 Ion 和 Ioff 的影响。这将为后续第 3 章“器件参数提取”奠定网格和模型基础。