A separable shadow Hamiltonian hybrid Monte Carlo method
J. Chem. Phys. 131, 174106 (2009); doi:10.1063/1.3253687
Published 3 November 2009
You are not logged in to this journal. Log in
Hybrid Monte Carlo (HMC) is a rigorous sampling method that uses molecular dynamics (MD) as a global Monte Carlo move. The acceptance rate of HMC decays exponentially with system size. The shadow hybrid Monte Carlo (SHMC) was previously introduced to reduce this performance degradation by sampling instead from the shadow Hamiltonian defined for MD when using a symplectic integrator. SHMC's performance is limited by the need to generate momenta for the MD step from a nonseparable shadow Hamiltonian. We introduce the separable shadow Hamiltonian hybrid Monte Carlo (S2HMC) method based on a formulation of the leapfrog/Verlet integrator that corresponds to a separable shadow Hamiltonian, which allows efficient generation of momenta. S2HMC gives the acceptance rate of a fourth order integrator at the cost of a second-order integrator. Through numerical experiments we show that S2HMC consistently gives a speedup greater than two over HMC for systems with more than 4000 atoms for the same variance. By comparison, SHMC gave a maximum speedup of only 1.6 over HMC. S2HMC has the additional advantage of not requiring any user parameters beyond those of HMC. S2HMC is available in the program PROTOMOL 2.1. A Python version, adequate for didactic purposes, is also in MDL (http://mdlab.sourceforge.net/s2hmc).
©2009 American Institute of Physics
| History: | Received 5 May 2009; accepted 2 October 2009; published 3 November 2009 |
| Permalink: |
http://link.aip.org/link/?JCPSA6/131/174106/1 |
REFERENCES (24)
For access to fully linked references, you need to log in.
For access to fully linked references, you need to Log in.
- S. Duane, A. D. Kennedy, B. J. Pendleton, and D. Roweth,
Phys. Lett. B 195, 216 (1987) . - A. Brass, B. J. Pendleton, Y. Chen, and B. Robson,
Biopolymers 33, 1307 (1993) . - B. Mehlig, D. W. Heermann, and B. M. Forrest, Phys. Rev. B 45, 679 (1992).
- M. Creutz, Phys. Rev. D 38, 1228 (1988).
- A. D. Kennedy and B. Pendleton,
Nucl. Phys. B, Proc. Suppl. 20, 118 (1991) . - M. López-Marcos, J. M. Sanz-Serna, and R. D. Skeel, Numerical Analysis 1995 (Longmans, Green, New York, 1996), pp. 107–122.
- E. Hairer, C. Lubich, and G. Wanner, Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations, Springer Series in Computational Mathematics Vol. 31 (Springer-Verlag, Berlin, 2002).
- J. A. Izaguirre and S. S. Hampton,
J. Comput. Phys. 200, 581 (2004) . - E. Akhmatskaya and S. Reich, New Algorithms for Macromolecular Simulation (Springer-Verlag, Berlin, 2006), Vol. 49, pp. 141–151.
- C. L. Wee, M. S. Sansom, S. Reich, and E. Akhmatskaya,
J. Phys. Chem. B 112, 5710 (2008) . - E. Akhmatskaya and S. Reich,
J. Comput. Phys. 227, 4934 (2008) . - E. Akhmatskaya, N. Bou-Rabee, and S. Reich,
J. Comput. Phys. 228, 2256 (2009) . - T. Matthey, T. Cickovski, S. S. Hampton, A. Ko, Q. Ma, M. Nyerges, T. Raeder, T. Slabach, and J. A. Izaguirre,
ACM Trans. Math. Softw. 30, 237 (2004) . - L. Verlet,
Phys. Rev. 159, 98 (1967) . - T. Cickovski, C. Sweet, and J. A. Izaguirre, Proceedings of the 40th Annual Simulation Symposium, 2007, pp. 256–266.
- R. D. Skeel and D. J. Hardy,
SIAM J. Sci. Comput. (USA) 23, 1172 (2001) . - R. D. Engle, R. D. Skeel, and M. Drees,
J. Comput. Phys. 206, 432 (2005) . - S. Blanes, S. Casas, and A. Murua,
SIAM (Soc. Ind. Appl. Math.) J. Numer. Anal. 42, 531 (2004) . - T. Littell, R. Skeel, and M. Zhang,
SIAM (Soc. Ind. Appl. Math.) J. Numer. Anal. 34, 1792 (1997) . - P. R. Bevington and D. K. Robinson, Data Reduction and Error Analysis for the Physical Sciences (McGraw-Hill, New York, 2003).
- M. E. Nelson and J. M. Bower,
Trends Neurosci. 13, 403 (1990) . - A. D. MacKerell, Jr., D. Bashford, M. Bellott, R. L. Dunbrack, Jr., J. Evanseck, M. J. Field, S. Fischer, J. Gao, H. Guo, S. Ha, D. Joseph, L. Kuchnir, K. Kuczera, F. T. K. Lau, C. Mattos, S. Michnick, T. Ngo, D. T. Nguyen, B. Prodhom, B. Roux, M. Schlenkrich, J. Smith, R. Stote, J. Straub, M. Watanabe, J. Wiorkiewicz-Kuczera, D. Yin, and M. Karplus, FASEB J. A143, 6 (1992).
- A. D. MacKerell, Jr., D. Bashford, M. Bellott, R. L. Dunbrack, Jr., J. Evanseck, M. J. Field, S. Fischer, J. Gao, H. Guo, S. Ha, D. Joseph, L. Kuchnir, K. Kuczera, F. T. K. Lau, C. Mattos, S. Michnick, T. Ngo, D. T. Nguyen, B. Prodhom, I. W. E. Reiher, B. Roux, M. Schlenkrich, J. Smith, R. Stote, J. Straub, M. Watanabe, J. Wiorkiewicz-Kuczera, D. Yin, and M. Karplus,
J. Phys. Chem. B 102, 3586 (1998) . - W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys. 79, 926 (1983).








