4. MOSFET工艺仿真中的网格策略

在Sentaurus工艺仿真(Sentaurus Process, SProcess)中,网格不是“越密越好”的附属设置,而是决定仿真能否收敛、结果是否可信、计算成本是否可接受的核心数值参数。工艺仿真求解的是扩散-反应-应力耦合的偏微分方程组,网格节点承载掺杂浓度、缺陷浓度、应力张量等物理量,网格的空间分布直接决定了离散误差的大小和线性系统的条件数。本节围绕一维/二维/三维网格生成、自适应网格细化、关键区域网格密度控制,以及网格质量对收敛与精度的影响展开。

一、网格的维度选择与生成方式

Sentaurus Process的网格维度由仿真任务决定。一维网格用于快速标定扩散模型参数,例如仅需提取结深和薄层电阻时,一维仿真可在秒级完成;二维网格是MOSFET工艺仿真的主力,用于刻画源漏结、LDD浅结、侧墙下的横向扩散和氧化层台阶;三维网格仅在需要分析窄沟道效应、STI角效应或三维应力分布时启用,其节点数通常是二维的10²~10³倍。

网格生成有两种基本途径:结构化网格和自适应网格。结构化网格由line、region等命令显式定义,适合几何规则、物理量梯度已知的区域;自适应网格则由求解器根据解的变化自动加密或粗化。实际工程中通常采用“结构化打底 + 自适应细化”的混合策略:先用较粗的结构化网格保证几何拓扑正确,再让自适应机制在掺杂梯度大、缺陷团簇活跃的区域自动加密。

MOSFET工艺仿真网格策略:维度—细化—密度—质量闭环 维度选择 1D:模型标定、结深/方块电阻 2D:源漏结、LDD、侧墙、氧化台阶 3D:窄沟道、STI角、三维应力 节点规模:1D ≪ 2D ≪ 3D 2D→3D 节点数约 ×10²~10³ 按分析目标选最低必要维度 自适应网格细化 触发量:掺杂梯度、缺陷通量、应力 加密:梯度超阈值 → 局部细分 粗化:梯度回落 → 合并节点 RefineBox 限定细化区域 避免全局加密导致规模爆炸 每步检查节点数增长趋势 关键区域网格密度控制 沟道/结区:0.5~2 nm 级加密 氧化层界面:沿法向 0.1~0.5 nm 源漏延伸区:横向 1~3 nm 衬底/体区:10~50 nm 粗网格 原则:梯度大处密,梯度小处疏 用 MaxLen / MinLen 约束步长 网格质量对收敛与精度的影响及反馈调整 收敛性 节点过密→矩阵病态 精度 网格过粗→结深偏差 计算成本 节点数↑→时间/内存↑ 鲁棒性 畸变单元→求解失败 网格收敛性验证:加密后结果变化 < 1%

二、自适应网格细化的触发与执行

Sentaurus Process的自适应细化由求解器在每一时间步后评估解的局部变化率。当某单元内掺杂浓度梯度、缺陷浓度梯度或应力梯度超过设定阈值时,该单元被标记并细分;当梯度回落到粗化阈值以下时,相邻单元合并。关键控制参数包括RefineBox(限定细化区域)、MaxLen(最大步长)、MinLen(最小步长)和RefineFactor(细化倍率)。

一个典型的二维NMOS工艺仿真中,可在沟道和源漏结区域设置RefineBox,将MaxLen限制在2 nm以内,而在衬底深处允许MaxLen达到50 nm。这样既保证了结区精度,又避免了全片加密带来的节点爆炸。自适应细化的执行顺序是:先由line命令生成初始结构化网格,再在diffuse或anneal步骤中启用自适应选项,求解器根据当前解场动态调整。

三、关键区域网格密度控制

MOSFET工艺仿真中,以下区域必须重点控制网格密度:

  • 栅氧化层及界面:氧化层厚度通常为1~3 nm,沿厚度方向至少需要3~5个节点才能准确求解氧化生长和界面态。建议在氧化层内设置MaxLen = 0.3 nm,并在硅/氧化层界面处额外加密。
  • 源漏结与LDD浅结:结深通常在20~100 nm,结区掺杂浓度在几个数量级内变化。横向和纵向MaxLen建议为1~3 nm,确保结深提取误差小于5%。
  • 侧墙下方与沟道边缘:横向扩散和短沟道效应在此区域最敏感,建议MaxLen = 1~2 nm。
  • 衬底体区:掺杂梯度平缓,MaxLen可放宽至20~50 nm,以控制总节点数。

密度控制的核心原则是“梯度大处密,梯度小处疏”。实际操作中,可先用较粗网格跑一遍,观察掺杂分布和结深,再在梯度大的区域逐步加密,直到结果变化小于1%。

四、网格质量对收敛与精度的影响

网格质量对仿真结果的影响体现在两个层面。第一是收敛性:节点过密会导致线性方程组规模增大、矩阵条件数恶化,Newton迭代可能不收敛;节点过疏则会使离散误差增大,甚至导致求解器在强非线性区域振荡。第二是精度:网格过粗时,结深、薄层电阻、阈值电压等关键参数会出现系统性偏差。例如,当结区MaxLen从1 nm放宽到5 nm时,仿真得到的结深可能偏浅3%~8%,阈值电压偏差可达10~20 mV。

判断网格是否足够细的实用方法是网格收敛性验证:将MaxLen减半,重新仿真,比较关键参数的变化。若结深、阈值电压等变化小于1%,则认为网格已收敛。若变化仍大于1%,需继续加密或检查是否存在畸变单元。此外,应避免长宽比过大的单元(建议控制在10:1以内),因为畸变单元会显著降低求解精度和稳定性。

实用建议:在Sentaurus Process中,可先用RefineBox在关键区域设置较细的MaxLen,再在diffuse步骤中启用自适应细化。每次加密后检查节点数和仿真时间,确保在可接受的计算成本内达到收敛。对于三维仿真,优先使用各向异性网格,在垂直方向加密、横向适当放宽,以控制总节点数。

五、网格策略的工程化流程

综合以上要点,MOSFET工艺仿真的网格策略可归纳为以下工程化流程:首先根据分析目标确定最低必要维度(1D/2D/3D);然后用结构化网格建立初始几何,在关键区域设置RefineBox和MaxLen;接着启用自适应细化,让求解器在梯度大的区域自动加密;最后进行网格收敛性验证,通过加密对比确认关键参数变化小于1%。整个过程中,需持续监控节点数、仿真时间和收敛状态,在精度与成本之间取得平衡。

网格策略不是一次设定就固定不变的,而是随着工艺步骤推进动态调整的。例如,在氧化生长步骤中,氧化层界面需要加密;在退火步骤中,缺陷团簇区域需要加密;在刻蚀步骤后,新暴露的表面需要重新生成网格。理解并掌握这些动态调整方法,是MOSFET工艺仿真从“能跑”到“跑得准、跑得快”的关键。