分子动力学模拟 (GROMACS)

分子动力学模拟 (GROMACS)：一次分子模拟的完整排练

「GROMACS」在许多计算组的组会上被简化成一个动词——"这个体系你 GROMACS 了多久"。能被当成动词用的软件不多，它算一个：把一团原子的坐标与速度交给它，它按牛顿方程一步一步往前推，推到足够长的时间，你就能从这堆轨迹里读出扩散系数、构象分布、自由能面这些实验上很难直接测量的量。本册要带你去做的，是把"跑一次分子模拟"当成一场完整的排练来学：从选剧本、搭舞台，到走台、正式演出，再到演出后的复盘与差评处理，每一环都亲手过一遍。

为什么用"排练"来讲分子模拟

分子模拟的学习曲线陡，不是因为它难在数学，而是难在环节多、耦合深。力场选错，后面所有参数调得再好都是白跑；能量最小化没做够，升一步温度体系就炸；平衡没判准，成品轨迹拿去分析，结论整段作废。这些坑几乎没有人能靠读理论书躲开——它们只在"你自己跑崩过"之后才真正长进脑子里。

所以本册把知识全部挂到一条主线上：你手里有一个想跑的体系，从零开始把它推到能出报告的那一天。 主线案例刻意避开了教程圈用滥的"蛋白质水盒子"，改用两类在工业界与学术界同样高频、却很少被系统讲解的体系——聚环氧乙烷（PEO）聚合物熔体与 DPPC 磷脂双分子层。蛋白体系只在需要对照的地方出场。为什么这样选？因为聚合物让你必须直面拓扑构建、长链平衡慢、尺度外推这些硬问题；脂膜则把面积压强、相变温度、序参数这些"二维物理"摆上台面。学会这两类，蛋白体系反而是降维。

全册知识地图

学习顺序基本就是排练顺序，但有几条捷径与几处必须绕开的坑，先在图里看一遍全局。

图 0-1 全册知识地图：从剧本到复盘的排练路线

图上有两条容易走错的依赖线值得点出：分析技术（第5章）名义上排在后面，但"先想清楚要分析什么，再去定模拟方案"是老手的工作方式——RDF 要看什么对、氢键判据怎么定，都会反过来改你的输出频率设置；第6章的增强采样则依赖第5章的轨迹预处理，否则 PMF 算出来是碎的。

各章在排练里扮演什么角色

第1章是剧本围读：搞清 MD 到底在算什么，势能函数的骨架长什么样，软件怎么装、第一条命令怎么敲。第2章是搭台布景：选力场、生成拓扑、装盒子、填溶剂、加离子，把一个能跑的物理体系立起来。第3章是走台：mdp 参数逐项过，能量最小化把结构里的小应力卸掉，NVT 与 NPT 两段平衡把温度、压强、密度带到目标状态。第4章是正式演出：grompp 编译、mdrun 推进、检查点续算，以及演出中途灯灭了怎么排查。第5章是复盘：轨迹预处理、RMSD 与 RMSF 判读、RDF 与氢键分析，把轨迹变成结论。第6章是加演场：普通模拟够不着的自由能与罕见事件，用 TI、伞形采样、metadynamics 去够；顺带把并行性能调到不浪费机时。第7章是汇演质检：两个完整案例从头到尾复盘，附一张可以直接抄进组内规范的质控清单。

三类读者怎么用这本册子

赶进度的同学走"主线速通"：导读、1.1、2.1 到 2.3、3.1 与 3.5、3.6、4.1、4.3、5.1 与 5.2，先跑通一个最小流程，其余章节回来补。带课题的研究生建议按顺序全读，尤其是第3章的六节——参数、积分、耦合、静电、最小化、平衡，这六章几乎所有的"跑崩"都出在这里。被别人留下的崩溃轨迹接手烂摊子的同学，直接翻 4.3 的排查手册和 7.2 的质控清单，对号入座后回到对应章节补课。

全册的数字口径是统一的：时间步长 2 fs，截断半径 1.0 到 1.2 nm，最小化收敛判据最大力低于 1000 kJ/mol/nm，温度耦合弛豫时间 0.5 ps，压强耦合弛豫时间 5 ps。你在自己的课题里当然可以改，但改之前要知道默认值为什么长这样——这正是每一节"演算"环节要做的事。书里的每段命令行会话和 Python 小演算都为"能跑通"而写：没有 GROMACS 环境时，纯 Python 的那部分照样能让你把算法原理摸透。