Accelerating Molecular Dynamics Simulations with Population Annealing.

Population annealing is a powerful tool for large-scale Monte Carlo simulations. We adapt this method to molecular dynamics simulations and demonstrate its excellent accelerating effect by simulating the folding of a short peptide commonly used to gauge the performance of algorithms. The method is compared to the well established parallel tempering approach and is found to yield similar performance for the same computational resources. In contrast to other methods, however, population annealing scales to a nearly arbitrary number of parallel processors, and it is thus a unique tool that enables molecular dynamics to tap into the massively parallel computing power available in supercomputers that is so much needed for a range of difficult computational problems.

[1]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[2]  Yukito IBA,et al.  Population Monte Carlo algorithms , 2000, cond-mat/0008226.

[3]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[4]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[5]  Wolfhard Janke,et al.  GPU accelerated population annealing algorithm , 2017, Comput. Phys. Commun..

[6]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[7]  H. Katzgraber,et al.  Population annealing: Theory and application in spin glasses. , 2015, Physical review. E, Statistical, nonlinear, and soft matter physics.

[8]  J. Machta,et al.  Population annealing with weighted averages: a Monte Carlo method for rough free-energy landscapes. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  C. Geyer Markov Chain Monte Carlo Maximum Likelihood , 1991 .

[10]  Vijay S. Pande,et al.  OpenMM 7: Rapid development of high performance algorithms for molecular dynamics , 2016, bioRxiv.

[11]  Wolfhard Janke,et al.  Exploring first-order phase transitions with population annealing , 2017 .

[12]  E. LESTER SMITH,et al.  AND OTHERS , 2005 .

[13]  C. Lowe,et al.  An alternative approach to dissipative particle dynamics , 1999 .

[14]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules J. Am. Chem. Soc. 1995, 117, 5179−5197 , 1996 .

[15]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[16]  U. Hansmann Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[17]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[18]  Jonathan Machta,et al.  Population annealing simulations of a binary hard-sphere mixture. , 2017, Physical review. E.

[19]  Koji Hukushima,et al.  Population Annealing and Its Application to a Spin Glass , 2003 .

[20]  Patrice Koehl,et al.  Faculty Opinions recommendation of OpenMM 7: Rapid development of high performance algorithms for molecular dynamics. , 2018, Faculty Opinions – Post-Publication Peer Review of the Biomedical Literature.

[21]  W. Janke,et al.  Dissipative Dynamics of a Single Polymer in Solution: A Lowe-Andersen Approach , 2018, Journal of Physics: Conference Series.

[22]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[23]  P. Grassberger Go with the Winners: a General Monte Carlo Strategy , 2002, cond-mat/0201313.

[24]  KumarShankar,et al.  The weighted histogram analysis method for free-energy calculations on biomolecules. I , 1992 .

[25]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[26]  W. Janke,et al.  Population annealing molecular dynamics with adaptive temperature steps , 2019, Journal of Physics: Conference Series.

[27]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical review letters.