1. 从零上手 Silvaco Atlas材料参数与模型定义到底在做什么刚接触 Silvaco Atlas 的人十有八九会在第二个星期卡在同一个地方前面的网格和电极都建好了代码跑起来却直接报错或者结果明显不对。问题往往不在语法而在于材料参数和物理模型没有定义清楚。Atlas 是一个基于物理的半导体器件仿真工具它不像某些高层仿真器那样给你一堆默认值就能跑它要求你明确告诉它这个器件用什么材料、材料在什么条件下表现如何、载流子怎么运动、复合和产生遵循什么规律。这些信息全部通过材料参数和模型语句来传递。我刚开始用 Atlas 的时候习惯性地把材料定义当成“填表”——把硅的禁带宽度、迁移率、介电常数填进去就完事。结果第一次跑 MOSFET 的 I-V 曲线阈值电压偏了将近 0.5V漏电流大了两个数量级。后来逐行排查才发现问题出在迁移率模型没有匹配掺杂浓度范围以及没有启用与浓度相关的复合模型。这件事让我意识到材料参数和模型定义不是可选项而是决定仿真结果是否可信的核心环节。这篇文章面向的是已经跑通过 Atlas 基础示例、能独立搭建简单器件结构、但对材料参数和物理模型还停留在“照抄示例”阶段的读者。我会把材料参数的定义逻辑、模型的选择依据、参数之间的耦合关系以及实际调试中踩过的坑按照我自己的操作顺序完整梳理一遍。你不需要有半导体物理的博士学位但需要知道什么是载流子、什么是掺杂、什么是能带。如果这些概念还比较模糊建议先补一下基础否则后面的参数调整会变成盲人摸象。Atlas 中的材料参数和模型定义本质上是在回答三个问题第一器件由哪些材料构成每种材料的基本电学和热学属性是什么第二在给定的偏置和温度下载流子如何输运、如何产生和复合第三这些行为如何随位置、浓度、电场和温度变化。把这三个问题回答清楚仿真结果才有物理意义。接下来的内容会围绕这三个问题展开每一步都给出我实际使用的语句和参数值并解释为什么这样选。2. 材料参数定义从基础属性到高级参数2.1 材料定义的基本语句结构Atlas 中定义材料参数的标准方式是使用material语句。它的基本结构是material materialSilicon parametervalue ...或者针对特定区域material region2 materialSilicon parametervalue ...前者修改的是全局材料属性后者只修改指定区域的属性。我通常的做法是先用全局语句设定材料的基础参数再针对关键区域如沟道、漂移区用区域语句做局部调整。这样做的好处是代码结构清晰后续修改时不容易误伤其他区域。一个典型的硅材料定义片段如下material materialSilicon EG3001.12 PERMITTIVITY11.7 \ MUN1400 MUP450 NC3002.8e19 NV3001.04e19 \ AFFINITY4.05这里每个参数都有明确的物理含义。EG300是 300K 时的禁带宽度单位是电子伏特。PERMITTIVITY是相对介电常数。MUN和MUP分别是低场电子和空穴迁移率单位是 cm²/(V·s)。NC300和NV300是导带和价带的有效状态密度。AFFINITY是电子亲和能。这些参数共同决定了材料在平衡状态下的基本行为。注意Atlas 中参数的单位体系是 cm、V、A、s 的组合迁移率用 cm²/(V·s)浓度用 cm⁻³电场用 V/cm。单位搞错是新手最常见的错误之一而且往往不会直接报错只会让结果悄悄偏掉。2.2 禁带宽度与温度的关系禁带宽度不是常数它随温度变化。Atlas 提供了两种方式来处理这种依赖关系。第一种是直接给出不同温度下的值第二种是使用内置的温度模型。我一般用EGALPHA和EGBETA两个参数来定义 Varshni 关系material materialSilicon EG3001.12 EGALPHA4.73e-4 EGBETA636Varshni 公式的形式是Eg(T) Eg(0) - alpha * T² / (T beta)其中Eg(0)是 0K 时的禁带宽度alpha和beta是材料常数。对于硅alpha约为 4.73×10⁻⁴ eV/Kbeta约为 636K。这个关系在温度仿真中非常重要尤其是做功率器件或高温应用时。我曾经在一个 SiC 器件的仿真中忽略了温度对禁带宽度的影响结果在 200°C 时计算出的本征载流子浓度比实际低了近一个数量级导致漏电流严重低估。2.3 迁移率模型与掺杂浓度的耦合迁移率是材料参数中最容易出问题的地方。Atlas 默认使用常数迁移率但实际材料中迁移率会随掺杂浓度、电场、温度变化。我通常至少启用与掺杂浓度相关的迁移率模型models materialSilicon CONMOBCONMOB模型使用 Masetti 等人提出的经验公式描述迁移率随掺杂浓度的下降。对于硅电子和空穴的迁移率在轻掺杂时分别约为 1400 和 450 cm²/(V·s)在重掺杂时会降到 100 以下。如果你做的是 MOSFET 的沟道仿真沟道掺杂通常在 10¹⁷ cm⁻³ 量级此时迁移率已经明显低于低场值。不启用这个模型计算出的驱动电流会偏大。更高阶的迁移率模型包括FLDMOB场相关迁移率、TMPMOB温度相关迁移率和SURFMOB表面迁移率。我的建议是先启用CONMOB和FLDMOB这两个覆盖了大多数情况。如果做的是深亚微米器件再考虑SURFMOB因为它考虑了表面粗糙度散射和库仑散射。models materialSilicon CONMOB FLDMOB实操心得启用多个迁移率模型时Atlas 会按照内置的优先级进行组合。不同模型之间的参数可能会相互影响建议每次只增加一个模型观察结果变化确认合理后再继续叠加。2.4 有效状态密度与温度有效状态密度NC300和NV300同样随温度变化关系是NC(T) NC300 * (T/300)^1.5Atlas 默认使用这个关系不需要额外定义。但如果你有实验数据表明材料在特定温度下偏离这个关系可以用NCTAU等参数做修正。对于大多数硅基器件默认值已经足够。对于化合物半导体比如 GaAs 或 GaN有效状态密度与温度的关系可能更复杂需要查阅文献确认。2.5 其他常用材料参数除了上述核心参数还有一些参数在特定场景下必须定义。比如做光电器件时需要定义折射率和吸收系数material materialSilicon REFRACTIVE3.5 \ ABSORPTION1e4做热仿真时需要定义热导率和热容material materialSilicon TC1.5 RHOCP2.3TC是热导率单位是 W/(cm·K)。RHOCP是体积热容单位是 J/(cm³·K)。这两个参数在电热耦合仿真中必不可少。我曾经在一个功率器件的仿真中只定义了电学参数忘了热导率结果 Atlas 使用了默认值导致结温计算完全不可信。3. 物理模型选择从基础输运到高级效应3.1 载流子输运模型Atlas 提供了多个层次的载流子输运模型从最简单的漂移-扩散模型到复杂的蒙特卡洛模型。对于大多数器件仿真漂移-扩散模型已经足够。它是 Atlas 的默认模型不需要显式启用。但如果你做的是亚微米或纳米尺度器件可能需要考虑能量平衡模型models materialSilicon ENERGYBALANCE能量平衡模型引入了载流子温度作为独立变量能够捕捉非局部输运效应比如速度过冲。代价是计算量显著增加收敛难度也更大。我的经验是沟道长度大于 0.5 微米时漂移-扩散模型足够小于 0.25 微米时能量平衡模型开始显现优势小于 0.1 微米时可能需要考虑蒙特卡洛方法但那已经超出 Atlas 的常规使用范围。3.2 产生-复合模型产生-复合模型决定了载流子的寿命和复合速率。最基础的是 Shockley-Read-Hall 复合models materialSilicon SRHSRH 复合通过缺陷能级捕获载流子需要定义载流子寿命TAUN和TAUPmaterial materialSilicon TAUN1e-7 TAUP1e-7寿命值对结果影响很大。硅中电子和空穴的寿命通常在微秒到毫秒量级但重掺杂区域会显著降低。我通常根据掺杂浓度设置寿命material materialSilicon TAUN1e-7 TAUP1e-7 \ NSRHN1e17 NSRHP1e17NSRHN和NSRHP是 SRH 复合的浓度参考值用于描述寿命随掺杂浓度的下降。除了 SRH 复合还有辐射复合和俄歇复合。辐射复合在硅中很弱通常可以忽略。俄歇复合在重掺杂区域很重要models materialSilicon AUGER俄歇复合需要定义AUGN和AUGP系数硅中典型值约为 2.8×10⁻³¹ cm⁶/s 和 9.9×10⁻³² cm⁶/s。在掺杂浓度超过 10¹⁸ cm⁻³ 时俄歇复合会成为主要复合机制。3.3 碰撞电离模型做击穿仿真时碰撞电离模型必不可少。Atlas 提供了多种碰撞电离模型我常用的是 Selberherr 模型models materialSilicon IMPACT SELB碰撞电离需要定义电离系数ALPHA和BETA以及临界电场Ecrit。这些参数对击穿电压的计算影响很大。我曾经在一个 LDMOS 的仿真中使用了默认的硅参数结果击穿电压比实验值高了 15%。后来查阅文献发现默认参数是针对均匀电场条件的而 LDMOS 的漂移区电场分布并不均匀需要根据实际电场分布调整参数。3.4 隧道模型当器件中存在薄氧化层或高电场区域时隧穿效应不可忽略。Atlas 提供了直接隧穿和 Fowler-Nordheim 隧穿模型models materialSilicon FNTUNNELFowler-Nordheim 隧穿需要定义势垒高度和有效质量。对于硅/二氧化硅界面势垒高度约为 3.15 eV电子有效质量约为 0.5m₀。这些参数在做 Flash 存储器或薄栅氧 MOSFET 仿真时非常关键。3.5 模型组合的实际案例以一个 0.18 微米 NMOS 的仿真为例我使用的模型组合是models materialSilicon CONMOB FLDMOB SRH AUGER \ IMPACT SELB FNTUNNEL这个组合覆盖了掺杂相关迁移率、场相关迁移率、SRH 复合、俄歇复合、碰撞电离和 Fowler-Nordheim 隧穿。对于这个尺寸的器件漂移-扩散模型已经足够不需要能量平衡模型。实际跑下来I-V 曲线与实验数据的偏差在 10% 以内对于工程仿真来说是可以接受的。注意事项模型不是越多越好。每增加一个模型计算量增加收敛难度也增加。我的原则是先跑通基础模型再根据与实验数据的偏差逐步添加模型。如果添加某个模型后结果没有明显改善就把它去掉。4. 参数调试与结果验证从报错到收敛4.1 常见报错与排查思路Atlas 的报错信息有时候比较隐晦但大多数与材料参数和模型相关的问题都有规律可循。下面是我整理的一份常见问题速查表报错信息可能原因解决方法Material not found材料名称拼写错误或未定义检查material语句中的材料名称Parameter out of range参数值超出物理合理范围检查参数单位确认数值量级Convergence failed模型组合过于复杂或参数不合理减少模型数量逐步添加Negative carrier concentration迁移率或复合参数设置不当检查MUN、MUP、TAUN、TAUPSingular matrix网格或材料参数导致方程病态检查网格质量调整材料参数我曾经遇到过一个典型的收敛问题一个 GaN HEMT 的仿真在栅压接近阈值时无法收敛。排查后发现是CONMOB模型中的MUN参数设置过大导致载流子浓度出现负值。将MUN从 2000 降到 1500 后收敛问题解决。这个案例说明材料参数不仅要物理合理还要与模型和网格匹配。4.2 参数敏感性分析材料参数中有些对结果影响很大有些则不太敏感。根据我的经验敏感性从高到低大致是迁移率参数 复合寿命 禁带宽度 有效状态密度 介电常数。做参数调试时应该优先关注迁移率和复合寿命。我通常用简单的扫描来评估参数敏感性。比如固定其他参数只改变MUN观察漏电流的变化solve init log outfmun_sweep.log solve vdrain0.1 solve vdrain0.5 solve vdrain1.0然后修改MUN值重复上述过程。通过对比不同MUN下的 I-V 曲线可以快速判断这个参数对结果的影响程度。如果影响很小就不需要花太多时间精确标定。4.3 与实验数据的对标材料参数和模型的最终检验标准是与实验数据对标。我通常分三步走第一步对标平衡状态下的电容-电压特性验证材料的基本参数第二步对标线性区的 I-V 特性验证迁移率模型第三步对标饱和区和击穿特性验证高场模型和碰撞电离模型。对标过程中不要试图一次性调整所有参数。我的做法是先固定模型只调参数参数调好后再考虑是否增加模型。每次只调一个参数观察结果变化记录在案。这样虽然慢但能保证每一步都有据可循。实操心得Atlas 的tonyplot工具可以方便地查看仿真结果。我习惯把不同参数下的 I-V 曲线画在同一张图上直观对比。另外log文件中的收敛信息也很有价值如果某个偏置点迭代次数异常多说明该处的模型或参数可能有问题。4.4 收敛性优化技巧收敛问题是 Atlas 用户最常遇到的难题。除了减少模型数量和调整参数还有一些技巧可以提高收敛性。第一使用method语句调整数值方法method newton trapnewton是默认方法trap是梯形积分法对某些问题更稳定。第二使用solve语句的ramp功能逐步增加偏置solve vdrain0.0 vstep0.1 vfinal1.0 namedrain第三调整网格密度。网格太疏会导致精度不足太密会导致计算量过大和收敛困难。我通常在高场区域和结附近加密网格在其他区域适当放疏。5. 进阶话题从材料参数到器件物理的深度理解5.1 材料参数的温度依赖性温度对材料参数的影响是多方面的。除了前面提到的禁带宽度和有效状态密度迁移率也随温度变化。Atlas 提供了TMPMOB模型models materialSilicon TMPMOB这个模型使用幂律关系描述迁移率随温度的变化mu(T) mu300 * (T/300)^(-gamma)对于硅电子的gamma约为 2.4空穴约为 2.2。这意味着温度从 300K 升到 400K 时电子迁移率会下降约 30%。做功率器件仿真时这个效应必须考虑否则结温升高后的电流会被高估。5.2 异质结与能带偏移做异质结器件如 HEMT、HBT时能带偏移是关键参数。Atlas 通过material语句中的AFFINITY和EG300来计算异质结的能带排列。对于 AlGaAs/GaAs 异质结导带偏移约为 0.65 倍的禁带宽度差。我通常这样定义material materialAlGaAs EG3001.79 AFFINITY3.5 material materialGaAs EG3001.42 AFFINITY4.07导带偏移就是两个材料电子亲和能的差4.07 - 3.5 0.57 eV。这个值决定了异质结处的势垒高度对二维电子气的形成至关重要。5.3 缺陷与陷阱模型实际材料中存在的缺陷和陷阱会显著影响器件性能。Atlas 提供了trap语句来定义陷阱trap region2 e.level0.5 density1e17 \ signacceptor lifetime1e-7这个语句在区域 2 中定义了一个能量为 0.5 eV、密度为 10¹⁷ cm⁻³ 的受主陷阱。陷阱会影响 SRH 复合速率和载流子浓度。做可靠性仿真如 NBTI、HCI时陷阱模型是核心。5.4 材料参数与模型选择的工程权衡在实际工程中材料参数和模型的选择需要在精度和效率之间权衡。高精度意味着更多的模型、更细的网格、更长的计算时间。我通常根据仿真目的来决定如果只是做趋势分析用基础模型和默认参数即可如果需要与实验精确对标就需要精细标定参数和模型。一个实用的建议是建立自己的材料参数库。把常用材料硅、二氧化硅、氮化硅、铝、多晶硅等的参数整理成独立的文件每次仿真时直接调用。这样既能保证一致性又能节省时间。我自己的参数库已经积累了十几种材料覆盖了大多数常见器件类型。最后分享一个小技巧Atlas 的示例文件是学习材料参数和模型定义的最好资源。安装目录下的examples文件夹中有大量经过验证的示例涵盖了从简单 PN 结到复杂 HEMT 的各种器件。我建议新手先通读这些示例理解每个参数和模型的设置逻辑然后再自己动手修改。这比从头摸索要快得多。