4. TCAD 仿真环境搭建：工具链、项目结构与流程概览
本章是进入实际仿真的第一步。我们将系统梳理 GaN 器件工艺角仿真所需的 TCAD 工具链（以 Sentaurus 和 Silvaco 为例）、标准项目文件结构以及从结构定义到数据提取的完整仿真流程。掌握这些基础，后续的工艺角扫描与 C-V 曲线拟合才能高效落地。
核心目标：理解 TCAD 仿真不是“一键运行”，而是由多个工具模块、输入文件与脚本协同完成的工程过程。本章结束后，你应能独立搭建一个可复用的 GaN HEMT 仿真项目文件夹。
4.1 TCAD 工具链概览
主流商用 TCAD 工具包括 Synopsys Sentaurus 与 Silvaco Victory/Atlas。两者在物理模型精度与网格算法上各有侧重，但核心工作流高度一致：结构定义 → 网格生成 → 物理求解 → 数据提取。下表对比了两套工具链的关键模块与对应功能。
| 功能阶段 | Sentaurus 工具模块 | Silvaco 工具模块 | 说明 |
|---|---|---|---|
| 结构/工艺定义 | Sentaurus Structure Editor (SDE) | Victory Process / Atlas Syntax | 定义器件几何、掺杂、材料层 |
| 网格生成 | Sentaurus Mesh (SnMesh) | Victory Mesh / Atlas Mesh | 生成非结构化/结构化网格 |
| 物理求解 | Sentaurus Device (SDevice) | Atlas / Victory Device | 求解 Poisson、Drift-Diffusion 等方程 |
| 数据可视化 | Sentaurus Visual (SVISUAL) | TonyPlot / DeckBuild | 绘制 C-V、I-V 曲线与电场分布 |
| 参数提取 | Inspect / 自定义 Python 脚本 | TonyPlot Extract / Python | 提取阈值电压、电容值等关键参数 |
选型建议：本课程以 Sentaurus 为例进行讲解，但所有流程概念与脚本逻辑均可直接迁移至 Silvaco。如果你所在课题组使用 Silvaco，只需将 SDE 命令替换为 Atlas 语法，SDevice 替换为 Atlas 求解语句即可。
4.2 项目文件结构规范
一个规范的 TCAD 项目文件夹能大幅减少调试时间。以下是针对 GaN 工艺角仿真推荐的目录结构，每个文件夹承担明确职责。
GaN_HEMT_Project/
├── 00_tech/ # 工艺技术文件（.par, .lib）
│ ├── GaN_material.par # GaN/AlGaN 材料参数
│ └── polarization.par # 极化电荷模型参数
├── 01_sde/ # 结构定义脚本
│ ├── gan_hemt_sde.cmd # SDE 命令流
│ └── gan_hemt_sde_dvs.cmd # 掺杂/电极定义
├── 02_mesh/ # 网格生成脚本
│ ├── gan_hemt_mesh.cmd # SnMesh 命令
│ └── gan_hemt_mesh.grd # 生成的网格文件
├── 03_sdevice/ # 物理求解脚本
│ ├── CV_sweep.cmd # C-V 扫描求解命令
│ ├── CV_sweep_plot.cmd # 输出控制
│ └── output/ # 仿真结果（.plt, .tdr）
├── 04_analysis/ # 后处理与拟合
│ ├── extract_cv.py # Python 提取 C-V 数据
│ └── cv_fitting.py # 曲线拟合脚本
├── 05_corner/ # 工艺角批量脚本
│ ├── run_corner.sh # 批量运行脚本
│ └── corner_table.csv # 工艺角参数表
└── README.md # 项目说明
常见错误：初学者常将所有文件放在同一目录下，导致版本混乱。请严格按上述结构组织，尤其是将
00_tech/ 中的材料参数文件与求解脚本分离，便于后续更换材料模型。
4.3 仿真流程概览
一个完整的 GaN 器件 C-V 仿真流程包含 5 个核心步骤。下面的 SVG 流程图直观展示了各步骤之间的数据传递关系。
4.4 各步骤详细说明
4.4.1 结构定义（SDE）
使用 Sentaurus Structure Editor 绘制 GaN HEMT 的二维截面。关键操作包括：
- 材料层堆叠：定义衬底（SiC/Si）、GaN 缓冲层、AlGaN 势垒层、GaN 帽层。
- 极化电荷设置：在 AlGaN/GaN 界面设置固定极化电荷（正值 1e13 cm⁻² 量级）。
- 电极定义：指定源极、漏极、栅极的金属接触区域与功函数。
4.4.2 网格生成（SnMesh）
网格质量直接影响仿真收敛性与精度。针对 GaN 器件，需在以下区域加密：
- AlGaN/GaN 异质结界面（2DEG 沟道区域）
- 栅极边缘（电场集中区域）
- 源/漏接触边缘
典型网格尺寸：沟道垂直方向 0.5 nm，水平方向 10 nm。
4.4.3 物理求解（SDevice）
在 SDevice 命令文件中，需指定：
- 物理模型：Poisson 方程、Drift-Diffusion 传输、高场迁移率退化、极化模型。
- 偏置条件：C-V 扫描时，栅极电压从 -10 V 扫描至 +2 V，步长 0.1 V；源漏短接（0 V）。
- 输出控制：记录每个偏置点下的栅极总电荷，用于后续电容计算。
4.4.4 数据提取与拟合
仿真完成后，使用 Python 脚本读取 .plt 文件中的电荷-电压数据，通过数值微分计算电容：
import numpy as np
import matplotlib.pyplot as plt
# 读取仿真数据（示例格式）
data = np.loadtxt('CV_sweep.plt', skiprows=1)
Vg = data[:, 0] # 栅极电压
Qg = data[:, 1] # 栅极总电荷
# 计算微分电容 C = dQ/dV
C = np.gradient(Qg, Vg)
# 绘制 C-V 曲线
plt.plot(Vg, C * 1e15, 'b-', linewidth=2) # 单位转换为 fF
plt.xlabel('Gate Voltage (V)')
plt.ylabel('Capacitance (fF)')
plt.grid(True)
plt.show()
效率提升：在批量工艺角仿真中，建议将上述 Python 脚本封装为函数，并配合
run_corner.sh 循环调用 SDevice，实现全自动扫描与拟合。
4.5 本章小结
本章建立了 TCAD 仿真环境的完整认知框架：
- 工具链：Sentaurus 与 Silvaco 在模块功能上高度对应，选择取决于课题组现有许可。
- 项目结构：按
tech → sde → mesh → sdevice → analysis分层组织，确保可复现性。 - 流程闭环：从结构定义到拟合分析形成完整数据流，反馈回路支持工艺角参数迭代优化。