3. 工艺仿真中的网格与结构表示

在Sentaurus工艺仿真（Sentaurus Process）中，网格不只是数值计算的“背景”，而是决定仿真能否收敛、结果是否可信的核心数据结构。工艺仿真与器件仿真最大的不同在于：氧化、扩散、注入、刻蚀、沉积等步骤会不断改变材料边界和掺杂分布，因此网格必须随结构演化而动态调整。理解一维、二维、三维网格的生成策略、自适应细化机制以及结构文件的组织方式，是搭建可靠工艺仿真流程的基础。

一、网格的维度与生成策略

Sentaurus Process支持一维、二维和三维结构表示。一维网格通常用于快速标定扩散模型或注入损伤模型，计算量小、迭代快，适合在正式二维/三维仿真前验证物理参数。二维网格是工艺仿真最常用的形式，能够描述MOSFET、二极管等器件的截面结构演化。三维网格则用于FinFET、纳米片、GAA等具有强三维几何效应的结构，但网格量和计算成本显著上升。

网格生成的核心是“在关键区域加密，在无关区域稀疏”。以二维MOS结构为例，沟道、源漏结、栅氧化层附近掺杂梯度大，需要细网格；衬底深处掺杂变化平缓，可用粗网格。Sentaurus Process通过grid命令或mgoals（网格目标）控制局部细化。例如：

grid set.min.normal.size=1e-3
grid set.max.normal.size=0.05
mgoals on location= [0 0] spacing= 0.002
mgoals on location= [0 0.05] spacing= 0.001
mgoals on location= [0 0.2] spacing= 0.01

上述示例中，set.min.normal.size和set.max.normal.size限定网格法向尺寸范围，mgoals在指定位置施加局部间距目标。实际工程中，网格间距应小于特征长度的1/3～1/5，例如栅氧化层厚度为2 nm时，氧化层内至少布置3～5层网格，间距约0.4～0.7 nm。

二、自适应网格细化

工艺仿真中，掺杂浓度、应力、氧化剂浓度等物理量在界面附近变化剧烈。自适应网格细化（Adaptive Mesh Refinement, AMR）根据解的梯度或误差估计自动加密网格，在保证精度的同时控制网格规模。Sentaurus Process的AMR通常由refinebox和adaptive相关参数触发。

常见的触发判据包括：掺杂浓度梯度超过阈值、氧化层厚度变化率过大、界面曲率过高、应力集中等。例如，在注入后的退火步骤中，结深附近掺杂浓度从1e20 cm⁻³骤降到1e15 cm⁻³，若网格过粗，结深位置会产生明显数值误差。此时可设置：

refinebox name= "junction" min= [0 0] max= [0.2 0.3]
refinebox add= "junction" spacing= 0.002
adaptive on max.level= 4 min.level= 1

其中max.level限制最大细化层级，避免网格无限加密导致内存溢出。实际使用中，建议先以较粗网格试算，观察收敛性和关键物理量（如结深、氧化层厚度）的变化，再逐步加密，直到结果变化小于设定容差（例如结深变化小于1%）。

三、结构文件格式与可视化

Sentaurus Process生成的结构文件主要有两种格式：.tif（TDR格式）和.grd（旧版网格格式）。TDR（Taurus Data Repository）是Sentaurus系列工具的标准结构格式，包含网格节点、单元连接、材料区域、掺杂分布、应力场等信息。TDR文件可由Sentaurus Visual、Sentaurus Workbench或Tecplot SV等工具打开。

结构文件的核心内容包括：节点坐标、单元拓扑（三角形/四面体/六面体）、材料标识、掺杂数据、缺陷数据、应力张量等。以二维TDR为例，其内部按“区域—材料—数据字段”组织，每个数据字段对应一个物理量。可视化时，常用切面图、等值线图、矢量图展示掺杂浓度、电场、应力分布。

在Sentaurus Visual中，可通过svisual命令加载TDR文件并绘制掺杂等值线：

svisual -b -p "load.tdr; plot.contour Boron; export.png"

实际工程中，建议在每步关键工艺后保存TDR文件，便于回溯和对比。例如，注入后、退火后、氧化后分别保存struct_after_implant.tdr、struct_after_anneal.tdr、struct_after_oxidation.tdr，通过对比结深、氧化层厚度和掺杂分布，判断工艺步骤是否符合预期。

四、网格质量与仿真收敛的关系

网格质量直接影响求解器的收敛性。常见问题包括：网格畸变（单元内角过大或过小）、长宽比过高、界面处网格不连续等。在Sentaurus Process中，可通过grid check命令检查网格质量，输出最小角、最大长宽比等指标。一般要求三角形单元最小角大于20°，长宽比小于10。

当仿真出现“求解不收敛”或“掺杂分布异常”时，应优先检查网格：关键区域是否足够细、界面处是否对齐、自适应细化是否触发。必要时可手动插入refinebox或调整mgoals间距。对于三维结构，还需注意网格数量与内存的平衡，通常三维网格节点数控制在50万～200万以内较为可行。

综上，工艺仿真中的网格与结构表示是一个“生成—细化—存储—可视化—验证”的闭环。一维/二维/三维网格的选择取决于器件结构和精度需求；自适应细化是控制精度与效率平衡的关键手段；TDR结构文件承载了完整的物理场信息，是后续器件仿真和结果分析的桥梁。实际工程中，建议在每步关键工艺后保存结构文件，并通过网格质量检查和收敛性验证，确保仿真结果可信。