简介这份文档面向微电子、纳米电子与光电子器件仿真学习者及科研人员围绕Sentaurus常用语句与典型仿真流程展开可用于快速查阅和复用器件电学特性分析配置。压缩包内仅含1个docx文件约16KB篇幅精简但覆盖I-V曲线、C-V曲线、辐照损伤与重粒子入射、瞬态仿真等模块并串联File、Electrode、Physics、Math、Plot、Solve等常用段落。内容涉及电流电压定义与输出、迁移率和复合模型设置以及电容混合模式AC耦合仿真、重离子入射的LET与时间参数、激光光照瞬态求解等关键语句并展示重离子撞击与激光瞬态的求解配置思路可帮助读者理解参数之间的对应关系。已有2036人学习下载适合需要对照Sentaurus命令结构、排查语法与组织参数的读者作为案头参考也可用于搭建二极管I-V/C-V及瞬态实验的语句模板减少重复配置时间为后续仿真实验提供可借鉴的语句组织方式。1. 一份 docx 装下五种仿真Sentaurus 语句笔记能解决的问题我第一次拆这种资料时打开文件就想关掉全是File{}、Physics{}、Solve{}这样的命令块堆在一起没有目录没有配套网格也没有一句注释说明谁在前谁在后。但只要真正跑过几个器件就会明白Sentaurus 这类 TCAD 仿真的难点并不在模型写成什么样而在于分寸——网格里的接触叫什么名字、迁移率模型开哪几个、准静态扫描的步长给多大才不发散。这些在手册里能写三百页在一份实际能跑通的输入文件里往往只有五十行。这份资源把 I-V 曲线、C-V 曲线、重离子撞击、激光注入和辐照陷阱五类 SDevice 输入整理到了一起覆盖二极管类器件的直流、小信号和瞬态场景。适合手上有.tdr网格、准备开始跑、并且已经被“不收敛”折腾过几轮的人。它的价值也就在这里——照着改一遍比从空文件开始试要省下不少时间。2. 命令骨架拆解五类仿真共用的四段结构2.1 File / Electrode / Physics / Solve谁在前谁依赖谁无论是 I-V 还是瞬态SDevice 输入文件的骨架都是由同样的四段拼起来的。先看一份可以直接跑的最简骨架File { Grid Diode_msh.tdr # 网格文件由前端的网格生成器导出 Current Diode_IV.plt # 电学曲线终端电压/电流逐点写入 Plot Diode_IV.tdr # 场量快照体积远大于 plt Output Diode_IV.log # 求解日志收敛问题全在这里 } Electrode { { NameElectrodeN Voltage0.0 } # 名字必须和网格里的接触一一对应 } Physics { AreaFactor 1 Temperature 300 Mobility(DopingDependence eHighFieldSaturation hHighFieldSaturation Enormal CarrierCarrierScattering) Recombination(SRH(DopingDependence) Auger(withGeneration) Avalanche(UniBo Eparallel) Band2Band(Hurkx) EffectiveIntrinsicDensity(OldSlotboom)) }File段的先后顺序其实是有讲究的。Grid是唯一一个“缺了就什么都跑不了”的条目名字对不上会直接报找不到文件连日志都写不完整Output反而是最不该省的一项因为后面所有的调整参数、降低步长、切换求解器全部要靠这个 log 文件里的残差曲线来判断。Current和Plot的区别经常被新手混淆.plt是标量时间序列几 MB画 I-V 就靠它.tdr是场分布快照每个输出点都会存一份完整网格数据跑一次瞬态能堆到几个 GB。Electrode段里的Voltage只是初值真正的扫描目标在Solve里定义。这个“初值”非常重要所有耦合求解的第一次迭代都从这个偏置出发如果给了一个离平衡态很远的初值第一步就容易出现“not converged”的报错。Physics段是整份文件里信息密度最高的地方这里逐项看一下AreaFactor 1二维切片仿真默认按单位深度计算。如果器件实际宽度是 100 µm而网格只画了 1 µm 宽的截面电流结果必须乘以这个因子才能对应真实器件。这个数值写错所有电流都会差几个数量级而曲线形状却完全正确很难发现。Temperature 300单位是开尔文。我见过写成27的跑出来的曲线平得像一条直线因为 27 K 下本征载流子浓度几乎为零。Mobility(...)DopingDependence是电离杂质散射eHighFieldSaturation/hHighFieldSaturation分别是电子和空穴的强场速度饱和Enormal是横向电场依赖反型层里不开这个模型沟道电流会虚高CarrierCarrierScattering是高注入下才需要补的载流子间散射。Recombination(...)SRH(DopingDependence)是带掺杂依赖的肖克利-里德-霍尔复合Auger(withGeneration)里的withGeneration表示反向偏置下也计算产生项漏电分析必须打开Avalanche(UniBo Eparallel)用的是 UniBo 雪崩模型Eparallel表示按电场平行分量做判据击穿分析离不开它Band2Band(Hurkx)是带间隧穿EffectiveIntrinsicDensity(OldSlotboom)是禁带变窄模型。这五行基本就是“一个硅器件该有的物理”。少开一项曲线在某个区间就会不对劲多开一项收敛难度就会立刻上升。我一般的习惯是先把雪崩和带间隧穿关掉跑通低偏置段确认接触和掺杂没错再把它们加回来。2.2 五类仿真分别落在哪些段这五种仿真看着差别很大位置实际上就集中在两三个段里仿真类型新增或改动的段关键语句主要产物I-V 直流SolveQuasiStationaryGoal{Name...}.plt曲线C-V 小信号Device/System/SolveVsource_pset、ACCoupledACExtract文件重离子瞬态PhysicsSolveHeavyIon(...)、Transienttrans_前缀曲线激光注入OpticsSolveOpticalGeneration、WaveTime光生载流子分布辐照陷阱Physics(material...)Traps(...)、Charge(...)陷阱占据、复合率把表竖着看会发现一个规律File和Electrode几乎不动Solve才是每类仿真真正的分化点。这也解释了为什么这些语句可以被整理进同一个文档——它们共享同一个底座差异全是“怎么把偏置推上去”。2.3 单位与坐标最不值得踩却最多人踩的地方SDevice 的单位体系并不统一它是按物理量分治的温度用 K长度用 µm浓度用 cm⁻³时间用 s电流用 A。这四套单位混在一个文件里写的时候根本看不出来。坐标是另一处隐患。Location (200, 0)这样的写法必须和网格坐标系严格一致——网格生成器导出时如果做过平移或缩放这个点可能就跑到器件外面去了结果是“仿真正常结束但物理事件没有任何效果”。提示重离子和激光窗口都用绝对坐标定位改网格尺寸后这两个坐标必须同步改否则事件会打在空气里。3. I-V 与 C-V 两条链路从准静态扫描到 AC 小信号3.1 准静态扫 I-V为什么要先解三次先看 I-V 的求解段这是整份资料里最完整的一段Solve { # 三步暖启动Poisson 单独解 → 加空穴 → 全耦合 Coupled(Iterations50) { Poisson } Coupled(Iterations15) { Hole Poisson } Coupled(Iterations15) { Electron Hole Poisson } QuasiStationary( InitialStep 1e-6 MaxStep 0.01 MinStep 1e-9 Goal { NameElectrodeN Voltage1000 } # 扫描终点 Plot { Range(0 1) Intervals2 } # 归一化扫描进度上出快照 ) { Coupled { Hole Electron Poisson } } Plot(FilePrefixIV_ Time(0.01;0.05;0.1;0.5) NoOverwrite) }开头那三次Coupled是整套流程里最容易被删掉、删掉后又最容易出问题的部分。器件初始状态默认是“中性近似”直接上全耦合电子和空穴的初值都偏离自洽解牛顿迭代的雅可比矩阵条件数很差第一步就可能失败。先单独解Poisson只有电势一个未知量收敛很快再依次引入空穴和电子等于给后续迭代铺了三级台阶。QuasiStationary的三个步长参数都表示比例不是电压值。MaxStep 0.01的意思是单步最多推进总扫描范围的 1%也就是至少 100 步走完全程。MinStep 1e-9是底线自适应步长一旦压到这个值以下还无法收敛求解器就会报错退出。很多“扫到某个电压就断”的情况根源并不在物理模型而是MinStep设得太小、求解器在硬撑。Goal { NameElectrodeN Voltage1000 }里的 1000 V 是演示用的量级。真实硅二极管正向 1 V 左右电流就已经很大了直接扫到 1000 V 会在中间某处彻底崩溃。我一般先给 5 V 看整体形状确认没有出现异常再往上加。Plot { Range(0 1) Intervals2 }里的Range是归一化扫描进度从 0 到 1Intervals2表示切成两段、输出三个快照。而末尾的Plot(FilePrefixIV_ Time(0.01;0.05;0.1;0.5))在准静态语境下Time同样指归一化进度不是真实的物理时间。这两处写法不同、含义相近比较容易把人绕晕。Math段还有一条容易被忽略的提前停止条件Math { Method Pardiso Number_of_threads 4 Extrapolate Derivatives AvalDerivatives RelErrControl Iterations 15 Notdamped 60 BreakCriteria { Current(ContactElectrodeN maxval1e-8) } }BreakCriteria的意思是当该接触的电流超过 1e-8 时就停掉扫描。这在雪崩区特别有用——电流一旦进入指数上升段再往下扫只会浪费时间而且大概率不收敛。Notdamped 60表示前 60 次迭代不做阻尼让求解器先尝试完整的牛顿步之后才引入阻尼因子这个值调小了收敛慢调大了容易过冲。3.2 Mixed-Mode ACCoupledC-V 为什么要两段式C-V 不能靠直流扫描算必须走小信号路径。资料里用的是 Mixed-Mode 写法把器件和电路分开描述Device Diode { Electrode { { NameElectrodeN Voltage0 } { NameElectrodeP Voltage0 } } File { Grid Diode_msh.tdr Current Diode_CV_1kHz.plt Plot Diode_CV_1kHz.tdr } Physics { /* 与 I-V 相同的一套模型 */ } Plot { eDensity hDensity Potential SpaceCharge } } System { Diode diodesystem (ElectrodeNfront ElectrodeP0) Vsource_pset vn (front 0) { dc0 } } 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( InitialStep 1e-6 MaxStep 1e-2 MinStep 1e-7 Increment 2 Decrement 4 Goal { Parameter vn.dc Voltage 100 } ) { ACCoupled( StartFrequency 1e3 EndFrequency 1e3 NumberOfPoints 1 Decade Iterations 15 Node(front) ACMethod Blocked ACSubMethod(diodesystem) ParDiSo ) { Poisson Electron Hole Contact Circuit } } }Device段描述器件本体网格、物理、输出System段描述外部电路。(ElectrodeNfront ElectrodeP0)把两个电极分别绑定到电路节点front和地。Vsource_pset vn (front 0)是一个分段线性电压源dc0是初值后面Goal { Parameter vn.dc Voltage 100 }扫的就是这个源而不是电极本身——这是 Mixed-Mode 和纯器件仿真的关键区别。ACCoupled才是真正的电容提取动作。StartFrequency EndFrequency 1e3配NumberOfPoints 1意思就是只算 1 kHz 这一个频点。做 C-V 单频扫描足够了如果要做频率色散分析把EndFrequency拉开、让NumberOfPoints覆盖几个十倍频程就行。ACMethod Blocked加ACSubMethod(diodesystem) ParDiSo是把交流求解和直流求解拆开让电路子系统单独走并行对角求解。器件规模一大、电路节点一多这个选项能显著缩短单点求解时间属于规模上来之后必须知道的一个开关。Node(front)指定小信号加在哪个节点上。漏了这行求解器不知道该往哪注入扰动输出的导纳矩阵就无从谈起。还有两个参数值得单独说Increment 2和Decrement 4。前者表示收敛顺利时步长放大 2 倍后者表示不收敛时收缩到 1/4。这两个值决定自适应步长的激进程度反偏段电容变化平缓可以放宽一点正偏接近开启时变化剧烈收紧一点更稳。3.3 结果文件怎么换算ACExtract输出的是复导纳不是电容。这点非常关键文件里的值需要按 Y G jωC 反解C Im(Y) / ω其中 ω 2πf。单频 1 kHz 时 ω ≈ 6283.2忘了除以 ω 直接当电容用数值会差三四个数量级。.plt文件是纯文本表头会写明每一列对应的物理量前几列一般是电极电压/电流或时间。用脚本提取时不要靠列序硬编码先读表头再定位换一台机器或换一个版本都能少出问题。4. 瞬态与辐照重离子、激光、陷阱三套物理4.1 HeavyIon时间轴必须和 Transient 对齐重离子撞击是这几个场景里最容易“跑完但什么也没发生”的一个。物理段写法如下Physics { /* 迁移率、复合模型同前 */ HeavyIon ( Direction (0,1) # 入射方向沿 y 轴进入器件 Location (200, 0) # 入射点必须落在网格范围内 Time 1e-9 # 撞击时刻 Length [ ... ] # 径迹长度 Wt_hi [ ... ] # 径迹特征宽度 LET_f [0 8.7e-6 8.7e-6 0] # LET 随时间的四点包络 Gaussian PicoCoulomb ) }Location和Direction决定径迹在器件里的几何位置Time决定它什么时候发生。LET_f的四个值是这条径迹的强度包络零 → 上升 → 平台 → 归零配合Gaussian的时间展宽形成一个脉冲。PicoCoulomb切换的是 LET 的单位体系——这时候LET_f按 pC/µm 给而不是常见的 MeV·cm²/mg。这两种单位的换算跟材料密度和电子空穴对生成能有关量级差得很远直接照抄数值基本都会出错。Wt_hi需要和网格尺度对齐。如果网格在入射点附近的单元尺寸远大于这个特征宽度沉积的能量会被均摊到一个形状不规则的体积上得到的瞬态电流波形会明显失真——峰值偏低、拖尾偏长。经验做法是先在入射路径上做局部加密把单元尺寸压到Wt_hi的几分之一。配套的求解段Solve { NewCurrentPrefix trans_ Transient ( InitialTime 0 FinalTime 35e-9 MinStep 1e-17 MaxStep 1e-10 ) { Coupled { Poisson Electron Hole Circuit } Plot (FilePrefixTransHI_ Time(0.5e-9;1e-9;2e-9;5e-9;10e-9) NoOverwrite) } }MinStep 1e-17看起来夸张但径迹中心附近的载流子梯度确实能在几个飞秒内变化好几个数量级步长压不下去牛顿迭代就会直接失败。MaxStep 1e-10保证在 1 ns 到 2 ns 这段关键窗口里有足够的采样点不会把脉冲峰值跳过去。NewCurrentPrefix trans_给本次运行的所有曲线加前缀。这是重跑实验时的一个实用习惯——不改前缀第二次跑的结果会直接追加到第一次的.plt里画出来是两条叠在一起的曲线很容易误判成“物理异常”。Plot(...Time(...))里的Time在瞬态里是真实物理时间和准静态里的归一化含义不同这一点切场景的时候特别容易混淆。4.2 激光注入Optics 段和它依赖的层栈激光激励走的是完全不同的入口不在Physics里而是独立的Optics段Optics ( 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 (Dx 10) Theta 0 * Angle from positive y-axis OpticalSolver ( OptBeam ( LayerStackExtraction ( WindowName L1 WindowPosition Center Mode ElementWise ComplexRefractiveIndex (WavelengthDep (real imag)) ) ) ) ) ) ) )ComputeFromMonochromaticSource表示光源是单色光光强分布由后面的Excitation给出。Wavelength 0.8 * um落在近红外硅在这个波段吸收深度是几十微米量级所以窗口深度要给够。TimeDependence决定脉冲形状。WaveTime (1e-9 1e-9)是脉冲起止时刻两个值相同表示这是一个极窄的窗口WaveTSigma 50e-12给高斯展宽。这套组合的整体效果是脉冲中心在 1 ns 附近宽度约 50 ps。Scaling 0表示不额外缩放光强全部按Intensity给的值来。这一项写成非零值光强会被乘上去结果和预期差很多倍而且日志里不会有任何提示。LayerStackExtraction是容易被跳过但一旦跳过就出错的部分。它按窗口从器件表面往下逐层拆分光在每一层里按复折射率衰减。ComplexRefractiveIndex(WavelengthDep(real imag))说明折射率的实部和虚部都要从材料参数里读——如果材料库里的光学参数不全这一步会退回默认值算出来的光生分布和平面的均匀分布差不多看起来“能跑”其实完全不对。对应的瞬态段比重离子长得多NewCurrentPrefix opt_ Transient (InitialTime0 FinalTime60e-9 MinStep1e-17 MaxStep1e-10) { Coupled { Poisson Electron Hole Circuit } }FinalTime 60e-9是为了让光生载流子有足够时间被电场扫走、电流回到基线。只跑到 10 ns 的话看得到的都是上升沿恢复过程完全缺失。4.3 Traps陷阱浓度是从注量折算出来的辐照损伤的核心在Physics的材料限定写法里Physics (materialSilicon) { Traps ( (Acceptor Level EnergyMid0.42 fromCondBand Conc2.3226E15 Randomize0.29 eXsection9.5E-15 hXsection9.5E-14) (Acceptor Level EnergyMid0.46 fromCondBand Conc1.8E15 Randomize0.23 eXsection5E-15 hXsection5E-14) (Donor Level EnergyMid0.36 fromValBand Conc1.8E15 Randomize0.31 eXsection3.23E-13 hXsection3.23E-14) ) } Physics (MaterialInterfaceOxide/Silicon) { Charge (Conc1.5e11) }Acceptor/Donor在这里指陷阱类型不是掺杂类型——它们决定陷阱带电时是捕获电子还是空穴。EnergyMid是能级位置fromCondBand/fromValBand决定这个位置是相对导带底还是价带顶量起。0.42 eV 和 0.46 eV 两个受主能级加上 0.36 eV 的施主能级构成一套典型的辐照缺陷能谱。Conc的值不是凭直觉定的。原文件里留了一条注释Conc Fluence * 1.1613。也就是说2.3226e15 是从 2e15 的注量乘上一个 1.1613 的折算系数得来的。这个系数跟离子种类、能量、材料都有关直接套用会系统性地偏但至少它把“注量 → 缺陷浓度”这条链路显式写出来了比直接填一个数字好得多。Randomize是能级的随机展宽用来模拟非晶化之后能级的弥散。这个值不为零时实际参与计算的不是单一能级而是一个以EnergyMid为中心的分布。eXsection/hXsection是电子和空穴的俘获截面单位是 cm²。注意第三组施主陷阱的电子截面 3.23e-13 比前两组大了将近两个数量级这会显著影响少子寿命也是这类模型里最容易调出不同结果的地方。界面那行的Charge(Conc1.5e11)是固定氧化物电荷单位是 cm⁻²。它和体陷阱不是一回事体陷阱参与 SRH 复合动力学界面电荷只是给沟道加了一个固定偏置不随偏置变化。把materialSilicon这个限定去掉陷阱会被应用到所有材料上包括氧化层结果就是栅氧里凭空多出一堆复合中心C-V 曲线的积累区会出现解释不了的凹陷。5. 避坑与排查五条血泪记录5.1 温度写成摄氏度曲线全是平的现象I-V 曲线从头到尾是一条水平线电流量级在 1e-30 附近。原因Temperature 27被当成了开尔文。27 K 下本征载流子浓度比 300 K 低了十几个数量级。解决SDevice 里的温度一律用 K。写完之后横向扫一眼整个文件凡是出现 27、25、85 这些数字的地方都要回头核对单位425 K 这种“看着像温度”的值反而容易忽略。5.2 电极名和网格对不上现象求解能跑起来但所有电极电流都是零或者日志里出现找不到接触的警告。原因Electrode { NameElectrodeN }里的名字和.tdr网格里的接触名不完全一致。大小写、下划线、数字后缀都可能差一个字符。解决先打开网格文件把接触名列表打印出来对着抄一遍别凭记忆写。改过网格之后这一项必须重新核对这也是唯一一个“网格换了就必须同步改”的地方。5.3 重离子跑完波形一点反应都没有现象瞬态正常结束电流曲线是一条平滑的基线看不到任何脉冲。原因HeavyIon的Time落在Transient的[InitialTime, FinalTime]区间之外或者MaxStep太大求解器直接跳过了撞击时刻附近那一段。解决先把FinalTime拉到Time之后足够远再把MaxStep临时压到 1e-11 做一次验证确认波形出现之后再放宽。另外Location也检查一遍坐标落在网格外面同样不会有任何效果而且不会报错。5.4 扫到某个电压就断日志里全是步长告警现象QuasiStationary到某个偏置电压时退出日志末尾提示步长已经到下限。原因三种可能混在一起——MaxStep给得太小导致总步数不够雪崩模型在电流还很低的区间被过早打开导致矩阵条件数变差BreakCriteria没设求解器在雪崩区强行往下推进。解决先加上BreakCriteria { Current(ContactElectrodeN maxval1e-8) }让它在该停的地方停下来再把Avalanche关掉跑一遍确认低偏置段是否本来就能收敛最后才考虑把MaxStep放大到 0.05 减少总步数。5.5 输出文件越跑越大数据还叠在一起现象第二次运行时.plt里的数据点数是第一次的两倍画出来两条曲线重合在一起或.tdr单个文件几百 MB。原因Current文件默认是追加写入的同名文件不会自动清空Plot段如果开了全部场量每个输出点都会存一份完整数据。解决每次重跑前用NewCurrentPrefix换个前缀或者手动删掉旧文件Plot只保留分析真正需要的量通常eDensity、hDensity、Potential、SpaceCharge加一两个电流矢量就够其余全删。NoOverwrite能挡住同名快照被覆盖但挡不住追加两个机制要配合用。6. 二次开发把散装语句变成参数化模板这份资料的原始形态是一堆并列的语句片段直接复制粘贴能跑但没法批量做事。真正想拿它做点东西第一步是抽公共部分第二步是参数化第三步是建立回归基线。抽公共部分很直接把File的四个路径、Physics的迁移率和复合模型、Math的求解器配置抽到一个模板文件里用占位符标出每次会变的地方。这里有个必须注意的细节——不能用{ }做占位符因为 SDevice 语法本身就大量使用大括号替换脚本会把结构一起改坏。用NAME这种双符号定界更安全from pathlib import Path TPL Path(sdevice_template.cmd).read_text(encodingutf-8) def render(**kw): out TPL for k, v in kw.items(): out out.replace(%s % k, str(v)) return out # 同一套物理扫三档 LET输出三份独立命令文件 for let in (8.7e-6, 4.0e-5, 1.0e-4): cmd render( GRIDDiode_msh.tdr, VSTOP5.0, LETlet, PREFIXhi_%.1e % let, ) Path(run_%.1e.cmd % let).write_text(cmd, encodingutf-8)模板里对应位置写成Grid GRID、LET_f [0 LET LET 0]、NewCurrentPrefix PREFIX。这样换网格、换注量、换扫描终点都不用再动模板本体。参数化之后一定要有回归基线否则改一个参数跑出来的结果对不对自己心里没底。我一般的做法是固定三个检查点正向 0.8 V 的电流量级、反向 5 V 的漏电流、零偏下的结电容。三个都对得上才认为这轮改动没有破坏基础物理。#!/usr/bin/env bash # 粗筛脚本只看关键点的量级不看完整曲线 set -euo pipefail sdevice run_base.cmd /dev/null 21 # 列索引要按 .plt 表头来定这里假设第 1 列电压、第 2 列电流 awk NF2 $100.8 {printf fwd0.8V %.3e\n, $2} Diode_IV.plt awk NF2 $10-5.0 {printf rev-5V %.3e\n, $2} Diode_IV.pltawk里那个列索引是最容易骗人的地方——不同版本、不同输出选项下列的顺序会变所以脚本开头必须先去.plt表头确认一遍别硬编码。再往上一层就是批量参数扫描。把render()外面套一层循环遍历LET、Fluence、Wavelength三组变量每次生成独立目录、独立前缀跑完统一收集.plt里的峰值电流和恢复时间。这时候NewCurrentPrefix和NoOverwrite的价值就体现出来了——没有它们几十个任务的结果会互相覆盖事后根本分不清哪条曲线对应哪组参数。最后说一个我自己吃过亏的地方。早期我总想一次把所有物理都开满——雪崩、带间隧穿、陷阱、界面电荷全开——然后对着不收敛的日志反复调MaxStep和Notdamped。后来改成先把模型一层层剥掉跑通最简版本再逐个加回来每次加一项就跑一遍基线问题出在哪一层立刻就能定位。从那以后我每次改Physics段都强制走一遍“关掉一半 → 跑通 → 加回一项 → 再跑通”的流程比对着日志猜要快得多。希望这个习惯也能帮到你。本文还有配套的精品资源点击获取
简介这份文档面向微电子、纳米电子与光电子器件仿真学习者及科研人员围绕Sentaurus常用语句与典型仿真流程展开可用于快速查阅和复用器件电学特性分析配置。压缩包内仅含1个docx文件约16KB篇幅精简但覆盖I-V曲线、C-V曲线、辐照损伤与重粒子入射、瞬态仿真等模块并串联File、Electrode、Physics、Math、Plot、Solve等常用段落。内容涉及电流电压定义与输出、迁移率和复合模型设置以及电容混合模式AC耦合仿真、重离子入射的LET与时间参数、激光光照瞬态求解等关键语句并展示重离子撞击与激光瞬态的求解配置思路可帮助读者理解参数之间的对应关系。已有2036人学习下载适合需要对照Sentaurus命令结构、排查语法与组织参数的读者作为案头参考也可用于搭建二极管I-V/C-V及瞬态实验的语句模板减少重复配置时间为后续仿真实验提供可借鉴的语句组织方式。1. 一份 docx 装下五种仿真Sentaurus 语句笔记能解决的问题我第一次拆这种资料时打开文件就想关掉全是File{}、Physics{}、Solve{}这样的命令块堆在一起没有目录没有配套网格也没有一句注释说明谁在前谁在后。但只要真正跑过几个器件就会明白Sentaurus 这类 TCAD 仿真的难点并不在模型写成什么样而在于分寸——网格里的接触叫什么名字、迁移率模型开哪几个、准静态扫描的步长给多大才不发散。这些在手册里能写三百页在一份实际能跑通的输入文件里往往只有五十行。这份资源把 I-V 曲线、C-V 曲线、重离子撞击、激光注入和辐照陷阱五类 SDevice 输入整理到了一起覆盖二极管类器件的直流、小信号和瞬态场景。适合手上有.tdr网格、准备开始跑、并且已经被“不收敛”折腾过几轮的人。它的价值也就在这里——照着改一遍比从空文件开始试要省下不少时间。2. 命令骨架拆解五类仿真共用的四段结构2.1 File / Electrode / Physics / Solve谁在前谁依赖谁无论是 I-V 还是瞬态SDevice 输入文件的骨架都是由同样的四段拼起来的。先看一份可以直接跑的最简骨架File { Grid Diode_msh.tdr # 网格文件由前端的网格生成器导出 Current Diode_IV.plt # 电学曲线终端电压/电流逐点写入 Plot Diode_IV.tdr # 场量快照体积远大于 plt Output Diode_IV.log # 求解日志收敛问题全在这里 } Electrode { { NameElectrodeN Voltage0.0 } # 名字必须和网格里的接触一一对应 } Physics { AreaFactor 1 Temperature 300 Mobility(DopingDependence eHighFieldSaturation hHighFieldSaturation Enormal CarrierCarrierScattering) Recombination(SRH(DopingDependence) Auger(withGeneration) Avalanche(UniBo Eparallel) Band2Band(Hurkx) EffectiveIntrinsicDensity(OldSlotboom)) }File段的先后顺序其实是有讲究的。Grid是唯一一个“缺了就什么都跑不了”的条目名字对不上会直接报找不到文件连日志都写不完整Output反而是最不该省的一项因为后面所有的调整参数、降低步长、切换求解器全部要靠这个 log 文件里的残差曲线来判断。Current和Plot的区别经常被新手混淆.plt是标量时间序列几 MB画 I-V 就靠它.tdr是场分布快照每个输出点都会存一份完整网格数据跑一次瞬态能堆到几个 GB。Electrode段里的Voltage只是初值真正的扫描目标在Solve里定义。这个“初值”非常重要所有耦合求解的第一次迭代都从这个偏置出发如果给了一个离平衡态很远的初值第一步就容易出现“not converged”的报错。Physics段是整份文件里信息密度最高的地方这里逐项看一下AreaFactor 1二维切片仿真默认按单位深度计算。如果器件实际宽度是 100 µm而网格只画了 1 µm 宽的截面电流结果必须乘以这个因子才能对应真实器件。这个数值写错所有电流都会差几个数量级而曲线形状却完全正确很难发现。Temperature 300单位是开尔文。我见过写成27的跑出来的曲线平得像一条直线因为 27 K 下本征载流子浓度几乎为零。Mobility(...)DopingDependence是电离杂质散射eHighFieldSaturation/hHighFieldSaturation分别是电子和空穴的强场速度饱和Enormal是横向电场依赖反型层里不开这个模型沟道电流会虚高CarrierCarrierScattering是高注入下才需要补的载流子间散射。Recombination(...)SRH(DopingDependence)是带掺杂依赖的肖克利-里德-霍尔复合Auger(withGeneration)里的withGeneration表示反向偏置下也计算产生项漏电分析必须打开Avalanche(UniBo Eparallel)用的是 UniBo 雪崩模型Eparallel表示按电场平行分量做判据击穿分析离不开它Band2Band(Hurkx)是带间隧穿EffectiveIntrinsicDensity(OldSlotboom)是禁带变窄模型。这五行基本就是“一个硅器件该有的物理”。少开一项曲线在某个区间就会不对劲多开一项收敛难度就会立刻上升。我一般的习惯是先把雪崩和带间隧穿关掉跑通低偏置段确认接触和掺杂没错再把它们加回来。2.2 五类仿真分别落在哪些段这五种仿真看着差别很大位置实际上就集中在两三个段里仿真类型新增或改动的段关键语句主要产物I-V 直流SolveQuasiStationaryGoal{Name...}.plt曲线C-V 小信号Device/System/SolveVsource_pset、ACCoupledACExtract文件重离子瞬态PhysicsSolveHeavyIon(...)、Transienttrans_前缀曲线激光注入OpticsSolveOpticalGeneration、WaveTime光生载流子分布辐照陷阱Physics(material...)Traps(...)、Charge(...)陷阱占据、复合率把表竖着看会发现一个规律File和Electrode几乎不动Solve才是每类仿真真正的分化点。这也解释了为什么这些语句可以被整理进同一个文档——它们共享同一个底座差异全是“怎么把偏置推上去”。2.3 单位与坐标最不值得踩却最多人踩的地方SDevice 的单位体系并不统一它是按物理量分治的温度用 K长度用 µm浓度用 cm⁻³时间用 s电流用 A。这四套单位混在一个文件里写的时候根本看不出来。坐标是另一处隐患。Location (200, 0)这样的写法必须和网格坐标系严格一致——网格生成器导出时如果做过平移或缩放这个点可能就跑到器件外面去了结果是“仿真正常结束但物理事件没有任何效果”。提示重离子和激光窗口都用绝对坐标定位改网格尺寸后这两个坐标必须同步改否则事件会打在空气里。3. I-V 与 C-V 两条链路从准静态扫描到 AC 小信号3.1 准静态扫 I-V为什么要先解三次先看 I-V 的求解段这是整份资料里最完整的一段Solve { # 三步暖启动Poisson 单独解 → 加空穴 → 全耦合 Coupled(Iterations50) { Poisson } Coupled(Iterations15) { Hole Poisson } Coupled(Iterations15) { Electron Hole Poisson } QuasiStationary( InitialStep 1e-6 MaxStep 0.01 MinStep 1e-9 Goal { NameElectrodeN Voltage1000 } # 扫描终点 Plot { Range(0 1) Intervals2 } # 归一化扫描进度上出快照 ) { Coupled { Hole Electron Poisson } } Plot(FilePrefixIV_ Time(0.01;0.05;0.1;0.5) NoOverwrite) }开头那三次Coupled是整套流程里最容易被删掉、删掉后又最容易出问题的部分。器件初始状态默认是“中性近似”直接上全耦合电子和空穴的初值都偏离自洽解牛顿迭代的雅可比矩阵条件数很差第一步就可能失败。先单独解Poisson只有电势一个未知量收敛很快再依次引入空穴和电子等于给后续迭代铺了三级台阶。QuasiStationary的三个步长参数都表示比例不是电压值。MaxStep 0.01的意思是单步最多推进总扫描范围的 1%也就是至少 100 步走完全程。MinStep 1e-9是底线自适应步长一旦压到这个值以下还无法收敛求解器就会报错退出。很多“扫到某个电压就断”的情况根源并不在物理模型而是MinStep设得太小、求解器在硬撑。Goal { NameElectrodeN Voltage1000 }里的 1000 V 是演示用的量级。真实硅二极管正向 1 V 左右电流就已经很大了直接扫到 1000 V 会在中间某处彻底崩溃。我一般先给 5 V 看整体形状确认没有出现异常再往上加。Plot { Range(0 1) Intervals2 }里的Range是归一化扫描进度从 0 到 1Intervals2表示切成两段、输出三个快照。而末尾的Plot(FilePrefixIV_ Time(0.01;0.05;0.1;0.5))在准静态语境下Time同样指归一化进度不是真实的物理时间。这两处写法不同、含义相近比较容易把人绕晕。Math段还有一条容易被忽略的提前停止条件Math { Method Pardiso Number_of_threads 4 Extrapolate Derivatives AvalDerivatives RelErrControl Iterations 15 Notdamped 60 BreakCriteria { Current(ContactElectrodeN maxval1e-8) } }BreakCriteria的意思是当该接触的电流超过 1e-8 时就停掉扫描。这在雪崩区特别有用——电流一旦进入指数上升段再往下扫只会浪费时间而且大概率不收敛。Notdamped 60表示前 60 次迭代不做阻尼让求解器先尝试完整的牛顿步之后才引入阻尼因子这个值调小了收敛慢调大了容易过冲。3.2 Mixed-Mode ACCoupledC-V 为什么要两段式C-V 不能靠直流扫描算必须走小信号路径。资料里用的是 Mixed-Mode 写法把器件和电路分开描述Device Diode { Electrode { { NameElectrodeN Voltage0 } { NameElectrodeP Voltage0 } } File { Grid Diode_msh.tdr Current Diode_CV_1kHz.plt Plot Diode_CV_1kHz.tdr } Physics { /* 与 I-V 相同的一套模型 */ } Plot { eDensity hDensity Potential SpaceCharge } } System { Diode diodesystem (ElectrodeNfront ElectrodeP0) Vsource_pset vn (front 0) { dc0 } } 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( InitialStep 1e-6 MaxStep 1e-2 MinStep 1e-7 Increment 2 Decrement 4 Goal { Parameter vn.dc Voltage 100 } ) { ACCoupled( StartFrequency 1e3 EndFrequency 1e3 NumberOfPoints 1 Decade Iterations 15 Node(front) ACMethod Blocked ACSubMethod(diodesystem) ParDiSo ) { Poisson Electron Hole Contact Circuit } } }Device段描述器件本体网格、物理、输出System段描述外部电路。(ElectrodeNfront ElectrodeP0)把两个电极分别绑定到电路节点front和地。Vsource_pset vn (front 0)是一个分段线性电压源dc0是初值后面Goal { Parameter vn.dc Voltage 100 }扫的就是这个源而不是电极本身——这是 Mixed-Mode 和纯器件仿真的关键区别。ACCoupled才是真正的电容提取动作。StartFrequency EndFrequency 1e3配NumberOfPoints 1意思就是只算 1 kHz 这一个频点。做 C-V 单频扫描足够了如果要做频率色散分析把EndFrequency拉开、让NumberOfPoints覆盖几个十倍频程就行。ACMethod Blocked加ACSubMethod(diodesystem) ParDiSo是把交流求解和直流求解拆开让电路子系统单独走并行对角求解。器件规模一大、电路节点一多这个选项能显著缩短单点求解时间属于规模上来之后必须知道的一个开关。Node(front)指定小信号加在哪个节点上。漏了这行求解器不知道该往哪注入扰动输出的导纳矩阵就无从谈起。还有两个参数值得单独说Increment 2和Decrement 4。前者表示收敛顺利时步长放大 2 倍后者表示不收敛时收缩到 1/4。这两个值决定自适应步长的激进程度反偏段电容变化平缓可以放宽一点正偏接近开启时变化剧烈收紧一点更稳。3.3 结果文件怎么换算ACExtract输出的是复导纳不是电容。这点非常关键文件里的值需要按 Y G jωC 反解C Im(Y) / ω其中 ω 2πf。单频 1 kHz 时 ω ≈ 6283.2忘了除以 ω 直接当电容用数值会差三四个数量级。.plt文件是纯文本表头会写明每一列对应的物理量前几列一般是电极电压/电流或时间。用脚本提取时不要靠列序硬编码先读表头再定位换一台机器或换一个版本都能少出问题。4. 瞬态与辐照重离子、激光、陷阱三套物理4.1 HeavyIon时间轴必须和 Transient 对齐重离子撞击是这几个场景里最容易“跑完但什么也没发生”的一个。物理段写法如下Physics { /* 迁移率、复合模型同前 */ HeavyIon ( Direction (0,1) # 入射方向沿 y 轴进入器件 Location (200, 0) # 入射点必须落在网格范围内 Time 1e-9 # 撞击时刻 Length [ ... ] # 径迹长度 Wt_hi [ ... ] # 径迹特征宽度 LET_f [0 8.7e-6 8.7e-6 0] # LET 随时间的四点包络 Gaussian PicoCoulomb ) }Location和Direction决定径迹在器件里的几何位置Time决定它什么时候发生。LET_f的四个值是这条径迹的强度包络零 → 上升 → 平台 → 归零配合Gaussian的时间展宽形成一个脉冲。PicoCoulomb切换的是 LET 的单位体系——这时候LET_f按 pC/µm 给而不是常见的 MeV·cm²/mg。这两种单位的换算跟材料密度和电子空穴对生成能有关量级差得很远直接照抄数值基本都会出错。Wt_hi需要和网格尺度对齐。如果网格在入射点附近的单元尺寸远大于这个特征宽度沉积的能量会被均摊到一个形状不规则的体积上得到的瞬态电流波形会明显失真——峰值偏低、拖尾偏长。经验做法是先在入射路径上做局部加密把单元尺寸压到Wt_hi的几分之一。配套的求解段Solve { NewCurrentPrefix trans_ Transient ( InitialTime 0 FinalTime 35e-9 MinStep 1e-17 MaxStep 1e-10 ) { Coupled { Poisson Electron Hole Circuit } Plot (FilePrefixTransHI_ Time(0.5e-9;1e-9;2e-9;5e-9;10e-9) NoOverwrite) } }MinStep 1e-17看起来夸张但径迹中心附近的载流子梯度确实能在几个飞秒内变化好几个数量级步长压不下去牛顿迭代就会直接失败。MaxStep 1e-10保证在 1 ns 到 2 ns 这段关键窗口里有足够的采样点不会把脉冲峰值跳过去。NewCurrentPrefix trans_给本次运行的所有曲线加前缀。这是重跑实验时的一个实用习惯——不改前缀第二次跑的结果会直接追加到第一次的.plt里画出来是两条叠在一起的曲线很容易误判成“物理异常”。Plot(...Time(...))里的Time在瞬态里是真实物理时间和准静态里的归一化含义不同这一点切场景的时候特别容易混淆。4.2 激光注入Optics 段和它依赖的层栈激光激励走的是完全不同的入口不在Physics里而是独立的Optics段Optics ( 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 (Dx 10) Theta 0 * Angle from positive y-axis OpticalSolver ( OptBeam ( LayerStackExtraction ( WindowName L1 WindowPosition Center Mode ElementWise ComplexRefractiveIndex (WavelengthDep (real imag)) ) ) ) ) ) ) )ComputeFromMonochromaticSource表示光源是单色光光强分布由后面的Excitation给出。Wavelength 0.8 * um落在近红外硅在这个波段吸收深度是几十微米量级所以窗口深度要给够。TimeDependence决定脉冲形状。WaveTime (1e-9 1e-9)是脉冲起止时刻两个值相同表示这是一个极窄的窗口WaveTSigma 50e-12给高斯展宽。这套组合的整体效果是脉冲中心在 1 ns 附近宽度约 50 ps。Scaling 0表示不额外缩放光强全部按Intensity给的值来。这一项写成非零值光强会被乘上去结果和预期差很多倍而且日志里不会有任何提示。LayerStackExtraction是容易被跳过但一旦跳过就出错的部分。它按窗口从器件表面往下逐层拆分光在每一层里按复折射率衰减。ComplexRefractiveIndex(WavelengthDep(real imag))说明折射率的实部和虚部都要从材料参数里读——如果材料库里的光学参数不全这一步会退回默认值算出来的光生分布和平面的均匀分布差不多看起来“能跑”其实完全不对。对应的瞬态段比重离子长得多NewCurrentPrefix opt_ Transient (InitialTime0 FinalTime60e-9 MinStep1e-17 MaxStep1e-10) { Coupled { Poisson Electron Hole Circuit } }FinalTime 60e-9是为了让光生载流子有足够时间被电场扫走、电流回到基线。只跑到 10 ns 的话看得到的都是上升沿恢复过程完全缺失。4.3 Traps陷阱浓度是从注量折算出来的辐照损伤的核心在Physics的材料限定写法里Physics (materialSilicon) { Traps ( (Acceptor Level EnergyMid0.42 fromCondBand Conc2.3226E15 Randomize0.29 eXsection9.5E-15 hXsection9.5E-14) (Acceptor Level EnergyMid0.46 fromCondBand Conc1.8E15 Randomize0.23 eXsection5E-15 hXsection5E-14) (Donor Level EnergyMid0.36 fromValBand Conc1.8E15 Randomize0.31 eXsection3.23E-13 hXsection3.23E-14) ) } Physics (MaterialInterfaceOxide/Silicon) { Charge (Conc1.5e11) }Acceptor/Donor在这里指陷阱类型不是掺杂类型——它们决定陷阱带电时是捕获电子还是空穴。EnergyMid是能级位置fromCondBand/fromValBand决定这个位置是相对导带底还是价带顶量起。0.42 eV 和 0.46 eV 两个受主能级加上 0.36 eV 的施主能级构成一套典型的辐照缺陷能谱。Conc的值不是凭直觉定的。原文件里留了一条注释Conc Fluence * 1.1613。也就是说2.3226e15 是从 2e15 的注量乘上一个 1.1613 的折算系数得来的。这个系数跟离子种类、能量、材料都有关直接套用会系统性地偏但至少它把“注量 → 缺陷浓度”这条链路显式写出来了比直接填一个数字好得多。Randomize是能级的随机展宽用来模拟非晶化之后能级的弥散。这个值不为零时实际参与计算的不是单一能级而是一个以EnergyMid为中心的分布。eXsection/hXsection是电子和空穴的俘获截面单位是 cm²。注意第三组施主陷阱的电子截面 3.23e-13 比前两组大了将近两个数量级这会显著影响少子寿命也是这类模型里最容易调出不同结果的地方。界面那行的Charge(Conc1.5e11)是固定氧化物电荷单位是 cm⁻²。它和体陷阱不是一回事体陷阱参与 SRH 复合动力学界面电荷只是给沟道加了一个固定偏置不随偏置变化。把materialSilicon这个限定去掉陷阱会被应用到所有材料上包括氧化层结果就是栅氧里凭空多出一堆复合中心C-V 曲线的积累区会出现解释不了的凹陷。5. 避坑与排查五条血泪记录5.1 温度写成摄氏度曲线全是平的现象I-V 曲线从头到尾是一条水平线电流量级在 1e-30 附近。原因Temperature 27被当成了开尔文。27 K 下本征载流子浓度比 300 K 低了十几个数量级。解决SDevice 里的温度一律用 K。写完之后横向扫一眼整个文件凡是出现 27、25、85 这些数字的地方都要回头核对单位425 K 这种“看着像温度”的值反而容易忽略。5.2 电极名和网格对不上现象求解能跑起来但所有电极电流都是零或者日志里出现找不到接触的警告。原因Electrode { NameElectrodeN }里的名字和.tdr网格里的接触名不完全一致。大小写、下划线、数字后缀都可能差一个字符。解决先打开网格文件把接触名列表打印出来对着抄一遍别凭记忆写。改过网格之后这一项必须重新核对这也是唯一一个“网格换了就必须同步改”的地方。5.3 重离子跑完波形一点反应都没有现象瞬态正常结束电流曲线是一条平滑的基线看不到任何脉冲。原因HeavyIon的Time落在Transient的[InitialTime, FinalTime]区间之外或者MaxStep太大求解器直接跳过了撞击时刻附近那一段。解决先把FinalTime拉到Time之后足够远再把MaxStep临时压到 1e-11 做一次验证确认波形出现之后再放宽。另外Location也检查一遍坐标落在网格外面同样不会有任何效果而且不会报错。5.4 扫到某个电压就断日志里全是步长告警现象QuasiStationary到某个偏置电压时退出日志末尾提示步长已经到下限。原因三种可能混在一起——MaxStep给得太小导致总步数不够雪崩模型在电流还很低的区间被过早打开导致矩阵条件数变差BreakCriteria没设求解器在雪崩区强行往下推进。解决先加上BreakCriteria { Current(ContactElectrodeN maxval1e-8) }让它在该停的地方停下来再把Avalanche关掉跑一遍确认低偏置段是否本来就能收敛最后才考虑把MaxStep放大到 0.05 减少总步数。5.5 输出文件越跑越大数据还叠在一起现象第二次运行时.plt里的数据点数是第一次的两倍画出来两条曲线重合在一起或.tdr单个文件几百 MB。原因Current文件默认是追加写入的同名文件不会自动清空Plot段如果开了全部场量每个输出点都会存一份完整数据。解决每次重跑前用NewCurrentPrefix换个前缀或者手动删掉旧文件Plot只保留分析真正需要的量通常eDensity、hDensity、Potential、SpaceCharge加一两个电流矢量就够其余全删。NoOverwrite能挡住同名快照被覆盖但挡不住追加两个机制要配合用。6. 二次开发把散装语句变成参数化模板这份资料的原始形态是一堆并列的语句片段直接复制粘贴能跑但没法批量做事。真正想拿它做点东西第一步是抽公共部分第二步是参数化第三步是建立回归基线。抽公共部分很直接把File的四个路径、Physics的迁移率和复合模型、Math的求解器配置抽到一个模板文件里用占位符标出每次会变的地方。这里有个必须注意的细节——不能用{ }做占位符因为 SDevice 语法本身就大量使用大括号替换脚本会把结构一起改坏。用NAME这种双符号定界更安全from pathlib import Path TPL Path(sdevice_template.cmd).read_text(encodingutf-8) def render(**kw): out TPL for k, v in kw.items(): out out.replace(%s % k, str(v)) return out # 同一套物理扫三档 LET输出三份独立命令文件 for let in (8.7e-6, 4.0e-5, 1.0e-4): cmd render( GRIDDiode_msh.tdr, VSTOP5.0, LETlet, PREFIXhi_%.1e % let, ) Path(run_%.1e.cmd % let).write_text(cmd, encodingutf-8)模板里对应位置写成Grid GRID、LET_f [0 LET LET 0]、NewCurrentPrefix PREFIX。这样换网格、换注量、换扫描终点都不用再动模板本体。参数化之后一定要有回归基线否则改一个参数跑出来的结果对不对自己心里没底。我一般的做法是固定三个检查点正向 0.8 V 的电流量级、反向 5 V 的漏电流、零偏下的结电容。三个都对得上才认为这轮改动没有破坏基础物理。#!/usr/bin/env bash # 粗筛脚本只看关键点的量级不看完整曲线 set -euo pipefail sdevice run_base.cmd /dev/null 21 # 列索引要按 .plt 表头来定这里假设第 1 列电压、第 2 列电流 awk NF2 $100.8 {printf fwd0.8V %.3e\n, $2} Diode_IV.plt awk NF2 $10-5.0 {printf rev-5V %.3e\n, $2} Diode_IV.pltawk里那个列索引是最容易骗人的地方——不同版本、不同输出选项下列的顺序会变所以脚本开头必须先去.plt表头确认一遍别硬编码。再往上一层就是批量参数扫描。把render()外面套一层循环遍历LET、Fluence、Wavelength三组变量每次生成独立目录、独立前缀跑完统一收集.plt里的峰值电流和恢复时间。这时候NewCurrentPrefix和NoOverwrite的价值就体现出来了——没有它们几十个任务的结果会互相覆盖事后根本分不清哪条曲线对应哪组参数。最后说一个我自己吃过亏的地方。早期我总想一次把所有物理都开满——雪崩、带间隧穿、陷阱、界面电荷全开——然后对着不收敛的日志反复调MaxStep和Notdamped。后来改成先把模型一层层剥掉跑通最简版本再逐个加回来每次加一项就跑一遍基线问题出在哪一层立刻就能定位。从那以后我每次改Physics段都强制走一遍“关掉一半 → 跑通 → 加回一项 → 再跑通”的流程比对着日志猜要快得多。希望这个习惯也能帮到你。本文还有配套的精品资源点击获取
