第一章 原理与排练准备

第一章 · 原理与排练准备

章节摘要：任何一次模拟在按下回车之前，都要先回答三个问题——程序在解什么方程、力从哪张势能表来、我手里的软件靠什么把方程变成轨迹。本章按这三问展开：先用一次双人原子的演算把 MD 的本质走通，再拆开力场的数学骨架看清每根弹簧与每份电荷，最后把 GROMACS 装起来、认识它的命令家族与文件管线。读完本章，你应当能对"模拟到底在算什么"给出一段不含糊话的回答，并让 gmx 命令在你自己的终端里跑出第一行输出。

学习目标

读完本章并做完演算，你应当能够：

1. 不查资料地说出 MD 的输入与输出分别是什么，并解释为什么"轨迹"本身还不是结论，统计平均才是。
2. 从拉格朗日量或牛顿第二定律出发，复述积分器在每一步要做的动作，并说出时间步长受什么物理量限制。
3. 写出势能函数的六大件：键长项、键角项、二面角项、非正常二面角项、Lennard-Jones 项与静电项，并解释各自守着分子的哪种几何。
4. 给定两种原子的 LJ 参数，手算交叉相互作用的 sigma 与 epsilon（两种混合规则各算一遍）。
5. 判断自己的机器上哪条安装路线最省事，用 gmx --version 确认版本、线程与 GPU 支持，并说出每个 GROMACS 工具"读什么、吐什么"。

核心概念速览

分子模拟的全部戏剧都发生在一张势能面上：力场给出地形的等高线图，积分器负责让每个原子沿着坡往下滚，轨迹是滚动的录像，分析是录像的观后感。地形画错了，滚得再久也到不了真实世界。

本章的知识链条可以用一句话串起来：物理假设 → 势能函数 → 力 → 轨迹 → 统计量。缺任何一环，后面的章节都会悬空。

子章节导航

• 1.1 MD 模拟到底在算什么——不写公式堆砌，直接演算两个原子如何被牛顿方程推着走，顺便算清楚自由度、温度与动能的账。本节是全册的地基。
• 1.2 势能函数与力场的数学骨架——把六大势能项逐个拆开，弄懂每种函数形式在守什么、参数从哪来、混合规则怎么算。为第2章的力场选型做铺垫。
• 1.3 GROMACS 定位与第一个命令——软件从哪来、为什么快、怎么装、命令家族怎么分工，并把"文件 → tpr → 轨迹"的管线走一遍。为第4章的实操打底。

子章节之间的逻辑关系

三节的顺序不能颠倒：不先明白"算什么"（1.1），就没法理解力场为什么长成函数表的样子（1.2）；不先理解力场的输入地位，就看不懂 GROMACS 为什么要把结构、拓扑、参数编译成一个 tpr 才肯开工（1.3）。

[物理图景 1.1] ——势能从哪来——> [力场骨架 1.2] ——参数落进文件——> [软件与管线 1.3]
      |                                  |                                |
   牛顿方程                         六大势能项                     结构+拓扑+mdp
      |                                  |                                |
      +--------------- 三者合流 -------------------> 编译 tpr，准备搭台（第二章）
​

本章演算与脚本清单

关键数字速查

本章贡献给全册的"九九表"：室温热能 k_B T ≈ 2.48 kJ/mol；玻尔兹曼常数 0.008314 kJ/(mol·K)；C–C 键拉伸 0.01 nm 约 28 kJ/mol；O–H 振动周期约 10 fs（3.2 步长推导的源头）；经典模拟典型规模 10⁴ 至 10⁷ 原子、皮秒到微秒可及。这些数字会在后续每章反复出场——第一次遇到就手推一遍，后面引用时才有底气。

前置知识与后续延伸

前置：大学物理的牛顿力学与基础热力学；一点 Python（能读懂 numpy 数组操作即可，演算脚本逐行有注释）。不要求量子化学背景——经典 MD 恰恰是靠"把电子积分进参数里"才换来速度，这件事在 1.2 会展开讲。

延伸：1.1 末尾留下的"时间步长从哪来"问题在第3章 3.2 正面回答；1.2 的势能项拆解在第2章直接变成选型与参数核对的动手内容；1.3 的文件管线在第4章变成完整的生产运行。带着问题往下走，比把定义背下来有用得多。