2. 器件仿真基础：TCAD仿真工具链介绍
TCAD（Technology Computer-Aided Design）是半导体工艺与器件仿真的核心手段。本章系统介绍主流TCAD工具链的组成、仿真标准流程以及各环节的关键技术要点，为后续工艺角优化与器件性能分析奠定基础。
2.1 主流TCAD工具链对比
当前工业界与学术界广泛使用的TCAD工具主要包括Synopsys Sentaurus、Silvaco TCAD以及开源工具QucsStudio。三者在架构设计、物理模型库、求解器性能及用户接口上各有侧重。
| 特性维度 | Sentaurus | Silvaco | QucsStudio |
|---|---|---|---|
| 开发商 | Synopsys | Silvaco Inc. | 开源社区（Qucs衍生） |
| 核心架构 | 模块化工具链（SDE/SDevice/SVisual） | 集成化平台（DeckBuild + Atlas） | 单界面集成（原理图+仿真+后处理） |
| 物理模型库 | 最全面（量子效应、应力、热效应） | 丰富（迁移率、复合、隧穿模型） | 基础（漂移-扩散、SRH复合） |
| 求解器 | 耦合/非耦合Newton迭代 | Gummel/Newton混合 | 直接稀疏矩阵求解 |
| 网格引擎 | Noffset3D（各向异性自适应） | Mesh（张量积/三角形） | 内置网格生成器 |
| 典型应用 | 先进节点（7nm以下）器件优化 | 功率器件、存储器、光电器件 | 教学演示、基础器件验证 |
| 许可证 | 商业授权（高成本） | 商业授权（中等成本） | GPL开源（免费） |
2.2 标准仿真流程详解
无论使用哪种TCAD工具，器件仿真均遵循统一的五阶段流程。每个阶段的选择直接影响仿真精度与计算效率。
- 结构定义：描述器件的几何尺寸、掺杂分布、电极位置与材料属性。
- 网格划分：将连续器件区域离散化为有限元或有限体积网格节点。
- 物理模型：选择载流子迁移率、复合、隧穿、量子效应等物理方程。
- 求解器：采用Newton、Gummel或耦合迭代算法求解非线性偏微分方程组。
- 后处理：提取I-V曲线、C-V曲线、电场分布、载流子浓度分布等结果。
2.3 结构定义与网格划分
结构定义是仿真的起点，直接影响后续所有环节的准确性。以Sentaurus SDE为例，典型的结构定义包含以下要素：
- 几何区域：定义衬底、阱区、栅极、源漏区的矩形或多边形边界。
- 掺杂分布：指定高斯、误差函数或均匀掺杂的峰值浓度与结深。
- 电极定义：标记源极、漏极、栅极、衬底接触的边界段。
- 材料属性：设置硅、二氧化硅、氮化硅、高k介质等材料的介电常数与禁带宽度。
网格划分是数值仿真的核心瓶颈。过粗的网格导致精度不足，过细的网格造成计算爆炸。工程实践中采用自适应网格策略：
| 网格类型 | 适用场景 | 优点 | 缺点 |
|---|---|---|---|
| 均匀网格 | 简单结构、教学演示 | 实现简单，无畸变 | 无法局部加密，效率低 |
| 张量积网格 | MOSFET、二极管 | 正交性好，适合漂移-扩散 | 对弯曲边界适应性差 |
| 三角形/四面体网格 | 复杂几何（FinFET、SOI） | 几何适应性强 | 数值耗散较大 |
| 自适应网格 | 高精度工程仿真 | 按误差指示自动加密/粗化 | 迭代次数多，设置复杂 |
2.4 物理模型与求解器配置
物理模型的选择决定了仿真能否捕捉真实的器件行为。以MOSFET仿真为例，必须激活的核心模型包括：
- 迁移率模型：Lombardi表面粗糙度散射、Caughey-Thomas高场饱和、掺杂依赖迁移率。
- 复合模型：Shockley-Read-Hall（SRH）复合、Auger复合、辐射复合。
- 量子修正模型：密度梯度（Density Gradient）或van Dort模型，用于沟道反型层量子化。
- 隧穿模型：Band-to-Band隧穿（BTBT）、Fowler-Nordheim隧穿（栅漏电）。
求解器配置涉及迭代算法选择与收敛控制。以下为Sentaurus SDevice的典型求解器设置片段：
# Sentaurus SDevice 求解器配置示例
Solve {
# 初始解：零偏置泊松方程
Poisson
# 耦合求解：Newton迭代
Coupled { Poisson Electron Hole }
# 电压扫描：漏极电压从0V到1.2V，步长0.05V
Quasistationary (
InitialStep = 0.01
Increment = 1.5
MinStep = 1e-5
MaxStep = 0.05
Goal { Name = "drain" Voltage = 1.2 }
)
{ Coupled { Poisson Electron Hole } }
}2.5 后处理与结果提取
后处理是将原始求解数据转化为工程可读信息的关键步骤。常用提取内容包括：
| 提取类型 | 典型曲线/分布 | 工程意义 |
|---|---|---|
| 电学特性 | ID-VG、ID-VD、IG-VG | 阈值电压、驱动电流、亚阈值摆幅、漏电 |
| 电容特性 | C-V曲线（CGG、CGD、CGS） | 栅电容、米勒电容、开关速度 |
| 内部物理量 | 电势分布、载流子浓度、电场强度 | 击穿位置、热载流子效应、短沟道效应 |
| 可靠性指标 | 碰撞电离率、栅漏电流密度 | 热载流子注入（HCI）、负偏置温度不稳定性（NBTI） |
2.6 工程案例：MOSFET仿真快速入门
以下通过一个完整的MOSFET仿真案例，串联上述所有流程。使用QucsStudio平台（开源，适合教学），实现从结构定义到I-V曲线提取的全过程。
步骤1：结构定义 — 在QucsStudio原理图中放置NMOS器件，设置沟道长度L=0.18μm，氧化层厚度Tox=4nm，衬底掺杂NA=1×1017 cm-3。
步骤2：网格划分 — QucsStudio自动生成非均匀网格，在沟道与氧化层界面处加密（最小网格尺寸0.5nm）。
步骤3：物理模型 — 激活漂移-扩散传输、SRH复合、Lombardi迁移率模型。
步骤4：求解器 — 采用Newton耦合迭代，设置漏极电压扫描范围0~1.8V，步长0.1V。
步骤5：后处理 — 提取ID-VG曲线，计算阈值电压Vth=0.45V，亚阈值摆幅SS=72 mV/dec。
// QucsStudio 网表示例（NMOS直流扫描）
.netlist
VGS 1 0 DC 0
VDS 2 0 DC 0.1
M1 2 1 0 0 NMOS L=0.18u W=10u
.model NMOS NMOS (VTO=0.4 TOX=4e-9 NSUB=1e17 U0=350)
.dc VGS 0 1.2 0.05 VDS 0.05 1.8 0.1
.plot i(VDS)
.end