TY - JOUR
T1 - Simultaneous computation of dynamical and equilibrium information using a weighted ensemble of trajectories
AU - Suárez, Ernesto
AU - Lettieri, Steven
AU - Zwier, Matthew C.
AU - Stringer, Carsen A.
AU - Subramanian, Sundar Raman
AU - Chong, Lillian T.
AU - Zuckerman, Daniel M.
PY - 2014/7/8
Y1 - 2014/7/8
N2 - Equilibrium formally can be represented as an ensemble of uncoupled systems undergoing unbiased dynamics in which detailed balance is maintained. Many nonequilibrium processes can be described by suitable subsets of the equilibrium ensemble. Here, we employ the "weighted ensemble" (WE) simulation protocol [Huber and Kim, Biophys. J. 1996, 70, 97-110] to generate equilibrium trajectory ensembles and extract nonequilibrium subsets for computing kinetic quantities. States do not need to be chosen in advance. The procedure formally allows estimation of kinetic rates between arbitrary states chosen after the simulation, along with their equilibrium populations. We also describe a related history-dependent matrix procedure for estimating equilibrium and nonequilibrium observables when phase space has been divided into arbitrary non-Markovian regions, whether in WE or ordinary simulation. In this proof-of-principle study, these methods are successfully applied and validated on two molecular systems: explicitly solvated methane association and the implicitly solvated Ala4 peptide. We comment on challenges remaining in WE calculations.
AB - Equilibrium formally can be represented as an ensemble of uncoupled systems undergoing unbiased dynamics in which detailed balance is maintained. Many nonequilibrium processes can be described by suitable subsets of the equilibrium ensemble. Here, we employ the "weighted ensemble" (WE) simulation protocol [Huber and Kim, Biophys. J. 1996, 70, 97-110] to generate equilibrium trajectory ensembles and extract nonequilibrium subsets for computing kinetic quantities. States do not need to be chosen in advance. The procedure formally allows estimation of kinetic rates between arbitrary states chosen after the simulation, along with their equilibrium populations. We also describe a related history-dependent matrix procedure for estimating equilibrium and nonequilibrium observables when phase space has been divided into arbitrary non-Markovian regions, whether in WE or ordinary simulation. In this proof-of-principle study, these methods are successfully applied and validated on two molecular systems: explicitly solvated methane association and the implicitly solvated Ala4 peptide. We comment on challenges remaining in WE calculations.
UR - http://www.scopus.com/inward/record.url?scp=84904125478&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84904125478&partnerID=8YFLogxK
U2 - 10.1021/ct401065r
DO - 10.1021/ct401065r
M3 - Article
AN - SCOPUS:84904125478
SN - 1549-9618
VL - 10
SP - 2658
EP - 2667
JO - Journal of Chemical Theory and Computation
JF - Journal of Chemical Theory and Computation
IS - 7
ER -