2. 工艺波动来源与统计描述

在Sentaurus TCAD中做良率敏感性分析，第一步不是急着跑仿真，而是把“波动从哪里来、用什么统计量描述、彼此是否相关”这三件事定义清楚。工艺波动本质上是一组随机变量：每个器件、每个晶圆、每个批次上的某个工艺量都不是单一确定值，而是服从某种分布。TCAD中的统计仿真，就是把这些分布映射到器件电学参数（如Vth、Ion、Ioff、SS）的分布上。本节围绕五类主要波动源展开，并给出分布建模与相关性定义的方法。

一、五类主要工艺波动源

随机掺杂波动（RDF）是纳米尺度器件中最主要的波动源之一。当沟道体积小到只容纳几百到几千个掺杂原子时，离散的杂质位置和数量会显著影响耗尽层电荷与阈值电压。RDF的物理本质是“离散性”，不是连续浓度误差。在Sentaurus中通常用两种方式处理：一是解析模型（如随机阈值电压偏移模型），二是基于SDE/SPROCESS的原子级离散掺杂。RDF的统计描述常用标准差σ(ΔVth)随沟道面积缩小而增大的关系，近似满足σ ∝ 1/√(W·L)。

线边缘粗糙度（LER/LWR）来自光刻与刻蚀过程中光刻胶边缘的随机起伏。它表现为栅极或鳍片边缘偏离理想直线的随机偏移，典型幅度为1～3 nm，相关长度约10～30 nm。LER会改变有效沟道长度，进而影响短沟道效应和Vth。在Sentaurus中可用随机边缘轮廓生成器（如基于自相关函数的粗糙度模型）叠加到结构上，再逐次仿真。

栅氧厚度波动包括平均厚度偏差和局部厚度起伏。平均偏差属于全局（die-to-die或wafer-to-wafer）波动，局部起伏属于局部（device-to-device）波动。栅氧厚度直接影响栅电容Cox和阈值电压，其敏感度约为ΔVth/Δtox ≈ -Vth/tox量级。在统计仿真中通常把tox建模为高斯分布，σ/mean约为1%～3%。

刻蚀、沉积与退火波动属于工艺设备与工艺条件引起的波动。刻蚀波动表现为栅长Lg、侧壁角、刻蚀深度的偏差；沉积波动表现为薄膜厚度、应力的偏差；退火波动表现为温度、时间、升降温速率的偏差，进而影响掺杂激活与扩散长度。这类波动往往具有空间相关性（wafer map上的径向或条纹分布），不能简单当作独立同分布处理。

工艺参数分布建模与相关性定义是把上述物理波动转化为可仿真的统计模型。核心工作包括：为每个波动量选择分布类型（正态、对数正态、均匀、截断正态等），设定均值μ和标准差σ，并定义变量之间的相关系数矩阵ρ。例如Lg与LER可能相关，tox与退火温度可能通过氧化速率相关。Sentaurus中可用CorrelationMatrix或协方差矩阵描述这些关系，再通过拉丁超立方采样（LHS）或蒙特卡洛采样生成样本。

二、分布建模的具体做法

在Sentaurus中建立统计模型时，通常按以下顺序操作。首先确定每个波动量的分布类型。对于RDF引起的Vth偏移，常用正态分布；对于Ioff这类跨数量级变化的量，常用对数正态分布；对于Lg、tox等有明确工艺窗口的量，常用截断正态分布，避免出现负值或超出物理极限的样本。

其次设定均值和标准差。均值来自工艺目标值（nominal），标准差来自工艺能力数据或文献。例如某28 nm工艺中，Lg目标值30 nm，σ(Lg)=1.2 nm；tox目标值1.8 nm，σ(tox)=0.03 nm；RDF引起的σ(ΔVth)在W·L=0.1 μm²时约为25 mV。这些数值属于教学示例，实际项目应使用自己的工艺数据。

然后定义相关性。相关系数矩阵必须是对称正定的，对角线为1。例如三个变量Lg、LER、tox的相关系数矩阵可写为：

ρ = [ 1.00   0.35   0.10
      0.35   1.00   0.05
      0.10   0.05   1.00 ]

在Sentaurus中可通过CorrelationMatrix语句或外部脚本生成相关样本。若使用LHS采样，建议样本数N至少为变量数的10倍以上；对于5个变量，N≥100组；若需要估计σ(Vth)的置信区间，N建议取200～500组。

三、从波动源到良率指标的映射

统计仿真的输出不是单一曲线，而是一组器件特性分布。以Vth为例，仿真完成后可得到N个Vth值，计算其均值μ(Vth)和标准差σ(Vth)。若规格要求Vth在[0.25 V, 0.35 V]之间，则可统计超出规格的样本比例，即失效率。进一步可计算Cpk = min[(USL-μ)/(3σ), (μ-LSL)/(3σ)]，Cpk越大表示工艺窗口越充裕。

敏感性分析则回答“哪个波动源贡献最大”。常用方法是计算每个输入变量与输出Vth的相关系数或回归系数。例如若ρ(Lg, Vth) = -0.72，ρ(RDF, Vth) = 0.45，ρ(tox, Vth) = -0.30，则说明Lg波动对Vth影响最显著，应优先控制Lg的工艺稳定性。Sentaurus的统计模块可直接输出这些相关系数，也可通过外部脚本对仿真结果做回归分析。