刘建川・GROMACS MDP 实战 第 7 篇｜输出设置与常见错误修正
本篇是 MDP 参数实战系列收官篇，重点讲解模拟输出参数配置与新手高频报错排查，帮你实现模拟稳定运行、文件规范输出、报错快速处理。
【14 年实战经验】输出参数并非越多越好，合理配置既能满足后续分析需求，又能节省硬盘空间；新手 90% 的报错都源于参数不匹配、运行顺序错误或文件缺失，按本篇方法可快速定位解决。
一、核心输出参数
nstxout = 5000 ; 每5000步保存一次全精度坐标（.trr文件）nstvout = 5000 ; 每5000步保存一次全精度速度（.trr文件）nstxtcout = 5000 ; 每5000步输出一次压缩轨迹（.xtc文件）nstenergy = 5000 ; 每5000步保存一次能量数据（.edr文件）nstlog = 5000 ; 每5000步更新一次日志文件（.log文件）
nstxout：控制将全精度坐标写入输出轨迹文件（通常为 .trr 文件）的频率。例如，nstxout = 5000 表示每隔 5000 步将当前所有原子的坐标写入一次 .trr 文件。nstvout：控制将全精度速度写入 .trr 文件的频率。同样，nstvout = 5000 表示每隔 5000 步输出一次速度信息，时间间隔由步长决定。
nstxtcout：控制将压缩轨迹（通常为 .xtc 文件）写入的频率。.xtc 文件是对坐标进行压缩存储，占用空间小，常用于后续分析。nstxtcout = 5000 即每隔 5000 步输出一次压缩后的坐标。
nstenergy：控制将能量信息（如温度、压力、势能等）写入能量文件（.edr）和日志文件的频率。nstenergy = 5000 表示每隔 5000 步计算并记录一次能量项。
nstlog：控制将模拟的日志信息（如步数、时间、能量等）写入 .log 文件的频率。nstlog = 5000 即每隔 5000 步更新一次日志文件。
✅ 关键说明：
轨迹文件：优先使用压缩轨迹
xtc，文件体积小、读写快，满足绝大多数分析需求；步频设置：常规蛋白 / 小分子体系，直接使用上述参数，无需随意修改；
文件作用：轨迹用于构象分析，能量文件用于判断体系是否平衡。
二、新手 3 类高频报错及解决方案
温度 / 压力急剧飙升，模拟直接崩溃
原因：时间步 dt 设置过大、跳过平衡阶段、控温控压参数错误
解决：固定 dt=0.002 ps，严格按照「能量最小化→NVT→NPT→生产模拟」顺序运行
MDP 参数报错，模拟无法启动
原因：控温分组
tc-grps、弛豫时间tau_t、参考温度ref_t数量不对应解决：分组数量 = tau_t 数值个数 = ref_t 数值个数，保持一一匹配
无轨迹 / 无能量文件输出
原因：输出参数设为 0 或关键字书写错误
解决：直接套用本篇输出参数，检查关键字拼写无误
三、GROMACS MDP 全套参数极简总结
积分器：能量最小化用 steep，动力学模拟用md时间步：无氢约束时固定 dt=0.002 ps温控： V-rescale+Protein SOL+0.1 0.1+300 300控压： Parrinello-Rahman+isotropic+ 标准大气压边界条件： pbc=xyz输出参数：使用本篇压缩输出配置，稳定高效
使用小结
输出优先用压缩轨迹，步频按本篇设置即可，不盲目追求高频输出； 模拟报错先查参数匹配、运行顺序、关键字拼写； 整套 MDP 参数按系列文章固定配置，可满足 90% 以上常规体系模拟。
