说到Silvaco TCAD做功率半导体器件的工程师和微电子方向的学生应该都不陌生。它几乎是器件仿真的入门标配工具而肖特基二极管反向击穿特性这个案例又是Silvaco TCAD仿真例库里非常经典的一类结构简单、物理机制清晰、仿真时间短但涉及的知识点却很全——从金属-半导体接触原理、碰撞电离模型到网格设计、数值收敛、参数扫描一套走下来你对TCAD仿真和功率器件反向特性的理解都会上一个台阶。这篇内容我想围绕“肖特基二极管反向击穿特性”这个具体案例把整套仿真思路、关键设置、常见坑位都掰开揉碎讲清楚。适合刚开始接触Silvaco TCAD的朋友也适合已经能跑通简单例程、但一遇到反向击穿仿真就发散或结果不合理的人。我会尽量用实际操练过的口吻来讲哪些地方容易翻车哪些参数值得较真希望你看完能直接照着搭出自己的仿真流程。1. 为什么拿肖特基二极管反向击穿做TCAD仿真的入门案例1.1 肖特基结与PN结反向特性的根本差异肖特基二极管的核心是金属-半导体接触。金属和n型半导体接触时如果金属功函数大于半导体电子亲和势半导体一侧就会形成耗尽区和势垒这就是肖特基势垒。和PN结不同肖特基二极管是多数载流子器件正向导通时电子从半导体注入金属没有少子存储效应所以反向恢复速度极快正向压降也低非常适合做高频整流和电源输出级的续流管。但反向特性上肖特基结恰恰比PN结“软”得多。PN结反偏时漏电主要由耗尽区内热产生载流子的复合决定电流很小而肖特基结反向漏电主要来自热电子发射电子越过势垒从金属流向半导体这个电流和势垒高度直接相关。更麻烦的是反向偏压增大时镜像力会降低有效势垒高度电场越强降低越多漏电流随电压呈指数上升。这就是为什么肖特基管的反向特性曲线看起来有一种“软击穿”的味道而不是像PN结那样在击穿前漏电平台基本不动。1.2 反向击穿的两条物理路径雪崩与隧穿怎么“抢戏”反向击穿机制一般分两类。一类是雪崩击穿耗尽区电场升高到一定程度后漂移过来的载流子在强电场中被加速碰撞电离产生新的电子-空穴对新载流子再被加速、再碰撞形成倍增效应电流迅速增大。另一类是隧穿击穿电场足够强时价带电子直接量子隧穿到导带形成大电流。对硅材料来说掺杂浓度低的时候击穿电压主要由雪崩机制决定掺杂浓度高到一定程度隧穿效应会提前接管。肖特基二极管大多做在低掺杂的n型漂移区上所以设计上一般按雪崩击穿来估算。但这不意味着仿真时可以不开隧穿模型。实际器件的表面、边缘以及高场局部区域电场可能比平行平面近似高不少隧穿很可能在局部提前出现。仿真时把两条路径都打开才能看到更接近真实情况的软击穿特征也才能解释为什么实测击穿电压往往比理想平行平面公式算出来的低。1.3 仿真能帮我们看清什么从伏安曲线到电场分布反向击穿特性的TCAD仿真输出物至少有两类。第一类是反向I-V曲线也就是反向电压从0扫描到击穿点观察电流密度随电压的变化这决定了器件的漏电水平和击穿电压。第二类是器件内部物理量的空间分布比如电场强度、碰撞电离率、载流子浓度的二维或一维分布图这可以帮助我们定位击穿发生在哪里、电场峰值在哪、哪种击穿机制先被触发。只看I-V曲线容易陷入“击穿电压多少V”这个单点结论但TCAD仿真真正的价值在于把曲线背后的机理可视化。比如电场峰值如果落在表面而不是体内说明表面终端处理是关键比如碰撞电离积分达到1的电压和I-V曲线上电流突增的电压是否吻合可以直接判断击穿是不是雪崩主导。这些分析在测试台上很难直接测到但TCAD可以这也是这类案例值得反复做、反复看的原因。2. 动手前理清思路器件结构、材料参数与物理模型选型2.1 结构建在Athena还是Atlas这个案例怎么选Silvaco TCAD两个主力模块Athena负责工艺仿真Atlas负责器件电学特性仿真。常见的流程是先用Athena把离子注入、扩散、氧化、刻蚀这些工艺步骤模拟一遍生成掺杂分布和结构文件再交给Atlas算电学特性。这种流程适合研究工艺变化对器件性能的影响比如注入剂量偏差对击穿电压的影响。但如果你的关注点就是“反向击穿特性”这个电学问题直接拿Atlas构建理想化结构反而更高效。用region语句定义材料区域用doping语句给一个均匀掺杂的漂移区用contact语句设定肖特基接触的功函数整个过程几行代码搞定不需要工艺仿真也不需要等Athena跑完。这样做的好处是参数完全可控比如漂移区掺杂浓度设成1e15 cm^-3就是1e15不会有工艺仿真里杂质再分布带来的不确定因素。对于学习案例来说从Atlas直接建结构把注意力放在物理模型和电学求解上是最合适的路径。2.2 关键材料参数与势垒高度计算别在小事上翻车仿真之前必须把材料参数想清楚。硅的禁带宽度在300K下约1.12 eV电子亲和势4.05 eV相对介电常数11.9。肖特基势垒高度由金属功函数和半导体电子亲和势决定φB φM - χ。举例来说铂的功函数约5.65 eVφB约1.6 eV镍的功函数约5.15 eVφB约1.1 eV。势垒越高反向漏电越小但正向压降也会变大这是一个经典的折中关系。掺杂浓度也会通过镜像力降低效应影响有效势垒。反向偏压下金属-半导体界面处电场很强势垒整体被拉低降低量大概正比于电场的平方根。很多初学者把φB当成固定值算出来的反向漏电比实际小好几个数量级就是这个原因。所以Atlas里定义肖特基接触时除了设置功函数还要确保模型里开启了barrier lowering的修正项否则反向特性在低电压段就会和理论对不上。2.3 物理模型开关怎么打碰撞电离、迁移率、复合物理模型是整个仿真的灵魂。只算常温小偏压特性默认模型可能够用但反向扫描到大电场必须主动把高场相关模型打开。碰撞电离模型是击穿仿真的核心。Silvaco里最常用的是Selberherr模型它把碰撞电离系数描述成电场强度E的关系式里面有a、b、β三个参数。硅材料下电子和空穴的参数不同不同文献略有差异但基本量级一致。仿真击穿时碰撞电离率在耗尽区里积分积分值达到1时对应雪崩击穿。这个模型必须开否则电压加再高电流也不会陡增。迁移率模型在高场下也很关键。浓度依赖迁移率模型conmob描述掺杂浓度对迁移率的影响电场依赖迁移率模型fldmob描述载流子速度饱和效应。反向大偏压下耗尽区电场很高载流子漂移速度趋于饱和这一项直接决定碰撞电离过程中载流子的运动和倍增行为不能省。复合模型方面SRH复合和俄歇复合对反向漏电有影响。低掺杂漂移区里耗尽区热产生电流通常不占主导但为了让曲线更真实一般还是打开。禁带变窄效应bgn在重掺杂衬底区有意义如果阴极附近的n衬底浓度很高建议打开。温度模型如果做常温仿真可以不额外处理但后面要扫温度就必须把温度相关参数一并考虑。2.4 网格布局是最便宜的精度投资TCAD仿真跑出来结果离谱十个里有八个是网格问题。网格太粗电场峰值会被平滑掉击穿电压偏高网格太密计算量翻倍甚至不收敛。肖特基二极管反向击穿关键是在金属-半导体接触界面附近那里势垒最窄、电场最强y方向上必须加密。我的经验是从界面往里0.1 μm范围内网格间距控制在0.001到0.01 μm量级往外到漂移区中部可以逐渐放宽靠近欧姆接触的底部衬底区网格不需要太密。x方向反而简单因为是矩形结构均匀分布就行。但要注意网格长宽比Silvaco对长宽比过大的网格很敏感容易导致矩阵求解病态。加密时最好保持纵横向网格长度接近别在一维方向上疯狂加密。网格这个东西前期多花几分钟想清楚后面收敛性会省很多事。3. 案例实操从零搭出一个可收敛的反向击穿仿真3.1 结构定义与网格生成下面我直接给一套可以跑通的结构定义硅材料宽度2 μm厚度10 μm顶部是肖特基阳极底部是n衬底加阴极欧姆接触。网格的关键点已经写进注释里了。# 定义网格 # x方向均匀分布结构宽度2um mesh x.mesh x.mesh location0.0 spacing0.5 x.mesh location2.0 spacing0.5 # y方向顶部金属-半导体界面处加密底部衬底稀疏 y.mesh location0.0 spacing0.005 y.mesh location0.05 spacing0.005 y.mesh location0.2 spacing0.01 y.mesh location0.5 spacing0.02 y.mesh location1.0 spacing0.05 y.mesh location3.0 spacing0.2 y.mesh location5.0 spacing0.5 y.mesh location10.0 spacing0.5 # 区域定义1为漂移区2为重掺杂衬底 region num1 y.min0.0 y.max5.0 silicon region num2 y.min5.0 y.max10.0 silicon # 电极定义 electrode nameanode x.min0.0 x.max2.0 y.min0.0 y.max0.0 electrode namecathode x.min0.0 x.max2.0 y.min10.0 y.max10.0这段代码里y方向上0到0.05 μm的网格间距设成0.005 μm就是为了捕捉肖特基结附近陡峭的电场变化。如果你的仿真发散或者算出来的击穿电压离谱偏高回头先把这段网格加密再试。底部衬底区域我单独给了region num2后面掺杂时可以区分不同区域。3.2 掺杂、接触与电极设置结构建好之后定义漂移区和衬底的掺杂。漂移区用均匀n型掺杂浓度从1e14到1e16 cm^-3都可以试这里先用1e15 cm^-3。衬底用高浓度n型掺杂1e19 cm^-3模拟重掺杂衬底的低阻特性。# 掺杂定义 doping uniform conc1e15 n.type region1 doping uniform conc1e19 n.type region2接触设置方面阳极是肖特基接触通过指定功函数来实现。我用镍的功函数5.15 eV做示例这样势垒约1.1 eV反向漏电不会太小也不会太夸张。阴极默认是欧姆接触不需要额外指定功函数Atlas会自动在metal-semiconductor界面按欧姆接触处理。# 阳极肖特基接触指定金属功函数 contact nameanode workfunction5.15需要提醒的是contact语句里workfunction的单位是eV别手滑填成伏特。另外如果你用的是铂之类的金属功函数填5.65势垒高度更高反向漏电会小好几个数量级击穿前曲线的“软”程度也会明显不同。这个参数值得自己做几组对比感受一下势垒高度对反向特性的影响。3.3 物理模型与数值求解设置模型配置直接决定反向击穿能不能算出来。针对这个案例我建议开启以下核心模型# 物理模型 models srh auger bgn conmob fldmob impact selb temp300 # 数值方法 method newton trapsrh和auger是复合模型bgn是禁带变窄conmob处理浓度依赖迁移率fldmob处理电场依赖迁移率impact selb就是打开Selberherr碰撞电离模型。这里的temp300表示晶格温度300 K。如果你要扫温度这个值要改成你想要的温度点。数值方法上我习惯用newton trap。newton是全耦合牛顿法trap是遇到不收敛时自动切换到步长减半的重启策略对反向大电压扫描这种强非线性问题特别有帮助。如果你发现发散频繁还可以在method里加itl50之类的参数提高迭代次数上限但别把迭代次数无脑调大先找找网格和模型问题才是正道。3.4 扫描策略与击穿电压提取反向击穿扫描必须从零偏压开始逐步往下加负电压千万别一上来就给定一个大反向偏压硬解。建议先做一个init求解零偏初始状态# 初始求解 solve init # 保存结构 save outfschottky_0.str # 开启日志文件 log outfschottky_reverse.log然后从0开始按0.1 V的步长向负方向扫描。第一遍可以扫到-150 V左右看击穿发生在哪里再根据实际情况调整扫描范围# 反向扫描 solve vanode0 vstep-0.1 vfinal-150 nameanode这里vanode是阳极电压vstep是步长vfinal是终止电压nameanode指定哪个电极在扫。注意如果设在-150 V已经击穿日志里电流会直接冲上去tonyplot打开日志文件就能看到完整的反向I-V曲线。击穿电压的判据怎么定工程上常用反向电流密度达到某个阈值来定义。二维仿真里电流单位是A/μm如果你定义的结构宽度是2 μm电流阈值可以取1e-6 A/μm量级平滑处理一下就能得到比较稳定的击穿电压。也可以在tonyplot里看I-V曲线的拐点电流斜率突变的位置就是击穿点。更定量一点可以用extract语句# 提取击穿电压电流达到1e-6 A/um对应的电压 extract namebv x.val from curve(v.anode, i.anode) where y.val1e-6extract语法不同版本略有差异如果报错直接在tonyplot里用鼠标取点也行。我比较推荐养成用log记录曲线、再统一提取参数的习惯这样后面做参数扫描时可以用脚本一次性提取所有结构的击穿电压效率高很多。4. 结果分析与参数扫描让仿真参与设计决策4.1 tonyplot里怎么看电场与碰撞电离仿真跑完用tonyplot打开最终的结构文件和日志文件能看到两类信息。I-V曲线直接反映了击穿行为但要知道击穿发生在哪、是什么机制必须结合电场分布和碰撞电离率分布一起看。在击穿临界点附近把电场分布图打出来你会发现电场峰值基本贴着金属-半导体的界面沿y方向向体内衰减。这个峰值电场如果超过了硅的临界电场量级比如在5e5 V/cm上下那基本可以判断击穿已经接近临界。再看碰撞电离率分布碰撞电离率的数值越大说明该位置载流子倍增越剧烈。通常它和电场峰值位置重叠。你还可以在tonyplot里沿着结构中心切一条线把电场强度和碰撞电离率画成一条一维曲线观察它们沿深度方向的变化趋势。这是诊断器件设计的利器如果电场峰值贴在表面说明表面状态和终端对击穿的影响巨大实际器件大概率会在表面提前击穿如果电场峰值在体内且宽度较宽平行平面近似可能更接近实际。4.2 漂移区掺杂浓度对击穿电压的灵敏度击穿电压和漂移区掺杂浓度是强相关的。平行平面突变结的击穿电压近似满足VB ≈ εs * Ec^2 / (2 * q * ND)掺杂浓度越高临界电场略高但VB整体下降。用TCAD把不同掺杂浓度的反向特性都跑一遍可以得到一条设计曲线这个比单点计算结果可靠得多。我建议直接把掺杂浓度设成一组扫描值1e14、3e14、1e15、3e15、1e16 cm^-3分别跑反向扫描再把击穿电压提取出来汇总。结果会呈现一个明显规律浓度从1e14提高到1e16击穿电压可能从几百伏快速掉到几十伏以下。这个规律直接决定耐压和导通电阻之间的折中是设计肖特基二极管最核心的权衡。漂移区掺杂浓度 ND (cm^-3)估算击穿电压范围仿真趋势说明1e14数百V量级耗尽区宽电场平缓击穿电压高1e15100 V左右常用中压功率肖特基管范围1e1630 V以下击穿电压低适合低压大电流场景实际仿出来的数值会受网格和模型参数影响但相对趋势和理论公式吻合。做这种扫描时最好的做法是用脚本自动改掺杂浓度、自动跑仿真、自动提取击穿电压别一个一个手动改浪费的时间足够写个循环脚本了。4.3 温度系数与反向漏电容易被忽略的“第二条曲线”功率器件设计不能只看常温特性。温度升高时晶格振动加剧载流子碰撞电离的阈值条件变难雪崩击穿电压通常会略微上升这是击穿电压的正温度系数但反向漏电流就完全相反热电子发射电流随温度指数上升漏电会显著增大。用TCAD做温度扫描也很简单把models里的temp改成需要的温度比如350 K、400 K重新跑反向扫描。你会发现一个很有意思的现象低温时漏电很小击穿比较“硬”高温时漏电基数被抬高I-V曲线在击穿前就明显翘起来击穿点反而不太好定义。这也是为什么功率肖特基管在高温应用里要特别关注反向漏电的原因仿真能提前把这个风险暴露出来。在工程实践中高温反向漏电往往是器件失效的伏笔。TCAD在这里的价值不是告诉你“漏电变大了”这种废话而是帮你定量评估漏电增加多少、击穿电压变化多少以及势垒高度和掺杂浓度如何影响这种温度敏感性。5. 常见问题与排查技巧实录5.1 仿真不收敛先检查网格而不是调参数我见过太多人一遇到“convergence failure”就去翻物理模型、把迭代次数调大结果越调越乱。扪心自问我自己早期也这么干过后来发现至少一半的发散问题根源在网格。肖特基二极管反向击穿仿真最典型的发散场景发生在电压扫描接近击穿点时。击穿前电流急剧变化如果网格在耗尽区边界附近太粗数值上无法准确表达这种突变牛顿迭代就会崩掉。另一种常见情况是网格长宽比太大矩阵条件数恶化求解器直接报错。处理办法很直接把关键区域网格加密控制长宽比把vstep步长减小比如从0.1 V改成0.01 V如果还不行把method设置改成newton trap让求解器在发散时自动减半步长重新尝试。如果一定要给个排查顺序我的习惯是先细化界面附近网格再把扫描步长减半最后才考虑改模型和迭代参数。这个顺序能解决绝大多数问题而且不会引入新的物理偏差。5.2 击穿电压偏离理论值太多按这个顺序排查仿真跑通了但击穿电压和理论公式差很多先从这几个方向查。第一是碰撞电离模型参数。Selberherr模型在Silvaco里有内置参数但不同工艺条件下你可能会看到不同版本的推荐值。如果你自定义了参数务必确认a、b、β数值的单位和来源是合理的。第二是网格粗细。网格太粗会把电场峰值磨平导致击穿电压偏高网格过分细密且某处畸变严重也可能导致局部电场异常击穿电压偏低。把界面附近网格加密再看击穿电压变化如果变化超过10%说明网格还没收敛。第三是结构是否穿通。如果漂移区厚度不够反向偏压还没到雪崩击穿耗尽区就已经穿透到衬底击穿特性会变成“穿通型”电压-电流曲线呈现很软的特征。这时候不是模型问题是结构设计问题增加漂移区厚度就能解决。第四是接触定义。阳极如果用成了欧姆接触而不是肖特基接触反向特性会完全错掉。检查contact语句里的workfunction有没有生效最简单的方法是用tonyplot看零偏时的能带图如果界面处有势垒说明肖特基接触定义正确。5.3 反向漏电流异常物理泄漏还是数值泄漏反向漏电流是我们这个案例里最敏感的物理量之一因为它对势垒高度、温度、复合模型都很敏感。但有时候漏电流大不是物理问题而是数值泄漏。一个容易忽略的点是载流子浓度下限。反向偏压下耗尽区载流子浓度极低数值求解时如果下限设置过高载流子无法被精确求解会出现虚假电流。另一个问题是网格边界效应器件左右边界如果处理成普通的欧姆边界可能会有载流子从边界“漏”进来导致反向漏电虚高。这个案例里左右边界没有电极默认是反射边界一般没问题但如果你自己改了边界条件要小心。区分物理泄漏和数值泄漏可以做一组快速对比把势垒高度提高0.2 eV看看漏电是不是明显下降。如果漏电几乎不变大概率是数值问题如果漏电跟着势垒高度指数变化说明是物理热电子发射主导仿真状态是正常的。5.4 问题速查表最后整理一个速查表方便你实操时对照排查。现象可能原因处理方式扫描到高压段不收敛网格太粗、步长过大加密界面附近网格vstep减半开启trap策略击穿电压严重偏高网格平滑了电场峰值细化y方向网格尤其是界面区域击穿电压严重偏低结构穿通、接触定义错、模型参数错增加漂移区厚度检查电极和功函数反向漏电异常大势垒太低、模型缺失、数值泄漏检查功函数确认barrier lowering检查载流子下限漏电对势垒高度不敏感数值泄漏或接触定义错误检查边界条件确认阳极是肖特基接触温度扫描结果怪异温度参数未正确传入模型检查models中的temp值确认材料参数随温度变化开启这套速查表我贴在手边用了挺长时间。每次仿真出问题照着表走一遍大多数情况能在十几分钟内定位到根因比自己瞎试省时间得多。最后说点我的个人体会这个案例看起来小但对理解TCAD仿真和功率器件设计有很强的牵引作用。我做了很多次反向击穿仿真之后印象最深的不是某条I-V曲线而是一个很朴素的道理TCAD仿真的核心价值在于趋势和机理而不是单点数值。击穿电压算出来是101 V还是105 V意义远不如搞明白“掺杂浓度上升后击穿电压按什么规律下降”“电场峰值为什么会贴在表面”这些问题来得重要。如果刚开始练习建议先跑通一套默认结构、默认模型的流程再逐步改掺杂浓度、改金属功函数、改温度和网格每次只动一个变量观察结果怎么变。这个过程会让你对仿真工具和器件物理建立非常扎实的直觉。还有一个实用的小技巧充分利用tonyplot的cutline功能沿着结构深度方向切一条线把电场、碰撞电离率、能带画成二维曲线很多看起来莫名其妙的结果一画线就真相大白了。这个案例后面还能往很多方向扩展比如加场限环研究终端耐压设计、做JBS结构对比正向反向折中、把工艺仿真导入研究实际掺杂分布的影响。先把反向击穿这个基础动作练扎实后面这些东西都会顺很多。