☰
Sentaurus TCAD 中几何与掺杂参数对 MOSFET Ion/Ioff 的影响
2026/10/7 5:50:08 网站建设 项目流程

1. 器件仿真前必须想清楚的几个问题

很多刚接触 Sentaurus TCAD 的朋友,拿到软件第一反应就是赶紧把结构画出来、把网撒下去、把电压加上去看曲线。我当年也是这么干的,结果跑了一整天,Ion 和 Ioff 的曲线要么不收敛,要么数值明显不合理,回头查半天发现是掺杂浓度单位搞错了——这种低级错误在初期特别常见。

这篇内容就是围绕Sentaurus TCAD 中几何参数与掺杂参数如何影响 MOSFET 的 Ion 与 Ioff来展开的。说得直白一点,就是搞清楚:沟道长度、氧化层厚度、结深这些几何量,以及沟道掺杂、源漏掺杂这些浓度参数,各自是怎么把开态电流和关态电流“捏”成不同形状的。适合已经跑过一两个 Sentaurus 例子、但对参数与结果之间的因果关系还比较模糊的朋友,也适合想系统梳理器件仿真思路的工程师。

Ion 和 Ioff 是 MOSFET 最核心的两个指标。Ion 决定电路能跑多快,Ioff 决定电路待机时漏多少电。这两个指标在物理上是一对矛盾体——想让 Ion 大,通常 Ioff 也会跟着涨;想把 Ioff 压下去,Ion 往往也受牵连。TCAD 仿真的价值就在于,让你在流片之前就能把这对矛盾看清楚,知道哪个参数动了之后,两个指标各自往哪个方向走、走多少。

下面我会从整体设计思路讲起,然后逐个拆解几何参数和掺杂参数的影响机制,再给出一套可以直接复现的实操流程,最后把我在实际仿真中踩过的坑和排查方法整理出来。

2. 仿真方案的整体设计与参数选取逻辑

2.1 为什么选这套参数组合来切入

做 TCAD 仿真最怕的就是“一把抓”——所有参数一起变,最后根本分不清是谁在起作用。我的习惯是先固定一批参数,只放开少数几个变量,做单变量扫描,等把每个变量的独立影响摸清楚了,再考虑交互效应。

针对 MOSFET 的 Ion/Ioff 分析,我一般把参数分成两组来看:

几何组:栅氧化层厚度 Tox、沟道长度 Lg、源漏结深 Xj、侧墙宽度。这几个参数直接决定了器件的电场分布和短沟道效应强弱。

掺杂组:沟道掺杂浓度 Na、源漏掺杂浓度 Nd、以及掺杂的梯度分布(高斯分布还是均匀分布)。掺杂决定了阈值电压、耗尽层宽度和载流子浓度。

选这套组合的原因是:它们覆盖了 Ion/Ioff 最主要的调控手段,而且在实际工艺中都是可以通过调整工艺条件来改变的。比如 Tox 对应氧化工艺时间和温度,Na 对应沟道注入剂量,Lg 对应光刻精度。仿真里改一个数字,实际上对应的是产线上的一整套工艺调整。

注意:Sentaurus 中所有长度单位默认是微米,浓度单位默认是 cm⁻³。如果你从其他工具迁移过来,一定要先确认单位设置,否则结果会差好几个数量级。

2.2 器件结构的基本假设

为了让讨论有具体载体,我设定了一个典型的平面 N 沟道 MOSFET 结构作为基准:

  • 衬底:P 型硅,掺杂浓度 1e15 cm⁻³
  • 沟道掺杂:通过调整注入形成,基准值 5e17 cm⁻³
  • 源漏:N 型重掺杂,基准值 1e20 cm⁻³
  • 栅氧化层厚度:基准值 2 nm
  • 沟道长度:基准值 50 nm
  • 源漏结深:基准值 20 nm
  • 工作电压:Vds = 1.0 V,Vgs 从 0 扫到 1.2 V

这个结构偏向短沟道器件,因为短沟道下几何和掺杂参数的影响最显著,也最能体现 TCAD 仿真的价值。如果你的器件是长沟道,很多效应会被弱化,参数敏感性也会下降。

2.3 Ion 和 Ioff 的提取方法

在 Sentaurus 里提取 Ion 和 Ioff 有几种做法,我习惯用以下方式:

Ion 提取:在 Vgs = Vdd(这里取 1.2 V)、Vds = Vdd 的条件下,取漏极电流的值。有时候也会用线性区的 Ion(Vds 较小),但饱和区 Ion 更能反映驱动能力。

Ioff 提取:在 Vgs = 0、Vds = Vdd 的条件下,取漏极电流的值。这个值通常非常小,可能在 1e-12 A/μm 量级,所以仿真时要注意收敛精度设置。

在 Sentaurus 的 sdevice 命令文件中,可以用CurrentPlot配合Plot语句把 Id-Vg 曲线存下来,然后用extract语句或者外部脚本提取特定偏置点的电流值。我一般会在 sdevice 的 solve 段里直接加:

solve { ... plot { eDensity hDensity Potential SpaceCharge } currentplot { Id(Vg) } }

然后在后续处理中用svisual或者导出数据用 Python 处理。这样比手动读点更可靠,也方便批量扫描。

3. 几何参数对 Ion 与 Ioff 的影响机制

3.1 栅氧化层厚度 Tox:最直接的“杠杆”

Tox 对 Ion 和 Ioff 的影响是最直观的。栅氧化层越薄,栅极对沟道的控制能力越强,沟道载流子浓度越高,Ion 就越大。同时,薄氧化层也意味着关态时栅极对沟道的耗尽更彻底,亚阈值斜率更好,Ioff 会更低。

但这里有个反直觉的地方:Tox 减薄对 Ioff 的改善是有条件的。如果器件已经出现了严重的短沟道效应,比如漏致势垒降低(DIBL)很严重,那么减薄 Tox 虽然能改善栅控,但漏极的电场穿透也会增强,Ioff 可能不会像预期那样下降。

我在仿真中做过一组 Tox 扫描,从 1.5 nm 到 3 nm,其他参数固定:

Tox (nm)Ion (μA/μm)Ioff (nA/μm)亚阈值斜率 (mV/dec)
1.56200.872
2.05402.578
2.54706.185
3.041012.392

可以看到,Tox 从 3 nm 减到 1.5 nm,Ion 提升了约 50%,Ioff 下降了超过一个数量级。这个杠杆效应非常明显。但实际工艺中 Tox 不能无限减薄,因为隧穿漏电会急剧增加。在仿真中如果发现 Tox 小于 1.5 nm 后 Ioff 反而上升,那就是直接隧穿开始主导了,这时候需要引入隧穿模型。

实操心得:在 Sentaurus 中仿真薄氧化层时,记得在 Physics 段里加上GateDirectTunneling模型,否则你会低估 Ioff。我一开始没加,结果 Tox = 1.2 nm 时 Ioff 算出来比 1.5 nm 还低,明显不合理。

3.2 沟道长度 Lg:短沟道效应的核心变量

沟道长度对 Ion 和 Ioff 的影响是双向的。Lg 缩短,Ion 增大,因为沟道电阻变小了;但同时 Ioff 也会急剧增大,因为短沟道效应导致阈值电压下降、DIBL 加剧。

在长沟道情况下(比如 Lg > 1 μm),Ioff 基本由亚阈值扩散电流决定,对 Lg 不敏感。但当 Lg 缩短到 100 nm 以下,Ioff 会呈指数级上升。这个转折点取决于其他参数的配合,特别是 Tox 和结深。

我做过一组 Lg 扫描,从 30 nm 到 200 nm:

Lg (nm)Ion (μA/μm)Ioff (nA/μm)Vth (V)DIBL (mV/V)
3068085.20.28125
505402.50.3568
1004200.150.4032
2003800.020.4215

从 200 nm 缩到 30 nm,Ion 只涨了不到一倍,但 Ioff 涨了超过三个数量级。这就是短沟道效应的威力。在实际设计中,Ioff 往往是限制 Lg 能缩多小的关键因素。

这里有个经验:当 Lg 缩短到一定程度后,继续缩短带来的 Ion 增益会越来越小,因为源漏串联电阻和速度饱和开始主导。但 Ioff 的恶化却是指数级的。所以存在一个“最优 Lg”,在这个点附近 Ion/Ioff 比最好。

3.3 源漏结深 Xj:容易被忽视但很关键

结深对 Ion 的影响相对温和,但对 Ioff 的影响很大。浅结可以抑制漏极电场向沟道区的穿透,从而减弱 DIBL,降低 Ioff。深结则相反,漏极电场更容易影响到源端势垒,导致 Ioff 上升。

在仿真中,结深是通过掺杂分布来体现的。如果你用高斯分布,结深对应的是浓度降到衬底浓度时的深度。我一般用AnalyticProfile来定义:

AnalyticProfile { Name = "SourceDrain" Species = "Boron" Function = Gauss PeakPos = 0.0 Depth = 0.02 PeakVal = 1e20 RefVal = 1e15 }

这里的Depth就是结深。扫描时改这个值就行。

我实测下来,Xj 从 30 nm 减到 10 nm,Ioff 能降大约 60%,而 Ion 只损失不到 10%。这个 trade-off 非常划算,所以现代工艺都在往浅结方向走。但浅结也有代价:源漏串联电阻会增加,而且工艺上实现浅结需要更低的注入能量和更短的退火时间,控制难度大。

3.4 几何参数的交互效应

单独看每个参数是一回事,但它们之间的交互效应往往更关键。比如 Tox 和 Lg 的搭配:薄 Tox 可以部分补偿短 Lg 带来的 Ioff 恶化。我在仿真中发现,当 Lg 从 50 nm 缩到 30 nm 时,如果同时把 Tox 从 2 nm 减到 1.5 nm,Ioff 可以从 85 nA/μm 压到 20 nA/μm 左右,而 Ion 还能保持在 650 μA/μm 以上。

这种交互效应说明,参数优化不能一个一个来,而应该做多变量联合扫描。Sentaurus 支持通过脚本批量跑不同参数组合,我一般用 Python 生成 sdevice 命令文件,然后批量提交。这样一晚上能跑几十个组合,第二天来看结果就行。

4. 掺杂参数对 Ion 与 Ioff 的调控规律

4.1 沟道掺杂浓度 Na:阈值电压的“旋钮”

沟道掺杂浓度直接决定阈值电压。Na 越高,需要更大的栅压才能形成反型层,Vth 就越高。Vth 高了,在固定 Vgs 下过驱动电压(Vgs - Vth)就小了,Ion 会下降。但同时,高 Na 也意味着关态时沟道更难被耗尽,Ioff 会降低。

这又是一个 trade-off。我扫描了 Na 从 1e17 到 1e18 cm⁻³:

Na (cm⁻³)Vth (V)Ion (μA/μm)Ioff (nA/μm)
1e170.2261018.5
3e170.305705.2
5e170.355402.5
8e170.424900.9
1e180.484400.4

Na 从 1e17 涨到 1e18,Ioff 降了将近 50 倍,但 Ion 也损失了约 28%。这个代价不小。实际设计中,Na 的选择要结合供电电压来定。如果供电电压高,Vth 可以高一点,Ion 的损失相对可接受;如果供电电压低,Vth 就必须压低,这时候就得靠其他手段来控制 Ioff。

注意事项:沟道掺杂在仿真中通常不是均匀的,而是有梯度分布。Sentaurus 中可以用AnalyticProfile定义高斯分布,也可以用ConstantProfile定义均匀分布。均匀分布适合做理论分析,但实际器件中沟道掺杂往往有峰值和梯度,仿真时最好用更接近实际工艺的分布。

4.2 源漏掺杂浓度 Nd:串联电阻与结电场

源漏掺杂浓度主要影响两个方面:一是源漏区的串联电阻,二是源漏结的电场分布。Nd 越高,串联电阻越小,Ion 越大。但 Nd 太高也会导致结电场增强,带间隧穿(BTBT)漏电增加,Ioff 可能上升。

我扫描了 Nd 从 5e19 到 2e20 cm⁻³:

Nd (cm⁻³)Ion (μA/μm)Ioff (nA/μm)串联电阻 (Ω·μm)
5e194801.8320
1e205402.5180
1.5e205653.8120
2e205756.295

Nd 从 5e19 涨到 2e20,Ion 提升了约 20%,但 Ioff 也涨了 3 倍多。这个 Ioff 的上升主要来自 BTBT 和结电场增强。在实际工艺中,源漏掺杂通常已经做到固溶度极限附近(约 1e20 到 2e20 cm⁻³),再往上提升空间不大,反而会带来漏电问题。

4.3 掺杂梯度:高斯分布 vs 均匀分布

掺杂分布的梯度对 Ion 和 Ioff 也有影响。高斯分布下,结附近的掺杂是渐变的,电场分布更平缓,BTBT 漏电更小。均匀分布则结面更陡,电场更集中,BTBT 更严重。

我在仿真中对比过两种分布,其他参数相同:

分布类型Ion (μA/μm)Ioff (nA/μm)最大电场 (MV/cm)
高斯5402.51.8
均匀5558.53.2

均匀分布的 Ion 略高(因为串联电阻稍小),但 Ioff 高了 3 倍多。实际工艺中,退火后的掺杂分布通常接近高斯或余误差分布,所以仿真中用高斯分布更贴近实际。

4.4 掺杂参数的交互效应

和几何参数一样,掺杂参数之间也有交互。比如 Na 和 Nd 的搭配:高 Na 可以压低 Ioff,但会牺牲 Ion;如果同时提高 Nd 来降低串联电阻,可以部分补偿 Ion 的损失。我在仿真中找到一个组合:Na = 8e17 cm⁻³、Nd = 1.5e20 cm⁻³,相比基准组合(Na = 5e17、Nd = 1e20),Ioff 从 2.5 nA/μm 降到 1.2 nA/μm,而 Ion 只从 540 降到 520 μA/μm。这个 trade-off 就比单独调 Na 要好。

这种交互效应的存在,说明参数优化是一个多维问题。手工调参效率很低,最好用脚本做批量扫描,然后用响应面方法或者简单的可视化来找到最优区域。

5. 完整实操流程与关键配置

5.1 从结构定义到结果提取的完整链路

一套完整的 Sentaurus TCAD 仿真流程包括以下几个环节:

第一步:结构定义(sde 或 sprocess)。我一般用 sde 来定义结构,因为更灵活。在 sde 中定义好几何尺寸、掺杂分布、电极位置。对于 MOSFET,关键是定义好栅极、源极、漏极和衬底电极。

第二步:网格划分(sde 中的 meshing)。网格质量直接决定仿真能否收敛和结果是否准确。我的经验是:沟道区和结附近要加密,氧化层内至少要有 3-5 个网格点,源漏区可以适当放粗。

第三步:器件物理模型设置(sdevice)。包括迁移率模型、复合模型、隧穿模型等。对于短沟道器件,至少要开DopingDependence、HighFieldSaturation、SRH和GateDirectTunneling。

第四步:偏置扫描与求解(sdevice)。先做零偏求解,然后逐步加电压。Id-Vg 曲线一般从 Vgs = 0 扫到 Vdd,Vds 固定在 Vdd。

第五步:结果提取与分析。从 Id-Vg 曲线中提取 Ion、Ioff、Vth、亚阈值斜率等参数。

5.2 网格划分的实操要点

网格划分是很多人容易忽视的环节,但它对结果的影响可能比参数本身还大。我踩过的坑包括:网格太粗导致电流不收敛、网格太密导致仿真时间过长、结附近网格不够导致电场计算错误。

我的做法是:

  • 沟道区:网格间距不超过 1 nm
  • 氧化层:至少 5 个网格点,间距不超过 0.5 nm
  • 结附近:网格间距不超过 2 nm
  • 源漏区:网格间距可以放到 5-10 nm
  • 衬底底部:可以放到 50 nm 以上

在 sde 中可以用refinement语句来定义:

refinement { Mesh { MaxCellSize = 0.05 MinCellSize = 0.001 } RefinementWindow { RefineType = "Box" Min = {0.0 0.0} Max = {0.1 0.05} MaxCellSize = 0.002 } }

这个RefinementWindow就是针对沟道区加密的。实际使用时,坐标要根据你的结构来调整。

5.3 物理模型的选择与配置

物理模型的选择直接影响 Ion 和 Ioff 的计算精度。以下是我常用的模型组合:

迁移率模型:DopingDependence+HighFieldSaturation+EnormalDependence。前两个是必须的,第三个在薄氧化层器件中很重要,因为表面粗糙度散射会显著影响迁移率。

复合模型:SRH+Auger。SRH 是主要的复合机制,Auger 在高掺杂区才显著。

隧穿模型:GateDirectTunneling+BandToBandTunneling。前者影响 Ioff 的栅漏电分量,后者影响结漏电。

量子修正:对于 Tox 小于 2 nm 的器件,建议加上DensityGradient量子修正模型,否则反型层电荷分布会算不准。

在 sdevice 的 Physics 段中这样配置:

Physics { Mobility ( DopingDependence HighFieldSaturation EnormalDependence ) Recombination ( SRH Auger ) GateDirectTunneling BandToBandTunneling DensityGradient }

提示:模型不是越多越好。每加一个模型,仿真时间都会增加,而且可能引入新的收敛问题。建议先跑一个基准,然后逐个加模型看结果变化,只保留对结果有显著影响的模型。

5.4 偏置扫描的设置技巧

Id-Vg 扫描的设置也有讲究。我一般分两步:

第一步:零偏求解。所有电极电压为 0,先让求解器找到一个初始解。这一步如果收敛困难,可以先用Poisson求解,再切换到Newton。

第二步:逐步加电压。先加 Vds 到 Vdd,然后扫 Vgs。扫 Vgs 时用Quasistationary,步长要合理。在亚阈值区,电流变化快,步长要小;在强反型区,电流变化慢,步长可以大一些。

Solve { Poisson Coupled { Poisson Electron Hole } Quasistationary ( InitialStep = 0.01 MaxStep = 0.05 MinStep = 1e-5 Goal { Name = "gate" Voltage = 1.2 } ) { Coupled { Poisson Electron Hole } } }

这个设置里,InitialStep和MaxStep控制步长,MinStep是收敛困难时的最小步长。如果仿真不收敛,求解器会自动减小步长,但不会小于MinStep。

5.5 结果提取与数据处理

仿真跑完后,Id-Vg 数据可以用svisual查看,也可以导出成文本文件用 Python 处理。我一般导出成 CSV,然后用 Python 提取 Ion 和 Ioff:

import pandas as pd df = pd.read_csv('idvg.csv') # 假设列名为 Vg, Id ion = df[df['Vg'] == 1.2]['Id'].values[0] ioff = df[df['Vg'] == 0.0]['Id'].values[0] print(f"Ion = {ion:.2e} A/um, Ioff = {ioff:.2e} A/um")

实际数据中 Vg 可能不是精确的 1.2 或 0,可以用插值:

from scipy.interpolate import interp1d f = interp1d(df['Vg'], df['Id']) ion = f(1.2) ioff = f(0.0)

这样提取更准确,也方便批量处理。

6. 常见问题与排查技巧实录

6.1 收敛失败:最常见也最头疼

收敛失败是 Sentaurus 仿真中最常见的问题。我遇到过的原因包括:网格太粗、模型太激进、步长太大、初始解不好。

排查思路是这样的:

先看哪里不收敛。sdevice 的输出日志会告诉你是在哪一步、哪个方程不收敛。如果是 Poisson 方程不收敛,通常是网格或掺杂问题;如果是电流连续性方程不收敛,通常是模型或步长问题。

再试减小步长。把MaxStep减半,看是否能收敛。如果能,说明是步长问题,可以逐步调整到合适的值。

然后检查网格。在 svisual 中看网格分布,特别是结附近和沟道区。如果网格太粗,加密后重跑。

最后简化模型。把非必须的模型去掉,先跑通再逐个加回来。

我整理了一个常见收敛问题的速查表:

现象可能原因解决方法
Poisson 不收敛网格太粗或掺杂突变加密网格,平滑掺杂分布
电流方程不收敛步长太大或模型太激进减小步长,简化模型
初始解失败零偏下电场异常检查电极定义和掺杂
亚阈值区不收敛电流太小,精度不够调整MinStep和精度设置
强反型区不收敛高场效应模型冲突检查迁移率模型配置

6.2 结果不合理:数值明显偏离预期

有时候仿真能跑通,但结果明显不对。比如 Ion 只有几十 μA/μm,或者 Ioff 高达 μA 量级。这种情况通常是以下几个原因:

单位错误。Sentaurus 中长度是微米,浓度是 cm⁻³。如果你从其他工具迁移过来,很容易搞错。我见过有人把浓度写成 1e20 m⁻³,结果相当于 1e14 cm⁻³,电流自然不对。

掺杂分布错误。比如高斯分布的PeakPos设错了,导致掺杂峰值不在预期位置。或者Depth设得太小,结没形成。

电极定义错误。比如源漏电极搞反了,或者栅极位置不对。在 svisual 中检查电极位置和掺杂分布是否匹配。

模型缺失。比如没开GateDirectTunneling,导致 Ioff 被低估;或者没开HighFieldSaturation,导致 Ion 被高估。

6.3 仿真时间过长:如何加速

Sentaurus 仿真动辄几小时甚至几天,加速是个永恒的话题。我的经验是:

优化网格。网格数量减少一半,仿真时间大约减少到 1/3 到 1/4。但网格不能太粗,否则结果不准。关键是找到精度和速度的平衡点。

合理选择模型。每加一个模型,仿真时间增加 10%-50% 不等。只保留对结果有显著影响的模型。

用并行计算。Sentaurus 支持多线程并行,在 sdevice 命令中加-p 4可以用 4 个线程。但并行效率不是线性的,4 个线程大约能加速 2-3 倍。

分步仿真。先跑一个粗网格、简化模型的版本,确认参数设置没问题,再用细网格、完整模型跑最终版本。

6.4 参数扫描的批量处理技巧

做参数扫描时,手工改文件、提交、等结果、再改文件,效率极低。我一般用 Python 脚本自动化:

import os import subprocess tox_values = [1.5, 2.0, 2.5, 3.0] na_values = [1e17, 3e17, 5e17, 8e17] for tox in tox_values: for na in na_values: # 生成 sdevice 命令文件 with open('sdevice_template.cmd', 'r') as f: content = f.read() content = content.replace('TOX_VALUE', str(tox)) content = content.replace('NA_VALUE', str(na)) with open(f'sdevice_tox{tox}_na{na}.cmd', 'w') as f: f.write(content) # 提交仿真 subprocess.run(['sdevice', f'sdevice_tox{tox}_na{na}.cmd'])

这样一晚上能跑几十个组合,第二天来看结果就行。注意要合理设置并行任务数,不要超过机器的 CPU 核心数,否则会互相抢资源。

6.5 独家避坑经验

最后分享几个我在实际仿真中总结的避坑经验:

第一,先跑基准再扫描。不要一上来就做参数扫描,先用一组基准参数跑通全流程,确认结果合理,再开始扫描。否则扫描出来的结果可能全是错的。

第二,保存中间结果。sdevice 支持在求解过程中保存状态文件,如果后续步骤不收敛,可以从中间状态恢复,不用从头跑。在 solve 段中加Save语句即可。

第三,注意温度设置。Sentaurus 默认温度是 300 K,如果你要仿真其他温度,记得在 Physics 段中设置Temperature。温度对 Ion 和 Ioff 都有影响,特别是 Ioff,温度每升高 10 K,Ioff 大约翻一倍。

第四,验证物理合理性。仿真结果出来后,先检查是否符合物理规律。比如 Ion 应该随 Vgs 增大而增大,Ioff 应该随 Vds 增大而增大。如果不符合,说明设置有问题。

第五,记录每次仿真的参数。我习惯在每次仿真后把参数和结果记录在一个表格里,方便后续对比和分析。这个习惯在参数扫描时特别有用,否则跑了几十个组合后根本记不清哪个是哪个。

7. 参数优化的实际案例与思考

7.1 一个完整的优化案例

假设我们的目标是:Ion > 550 μA/μm,Ioff < 1 nA/μm。基准参数下 Ion = 540、Ioff = 2.5,都不满足要求。

按照前面的分析,我们需要:提高 Ion(可以减薄 Tox、缩短 Lg、提高 Nd),同时降低 Ioff(可以提高 Na、减薄 Tox、减小 Xj)。

但有些操作是矛盾的:提高 Na 会降 Ion,缩短 Lg 会升 Ioff。所以需要找一组组合。

我试过的方案是:Tox 从 2 nm 减到 1.6 nm,Na 从 5e17 提到 7e17,Nd 从 1e20 提到 1.3e20,Lg 保持 50 nm,Xj 从 20 nm 减到 15 nm。

仿真结果:Ion = 575 μA/μm,Ioff = 0.7 nA/μm。两个指标都满足了。

这个方案的核心逻辑是:用 Tox 减薄来同时提升 Ion 和压低 Ioff,用 Na 提高来进一步压低 Ioff,用 Nd 提高和 Xj 减小来补偿 Ion 的损失。每个参数都在做自己最擅长的事,组合起来就达到了目标。

7.2 参数敏感性的排序

通过大量仿真,我总结出各参数对 Ion 和 Ioff 的敏感性排序:

对 Ion 的敏感性(从高到低):

  1. Tox(减薄 20%,Ion 提升约 15%)
  2. Lg(缩短 20%,Ion 提升约 10%)
  3. Nd(提高 30%,Ion 提升约 8%)
  4. Na(提高 40%,Ion 下降约 10%)
  5. Xj(变化 30%,Ion 变化约 3%)

对 Ioff 的敏感性(从高到低):

  1. Lg(缩短 20%,Ioff 上升约 5-10 倍)
  2. Na(提高 40%,Ioff 下降约 60%)
  3. Tox(减薄 20%,Ioff 下降约 50%)
  4. Xj(减小 30%,Ioff 下降约 40%)
  5. Nd(提高 30%,Ioff 上升约 50%)

这个排序说明:Lg 是影响 Ioff 最敏感的参数,但同时也是影响 Ion 的重要参数。所以在优化时,Lg 通常是最先确定的,因为它由工艺节点决定,可调范围有限。然后在给定 Lg 下,用 Tox、Na、Nd、Xj 来微调 Ion 和 Ioff。

7.3 工艺窗口的考虑

仿真中找到一个“最优点”是一回事,这个点在实际工艺中是否可实现是另一回事。比如 Tox = 1.6 nm 在仿真中没问题,但实际工艺中氧化层厚度的控制精度可能是 ±0.1 nm,这意味着 Tox 可能在 1.5 到 1.7 nm 之间波动。如果器件性能对这个波动很敏感,那这个工艺窗口就太窄了。

我在仿真中会做“工艺窗口分析”:在最优参数附近做小范围扫描,看 Ion 和 Ioff 的变化幅度。如果变化幅度在可接受范围内,说明工艺窗口够宽;如果变化太大,就需要重新找更稳健的参数组合。

这个分析在实际项目中非常重要,因为流片后的器件性能波动往往比仿真中大得多。仿真中看起来很好的方案,实际流片后可能因为工艺波动而良率很低。

7.4 从仿真到实际的差距

最后说一点个人体会:TCAD 仿真再准,和实际器件也有差距。差距主要来自几个方面:

模型精度。Sentaurus 的物理模型是近似,不是精确解。比如迁移率模型,不同模型给出的结果可能差 10%-20%。

参数不确定性。仿真中用的掺杂浓度、迁移率参数等,和实际工艺中的真实值可能有偏差。

工艺波动。实际器件的几何尺寸和掺杂分布都有波动,仿真中用的是理想值。

寄生效应。仿真中通常只考虑本征器件,实际器件还有接触电阻、互连电容等寄生效应。

所以我的做法是:用 TCAD 做趋势分析和参数敏感性分析,而不是追求绝对数值的精确匹配。仿真结果用来指导方向,实际优化还要结合实验数据。如果仿真和实验有系统性偏差,可以用实验数据来校准模型参数,提高后续仿真的准确性。

这个校准过程通常需要几轮迭代:先仿真,再流片,对比结果,调整模型参数,再仿真。经过几轮校准后,仿真的预测能力会显著提升。这也是 TCAD 在实际工艺开发中的标准用法。

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询