[论文解读] The hydration free energy of water

Gerhard Hummer, Lawrence R. Pratt+1 more|ArXiv.org|Jul 20, 1995

Fluid dynamics and aerodynamics studies被引用 14

一句话总结

本研究利用涨落理论与多种模拟方法，计算了水的溶剂化自由能随电荷的变化关系，发现其显著偏离连续介质模型所预测的二次函数行为。作者在重叠直方图、接受率和几何插值方法之间均得到一致结果，验证了Rick和Berne的热力学积分结果，同时揭示了非高斯的静电涨落，从而否定了简单线性模型的有效性。

ABSTRACT

We study the chemical potential of water as a function of charge based on perturbation theory. By calculating the electrostatic-energy fluctuations of two states (fully charged and uncharged) we are able to determine accurate values for the dependence of the chemical potential on charge. We find identical results for the chemical-potential difference of fully charged and uncharged water from overlapping-histogram and acceptance-ratio methods and by smoothly connecting the curves of direct exponential averages. Our results agree with those of Rick and Berne (J. Am. Chem. Soc., 1994, 116, 3949) with respect to both the chemical-potential difference and its dependence on the charge coupling parameter. We observe significant deviations from simple Gaussian-fluctuation statistics. The dependence on the coupling parameter is not quadratic, as would be inferred from linear continuum models of electrostatics.

研究动机与目标

• 通过平衡模拟，准确确定水的电化学势随电荷变化的静电贡献。
• 检验二级微扰理论和高斯涨落统计在模拟水的溶剂化自由能中的有效性。
• 评估在水的电静电模拟中有限尺寸效应和自相互作用校正的影响。
• 在相同体系上比较多种自由能估算技术——重叠直方图法、接受率法和几何插值法。
• 与Rick和Berne（1994年）的热力学积分数据进行结果验证。

提出的方法

• 使用微扰理论，将带电与不带电水之间的化学势差与静电能差的指数平均值关联起来。
• 应用电势分布定理，通过固定电荷态下的构型平均计算过剩化学势。
• 采用Ewald求和法与广义反应场（GRF）方法处理长程静电作用，并一致地包含自相互作用校正。
• 使用Bennett的重叠直方图法与接受率法，从两个电荷态下的独立模拟中估算自由能差。
• 采用一种几何方法，通过电荷耦合参数λ平滑连接不同电荷态下的直接指数平均值。
• 在不同尺寸的系统（N = 64 和 N = 256）上进行模拟，以评估有限尺寸效应并验证收敛性。

实验结果

研究问题

• RQ1 水的溶剂化自由能是否如Born模型所预测的那样，对电荷耦合参数呈现二次依赖关系？
• RQ2 基于涨落的方法（如重叠直方图法与接受率法）在估算带电与不带电水之间化学势差时的准确性如何？
• RQ3 有限尺寸效应以及自相互作用项的缺失在多大程度上扭曲了水模拟中的自由能估算？
• RQ4 当包含自相互作用校正时，不同静电模型（Ewald vs. GRF）能否获得一致的结果？
• RQ5 静电能中的非高斯涨落如何影响二级微扰理论在水溶液体系中的有效性？

主要发现

• 完全带电与不带电水之间的化学势差为 Δμ^ex = 35.56 kJ mol⁻¹，与Rick和Berne的热力学积分结果（35.52 kJ mol⁻¹）高度一致。
• 重叠直方图法、接受率法与几何插值法对化学势差的估算结果完全一致，证实了方法的一致性。
• 观察到显著偏离电荷耦合参数λ的二次依赖关系，否定了简单线性连续介质电静电模型的有效性。
• 包含自相互作用项可减少系统尺寸依赖性，并改善Ewald与GRF模拟结果的一致性，尤其在小系统（N = 64）中效果显著。
• 若不包含自相互作用校正，Ewald与GRF结果之间的差异可达约1.7 kJ mol⁻¹，凸显此类校正的必要性。
• 本研究证明，在两个采样充分的状态之间进行插值，比高阶累积量展开更为准确，后者因分布尾部采样不足而存在缺陷。