Weighted least square analysis method for free energy calculations

Dan Hu, Xiaoqing Guan, Yukun Wang

Research output: Contribution to journalArticle

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.

Original languageEnglish (US)
Pages (from-to)2397-2404
Number of pages8
JournalJournal of Computational Chemistry
Volume39
Issue number28
DOIs
StatePublished - Oct 30 2018

Fingerprint

Weighted Least Squares
Bins
Free energy
Free Energy
Segmentation
Path
Rare Events
Calculate
Sampling
Sample point
Equilibrium State
Sufficient
Series
Simulation

Keywords

  • free energy
  • rare event
  • umbrella sampling
  • weighted least square
  • Welsam

ASJC Scopus subject areas

  • Chemistry(all)
  • Computational Mathematics

Cite this

Weighted least square analysis method for free energy calculations. / Hu, Dan; Guan, Xiaoqing; Wang, Yukun.

In: Journal of Computational Chemistry, Vol. 39, No. 28, 30.10.2018, p. 2397-2404.

Research output: Contribution to journalArticle

Hu, Dan ; Guan, Xiaoqing ; Wang, Yukun. / Weighted least square analysis method for free energy calculations. In: Journal of Computational Chemistry. 2018 ; Vol. 39, No. 28. pp. 2397-2404.
@article{6ff8c84da5d5401ab1a46bcbd48f849a,
title = "Weighted least square analysis method for free energy calculations",
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.",
keywords = "free energy, rare event, umbrella sampling, weighted least square, Welsam",
author = "Dan Hu and Xiaoqing Guan and Yukun Wang",
year = "2018",
month = "10",
day = "30",
doi = "10.1002/jcc.25580",
language = "English (US)",
volume = "39",
pages = "2397--2404",
journal = "Journal of Computational Chemistry",
issn = "0192-8651",
publisher = "John Wiley and Sons Inc.",
number = "28",

}

TY - JOUR

T1 - Weighted least square analysis method for free energy calculations

AU - Hu, Dan

AU - Guan, Xiaoqing

AU - Wang, Yukun

PY - 2018/10/30

Y1 - 2018/10/30

N2 - 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.

AB - 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.

KW - free energy

KW - rare event

KW - umbrella sampling

KW - weighted least square

KW - Welsam

UR - http://www.scopus.com/inward/record.url?scp=85055696748&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85055696748&partnerID=8YFLogxK

U2 - 10.1002/jcc.25580

DO - 10.1002/jcc.25580

M3 - Article

C2 - 30374999

AN - SCOPUS:85055696748

VL - 39

SP - 2397

EP - 2404

JO - Journal of Computational Chemistry

JF - Journal of Computational Chemistry

SN - 0192-8651

IS - 28

ER -