[论文解读] Use of Polynomial Interpolation to Reduce Bias and Uncertainty of Free Energy Estimates via Thermodynamic Integration

Conrad Shyu, F. Marty Ytreberg|arXiv (Cornell University)|Sep 5, 2008

Phase Equilibria and Thermodynamics参考文献 24被引用 2

一句话总结

本文提出使用非等距的λ值进行多项式插值，以减少通过热力学积分法估算自由能差时的偏差和不确定性。该方法在不需使用软核势或额外模拟的情况下，提高了溶剂化自由能计算的准确性和精确度，已在可解析求解的系统及小分子上得到验证。

ABSTRACT

This report details an approach to improve the accuracy and precision of free energy difference estimates using thermodynamic integration data (slope of the free energy with respect to the switching variable lambda) and its application to calculating solvation free energy. The central idea is to utilize polynomial fitting schemes to approximate the thermodynamic integration data to improve the accuracy and precision of the free energy difference estimates. In this report we introduce polynomial and spline interpolation techniques. Two systems with analytically solvable relative free energies are used to test the accuracy and precision of the interpolation approach (Shyu and Ytreberg, J Comput Chem 30: 2297-2304, 2009). We also use both interpolation and extrapolation methods to determine a small molecule salvation free energy. Our simulations show that, using such polynomial techniques and non-equidistant lambda values, the solvation free energy can be estimated with high accuracy without using soft-core scaling and separate simulations for Lennard-Jones and partial charges. The results from our study suggest these polynomial techniques, especially with use of non-equidistant lambda values, improve the accuracy and precision for dF estimates without demanding additional simulations. To allow researchers to immediately utilize these methods, free software and documentation is provided via this http URL

研究动机与目标

• 解决从热力学积分数据中估算自由能差时存在的偏差和高不确定性问题。
• 开发一种方法，以在不依赖软核缩放或额外模拟的情况下，提高溶剂化自由能计算的准确性和精确度。
• 在具有已知解析解的系统上，评估多项式与样条插值技术的性能。
• 证明非等距λ间距在提升估算质量方面的有效性。
• 提供开源软件，以支持研究社区立即采用该方法。

提出的方法

• 对热力学积分数据应用多项式拟合，其中自由能是切换参数λ的函数。
• 使用非等距的λ值以更好地捕捉自由能曲率较大的区域，从而提高插值精度。
• 实现插值与外推技术，以从模拟数据中估算自由能差。
• 在两个具有已知解析相对自由能的系统上验证该方法，以评估偏差和精确度。
• 与标准方法进行比较，量化在准确度提升和不确定性降低方面的改进。
• 发布免费软件及文档，以支持可重现性及社区采纳。

实验结果

研究问题

• RQ1 多项式插值能否显著减少从热力学积分中获得的自由能差估算的偏差和不确定性？
• RQ2 与等距间距相比，使用非等距λ值如何影响自由能估算的准确度和精确度？
• RQ3 该方法在多大程度上可消除对Lennard-Jones项和库仑项使用软核势及单独模拟的需求？
• RQ4 该插值方法在具有已知解析解的系统上的表现如何？
• RQ5 该方法能否可靠地应用于真实世界的溶剂化自由能计算，例如小分子？

主要发现

• 使用非等距λ值的多项式插值显著降低了自由能差估算的偏差和不确定性。
• 该方法在无需软核缩放或对Lennard-Jones项与部分电荷项进行单独模拟的情况下，实现了高精度的溶剂化自由能估算。
• 在具有已知解析自由能差的测试系统上，该方法优于标准方法。
• 非等距λ间距提高了插值的保真度，特别是在自由能快速变化的区域。
• 该技术仅通过一组经过优化λ间距的模拟，即可实现对溶剂化自由能的精确估算。
• 作者提供了免费的软件，使该方法能够立即被实现并广泛采用。