1. 从一道反向击穿曲线说起为什么肖特基二极管值得单独建一个仿真案例做功率器件或者电源设计的朋友对肖特基二极管应该都不陌生。它的正向压降低、开关速度快在开关电源、续流回路、保护电路里几乎是标配。但很多人对它的认识停留在正向导通0.3V左右这个层面一到反向特性就含糊了——尤其是反向击穿这一段数据手册上往往只给一个击穿电压的典型值曲线画得也很粗略真正设计时心里没底。我在实际项目里就吃过这个亏。早些年做一款低压大电流的同步整流方案选了一颗标称60V的肖特基理论上留了将近一倍的余量结果在高温满载的工况下反复出现漏电偏大、效率掉点的问题最后定位到是反向漏电流随温度急剧上升工作点已经逼近了击穿前的软退化区。那次之后我才意识到光看手册上的击穿电压数字远远不够必须把反向击穿这一段特性真正算清楚、看明白。这就是这个仿真案例的价值所在。用Silvaco TCAD把肖特基二极管的反向击穿特性完整跑一遍你能看到电场在结区的分布、耗尽层怎么随反偏电压扩展、击穿到底是雪崩机制还是隧穿机制主导、温度升高时击穿电压往哪个方向漂。这些信息在数据手册里是拿不到的但对选型和可靠性评估又极其关键。这个案例适合几类人一是做功率器件设计和仿真的工程师想系统掌握击穿仿真的建模方法二是电源、电机驱动方向的硬件工程师需要评估二极管在极端工况下的安全边界三是高校里做半导体器件方向的学生想找一个能跑通、能出图的完整案例练手。下面我就把这个案例从建模思路到参数设置再到结果解读完整拆一遍。2. 肖特基二极管反向击穿的物理图像先想明白再动手建模2.1 肖特基结和PN结在击穿行为上的本质差异很多人习惯用PN结的那套经验去套肖特基二极管这是第一个容易踩的坑。肖特基结是金属和半导体接触形成的势垒多数载流子导电没有少子存储效应所以正向压降低、反向恢复快。但到了反向击穿这一段两者的差异更值得注意。PN结的反向击穿通常由雪崩倍增主导击穿电压对掺杂浓度和结深非常敏感而且击穿曲线比较硬一旦到了击穿点电流上升很陡。肖特基结则不同它的反向特性里漏电流本身就比PN结大好几个数量级而且击穿机制可能是雪崩也可能是势垒降低导致的隧穿或者热电子发射增强具体哪种主导取决于半导体的掺杂浓度和势垒高度。在Silvaco里建模时这个差异直接决定了你要不要开启某些物理模型。比如做PN结击穿雪崩模型基本是必开的做肖特基除了雪崩还得考虑势垒降低效应和镜像力降低否则算出来的反向漏电流会明显偏小击穿点也会偏。2.2 击穿电压到底由什么决定从器件物理的角度肖特基二极管的反向击穿电压主要受三个因素控制半导体外延层的掺杂浓度、外延层厚度、以及金属-semiconductor界面的势垒高度。掺杂浓度决定了耗尽层在给定反偏电压下能扩展多宽也决定了峰值电场能达到多高。外延层厚度则决定了耗尽层能不能充分扩展——如果外延层太薄器件会在耗尽层还没完全展开时就发生穿通击穿这时候击穿电压由厚度决定而不是掺杂决定。势垒高度影响的是反向漏电流的大小势垒越低漏电越大但击穿电压本身对势垒高度不那么敏感。这里有个经验公式可以帮你快速估算对于突变结近似击穿电压和掺杂浓度的关系大致是$$V_{BR} \propto \frac{1}{N_D}$$也就是说掺杂浓度翻倍击穿电压大致减半。这个关系在初步选参数时非常好用。比如你要设计一个100V的肖特基外延层掺杂大概在1e16 cm^-3量级要设计600V的掺杂就得降到1e15甚至更低。2.3 为什么必须用TCAD而不是解析公式解析公式能给你一个数量级的估算但到了实际设计就不够用了。原因有几个一是实际器件的电场分布不是理想的突变结或者线性缓变结边缘效应、场板结构、结曲率都会让局部电场远高于一维估算值二是温度效应、载流子迁移率随电场的退化、碰撞电离率的各向异性这些解析公式都处理不了三是肖特基结的势垒降低和隧穿效应本身就是量子层面的必须数值求解。TCAD的价值就在于把这些效应都放进一个自洽的数值框架里求解。你在Silvaco里设好结构、掺杂、模型它会把泊松方程、连续性方程、电流密度方程联立起来迭代求解最后给你一个物理上自洽的电场和载流子分布。这也是为什么做击穿仿真TCAD几乎是绕不开的工具。3. Silvaco Atlas仿真环境的搭建与结构定义3.1 网格划分击穿仿真里最容易被低估的一步击穿仿真的收敛性和精度八成取决于网格划得好不好。我见过太多人结构定义写得漂漂亮亮一到击穿附近就发散最后发现是结区网格太粗电场峰值根本没被捕捉到。肖特基二极管的反向击穿电场峰值出现在金属-半导体界面附近而且随着反偏电压升高峰值位置会略微向半导体内部移动。所以网格策略是界面附近要非常密往外延层方向逐渐放粗。一个可用的网格设置思路是这样的# 定义Y方向网格深度方向 # 界面附近0.01um步长逐渐过渡到0.1um mesh space.mult1.0 x.mesh loc0.0 spacing0.5 x.mesh loc5.0 spacing0.5 y.mesh loc0.0 spacing0.005 y.mesh loc0.05 spacing0.005 y.mesh loc0.2 spacing0.02 y.mesh loc1.0 spacing0.1 y.mesh loc5.0 spacing0.5注意这里的spacing单位是微米。界面附近0.005um也就是5nm的步长对于捕捉势垒区的电场变化是必要的。外延层深处可以放到0.1um甚至更粗因为那里的电场梯度小。提示网格不是越密越好。网格太密会导致矩阵规模爆炸迭代时间成倍增加而且浮点误差累积反而可能让收敛变差。我的经验是界面附近5nm、过渡区20nm、衬底区100nm这个梯度比较平衡。3.2 掺杂与材料参数从实际外延片参数反推结构定义里第二个关键是掺杂分布。肖特基二极管通常是外延结构重掺杂衬底上生长一层轻掺杂外延层金属接触做在外延层表面。假设我们要仿真一颗100V级别的硅肖特基典型参数是区域厚度掺杂类型浓度外延层5 umN型1e16 cm^-3衬底2 umN型1e19 cm^-3金属接触--势垒高度0.75 eV在Silvaco里定义掺杂可以这样写doping region1 uniform conc1e16 n.type doping region2 uniform conc1e19 n.type材料参数方面硅的禁带宽度、迁移率模型、碰撞电离系数这些Atlas都有内置但要注意碰撞电离系数对温度敏感做变温仿真时得确认模型是否包含温度依赖。3.3 电极定义与肖特基接触的设置要点电极定义看起来简单但肖特基接触的设置有个细节很容易出错。在Atlas里肖特基接触是通过在电极上指定workfunction或者barrier height来实现的。electrode nameanode top workfunction4.8 electrode namecathode bottom这里workfunction4.8 eV对应的是势垒高度约0.75 eV硅的电子亲和能约4.05 eV4.8-4.050.75。如果你直接用barrier height指定也可以electrode nameanode top barrier0.75注意势垒高度这个参数对反向漏电流影响极大。0.75 eV和0.70 eV算出来的漏电流可能差一个数量级。实际建模时最好用CV或者IV测试反推出来的势垒高度不要随便填一个教科书值。4. 物理模型的选择击穿仿真能不能收敛全看这里4.1 迁移率模型从恒定迁移率到场相关迁移率新手最容易犯的错是用默认的恒定迁移率模型去跑击穿。恒定迁移率在低场下没问题但击穿附近电场能到3e5 V/cm以上载流子速度早就饱和了这时候必须用场相关迁移率模型。Atlas里常用的组合是models conmob fldmob srh augerconmob浓度相关迁移率考虑掺杂对散射的影响fldmob场相关迁移率处理速度饱和srh肖克利-雷德-霍尔复合auger俄歇复合高掺杂区重要如果做的是宽禁带半导体比如SiC还得加上各向异性迁移率模型因为SiC的迁移率在不同晶向上差异明显。4.2 碰撞电离模型雪崩击穿的核心碰撞电离是雪崩击穿的物理基础这个模型必须开。Atlas里用impact关键字models impact selbselb是Selberherr碰撞电离模型它把电离率表示成电场的函数$$\alpha \alpha_\infty \exp\left(-\frac{E_{crit}}{E}\right)$$其中alpha_inf和E_crit是材料相关的常数。对硅来说电子的E_crit约1.1e6 V/cm空穴约1.8e6 V/cm。这里有个实操经验如果仿真出来的击穿电压比预期高很多先检查碰撞电离模型开没开。我见过有人跑了一下午击穿电压算出来300V结果发现impact模型忘了加实际器件只有100V。4.3 肖特基特有的模型势垒降低与隧穿肖特基结的反向电流里除了漂移扩散还有两个机制不能忽略镜像力势垒降低和热电子场发射。镜像力降低是因为电子靠近金属表面时会在金属里感应出正电荷这个感应电荷对电子的吸引会降低有效势垒。Atlas里通过models barrier来开启。隧穿则用models tunn对于重掺杂的肖特基隧穿电流可能成为反向漏电的主要成分这时候tunn模型必须开。4.4 数值方法Newton法、Gummel法和收敛参数模型选对了数值方法不对照样发散。Atlas默认用Newton法解耦合方程收敛快但对初值敏感。击穿仿真因为电流变化范围大从pA到mA跨9个数量级建议用以下策略method newton trap maxtrap10 solve init solve vstep0.1 vfinal1.0 nameanode solve vstep0.5 vfinal50 nameanode solve vstep1.0 vfinal100 nameanode分段加电压先小步后大步。在击穿点附近再减小步长solve vstep0.1 vfinal110 nameanode提示如果还是发散可以试试method gummel或者method block。Gummel法收敛慢但更稳适合初值很差的情况。另外climit参数可以限制每次迭代的电流变化量对击穿附近的收敛很有帮助。5. 反向击穿仿真的完整跑通流程与结果解读5.1 从0V到击穿的分段扫描策略把前面所有设置串起来一个完整的反向扫描流程大概是这样# 初始化 solve init # 正向小电压确认接触正常 solve vanode0.1 solve vanode0.2 # 反向扫描 solve vanode-0.5 solve vanode-1.0 solve vanode-5.0 solve vanode-10.0 # 接近击穿区分段 solve vanode-20.0 solve vanode-50.0 solve vanode-80.0 # 击穿附近细扫 solve vanode-90.0 solve vanode-95.0 solve vanode-98.0 solve vanode-100.0 solve vanode-102.0每一步都会输出一个解最后用TonyPlot把IV曲线画出来。击穿点就是电流突然上升的那个拐点。5.2 击穿电压的提取与温度系数的观察跑完一遍之后提取击穿电压的标准做法是在IV曲线上找电流达到某个判据值比如1e-6 A或者1mA对应的电压。不同判据得到的击穿电压会有差异报告数据时要注明判据。温度系数是肖特基击穿仿真的一个重点。硅肖特基的击穿电压通常是正温度系数也就是温度升高击穿电压升高。这跟PN结类似原因是温度升高晶格散射增强载流子平均自由程缩短需要更高的电场才能达到碰撞电离的阈值。做变温仿真时在Atlas里用models temperature300然后分别跑300K、350K、400K对比击穿电压。我实测下来硅肖特基在300K到400K之间击穿电压大概有5%到8%的上升具体数值跟掺杂浓度有关。5.3 电场分布和耗尽层扩展的可视化IV曲线只是结果真正有价值的是内部的物理图像。用TonyPlot看电场分布你能看到低反偏时电场峰值在界面处耗尽层很窄反偏升高耗尽层向半导体内部扩展电场峰值略微下降但积分面积增大接近击穿时电场峰值达到临界值碰撞电离率急剧上升这些图对于理解为什么击穿电压是那个值、为什么提高掺杂会降低击穿电压比任何公式都直观。5.4 仿真结果和实测数据的对标方法仿真跑通了不代表就对了必须和实测对标。对标时重点看三个量反向漏电流在低反偏下的量级、击穿电压的位置、击穿曲线的陡度。如果漏电流仿真值比实测小很多可能是势垒高度设高了或者隧穿模型没开。如果击穿电压仿真值偏高可能是碰撞电离系数不对或者网格太粗没捕捉到峰值电场。如果击穿曲线太陡可能是串联电阻没考虑实际器件的外延层和衬底都有电阻会让击穿后的曲线变缓。对标是个迭代过程通常要调两三轮才能把仿真和实测对上。但一旦对上了这个模型就可以用来做外推和优化价值就出来了。6. 击穿仿真发散时的排查链路我踩过的那些坑6.1 发散的第一反应先看网格再看模型仿真发散是击穿仿真最常见的报错。我的排查顺序是网格→模型→数值方法→电压步长。网格问题最隐蔽。表面上看网格已经够密了但电场峰值可能出现在你没想到的位置。比如肖特基结的电场峰值在界面但如果外延层和衬底界面处掺杂梯度大那里也可能出现电场尖峰。解决办法是先用比较粗的网格跑一个低电压的解看电场分布在哪里集中再针对性地加密。模型问题通常是漏开或者冲突。比如同时开了fldmob和某些高场模型参数不兼容就会发散。这时候逐个关掉模型试能快速定位。6.2 电流跨度太大导致的收敛困难反向漏电流在低反偏时可能只有1e-12 A击穿时到1e-3 A跨了9个数量级。Newton法在这种大跨度下很容易失去收敛。我的应对办法是分两段跑先跑到击穿前比如击穿电压的80%保存解然后从这个解重启用更小的步长跑到击穿。这样初值好收敛概率高很多。# 第一段 solve vanode-80 save outfpre_breakdown.str # 重启 quit然后新开一个deck用mesh infilepre_breakdown.str读入之前的解作为初值继续往下扫。6.3 碰撞电离模型参数不匹配的典型症状碰撞电离模型的默认参数是针对硅的如果你仿的是SiC或者GaN直接用默认参数会得到完全错误的击穿电压。SiC的临界电场比硅高将近一个数量级击穿电压也高得多。症状是击穿电压算出来比预期低很多或者电场还没到临界值就击穿了。这时候要检查材料文件里的碰撞电离系数是不是对应你用的材料。6.4 从发散日志里读出有用信息Atlas发散时会打印一堆信息很多人直接忽略。其实这些日志里有金矿。比如singular matrix通常意味着某个节点没有连接好或者掺杂定义有误voltage step too large就是字面意思减小步长convergence failure at bias xxx告诉你具体在哪个电压点失败可以针对性地在那个点附近加密步长养成看日志的习惯能省很多瞎试的时间。7. 把这个案例用起来从跑通到真正指导设计7.1 用仿真做掺杂-击穿电压的灵敏度分析跑通一个点之后最有价值的用法是做参数扫描。固定结构只变外延层掺杂浓度从5e15到5e16扫一遍看击穿电压怎么变。这个曲线能直接指导你的外延片选型。在Atlas里可以用循环实现loop steps5 assign namend n.value5e15*(1.5^$loop) doping region1 uniform conc$nd n.type solve init solve vanode-100 extract namebv max(v.anode) end这样一次跑完就能得到掺杂和击穿电压的关系。7.2 温度对击穿特性的影响怎么评估前面提到击穿电压是正温度系数但具体到你的器件这个系数是多少得自己算。做300K、350K、400K三个温度点提取击穿电压线性拟合一下就能得到温度系数。这个数据在可靠性评估里很有用。比如你的器件工作温度范围是-40到150度击穿电压的温度系数是正的话低温下击穿电压反而低设计余量要按低温来留。7.3 结构优化场板、结终端对击穿电压的提升平面肖特基二极管的一个问题是结边缘电场集中实际击穿电压往往远低于一维理论值。解决办法是加场板或者结终端结构。在Silvaco里可以很方便地加一个场板在氧化层上再定义一层金属延伸到结外。仿真对比加场板前后的电场分布你能直观看到边缘电场峰值被压下去了击穿电压提上来了。这个优化过程用TCAD做比流片试错成本低太多了。7.4 从单管到模块仿真模型的复用思路单管模型跑通之后可以往上扩展到模块级。比如把肖特基的SPICE模型提取出来放到系统仿真里看开关过程中的反向恢复和电压尖峰。TCAD和SPICE的联合仿真在电源设计里是很实用的技能。提取SPICE模型时重点参数是正向压降、反向漏电、结电容、串联电阻。这些都能从TCAD的IV和CV仿真里提取出来。8. 一些不在手册里的实操心得仿真跑多了有些经验是文档里不会写的。比如网格的spacing值不要用整数用0.005、0.02这种因为Atlas内部对网格位置有舍入整数有时候会导致相邻网格重合。再比如solve的时候电压步长不要用0.1、0.2这种规整值用0.13、0.27这种不规则的值能减少数值共振导致的收敛问题。还有一点TonyPlot看电场分布时默认的色标范围可能把峰值淹没了。手动把色标上限调到峰值电场的1.2倍细节才看得清楚。这些小事看起来不起眼但能让你少走很多弯路。肖特基二极管的反向击穿仿真说到底是一个把物理图像、数值方法和工程经验揉在一起的事情。物理图像帮你判断结果合不合理数值方法帮你把结果算出来工程经验帮你把结果用起来。这个案例跑通一遍你对肖特基器件的理解会上一个台阶后面做SiC、GaN的击穿仿真方法论是相通的。