Silvaco Atlas仿真结果解析与常见报错排查实战指南

1. 半导体工程师必备：Silvaco Atlas仿真结果深度解析与常见报错排查指南

做半导体器件仿真这行十几年，Silvaco Atlas 是我用得最顺手也最让人头疼的工具之一。说它顺手，是因为从 PN 结到 HEMT、从击穿特性到电容电压曲线，它几乎能覆盖器件物理的方方面面；说它头疼，是因为它的报错信息往往只有一行，但背后可能牵扯到网格、材料参数、物理模型、数值收敛等一大堆问题。很多刚入行的朋友拿到仿真结果，要么看不懂曲线为什么长这样，要么遇到报错就卡住，翻手册也找不到答案。这篇内容就是把我这些年踩过的坑、总结出来的排查思路和结果解读方法，系统地梳理一遍。不管你是刚接触 TCAD 的在校研究生，还是已经工作几年但主要靠试错推进的工程师，应该都能从中找到能直接用的东西。

Silvaco Atlas 本质上是一个基于漂移扩散模型和泊松方程耦合求解的二维/三维器件仿真器。你给它结构、掺杂、电极和物理模型，它给你电势分布、载流子浓度、电流密度、I-V 曲线、C-V 曲线这些结果。听起来简单，但实际用起来，结果合不合理、报错怎么解，才是真正拉开水平的地方。下面我按“设计思路—核心细节—实操过程—问题排查”这条线，把整个流程拆开讲。

2. 仿真整体设计与思路拆解

2.1 为什么仿真流程要分阶段推进

很多人拿到一个器件结构，第一反应是直接把所有物理模型打开，网格铺得很密，然后跑一个完整的 I-V 扫描。结果往往是跑了一晚上，要么不收敛，要么结果明显不对。我的习惯是把仿真分成三个阶段：结构搭建与基础验证、物理模型逐步加载、目标特性精细扫描。

第一阶段只做结构定义和掺杂，用最简单的漂移扩散模型跑一个零偏压求解，看看电势和载流子分布是否符合预期。这一步的目的是确认几何结构和掺杂没有低级错误。第二阶段逐步加入迁移率模型、复合模型、隧穿模型等，每加一个模型就重新跑一次基础偏压，观察结果变化是否合理。第三阶段才是针对目标特性做精细的电压或频率扫描。

这样做的逻辑很简单：TCAD 仿真的变量太多，如果一次性全打开，出了问题你根本不知道是哪个环节导致的。分阶段推进虽然看起来慢，但实际上省时间。我见过太多人因为一次性加载所有模型导致不收敛，然后花几天时间逐个排除，还不如一开始就一步步来。

2.2 网格策略：什么时候该密，什么时候该疏

网格是 Atlas 仿真里最容易被忽视但又最影响结果和收敛性的因素。新手常见的两个极端：要么全局均匀密网格，导致计算量爆炸；要么网格太疏，关键区域根本解析不出来。

我的经验是，网格密度应该跟物理量的梯度挂钩。PN 结耗尽区、异质结界面、沟道区域、碰撞电离发生区，这些地方电势和载流子浓度变化剧烈，网格必须密。而衬底底部、电极边缘外侧这些地方，物理量变化平缓，网格可以适当放粗。

具体操作上，Atlas 提供了x.mesh和y.mesh命令，可以通过spacing参数控制局部网格密度。比如在 PN 结附近，我会把 spacing 设到 0.5 nm 到 1 nm 量级，而在衬底区域可以放到 50 nm 甚至更大。另外，electrode定义处的网格也需要加密，因为电极边缘的电场集中效应很容易导致收敛困难。

注意：网格加密不是越多越好。当网格小到一定程度，数值误差反而可能增大，而且计算时间会急剧上升。一般原则是，关键区域的网格尺寸应该小于德拜长度的三分之一。

2.3 物理模型选型：别把该开的漏了，也别把不该开的硬塞

Atlas 的物理模型非常多，从基本的漂移扩散到复杂的量子修正、热电子发射、Fowler-Nordheim 隧穿等。选型的原则是：根据器件类型和工作条件，选择必要的模型，而不是把所有模型都打开。

举个例子，做硅基 PN 结二极管的反向击穿仿真，你必须打开碰撞电离模型（impact selb或impact grant），否则看不到击穿。但如果你做的是低电压 MOS 管的输出特性，碰撞电离可能不是主导，打开反而增加收敛难度。再比如，做 GaN HEMT 仿真时，极化效应模型（polarization）是必须的，否则二维电子气根本出不来。而做硅器件时，这个模型完全不需要。

我一般会参考以下几个维度来选模型：材料体系（硅、锗、GaAs、GaN、SiC 等）、器件类型（二极管、MOSFET、HEMT、HBT 等）、工作条件（低场、高场、高频、高温等）、目标特性（I-V、C-V、击穿、噪声等）。每个维度对应一组常用模型，具体可以查 Atlas 手册的模型章节，但更重要的是靠实际跑一遍看结果是否合理。

3. 核心细节解析与实操要点

3.1 掺杂定义：均匀掺杂、高斯掺杂和解析掺杂怎么选

Atlas 支持多种掺杂定义方式，最常用的是doping uniform、doping gaussian和doping analytic。均匀掺杂适合衬底和简单结构，高斯掺杂适合离子注入后的杂质分布，解析掺杂则适合需要精确控制分布形状的场景。

高斯掺杂的关键参数是peak（峰值浓度）、char（特征长度）和junction（结深位置）。这里有个容易搞错的地方：char不是标准差，而是高斯分布中浓度降到峰值 1/e 处的距离。如果你用注入能量和剂量来估算，需要先通过蒙特卡洛工具（比如 SRIM）得到射程和横向散布，再换算成 Atlas 的参数。

解析掺杂用doping analytic配合function参数，可以定义 erfc、高斯、指数等分布形式。我一般在高斯掺杂无法很好拟合实际 SIMS 曲线时，才会用解析掺杂做精细调整。

实操心得：掺杂定义完成后，一定要用tonyplot把掺杂分布画出来看一眼。我遇到过不止一次因为junction参数设错，导致结位置偏了几十纳米，结果整个仿真全废的情况。

3.2 电极定义与边界条件：欧姆接触和肖特基接触的区别

Atlas 里的电极分为欧姆接触和肖特基接触两种。欧姆接触用electrode name=...定义，默认就是欧姆接触。肖特基接触需要额外指定workfunction参数。

欧姆接触的物理含义是电极处载流子浓度等于平衡态浓度，电势等于外加偏压。肖特基接触则需要定义金属功函数，Atlas 会根据功函数和半导体电子亲和能计算势垒高度。

这里有个常见误区：很多人以为定义了electrode就完事了，但实际上电极处的网格密度、电极覆盖范围、电极边缘的曲率都会影响结果。特别是肖特基接触，如果网格不够密，势垒区的电场分布会严重失真。

3.3 求解方法选择：Newton、Gummel 和 Block 的适用场景

Atlas 提供了多种非线性求解方法，最常用的是 Newton 法和 Gummel 法。Newton 法收敛快，但对初值敏感，适合结构简单、物理模型不太复杂的情况。Gummel 法收敛慢但更稳定，适合多物理场耦合、模型复杂的场景。

我的习惯是：先用 Gummel 法跑一个零偏压解，作为初值，然后切换到 Newton 法做偏压扫描。Atlas 里可以通过method gummel newton这样的组合来指定。另外，block方法适合求解泊松方程和连续性方程分离的情况，但一般器件仿真用得少。

还有一个关键参数是maxtrap，控制数值陷阱的次数。默认值通常是 4，如果收敛困难可以适当增大，但太大会导致结果不准确。

4. 实操过程与核心环节实现

4.1 从零搭建一个 PN 结二极管仿真

下面我以一个硅基 PN 结二极管为例，把完整流程走一遍。这个例子虽然简单，但涵盖了 Atlas 仿真的核心步骤，适合用来理解整个框架。

首先是结构定义。假设我们做一个 1 μm × 1 μm 的二极管，P 区在左，N 区在右，结面在 x = 0.5 μm 处。

go atlas mesh space.mult=1.0 x.mesh loc=0.0 spacing=0.05 x.mesh loc=0.5 spacing=0.005 x.mesh loc=1.0 spacing=0.05 y.mesh loc=0.0 spacing=0.005 y.mesh loc=0.5 spacing=0.05 y.mesh loc=1.0 spacing=0.1

这里space.mult是全局网格倍率，一般设为 1.0。x 方向在结面附近加密到 5 nm，y 方向在表面附近加密到 5 nm。

然后是区域和电极定义。

region num=1 material=silicon electrode name=anode top electrode name=cathode bottom

接着是掺杂。

doping uniform conc=1e18 p.type region=1 doping gaussian conc=1e20 peak=0.5 char=0.05 n.type

这里 P 区均匀掺杂 1e18 cm⁻³，N 区用高斯掺杂，峰值 1e20 cm⁻³，特征长度 0.05 μm。

物理模型和求解方法。

models srh auger bgn method gummel newton

SRH 和 Auger 复合模型是基础，bgn 是带隙窄化模型。求解方法先用 Gummel 再用 Newton。

最后是求解和结果输出。

solve init save outf=pn_init.str log outf=pn_iv.log solve vanode=0.0 vstep=0.1 vfinal=1.0 name=anode

这段代码先求解零偏压，保存结构文件，然后从 0 V 扫描到 1 V，步长 0.1 V，结果保存到 log 文件。

4.2 结果文件解读：log 文件、str 文件和 tonyplot 可视化

Atlas 的输出文件主要有三类：log 文件记录电学特性，str 文件保存结构信息，tonyplot 用于可视化。

log 文件里最重要的是 I-V 数据。你可以用 tonyplot 直接打开 log 文件，它会自动画出电流随电压的变化曲线。但要注意，log 文件里的电流是总电流，如果你需要分区域的电流密度，需要在仿真时用save命令保存特定时刻的结构文件，然后在 tonyplot 里查看电流密度分布。

str 文件保存的是网格、掺杂、电势、载流子浓度等信息。用 tonyplot 打开后，可以画出电势分布、载流子浓度分布、电场分布等。我一般会重点看几个东西：结区的耗尽宽度是否合理、峰值电场位置是否正确、载流子浓度在电极处是否达到欧姆接触条件。

注意：tonyplot 默认的配色方案有时候不太直观，建议在显示设置里把等值线调密一些，特别是看电势和电场分布时。

4.3 关键参数计算：如何从仿真结果提取器件参数

仿真跑完之后，真正有价值的是从结果里提取器件参数。以 PN 结二极管为例，几个关键参数包括：

• 开启电压：从 I-V 曲线里找电流达到某个阈值（比如 1 A/cm²）对应的电压。硅二极管一般在 0.6-0.7 V 左右。
• 反向饱和电流：反向偏压下的电流密度，跟掺杂浓度和结面积有关。
• 击穿电压：反向电流急剧增大时的电压，需要碰撞电离模型才能看到。
• 理想因子：从 I-V 曲线的斜率提取，反映复合电流和扩散电流的比例。

提取这些参数时，我一般会用 tonyplot 的曲线拟合功能，或者把 log 文件导出到 Excel 里做后处理。Atlas 本身也提供了一些提取命令，比如extract可以计算特定偏压下的电流、电容等。

5. 常见报错与排查技巧实录

5.1 收敛失败：最常见但也最让人抓狂的问题

收敛失败是 Atlas 仿真里出现频率最高的问题。报错信息通常是 “Newton iteration failed to converge” 或者 “Gummel iteration failed to converge”。原因可能有很多，我按排查优先级列一下：

我遇到最多的情况是网格问题。特别是做击穿仿真时，碰撞电离率对电场极其敏感，如果网格不够密，电场计算误差会导致电离率偏差几个数量级，直接导致不收敛。这时候把结区网格加密到 0.1 nm 量级往往能解决问题。

另一个常见原因是偏压步长太大。Atlas 默认的vstep是 0.1 V，但在击穿附近，电流变化剧烈，0.1 V 的步长可能太大。这时候可以用solve vanode=... vstep=0.01把步长减小，或者用ramp命令让 Atlas 自动调整步长。

5.2 结果异常：曲线形状不对、数值离谱

有时候仿真能跑完，但结果明显不对。比如 I-V 曲线在正向偏压下电流不增反降，或者 C-V 曲线在某个电压下突然跳变。这类问题通常跟物理模型或参数设置有关。

一个经典问题是迁移率模型选错。硅器件在低场下用常数迁移率模型就够了，但在高场下必须用场相关迁移率模型（比如fldmob）。如果该开的时候没开，电流会被高估。反过来，如果做低场仿真时开了高场模型，可能会引入不必要的数值噪声。

另一个常见问题是复合模型。SRH 复合需要定义载流子寿命，Atlas 默认值是 1e-7 s，但实际器件的寿命可能从纳秒到毫秒不等。如果寿命设得不对，反向饱和电流会差好几个数量级。

实操心得：每次修改物理模型或参数后，先跑一个简单的 I-V 曲线，跟之前的版本对比。如果曲线形状发生剧烈变化，说明修改有问题，需要回退检查。

5.3 报错信息速查表

下面这张表是我这些年积累的常见报错和对应解法，可以直接拿来查。

5.4 独家避坑技巧：那些手册上不会写的东西

第一个技巧：保存中间状态。Atlas 的save命令可以在任意求解步骤保存结构文件。我习惯在每个关键偏压点都保存一个 str 文件，这样如果后续仿真出问题，可以从最近的保存点重新开始，不用从头跑。

第二个技巧：用tonyplot做动画。如果你保存了多个偏压下的 str 文件，可以用 tonyplot 的动画功能看电势和载流子分布随偏压的变化。这对理解器件物理非常有帮助，特别是看耗尽区扩展和沟道形成过程。

第三个技巧：日志文件要分段看。Atlas 的 log 文件会记录每一步求解的收敛情况。如果仿真中途不收敛，log 文件里会显示是哪一步、哪个变量出了问题。我一般会用文本编辑器打开 log 文件，搜索 “failed” 或 “error”，快速定位问题。

第四个技巧：别忽视space.mult参数。这个参数控制全局网格倍率，默认是 1.0。如果你发现网格太密导致计算太慢，可以把它设成 0.5 或 0.8，整体网格会变疏。但要注意，这会影响所有区域的网格密度，关键区域可能因此变得太疏。

第五个技巧：材料参数要核对。Atlas 内置的材料参数（如硅的禁带宽度、电子亲和能、迁移率等）是默认值，但实际器件可能因为工艺差异而不同。我一般会查最新文献，把关键参数手动覆盖。比如 GaN 的极化系数，不同文献给出的值差异很大，选错了结果会差很多。

6. 仿真效率优化与结果验证

6.1 如何缩短仿真时间

Atlas 仿真的计算量主要取决于网格数量和物理模型复杂度。缩短仿真时间的核心思路是：在保证精度的前提下，尽量减少网格数量和模型数量。

一个有效的方法是使用非均匀网格。在关键区域加密，在非关键区域放粗。我一般会把衬底底部的网格放到 100 nm 甚至更大，因为那里的物理量变化很平缓。另一个方法是使用method命令里的itlimit参数限制迭代次数，避免在不收敛的情况下无限循环。

还有一个技巧是使用solve命令的previous参数，让 Atlas 用上一次求解的结果作为初值。这在做偏压扫描时特别有用，可以显著减少迭代次数。

6.2 结果验证：怎么知道仿真结果是对的

仿真结果对不对，不能只看曲线好不好看。我一般会从三个层面验证：

第一，物理合理性检查。比如 PN 结的开启电压应该在 0.6-0.7 V 左右，如果仿真出来是 0.3 V 或 1.0 V，那肯定有问题。再比如，MOSFET 的阈值电压应该跟氧化层厚度和掺杂浓度匹配，如果偏差太大，说明模型或参数有问题。

第二，与解析模型对比。对于简单结构，可以用解析公式估算关键参数，然后跟仿真结果对比。比如 PN 结的耗尽宽度可以用W = sqrt(2*eps*(Vbi-V)/(q*N))估算，如果仿真结果跟这个公式差很多，说明掺杂或网格有问题。

第三，与实验数据对比。如果有实验数据，那是最好的验证。但要注意，实验数据可能受到工艺偏差、测量误差等因素影响，对比时要留出合理的误差范围。

注意：仿真结果跟实验数据对不上，不一定是仿真错了。我遇到过很多次是实验数据本身有问题，比如接触电阻没扣掉、测量设备校准不准等。所以对比之前，先确认实验数据的可靠性。

6.3 从仿真到设计：如何用结果指导器件优化

仿真的最终目的是指导设计。以 MOSFET 为例，如果你发现仿真出来的阈值电压偏高，可以通过调整沟道掺杂浓度或氧化层厚度来优化。如果发现击穿电压不够，可以调整漂移区浓度和厚度。

我一般会做参数扫描，比如把沟道掺杂浓度从 1e17 扫到 1e18，看阈值电压怎么变。Atlas 支持用loop命令做参数扫描，可以自动跑多个仿真并保存结果。这样你就能得到一条设计曲线，直接用来选参数。

另外，Atlas 还支持与优化工具联动，比如用optimizer做自动优化。不过我个人更倾向于手动扫描，因为自动优化有时候会陷入局部最优，而且不容易理解背后的物理机制。

7. 一些个人体会

做 TCAD 仿真这些年，我最大的体会是：仿真结果的可信度，取决于你对器件物理的理解程度。如果你不理解 PN 结的耗尽区怎么形成、MOSFET 的沟道怎么导通，那仿真对你来说就是一个黑盒子，跑出什么结果你也不知道对不对。反过来，如果你对物理机制很清楚，仿真就是一个强大的验证和优化工具。

另外，别怕报错。Atlas 的报错信息虽然有时候很模糊，但每一个报错背后都有明确的物理或数值原因。我遇到报错的时候，第一反应不是去搜解决方案，而是先想：这个报错可能跟哪个物理过程有关？是网格问题、模型问题还是参数问题？想清楚了再动手改，往往比盲目试错快得多。

最后，多保存、多对比、多记录。我习惯给每个仿真项目建一个文件夹，里面按日期和版本号保存 inp 文件、log 文件和 str 文件。每次修改参数或模型，都新建一个版本，而不是直接覆盖。这样如果新版本结果不对，可以随时回退到旧版本对比。这个习惯帮我省了很多时间，也让我能清楚地看到每个参数修改对结果的影响。