2013, 2013(special): 171-181. doi: 10.3934/proc.2013.2013.171

Analysis of the accelerated weighted ensemble methodology

1. 

Mechanical Engineering Department, Stanford University, CA, United States

2. 

Computer Science and Engineering, University of Notre Dame, IN, United States, United States

3. 

Mechanical Engineering Department and Institute for Computational and Mathematical Engineering, Stanford University, CA, United States

Received  September 2012 Revised  September 2013 Published  November 2013

The main issue addressed in this note is the study of an algorithm to accelerate the computation of kinetic rates in the context of molecular dynamics (MD). It is based on parallel simulations of short-time trajectories and its main components are: a decomposition of phase space into macrostates or cells, a resampling procedure that ensures that the number of parallel replicas (MD simulations) in each macro-state remains constant, the use of multiple populations (colored replicas) to compute multiple rates (e.g., forward and backward rates) in one simulation. The method leads to enhancing the sampling of macro-states associated to the transition states, since in traditional MD these are likely to be depleted even after short-time simulations. By allowing parallel replicas to carry different probabilistic weights, the number of replicas within each macro-state can be maintained constant without introducing any bias. The empirical variance of the estimated reaction rate, defined as a probability flux, is expectedly diminished. This note is a first attempt towards a better mathematical and numerical understanding of this method. It provides first a mathematical formalization of the notion of colors. Then, the link between this algorithm and a set of closely related methods having been proposed in the literature within the past few years is discussed. Lastly, numerical results are provided that illustrate the efficiency of the method.
Citation: Ronan Costaouec, Haoyun Feng, Jesús Izaguirre, Eric Darve. Analysis of the accelerated weighted ensemble methodology. Conference Publications, 2013, 2013 (special) : 171-181. doi: 10.3934/proc.2013.2013.171
References:
[1]

D. Bhatt and D. M. Zuckerman, Heterogeneous path ensembles for conformational transitions in semiatomistic models of adenylate kinase,, Journal of Chemical Theory and Computation, 6 (2010), 3527.   Google Scholar

[2]

D. Bhatt, B. W. Zhang, and D. M. Zuckerman, Steady-state simulations using weighted ensemble path sampling,, The Journal of Chemical Physics, 133 (2010).   Google Scholar

[3]

D. Bhatt and D. M. Zuckerman, Symmetry of forward and reverse path populations,, arXiv:1002.2402 (physics.comp-ph), (2010).   Google Scholar

[4]

E. Darve and E. Ryu, Computing reaction rates in bio-molecular systems using discrete macro-states,, Innovations in Biomolecular Modeling and Simulations, 1 (2012), 138.   Google Scholar

[5]

A. Dickson and A. R. Dinner., Enhanced Sampling of Nonequilibrium Steady States,, Annu. Rev. Phys. Chem., 61 (2010), 441.   Google Scholar

[6]

A. Dickson, A. Warmflash and A. R. Dinner., Separating forward and backward pathways in nonequilibrium umbrella sampling,, J. Chem. Phys., 131 (2009).   Google Scholar

[7]

A. Dickson, A. Warmflash and A. R. Dinner., Nonequilibrium umbrella sampling in spaces of many order parameters,, J. Chem. Phys., 130 (2009).   Google Scholar

[8]

W. E and E. Vanden-Eijnden, Towards a Theory of Transition Paths,, Journal of Statistical Physics, 123 (2006), 503.   Google Scholar

[9]

G. A. Huber and S. Kim, Weighted ensemble Brownian dynamics simulations for protein association reactions,, Biophysical journal, 70 (1996).   Google Scholar

[10]

P. Metzner, C. Schütte and E. Vanden-Eijnden, Illustration of transition path theory on a collection of simple examples,, The Journal of Chemical Physics, 125 (2006).   Google Scholar

[11]

V. S. Pande, K. Beauchamp and G. R. Bowman, Everything you wanted to know about Markov state models but were afraid to ask,, Methods, 52 (2010).   Google Scholar

[12]

N. Singhal and V. S. Pande, Error analysis and efficient sampling in Markovian state models for molecular dynamics,, The Journal of Chemical Physics, 123 (2005).   Google Scholar

[13]

W. C. Swope, J. W. Pitera and F. Suits, Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 1. Theory,, The Journal of Physical Chemistry B, 108 (2004).   Google Scholar

[14]

E. Vanden-Eijnden, M. Venturoli, G. Ciccotti and R. Elber, On the assumptions underlying milestoning,, The Journal of Chemical Physics, 129 (2008).   Google Scholar

[15]

E. Vanden-Eijnden and M. Venturoli., Exact rate calculations by trajectory parallelization and tilting,, J. Chem. Phys., 131 (2009).   Google Scholar

[16]

A. Warmflash, P. Bhimalapuram and A. R. Dinner., Umbrella sampling for nonequilibrium processes,, J. Chem. Phys., 127 (2007).   Google Scholar

[17]

B. W Zhang, D. Jasnow and D. M. Zuckerman, "Weighted Ensemble Path Sampling for Multiple Reaction Channels,", arXiv.org, (2009).   Google Scholar

[18]

B. W. Zhang, D. Jasnow and D. M. Zuckerman, The "weighted ensemble" path sampling method is statistically exact for a broad class of stochastic processes and binning procedures,, The Journal of Chemical Physics, 132 (2010).   Google Scholar

show all references

References:
[1]

D. Bhatt and D. M. Zuckerman, Heterogeneous path ensembles for conformational transitions in semiatomistic models of adenylate kinase,, Journal of Chemical Theory and Computation, 6 (2010), 3527.   Google Scholar

[2]

D. Bhatt, B. W. Zhang, and D. M. Zuckerman, Steady-state simulations using weighted ensemble path sampling,, The Journal of Chemical Physics, 133 (2010).   Google Scholar

[3]

D. Bhatt and D. M. Zuckerman, Symmetry of forward and reverse path populations,, arXiv:1002.2402 (physics.comp-ph), (2010).   Google Scholar

[4]

E. Darve and E. Ryu, Computing reaction rates in bio-molecular systems using discrete macro-states,, Innovations in Biomolecular Modeling and Simulations, 1 (2012), 138.   Google Scholar

[5]

A. Dickson and A. R. Dinner., Enhanced Sampling of Nonequilibrium Steady States,, Annu. Rev. Phys. Chem., 61 (2010), 441.   Google Scholar

[6]

A. Dickson, A. Warmflash and A. R. Dinner., Separating forward and backward pathways in nonequilibrium umbrella sampling,, J. Chem. Phys., 131 (2009).   Google Scholar

[7]

A. Dickson, A. Warmflash and A. R. Dinner., Nonequilibrium umbrella sampling in spaces of many order parameters,, J. Chem. Phys., 130 (2009).   Google Scholar

[8]

W. E and E. Vanden-Eijnden, Towards a Theory of Transition Paths,, Journal of Statistical Physics, 123 (2006), 503.   Google Scholar

[9]

G. A. Huber and S. Kim, Weighted ensemble Brownian dynamics simulations for protein association reactions,, Biophysical journal, 70 (1996).   Google Scholar

[10]

P. Metzner, C. Schütte and E. Vanden-Eijnden, Illustration of transition path theory on a collection of simple examples,, The Journal of Chemical Physics, 125 (2006).   Google Scholar

[11]

V. S. Pande, K. Beauchamp and G. R. Bowman, Everything you wanted to know about Markov state models but were afraid to ask,, Methods, 52 (2010).   Google Scholar

[12]

N. Singhal and V. S. Pande, Error analysis and efficient sampling in Markovian state models for molecular dynamics,, The Journal of Chemical Physics, 123 (2005).   Google Scholar

[13]

W. C. Swope, J. W. Pitera and F. Suits, Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 1. Theory,, The Journal of Physical Chemistry B, 108 (2004).   Google Scholar

[14]

E. Vanden-Eijnden, M. Venturoli, G. Ciccotti and R. Elber, On the assumptions underlying milestoning,, The Journal of Chemical Physics, 129 (2008).   Google Scholar

[15]

E. Vanden-Eijnden and M. Venturoli., Exact rate calculations by trajectory parallelization and tilting,, J. Chem. Phys., 131 (2009).   Google Scholar

[16]

A. Warmflash, P. Bhimalapuram and A. R. Dinner., Umbrella sampling for nonequilibrium processes,, J. Chem. Phys., 127 (2007).   Google Scholar

[17]

B. W Zhang, D. Jasnow and D. M. Zuckerman, "Weighted Ensemble Path Sampling for Multiple Reaction Channels,", arXiv.org, (2009).   Google Scholar

[18]

B. W. Zhang, D. Jasnow and D. M. Zuckerman, The "weighted ensemble" path sampling method is statistically exact for a broad class of stochastic processes and binning procedures,, The Journal of Chemical Physics, 132 (2010).   Google Scholar

[1]

Huiying Fan, Tao Ma. Parabolic equations involving Laguerre operators and weighted mixed-norm estimates. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5487-5508. doi: 10.3934/cpaa.2020249

[2]

Håkon Hoel, Gaukhar Shaimerdenova, Raúl Tempone. Multilevel Ensemble Kalman Filtering based on a sample average of independent EnKF estimators. Foundations of Data Science, 2020  doi: 10.3934/fods.2020017

[3]

Zhiyan Ding, Qin Li, Jianfeng Lu. Ensemble Kalman Inversion for nonlinear problems: Weights, consistency, and variance bounds. Foundations of Data Science, 2020  doi: 10.3934/fods.2020018

[4]

Anna Abbatiello, Eduard Feireisl, Antoní Novotný. Generalized solutions to models of compressible viscous fluids. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 1-28. doi: 10.3934/dcds.2020345

[5]

Soniya Singh, Sumit Arora, Manil T. Mohan, Jaydev Dabas. Approximate controllability of second order impulsive systems with state-dependent delay in Banach spaces. Evolution Equations & Control Theory, 2020  doi: 10.3934/eect.2020103

[6]

Wei Feng, Michael Freeze, Xin Lu. On competition models under allee effect: Asymptotic behavior and traveling waves. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5609-5626. doi: 10.3934/cpaa.2020256

[7]

Anna Canale, Francesco Pappalardo, Ciro Tarantino. Weighted multipolar Hardy inequalities and evolution problems with Kolmogorov operators perturbed by singular potentials. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2020274

[8]

Weiwei Liu, Jinliang Wang, Yuming Chen. Threshold dynamics of a delayed nonlocal reaction-diffusion cholera model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020316

[9]

Abdelghafour Atlas, Mostafa Bendahmane, Fahd Karami, Driss Meskine, Omar Oubbih. A nonlinear fractional reaction-diffusion system applied to image denoising and decomposition. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020321

[10]

Serena Dipierro, Benedetta Pellacci, Enrico Valdinoci, Gianmaria Verzini. Time-fractional equations with reaction terms: Fundamental solutions and asymptotics. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 257-275. doi: 10.3934/dcds.2020137

[11]

Yongge Tian, Pengyang Xie. Simultaneous optimal predictions under two seemingly unrelated linear random-effects models. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2020168

[12]

Annegret Glitzky, Matthias Liero, Grigor Nika. Dimension reduction of thermistor models for large-area organic light-emitting diodes. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020460

[13]

Leilei Wei, Yinnian He. A fully discrete local discontinuous Galerkin method with the generalized numerical flux to solve the tempered fractional reaction-diffusion equation. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020319

[14]

H. M. Srivastava, H. I. Abdel-Gawad, Khaled Mohammed Saad. Oscillatory states and patterns formation in a two-cell cubic autocatalytic reaction-diffusion model subjected to the Dirichlet conditions. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020433

[15]

Lin Shi, Xuemin Wang, Dingshi Li. Limiting behavior of non-autonomous stochastic reaction-diffusion equations with colored noise on unbounded thin domains. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5367-5386. doi: 10.3934/cpaa.2020242

 Impact Factor: 

Metrics

  • PDF downloads (51)
  • HTML views (0)
  • Cited by (0)

[Back to Top]