Weighted least square analysis method for free energy calculations. Issue 28 (30th October 2018)

• Record Type: Journal Article
• Authors: Hu, Dan; Guan, Xiaoqing; Wang, Yukun
• Is Part Of: Journal of computational chemistry. Volume 39:Issue 28(2018)
• Page Start: 2397
• Page End: 2404
• Publication Date: 2018-10-30
• Subjects: weighted least square -- free energy -- rare event -- umbrella sampling -- Welsam; Chemistry -- Data processing -- Periodicals; 542.85
• DOI: 10.1002/jcc.25580
• Languages: English
• ISSNs: 0192-8651

Abstract

Abstract : Free energy calculation is an efficient way for studying rare event dynamics. For a complex rare event dynamics, multiple reaction coordinates may be required to describe the transition path between equilibrium states. Theoretically, a one‐dimensional sampling along the transition path can provide sufficient information to calculate the potential of mean force (PMF) along the transition path. In the widely used free energy analysis method Wham, the sample data are divided into a series of bins to calculate PMF. However, bin segmentation in Wham is coupled with the umbrella potentials applied in each window, because each umbrella potential is assumed to have a close value for all sample points in each bin. This coupling makes it difficult to perform one‐dimensional bin segmentation along the transition path when multivariable umbrella potentials are used in sampling. Here, we develop a weighted least square analysis method (Welsam) to take the place of Wham for free energy analysis. In the new method Welsam, bin segmentation is decoupled from application of umbrella potentials. As a result, it becomes very convenient to perform one‐dimensional bin segmentation and calculate one‐dimensional PMF along the transition path. Our simulation results suggest that Welsam has a comparable statistical error with Wham. Furthermore, Welsam can be used to reduce waste of sample data obtained during exploration of reaction coordinates. © 2018 Wiley Periodicals, Inc. Abstract : A free energy analysis method based on weighted least squares is developed for free energy calculations. The new method has a comparable statistical error with Wham while it is more convenient for bin‐segmentation the Wham. In the figure, two‐dimensional potentials are used to perform umbrella sampling along the transition path shown by the red line. Bin‐segmentation along the transition path is then used to calculate the potential of mean forces.