第三章 参数最小化与平衡

第三章 · 参数、最小化与平衡

章节摘要：走台阶段决定模拟的成败一半以上：mdp 文件里的每一行参数都是一个物理承诺，能量最小化把构建期的高能关节卸掉，两段系综平衡把体系带到目标温压状态。本章把最常被含糊带过的六个主题一次讲透——参数逐项精读、积分算法与步长上限的推导、恒温恒压算法的取舍对照、PME 静电的空间分工、最小化收敛判据、平衡曲线的判读与放行标准。每节都有演算或脚本，读完你应当能对"我的参数为什么这么设"给出有依据的回答，而不是"教程里就这么写的"。

学习目标

1. 拿到一份 mdp 文件能逐行说出每个参数的物理含义与默认值的来历，并能用脚本对比两份 mdp 的差异。
2. 推导积分步长的稳定性上限：从最快振动周期（如 O–H 键约 10 fs）出发，说出为什么约束后可以放心用 2 fs。
3. 说出三类恒温器与三类恒压器各自的可逆性、涨落保真度与适用阶段，并解释为什么 Berendsen 只能用于平衡前段。
4. 手算一个盒子在给定网格间距下的 PME 格点数，解释实空间与倒易空间的误差分工。
5. 用最大力判据（低于 1000 kJ/mol/nm）放行能量最小化，用密度、温度、能量三条曲线的平稳性放行 NVT 与 NPT 平衡——每条判据都知道它在防什么病。

核心概念速览

走台的全部逻辑是一条因果链：参数定物理（mdp）→ 步长定稳定（积分器）→ 耦合定系综（恒温恒压）→ 静电定精度（PME）→ 最小化定起点 → 平衡定放行。任何一环偷懒，第4章的正式演出都会以崩溃或坏数据的方式讨回来。

子章节导航

• 3.1 mdp 参数逐项精读——一份带注释的完整 mdp 文件加一个差异对比脚本，把参数从"咒语"还原成"决策"。
• 3.2 积分算法与步长选择——Velocity Verlet 与 leapfrog 的同构关系、稳定性上限的推导、约束为什么能救步长，附 Euler 与 Verlet 的能量漂移对照演算。
• 3.3 温度与压力耦合对照——六种耦合算法一张表，涨落保真度的含义，热容演算，以及"平衡期用松耦合、生产期用严格系综"的阶段纪律。
• 3.4 长程静电与 PME——为什么截断静电是罪、Ewald 怎么拆、PME 网格参数怎么配、格点数怎么算。
• 3.5 能量最小化实操——最速下降的流程会话、最大力判据的含义、收敛曲线判读与收敛检查脚本。
• 3.6 NVT 与 NPT 平衡实操——两段平衡各管什么、位置限制怎么开、放行判据怎么量化，附一段收敛性检查脚本。

子章节之间的逻辑关系

六节的依赖是单向的：参数文件里的 integrator 与 dt 指向 3.2 的算法选择；tcoupl 与 pcoupl 指向 3.3 的系综设计；coulombtype 指向 3.4 的静电方案；然后才谈得上 3.5 的最小化（它用的正是前四节定下的参数）与 3.6 的两段平衡（它检查的曲线来自 3.3 选的耦合算法）。学的时候可以跳跃，用的时候必须成链。

[3.1 参数] --integrator/dt--> [3.2 积分]      [3.4 PME]
     |---tcoupl/pcoupl---->  [3.3 耦合]           |
     +---coulombtype---------+---------------------+
                  v
          [3.5 最小化] --> [3.6 NVT/NPT] --> 放行判据 --> [第4章]
​

本章演算与脚本清单

关键数字速查

• dt = 0.002 ps，前提 constraints = h-bonds 全约束。
• emtol = 1000 kJ/mol/nm：EM 放行线。
• tau_t = 0.5 ps、tau_p = 5 ps、4.5e-5 bar⁻¹：耦合三件套。
• rcoulomb = 1.0、fourierspacing = 0.12、pme_order = 4：PME 三件套。
• 密度漂移 < 1%、漂移涨落比 < 0.1：NPT 放行判据。

前置知识与后续延伸

前置：1.1 的自由度与温度换算（3.2、3.3 直接引用）、1.2 的静电项（3.4 的主角）、2.3 装配完成的体系（本章一切操作的输入）。

延伸：3.5 的收敛判据在 4.3 排错手册里是第一道分诊线；3.6 的放行标准在 7.2 的质控清单里被扩展成完整表格；3.3 的涨落保真度概念在 6.1 的自由能计算里再次出现——那边对系综正确性的要求更苛刻。