从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 中定义解析掺杂分布,跳过工艺仿真以加速模型卡开发。
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,否则可能导致矩阵病态,收敛困难。
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 模型选择决策流程
以下决策树帮助工程师根据器件特征快速锁定模型组合:
- 判断沟道长度:Lg > 1μm → 漂移-扩散模型;Lg < 100nm → 必须启用流体动力学或密度梯度。
- 判断工作电压:Vds > 击穿电压的 80% → 必须开启碰撞电离模型。
- 判断衬底类型:SOI 或超薄体 → 开启热力学模型以捕捉自热效应。
- 判断材料体系:SiGe/SiC/GaN → 使用专用迁移率模型(如 Caughey-Thomas 参数化)。
- 验证收敛性:若仿真发散,先关闭非必要模型(如量子修正),从纯 DD 模型开始逐步添加。
模型校准黄金法则:在开始 SPICE 模型卡提取前,务必使用至少一组实测 I-V 数据对 TCAD 物理模型参数(如迁移率系数、SRH 寿命)进行校准。未校准的 TCAD 仿真输出不能直接用于模型卡生成,否则误差会传递到电路级。
2.4 本章小结
- Sentaurus 套件以 SWB 为集成平台,SProcess/SDE/SDevice 分别负责工艺、结构、电学仿真,SVISUAL 用于结果分析。
- 网格生成遵循“全局粗网格 + 局部加密”原则,关键结区和沟道区最小边长建议 ≤ 0.01μm,并检查长宽比和内角。
- 物理模型选择需根据沟道长度、工作电压、材料体系分层决策,从漂移-扩散基础模型逐步添加流体动力学、碰撞电离、量子修正等高级模型。
- 模型校准是 TCAD 到 SPICE 流程的关键质量门,未经实测数据校准的仿真结果不能用于模型卡生成。
下一步实践:在 Sentaurus Workbench 中创建一个简单的 2D MOSFET 工程,分别使用 DD 和 HD 模型仿真转移特性曲线,观察迁移率退化模型对 Ion 和 Ioff 的影响。这将为后续第 3 章“器件参数提取”奠定网格和模型基础。