TY - JOUR

T1 - Estimation of free-energy differences from computed work distributions

T2 - An application of Jarzynski's equality

AU - Echeverria, Ignacia

AU - Amzel, L. Mario

N1 - Copyright:
Copyright 2015 Elsevier B.V., All rights reserved.

PY - 2012/9/13

Y1 - 2012/9/13

N2 - Equilibrium free-energy differences can be computed from nonequilibrium molecular dynamics (MD) simulations using Jarzynski's equality (Jarzynski, C. Phys. Rev. Lett.1997, 78, 2690) by combining a large set of independent trajectories (path ensemble). Here we present the multistep trajectory combination (MSTC) method to compute free-energy differences, which by combining trajectories significantly reduces the number of trajectories necessary to generate a representative path ensemble. This method generates well-sampled work distributions, even for large systems, by combining parts of a relatively small number of trajectories carried out in steps. To assess the efficiency of the MSTC method, we derived analytical expressions and used them to compute the bias and the variance of the free-energy estimates along with numerically calculated values. We show that the MSTC method significantly reduces both the bias and variance of the free-energy estimates compared to the estimates obtained using single-step trajectories. In addition, because in the MSTC method the process is divided into steps, it is feasible to compute the reverse transition. By combining the forward and reverse processes, the free-energy difference can be computed using the Crooks' fluctuation theorem (Crooks, G. E. J. Stat. Phys.1998, 90, 1481 and Crooks, G. E. Phys. Rev. E2000, 61, 2361) or Bennett's acceptance ratio (Bennett, C. H. J. Comput. Phys. 1976, 22, 245), which further reduces the bias and variance of the estimates.

AB - Equilibrium free-energy differences can be computed from nonequilibrium molecular dynamics (MD) simulations using Jarzynski's equality (Jarzynski, C. Phys. Rev. Lett.1997, 78, 2690) by combining a large set of independent trajectories (path ensemble). Here we present the multistep trajectory combination (MSTC) method to compute free-energy differences, which by combining trajectories significantly reduces the number of trajectories necessary to generate a representative path ensemble. This method generates well-sampled work distributions, even for large systems, by combining parts of a relatively small number of trajectories carried out in steps. To assess the efficiency of the MSTC method, we derived analytical expressions and used them to compute the bias and the variance of the free-energy estimates along with numerically calculated values. We show that the MSTC method significantly reduces both the bias and variance of the free-energy estimates compared to the estimates obtained using single-step trajectories. In addition, because in the MSTC method the process is divided into steps, it is feasible to compute the reverse transition. By combining the forward and reverse processes, the free-energy difference can be computed using the Crooks' fluctuation theorem (Crooks, G. E. J. Stat. Phys.1998, 90, 1481 and Crooks, G. E. Phys. Rev. E2000, 61, 2361) or Bennett's acceptance ratio (Bennett, C. H. J. Comput. Phys. 1976, 22, 245), which further reduces the bias and variance of the estimates.

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

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

U2 - 10.1021/jp300527q

DO - 10.1021/jp300527q

M3 - Article

C2 - 22849340

AN - SCOPUS:84866395823

VL - 116

SP - 10986

EP - 10995

JO - Journal of Physical Chemistry B

JF - Journal of Physical Chemistry B

SN - 1520-6106

IS - 36

ER -