1. 半导体仿真入门为什么选择Silvaco Atlas1.1 从“跑一个PN结”说起如果你正在读微电子、集成电路或者半导体物理相关专业大概率逃不过“器件仿真”这一关。老师布置一个大作业让你画出一个PN结的电场分布或者I-V曲线你打开搜索引擎发现满屏都是Sentaurus、Silvaco、COMSOL这些名字。其中Silvaco Atlas是很多高校和初创公司用得最多的TCAD工具之一原因很简单它有一套完整的生态从工艺仿真Athena到器件仿真Atlas再到可视化TonyPlot全部打通而且网上能找到的教程和示例相对丰富。但问题也恰恰出在这里。很多新手拿到软件之后面对DeckBuild那个黑底白字的界面完全不知道从哪下手。菜单栏里几十个选项示例文件夹里几百个.in文件打开一个看看全是密密麻麻的语句像天书一样。我见过太多人卡在第一步——软件装好了License配好了但就是跑不出第一个能看的曲线。这篇内容就是帮你跨过这道坎的。我会从最基础的安装注意事项讲起然后手把手带你写第一个PN结的仿真输入文件解释每一行代码在干什么最后用TonyPlot把结果画出来。整个过程不需要你事先精通半导体物理只要知道PN结是什么、能看懂基本的能带图就行。目标很明确让你在半天之内从零到一跑通一个完整的半导体器件仿真流程。1.2 Atlas到底能做什么不能做什么在动手之前有必要先搞清楚Atlas的定位。Atlas是一个基于物理的二维/三维器件仿真器它通过求解泊松方程、连续性方程和漂移-扩散模型来计算半导体器件内部的电势、载流子浓度、电流密度等物理量。你可以把它理解成一个“虚拟实验台”给定器件的结构、掺杂分布和电极偏压它告诉你器件会表现出什么样的电学特性。它能做的事情很多PN结二极管、MOSFET、BJT、太阳能电池、LED、功率器件甚至一些新型器件如隧穿场效应晶体管都可以用Atlas来仿真。但它也有明确的边界。Atlas不负责电路级仿真那是SPICE工具的事情它也不直接处理光刻、刻蚀这些工艺步骤的几何形貌变化那是Athena的范畴。所以如果你的目标是“设计一个放大器电路”Atlas不是正确的工具但如果你想知道“这个MOSFET的阈值电压随沟道掺杂浓度怎么变”那Atlas就是你要找的东西。对于新手来说最友好的切入点是PN结。结构简单物理清晰仿真时间短而且结果容易验证——毕竟PN结的I-V特性在教科书上已经讲得很透了。跑通PN结之后你再去看MOSFET或者HBT的示例会发现很多语句是相通的只是结构定义和物理模型更复杂一些。1.3 学习路径与心理预期我建议的学习路径是这样的先跑通一个最小可用的PN结示例理解DeckBuild的基本操作和Atlas的语句结构然后尝试修改掺杂浓度、结深这些参数观察结果的变化接着学习如何添加物理模型如SRH复合、迁移率模型最后再挑战更复杂的器件结构。心理预期方面不要指望第一天就能仿真出和文献一模一样的曲线。半导体仿真涉及大量的物理模型和数值参数默认设置往往只能给出定性正确的结果。定量准确需要你根据具体器件调整模型参数甚至需要校准。但好消息是对于学习目的来说定性正确已经足够让你理解器件的工作原理了。另外Atlas的报错信息有时候比较晦涩一个“Convergence failure”可能对应十几种不同的原因。这时候不要慌先检查结构定义有没有问题再检查掺杂和电极设置最后看物理模型和数值参数。大部分新手遇到的问题都集中在结构定义和电极接触上真正因为物理模型选错而跑不通的情况反而比较少。2. 安装与环境配置少走弯路的几个关键点2.1 安装前的准备工作Silvaco的安装包通常是一个ISO镜像或者压缩包里面包含主程序、示例库和License管理工具。在开始安装之前有几件事需要提前确认。首先是操作系统。Silvaco官方支持Windows和Linux两个平台但根据我的经验Linux版本在稳定性和性能上明显更好尤其是当你需要跑大规模仿真或者使用多核并行的时候。如果你用的是Windows建议至少是Windows 10 64位内存16GB起步32GB会更从容。硬盘空间方面完整安装大约需要20-30GB但如果你只安装Atlas和TonyPlot可以精简到10GB左右。其次是License。Silvaco采用FlexLM授权管理你需要从学校或公司获取License文件。这个文件通常是一个.dat或者.lic文件里面包含了授权的主机ID和功能模块。注意License是和网卡MAC地址绑定的如果你换了网卡或者用了虚拟机MAC地址变了License就会失效。我见过不少人在虚拟机里装好之后因为网络适配器设置变动导致License报错折腾半天才发现是MAC地址的问题。提示安装之前先把License文件放到一个固定的、路径中不含中文和空格的目录下。后续配置环境变量时会反复用到这个路径路径里有空格或者特殊字符容易出问题。2.2 安装步骤与常见报错处理安装过程本身不算复杂双击安装程序一路Next选择安装路径等待进度条走完。但有几个地方容易出问题。第一个是安装路径。默认路径通常是C:\Silvaco或者/usr/local/silvaco建议保持默认不要自己改成带中文或者空格的路径。我试过把安装目录设在“D:\我的软件\Silvaco”下面结果DeckBuild启动时报了一堆找不到库文件的错误改回默认路径就正常了。第二个是环境变量。安装完成后需要设置几个环境变量才能让DeckBuild和Atlas正常工作。在Windows上主要是设置SILVACO变量指向安装根目录然后把%SILVACO%\bin加到PATH里面。在Linux上通常安装脚本会自动在.bashrc或者.cshrc里添加这些设置但如果你用的是自定义shell可能需要手动添加。第三个是License配置。你需要设置LM_LICENSE_FILE或者SILVACO_LICENSE_FILE环境变量指向License文件的完整路径。如果设置正确运行lmstat命令应该能看到License服务器的状态和可用模块。如果报错“Cannot find license file”先检查环境变量有没有生效再检查License文件里的主机名和MAC地址是否和当前机器匹配。常见的一个坑是安装程序自带的License管理工具在某些Windows版本上会闪退。这时候可以尝试用管理员权限运行或者直接手动配置环境变量跳过图形化工具。另外如果你在Linux上安装注意安装脚本可能需要root权限来创建符号链接和设置文件权限但实际使用时不建议用root账户跑仿真容易产生权限混乱。2.3 验证安装是否成功安装完成后怎么确认一切正常最直接的方法是跑一个自带的示例。在安装目录下找到examples文件夹里面按器件类型分了子目录比如diode、mosfet、bjt等。进入diode目录找一个最简单的.in文件比如pn_diode.in然后在DeckBuild里打开它点击运行按钮。如果一切正常你会看到DeckBuild的输出窗口开始滚动日志信息最后显示“Simulation completed successfully”。同时TonyPlot会自动弹出来显示仿真结果。如果报错先看错误信息是什么。常见的错误包括“Cannot open input file”文件路径不对、“Unknown statement”语句拼写错误、“Convergence failure”数值不收敛。前两个通常是配置问题第三个可能是示例本身需要调整参数但自带的示例一般都是调好的所以如果连示例都跑不通大概率是安装或License的问题。另一个验证方法是直接在命令行运行Atlas。在Windows上打开CMD输入atlas -v应该会显示版本信息。在Linux上类似。如果提示“command not found”说明PATH环境变量没设置对。3. DeckBuild与Atlas基础理解输入文件的结构3.1 DeckBuild界面速览第一次打开DeckBuild你可能会觉得这个界面有点“复古”。顶部是菜单栏左边是文件树中间是编辑区底部是输出窗口。菜单栏里的选项很多但新手常用的就那么几个File菜单用来打开和保存文件Commands菜单用来插入常用的Atlas语句Run菜单用来执行仿真Tools菜单里可以启动TonyPlot。编辑区是你写代码的地方。Atlas的输入文件通常以.in为扩展名内容是一系列按顺序执行的语句。这些语句可以分为几大类结构定义定义器件几何形状和掺杂、物理模型指定使用哪些物理模型、数值方法控制求解器的行为、求解命令指定要计算什么比如I-V曲线、结果输出保存哪些数据画哪些图。输出窗口会实时显示仿真过程中的日志信息包括每一步的求解状态、收敛情况、耗时等。如果出错错误信息也会显示在这里。养成看输出窗口的习惯很多问题从日志里就能找到线索。3.2 Atlas输入文件的“三段式”结构一个典型的Atlas输入文件可以分成三个部分结构定义、求解设置、结果输出。这三部分有严格的顺序要求不能颠倒。结构定义部分用mesh、region、electrode、doping这些语句来描述器件的几何形状、材料分布、电极位置和掺杂浓度。这部分相当于“画图纸”告诉Atlas你要仿真的器件长什么样。求解设置部分用models、method、solve这些语句来指定物理模型、数值方法和求解条件。这部分相当于“定规则”告诉Atlas用什么物理方程、用什么数值算法、在什么偏压下求解。结果输出部分用save、log、extract、tonyplot这些语句来保存数据、提取参数、绘制曲线。这部分相当于“出报告”把仿真结果以你需要的格式输出出来。新手最容易犯的错误是把顺序搞乱比如在定义结构之前就写了solve语句或者在求解之后才定义掺杂。Atlas是顺序执行的前面的语句会影响后面的语句所以顺序很重要。3.3 常用语句速查与含义解释下面这张表列出了新手最常用的Atlas语句及其含义你可以把它当作一个速查手册。语句作用常用参数示例mesh定义网格mesh space.mult1.0x.mesh定义X方向网格线x.mesh loc0.0 spacing0.1y.mesh定义Y方向网格线y.mesh loc0.0 spacing0.01region定义材料区域region num1 materialSiliconelectrode定义电极electrode nameanode topdoping定义掺杂doping uniform conc1e16 n.typemodels指定物理模型models srh conmob fldmobmethod指定数值方法method newton trapsolve求解solve init/solve vanode0.1log记录日志log outfpn.logsave保存结构save outfpn.strtonyplot绘图tonyplot pn.log这些语句的具体用法会在后面的实操部分详细展开。现在你只需要有个印象Atlas的输入文件就是由这些语句按顺序排列组成的。4. 第一个PN结仿真从零到一完整实操4.1 器件结构定义网格、区域与电极我们从一个最简单的PN结开始。结构很简单一块硅材料左边是P型掺杂右边是N型掺杂两端各有一个电极。为了简化我们做一维仿真但Atlas是二维仿真器所以我们需要定义一个很薄的二维结构然后在X方向上看结果。先定义网格。网格的疏密直接影响仿真精度和速度。太疏了结果不准太密了跑得慢。对于PN结结附近需要密一些远离结的地方可以疏一些。go atlas mesh space.mult1.0 x.mesh loc0.0 spacing0.1 x.mesh loc1.0 spacing0.1 x.mesh loc2.0 spacing0.1 y.mesh loc0.0 spacing0.01 y.mesh loc0.5 spacing0.01这里X方向从0到2微米Y方向从0到0.5微米。X方向的网格间距是0.1微米Y方向是0.01微米。为什么Y方向要这么密因为我们要在Y方向上加电极如果Y方向网格太疏电极接触会不连续导致电流分布不均匀。实际上对于一维物理问题Y方向只需要一层网格就够了但Atlas要求至少有两层所以设成0.5微米间距0.01微米这样有50个网格点足够了。接下来定义区域和材料。region num1 materialSilicon这一行定义了一个硅区域编号为1。因为整个结构都是硅所以只需要一个区域。然后定义电极。electrode nameanode top electrode namecathode bottomanode在顶部cathode在底部。注意电极的名字可以随便取但后面加偏压的时候要用到。这里我用了anode和cathode你也可以用p和n或者left和right只要前后一致就行。4.2 掺杂分布设置均匀掺杂与高斯掺杂掺杂是PN结的核心。我们做最简单的均匀掺杂左边P型右边N型。doping uniform conc1e17 n.type doping uniform conc1e17 p.type等等这样写不对。doping uniform会把整个区域都掺杂成同一种类型我们需要分区域掺杂。正确的方法是结合region或者x.min、x.max来限定掺杂范围。doping region1 conc1e17 n.type x.min0.0 x.max1.0 doping region1 conc1e17 p.type x.min1.0 x.max2.0这样X从0到1微米是N型掺杂浓度1e17 cm^-3X从1到2微米是P型掺杂浓度也是1e17 cm^-3。结面在X1微米处。如果你想让掺杂更接近实际工艺可以用高斯分布。比如doping gauss conc1e20 n.type x.min0.0 x.max1.0 peak0.0 char0.1这表示在X0处有峰值浓度1e20的N型高斯掺杂特征长度0.1微米。高斯掺杂更接近离子注入后的实际分布但计算量也更大。对于第一个示例均匀掺杂就够了。注意掺杂浓度的单位是cm^-3不是m^-3。这是半导体领域的惯例Atlas也遵循这个惯例。如果你不小心写了1e23那相当于1e17 cm^-3乘以1e6结果会完全不对。4.3 物理模型与数值方法选择物理模型决定了Atlas用什么方程来描述载流子的行为。对于硅基PN结最基本的模型包括srhShockley-Read-Hall复合模型描述通过缺陷能级的复合。conmob浓度依赖的迁移率模型考虑掺杂浓度对迁移率的影响。fldmob场依赖的迁移率模型考虑电场对迁移率的影响。boltz玻尔兹曼统计适用于非简并半导体。models srh conmob fldmob boltz数值方法控制求解器的行为。对于PN结newton方法通常收敛得比较快配合trap可以处理一些数值上的刚性。method newton trapnewton是牛顿迭代法trap是梯形积分法。这两个组合是Atlas中最常用的数值方法。如果遇到收敛问题可以尝试改成method gummel或者调整迭代次数和误差容限。4.4 求解设置与I-V曲线扫描求解分两步先求初始解再扫描偏压。solve initsolve init让Atlas在零偏压下求解一次得到初始的载流子浓度和电势分布。这一步相当于“热身”为后续的偏压扫描提供初始猜测。然后扫描阳极电压log outfpn_diode.log solve vanode0.0 solve vanode0.1 solve vanode0.2 ...手动一行行写太麻烦可以用循环log outfpn_diode.log solve vanode0.0 solve vanode0.1 vstep0.1 vfinal1.0 nameanodesolve vanode0.1 vstep0.1 vfinal1.0 nameanode表示从0.1V开始每次增加0.1V一直加到1.0V。nameanode指定了扫描的电极。log outfpn_diode.log把每次求解的结果记录到日志文件里后面用TonyPlot画I-V曲线。如果你想同时扫描负偏压可以写成solve vanode-1.0 vstep0.1 vfinal1.0 nameanode这样从-1.0V扫到1.0V覆盖了反向偏压和正向偏压。4.5 结果输出与TonyPlot可视化求解完成后保存结构文件和日志文件save outfpn_diode.str tonyplot pn_diode.logsave把最终的结构和物理量保存成.str文件可以用TonyPlot打开查看电场、载流子浓度等分布。tonyplot直接打开日志文件画出I-V曲线。TonyPlot的界面比较直观。打开日志文件后默认会画出阳极电流随阳极电压的变化。你可以通过菜单选择不同的物理量比如电场、电势、电子浓度、空穴浓度等。对于PN结我建议重点看几个图I-V曲线验证整流特性、电场分布看结区电场峰值、载流子浓度分布看耗尽区宽度。提示TonyPlot支持导出图片和数据。如果你要写报告可以直接把图导出成PNG或者EPS格式。数据可以导出成CSV方便用其他工具进一步处理。5. 常见问题与排查技巧实录5.1 收敛失败最常见也最头疼的问题“Convergence failure”是Atlas新手遇到最多的报错。它的原因很多但排查思路是有章可循的。首先检查结构定义。如果网格太疏尤其是结附近网格太疏会导致电场变化剧烈求解器难以收敛。解决办法是在结附近加密网格。比如在X1微米附近加几条网格线x.mesh loc0.9 spacing0.01 x.mesh loc1.0 spacing0.01 x.mesh loc1.1 spacing0.01其次检查掺杂。如果掺杂浓度突变太剧烈比如从1e20突然跳到1e15也会导致收敛困难。可以尝试用doping gauss代替doping uniform让掺杂有一个平滑的过渡。然后检查物理模型。有些模型在特定条件下会导致方程刚性增强比如fldmob在高电场下。如果收敛不了可以先去掉fldmob只保留srh和conmob跑通了再加回来。最后检查数值方法。newton方法对初始猜测比较敏感如果初始解不好容易发散。可以尝试先用gummel方法跑几步再用newton。或者调整method的参数method newton trap itlimit50itlimit50把每次求解的最大迭代次数增加到50给求解器更多机会收敛。5.2 结果异常曲线不对、电流为零有时候仿真跑通了但结果明显不对。比如I-V曲线是一条直线或者电流始终为零。这种情况通常不是数值问题而是设置问题。电流为零最常见的原因是电极没接触好。检查electrode语句确保电极定义在正确的边界上。如果电极定义在top但你的结构在Y方向只有一层网格电极可能没有真正接触到半导体。解决办法是确保Y方向至少有两层网格并且电极定义在最外层。另一个原因是掺杂类型搞反了。PN结的整流特性依赖于P型和N型的正确分布。如果两边都是N型那就是一个电阻I-V曲线当然是直线。检查doping语句里的n.type和p.type有没有写错。还有一种情况是偏压加错了电极。solve vanode0.1 nameanode表示在anode上加0.1Vcathode接地。如果你把name写成了cathode那相当于在cathode上加正压anode接地结果正好反过来。5.3 性能优化让仿真跑得更快Atlas的仿真速度取决于网格数量、物理模型复杂度和求解步长。对于PN结这种简单器件默认设置通常几秒钟就能跑完。但如果你做的是三维仿真或者复杂器件速度可能就成了问题。减少网格数量是最直接的办法。在保证精度的前提下尽量用稀疏的网格。对于PN结结附近密一些远离结的地方可以很疏。我通常会在结附近用0.01微米的间距远离结的地方用0.1微米甚至0.5微米。简化物理模型也能提速。如果只是看定性趋势可以去掉fldmob和srh只保留最基本的漂移-扩散模型。当然结果的精度会下降但对于理解器件原理来说足够了。调整求解步长也很重要。vstep0.1比vstep0.01快十倍但曲线会粗糙一些。如果只是看趋势0.1V的步长足够了。如果需要精确提取某个参数比如开启电压可以先用大步长扫一遍找到大致范围后再用小步长细扫。5.4 常见问题速查表问题现象可能原因解决办法Convergence failure网格太疏、掺杂突变、模型太复杂加密结区网格、用高斯掺杂、简化模型电流始终为零电极未接触、掺杂类型错误检查electrode和doping语句I-V曲线是直线两边同型掺杂、偏压加错电极检查掺杂类型和solve的name参数仿真速度慢网格太密、模型太多、步长太小减少网格、简化模型、增大步长TonyPlot打不开文件路径错误、文件损坏检查路径、重新运行仿真License报错环境变量未设置、MAC地址变化检查LM_LICENSE_FILE、重新获取License6. 从PN结到更复杂的器件下一步怎么走6.1 参数扫描与结果分析跑通第一个PN结之后下一步是学会做参数扫描。比如改变掺杂浓度看I-V曲线怎么变改变结深看电场分布怎么变。Atlas支持用循环语句做参数扫描但更简单的方法是用DeckBuild的“参数化仿真”功能或者手动修改输入文件然后批量运行。我个人的习惯是写一个模板文件把需要扫描的参数用变量代替然后用脚本批量生成和运行。比如for conc in 1e16 1e17 1e18; do sed s/CONC/$conc/g pn_template.in pn_$conc.in deckbuild -run pn_$conc.in done这样可以在命令行批量跑多个仿真省去手动修改的麻烦。当然前提是你已经配置好了命令行环境。结果分析方面除了看I-V曲线还可以提取一些关键参数比如开启电压、反向饱和电流、理想因子等。Atlas提供了extract语句可以自动提取这些参数extract nameVturnon x.val from curve(v.anode, i.anode) where y.val1e-6这行语句提取阳极电流达到1e-6 A时的阳极电压也就是开启电压。类似的提取语句还有很多具体用法可以参考Atlas手册。6.2 学习资源与进阶方向Atlas的学习资源主要有三个来源官方手册、示例库和社区论坛。官方手册是最权威的但内容很多新手容易迷失。我的建议是先看“Atlas Users Manual”的前几章了解基本概念和语句结构然后直接去示例库里找对应的例子。示例库里的文件都是调好的跑通之后修改参数观察结果变化这是最快的学习方式。进阶方向取决于你的具体需求。如果你做的是功率器件需要学习碰撞电离模型和热效应仿真如果你做的是光电器件需要学习光吸收和光电流模型如果你做的是射频器件需要学习小信号AC仿真。每个方向都有对应的物理模型和求解设置但基础的结构定义和求解流程是相通的。6.3 一些个人体会最后分享几个我在使用Atlas过程中总结的小经验。第一养成备份输入文件的习惯。每次修改之前先另存一份尤其是当仿真跑通之后。我吃过亏调参数调了半天结果不小心把跑通的版本覆盖了只能从头再来。第二善用TonyPlot的“比较”功能。TonyPlot可以同时打开多个日志文件把不同参数下的I-V曲线画在同一张图上。这对于分析参数影响非常有用。第三不要迷信默认参数。Atlas的默认设置是为了通用性不一定适合你的具体器件。多尝试不同的模型组合和数值参数找到最适合你问题的配置。第四遇到问题先看日志。DeckBuild的输出窗口会显示详细的求解过程包括每一步的迭代次数、误差、收敛状态。很多问题从日志里就能看出端倪比盲目试错高效得多。第五保持耐心。半导体仿真是一个需要积累的过程刚开始可能会被各种报错折磨但跑通几个例子之后你会发现很多问题是重复的解决起来也越来越快。我当初学Atlas的时候第一个PN结跑了整整两天才跑通但现在回头看那两天学到的东西比后面两个月都多。 Silvaco Atlas半导体仿真入门：从零跑通第一个PN结
1. 半导体仿真入门为什么选择Silvaco Atlas1.1 从“跑一个PN结”说起如果你正在读微电子、集成电路或者半导体物理相关专业大概率逃不过“器件仿真”这一关。老师布置一个大作业让你画出一个PN结的电场分布或者I-V曲线你打开搜索引擎发现满屏都是Sentaurus、Silvaco、COMSOL这些名字。其中Silvaco Atlas是很多高校和初创公司用得最多的TCAD工具之一原因很简单它有一套完整的生态从工艺仿真Athena到器件仿真Atlas再到可视化TonyPlot全部打通而且网上能找到的教程和示例相对丰富。但问题也恰恰出在这里。很多新手拿到软件之后面对DeckBuild那个黑底白字的界面完全不知道从哪下手。菜单栏里几十个选项示例文件夹里几百个.in文件打开一个看看全是密密麻麻的语句像天书一样。我见过太多人卡在第一步——软件装好了License配好了但就是跑不出第一个能看的曲线。这篇内容就是帮你跨过这道坎的。我会从最基础的安装注意事项讲起然后手把手带你写第一个PN结的仿真输入文件解释每一行代码在干什么最后用TonyPlot把结果画出来。整个过程不需要你事先精通半导体物理只要知道PN结是什么、能看懂基本的能带图就行。目标很明确让你在半天之内从零到一跑通一个完整的半导体器件仿真流程。1.2 Atlas到底能做什么不能做什么在动手之前有必要先搞清楚Atlas的定位。Atlas是一个基于物理的二维/三维器件仿真器它通过求解泊松方程、连续性方程和漂移-扩散模型来计算半导体器件内部的电势、载流子浓度、电流密度等物理量。你可以把它理解成一个“虚拟实验台”给定器件的结构、掺杂分布和电极偏压它告诉你器件会表现出什么样的电学特性。它能做的事情很多PN结二极管、MOSFET、BJT、太阳能电池、LED、功率器件甚至一些新型器件如隧穿场效应晶体管都可以用Atlas来仿真。但它也有明确的边界。Atlas不负责电路级仿真那是SPICE工具的事情它也不直接处理光刻、刻蚀这些工艺步骤的几何形貌变化那是Athena的范畴。所以如果你的目标是“设计一个放大器电路”Atlas不是正确的工具但如果你想知道“这个MOSFET的阈值电压随沟道掺杂浓度怎么变”那Atlas就是你要找的东西。对于新手来说最友好的切入点是PN结。结构简单物理清晰仿真时间短而且结果容易验证——毕竟PN结的I-V特性在教科书上已经讲得很透了。跑通PN结之后你再去看MOSFET或者HBT的示例会发现很多语句是相通的只是结构定义和物理模型更复杂一些。1.3 学习路径与心理预期我建议的学习路径是这样的先跑通一个最小可用的PN结示例理解DeckBuild的基本操作和Atlas的语句结构然后尝试修改掺杂浓度、结深这些参数观察结果的变化接着学习如何添加物理模型如SRH复合、迁移率模型最后再挑战更复杂的器件结构。心理预期方面不要指望第一天就能仿真出和文献一模一样的曲线。半导体仿真涉及大量的物理模型和数值参数默认设置往往只能给出定性正确的结果。定量准确需要你根据具体器件调整模型参数甚至需要校准。但好消息是对于学习目的来说定性正确已经足够让你理解器件的工作原理了。另外Atlas的报错信息有时候比较晦涩一个“Convergence failure”可能对应十几种不同的原因。这时候不要慌先检查结构定义有没有问题再检查掺杂和电极设置最后看物理模型和数值参数。大部分新手遇到的问题都集中在结构定义和电极接触上真正因为物理模型选错而跑不通的情况反而比较少。2. 安装与环境配置少走弯路的几个关键点2.1 安装前的准备工作Silvaco的安装包通常是一个ISO镜像或者压缩包里面包含主程序、示例库和License管理工具。在开始安装之前有几件事需要提前确认。首先是操作系统。Silvaco官方支持Windows和Linux两个平台但根据我的经验Linux版本在稳定性和性能上明显更好尤其是当你需要跑大规模仿真或者使用多核并行的时候。如果你用的是Windows建议至少是Windows 10 64位内存16GB起步32GB会更从容。硬盘空间方面完整安装大约需要20-30GB但如果你只安装Atlas和TonyPlot可以精简到10GB左右。其次是License。Silvaco采用FlexLM授权管理你需要从学校或公司获取License文件。这个文件通常是一个.dat或者.lic文件里面包含了授权的主机ID和功能模块。注意License是和网卡MAC地址绑定的如果你换了网卡或者用了虚拟机MAC地址变了License就会失效。我见过不少人在虚拟机里装好之后因为网络适配器设置变动导致License报错折腾半天才发现是MAC地址的问题。提示安装之前先把License文件放到一个固定的、路径中不含中文和空格的目录下。后续配置环境变量时会反复用到这个路径路径里有空格或者特殊字符容易出问题。2.2 安装步骤与常见报错处理安装过程本身不算复杂双击安装程序一路Next选择安装路径等待进度条走完。但有几个地方容易出问题。第一个是安装路径。默认路径通常是C:\Silvaco或者/usr/local/silvaco建议保持默认不要自己改成带中文或者空格的路径。我试过把安装目录设在“D:\我的软件\Silvaco”下面结果DeckBuild启动时报了一堆找不到库文件的错误改回默认路径就正常了。第二个是环境变量。安装完成后需要设置几个环境变量才能让DeckBuild和Atlas正常工作。在Windows上主要是设置SILVACO变量指向安装根目录然后把%SILVACO%\bin加到PATH里面。在Linux上通常安装脚本会自动在.bashrc或者.cshrc里添加这些设置但如果你用的是自定义shell可能需要手动添加。第三个是License配置。你需要设置LM_LICENSE_FILE或者SILVACO_LICENSE_FILE环境变量指向License文件的完整路径。如果设置正确运行lmstat命令应该能看到License服务器的状态和可用模块。如果报错“Cannot find license file”先检查环境变量有没有生效再检查License文件里的主机名和MAC地址是否和当前机器匹配。常见的一个坑是安装程序自带的License管理工具在某些Windows版本上会闪退。这时候可以尝试用管理员权限运行或者直接手动配置环境变量跳过图形化工具。另外如果你在Linux上安装注意安装脚本可能需要root权限来创建符号链接和设置文件权限但实际使用时不建议用root账户跑仿真容易产生权限混乱。2.3 验证安装是否成功安装完成后怎么确认一切正常最直接的方法是跑一个自带的示例。在安装目录下找到examples文件夹里面按器件类型分了子目录比如diode、mosfet、bjt等。进入diode目录找一个最简单的.in文件比如pn_diode.in然后在DeckBuild里打开它点击运行按钮。如果一切正常你会看到DeckBuild的输出窗口开始滚动日志信息最后显示“Simulation completed successfully”。同时TonyPlot会自动弹出来显示仿真结果。如果报错先看错误信息是什么。常见的错误包括“Cannot open input file”文件路径不对、“Unknown statement”语句拼写错误、“Convergence failure”数值不收敛。前两个通常是配置问题第三个可能是示例本身需要调整参数但自带的示例一般都是调好的所以如果连示例都跑不通大概率是安装或License的问题。另一个验证方法是直接在命令行运行Atlas。在Windows上打开CMD输入atlas -v应该会显示版本信息。在Linux上类似。如果提示“command not found”说明PATH环境变量没设置对。3. DeckBuild与Atlas基础理解输入文件的结构3.1 DeckBuild界面速览第一次打开DeckBuild你可能会觉得这个界面有点“复古”。顶部是菜单栏左边是文件树中间是编辑区底部是输出窗口。菜单栏里的选项很多但新手常用的就那么几个File菜单用来打开和保存文件Commands菜单用来插入常用的Atlas语句Run菜单用来执行仿真Tools菜单里可以启动TonyPlot。编辑区是你写代码的地方。Atlas的输入文件通常以.in为扩展名内容是一系列按顺序执行的语句。这些语句可以分为几大类结构定义定义器件几何形状和掺杂、物理模型指定使用哪些物理模型、数值方法控制求解器的行为、求解命令指定要计算什么比如I-V曲线、结果输出保存哪些数据画哪些图。输出窗口会实时显示仿真过程中的日志信息包括每一步的求解状态、收敛情况、耗时等。如果出错错误信息也会显示在这里。养成看输出窗口的习惯很多问题从日志里就能找到线索。3.2 Atlas输入文件的“三段式”结构一个典型的Atlas输入文件可以分成三个部分结构定义、求解设置、结果输出。这三部分有严格的顺序要求不能颠倒。结构定义部分用mesh、region、electrode、doping这些语句来描述器件的几何形状、材料分布、电极位置和掺杂浓度。这部分相当于“画图纸”告诉Atlas你要仿真的器件长什么样。求解设置部分用models、method、solve这些语句来指定物理模型、数值方法和求解条件。这部分相当于“定规则”告诉Atlas用什么物理方程、用什么数值算法、在什么偏压下求解。结果输出部分用save、log、extract、tonyplot这些语句来保存数据、提取参数、绘制曲线。这部分相当于“出报告”把仿真结果以你需要的格式输出出来。新手最容易犯的错误是把顺序搞乱比如在定义结构之前就写了solve语句或者在求解之后才定义掺杂。Atlas是顺序执行的前面的语句会影响后面的语句所以顺序很重要。3.3 常用语句速查与含义解释下面这张表列出了新手最常用的Atlas语句及其含义你可以把它当作一个速查手册。语句作用常用参数示例mesh定义网格mesh space.mult1.0x.mesh定义X方向网格线x.mesh loc0.0 spacing0.1y.mesh定义Y方向网格线y.mesh loc0.0 spacing0.01region定义材料区域region num1 materialSiliconelectrode定义电极electrode nameanode topdoping定义掺杂doping uniform conc1e16 n.typemodels指定物理模型models srh conmob fldmobmethod指定数值方法method newton trapsolve求解solve init/solve vanode0.1log记录日志log outfpn.logsave保存结构save outfpn.strtonyplot绘图tonyplot pn.log这些语句的具体用法会在后面的实操部分详细展开。现在你只需要有个印象Atlas的输入文件就是由这些语句按顺序排列组成的。4. 第一个PN结仿真从零到一完整实操4.1 器件结构定义网格、区域与电极我们从一个最简单的PN结开始。结构很简单一块硅材料左边是P型掺杂右边是N型掺杂两端各有一个电极。为了简化我们做一维仿真但Atlas是二维仿真器所以我们需要定义一个很薄的二维结构然后在X方向上看结果。先定义网格。网格的疏密直接影响仿真精度和速度。太疏了结果不准太密了跑得慢。对于PN结结附近需要密一些远离结的地方可以疏一些。go atlas mesh space.mult1.0 x.mesh loc0.0 spacing0.1 x.mesh loc1.0 spacing0.1 x.mesh loc2.0 spacing0.1 y.mesh loc0.0 spacing0.01 y.mesh loc0.5 spacing0.01这里X方向从0到2微米Y方向从0到0.5微米。X方向的网格间距是0.1微米Y方向是0.01微米。为什么Y方向要这么密因为我们要在Y方向上加电极如果Y方向网格太疏电极接触会不连续导致电流分布不均匀。实际上对于一维物理问题Y方向只需要一层网格就够了但Atlas要求至少有两层所以设成0.5微米间距0.01微米这样有50个网格点足够了。接下来定义区域和材料。region num1 materialSilicon这一行定义了一个硅区域编号为1。因为整个结构都是硅所以只需要一个区域。然后定义电极。electrode nameanode top electrode namecathode bottomanode在顶部cathode在底部。注意电极的名字可以随便取但后面加偏压的时候要用到。这里我用了anode和cathode你也可以用p和n或者left和right只要前后一致就行。4.2 掺杂分布设置均匀掺杂与高斯掺杂掺杂是PN结的核心。我们做最简单的均匀掺杂左边P型右边N型。doping uniform conc1e17 n.type doping uniform conc1e17 p.type等等这样写不对。doping uniform会把整个区域都掺杂成同一种类型我们需要分区域掺杂。正确的方法是结合region或者x.min、x.max来限定掺杂范围。doping region1 conc1e17 n.type x.min0.0 x.max1.0 doping region1 conc1e17 p.type x.min1.0 x.max2.0这样X从0到1微米是N型掺杂浓度1e17 cm^-3X从1到2微米是P型掺杂浓度也是1e17 cm^-3。结面在X1微米处。如果你想让掺杂更接近实际工艺可以用高斯分布。比如doping gauss conc1e20 n.type x.min0.0 x.max1.0 peak0.0 char0.1这表示在X0处有峰值浓度1e20的N型高斯掺杂特征长度0.1微米。高斯掺杂更接近离子注入后的实际分布但计算量也更大。对于第一个示例均匀掺杂就够了。注意掺杂浓度的单位是cm^-3不是m^-3。这是半导体领域的惯例Atlas也遵循这个惯例。如果你不小心写了1e23那相当于1e17 cm^-3乘以1e6结果会完全不对。4.3 物理模型与数值方法选择物理模型决定了Atlas用什么方程来描述载流子的行为。对于硅基PN结最基本的模型包括srhShockley-Read-Hall复合模型描述通过缺陷能级的复合。conmob浓度依赖的迁移率模型考虑掺杂浓度对迁移率的影响。fldmob场依赖的迁移率模型考虑电场对迁移率的影响。boltz玻尔兹曼统计适用于非简并半导体。models srh conmob fldmob boltz数值方法控制求解器的行为。对于PN结newton方法通常收敛得比较快配合trap可以处理一些数值上的刚性。method newton trapnewton是牛顿迭代法trap是梯形积分法。这两个组合是Atlas中最常用的数值方法。如果遇到收敛问题可以尝试改成method gummel或者调整迭代次数和误差容限。4.4 求解设置与I-V曲线扫描求解分两步先求初始解再扫描偏压。solve initsolve init让Atlas在零偏压下求解一次得到初始的载流子浓度和电势分布。这一步相当于“热身”为后续的偏压扫描提供初始猜测。然后扫描阳极电压log outfpn_diode.log solve vanode0.0 solve vanode0.1 solve vanode0.2 ...手动一行行写太麻烦可以用循环log outfpn_diode.log solve vanode0.0 solve vanode0.1 vstep0.1 vfinal1.0 nameanodesolve vanode0.1 vstep0.1 vfinal1.0 nameanode表示从0.1V开始每次增加0.1V一直加到1.0V。nameanode指定了扫描的电极。log outfpn_diode.log把每次求解的结果记录到日志文件里后面用TonyPlot画I-V曲线。如果你想同时扫描负偏压可以写成solve vanode-1.0 vstep0.1 vfinal1.0 nameanode这样从-1.0V扫到1.0V覆盖了反向偏压和正向偏压。4.5 结果输出与TonyPlot可视化求解完成后保存结构文件和日志文件save outfpn_diode.str tonyplot pn_diode.logsave把最终的结构和物理量保存成.str文件可以用TonyPlot打开查看电场、载流子浓度等分布。tonyplot直接打开日志文件画出I-V曲线。TonyPlot的界面比较直观。打开日志文件后默认会画出阳极电流随阳极电压的变化。你可以通过菜单选择不同的物理量比如电场、电势、电子浓度、空穴浓度等。对于PN结我建议重点看几个图I-V曲线验证整流特性、电场分布看结区电场峰值、载流子浓度分布看耗尽区宽度。提示TonyPlot支持导出图片和数据。如果你要写报告可以直接把图导出成PNG或者EPS格式。数据可以导出成CSV方便用其他工具进一步处理。5. 常见问题与排查技巧实录5.1 收敛失败最常见也最头疼的问题“Convergence failure”是Atlas新手遇到最多的报错。它的原因很多但排查思路是有章可循的。首先检查结构定义。如果网格太疏尤其是结附近网格太疏会导致电场变化剧烈求解器难以收敛。解决办法是在结附近加密网格。比如在X1微米附近加几条网格线x.mesh loc0.9 spacing0.01 x.mesh loc1.0 spacing0.01 x.mesh loc1.1 spacing0.01其次检查掺杂。如果掺杂浓度突变太剧烈比如从1e20突然跳到1e15也会导致收敛困难。可以尝试用doping gauss代替doping uniform让掺杂有一个平滑的过渡。然后检查物理模型。有些模型在特定条件下会导致方程刚性增强比如fldmob在高电场下。如果收敛不了可以先去掉fldmob只保留srh和conmob跑通了再加回来。最后检查数值方法。newton方法对初始猜测比较敏感如果初始解不好容易发散。可以尝试先用gummel方法跑几步再用newton。或者调整method的参数method newton trap itlimit50itlimit50把每次求解的最大迭代次数增加到50给求解器更多机会收敛。5.2 结果异常曲线不对、电流为零有时候仿真跑通了但结果明显不对。比如I-V曲线是一条直线或者电流始终为零。这种情况通常不是数值问题而是设置问题。电流为零最常见的原因是电极没接触好。检查electrode语句确保电极定义在正确的边界上。如果电极定义在top但你的结构在Y方向只有一层网格电极可能没有真正接触到半导体。解决办法是确保Y方向至少有两层网格并且电极定义在最外层。另一个原因是掺杂类型搞反了。PN结的整流特性依赖于P型和N型的正确分布。如果两边都是N型那就是一个电阻I-V曲线当然是直线。检查doping语句里的n.type和p.type有没有写错。还有一种情况是偏压加错了电极。solve vanode0.1 nameanode表示在anode上加0.1Vcathode接地。如果你把name写成了cathode那相当于在cathode上加正压anode接地结果正好反过来。5.3 性能优化让仿真跑得更快Atlas的仿真速度取决于网格数量、物理模型复杂度和求解步长。对于PN结这种简单器件默认设置通常几秒钟就能跑完。但如果你做的是三维仿真或者复杂器件速度可能就成了问题。减少网格数量是最直接的办法。在保证精度的前提下尽量用稀疏的网格。对于PN结结附近密一些远离结的地方可以很疏。我通常会在结附近用0.01微米的间距远离结的地方用0.1微米甚至0.5微米。简化物理模型也能提速。如果只是看定性趋势可以去掉fldmob和srh只保留最基本的漂移-扩散模型。当然结果的精度会下降但对于理解器件原理来说足够了。调整求解步长也很重要。vstep0.1比vstep0.01快十倍但曲线会粗糙一些。如果只是看趋势0.1V的步长足够了。如果需要精确提取某个参数比如开启电压可以先用大步长扫一遍找到大致范围后再用小步长细扫。5.4 常见问题速查表问题现象可能原因解决办法Convergence failure网格太疏、掺杂突变、模型太复杂加密结区网格、用高斯掺杂、简化模型电流始终为零电极未接触、掺杂类型错误检查electrode和doping语句I-V曲线是直线两边同型掺杂、偏压加错电极检查掺杂类型和solve的name参数仿真速度慢网格太密、模型太多、步长太小减少网格、简化模型、增大步长TonyPlot打不开文件路径错误、文件损坏检查路径、重新运行仿真License报错环境变量未设置、MAC地址变化检查LM_LICENSE_FILE、重新获取License6. 从PN结到更复杂的器件下一步怎么走6.1 参数扫描与结果分析跑通第一个PN结之后下一步是学会做参数扫描。比如改变掺杂浓度看I-V曲线怎么变改变结深看电场分布怎么变。Atlas支持用循环语句做参数扫描但更简单的方法是用DeckBuild的“参数化仿真”功能或者手动修改输入文件然后批量运行。我个人的习惯是写一个模板文件把需要扫描的参数用变量代替然后用脚本批量生成和运行。比如for conc in 1e16 1e17 1e18; do sed s/CONC/$conc/g pn_template.in pn_$conc.in deckbuild -run pn_$conc.in done这样可以在命令行批量跑多个仿真省去手动修改的麻烦。当然前提是你已经配置好了命令行环境。结果分析方面除了看I-V曲线还可以提取一些关键参数比如开启电压、反向饱和电流、理想因子等。Atlas提供了extract语句可以自动提取这些参数extract nameVturnon x.val from curve(v.anode, i.anode) where y.val1e-6这行语句提取阳极电流达到1e-6 A时的阳极电压也就是开启电压。类似的提取语句还有很多具体用法可以参考Atlas手册。6.2 学习资源与进阶方向Atlas的学习资源主要有三个来源官方手册、示例库和社区论坛。官方手册是最权威的但内容很多新手容易迷失。我的建议是先看“Atlas Users Manual”的前几章了解基本概念和语句结构然后直接去示例库里找对应的例子。示例库里的文件都是调好的跑通之后修改参数观察结果变化这是最快的学习方式。进阶方向取决于你的具体需求。如果你做的是功率器件需要学习碰撞电离模型和热效应仿真如果你做的是光电器件需要学习光吸收和光电流模型如果你做的是射频器件需要学习小信号AC仿真。每个方向都有对应的物理模型和求解设置但基础的结构定义和求解流程是相通的。6.3 一些个人体会最后分享几个我在使用Atlas过程中总结的小经验。第一养成备份输入文件的习惯。每次修改之前先另存一份尤其是当仿真跑通之后。我吃过亏调参数调了半天结果不小心把跑通的版本覆盖了只能从头再来。第二善用TonyPlot的“比较”功能。TonyPlot可以同时打开多个日志文件把不同参数下的I-V曲线画在同一张图上。这对于分析参数影响非常有用。第三不要迷信默认参数。Atlas的默认设置是为了通用性不一定适合你的具体器件。多尝试不同的模型组合和数值参数找到最适合你问题的配置。第四遇到问题先看日志。DeckBuild的输出窗口会显示详细的求解过程包括每一步的迭代次数、误差、收敛状态。很多问题从日志里就能看出端倪比盲目试错高效得多。第五保持耐心。半导体仿真是一个需要积累的过程刚开始可能会被各种报错折磨但跑通几个例子之后你会发现很多问题是重复的解决起来也越来越快。我当初学Atlas的时候第一个PN结跑了整整两天才跑通但现在回头看那两天学到的东西比后面两个月都多。