1. GROMACS分子动力学模拟工程导论
分子动力学（Molecular Dynamics, MD）模拟的本质，是在给定力场与边界条件下，对由成千上万乃至上亿个原子组成的体系求解牛顿运动方程，从而获得体系随时间演化的轨迹。GROMACS 正是这一领域应用最广、性能最成熟的开源引擎之一。它最初为蛋白质与水溶液的模拟而设计，经过二十余年发展，已能覆盖生物大分子、膜体系、聚合物、无机材料、离子液体乃至粗粒化模型。理解 GROMACS 的定位，关键不在于记住命令，而在于认清它在整个“分子模拟工程”链条中扮演的角色：它是计算内核，而不是建模工具、不是可视化工具、也不是数据分析工具。
从工程视角看，一次完整的模拟任务可以拆成四个相互衔接的环节：体系构建、能量与平衡、生产模拟、轨迹分析。GROMACS 主要承担后三个环节，尤其是“平衡 + 生产”这一计算密集型部分。体系构建通常由 pdb2gmx、editconf、genbox、genion、solvate 等 GROMACS 自带工具完成，但拓扑生成、力场选择、质子化状态判断往往需要外部工具（如 pdb2pqr、CHARMM-GUI、Packmol）配合。因此，一个合格的 GROMACS 使用者，必须同时具备“力场与拓扑”的知识和“命令行与文件格式”的工程能力。
.gro（坐标）、.top（拓扑）、.mdp（运行参数）三类文件，输出是 .trr/.xtc（轨迹）、.edr（能量）、.log（日志）。任何模拟问题，最终都可以归结为这三类输入是否自洽、输出是否合理。
GROMACS 的核心能力可以归纳为五点。第一，高效的非键计算：通过邻居列表（neighbor list）、Verlet 截断、PME（Particle Mesh Ewald）静电求和，把 O(N²) 的相互作用降到接近 O(N log N)。第二，丰富的积分与控温控压算法：leap-frog 积分、velocity-rescale 与 Nosé-Hoover 控温、Berendsen 与 Parrinello-Rahman 控压，可适配从平衡到生产的多种需求。第三，多平台并行：支持 OpenMP、MPI、CUDA、HIP，可在单机多卡、CPU 集群、GPU 集群上运行。第四，约束与自由能方法：LINCS/SHAKE 约束、umbrella sampling、自由能微扰（FEP）、热力学积分（TI），支撑结合自由能计算。第五，可复现的工程化流程：所有参数写入 .mdp，所有拓扑写入 .top，配合 .tpr 二进制运行文件，保证同一输入在不同机器上得到一致结果。
典型工作流可以概括为“建模—溶剂化—加离子—能量最小化—NVT/NPT 平衡—生产模拟—分析”。以溶菌酶在水中的模拟为例：先用 pdb2gmx 选择力场（如 AMBER99SB-ILDN）并生成拓扑；用 editconf 设置 1.0 nm 的立方盒子；用 solvate 填充 SPC/E 水；用 genion 加入 0.15 mol/L NaCl 中和电荷；随后依次运行能量最小化（steep 最速下降，最多 50000 步）、NVT 平衡（100 ps，V-rescale，300 K）、NPT 平衡（100 ps，Parrinello-Rahman，1 bar）、生产模拟（100 ns，2 fs 步长，LINCS 约束氢键）。每一步的 .mdp 参数都直接决定后续轨迹是否物理合理。
学习路径上，建议按“文件格式 → 力场与拓扑 → 平衡策略 → 生产参数 → 分析工具”的顺序推进。初学者最容易犯的错误，是跳过能量最小化直接跑生产模拟，导致体系因初始构象冲突而崩溃；或者把 NVT 和 NPT 的顺序颠倒，造成盒子尺寸与温度耦合异常。另一个常见误区是忽视 .mdp 中 tcoupl、pcoupl、constraints 三者的相互影响：例如使用 2 fs 步长时，必须对含氢键的体系启用 LINCS 约束，否则高频振动会导致积分发散。
从工程能力角度看，GROMACS 使用者需要掌握三类技能。第一类是文件操作能力：能读懂 .gro 的固定格式（标题行、原子数、坐标行、盒子向量），能检查 .top 中 [ atoms ]、[ bonds ]、[ angles ] 的原子编号是否与坐标一致。第二类是参数调试能力：能根据体系大小和硬件条件，合理设置 cutoff-scheme、coulombtype、nstlist、nsteps。第三类是结果判读能力：能从 .log 中读取平均温度、压力、密度、势能，从 .edr 中提取时间序列，判断模拟是否达到平衡。