BAR-based optimum adaptive sampling regime for variance minimization in alchemical transformation

The nonequilibrium stratification

Xiaohui Wang, Xingzhao Tu, John Zhang, Zhaoxi Sun

Research output: Contribution to journalArticle

Abstract

Following the previously proposed equilibrate-state sampling based adaptive sampling regime Optimum Bennett Acceptance Ratio (OBAR), we introduce its nonequilibrium extension, the Optimum Crooks' Equation (OCE) in the current work. The efficiency of the NonEquilibrium Work (NEW) stratification is improved by adaptively manipulating the significance of each nonequilibrium realization followed by importance sampling. As is exhibited in the equilibrium case, the nonequilibrium extension outperforms the simple equal time rule used in nonequilibrium stratification in the sense of minimizing the total variance of the free energy estimate. The speedup of this non-equal time rule is more than 1-fold. The Time Derivative of total Variance (TDV) proposed for the OBAR criterion is extended to determine the importance of each nonequilibrium transformation, which is linearly dependent on the variance. The TDV in the nonequilibrium case gives a totally different importance rank from the standard errors of the free energy differences and OBAR TDV due to the duration of nonequilibrium pulling being added into the OCE equation. The performance of the OCE workflow is demonstrated in the solvation of several small molecules with a series of lambda increments and relaxation times between successive perturbations. To the best of our knowledge, such a nonequilibrium adaptive sampling regime in alchemical transformation is unprecedented.

Original languageEnglish (US)
Pages (from-to)2009-2021
Number of pages13
JournalPhysical Chemistry Chemical Physics
Volume20
Issue number3
DOIs
StatePublished - Jan 1 2018

Fingerprint

stratification
sampling
Sampling
Derivatives
Free energy
optimization
Importance sampling
acceptability
Solvation
Relaxation time
free energy
Molecules
pulling
solvation
relaxation time
perturbation
estimates
molecules

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

BAR-based optimum adaptive sampling regime for variance minimization in alchemical transformation : The nonequilibrium stratification. / Wang, Xiaohui; Tu, Xingzhao; Zhang, John; Sun, Zhaoxi.

In: Physical Chemistry Chemical Physics, Vol. 20, No. 3, 01.01.2018, p. 2009-2021.

Research output: Contribution to journalArticle

@article{cbc0653feebb48ac9560b1e882f55a5d,
title = "BAR-based optimum adaptive sampling regime for variance minimization in alchemical transformation: The nonequilibrium stratification",
abstract = "Following the previously proposed equilibrate-state sampling based adaptive sampling regime Optimum Bennett Acceptance Ratio (OBAR), we introduce its nonequilibrium extension, the Optimum Crooks' Equation (OCE) in the current work. The efficiency of the NonEquilibrium Work (NEW) stratification is improved by adaptively manipulating the significance of each nonequilibrium realization followed by importance sampling. As is exhibited in the equilibrium case, the nonequilibrium extension outperforms the simple equal time rule used in nonequilibrium stratification in the sense of minimizing the total variance of the free energy estimate. The speedup of this non-equal time rule is more than 1-fold. The Time Derivative of total Variance (TDV) proposed for the OBAR criterion is extended to determine the importance of each nonequilibrium transformation, which is linearly dependent on the variance. The TDV in the nonequilibrium case gives a totally different importance rank from the standard errors of the free energy differences and OBAR TDV due to the duration of nonequilibrium pulling being added into the OCE equation. The performance of the OCE workflow is demonstrated in the solvation of several small molecules with a series of lambda increments and relaxation times between successive perturbations. To the best of our knowledge, such a nonequilibrium adaptive sampling regime in alchemical transformation is unprecedented.",
author = "Xiaohui Wang and Xingzhao Tu and John Zhang and Zhaoxi Sun",
year = "2018",
month = "1",
day = "1",
doi = "10.1039/c7cp07573a",
language = "English (US)",
volume = "20",
pages = "2009--2021",
journal = "Physical Chemistry Chemical Physics",
issn = "1463-9076",
publisher = "Royal Society of Chemistry",
number = "3",

}

TY - JOUR

T1 - BAR-based optimum adaptive sampling regime for variance minimization in alchemical transformation

T2 - The nonequilibrium stratification

AU - Wang, Xiaohui

AU - Tu, Xingzhao

AU - Zhang, John

AU - Sun, Zhaoxi

PY - 2018/1/1

Y1 - 2018/1/1

N2 - Following the previously proposed equilibrate-state sampling based adaptive sampling regime Optimum Bennett Acceptance Ratio (OBAR), we introduce its nonequilibrium extension, the Optimum Crooks' Equation (OCE) in the current work. The efficiency of the NonEquilibrium Work (NEW) stratification is improved by adaptively manipulating the significance of each nonequilibrium realization followed by importance sampling. As is exhibited in the equilibrium case, the nonequilibrium extension outperforms the simple equal time rule used in nonequilibrium stratification in the sense of minimizing the total variance of the free energy estimate. The speedup of this non-equal time rule is more than 1-fold. The Time Derivative of total Variance (TDV) proposed for the OBAR criterion is extended to determine the importance of each nonequilibrium transformation, which is linearly dependent on the variance. The TDV in the nonequilibrium case gives a totally different importance rank from the standard errors of the free energy differences and OBAR TDV due to the duration of nonequilibrium pulling being added into the OCE equation. The performance of the OCE workflow is demonstrated in the solvation of several small molecules with a series of lambda increments and relaxation times between successive perturbations. To the best of our knowledge, such a nonequilibrium adaptive sampling regime in alchemical transformation is unprecedented.

AB - Following the previously proposed equilibrate-state sampling based adaptive sampling regime Optimum Bennett Acceptance Ratio (OBAR), we introduce its nonequilibrium extension, the Optimum Crooks' Equation (OCE) in the current work. The efficiency of the NonEquilibrium Work (NEW) stratification is improved by adaptively manipulating the significance of each nonequilibrium realization followed by importance sampling. As is exhibited in the equilibrium case, the nonequilibrium extension outperforms the simple equal time rule used in nonequilibrium stratification in the sense of minimizing the total variance of the free energy estimate. The speedup of this non-equal time rule is more than 1-fold. The Time Derivative of total Variance (TDV) proposed for the OBAR criterion is extended to determine the importance of each nonequilibrium transformation, which is linearly dependent on the variance. The TDV in the nonequilibrium case gives a totally different importance rank from the standard errors of the free energy differences and OBAR TDV due to the duration of nonequilibrium pulling being added into the OCE equation. The performance of the OCE workflow is demonstrated in the solvation of several small molecules with a series of lambda increments and relaxation times between successive perturbations. To the best of our knowledge, such a nonequilibrium adaptive sampling regime in alchemical transformation is unprecedented.

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

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

U2 - 10.1039/c7cp07573a

DO - 10.1039/c7cp07573a

M3 - Article

VL - 20

SP - 2009

EP - 2021

JO - Physical Chemistry Chemical Physics

JF - Physical Chemistry Chemical Physics

SN - 1463-9076

IS - 3

ER -