简介：《分子动力学步骤说明》是一份面向分子模拟初学者与准备撰写相关论文的研究人员的系统导读资料。内容按模拟进行的线性顺序展开，从体系构建（加溶剂、带电离子）开始，依次讲解能量最小化、NVT预平衡、NPT预平衡、生产项模拟，再到数据分析中的重包装/重映像、RMSD、回旋半径、SASA等关键指标，并对ORCA量子化学包和GROMACS分子动力学软件做了简洁介绍。整个资源包仅含1个PDF文档，大小约6.87MB，结构清晰，适合对照阅读。内容中特别强调了数据分析阶段常见误区，例如不能仅凭RMSD曲线判断模拟收敛，同时以线性书写方式串联实验流程，对初次接触ORCA等专业名词的读者很友好。目前已有442人浏览学习，可帮助读者快速建立对MD全流程的整体认知，也能作为论文“方法”部分写作前的导读与查漏补缺参考。

1. 分子动力学模拟的完整链路：为什么要按这个顺序跑

拿到一个对接好的蛋白-配体复合物，导师只丢下一句“跑个50纳秒的分子动力学看看稳定性”，第一次接触MD的人通常会懵在原地：先加盒子还是先加离子？要不要做能量最小化？NVT和NPT能不能换顺序？这份《分子动力学步骤说明》PDF解决的问题就是把这些事按线性流程理顺：体系构建、能量最小化、NVT预平衡、NPT预平衡、生产项模拟、数据分析，一环扣一环，每一步在GROMACS里对应什么命令、失败时该看哪个指标、论文里该怎么写，都对照着给了示例。我按这个顺序跑完一个体系后发现，真正省时间的不是某个命令技巧，而是它把“为什么NVT必须先于NPT”“为什么RMSD平台不等于收敛”这些概念讲透了，后续做生产模拟基本没走回头路。适合准备在论文里加MD验证的研究生，也适合第一次完整搭建GROMACS模拟流程的从业者。

2. 体系构建与能量最小化：pdb2gmx、solvate、genion与50000步细节

2.1 输入结构准备：质子化状态先判断

体系构建的第一步不是敲命令，而是决定用哪个结构文件。通常会把对接得到的复合物PDB作为初始结构。这里有个容易被忽略的问题：PDB文件里的氢原子基本是缺失的，pdb2gmx会按预设几何构型补氢，但补成什么样取决于力场的氨基酸模板，而不是结构生物学家的pH判断。也就是说，组氨酸的HIE/HID状态、半胱氨酸是否成二硫键，这些在输入前就要确认好，pdb2gmx只是执行者，不是决策者。

常见做法是先打开PDB检查有没有非标准残基和原子名异常，再执行补氢与位置重构：

gmx pdb2gmx -f complex.pdb -o complex_processed.gro -ff amber14sb -water tip3p -ignh

这条命令里，-ff指定蛋白力场为amber14sb，它和配体端常用的gaff2力场组合很常见，兼容性也好；-water tip3p指定水模型为TIP3P；-ignh表示忽略PDB中已有的氢，统一按力场模板重新生成，避免PDB氢原子位置和模板键长冲突导致后续几何crash。跑完会生成topol.top和posre.itp，topol.top之后需要手工维护。

如果命令直接报错，八成是三类问题：PDB缺少必需的主链原子、残基原子名与力场模板对不上、存在非标准残基。这时候不要硬调PDB，正确做法是回到对接软件重新导出最完整的结构，再把配体单独摘出来走量化优化流程。

2.2 加溶剂：盒子形状与缓冲距离

补完氢、生成拓扑后，下一步是把体系放进周期性盒子。盒子形状我一般优先选十二面体，同样体积下比立方体少约30%的水分子，省计算量，而且周期性镜像距离更均匀，对后续PME静电处理友好。盒子尺寸不是越大越好，留足缓冲距离即可：

gmx editconf -f complex_processed.gro -o complex_box.gro -c -d 1.0 -bt dodecahedron gmx solvate -cp complex_box.gro -cs spc216.gro -o complex_solv.gro -p topol.top

-d 1.0表示溶质任意原子到盒子边界的最小距离为1.0nm。缓冲太小，周期性镜像里的分子会离自己太近，产生不存在的分子间作用；缓冲太大则纯浪费算力。-cs spc216.gro是预平衡好的水构型，用在TIP3P水模型上没问题，因为这里只需要水分子的初始坐标，水参数和电荷由topol.top里-water tip3p决定。

solvate跑完后topol.top里会多出SOL条目，水分子数和总原子数在输出log里可以查到。这时候我一般先打开topol.top确认蛋白链、配体、水的顺序没乱，因为后续genion依赖这个顺序来替换残基。

提示：solvate和genion都会改写topol.top，每跑一步前先备份一份，出问题有后悔药。

2.3 加离子：电中性是硬约束

生物分子模拟的体系必须电中性，否则PME的周期性求和会在长程尺度上产生不真实的静电势漂移，静电相互作用完全失真。加离子的做法是先定义一个只有离子替换规则的mdp文件，grompp生成临时tpr，再用genion把指定数量的水分子替换成钠离子或氯离子：

gmx grompp -f ions.mdp -c complex_solv.gro -p topol.top -o ions.tpr gmx genion -s ions.tpr -o complex_ions.gro -p topol.top -pname NA -nname CL -neutral

ions.mdp内容很简单，一个合法的cutoff设置加integrator=steep就够，因为grompp需要完整参数才能生成tpr。genion的-neutral参数会自动计算体系净电荷并加入对应数量的抗衡离子，不用自己数电荷。蛋白带负电时加Na+，带正电加Cl-，具体加哪种不会改变水结构，只影响论文里的描述措辞。

这里有个常见翻车点：genion随机挑水分子替换，不关心替换位置离蛋白残基多近。加完离子后我习惯用VMD看一眼离子分布，如果在口袋内部或两个连续残基之间出现明显不合理的离子，就换个随机种子重新跑一次，避免后续平衡阶段出现局部相互作用异常。

2.4 能量最小化：最陡下降法与收敛判断

构建好的体系不能直接跑生产模拟。solvate在格点上放水分子时，氢键取向是统计学的，genion替换离子是随机的，原子之间完全可能出现距离过小、范德华力巨大的情况，严重时直接crash。能量最小化的目的就是沿着原子受力方向移动位置，直到体系达到一个几何合理、能量较低的起始构型。

GROMACS里最常用的EM配置是最陡下降法：

integrator = steep nsteps = 50000 emtol = 1000.0 emstep = 0.01 nstlist = 10 cutoff-scheme = Verlet coulombtype = PME rcoulomb = 1.4 rvdw = 1.4

运行命令：

gmx grompp -f em.mdp -c complex_ions.gro -p topol.top -o em.tpr gmx mdrun -deffnm em -v

integrator=steep即最陡下降法；nsteps=50000是最大步数，并不是必须跑满，只要能量梯度低于emtol的1000.0 kJ/mol/nm或最大力收敛，mdrun会提前结束；cutoff-scheme=Verlet是GROMACS默认邻居搜索方案，配合PME处理长程静电，rcoulomb和rvdw都取1.4nm。原文里还提到cg算法，即共轭梯度，适合在steep初步收敛后做二次精细化，多数时候steep一步到位就够了。收敛后看em.log里的梯度信息，如果50000步跑满还没收敛，说明初始结构里有严重冲突，回到solvate和genion重做，而不是加步数死磕。

收敛不等于能量越低越好。EM的目的从来不是找全局能量最小点，而是找一个合理的起始构型。原文有句话我印象很深：“我们不可能知道体系是否处于全局能量最小点，这也不是能量最小化的目的。”每次跑到这一步，我会用gmx energy看一眼EM阶段的总能量趋势，确认它在持续下降且梯度收敛，再进预平衡。

3. NVT与NPT预平衡：为什么温度要先于压力稳定

3.1 NVT预平衡：让温度先把溶剂焐热

EM之后体系几何上合理了，但溶剂和离子的速度分布还是人为给定的，溶质取向也没有在目标温度下弛豫过。如果跳过预平衡直接跑无限制MD，最初的几皮秒内，速度产生与坐标缩放同时作用在一个远离平衡的体系上，温度会剧烈震荡，局部速度过大时LJ势直接爆炸。所以要做NVT预平衡，也就是在粒子数N、体积V和平均温度T都保持恒定的系综下，让溶剂和离子在蛋白周围把取向调整好。

原文里有个判断我很认同：在加恒压器之前先做NVT，通常比直接上NPT更稳健。恒压器会对坐标做整体缩放，速度产生和坐标缩放叠加，会放大一个尚未平衡体系的不稳定因素。先稳定温度，再让体积去适应压力，分步处理比一步到位可靠得多。

NVT的mdp关键参数如下：

integrator = md nsteps = 50000 dt = 0.002 tcoupl = v-rescale tc-grps = Protein Non-Protein tau_t = 0.1 0.1 ref_t = 300 300 pcoupl = no constraints = h-bonds

运行命令：

gmx grompp -f nvt.mdp -c em.gro -r em.gro -p topol.top -o nvt.tpr gmx mdrun -deffnm nvt -v

tcoupl=v-rescale是比Berendsen更温和的恒温方式，它会在目标温度附近加一个随机项，避免Berendsen那种温度“死寂”的现象；tc-grps把蛋白和非蛋白分开耦合，让溶剂弛豫速度不被蛋白拖慢；tau_t=0.1表示耦合时间常数0.1ps，值越小温度弛豫越快，但过小会让温度波动不自然；pcoupl=no明确这个阶段不控压。-r em.gro的作用是写入参考结构，给预平衡阶段的位置限制提供坐标基准，把蛋白重原子约束在初始位置附近，允许柔性侧链和溶剂自由运动。

3.2 NPT预平衡：压力和密度一起稳定

NVT阶段的目标是让温度稳定。但大多数实验室条件是恒压而非恒容，所以真正的成品模拟通常在NPT系综（等温等压系综）下进行，盒子尺寸随压力波动自动调整，直到密度达到目标值。先做NVT再做NPT，相当于先把温度这个自由度稳定下来，再让体积去适应压力。

NPT阶段需要在NVT基础上打开恒压器：

integrator = md nsteps = 50000 dt = 0.002 tcoupl = v-rescale tc-grps = Protein Non-Protein tau_t = 0.1 0.1 ref_t = 300 300 pcoupl = Parrinello-Rahman pcoupltype = isotropic tau_p = 2.0 ref_p = 1.0 compressibility = 4.5e-5 constraints = h-bonds

运行命令：

gmx grompp -f npt.mdp -c nvt.gro -r nvt.gro -p topol.top -o npt.tpr gmx mdrun -deffnm npt -v

pcoupl=Parrinello-Rahman是GROMACS在高精度模拟里常用的恒压器，tau_p=2.0是压力耦合时间常数，ref_p=1.0对应1个标准大气压，compressibility=4.5e-5是水的等温压缩系数，这套数值用于水溶液体系很稳。如果跑的是膜蛋白或其他方向性体系，pcoupltype要改成semiisotropic或anisotropic，各向同性水溶液用isotropic即可。

3.3 grompp与mdrun：三步命令与检查点

把EM、NVT、NPT三个阶段串起来看，用的是同一套命令模板，改的就是integrator、tcoupl、pcoupl三个族参数。我习惯在每个阶段之间加一个检查点：EM看em.log梯度是否收敛，NVT看温度是否稳定在300K附近，NPT看密度是否收敛到合理区间。

gmx energy -f nvt.edr -o temperature.xvg gmx energy -f npt.edr -o density.xvg

在交互界面选Temperature或Density对应的序号，按回车就能导出时间序列曲线。判断标准：NVT末段温度均值与300K偏差不超过3K，NPT末段密度在1000 kg/m^3附近波动，如果波动持续很大，说明盒子尺寸还在调整，把nsteps加大到100000再平衡。

3.4 预平衡阶段最容易翻车的三个参数

第一个是dt。含氢键约束的体系常规用2fs（0.002ps），但如果体系里有高频振动的小分子或特殊金属中心，建议降到1fs，否则能量会在几步内爆炸。

第二个是位置限制。NVT和NPT阶段不加位置限制的话，蛋白整体会在预平衡期间发生平移或旋转，导致生产模拟从错误的局部几何开始。常规做法是在mdp里加define=-DPOSRES让posre.itp生效，生产模拟阶段再把这个宏去掉。

第三个是温度耦合组。tc-grps只写一个System会让蛋白和溶剂共享一个温度耦合，蛋白内部温度弛豫被抹平，应该分开Protein和非Protein两组。原文流程里提到的预平衡顺序是等压模拟在前、等温模拟在后，我在实操中按NVT→NPT执行，对应等温→等压，顺序不能反：让未平衡温度体系过早接触恒压器，坐标缩放会把温度模型连带破坏。

提示：NVT和NPT的mdp可以共用一份模板，只改恒压器相关字段，不要每次从零重写。

4. 数据分析避坑：RMSD、RMSF、Rg、SASA的判读红线

4.1 先做PBC重包装：轨迹里的“破碎”分子是假象

拿到生产模拟轨迹后，第一件事不是画RMSD，而是处理周期性边界条件。由于PBC，轨迹帧里的蛋白可能看起来“破碎”——一部分原子在盒子一侧，另一部分在另一侧，甚至整体像跑出了盒子。这完全来自坐标记录与可视化方式，分子实际相互作用始终按拓扑结构计算。mdrun输出“破碎”轨迹是为了计算效率，分子没有真的跑出去，因为对无限体系而言没有所谓的外部。

分析前必须先做重包装：

gmx trjconv -s npt.tpr -f md.xtc -o md_pbc1.xtc -pbc whole gmx trjconv -s npt.tpr -f md_pbc1.xtc -o md_pbc2.xtc -pbc mol

-pbc whole先把多聚体按完整分子编号处理，让蛋白分成不碎；-pbc mol再把每个分子的中心重新映射回盒子里，恢复完整构象。如果后续要算RMSD，建议再做一步-fit rot+trans，把所有帧叠合到参考结构上，否则蛋白整体的平动和转动分量会污染RMSD数值。很多人拿到轨迹直接算RMSD，曲线飘得离谱，缺的就是这步处理。

4.2 RMSD：能看平衡，但别迷信平台值

RMSD是将结构与参考结构做最小二乘叠合后逐帧计算的均方根偏差，GROMACS里对应gmx rms。参考结构一般选EM或NVT阶段的结构，也可以是实验晶体结构。原文示例里，蛋白和配体的RMSD在初期短暂扰动后，大约在8ns、24ns处出现波动，之后趋于约0.7nm的基线，上下误差不超过0.2nm。这里有两个常见误解必须破除：

第一，RMSD是简化度量，一个单一数值可以有多种数学解，完全不同的构象偏移可能算出同一个RMSD。所以RMSD本身不能区分两种构象，只能说明整体偏离参考量有多大。

第二，不能仅凭RMSD平台就断言模拟已经收敛。结构仍然可以在RMSD平稳的同时发生系统性构象漂移。原文写得很直接：“不能仅仅通过查看RMSD就认为模拟是‘稳定’的或‘收敛’的而终止模拟。”所以必须结合Rg、SASA、氢键数量等多个指标交叉验证。

4.3 RMSF：用0.3nm截止值锁定高柔性区域

RMSF衡量每个残基在模拟期间相对平均位置的波动大小，反映局部位点柔性，GROMACS里用gmx rmsf计算。原文示例以0.3nm波动为截止值，A链的2-11、36-94、105-107、244-257、482-508等区域，B链的2-14、26-35、44-101、245-263、464-530这些区段被标记为高波动。口袋附近的氨基酸分布，可以结合配体结合位点周围残基来看，如果高波动区域正好落在结合口袋附近，说明结合界面的柔性在参与调控配体停留。

我一般会把每个残基的RMSF数值导出来，和蛋白二级结构一起看。loop区域波动大是正常的，但α-helix中段出现大面积高波动，就要回查体系是否真的平衡了。

4.4 Rg与SASA：构象转变和密实度的交叉验证

回旋半径Rg衡量的是蛋白密实度。一个折叠良好的蛋白，Rg应基本稳定；如果去折叠或构象大幅转变，Rg会随明显变化。原文例子里7ns处有一个比较明显的扭转，24ns还有一次波动，正好对应RMSD的两个波折，这种RMSD与Rg同步波动的模式，比单独看RMSD可靠得多。

SASA是溶剂可及表面积。原文示例中0-50ns里蛋白质的SASA波动下降，表明蛋白逐步紧缩，结构趋于紧凑。SASA下降通常是结构紧凑化或疏水塌缩的信号，和Rg下降同时出现时，基本可以确认蛋白在向更紧密的构象转移。如果Rg和SASA一个降一个升，就要小心，那可能是构象变化方向不一致，比如蛋白表面扩张的同时核心收缩，这时需要DSSP确认二级结构变化。

4.5 五条高频踩坑记录

记录一：RMSD一平台就判定收敛 现象：跑到20ns看到RMSD曲线平台化，直接写“体系已达平衡”。 原因：RMSD是简化度量，单个值对应多种构象解，平台只代表整体偏离量稳定，不代表构象固定。 解决：结合Rg、SASA、氢键数量交叉验证，再在VMD里看几帧轨迹，确认没有系统性漂移。

记录二：xvg文件用Excel打不开 现象：双击rmsd.xvg弹出乱码或直接被当成未知格式。 原因：xvg是纯文本格式，带注释和数据类型说明，Excel不认。 解决：用文本编辑器打开看前几行，或转成csv再进Excel。原文配套的常规动力学分析汇总excel已经帮你把rmsd这类数据整理好了，省去手动转换。

记录三：轨迹里蛋白“跑出盒子” 现象：VMD里蛋白一半在盒子内一半在盒子外，甚至像“分子碎了”。 原因：周期性边界下分子跨边界是正常现象，坐标记录方式造成破碎感。 解决：分析前必须执行trjconv -pbc whole与-pbc mol，让分子恢复完整形态再做后续计算。

记录四：RMSD和MMPBSA残基分解结论对不上 现象：稳定构象的氢键分析显示某残基与配体成氢键，但MMPBSA残基能量分解显示该残基贡献为正。 原因：稳态分析只看平衡后的某一帧或短时间窗口，残基能量分解是对全模拟轨迹的平均统计，时间尺度不同。 解决：论文里分开表述——稳态分析描述“稳定构象中的相互作用”，能量分解描述“整个模拟过程中的平均贡献”，两者互为补充。

记录五：最小化没收敛就硬跑生产模拟 现象：NVT刚跑几步体系温度飙升，能量出现NaN，LJ势爆炸。 原因：EM阶段没有真正收敛，初始结构里仍存在原子冲突。 解决：回查em.log最终最大力是否低于emtol，没收敛就重新solvate/genion换随机种子，而不是加步数硬扛。

5. MMPBSA结合自由能：从TOTAL到残基分解的解读顺序

5.1 MMPBSA的原理：一份轨迹，多帧采样

MMPBSA是估算蛋白-配体结合自由能最常用的方法之一。它不直接算绝对自由能，而是从MD轨迹里提取大量构象，对每帧构象分别计算气相能（分子力学项）和溶剂化能（泊松-玻尔兹曼/广义伯恩表面积项），再在bound和unbound状态间做差。动力学模拟时间越长，构象采样越充分，能量项数值越稳定。我一般用gmx_MMPBSA工具跑，输入只需要tpr、轨迹索引和配体残基组定义。

5.2 能量项拆解：GGAS为负、GSOLV为正说明什么

拿到MMPBSA输出后，第一页先看TOTAL那一项。原文示例给出TOTAL平均值-19.28 kcal/mol，为负说明蛋白与配体存在稳定结合的倾向。但只看TOTAL不够，要拆开看各项贡献：

这张表的价值在于看清“对冲账”：GGAS为负说明相互作用本身有利，GSOLV为正说明溶剂化在惩罚结合，两项互相抵消，最终TOTAL的符号和量级决定结合能否成立。写论文时不要只报TOTAL，把GGAS和GSOLV的拆解也放进去，否则审稿人会问“自由能靠什么驱动”。

5.3 残基分解：为什么与稳态氢键分析对不上

残基能量分解把总自由能按残基归属，能定位哪些残基在能量上主要贡献结合。原文示例里，A链的ALA-1和ARG-21自由能贡献为正，属于干扰型；A-PRO-23、A-MET-233、A-PRO-234、B-LYS-329、B-SER-350、B-GLU-353、B-GLN-354为负，属于有利贡献。

这里我踩过坑：把残基能量分解结果和稳态截帧的氢键分析对比，发现结论对不上，一度以为计算错了。后来想明白——稳态分析是从平衡模拟的某一帧去找相互作用，看到的是“当前稳定构象里的接触关系”；残基分解是从整条轨迹上千帧构象里做平均，看到的是“整个模拟过程中这个残基对结合的净贡献”。一个侧重结果，一个侧重过程。比如某个残基在稳态构象里离配体很近形成氢键，但它在模拟前几十纳秒反复远离配体，平均下来净贡献可能反而是正的。所以论文里这两类分析要分开表述，各说各的时间尺度，不要强行让它们互相支持。

5.4 实操建议：取帧、选轨迹、看TOTAL

跑MMPBSA前，我会把轨迹做两步预处理：先用trjconv去掉水分子和离子，只保留蛋白-配体复合物拓扑，再按固定间隔抽帧，避免高度相关的相邻帧重复采样。常见做法是取平衡后段（比如最后20ns）的轨迹，每隔100-200ps取一帧，样本量到上百个，平均值和标准差才稳。如果TOTAL的标准差超过平均值，说明采样不充分或体系在构象盆地之间跳动，这时要先回头确认RMSD是否真的稳定，而不是急着下结论。

能量分解结果和稳态分析对不上很正常，这在方法论上是两回事，不是错误。我后来写论文时会把它们分开表述：“能量分解显示残基X在模拟全程对结合贡献显著，而稳态构象中观察到的主要相互作用残基为Y。”这样既保留了两类信息，也避免了自相矛盾的质疑。

6. 把MD写进论文：方法学段落与结果图的可抄模板

6.1 方法部分：一段话写满全部模拟条件

原文给了一个可以直接改到方法部分的模板，我按自己的表达习惯调整后是这样：

将对接结果选为初始结构，蛋白力场选择AMBER14SB，小分子配体在量化软件中于B3LYP/6-31G*基组条件下进行量子化学优化，涉及键长、键角、二面角的修正以及RESP 2.0固定电荷的计算，后选择GAFF2力场描述配体。采用TIP3P水模型对蛋白-配体复合物体系添加溶剂，建立水盒子并添加钠离子平衡体系电荷。短程相互作用采用Verlet截断方案，静电相互作用采用PME方法，使用最陡下降法进行最大步数50000步的能量最小化，库仑力截距与范德华半径截距均设为1.4nm。随后在NVT正则系综与NPT等温等压系综下分别预平衡体系，最后在常温常压下开展50ns的无限制分子动力学模拟。

写方法部分时最关键的是九个要素不能缺：力场、水模型、盒子、离子、静电处理、截距、最小化方法、系综顺序、模拟时长。缺一个，之后投稿就可能要补材料。

6.2 结果部分：从RMSD到HBNUM的标准表述

结果部分按“宏观稳定性→局域柔性→结构紧密度→表面极性→氢键→二级结构”的顺序组织。原文给了很具体的表述模板：

RMSD：体系中蛋白-配体复合物的RMSD值在经历7ns与24ns两次波动后趋于约0.7nm，上下偏差不超过0.2nm，表明体系在模拟后期达到稳定的构象区间。 RMSF：以0.3nm为截止值，A链2-11、36-94等区段与B链对应柔性区域波动显著，覆盖配体结合口袋附近残基，说明结合界面存在构象可塑性。 Rg：回转半径在7ns与24ns处各有一次波动，正好对应RMSD中的波折，表明体系在这两个时间点经历从多个不稳定状态的转变。 SASA：0-50ns内SASA波动下降，表明蛋白逐步紧缩，结构趋于紧凑。 HBNUM：模拟期间氢键数量基本平稳，在1-2个范围内波动，无重大变化。 DSSP：初始帧与末尾帧对比显示，coil卷曲减少，B-sheet与Bend增多，A-helix增多，二级结构总体发生变化。

这段文字我在跑完体系后直接复制过来，只改数值和蛋白名字，能省一晚上写结果的时间。

6.3 图表规格与引用习惯

图表方面，黑白图按600dpi提供，彩色图按300dpi提供，这是期刊投稿的通用要求。引用文献我按原文涉及的工具整理一份习惯用法：GROMACS引用官方推荐的方法学文献；PDB数据库引其数据库专辑；Pymol是开源软件，可不强制引用；量化软件Orca需要引方法学论文；HDock和gmx_MMPBSA分别引各自的服务器/方法学文献；POCASA和proteinplus这类口袋预测工具在体系准备部分给一条引文更稳妥。

从那以后我每次拿到新体系，都强制按这套顺序走一遍：EM看能量梯度、NVT看温度、NPT看密度，任何一环不稳就回头查参数，而不是急着跑50ns的成品模拟。这份PDF最实用的地方是把模拟流程和论文写法放到了一起，跑完一个体系，方法部分和结果部分也就有了底稿，少翻车是真的。希望帮到你。

本文还有配套的精品资源，点击获取