简介《各种语句sentaurus.docx》是一份面向微电子、纳米电子与光电子器件仿真学习者的 Sentaurus 语句参考文档适合已具备 TCAD 基础、需要查阅或复用仿真命令的读者。内容围绕 I-V、C-V 曲线计算辐照损伤与重粒子入射以及瞬态仿真展开涵盖电流电压文件、电极与物理模型、Math 数值方法、Solve 求解流程、AC 耦合电容提取、HeavyIon 重离子撞击和激光光照等典型语句片段。读者可据此快速定位关键参数设置、模型选择及输出控制方式理解器件仿真输入文件的组织逻辑并为搭建、调试和优化仿真 deck 提供可直接对照的模板减少从零编写命令时的试错成本。压缩包共 1 个 docx 文件约 16KB体积轻便便于检索与摘录。该文档已有 2038 人学习适合作为 Sentaurus 常用语句的速查资料和课程实验、科研建模的辅助参考。1. 一行 Physics 开关为什么曲线形状和收敛性会绑在一起在 Sentaurus 这类 TCAD 工具里最反直觉的一件事是物理模型和数值收敛并不是两件事。你在 Physics 段里多加一个Avalanche (UniBo Eparallel)曲线上会多出一条击穿分支但同一个文件里 Math 段的Notdamped、BreakCriteria和 Solve 段的Coupled括号顺序决定了这条分支是被算出来还是让整个仿真在偏置扫到一半时崩掉。很多人第一次遇到发散会本能地怀疑模型不对真正的原因却常是求解路径没铺好。这套语句真正值得拆开看的地方是它把“写哪个方程”“按什么顺序解”“结果往哪儿落”三件事拆成了独立层级每层都能单独调试。搞懂分层比背参数有用。2. 声明与装配Device/File/Electrode/Physics 为什么必须分层写拿到一份能跑的 Sentaurus 语句第一眼看到的是几个平级的块File、Electrode、Physics、Math、Solve。看起来像一份配置清单实际是一次自上而下的装配先声明数据从哪来、往哪去再声明边界条件再声明方程里带哪些物理项最后才声明怎么解。顺序不能乱因为后面的块依赖前面的声明。2.1 File 块为什么输入输出要分成四条独立的流一个二极管 I-V 求解的开头通常是这样File { Grid Diode_msh.tdr # 输入网格 材料区 掺杂分布由结构编辑器生成 Current Diode_IV.plt # 输出电极上的标量电压/电流序列一条曲线一个文件 Plot Diode_IV.tdr # 输出全场量快照电势、载流子密度、电场矢量 Output Diode_IV.log # 输出求解日志迭代次数、残差、步长回退记录 }这四条流分开写不是为了好看。Current存的是电极上的一维时间序列几百个点几十 KBPlot存的是全场量每个节点十几个物理量动辄几百 MB。把它们塞进一个文件代价是每次只想看一条曲线都要读全场数据。更重要的是第四条Output日志是唯一能解释“为什么没算完”的地方。求解器不会告诉你物理哪里错了它只会告诉你第几步、残差多少、步长被砍到多少。把它当成黑匣子里的飞行记录仪出问题第一件事就是翻它而不是改模型。工程实践中还有一个细节Plot里的量由Plot段二次筛选File只负责命名目的地。两者分离的设计让你可以重跑后处理而不重算——前提是NoOverwrite这类开关用对否则文件名被覆盖后悔药就没了。2.2 Physics 块每个模型开关都是对方程组的一次修改这是整套语句里信息密度最高的地方Physics { AreaFactor 1 # 二维截面结果折算为等效三维时的宽度因子 Temperature 300 # 晶格温度单位 K不是载流子温度 Mobility ( DopingDependence # 迁移率随掺杂浓度退化 eHighFieldSaturation # 电子高场速度饱和 hHighFieldSaturation # 空穴高场速度饱和 EnormalCarrier # 表面法向电场引起的额外退化 CarrierCarrierScattering # 载流子-载流子散射 ) Recombination ( SRH (DopingDependence) # 深能级复合寿命随掺杂下降 Auger (withGeneration) # 三载流子过程高注入下主导 Avalanche (UniBo Eparallel) # 碰撞电离按平行电场驱动 Band2Band (Hurkx) # 带间隧穿折算成等效产生率 EffectiveIntrinsicDensity (OldSlotboom) # 重掺杂下的禁带变窄 ) }关键在于这些括号不是“精度档位”每一个都往漂移-扩散方程组里加一项加项就改雅可比矩阵的结构和条件数。SRH让复合率依赖寿命 τ(DopingDependence)再把 τ 写成掺杂浓度的函数于是高掺杂区的少子寿命被压到纳秒级反向饱和电流随之变化Auger (withGeneration)在 SRH 之外再加一个 ∝n²p 的复合通道withGeneration让它允许把多余能量转给另一个载流子去激发这在重掺杂发射区里不能省。Avalanche (UniBo Eparallel)更典型。碰撞电离率对电场是指数敏感的Eparallel指定驱动电场只取平行于电流方向的分量——如果取全电场模在结边缘或异质结界面上电场方向不纯会凭空触发大量电离击穿电压直接跑偏。这类模型细节在曲线上表现为“击穿点提前了 20%”不去拆开关根本定位不到。Band2Band (Hurkx)和EffectiveIntrinsicDensity (OldSlotboom)属于另一类前者把带间隧穿折算成局部产生率避免对隧穿路径做能量积分后者用经典模型描述禁带变窄让本征载流子浓度在重掺杂区抬升直接改变 PN 结的注入比。它们的共同点是——不加结果可能看起来“还算像”加上才和实测对得上但加了以后收敛变难。这就是物理精度与数值鲁棒性的第一次正面冲突。2.3 Electrode 与 AreaFactor二维结果怎么变成“等效三维”Electrode { NameElectrodeN Voltage 0.0 }只做两件事给接触起名给初值。名字非常重要——Plot段、BreakCriteria段、混合模式里的电路节点映射全都靠这个名字索引。名字写错语句照样能跑只是断点判据失效电流跑飞也没人拦。AreaFactor 1是另一个容易忽视的换算口。二维仿真的电流天然是“每单位深度”的电流A/μm真实器件的电流要乘以宽度。设成 1 意味着你在看“每微米宽度”的结果。忘了改这一项得到的电流可能和实测差两个数量级而曲线形状完全正确——这是最阴的一类错因为形状对人会先怀疑物理参数。提示二维结构里做校准之前先把 AreaFactor 按实际器件宽度设好再把实验结果换算到同一口径。顺序反了后面所有拟合都是白做。3. 求解序列Coupled 的括号顺序比模型精度更决定成败模型定了接下来是解。漂移-扩散是一组强非线性耦合方程泊松方程管静电势两条连续性方程管电子和空穴任何一条单独解都得不到有意义的电流。真正的难点不是“要不要耦合”而是“从什么状态开始耦合”。3.1 从 Poisson 到 Electron/Hole/Poisson 的渐进耦合Solve { Coupled (Iterations50) { Poisson } # 第一步只要静电势 Coupled (Iterations15) { Hole Poisson } # 加空穴连续性 Coupled (Iterations15) { Electron Hole Poisson } # 再加电子连续性 QuasiStationary ( InitialStep 1e-6 MaxStep 0.01 MinStep 1e-9 Goal { NameElectrodeN Voltage 1000 } Plot { Range (0 1) Intervals 2 } ) { Coupled { Hole Electron Poisson } # 每个偏置点内部仍是全耦合牛顿 Plot (FilePrefixIV_ Time(0.01;0.05;0.1;0.5) NoOverwrite) } }前三条Coupled是渐进式装配先解泊松拿到电势分布把它当作空穴方程的初值再拿两者当电子方程的初值。零初值下直接上三变量耦合载流子密度在数量级上跨越十几个量级牛顿迭代的线性系统条件数极差几乎必然发散。分步走相当于给每一步一个“差不多合理”的起点代价是多花几次迭代。注意每步的Iterations不一样泊松 50 次后面 15 次。泊松方程相对好解但初值最差给宽一点连续性方程已有前一步的场做引导15 次够用。这个数字不是固定的网格加密后经常要调大。QuasiStationary里的Goal是“扫到某个电压”求解器自己决定中间踩哪些点。InitialStep、MaxStep、MinStep控制自适应步长步长太大越过击穿拐点牛顿在无解区反复回退步长太小几百个偏置点白烧机时。MaxStep 0.01意味着单步不超过 0.01 V在 1000 V 的扫描范围里意味着至少 10 万个点——这显然不是常规用法更像是先给一个保守上限确保能跑完。真实项目里通常会先粗后细用大MaxStep跑通全程看形状再在小范围内收紧步长抓拐点。3.2 Math 块收敛失败的病根常在导数Math { # Cylindrical # 柱坐标对称注释掉说明当前是笛卡尔坐标 Method Pardiso # 稀疏直接法求解线性系统 Number_of_threads 4 Stacksize ... # 递归求解时防止栈溢出 ExtrapolateDerivatives # 导数用外推近似省一次函数求值 AvalDerivatives # 碰撞电离项参与雅可比必须给真导数 RelErrControl # 误差控制基于相对量 Iterations 15 # 单个偏置点允许的最大牛顿迭代 Notdamped 60 # 前 60 次迭代不做步长阻尼 BreakCriteria { Current (ContactElectrodeN maxval 1e-8) } }Pardiso是直接法分解矩阵一次解到底鲁棒但内存随节点数非线性增长。强非线性问题击穿、雪崩、辐照瞬态一般不用迭代法因为迭代法省下的内存会被“反复不收敛”加倍还回去。ExtrapolateDerivatives和AvalDerivatives是一对反向操作。牛顿法需要雅可比矩阵里的 ∂F/∂x有些项的解析导数代价很高外推近似能省时间但会破坏二次收敛碰撞电离系数是电场的指数函数导数错一点迭代方向就差很远所以必须用真实导数。这两个开关的取舍原则很朴素线性或弱非线性项可以外推指数型强非线性项不行。Notdamped 60和Iterations 15之间其实有个矛盾单点最多 15 次迭代可阻尼要到 60 次才启用。等于说在这套设置下阻尼形同虚设。这是血泪经验的产物——对某些结构早开阻尼会让迭代在局部极小值附近磨蹭反而不如放开步长让它冲过去。但这个值必须配合BreakCriteria使用否则冲过头就是发散。BreakCriteria是整套语句里最像后悔药的一行接触电流超过 1e-8 就停。它不改变物理只在数值失控时把进程从跑飞状态拉回来并留下一个明确的终止位置方便你回去看那一步的日志。4. 外部激励AC 小信号、重离子与激光为什么共用一套骨架I-V 是直流稳态C-V 是频率域的小扰动重离子和激光是时域的大扰动看起来是三套完全不同的东西。但它们在语句结构上高度同构都要先有一个收敛的稳态解再把外部激励作用上去再沿着某个自变量电压、频率、时间推进。理解这个同构就不会把它们当成三份互不相干的模板去背。4.1 混合模式与 ACCoupledC-V 为什么必须挂在直流解上System { Diode diodesystem (ElectrodeNfront ElectrodeP0) # 器件电极映射到电路节点 Vsource_pset vn (front 0) { dc 0 } # 电压源dc 是被扫描的量 } File { Output Diode_CV_1kHz.log ACExtract Diode_CV_AC_1kHz.plt } Solve { Coupled (Iterations50) { Poisson } Coupled (Iterations15) { Hole Poisson } Coupled (Iterations15) { Electron Hole Poisson } Coupled (Iterations15) { Electron Hole Poisson Contact } QuasiStationary (InitialStep1e-6 MaxStep1e-2 MinStep1e-7 Increment2 Decrement4 Goal { Parameter vn.dc Voltage 100 }) { ACCoupled (StartFrequency1e3 EndFrequency1e3 NumberOfPoints1 Decade Iterations15 Node(front) ACMethodBlocked ACSubMethod(diodesystem)ParDiSo) { Poisson Electron Hole Contact Circuit } } }小信号分析的数学前提是线性化把解写成 x x₀ Re(δx·e^{jωt})代入原方程后丢掉二次以上的小量得到一组复数线性方程。x₀ 就是直流工作点没有它就没有扰动可加。所以ACCoupled必须嵌在QuasiStationary内部每扫到一个偏置点先解直流再在这个点上加正弦扰动最后从复导纳里取出电容。这也是为什么ACExtract输出的是“电容-电压”而不是“电容-频率”——频率在这里只是求导的算子不是扫描变量。StartFrequency EndFrequency 1e3加NumberOfPoints 1等于只算 1 kHz 单点。这是 C-V 的常见用法固定频率扫偏置。若真想做频率扫描这两个值要拉开并给点数此时Decade才起作用。ACMethod Blocked把器件方程和电路方程分块求解ACSubMethod(diodesystem) ParDiSo指定该子系统的线性解法。分块的代价是丢掉器件-电路之间的一部分交叉项在“一个器件加一个理想源”的拓扑下几乎无影响一旦电路侧有多个强耦合元件就要谨慎评估。Node(front)指明扰动从哪个电路节点注入Contact出现在未知量列表里意味着接触被当作额外变量处理对接触电阻和边界通量的处理更稳。低频 1 kHz 这个选择本身也有物理含义频率足够低时深能级陷阱跟得上信号测到的是“含陷阱响应”的电容频率高到陷阱来不及响应电容会掉一截。用单频算 C-V等于默认陷阱响应可以忽略或已经包含在直流模型里。4.2 HeavyIonLET 波形怎么变成体电荷Physics { HeavyIon ( Direction (0, 1) # 入射方向 Location (200, 0) # 入射点 Time 1e-9 # 入射时刻 Length [ ] # 径迹长度区间 Wt_hi [ ] # 径迹横向分布的权重 LET_f [0 8.7e-6 8.7e-6 0] # 线性能量转移随时间的分段线性波形 Gaussian # 横向按高斯分布铺开 PicoCoulomb # 电荷量纲避免与能量量纲混用 ) }机制链条很短但每一步都可能出错LET(t) 给出单位长度上的能量沉积率除以产生一对电子空穴所需的平均能量硅约 3.6 eV得到单位长度、单位时间的载流子对产生率再按Gaussian的横向分布摊到径迹周围的网格节点上转成体产生率塞进连续性方程。LET_f那四个数就是波形的骨架起、升、平台、落。平台段的斜率决定了瞬态电流尖峰有多陡也直接决定后续Transient的MinStep要压到多小。Location和Direction决定径迹穿过哪些区域——打在结区还是衬底产生的电荷收集效率完全不同这是单粒子效应仿真里最敏感的几何参数。Gaussian不能省。径迹在物理上有横向宽度如果强行让所有电荷落在一列节点上结果会强烈依赖网格网格换一套峰值电流就变一次。加上横向分布后电荷铺到若干列节点上网格依赖性显著下降。注意PicoCoulomb是量纲声明不是单位换算开关。LET 数据、器件尺寸、体电荷密度三者口径不一致时得到的瞬态电流可能差好几个数量级而且曲线形状还是对的。4.3 激光与瞬态步长为什么 MinStep 要压到 1e-17Optics ( OpticalGeneration ( ComputeFromMonochromaticSource () TimeDependence (WaveTime (1e-9 1e-9) WaveTSigma 50e-12) # 脉冲起始与高斯展宽 Scaling 0 Excitation ( Wavelength 0.8 *um Intensity 0.06 *W/cm2 Window(L1) (Origin(200,0) XDirection(1,0,0) Line(Dx10) Theta0) OpticalSolver (OptBeam (LayerStackExtraction ( WindowName L1 WindowPosition Center Mode ElementWise ComplexRefractiveIndex (WavelengthDep (real imag))))) ) ) )光生载流子的产生率不是均匀的一个常数而是“吸收功率密度 ÷ 光子能量”。吸收功率密度由波长、入射强度和材料层堆叠的复折射率共同决定LayerStackExtraction加OptBeam就是沿着光束方向做一层一层的吸收计算。所以Wavelength决定的是分布形状短波长在表面就被吃掉长波长穿得深Intensity只决定总量。只调强度不改波长空间分布永远不对。时间维度上WaveTSigma 50e-12是一个 50 ps 量级的高斯脉冲。瞬态求解器的MaxStep 1e-10100 fs远小于这个宽度才能把脉冲的上升沿采样出来MinStep 1e-17是兜底防止在载流子密度剧烈变化的瞬间步长被砍到零而卡死。重离子那段的FinalTime 35e-9配MaxStep 1e-10同样是因为径迹沉积发生在皮秒量级时间步必须比物理事件的时间尺度小一到两个数量级否则瞬态电流尖峰会被直接抹平——曲线看起来平滑其实什么都没捕到。4.4 陷阱谱辐照损伤为什么用“注量 × 引入率”标定Physics (material Silicon) { Traps ( (Acceptor Level EnergyMid 0.42 fromCondBand Conc 2.3226E15 Randomize 0.29 eXsection 9.5E-15 hXsection 9.5E-14) # 上面这行 Conc Fluence * 1.1613 (Acceptor Level EnergyMid 0.46 fromCondBand Conc 1.8E15 Randomize 0.23 eXsection 5E-15 hXsection 5E-14) (Donor Level EnergyMid 0.36 fromValBand Conc 1.8E15 Randomize 0.31 eXsection 3.23E-13 hXsection 3.23E-14) ) } Physics (MaterialInterface Oxide/Silicon) { Charge (Conc 1.5e11) }辐照在禁带里引入深能级载流子的复合与产生改走这些能级等效于把 SRH 的能级参数换成了“辐照后的谱”。EnergyMid指定能级位置相对导带或价带Conc是体浓度注释里那行Fluence * 1.1613才是真正的机制说明陷阱浓度由辐照注量乘一个引入率换算而来引入率是材料与辐射类型的经验常数。这意味着整套模型的可迁移性完全取决于这个系数换个工艺条件就得重新标定。Randomize给能级一个随机展宽让同一批陷阱的能级不是一根线而是一个分布。这个细节的价值在于单一能级会让复合率在某个偏压区出现尖锐变化展宽之后曲线更接近实测的“弥散”过渡。eXsection和hXsection差一个数量级说明该陷阱对电子和空穴的俘获概率不对称直接影响少子寿命的估算。界面那行Charge (Conc 1.5e11)是面电荷单位是 cm⁻²不是体浓度。混淆这两个单位阈值电压会偏出好几伏。5. 避坑语句能跑通和结果可信之间差的那几步能跑出曲线不算成功能解释曲线才算。下面几条是这套语句最容易翻车的地方几乎都表现为“形状对、数值不对”。5.1 声明阶段的坑现象加了Avalanche之后击穿电压比预期高出一大截而且曲线拐点很软。原因通常是两条结区网格太粗电场峰值被数值扩散抹平同时没开AvalDerivatives碰撞电离项在雅可比里只是一阶近似迭代走不到真实的击穿点。解决办法是先把结区网格加密到德拜长度量级再确认AvalDerivatives已开最后把QuasiStationary的MaxStep降一个数量级重新逼近拐点。现象二维仿真的电流和实测差了两个数量级但曲线形状几乎重合。原因几乎总是AreaFactor 1也就是在看“每微米宽度”的电流而实测是整个器件的电流。解决办法是先把AreaFactor设成实际宽度与 1 μm 的比值再重新校准。顺序不能反否则前面所有拟合参数都是按错的口径定的。现象BreakCriteria写了但完全不生效电流跑飞也没被拦。原因是Contact名字和Electrode里定义的名字不一致判据挂在一个不存在的接触上。这种错不会报错只会安静地失效。5.2 求解阶段的坑现象QuasiStationary在某个偏置点反复回退MinStep触底后报错终止。原因多半是该点发生了强非线性跃变耗尽区穿通、击穿而MaxStep太大加上Notdamped太大牛顿一步跨过去落在无解区。解决办法是先把MaxStep降到原值的十分之一观察是否还回退仍然失败就把Notdamped调小让阻尼早点介入。现象Transient第一步就发散日志里残差暴涨。原因是初值不是稳态解直接把零偏置的初始场丢进瞬态。正确做法是先QuasiStationary把偏置扫到目标工作点再切换Transient让时域推进从一个自洽的状态出发。现象瞬态能跑完但峰值电流比预期低一大截波形被削平。原因通常是MaxStep大于物理事件的时间尺度脉冲上升沿被跨过去。检查方法很直接把MaxStep减半重跑看峰值是否变化还在变说明步长没收敛。5.3 输出与单位的坑现象想扫 C-V但输出文件里每个偏置点只有一组数看不出频率依赖。原因是StartFrequency和EndFrequency设成了同一个值且NumberOfPoints 1只算了一个频点。做 C-V 时这是正确的做频率扫描时才需要拉开区间并给点数。现象改了PicoCoulomb相关的量之后瞬态电流整体差了六个数量级。原因是把 LET 的能量量纲和电荷量纲混用而器件尺寸又是微米量级量纲缺口正好落在 10⁶ 附近。核对方式是拿一个已知解析解的简单结构均匀电场下的平板先跑一遍确认量级对得上再上真实结构。现象Plot段列了一堆物理量输出文件里却没有。原因是Plot文件只在Solve里带Plot语句的偏置点或时刻写出光在顶层Plot段声明不够。想固定输出就在QuasiStationary或Transient的求解体里显式给Plot的范围和前缀。6. 取舍模型开得多和结果可信之间没有必然关系把所有模型都打开是最容易做也最没用的选择。每开一个模型就引入至少一个未知参数SRH 给了掺杂依赖但寿命的掺杂曲线本身要标定Auger 系数有文献值但器件里实际的有效系数受工艺影响陷阱模型的能级、浓度、俘获截面加起来接近十个自由度。而一份 I-V 曲线能提供的信息量远不足以唯一确定这些参数。参数不可辨识的时候拟合得再好也只是把误差摊到了别处。可行的做法是分层校准一次只放开一组参数阶段固定住的部分用来拟合的数据主要标定对象1关掉雪崩、带间隧穿、陷阱正向 I-V迁移率模型与 SRH 寿命2继承阶段 1反向 I-V产生寿命、禁带变窄3继承阶段 2C-V掺杂剖面、界面电荷4继承阶段 3瞬态波形陷阱谱、LET 响应参数顺序本身有物理依据正向电流由串联电阻和注入决定对迁移率和寿命最敏感反向电流由产生过程决定电容对掺杂和界面电荷最敏感瞬态响应才对时间常数敏感。反过来做用瞬态数据去定迁移率参数之间会互相补偿拟合出花来也对不上另一组实验。数值参数同样要做收敛检查它们也是“模型”。MaxStep减半、网格整体加密一倍看关键结果击穿电压、峰值电流、电容平台值变化是否小于 1%如果还在变说明之前的结果是数值伪影而不是物理结论。这个检查花不了多少机时但能挡掉大部分“仿真和实测对不上”的争论。我自己的习惯是任何一份新的语句文件先关掉所有非必要模型跑一遍收敛性检查确认网格和步长都收敛再逐个把模型打开每开一个就重跑同一组数据对比增量。有一次跳过这步直接在一份继承来的文件上改LET_f结果连续三天都在怀疑物理模型最后发现是MinStep设得太大、峰值根本没被采到。从那以后我每次改瞬态参数都强制先做一次步长减半的对照跑再看曲线。希望这能帮你省掉那三天。本文还有配套的精品资源点击获取
简介《各种语句sentaurus.docx》是一份面向微电子、纳米电子与光电子器件仿真学习者的 Sentaurus 语句参考文档适合已具备 TCAD 基础、需要查阅或复用仿真命令的读者。内容围绕 I-V、C-V 曲线计算辐照损伤与重粒子入射以及瞬态仿真展开涵盖电流电压文件、电极与物理模型、Math 数值方法、Solve 求解流程、AC 耦合电容提取、HeavyIon 重离子撞击和激光光照等典型语句片段。读者可据此快速定位关键参数设置、模型选择及输出控制方式理解器件仿真输入文件的组织逻辑并为搭建、调试和优化仿真 deck 提供可直接对照的模板减少从零编写命令时的试错成本。压缩包共 1 个 docx 文件约 16KB体积轻便便于检索与摘录。该文档已有 2038 人学习适合作为 Sentaurus 常用语句的速查资料和课程实验、科研建模的辅助参考。1. 一行 Physics 开关为什么曲线形状和收敛性会绑在一起在 Sentaurus 这类 TCAD 工具里最反直觉的一件事是物理模型和数值收敛并不是两件事。你在 Physics 段里多加一个Avalanche (UniBo Eparallel)曲线上会多出一条击穿分支但同一个文件里 Math 段的Notdamped、BreakCriteria和 Solve 段的Coupled括号顺序决定了这条分支是被算出来还是让整个仿真在偏置扫到一半时崩掉。很多人第一次遇到发散会本能地怀疑模型不对真正的原因却常是求解路径没铺好。这套语句真正值得拆开看的地方是它把“写哪个方程”“按什么顺序解”“结果往哪儿落”三件事拆成了独立层级每层都能单独调试。搞懂分层比背参数有用。2. 声明与装配Device/File/Electrode/Physics 为什么必须分层写拿到一份能跑的 Sentaurus 语句第一眼看到的是几个平级的块File、Electrode、Physics、Math、Solve。看起来像一份配置清单实际是一次自上而下的装配先声明数据从哪来、往哪去再声明边界条件再声明方程里带哪些物理项最后才声明怎么解。顺序不能乱因为后面的块依赖前面的声明。2.1 File 块为什么输入输出要分成四条独立的流一个二极管 I-V 求解的开头通常是这样File { Grid Diode_msh.tdr # 输入网格 材料区 掺杂分布由结构编辑器生成 Current Diode_IV.plt # 输出电极上的标量电压/电流序列一条曲线一个文件 Plot Diode_IV.tdr # 输出全场量快照电势、载流子密度、电场矢量 Output Diode_IV.log # 输出求解日志迭代次数、残差、步长回退记录 }这四条流分开写不是为了好看。Current存的是电极上的一维时间序列几百个点几十 KBPlot存的是全场量每个节点十几个物理量动辄几百 MB。把它们塞进一个文件代价是每次只想看一条曲线都要读全场数据。更重要的是第四条Output日志是唯一能解释“为什么没算完”的地方。求解器不会告诉你物理哪里错了它只会告诉你第几步、残差多少、步长被砍到多少。把它当成黑匣子里的飞行记录仪出问题第一件事就是翻它而不是改模型。工程实践中还有一个细节Plot里的量由Plot段二次筛选File只负责命名目的地。两者分离的设计让你可以重跑后处理而不重算——前提是NoOverwrite这类开关用对否则文件名被覆盖后悔药就没了。2.2 Physics 块每个模型开关都是对方程组的一次修改这是整套语句里信息密度最高的地方Physics { AreaFactor 1 # 二维截面结果折算为等效三维时的宽度因子 Temperature 300 # 晶格温度单位 K不是载流子温度 Mobility ( DopingDependence # 迁移率随掺杂浓度退化 eHighFieldSaturation # 电子高场速度饱和 hHighFieldSaturation # 空穴高场速度饱和 EnormalCarrier # 表面法向电场引起的额外退化 CarrierCarrierScattering # 载流子-载流子散射 ) Recombination ( SRH (DopingDependence) # 深能级复合寿命随掺杂下降 Auger (withGeneration) # 三载流子过程高注入下主导 Avalanche (UniBo Eparallel) # 碰撞电离按平行电场驱动 Band2Band (Hurkx) # 带间隧穿折算成等效产生率 EffectiveIntrinsicDensity (OldSlotboom) # 重掺杂下的禁带变窄 ) }关键在于这些括号不是“精度档位”每一个都往漂移-扩散方程组里加一项加项就改雅可比矩阵的结构和条件数。SRH让复合率依赖寿命 τ(DopingDependence)再把 τ 写成掺杂浓度的函数于是高掺杂区的少子寿命被压到纳秒级反向饱和电流随之变化Auger (withGeneration)在 SRH 之外再加一个 ∝n²p 的复合通道withGeneration让它允许把多余能量转给另一个载流子去激发这在重掺杂发射区里不能省。Avalanche (UniBo Eparallel)更典型。碰撞电离率对电场是指数敏感的Eparallel指定驱动电场只取平行于电流方向的分量——如果取全电场模在结边缘或异质结界面上电场方向不纯会凭空触发大量电离击穿电压直接跑偏。这类模型细节在曲线上表现为“击穿点提前了 20%”不去拆开关根本定位不到。Band2Band (Hurkx)和EffectiveIntrinsicDensity (OldSlotboom)属于另一类前者把带间隧穿折算成局部产生率避免对隧穿路径做能量积分后者用经典模型描述禁带变窄让本征载流子浓度在重掺杂区抬升直接改变 PN 结的注入比。它们的共同点是——不加结果可能看起来“还算像”加上才和实测对得上但加了以后收敛变难。这就是物理精度与数值鲁棒性的第一次正面冲突。2.3 Electrode 与 AreaFactor二维结果怎么变成“等效三维”Electrode { NameElectrodeN Voltage 0.0 }只做两件事给接触起名给初值。名字非常重要——Plot段、BreakCriteria段、混合模式里的电路节点映射全都靠这个名字索引。名字写错语句照样能跑只是断点判据失效电流跑飞也没人拦。AreaFactor 1是另一个容易忽视的换算口。二维仿真的电流天然是“每单位深度”的电流A/μm真实器件的电流要乘以宽度。设成 1 意味着你在看“每微米宽度”的结果。忘了改这一项得到的电流可能和实测差两个数量级而曲线形状完全正确——这是最阴的一类错因为形状对人会先怀疑物理参数。提示二维结构里做校准之前先把 AreaFactor 按实际器件宽度设好再把实验结果换算到同一口径。顺序反了后面所有拟合都是白做。3. 求解序列Coupled 的括号顺序比模型精度更决定成败模型定了接下来是解。漂移-扩散是一组强非线性耦合方程泊松方程管静电势两条连续性方程管电子和空穴任何一条单独解都得不到有意义的电流。真正的难点不是“要不要耦合”而是“从什么状态开始耦合”。3.1 从 Poisson 到 Electron/Hole/Poisson 的渐进耦合Solve { Coupled (Iterations50) { Poisson } # 第一步只要静电势 Coupled (Iterations15) { Hole Poisson } # 加空穴连续性 Coupled (Iterations15) { Electron Hole Poisson } # 再加电子连续性 QuasiStationary ( InitialStep 1e-6 MaxStep 0.01 MinStep 1e-9 Goal { NameElectrodeN Voltage 1000 } Plot { Range (0 1) Intervals 2 } ) { Coupled { Hole Electron Poisson } # 每个偏置点内部仍是全耦合牛顿 Plot (FilePrefixIV_ Time(0.01;0.05;0.1;0.5) NoOverwrite) } }前三条Coupled是渐进式装配先解泊松拿到电势分布把它当作空穴方程的初值再拿两者当电子方程的初值。零初值下直接上三变量耦合载流子密度在数量级上跨越十几个量级牛顿迭代的线性系统条件数极差几乎必然发散。分步走相当于给每一步一个“差不多合理”的起点代价是多花几次迭代。注意每步的Iterations不一样泊松 50 次后面 15 次。泊松方程相对好解但初值最差给宽一点连续性方程已有前一步的场做引导15 次够用。这个数字不是固定的网格加密后经常要调大。QuasiStationary里的Goal是“扫到某个电压”求解器自己决定中间踩哪些点。InitialStep、MaxStep、MinStep控制自适应步长步长太大越过击穿拐点牛顿在无解区反复回退步长太小几百个偏置点白烧机时。MaxStep 0.01意味着单步不超过 0.01 V在 1000 V 的扫描范围里意味着至少 10 万个点——这显然不是常规用法更像是先给一个保守上限确保能跑完。真实项目里通常会先粗后细用大MaxStep跑通全程看形状再在小范围内收紧步长抓拐点。3.2 Math 块收敛失败的病根常在导数Math { # Cylindrical # 柱坐标对称注释掉说明当前是笛卡尔坐标 Method Pardiso # 稀疏直接法求解线性系统 Number_of_threads 4 Stacksize ... # 递归求解时防止栈溢出 ExtrapolateDerivatives # 导数用外推近似省一次函数求值 AvalDerivatives # 碰撞电离项参与雅可比必须给真导数 RelErrControl # 误差控制基于相对量 Iterations 15 # 单个偏置点允许的最大牛顿迭代 Notdamped 60 # 前 60 次迭代不做步长阻尼 BreakCriteria { Current (ContactElectrodeN maxval 1e-8) } }Pardiso是直接法分解矩阵一次解到底鲁棒但内存随节点数非线性增长。强非线性问题击穿、雪崩、辐照瞬态一般不用迭代法因为迭代法省下的内存会被“反复不收敛”加倍还回去。ExtrapolateDerivatives和AvalDerivatives是一对反向操作。牛顿法需要雅可比矩阵里的 ∂F/∂x有些项的解析导数代价很高外推近似能省时间但会破坏二次收敛碰撞电离系数是电场的指数函数导数错一点迭代方向就差很远所以必须用真实导数。这两个开关的取舍原则很朴素线性或弱非线性项可以外推指数型强非线性项不行。Notdamped 60和Iterations 15之间其实有个矛盾单点最多 15 次迭代可阻尼要到 60 次才启用。等于说在这套设置下阻尼形同虚设。这是血泪经验的产物——对某些结构早开阻尼会让迭代在局部极小值附近磨蹭反而不如放开步长让它冲过去。但这个值必须配合BreakCriteria使用否则冲过头就是发散。BreakCriteria是整套语句里最像后悔药的一行接触电流超过 1e-8 就停。它不改变物理只在数值失控时把进程从跑飞状态拉回来并留下一个明确的终止位置方便你回去看那一步的日志。4. 外部激励AC 小信号、重离子与激光为什么共用一套骨架I-V 是直流稳态C-V 是频率域的小扰动重离子和激光是时域的大扰动看起来是三套完全不同的东西。但它们在语句结构上高度同构都要先有一个收敛的稳态解再把外部激励作用上去再沿着某个自变量电压、频率、时间推进。理解这个同构就不会把它们当成三份互不相干的模板去背。4.1 混合模式与 ACCoupledC-V 为什么必须挂在直流解上System { Diode diodesystem (ElectrodeNfront ElectrodeP0) # 器件电极映射到电路节点 Vsource_pset vn (front 0) { dc 0 } # 电压源dc 是被扫描的量 } File { Output Diode_CV_1kHz.log ACExtract Diode_CV_AC_1kHz.plt } Solve { Coupled (Iterations50) { Poisson } Coupled (Iterations15) { Hole Poisson } Coupled (Iterations15) { Electron Hole Poisson } Coupled (Iterations15) { Electron Hole Poisson Contact } QuasiStationary (InitialStep1e-6 MaxStep1e-2 MinStep1e-7 Increment2 Decrement4 Goal { Parameter vn.dc Voltage 100 }) { ACCoupled (StartFrequency1e3 EndFrequency1e3 NumberOfPoints1 Decade Iterations15 Node(front) ACMethodBlocked ACSubMethod(diodesystem)ParDiSo) { Poisson Electron Hole Contact Circuit } } }小信号分析的数学前提是线性化把解写成 x x₀ Re(δx·e^{jωt})代入原方程后丢掉二次以上的小量得到一组复数线性方程。x₀ 就是直流工作点没有它就没有扰动可加。所以ACCoupled必须嵌在QuasiStationary内部每扫到一个偏置点先解直流再在这个点上加正弦扰动最后从复导纳里取出电容。这也是为什么ACExtract输出的是“电容-电压”而不是“电容-频率”——频率在这里只是求导的算子不是扫描变量。StartFrequency EndFrequency 1e3加NumberOfPoints 1等于只算 1 kHz 单点。这是 C-V 的常见用法固定频率扫偏置。若真想做频率扫描这两个值要拉开并给点数此时Decade才起作用。ACMethod Blocked把器件方程和电路方程分块求解ACSubMethod(diodesystem) ParDiSo指定该子系统的线性解法。分块的代价是丢掉器件-电路之间的一部分交叉项在“一个器件加一个理想源”的拓扑下几乎无影响一旦电路侧有多个强耦合元件就要谨慎评估。Node(front)指明扰动从哪个电路节点注入Contact出现在未知量列表里意味着接触被当作额外变量处理对接触电阻和边界通量的处理更稳。低频 1 kHz 这个选择本身也有物理含义频率足够低时深能级陷阱跟得上信号测到的是“含陷阱响应”的电容频率高到陷阱来不及响应电容会掉一截。用单频算 C-V等于默认陷阱响应可以忽略或已经包含在直流模型里。4.2 HeavyIonLET 波形怎么变成体电荷Physics { HeavyIon ( Direction (0, 1) # 入射方向 Location (200, 0) # 入射点 Time 1e-9 # 入射时刻 Length [ ] # 径迹长度区间 Wt_hi [ ] # 径迹横向分布的权重 LET_f [0 8.7e-6 8.7e-6 0] # 线性能量转移随时间的分段线性波形 Gaussian # 横向按高斯分布铺开 PicoCoulomb # 电荷量纲避免与能量量纲混用 ) }机制链条很短但每一步都可能出错LET(t) 给出单位长度上的能量沉积率除以产生一对电子空穴所需的平均能量硅约 3.6 eV得到单位长度、单位时间的载流子对产生率再按Gaussian的横向分布摊到径迹周围的网格节点上转成体产生率塞进连续性方程。LET_f那四个数就是波形的骨架起、升、平台、落。平台段的斜率决定了瞬态电流尖峰有多陡也直接决定后续Transient的MinStep要压到多小。Location和Direction决定径迹穿过哪些区域——打在结区还是衬底产生的电荷收集效率完全不同这是单粒子效应仿真里最敏感的几何参数。Gaussian不能省。径迹在物理上有横向宽度如果强行让所有电荷落在一列节点上结果会强烈依赖网格网格换一套峰值电流就变一次。加上横向分布后电荷铺到若干列节点上网格依赖性显著下降。注意PicoCoulomb是量纲声明不是单位换算开关。LET 数据、器件尺寸、体电荷密度三者口径不一致时得到的瞬态电流可能差好几个数量级而且曲线形状还是对的。4.3 激光与瞬态步长为什么 MinStep 要压到 1e-17Optics ( OpticalGeneration ( ComputeFromMonochromaticSource () TimeDependence (WaveTime (1e-9 1e-9) WaveTSigma 50e-12) # 脉冲起始与高斯展宽 Scaling 0 Excitation ( Wavelength 0.8 *um Intensity 0.06 *W/cm2 Window(L1) (Origin(200,0) XDirection(1,0,0) Line(Dx10) Theta0) OpticalSolver (OptBeam (LayerStackExtraction ( WindowName L1 WindowPosition Center Mode ElementWise ComplexRefractiveIndex (WavelengthDep (real imag))))) ) ) )光生载流子的产生率不是均匀的一个常数而是“吸收功率密度 ÷ 光子能量”。吸收功率密度由波长、入射强度和材料层堆叠的复折射率共同决定LayerStackExtraction加OptBeam就是沿着光束方向做一层一层的吸收计算。所以Wavelength决定的是分布形状短波长在表面就被吃掉长波长穿得深Intensity只决定总量。只调强度不改波长空间分布永远不对。时间维度上WaveTSigma 50e-12是一个 50 ps 量级的高斯脉冲。瞬态求解器的MaxStep 1e-10100 fs远小于这个宽度才能把脉冲的上升沿采样出来MinStep 1e-17是兜底防止在载流子密度剧烈变化的瞬间步长被砍到零而卡死。重离子那段的FinalTime 35e-9配MaxStep 1e-10同样是因为径迹沉积发生在皮秒量级时间步必须比物理事件的时间尺度小一到两个数量级否则瞬态电流尖峰会被直接抹平——曲线看起来平滑其实什么都没捕到。4.4 陷阱谱辐照损伤为什么用“注量 × 引入率”标定Physics (material Silicon) { Traps ( (Acceptor Level EnergyMid 0.42 fromCondBand Conc 2.3226E15 Randomize 0.29 eXsection 9.5E-15 hXsection 9.5E-14) # 上面这行 Conc Fluence * 1.1613 (Acceptor Level EnergyMid 0.46 fromCondBand Conc 1.8E15 Randomize 0.23 eXsection 5E-15 hXsection 5E-14) (Donor Level EnergyMid 0.36 fromValBand Conc 1.8E15 Randomize 0.31 eXsection 3.23E-13 hXsection 3.23E-14) ) } Physics (MaterialInterface Oxide/Silicon) { Charge (Conc 1.5e11) }辐照在禁带里引入深能级载流子的复合与产生改走这些能级等效于把 SRH 的能级参数换成了“辐照后的谱”。EnergyMid指定能级位置相对导带或价带Conc是体浓度注释里那行Fluence * 1.1613才是真正的机制说明陷阱浓度由辐照注量乘一个引入率换算而来引入率是材料与辐射类型的经验常数。这意味着整套模型的可迁移性完全取决于这个系数换个工艺条件就得重新标定。Randomize给能级一个随机展宽让同一批陷阱的能级不是一根线而是一个分布。这个细节的价值在于单一能级会让复合率在某个偏压区出现尖锐变化展宽之后曲线更接近实测的“弥散”过渡。eXsection和hXsection差一个数量级说明该陷阱对电子和空穴的俘获概率不对称直接影响少子寿命的估算。界面那行Charge (Conc 1.5e11)是面电荷单位是 cm⁻²不是体浓度。混淆这两个单位阈值电压会偏出好几伏。5. 避坑语句能跑通和结果可信之间差的那几步能跑出曲线不算成功能解释曲线才算。下面几条是这套语句最容易翻车的地方几乎都表现为“形状对、数值不对”。5.1 声明阶段的坑现象加了Avalanche之后击穿电压比预期高出一大截而且曲线拐点很软。原因通常是两条结区网格太粗电场峰值被数值扩散抹平同时没开AvalDerivatives碰撞电离项在雅可比里只是一阶近似迭代走不到真实的击穿点。解决办法是先把结区网格加密到德拜长度量级再确认AvalDerivatives已开最后把QuasiStationary的MaxStep降一个数量级重新逼近拐点。现象二维仿真的电流和实测差了两个数量级但曲线形状几乎重合。原因几乎总是AreaFactor 1也就是在看“每微米宽度”的电流而实测是整个器件的电流。解决办法是先把AreaFactor设成实际宽度与 1 μm 的比值再重新校准。顺序不能反否则前面所有拟合参数都是按错的口径定的。现象BreakCriteria写了但完全不生效电流跑飞也没被拦。原因是Contact名字和Electrode里定义的名字不一致判据挂在一个不存在的接触上。这种错不会报错只会安静地失效。5.2 求解阶段的坑现象QuasiStationary在某个偏置点反复回退MinStep触底后报错终止。原因多半是该点发生了强非线性跃变耗尽区穿通、击穿而MaxStep太大加上Notdamped太大牛顿一步跨过去落在无解区。解决办法是先把MaxStep降到原值的十分之一观察是否还回退仍然失败就把Notdamped调小让阻尼早点介入。现象Transient第一步就发散日志里残差暴涨。原因是初值不是稳态解直接把零偏置的初始场丢进瞬态。正确做法是先QuasiStationary把偏置扫到目标工作点再切换Transient让时域推进从一个自洽的状态出发。现象瞬态能跑完但峰值电流比预期低一大截波形被削平。原因通常是MaxStep大于物理事件的时间尺度脉冲上升沿被跨过去。检查方法很直接把MaxStep减半重跑看峰值是否变化还在变说明步长没收敛。5.3 输出与单位的坑现象想扫 C-V但输出文件里每个偏置点只有一组数看不出频率依赖。原因是StartFrequency和EndFrequency设成了同一个值且NumberOfPoints 1只算了一个频点。做 C-V 时这是正确的做频率扫描时才需要拉开区间并给点数。现象改了PicoCoulomb相关的量之后瞬态电流整体差了六个数量级。原因是把 LET 的能量量纲和电荷量纲混用而器件尺寸又是微米量级量纲缺口正好落在 10⁶ 附近。核对方式是拿一个已知解析解的简单结构均匀电场下的平板先跑一遍确认量级对得上再上真实结构。现象Plot段列了一堆物理量输出文件里却没有。原因是Plot文件只在Solve里带Plot语句的偏置点或时刻写出光在顶层Plot段声明不够。想固定输出就在QuasiStationary或Transient的求解体里显式给Plot的范围和前缀。6. 取舍模型开得多和结果可信之间没有必然关系把所有模型都打开是最容易做也最没用的选择。每开一个模型就引入至少一个未知参数SRH 给了掺杂依赖但寿命的掺杂曲线本身要标定Auger 系数有文献值但器件里实际的有效系数受工艺影响陷阱模型的能级、浓度、俘获截面加起来接近十个自由度。而一份 I-V 曲线能提供的信息量远不足以唯一确定这些参数。参数不可辨识的时候拟合得再好也只是把误差摊到了别处。可行的做法是分层校准一次只放开一组参数阶段固定住的部分用来拟合的数据主要标定对象1关掉雪崩、带间隧穿、陷阱正向 I-V迁移率模型与 SRH 寿命2继承阶段 1反向 I-V产生寿命、禁带变窄3继承阶段 2C-V掺杂剖面、界面电荷4继承阶段 3瞬态波形陷阱谱、LET 响应参数顺序本身有物理依据正向电流由串联电阻和注入决定对迁移率和寿命最敏感反向电流由产生过程决定电容对掺杂和界面电荷最敏感瞬态响应才对时间常数敏感。反过来做用瞬态数据去定迁移率参数之间会互相补偿拟合出花来也对不上另一组实验。数值参数同样要做收敛检查它们也是“模型”。MaxStep减半、网格整体加密一倍看关键结果击穿电压、峰值电流、电容平台值变化是否小于 1%如果还在变说明之前的结果是数值伪影而不是物理结论。这个检查花不了多少机时但能挡掉大部分“仿真和实测对不上”的争论。我自己的习惯是任何一份新的语句文件先关掉所有非必要模型跑一遍收敛性检查确认网格和步长都收敛再逐个把模型打开每开一个就重跑同一组数据对比增量。有一次跳过这步直接在一份继承来的文件上改LET_f结果连续三天都在怀疑物理模型最后发现是MinStep设得太大、峰值根本没被采到。从那以后我每次改瞬态参数都强制先做一次步长减半的对照跑再看曲线。希望这能帮你省掉那三天。本文还有配套的精品资源点击获取
 与 Runner loop_factory 实战)
