分子动力学模拟氢键

结合案例图分析，本文以这个含羟基体系氢键网络项目为例，讲氢键定义和 RDF 怎么用。这是一个含羟基体系氢键网络研究项目。客户研究水溶液或含羟基分子的氢键动力学，想了解平均氢键数、寿命和溶剂化结构。第一张案例图是 20 ns 内氢键数量随时间变化，平均 ⟨n_HB⟩≈6。第二张案例图是 O-H···O 的 RDF，第一峰约 1.9 Å，第二峰约 2.8 Å。

项目过程：含羟基体系氢键网络项目的执行链路

方法上氢键判据有几何法（距离+角度）和能量法（相互作用能阈值）两类。几何法快、可复现、易与文献对比，所以我固定用几何标准；能量法对力场敏感、阈值难定。我的标准：供体-受体距离 < 3.5 Å、D-H···A 角度 > 120°，这是文献最常用的组合，偏差可控。 我的流程：建系平衡后跑生产性 MD（2 fs 步长、298 K、总 20 ns），每隔 10 ps 帧用几何标准分析氢键并算 RDF。客户最初只要平均氢键数，我额外给了 RDF 和分布，因为单一个数说明不了结构。

案例分析：从 RDF 到溶剂化结构

图里读出的数要落地。第一步确认氢键数已收敛——第一张图 2 ns 达峰值后波动，我对后 15 ns 算均值 ⟨n_HB⟩≈6.0±0.5，标准差小才可信。我对 RDF 第一峰位置算统计误差，1.9±0.05 Å 与实验 X 射线值吻合，是判断力场质量的硬指标；若第一峰偏到 2.1 Å 说明力场氢键偏软。 第二步读第二峰 2.8 Å 为第二溶剂化壳，两峰间距反映四面体有序度。误差量级上，角度阈值 120° 比 150° 把更多弱作用算入，第一峰高度差约 15%。

总结教训：氢键定义与 RDF 的坑

做的过程中我踩过实打实的坑。第一坑是阈值：我画过平均氢键数随截断距离的曲线——距离阈值从 3.0 Å 放宽到 3.5 Å 时 ⟨n_HB⟩ 从 5.1 升到 6.0，3.8 Å 又跳到 6.8，说明阈值每动 0.1 Å 平均数量就变 0.3，必须写死标准。 第二坑是 RDF 第一峰位置，角度阈值 120° 固定并标注，否则不同文献没法比。第三坑是溶质参与，纯水 ⟨n_HB⟩≈3.4，含羟基溶质会升到 6，必须说明体系。

项目经验：氢键交付的可复用清单

作为可复用的项目经验，第一，RDF 第一峰位置和高度是力场质量判据，必须和实验比对。第二，要研究寿命需算时间相关函数 C_HB(t)，不只平均数量。第三，几何标准（距离 3.5 Å、角度 120°）必须写清并附误差带。客户据此理解溶质溶剂化结构。说到底，氢键项目最容易”同一个体系三家数”的，就是定义没对齐；几何标准写清楚，数据才进得了文献对比。

分子动力学模拟氢键

计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验

需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。