The adaptive buffered force QM/MM method in the CP2K and AMBER software packages
Abstract
We present the implementation and validation of the adaptive buffered force (AdBF) quantum-mechanics/molecular-mechanics (QM/MM) method in two popular packages, CP2K and AMBER. The implementations build on the existing QM/MM functionality in each code, extending it to allow for redefinition of the QM and MM regions during the simulation and reducing QM-MM interface errors by discarding forces near the boundary according to the buffered force-mixing approach. New adaptive thermostats, needed by force-mixing methods, are also implemented. Different variants of the method are benchmarked by simulating the structure of bulk water, water autoprotolysis in the presence of zinc and dimethyl-phosphate hydrolysis using various semiempirical Hamiltonians and density functional theory as the QM model. It is shown that with suitable parameters, based on force convergence tests, the AdBF QM/MM scheme can provide an accurate approximation of the structure in the dynamical QM region matching the corresponding fully QM simulations, as well as reproducing the correct energetics in all cases. Adaptive unbuffered force-mixing and adaptive conventional QM/MM methods also provide reasonable results for some systems, but are more likely to suffer from instabilities and inaccuracies.
- Authors:
-
- Engineering Department University of Cambridge Cambridge CB2 1PZ United Kingdom
- Institute for Condensed Matter and Complex Systems School of Physics and Astronomy University of Edinburgh Edinburgh EH9 3JZ United Kingdom
- San Diego Supercomputer Center University of California San Diego La Jolla California 92093
- Mathematical and Computational Sciences Department IBM Research–Zurich Säumerstrasse 4 8803 Rüschlikon Switzerland
- San Diego Supercomputer Center University of California San Diego La Jolla California 92093, Department of Chemistry and Biochemistry University of California San Diego La Jolla California 92093
- The Maxwell Institute and School of Mathematics, University of Edinburgh Edinburgh EH9 3JZ United Kingdom
- Center for Computational Material Science Naval Research Laboratory Washington DC 20375
- Publication Date:
- Research Org.:
- Univ. of California, San Diego, CA (United States)
- Sponsoring Org.:
- USDOE; National Institutes of Health (NIH); National Science Foundation (NSF)
- OSTI Identifier:
- 1345677
- Alternate Identifier(s):
- OSTI ID: 1343098; OSTI ID: 1345678
- Grant/Contract Number:
- DE‐AC36‐99GO‐10337; AC36-99GO10337; OCI-1148358; ACI-1053575; EP/G036136/1; EP/J01298X/1; GM100934
- Resource Type:
- Published Article
- Journal Name:
- Journal of Computational Chemistry
- Additional Journal Information:
- Journal Name: Journal of Computational Chemistry Journal Volume: 36 Journal Issue: 9; Journal ID: ISSN 0192-8651
- Publisher:
- Wiley Blackwell (John Wiley & Sons)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; quantum-mechanics/molecular-mechanics; adaptive quantum-mechanics/molecular-mechanics; force-mixing; multiscale
Citation Formats
Mones, Letif, Jones, Andrew, Götz, Andreas W., Laino, Teodoro, Walker, Ross C., Leimkuhler, Ben, Csányi, Gábor, and Bernstein, Noam. The adaptive buffered force QM/MM method in the CP2K and AMBER software packages. United States: N. p., 2015.
Web. doi:10.1002/jcc.23839.
Mones, Letif, Jones, Andrew, Götz, Andreas W., Laino, Teodoro, Walker, Ross C., Leimkuhler, Ben, Csányi, Gábor, & Bernstein, Noam. The adaptive buffered force QM/MM method in the CP2K and AMBER software packages. United States. https://doi.org/10.1002/jcc.23839
Mones, Letif, Jones, Andrew, Götz, Andreas W., Laino, Teodoro, Walker, Ross C., Leimkuhler, Ben, Csányi, Gábor, and Bernstein, Noam. Tue .
"The adaptive buffered force QM/MM method in the CP2K and AMBER software packages". United States. https://doi.org/10.1002/jcc.23839.
@article{osti_1345677,
title = {The adaptive buffered force QM/MM method in the CP2K and AMBER software packages},
author = {Mones, Letif and Jones, Andrew and Götz, Andreas W. and Laino, Teodoro and Walker, Ross C. and Leimkuhler, Ben and Csányi, Gábor and Bernstein, Noam},
abstractNote = {We present the implementation and validation of the adaptive buffered force (AdBF) quantum-mechanics/molecular-mechanics (QM/MM) method in two popular packages, CP2K and AMBER. The implementations build on the existing QM/MM functionality in each code, extending it to allow for redefinition of the QM and MM regions during the simulation and reducing QM-MM interface errors by discarding forces near the boundary according to the buffered force-mixing approach. New adaptive thermostats, needed by force-mixing methods, are also implemented. Different variants of the method are benchmarked by simulating the structure of bulk water, water autoprotolysis in the presence of zinc and dimethyl-phosphate hydrolysis using various semiempirical Hamiltonians and density functional theory as the QM model. It is shown that with suitable parameters, based on force convergence tests, the AdBF QM/MM scheme can provide an accurate approximation of the structure in the dynamical QM region matching the corresponding fully QM simulations, as well as reproducing the correct energetics in all cases. Adaptive unbuffered force-mixing and adaptive conventional QM/MM methods also provide reasonable results for some systems, but are more likely to suffer from instabilities and inaccuracies.},
doi = {10.1002/jcc.23839},
journal = {Journal of Computational Chemistry},
number = 9,
volume = 36,
place = {United States},
year = {Tue Feb 03 00:00:00 EST 2015},
month = {Tue Feb 03 00:00:00 EST 2015}
}
https://doi.org/10.1002/jcc.23839
Web of Science
Works referenced in this record:
Third-Generation Hydrogen-Bonding Corrections for Semiempirical QM Methods and Force Fields
journal, November 2010
- Korth, Martin
- Journal of Chemical Theory and Computation, Vol. 6, Issue 12
An Efficient Real Space Multigrid QM/MM Electrostatic Coupling
journal, November 2005
- Laino, Teodoro; Mohamed, Fawzi; Laio, Alessandro
- Journal of Chemical Theory and Computation, Vol. 1, Issue 6
Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme
journal, May 1976
- Warshel, A.; Levitt, M.
- Journal of Molecular Biology, Vol. 103, Issue 2
An Efficient Linear-Scaling Electrostatic Coupling for Treating Periodic Boundary Conditions in QM/MM Simulations
journal, August 2006
- Laino, Teodoro; Mohamed, Fawzi; Laio, Alessandro
- Journal of Chemical Theory and Computation, Vol. 2, Issue 5
Density-functional exchange-energy approximation with correct asymptotic behavior
journal, September 1988
- Becke, A. D.
- Physical Review A, Vol. 38, Issue 6
A combinedab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems: Applications to the CH3Cl + Cl? exchange reaction and gas phase protonation of polyethers
journal, December 1986
- Singh, U. Chandra; Kollman, Peter A.
- Journal of Computational Chemistry, Vol. 7, Issue 6
Hybrid Models for Combined Quantum Mechanical and Molecular Mechanical Approaches
journal, January 1996
- Bakowies, Dirk; Thiel, Walter
- The Journal of Physical Chemistry, Vol. 100, Issue 25
Results obtained with the correlation energy density functionals of becke and Lee, Yang and Parr
journal, May 1989
- Miehlich, Burkhard; Savin, Andreas; Stoll, Hermann
- Chemical Physics Letters, Vol. 157, Issue 3
Application of RESP charges to calculate conformational energies, hydrogen bond energies, and free energies of solvation
journal, October 1993
- Cornell, Wendy D.; Cieplak, Piotr; Bayly, Christopher I.
- Journal of the American Chemical Society, Vol. 115, Issue 21
Adaptive resolution molecular-dynamics simulation: Changing the degrees of freedom on the fly
journal, December 2005
- Praprotnik, Matej; Delle Site, Luigi; Kremer, Kurt
- The Journal of Chemical Physics, Vol. 123, Issue 22
Ab initio quantum mechanical charge field (QMCF) molecular dynamics: a QM/MM – MD procedure for accurate simulations of ions and complexes
journal, December 2005
- Rode, Bernd M.; Hofer, Thomas S.; Randolf, Bernhard R.
- Theoretical Chemistry Accounts, Vol. 115, Issue 2-3
Adaptive resolution scheme for efficient hybrid atomistic-mesoscale molecular dynamics simulations of dense liquids
journal, June 2006
- Praprotnik, Matej; Delle Site, Luigi; Kremer, Kurt
- Physical Review E, Vol. 73, Issue 6
Density, structure, and dynamics of water: The effect of van der Waals interactions
journal, January 2011
- Wang, Jue; Román-Pérez, G.; Soler, Jose M.
- The Journal of Chemical Physics, Vol. 134, Issue 2
ONIOM-XS: an extension of the ONIOM method for molecular simulation in condensed phase
journal, April 2002
- Kerdcharoen, Teerakiat; Morokuma, Keiji
- Chemical Physics Letters, Vol. 355, Issue 3-4
QM/MM Methods for Biomolecular Systems
journal, January 2009
- Senn, Hans Martin; Thiel, Walter
- Angewandte Chemie International Edition, Vol. 48, Issue 7
All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins †
journal, April 1998
- MacKerell, A. D.; Bashford, D.; Bellott, M.
- The Journal of Physical Chemistry B, Vol. 102, Issue 18
Why nature really chose phosphate
journal, January 2013
- Kamerlin, Shina C. L.; Sharma, Pankaz K.; Prasad, Ram B.
- Quarterly Reviews of Biophysics, Vol. 46, Issue 1
Coordination numbers as reaction coordinates in constrained molecular dynamics
journal, January 1998
- Sprik, Michiel
- Faraday Discussions, Vol. 110
Comparison of multiple Amber force fields and development of improved protein backbone parameters
journal, November 2006
- Hornak, Viktor; Abel, Robert; Okur, Asim
- Proteins: Structure, Function, and Bioinformatics, Vol. 65, Issue 3
Specific Reaction Parametrization of the AM1/d Hamiltonian for Phosphoryl Transfer Reactions: H, O, and P Atoms
journal, January 2007
- Nam, Kwangho; Cui, Qiang; Gao, Jiali
- Journal of Chemical Theory and Computation, Vol. 3, Issue 2
Extension of the MNDO formalism tod orbitals: Integral approximations and preliminary numerical results
journal, January 1992
- Thiel, Walter; Voityuk, Alexander A.
- Theoretica Chimica Acta, Vol. 81, Issue 6
Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model
journal, June 1985
- Dewar, Michael J. S.; Zoebisch, Eve G.; Healy, Eamonn F.
- Journal of the American Chemical Society, Vol. 107, Issue 13
Efficient molecular dynamics and hybrid Monte Carlo algorithms for path integrals
journal, August 1993
- Tuckerman, Mark E.; Berne, Bruce J.; Martyna, Glenn J.
- The Journal of Chemical Physics, Vol. 99, Issue 4
Developing ab initio quality force fields from condensed phase quantum-mechanics/molecular-mechanics calculations through the adaptive force matching method
journal, August 2008
- Akin-Ojo, Omololu; Song, Yang; Wang, Feng
- The Journal of Chemical Physics, Vol. 129, Issue 6
Bridging the gap between thermodynamic integration and umbrella sampling provides a novel analysis method: “Umbrella integration”
journal, October 2005
- Kästner, Johannes; Thiel, Walter
- The Journal of Chemical Physics, Vol. 123, Issue 14
A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: the RESP model
journal, October 1993
- Bayly, Christopher I.; Cieplak, Piotr; Cornell, Wendy
- The Journal of Physical Chemistry, Vol. 97, Issue 40
Dynamic QM/MM: A Hybrid Approach to Simulating Gas-Liquid Interactions
book, January 2011
- Yockel, Scott; Schatz, George C.
- Multiscale Molecular Methods in Applied Chemistry
Multiscale Modeling of Chemistry in Water: Are We There Yet?
journal, November 2013
- Bulo, Rosa E.; Michel, Carine; Fleurat-Lessard, Paul
- Journal of Chemical Theory and Computation, Vol. 9, Issue 12
A Generalized Hybrid Orbital (GHO) Method for the Treatment of Boundary Atoms in Combined QM/MM Calculations
journal, June 1998
- Gao, Jiali; Amara, Patricia; Alhambra, Cristobal
- The Journal of Physical Chemistry A, Vol. 102, Issue 24
RM1: A reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I
journal, July 2006
- Rocha, Gerd B.; Freire, Ricardo O.; Simas, Alfredo M.
- Journal of Computational Chemistry, Vol. 27, Issue 10
Energetics and Dynamics of Enzymatic Reactions
journal, August 2001
- Villà, Jordi; Warshel, Arieh
- The Journal of Physical Chemistry B, Vol. 105, Issue 33
The Solvation Structure of Na + and K + in Liquid Water Determined from High Level ab Initio Molecular Dynamics Simulations
journal, March 2012
- Rowley, Christopher N.; Roux, Benoı̂t
- Journal of Chemical Theory and Computation, Vol. 8, Issue 10
The quest for the best nonpolarizable water model from the adaptive force matching method
journal, August 2010
- Akin-Ojo, Omololu; Wang, Feng
- Journal of Computational Chemistry, Vol. 32, Issue 3
A pseudobond approach to combining quantum mechanical and molecular mechanical methods
journal, January 1999
- Zhang, Yingkai; Lee, Tai-Sung; Yang, Weitao
- The Journal of Chemical Physics, Vol. 110, Issue 1
An overview of the Amber biomolecular simulation package: Amber biomolecular simulation package
journal, September 2012
- Salomon-Ferrer, Romelia; Case, David A.; Walker, Ross C.
- Wiley Interdisciplinary Reviews: Computational Molecular Science, Vol. 3, Issue 2
Optimization of parameters for semiempirical methods I. Method
journal, March 1989
- Stewart, James J. P.
- Journal of Computational Chemistry, Vol. 10, Issue 2
On the Convergence of QM/MM Energies
journal, February 2011
- Hu, LiHong; Söderhjelm, Pär; Ryde, Ulf
- Journal of Chemical Theory and Computation, Vol. 7, Issue 3
A QM/MM simulation method applied to the solution of Li+ in liquid ammonia
journal, November 1996
- Kerdcharoen, Teerakiat; Liedl, Klaus R.; Rode, Bernd M.
- Chemical Physics, Vol. 211, Issue 1-3
Canonical sampling through velocity rescaling
journal, January 2007
- Bussi, Giovanni; Donadio, Davide; Parrinello, Michele
- The Journal of Chemical Physics, Vol. 126, Issue 1
Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density
journal, January 1988
- Lee, Chengteh; Yang, Weitao; Parr, Robert G.
- Physical Review B, Vol. 37, Issue 2
Application of Adaptive QM/MM Methods to Molecular Dynamics Simulations of Aqueous Systems
journal, July 2012
- Park, Kyoyeon; Götz, Andreas W.; Walker, Ross C.
- Journal of Chemical Theory and Computation, Vol. 8, Issue 8
Particle mesh Ewald: An N ⋅log( N ) method for Ewald sums in large systems
journal, June 1993
- Darden, Tom; York, Darrin; Pedersen, Lee
- The Journal of Chemical Physics, Vol. 98, Issue 12
Toward a Practical Method for Adaptive QM/MM Simulations
journal, July 2009
- Bulo, Rosa E.; Ensing, Bernd; Sikkema, Jetze
- Journal of Chemical Theory and Computation, Vol. 5, Issue 9
Generalized Langevin equation approach for atom/solid-surface scattering: General formulation for classical scattering off harmonic solids
journal, January 1976
- Adelman, S. A.
- The Journal of Chemical Physics, Vol. 64, Issue 6
A Gentle Stochastic Thermostat for Molecular Dynamics
journal, April 2009
- Leimkuhler, Ben; Noorizadeh, Emad; Theil, Florian
- Journal of Statistical Physics, Vol. 135, Issue 2
Force field design for metalloproteins
journal, October 1991
- Hoops, Stephen C.; Anderson, Kenneth W.; Merz, Kenneth M.
- Journal of the American Chemical Society, Vol. 113, Issue 22
Hybrid atomistic simulation methods for materials systems
journal, January 2009
- Bernstein, N.; Kermode, J. R.; Csányi, G.
- Reports on Progress in Physics, Vol. 72, Issue 2
Basic ideas for the correction of semiempirical methods describing H-bonded systems
journal, November 2000
- Bernal-Uruchurtu, M. I.; Ruiz-López, M. F.
- Chemical Physics Letters, Vol. 330, Issue 1-2
Adaptive stochastic methods for sampling driven molecular systems
journal, August 2011
- Jones, Andrew; Leimkuhler, Ben
- The Journal of Chemical Physics, Vol. 135, Issue 8
Separable dual-space Gaussian pseudopotentials
journal, July 1996
- Goedecker, S.; Teter, M.; Hutter, J.
- Physical Review B, Vol. 54, Issue 3
Semiempirical Self-Consistent Polarization Description of Bulk Water, the Liquid−Vapor Interface, and Cubic Ice
journal, June 2011
- Murdachaew, Garold; Mundy, Christopher J.; Schenter, Gregory K.
- The Journal of Physical Chemistry A, Vol. 115, Issue 23
Topologically Invariant Reaction Coordinates for Simulating Multistate Chemical Reactions
journal, December 2012
- Mones, Letif; Csányi, Gábor
- The Journal of Physical Chemistry B, Vol. 116, Issue 51
Energy Conservation in Adaptive Hybrid Atomistic/Coarse-Grain Molecular Dynamics
journal, March 2007
- Ensing, Bernd; Nielsen, Steven O.; Moore, Preston B.
- Journal of Chemical Theory and Computation, Vol. 3, Issue 3
QM/MM simulation of liquid water with an adaptive quantum region
journal, January 2012
- Bernstein, Noam; Várnai, Csilla; Solt, Iván
- Phys. Chem. Chem. Phys., Vol. 14, Issue 2
The implementation of a fast and accurate QM/MM potential method in Amber
journal, January 2008
- Walker, Ross C.; Crowley, Michael F.; Case, David A.
- Journal of Computational Chemistry, Vol. 29, Issue 7
Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements
journal, September 2007
- Stewart, James J. P.
- Journal of Molecular Modeling, Vol. 13, Issue 12
Semiempirical GGA-type density functional constructed with a long-range dispersion correction
journal, January 2006
- Grimme, Stefan
- Journal of Computational Chemistry, Vol. 27, Issue 15, p. 1787-1799
Ground states of molecules. 38. The MNDO method. Approximations and parameters
journal, June 1977
- Dewar, Michael J. S.; Thiel, Walter
- Journal of the American Chemical Society, Vol. 99, Issue 15
An Efficient Linear-Scaling Ewald Method for Long-Range Electrostatic Interactions in Combined QM/MM Calculations
journal, October 2004
- Nam, Kwangho; Gao, Jiali; York, Darrin M.
- Journal of Chemical Theory and Computation, Vol. 1, Issue 1
“Learn on the Fly”: A Hybrid Classical and Quantum-Mechanical Molecular Dynamics Simulation
journal, October 2004
- Csányi, Gabor; Albaret, T.; Payne, M. C.
- Physical Review Letters, Vol. 93, Issue 17
An extensible interface for QM/MM molecular dynamics simulations with AMBER
journal, October 2013
- Götz, Andreas W.; Clark, Matthew A.; Walker, Ross C.
- Journal of Computational Chemistry, Vol. 35, Issue 2
Adaptive Partitioning in Combined Quantum Mechanical and Molecular Mechanical Calculations of Potential Energy Functions for Multiscale Simulations
journal, March 2007
- Heyden, Andreas; Lin, Hai; Truhlar, Donald G.
- The Journal of Physical Chemistry B, Vol. 111, Issue 9
Hybrid Quantum Mechanics/Molecular Mechanics-Based Molecular Dynamics Simulation of Acid-Catalyzed Dehydration of Polyols in Liquid Water
journal, August 2011
- Caratzoulas, Stavros; Courtney, Timothy; Vlachos, Dionisios G.
- The Journal of Physical Chemistry A, Vol. 115, Issue 32
Notes on “Ewald summation of electrostatic multipole interactions up to quadrupolar level” [J. Chem. Phys. 119, 7471 (2003)]
journal, August 2008
- Laino, Teodoro; Hutter, Jürg
- The Journal of Chemical Physics, Vol. 129, Issue 7
Die Berechnung optischer und elektrostatischer Gitterpotentiale
journal, January 1921
- Ewald, P. P.
- Annalen der Physik, Vol. 369, Issue 3
Evaluating Boundary Dependent Errors in QM/MM Simulations
journal, April 2009
- Solt, Iván; Kulhánek, Petr; Simon, István
- The Journal of Physical Chemistry B, Vol. 113, Issue 17
Conservative Algorithm for an Adaptive Change of Resolution in Mixed Atomistic/Coarse-Grained Multiscale Simulations
journal, January 2008
- Heyden, Andreas; Truhlar, Donald G.
- Journal of Chemical Theory and Computation, Vol. 4, Issue 2
Tests of an Adaptive QM/MM Calculation on Free Energy Profiles of Chemical Reactions in Solution
journal, September 2013
- Várnai, Csilla; Bernstein, Noam; Mones, Letif
- The Journal of Physical Chemistry B, Vol. 117, Issue 40
Thermostats for “Slow” Configurational Modes
journal, July 2007
- Samoletov, Alex A.; Dettmann, Carl P.; Chaplain, Mark A. J.
- Journal of Statistical Physics, Vol. 128, Issue 6