1. 为什么这四种文件能决定你模拟的成败做分子动力学模拟的人几乎每天都和GROMACS打交道。但我发现很多刚接触GROMACS的朋友甚至一些跑了一两年模拟的同行对.tpr、.xtc、.edr、.cpt这四种文件的理解还停留在“知道名字、会用命令”的层面。等到真正分析数据或者中途续跑的时候问题就一个接一个冒出来为什么跑出来的轨迹坐标和能量对不上为什么续跑之后之前的进度全丢了为什么别人给的.tpr文件我这边跑不了这篇文章聚焦一个核心问题GROMACS模拟生命周期里这四类文件分别承担什么角色怎么生成、怎么使用、怎么避坑。我会按照实际模拟流程——从准备输入、运行模拟、到最后分析——把每条命令背后的逻辑和每个文件的正确打开方式讲清楚。适合刚入门MD模拟的学生也适合想系统梳理文件体系、排查历史疑难问题的研究者。先说一个我自己的惨痛经历。有一次跑一个蛋白质配体复合物的400 ns模拟跑到第280 ns的时候集群节点故障任务被杀了。我当时没在意以为有.cpt检查点文件就能无缝续跑。结果恢复的时候发现.cpt文件是损坏的不得不从200 ns的备份重新开始白白浪费了三天机时。从那以后我花了大量时间研究这四类文件的生成机制、写入频率和备份策略踩过的坑远比文档里写的要多。希望你读完这篇文章之后能比我当年少走一些弯路。2. 四种文件角色拆解谁是输入、谁是输出、谁管恢复2.1 .tpr文件整个模拟的“施工图纸”.tpr文件全称是Portable Run Input File中文可以理解为“便携式运行输入文件”。它是在运行mdrun之前由grompp命令把三样东西打包在一起生成的分子的拓扑结构也就是原子之间的连接关系、化学键、键角、二面角等模拟的坐标信息所有原子的三维坐标全部的运行参数温度、压力、步长、输出频率、约束算法等你可以把.tpr想象成一张完整的建筑图纸加施工方案。mdrun读入.tpr之后就按照这张图纸一步步计算每个原子之间的相互作用力、更新坐标和速度。在实际操作中.tpr文件有四个特别容易忽视的点第一.tpr文件应当在模拟开始时就用grompp生成而不是在模拟过程中反复修改。我之前见过有同学为了提高性能手动去编辑.tpr里的参数这非常危险。.tpr是二进制文件任何不当修改都会导致计算结果不可复现甚至直接报错。正确的做法是修改.mdp参数文件然后重新生成.tpr再重新跑模拟。如果只是想延续之前的模拟应该用的续跑机制是基于.cpt的后面细说而不是改.tpr。第二.tpr文件里的坐标是模拟起始的坐标不是结束时的坐标。很多人在分析结束构象时习惯性去.tpr里找最后的坐标这是完全错误的。mdrun运行过程中坐标更新后是写进.xtc或.trr轨迹文件里的.tpr只在启动时被读取一次。如果你需要从.tpr中提取起始构象用gmx convert-tpr如果需要提取模拟结束时的构象需要从.xtc轨迹里提取最后一帧。第三.tpr文件和GROMACS版本强绑定。用GROMACS 2021生成的.tpr文件GROMACS 2023基本可以读有向后兼容机制但反过来不行——旧版本读不了新版本生成的.tpr。更麻烦的是如果你用的集群上安装的GROMACS版本和你本地不完全一致生成的.tpr可能无法运行。所以在提交任务前建议用gmx check命令检查一下.tpr文件的版本和完整度。第四.tpr文件较小通常几百KB到几MB但它和.cpt文件配合使用才能实现真正的断点续跑。很多初学朋友误以为.tpr是唯一的运行输入续跑时只要再次指定.tpr就行。实际上如果没有.cpt检查点文件mdrun只能从初始坐标重新开始算之前算的全部作废。2.2 .xtc文件最常用的压缩轨迹.xtc是GROMACS默认输出的压缩轨迹文件只保存原子坐标信息不保存速度、力、能量等数据。它的核心卖点是“压缩率高、体积小”这也是为什么它成为构象分析的默认选择——一个包含10万原子的体系跑100 ns模拟每10 ps存一帧.xtc可能只有几百MB而如果全部存成未压缩的.trr轨迹文件体积会膨胀到几十GB甚至上百GB。但正因为.xtc是有损压缩它有一个实际使用中需要了解的细节坐标精度默认保留到小数点后3位单位是nm。对于大多数分析——比如计算RMSD、RMSF、氢键、自由能曲面——这个精度完全足够。但如果你要做高精度的能量分解或者极细微的构象差异分析建议同时输出.trr文件无损轨迹用于分析。我在跑配体结合自由能计算比如MM/PBSA时通常会让.xtc和.trr同时输出因为MM/PBSA需要速度和力信息.xtc提供不了。.xtc文件在使用过程中有几个高频问题轨迹时间轴和.edr对不上.xtc的写入步数和.edr的能量写入步数是独立控制的分别由.mdp中的nstxout-compressed和nstenergy决定。如果你把nstxout-compressed设为5000步nstenergy设为1000步那么它们的帧数、时间戳就完全不同。分析时如果直接把.xtc的第N帧对应.edr的第N行数据结果当然会错位。正确做法是各自按时间ns对齐或者统一设置输出频率。轨迹断帧模拟被中断后.xtc文件的最后一帧可能不完整用gmx check检测时会提示“last frame incomplete”。遇到这种情况要么用gmx trjconv截掉最后损坏的帧要么用.cpt续跑后让GROMACS自动补全。时间单位默认是ps.xtc内部记录的时间单位是皮秒ps很多第三方分析工具期望的是纳秒ns不换算直接分析会让你得到离谱的结果。我自己就吃过这个亏分析RMSF曲线发现横坐标到几百一开始还以为是正常现象后来才反应过来单位没换算。2.3 .edr文件能量数据的唯一来源.edr文件是GROMACS的能量文件以一个紧凑的二进制格式保存了模拟过程中所有热力学量包括势能、动能、总能量、温度、压力、密度、体积、盒长、各项相互作用能库仑、LJ等、氢键数、体系总偶极矩等。所有和“能量”相关的分析——能量漂移检查、体系平衡判断、自由能计算、温度压力波动分析——都要用到它。.edr最常见的坑有三个第一个坑.edr文件是二进制的不能直接less或cat查看。必须用gmx energy读取。这个命令会交互式地询问你想输出哪些能量项。你可以输入1 2 3 4这样的索引号也可以输入Temperature、Pressure之类的名字最后会生成一个.xvg文件再到Grace或Python里画图。我第一次用的时候不知道这个交互机制卡在那里不知所措。实际上gmx energy是支持管道输入的比如echo 1 2 3 4 0 | gmx energy -s topol.tpr -f prod.edr -o energy.xvg这样可以避免人工交互。第二个坑.edr文件的频率和轨迹文件不同步。前面已经提到nstenergy和nstxout-compressed是两个完全独立的参数。如果你在.mdp里设置nstenergy1000、nstxout-compressed5000那么能量每2 ps一个点而坐标每10 ps一帧。分析的时候如果想做“某一帧对应的能量”需要按时间查找不能按行号索引。第三个坑.edr文件的体积会持续增长。对于长时间的模拟比如微秒级.edr文件可以达到好几个GB。磁盘空间不够时mdrun不会自动清理而是直接报错退出。所以对于长时间模拟我建议设置合理的nstenergy——不是越小越好也不是越大越好。默认值是1000步相当于2 ps对于大部分体系足够精细。如果磁盘吃紧调到5000步10 ps也不会损失太多能量信息。2.4 .cpt文件断点续跑的救命稻草.cpt检查点文件可能是四者中利用率最低、却最关键的。它保存了模拟的完整状态——包括所有原子的坐标、速度、力的信息以及随机数生成器的状态、温度耦合和压力耦合的内部状态。一句话.cpt是把模拟“暂停”的完整快照。使用场景非常明确模拟中断后基于.cpt续跑而不是从头跑。命令也很简单mdrun -s topol.tpr -cpi state.cpt -deffnm prod-cpi参数告诉GROMACS“用这个检查点文件继续跑”。GROMACS会从state.cpt记录的精确状态接着算而不是从.tpr的初始坐标重来。这就意味着即使你的模拟已经在集群上跑了200 ns突然被系统杀掉只要最近的.cpt文件写入成功你也能从那个时间点继续而不是损失全部进度。但是.cpt文件有一个非常致命的特点它是定期写入的不是实时写入的。默认情况下mdrun每15分钟写一次检查点。如果你的模拟在14分59秒时崩溃那么你最多丢失将近15分钟的模拟进度。对于长模拟来说15分钟不算多但对于高精度短模拟比如自由能微扰15分钟的损失可能就要重跑很多步。如何调整检查点写入频率在.mdp文件里加一行nstcheckpoint 5000这里的5000是步数不是秒数。以2 fs步长算nstcheckpoint 5000相当于每10 ps写一次检查点。如果你的体系性能很高比如每天能跑50 ns那么10 ps一个检查点意味着每天写5000次每次检查点文件可能几十MB这个磁盘吞吐量可能变成瓶颈。所以我一般建议常规模拟nstcheckpoint 50000即100 ps折中安全和性能昂贵模拟比如增强采样、自由能计算nstcheckpoint 10000即20 ps磁盘非常充裕的集群可以更频繁但实测下来没必要另外一个容易忽略的点-cpi可以有多个检查点吗答案是肯定的。如果你在.mdp里设置了nstcheckpoint那么mdrun会在模拟目录下生成多个检查点文件默认命名如state.cpt、state1.cpt、state2.cpt...。使用-cpi state.cpt会从最近一次成功写入的检查点继续。但如果你希望回退到某个更早的状态比如怀疑某段模拟有问题可以指定-cpi state100.cpt。这个特性在排查模拟发散问题时特别有用——你可以把模拟切回到发散前的某个时间点修改参数后重新跑。3. 从生成到分析一条完整的时间线3.1 模拟准备阶段用grompp生成.tpr并检查模拟的第一步永远是准备拓扑和坐标然后生成.tpr。这一步最常用的命令gmx grompp -f prod.mdp -c npt.gro -t npt.cpt -p topol.top -n index.ndx -o prod.tpr -maxwarn 2这个命令里-f指定.mdp参数文件-c指定输入的坐标文件通常是前一步平衡完的.gro文件-t指定可选的检查点文件当你需要从某个检查点继续准备时用-p指定拓扑文件-n指定指数文件对于定义组非常关键-o是输出的.tpr文件。一个容易忽略的细节-maxwarn参数。GROMACS的grompp在遇到警告时默认直接拒绝生成.tpr这是为了保证物理合理性。但有些警告是无害的——比如“你的体系中含有孤对电子”或者“某些二面角参数缺失但用通用参数补齐了”。如果你确定警告不影响结果可以用-maxwarn来提高容忍度。但是永远不要无脑设-maxwarn 1000。每个警告都应该被逐条审查。曾经有个朋友因为忽略了“温度耦合参考温度与实际不符”的警告导致整个模拟温度一直飘最后数据全部作废。3.2 模拟运行阶段mdrun的参数选择和文件产出运行模拟的核心命令是mdrun。最常规的并行运行方式gmx mdrun -deffnm prod -v -ntomp 8 -nb gpu -pme gpu -bonded gpu -update gpu-deffnm prod意味着所有输出文件都以prod为前缀。于是你会看到以下文件依次生成prod.xtc压缩轨迹prod.edr能量文件prod.log日志文件prod.cpt检查点文件prod.trr如果.mdp中设置了nstxout、nstvout或nstfout会生成无压缩轨迹默认不生成-cpi续跑的场景gmx mdrun -s prod.tpr -cpi prod.cpt -deffnm prod -ntomp 8 -nb gpu -pme gpu这里有个重要提醒续跑时-s指定的.tpr必须和初始运行时完全一致否则GROMACS会报错。这是版本校验的一部分防止你在中途偷偷改了参数然后假装继续跑。mdrun运行过程中除了.cpt还有一个文件值得注意#prod.cpt.1#这类带数字后缀的临时文件。这些是mdrun在写入新检查点时产生的原子写文件目的是防止意外断电导致检查点文件损坏。你不能手动删除它们因为mdrun会在合适的时机自动清理。如果在目录里看到一堆临时文件说明上次运行没有正常退出或者磁盘几乎满了。3.3 数据后处理阶段xtc与edr的正确打开方式模拟完成后分析阶段涉及最多的两个文件就是.xtc和.edr。对.xtc做的最多的操作是去除周期性边界效应、对齐和截取。这些操作都用gmx trjconv完成gmx trjconv -s prod.tpr -f prod.xtc -o prod_noPBC.xtc -pbc mol -center-pbc mol表示以分子为单位处理周期性边界条件——对于分析跨越盒子边界的分子构象至关重要。如果不做这一步你会看到某些原子“跳”到盒子对面RMSD和氢键分析结果完全错误。gmx trjconv还有很多使用细节比如-dump可以提取某一帧-fit可以对轨迹做最小二乘拟合作叠合-skip可以间隔输出。这些都是日常高频操作。处理.edr使用gmx energy但实际工作中我很少手动交互。推荐一个更高效率的做法——直接用Python库比如MDAnalysis或mdanalysis的.edr解析读取或者用GROMACS自带的gmx energy配合管道输入批量提取。举个例子echo 10 11 12 0 | gmx energy -f prod.edr -o temperature_pressure_density.xvg这样一次性提取温度第10项、压力第11项、密度第12项生成一个.xvg文件再用Python绘制平衡曲线。3.4 数据转换与扩展工具链除了GROMACS自带的工具我自己在实际工作中还会用到两款生态工具在这里一并分享MDAnalysisPython库可以无缝读取.tpr、.xtc和.edr。它的MDAnalysis.coordinates.XTC模块可以直接读取.xtc配合MDAnalysis.analysis.rms、rmsf等模块可以写出比gmx rms更灵活的分析脚本。特别是当你需要对轨迹做自定义分析比如计算某个二面角的分布、分析特定残基的接触图MDAnalysis会非常顺手。VMD虽然VMD不能直接读.edr但它读.tpr和.xtc非常流畅适合快速查看轨迹和渲染分子图。VMD里有一个gromacs插件可以直接加载.tpr作为拓扑.xtc作为轨迹然后做可视化排查——比如检查模拟中是否出现原子重叠、配体是否脱离结合口袋等情况。每次模拟完成后我会用VMD快速“扫”一遍轨迹的关键帧确认没有问题后再开始批量分析。4. 常见报错与坑位排查实录这些年帮人排查模拟问题时我总结出了几类和这四种文件相关的典型报错。下面按现象整理成一个速查表后面再逐个展开。报错/现象可能原因解决方案Cannot find a usable .tpr file路径输错或.tpr未生成检查文件路径重新gromppFatal error: step ... too long模拟发散体系爆炸检查.edr能量用.cpt回退减小步长GROMACS does not have a version of tpr.tpr版本过新或过旧统一GROMACS版本或用高版本读低版本Last frame incomplete.xtc最后几帧写入不全用gmx trjconv截断或续跑补全Reading energy file failed.edr文件损坏用.cpt续跑生成新.edr检查磁盘空间Can not read checkpoint file.cpt与.tpr不匹配或损坏检查文件命名用备份.cpt回到最近完整检查点能量文件行数和轨迹帧数不一致nstenergy和nstxout-compressed不同按时间对齐统一输出频率轨迹中原子“飞舞”周期性边界条件未处理trjconv -pbc mol -center续跑后之前的进度丢失未使用-cpi或.cpt过期确认.cpt写入时间用-cpi续跑4.1 模拟发散的诊断从.edr和.log的异常行为找线索模拟发散是最常见的“大事故”之一。表现是mdrun突然报错退出日志文件里出现大量Step ... too long或者WARNING: Listed nonbonded interaction between particles ...。排查思路是这样的先用gmx energy查看.edr里最后几个能量点。如果势能在最后一个输出点急剧上升比如从-100000 kJ/mol跳到1e10说明体系已经发散发生的时间大约在.edr最后一个正常点到后续几个异常点之间。打开.log文件搜索Step和Energy关键字。mdrun会在每步计算能量时记录取决于nstcalcenergy看是否存在突然异常值。用.cpt回退到发散前的检查点修改参数比如减小步长dt、增加约束算法容差后重新跑。我遇到过一个案例体系里有两个氯离子靠得太近导致静电排斥能接近无穷大。当时.edr的温度曲线完全正常但势能在某一步突然暴涨。排查后发现是我在准备体系时不小心把离子放到了蛋白质内部而不是溶剂中。修正坐标后重新生成.tpr问题解决。4.2 .xtc轨迹撕裂问题周期性边界的经典翻车现场分析轨迹时最常见的“诡异”现象就是某几个原子在轨迹中突然从盒子一侧“瞬移”到另一侧整个分子看起来被撕成两半。这不是分子真的穿过了盒子而是因为GROMACS的坐标存储是基于最小镜像约定的分子跨过盒子边界时坐标表示会跳变。解决办法就是你在.mdp里设置pbc xyz然后在分析前用gmx trjconv -pbc mol -center处理轨迹。我甚至建议在grompp生成.tpr时就直接确认.mdp中的pbc设置因为后续很多问题的祸根都在这一步埋下。我在处理计算蛋白质-配体结合自由能时还养成了一个习惯先用gmx trjconv -dump提取模拟中间的一帧到VMD里快速查看配体是否还在结合口袋中。如果配体已经飘走那么这个模拟前期的轨迹基本不能用于结合模式分析最多只能用于计算结合自由能的热力学积分。4.3 续跑失败的常见原因.cpt和.tpr的版本、时空一致性续跑失败的报错信息五花八门但核心原因绝大多数是两类.tpr和.cpt不匹配比如你重新生成了.tpr然后尝试用旧的.cpt续跑。GROMACS会严格检查.tpr中的体系描述原子数、拓扑、盒子是否与.cpt一致。不一致就报错。.cpt文件过期比如模拟中断前最后一次成功写入检查点是在50 ns但.cpt文件本身是500 ns模拟中的第50 ns的你直接用错文件名导致从起点重跑。这里有一个被很多人忽略的命令gmx check -f prod.cpt它可以输出检查点文件的基本信息包括它对应的时间点。续跑之前先跑一下这个命令确认检查点的时间和预期一致能避免很多无效任务。4.4 分析结果异常从单位、频率和时间轴找原因如果你发现RMSD曲线明显偏高、氢键数目特别少、或者自由能曲线的形状怪异大概率不是GROMACS算错了而是你在分析阶段用错了文件或者忽略了单位/频率/时间轴问题。举个例子用gmx rms分析时如果你选择“将轨迹对初始结构进行拟合”这一步需要指定-s作为参考结构。这时-s指定的.tpr文件和-f指定的.xtc文件必须属于同一次模拟。如果你用A体系的.tpr去分析B体系的.xtcGROMACS会直接报错或者输出无意义结果因为它按.tpr里的原子顺序去读.xtc的坐标而原子顺序不一致时RMSD必然巨大。此外.xvg文件里的时间默认是ps。如果你要画“纳秒级的收敛曲线”记得除以1000。很多论文里的坐标轴单位错误就是这样来的——不是模拟算错了是脚本里忘了换算。5. 文件管理习惯长模拟能跑完也别忘了管好备份很多人模拟跑完了任务不报错就觉得万事大吉。等到分析阶段需要某个文件时才发现.edr被误删了、.xtc被截断了、.cpt因为磁盘空间不足根本没写进去。这些本来都是可以避免的。5.1 磁盘空间规划提前算好三种文件的体积增长长时间模拟前先用以下公式估算体积.xtc体积 ≈ 原子数 × 3 × 4字节 × 压缩率通常0.3~0.5× 帧数.edr体积 ≈ 能量项数 × 8字节 × 能量点数相对较小.cpt体积 ≈ 原子数 × (333若干) × 4字节 ≈ 和.trr的一帧差不多大.log体积随模拟时长线性增长每百万步约几百KB用nstxout-compressed 5000步长2fs时每10ps存一帧和nstenergy 5000作为参考配置一个50000原子的体系跑500ns.xtc大概在1~2GB左右.edr在几百MB范围。如果你输出三份轨迹文件.xtc、.trr、.cpt磁盘需求可能翻3~5倍。所以在提交任务前用df -h看一下集群挂载盘的剩余空间。5.2 备份策略定期归档检查点对于超长模拟几百ns以上我强烈建议你定期把关键时间点的.cpt文件复制到另一个磁盘目录或者远程存储。频率可以是每50ns复制一次或者每3天一次。这样一个检查点文件大约几十MB到几百MB占不了多少空间但能让你在模拟后期出现问题时回退到最近的健康状态。我的习惯是mkdir -p checkpoint_backup cp prod.cpt checkpoint_backup/prod_200ns.cpt这样即使模拟中断你损失的也只是上次备份到崩溃之间的那段时间的进度。5.3 记录参数把mdp和版本信息随文件一起归档最后一条小建议每次模拟结束后把.mdp参数文件、.top拓扑文件与.tpr、.xtc、.edr一起归档。同时在README里写清楚GROMACS版本号、使用的力场版本、是否用了GPU加速。很多人隔了半年再回头分析数据时面对一堆没有说明的文件完全无从下手。你现在的五分钟记录能省掉将来五小时的痛苦。6. 从一次真实模拟复盘看文件全流程的坑为了让这篇梳理更有现场感我把一个典型模拟项目的文件流转过程完整复盘一遍把上面提到的所有坑都串起来。假设我们要模拟一个膜蛋白与多肽配体的结合体系计划跑200 ns常规模拟。阶段一准备用pdb2gmx生成拓扑定义二硫键和质子化状态用editconf生成盒子solvate加溶剂genion加离子用grompp生成能量最小化的.tpr用mdrun做能量最小化此时主要产物是.gro和.log进入NVT和NPT平衡。NPT平衡结束时生成一个NPT平衡的.cpt文件这个文件在后续正式模拟中经常要用——因为正式模拟的初始坐标应该从NPT平衡结束的那一帧开始。这里就有一个常见的错误有人用-c npt.groNPT平衡的坐标但忘了带-t npt.cpt。这样虽然坐标对了但速度如果.dtp没写velocity信息和压力耦合状态没有继承NPT平衡的充分演化。在GROMACS中-c只给坐标-t才给完整的运行状态-cpi是续跑时用。三者的区别必须分清楚。阶段二正式模拟grompp -f prod.mdp -c npt.gro -t npt.cpt -p topol.top -o prod.tpr提交mdrun任务-deffnm prod设置nstxout-compressed1000020ps一帧、nstenergy500010ps一个能量点、nstcheckpoint50000100ps一个检查点跑了几天后集群提示维护任务被SIGTERM杀掉此时在目录里应该有prod.tpr prod.xtc 可能是0~120ns的内容 prod.edr prod.log prod.cpt 可能是维护前的最后一个检查点 #prod.cpt.1# 中间临时文件阶段三续跑查看prod.log的尾部确认最后完成的时间步。运行gmx check -f prod.cpt确认检查点对应的时间。用mdrun -s prod.tpr -cpi prod.cpt -deffnm prod续跑。GROMACS会打印“Restarting from checkpoint ... at step 60000000”。这一步有个关键点续跑后的输出文件是直接在原文件上追加还是生成新文件答案是追加。如果你指定相同的前缀-deffnm prodGROMACS会检测到已有同名的.xtc和.edr文件然后将新数据追加到末尾要求拓扑一致。这个机制很方便但也意味着如果你不小心用错.tpr追加的部分会和前面完全不兼容后期分析时各种奇异波形都可能出现。所以续跑前务必检查.tpr的哈希一致性用md5sum记录下最初的.tpr每次续跑前对比一下。阶段四分析gmx trjconv -s prod.tpr -f prod.xtc -o prod_noPBC.xtc -pbc mol -center去掉周期性边界条件gmx rms -s prod.tpr -f prod_noPBC.xtc -o rmsd.xvg计算RMSD收敛曲线gmx energy -f prod.edr -o temperature.xvg检查温度是否稳定在目标温度附近阶段五归档把所有关键文件压缩打包附上READMEtar -czf simulation_archive.tar.gz prod.tpr prod.xtc prod.edr prod.cpt README现在回头看看如果每一步都了解文件的角色和生命周期这些操作其实很流水线但是一旦忽略其中某一步返工的成本就是几天甚至几周的机时。这大概也是为什么圈子里总说GROMACS分析结果的可靠性从你准备.tpr的那一刻就决定了。文件用不对后面的分析再好看也是白搭。我个人的体会是最好的学习方式就是亲自踩坑。但既然这篇文章能帮你避开一部分坑也算没白写。以后跑模拟遇到这四种文件相关的疑难杂症不妨把本文的排查表打印出来贴在显示器旁边。