2. 分子动力学基础回顾与增强采样必要性
在进入增强采样方法之前，有必要把标准分子动力学（MD）的“工作边界”重新梳理清楚。增强采样并不是对 MD 的否定，而是对 MD 在有限时间内无法充分访问构象空间这一事实的工程化补救。理解力场、积分算法、系综、周期性边界、时间步长这些基础要素，才能准确判断一个体系究竟卡在哪里：是势能面本身有高垒，还是积分步长太小导致模拟时长不足，抑或是遍历性破缺让轨迹长期困在某个局部盆地。
力场与积分算法：决定“能走多远”的两条腿
力场给出原子间的势能函数 U(r)，常见形式包括键伸缩、键角弯曲、二面角扭转、非键范德华与静电项。以 AMBER/CHARMM 类全原子力场为例，键与角多用谐振子近似，二面角用周期函数，非键部分用 Lennard-Jones 加库仑项。力场决定了势能面的形状，也决定了势垒高度与鞍点位置；如果力场把某个二面角旋转势垒估得过高，那么即使采样算法再强，也会人为放大稀有事件难度。
积分算法把牛顿方程离散化。最常用的是 velocity Verlet：
v(t+Δt/2) = v(t) + (Δt/2m)·F(t)
r(t+Δt) = r(t) + Δt·v(t+Δt/2)
F(t+Δt) = -∇U(r(t+Δt))
v(t+Δt) = v(t+Δt/2) + (Δt/2m)·F(t+Δt)它时间可逆、辛结构好，能量漂移小。Δt 的选取受最高频振动限制：含氢键伸缩的体系通常取 1 fs；若用 SHAKE/LINCS 约束含氢键，可取 2 fs；配合氢质量重分配（HMR）可到 4 fs。Δt 过大，积分会失稳，表现为能量突增、原子飞散；Δt 过小，则同样机时只能覆盖更短物理时间，直接压缩采样量。
系综、周期性边界与采样极限
系综决定哪些热力学量被控制。NVE 能量守恒，适合检验积分稳定性；NVT 用 Nosé–Hoover 或 Langevin 热浴控制温度；NPT 再加 Parrinello–Rahman 或 Berendsen 控压，用于接近实验条件。系综选择错误会让结果不可比，例如在 NVT 下研究压力相关构象变化就不合适。
周期性边界（PBC）用最小镜像约定消除表面效应，但要求非键截断小于盒子半边长。若盒子过小，分子可能与自身镜像作用，产生虚假周期性；若截断处理不当，长程静电需用 PME 等格点求和。PBC 不改变采样极限，它只保证体相条件正确。
采样极限来自两个层面：一是时间尺度，二是遍历性。一个蛋白质在 300 K 下可能需要在微秒到毫秒量级才能完成一次大尺度构象转变，而常规全原子 MD 在 GPU 上通常只能跑到百纳秒到微秒。更关键的是遍历性破缺：轨迹可能长时间停留在某个亚稳态，无法按玻尔兹曼权重访问其他盆地。此时增加模拟时长收益递减，因为轨迹并没有“均匀”探索，而是被势垒和摩擦困住。
稀有事件与遍历性破缺的实际表现
稀有事件通常指在可观测时间内发生概率极低的构象跃迁，例如蛋白质别构口袋打开、配体从深结合位点逃逸、离子跨膜、肽段从 α 螺旋到 β 折叠的转变。它们的共同特征是：自由能垒远高于 kBT，跃迁等待时间远长于模拟时长。遍历性破缺则表现为轨迹在某个盆地内反复振荡，却始终不跨越势垒。判断依据可以看二面角时间序列、RMSD 分布、聚类占比：如果某个簇占据 90% 以上且长时间不衰减，就说明采样没有充分混合。
以一个含 20 个残基的肽段为例，假设在 300 K、NVT 下跑 500 ns，使用 2 fs 步长、SHAKE 约束，每 10 ps 保存一帧。若分析 φ/ψ 二面角分布，发现某残基长期停留在 α 区，而 β 区几乎无采样，则说明该二面角势垒未被跨越。此时继续延长到 1 μs 可能仍无改善，因为等待时间可能是指数级的。增强采样的必要性就在这里：通过加偏置、调温、交换副本或沿反应坐标加速，主动提高稀有事件发生率，再通过重加权恢复无偏玻尔兹曼分布。
工程判断：先诊断，再选增强策略
面对一个具体体系，建议先做基础诊断，而不是直接上增强采样。诊断顺序可以是：检查力场与体系是否匹配；确认 Δt 与约束设置是否合理；检查 PBC 与截断是否产生伪影；观察能量、温度、压力是否稳定；再看关键反应坐标的采样是否充分。若标准 MD 已能覆盖目标事件，则无需增强；若不能，再根据事件类型选择方法：沿已知反应坐标可用伞形采样或 metadynamics；无明确反应坐标可用副本交换或模拟退火；需要加速慢自由度可用加速分子动力学或高斯加速分子动力学。