第五章 轨迹分析

第五章 · 轨迹分析

章节摘要：轨迹是原料，结论要靠分析做出来。本章按分析工序讲三件事：先用 trjconv 把周期性边界切开的分子回卷、居中、去跳动（不做这步，后面全部白算）；再用 RMSD 与 RMSF 判断体系稳定了吗、哪些部位在动；最后用 RDF、氢键计数与取向序参数回答"结构长什么样"。每节配一段从零实现的 Python 演算（Kabsch 叠合、径向分布函数），让算法不被工具黑箱吞掉。读完你应当能把一条生轨迹加工成一组可写进论文的量，并知道每个量的适用边界。

学习目标

1. 说出 -pbc 的 mol、whole、nojump 三种模式各自修什么病，并能按分析目标选对组合。
2. 从零推导 RMSD 的计算流程（配准 + 逐原子距离统计），用 Kabsch 算法的小实现验证 gmx rms 的结果口径。
3. 区分 RMSD 与 RMSF 回答的两类问题，并识别 RMSD 曲线的"上升—平台—漂移"三段式含义。
4. 写出 RDF 的归一化公式并解释每一项的来历，读懂聚合物熔体与脂膜双层各自的特征峰。
5. 用氢键判据（距离加角度）算出氢键数目时间序列，并解释"平均氢键数"为什么不是强度的可靠指标。

核心概念速览

分析的第一原则是先修几何、再算物理：周期性边界让轨迹里的分子四分五裂，不先重装就统计，算出的每个量都带人为伪影。第二原则是每个量回答一个问题：稳定性问 RMSD、柔性分布问 RMSF、局域结构问 RDF、氢键网络问 hbond、取向序问序参数——不要指望一个量回答所有问题。

子章节导航

• 5.1 轨迹预处理与周期性边界——PBC 为什么会让分子"断成两截"，-pbc 三种模式与 -center 的组合拳，先修几何再谈分析。
• 5.2 RMSD 与 RMSF 判读——从 Kabsch 算法手写 RMSD 开始，讲到曲线三段式判读、参考结构与叠合组的选择陷阱。
• 5.3 RDF、氢键与取向分析——归一化公式的逐项推导，聚合物与脂膜的特征 RDF 峰，氢键判据的两要素，序参数在脂膜案例里的用法。

子章节之间的逻辑关系

三节是严格的流水线：5.1 的预处理产物是 5.2 与 5.3 的唯一合法输入（跳过预处理的分析结果应视为无效）；5.2 回答"整体稳不稳"这一前提问题，体系都不稳就谈不上结构统计；5.3 的结构量（RDF、氢键、序参数）则是 7.1 与 7.2 案例复盘的主角。第6章的增强采样同样建立在本章的预处理之上。

[5.1 修几何] --合法轨迹--> [5.2 稳定性] --通过--> [5.3 结构统计] --> [第6、7章]
    |                          |                       |
  -pbc 组合拳              RMSD/RMSF               RDF/氢键/序参数
    +---------- 跳过任何一步，后面的数字都不可信 ----------+
​

本章演算与脚本清单

关键数字速查

• whole → nojump → center：预处理的标准顺序，nojump 必须从第一帧连续处理。
• 0.001 nm：xtc 压缩精度，分析分辨率之上的安全线。
• Re/Rg ≈ 2.45：理想链环境的判据（熔体章会用）。
• 氢键两要素：距离阈值加角度阈值，缺一会把路过误判成键。

前置知识与后续延伸

前置：2.3 的周期性边界与最小镜像约定（本章一切操作的几何根源）、1.1 的时间平均思想、3.6 的平衡判据（分析只能对平衡后的轨迹做）。

延伸：5.2 的叠合思想在 6.2 的自由能面（构象聚类）里重现；5.3 的 RDF 是 7.1 聚合物案例的核心观测量，序参数是 7.2 脂膜案例的验收指标。学完本章，回头看 4.2 的"输出频率是分析设计的一部分"会有新的体会。