Silvaco Atlas进阶指南（四）——材料参数与物理模型的实战配置

1. 从“默认”到“定制”：为什么材料参数与物理模型是仿真的灵魂

如果你用过Silvaco Atlas做过一些基础的器件仿真，比如一个简单的硅基MOSFET，你可能会发现，有时候你只需要定义好几何结构、掺杂和网格，然后跑一个简单的MODELS MOS，仿真就能跑起来，并且给出一个看起来还不错的结果。这很容易给人一种错觉：材料参数和物理模型嘛，用默认的就好了，反正软件都内置好了。

我刚开始用Atlas的时候也是这么想的，直到我尝试仿真一个稍微复杂点的器件——比如一个GaN HEMT（高电子迁移率晶体管）。那次仿真简直是一场灾难：仿真要么不收敛，直接报错退出；要么勉强收敛了，但输出的电流-电压曲线和实验数据差了十万八千里，跨导的形状完全不对。我折腾了好几天，反复检查网格和掺杂，都无济于事。最后，还是一位有经验的同事点醒了我：“你用的迁移率模型还是硅的那套默认参数吧？GaN的电子饱和速度和临界电场跟硅能一样吗？”

这句话让我醍醐灌顶。原来，Atlas提供的那些“默认”参数，很大程度上是基于硅（Silicon）这种成熟材料优化好的。当你换一种材料，比如砷化镓（GaAs）、氮化镓（GaN），或者设计一个由多种材料组成的异质结器件时，这些默认值就不再适用了。它们就像一套为标准身材准备的“均码”衣服，穿在硅身上可能刚好，但穿在GaN或者InP身上，要么太紧（参数过高导致仿真失真），要么太松（参数过低无法体现材料特性）。

所以，材料参数与物理模型的配置，本质上是一个“量体裁衣”的过程。它的目标是将仿真软件的物理引擎，精准地“校准”到你所要研究的真实材料和器件物理上。这不是一个可选的“高级功能”，而是从“能仿真”迈向“仿真得准”的必经之路。一个配置得当的模型，能让仿真结果无限贴近现实；而一个粗糙的默认配置，则可能让你在错误的方向上越走越远，得出完全误导性的结论。接下来，我们就抛开那些笼统的概念，直接进入实战，看看怎么给我们的器件“量体”并“裁衣”。

2. 接触（CONTACT）定义：别小看电极，它是器件的“门户”

在仿真中，我们经常把注意力集中在器件有源区的那些复杂物理过程上，而容易忽略电极接触的定义。但实际上，电极是器件与外部电路交互的唯一“门户”。门户定义得不清楚，里面的情况再真实，传递出来的信号也是失真的。CONTACT语句就是用来精确定义这个门户属性的。

2.1 欧姆接触与肖特基接触：一字之差，天壤之别

这是最容易出错的地方之一。Atlas有一个默认规则：如果你不给接触指定功函数（WORKFUNCTION），它就认为这是一个理想的欧姆接触。欧姆接触意味着接触处的电压降为零，电流-电压呈线性关系，这是重掺杂半导体与金属接触的常见情况。比如，对于MOSFET的源漏区，我们通常就用欧姆接触。

但是，当你需要模拟一个肖特基二极管、一个HEMT的栅极，或者任何存在金属-半导体势垒的结构时，就必须明确地告诉Atlas：“这是个肖特基接触”。方法就是使用CONTACT语句的WORKFUNCTION参数。

# 错误示例：将肖特基栅极误设为欧姆接触 CONTACT NAME=gate # 软件会认为gate是欧姆接触，势垒高度为0，仿真结果完全错误 # 正确示例：定义一个功函数为4.8 eV的肖特基栅极 CONTACT NAME=gate WORKFUNCTION=4.8

这里有个非常实用的技巧：对于常用的金属，Atlas内置了材料名，你可以直接使用，而不用去记具体的功函数数值。比如，ALUMINUM代表铝（功函数约4.1 eV），N.POLYSILICON代表N型多晶硅（功函数约4.05 eV）。这样写更直观，也不容易出错：

CONTACT NAME=gate N.POLYSILICON # 使用N型多晶硅的功函数

我踩过的一个坑：曾经仿真一个肖特基势垒二极管，明明定义了功函数，但反向击穿电压总是比预期低很多。后来才发现，是忽略了肖特基势垒降低效应。在高电场下，肖特基势垒的高度会因镜像力降低而下降，导致电流提前增大。在Atlas中，可以用BARRIER和ALPHA参数来模拟这个效应：

CONTACT NAME=anode WORKFUNCTION=4.9 BARRIER ALPHA=1.0e-7

这条语句在定义功函数的同时，启用了势垒降低模型，并设置了降低系数。ALPHA=1.0e-7中的单位是厘米，代表1纳米，这是一个典型的尺度。加上这个模型后，反向特性立刻变得合理了。记住，当涉及肖特基接触时，务必在接触下方的半导体区域设置足够精细的Y方向网格，否则软件无法准确解析那个非常薄的耗尽区，什么高级模型都白搭。

2.2 电流边界与浮动节点：模拟特殊工作状态的利器

电压控制是我们最常用的方式。但有些器件特性，你必须用电流来控制电压才能得到。最典型的例子就是半导体器件的击穿特性，比如二极管的反向击穿、BJT的雪崩击穿。在击穿点附近，电压对电流极其敏感，微小的电压变化会导致电流急剧上升，甚至出现“负阻”现象（电流增大，电压反而下降）。如果用电压源去扫，求解器很难在这个陡峭的曲线上找到稳定解，极易发散。

这时，就需要把电极设置为电流控制：

CONTACT NAME=collector CURRENT

这样，你在SOLVE语句中设定的就是电流值，Atlas会去求解对应的电压是多少。这能极大地提升击穿仿真这类难题的收敛性。不过要注意，使用电流边界条件后，你必须选用BLOCK或NEWTON这类更强大的非线性求解器，简单的GUMMEL迭代是应付不来的。

另一个重要的概念是浮动接触。比如，在仿真一个EEPROM存储单元时，浮栅上没有物理电极连接，它的电位是由周围电容耦合决定的。模拟这种结构，就需要用到FLOATING参数：

CONTACT NAME=floating_gate FLOATING

这告诉Atlas，这个电极是浮空的，其上的净电流为零（电荷守恒边界条件）。对于直接做在半导体上又没有引线的浮动场板，则不能直接用FLOATING。我常用的变通方法是给它串联一个极大的电阻来模拟开路：

CONTACT NAME=field_plate RESISTANCE=1e20

这个电阻值（1e20欧姆）要足够大，使得流过它的电流远小于器件其他支路的电流，从而在电学上近似为断开。这个方法比直接删除电极或设置零电流边界更稳定，更容易收敛。

2.3 寄生元件与电极短接：连接真实世界的桥梁

真实的器件从来不是孤立的，它总是通过导线、键合线连接到测试电路或芯片内部，这些连接会引入寄生电阻、电容和电感。在仿真中忽略它们，尤其是在高频或大电流应用中，会导致严重误差。CONTACT语句可以方便地在电极上并联或串联集总元件。

CONTACT NAME=drain RESISTANCE=5 CAPACITANCE=1e-12 INDUCTANCE=1e-9

这条语句在漏极并联了一个5欧姆的电阻、1皮法的电容和1纳亨的电感。这对于模拟器件的开关瞬态、振荡行为或匹配特性至关重要。这里有一个关键的缩放问题需要注意：在2D仿真中，你定义的电阻值R，实际代表的是“方块电阻”或“每微米宽度的电阻”。如果仿真结构在第三维的宽度是W（默认为1微米），那么实际的电阻是R/W。如果你的器件宽度是10微米，想模拟一个50欧姆的串联电阻，这里应该填RESISTANCE=500（50欧姆 * 10微米）。

电极短接功能在仿真多指条晶体管、并联器件或某些特殊结构（如双基极晶体管）时非常有用。它能让多个物理上分开的电极保持电位相等。

ELECTRODE NAME=base1 ... ELECTRODE NAME=base2 ... ... CONTACT NAME=base1 COMMON=base SOLVE VBASE=0.5

这里，base1电极被短接到名为base的公共电极上。你只需要对base电极施加0.5V偏压，base1就会自动获得相同的电压。Atlas会分别计算两个电极的电流，如果你需要知道总电流，可以用SHORT语句将它们合并输出。这个功能在构建复杂器件版图对应的紧凑模型时，能节省大量时间。

3. 材料（MATERIAL）属性：为每一种物质注入“个性”

如果说几何结构是器件的“骨架”，那么材料属性就是它的“血肉”。MATERIAL语句让你能够精确地定义或覆盖仿真区域内任何一种材料的本征物理参数。

3.1 基础参数定义：从硅到宽禁带半导体

对于同质材料，比如整个衬底都是硅，你可以全局设置其参数：

MATERIAL MATERIAL=Silicon EG300=1.12 AFFINITY=4.05 NC300=2.8e19 NV300=1.04e19

EG300是300K下的带隙，AFFINITY是电子亲和势，NC300和NV300是导带和价带的有效态密度。这些是决定半导体能带结构的最核心参数。迁移率参数同样关键：

MATERIAL MATERIAL=Silicon MUN=1450 MUP=500 VSATN=1.0e7 VSATP=1.0e7

MUN和MUP是低场电子和空穴迁移率，VSATN和VSATP是饱和速度。对于硅，这些默认值可能够用。但当你切换到宽禁带半导体如碳化硅（SiC）或氮化镓（GaN）时，必须手动修改：

MATERIAL MATERIAL=GaN EG300=3.39 AFFINITY=4.1 NC300=2.3e18 NV300=4.6e19 MATERIAL MATERIAL=GaN MUN=1000 VSATN=2.5e7

你看，GaN的带隙（3.39 eV）远大于硅，饱和速度（2.5e7 cm/s）也高得多，直接沿用硅的参数会严重低估GaN器件的电流能力和高频性能。

更常见的情况是，同一个器件内不同区域的材料相同，但我们需要设置不同的局部属性。比如，一个双极晶体管（BJT）的基区载流子寿命和发射区就不同。这时可以用REGION或NAME参数来指定：

# 假设region 2是基区 MATERIAL REGION=2 TAUN0=2e-7 TAUP0=1e-6 # 设置SRH复合寿命 # 或者，如果之前用REGION NAME=base定义了基区 MATERIAL NAME=base TAUN0=2e-7

TAUN0和TAUP0分别是电子和空穴的SRH寿命。在基区设置较短的寿命，可以模拟复合对电流增益的影响。

3.2 异质结材料：能带对齐的奥秘

异质结器件，比如HBT（异质结双极晶体管）或HEMT，其卓越性能源于不同材料界面处的能带突变。在Atlas中模拟异质结，除了定义各自材料的参数，最关键的一步是正确设置能带偏移的分配。

当两种材料形成异质结时，它们的导带底和价带顶会发生相对移动，这个移动量就是能带偏移。总带隙差（ΔEg）有多少分给导带（ΔEc），多少分给价带（ΔEv），是由材料的电子亲和势和界面偶极子等因素决定的。在Atlas中，我们用ALIGN参数来定义这个分配比例。

ALIGN参数表示导带偏移占总带隙差的比例。它的值在0到1之间。

• ALIGN=0.7表示70%的带隙差分配给了导带，30%给了价带。这会形成一个较高的电子势垒和较低的空穴势垒。
• ALIGN=0.3则相反，形成低电子势垒和高空穴势垒。

举个例子，假设我们仿真一个AlGaAs/GaAs异质结，已知AlGaAs的带隙比GaAs大0.3 eV。如果我们从文献或实验知道，这个差值的约65%体现在导带上，那么：

MATERIAL MATERIAL=AlGaAs EG300=1.80 ALIGN=0.65 MATERIAL MATERIAL=GaAs EG300=1.50 ALIGN=0.65

这样，Atlas会自动计算：ΔEc = 0.3 eV * 0.65 = 0.195 eV， ΔEv = 0.3 eV * 0.35 = 0.105 eV。这个能带对齐方式直接决定了HEMT中二维电子气（2DEG）的 confinement 和 HBT中载流子的注入效率，如果ALIGN参数设错，整个器件的输运机制就全乱了。

对于更复杂的多元化合物半导体，如InGaAs、AlInGaN等，其参数（如带隙、晶格常数）会随组分（x, y）连续变化。Atlas支持通过MATERIAL语句定义这些组分相关的参数，通常需要配合查找表或解析公式。这是异质结仿真中最需要查阅具体材料数据库的部分。

4. 界面（INTERFACE）与物理模型（MODELS）：刻画微观物理的“画笔”

定义好材料和接触，相当于准备好了画布和颜料。而界面属性和物理模型，则是我们用来刻画器件内部复杂物理现象的“画笔”。用什么样的画笔，决定了你能画出多精细的图画。

4.1 界面属性：半导体表面的“疤痕”

半导体与绝缘体（如SiO₂）的界面，或者不同半导体之间的异质界面，从来都不是完美的。那里存在悬挂键、杂质等，导致界面处有电荷、能态，严重影响器件性能。INTERFACE语句就是用来描述这些“疤痕”的。

最常见的界面属性是固定界面电荷（Fixed Interface Charge），用QF参数表示。在硅基MOSFET中，这就是影响阈值电压的固定氧化物电荷。

INTERFACE QF=1e10

这条语句会在所有半导体-绝缘体界面施加1e10 cm⁻²的固定正电荷。但通常，我们只关心特定区域的界面，比如MOSFET的沟道区。这时可以用坐标范围来限定：

INTERFACE QF=1e10 X.MIN=0.5 X.MAX=1.5 Y.MIN=0.0 Y.MAX=0.05

这样，电荷只被添加到指定的矩形区域内的界面上，非常灵活。

除了固定电荷，界面复合速度对器件的漏电和少数载流子寿命至关重要。特别是在光伏器件或高灵敏度的光电探测器中，表面复合是主要的损失机制。

INTERFACE SN=1e3 SP=1e3

SN和SP分别代表电子和空穴的表面复合速度，单位是cm/s。数值越大，表示表面复合越强。

4.2 物理模型选择：从漂移扩散到能量平衡

MODELS语句是Atlas仿真的“大脑中枢”，它决定了软件用哪些物理方程来描述载流子的运动。模型的选择不是越高级越好，而是适合你的仿真目标和器件类型。

基础模型集：对于常规的、工作在较低电场下的硅基器件（如CMOS），一套标准的模型集通常就足够了：

MODELS MOS CONMOB SRH FERMI

• MOS: 启用适用于MOSFET的默认模型集（包括迁移率退化等）。
• CONMOB: 载流子浓度相关的迁移率模型（考虑掺杂浓度对迁移率的影响）。
• SRH: 肖克利-里德-霍尔复合模型。
• FERMI: 费米-狄拉克统计（对于高掺杂区域，必须用它来替代默认的玻尔兹曼统计，否则载流子浓度会算错）。

进阶模型：当你仿真功率器件、高频器件或者工作在强场下的情况时，就需要引入更精细的模型。

• 高场饱和与速度过冲：FLDMOB或FLMOB模型描述了迁移率随电场增强而下降直至饱和的效应。对于化合物半导体，其速度-电场曲线与硅不同，必须使用正确的参数：

  MODELS FLDMOB VSATN=2.0e7 ECRITN=4.0e3

  这里ECRITN是电子迁移率开始显著下降的临界电场。
• 碰撞电离与雪崩击穿：这是仿真二极管反向击穿、MOSFET漏致势垒降低（DIBL）等效应的关键。通过IMPACT语句启用：

  IMPACT SELB

  SELB是一个常用的碰撞电离系数模型。启用后，Atlas会计算高电场区域产生的电子-空穴对。
• 能量平衡模型（Energy Balance）：这是最高级的漂移-扩散类模型。在纳米尺度器件或强非平衡状态下（如亚微米MOSFET的沟道热载流子），载流子的温度可能远高于晶格温度，导致速度过冲等效应。标准的漂移扩散模型无法描述这些，必须启用能量平衡模型：

  MODELS HCTE.HO

  HCTE.HO同时启用电子和空穴的能量平衡输运。这个模型计算量巨大，但能更准确地模拟深亚微米器件的输运特性。

4.3 异质结器件的模型分区配置

这是材料参数与模型配置相结合的高级技巧。在一个AlGaN/GaN HEMT中，AlGaN势垒层和GaN沟道层是两种材料。GaN沟道中的电子迁移率很高，且受极化场影响；而AlGaN层中几乎无电流。因此，我们可以对不同的材料区域启用不同的迁移率模型：

# 在GaN沟道区使用高场迁移率模型和特定的饱和速度 MODEL MATERIAL=GaAs FLDMOB VSATN=2.5e7 ECRITN=1.5e5 # 在AlGaN势垒层，可以使用简单的常数迁移率模型，甚至不激活输运（如果它是完全耗尽的） MODEL MATERIAL=AlGaAs CONMOB MUN=300

这种按材料分区配置模型的能力，使得我们可以用最合适的计算资源去描述器件中不同部分的核心物理，在保证精度的同时提升仿真效率。

最后，分享一个我调试模型时的习惯：善用PRINT或OUTPUT语句。在MODELS语句后加上PRINT参数，Atlas会在日志文件中详细列出所有被激活的模型及其使用的参数值。这是验证你的配置是否按预期生效的最直接方法，能帮你快速定位是哪个参数设错了，还是想要的模型根本没被启用。模型配置这条路，坑不少，但每填平一个坑，你对器件物理和仿真工具的理解就会深一层。