DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Representation of the QM Subsystem for Long-Range Electrostatic Interaction in Non-Periodic Ab Initio QM/MM Calculations

Journal Article · · Molecules

In QM/MM calculations, it is essential to handle electrostatic interactions between the QM and MM subsystems accurately and efficiently. To achieve maximal efficiency, it is convenient to adopt a hybrid scheme, where the QM electron density is used explicitly in the evaluation of short-range QM/MM electrostatic interactions, while a multipolar representation for the QM electron density is employed to account for the long-range QM/MM electrostatic interactions. In order to avoid energy discontinuity at the cutoffs, which separate the short- and long-range QM/MM electrostatic interactions, a switching function should be utilized to ensure a smooth potential energy surface. In this study, we benchmarked the accuracy of such hybrid embedding schemes for QM/MM electrostatic interactions using different multipolar representations, switching functions and cutoff distances. For test systems (neutral and anionic oxyluciferin in MM (aqueous and enzyme) environments), the best accuracy was acquired with a combination of QM electrostatic potential (ESP) charges and dipoles and two switching functions (long-range electrostatic corrections (LREC) and Switch) in the treatment of long-range QM/MM electrostatics. It allowed us to apply a 10Å distance cutoff and still obtain QM/MM electrostatics/polarization energies within 0.1 kcal/mol and time-dependent density functional theory (TDDFT)/MM vertical excitation energies within 10−3 eV from theoretical reference values.

Research Organization:
Q-Chem, Inc., Pleasanton, CA (United States)
Sponsoring Organization:
USDOE; USDOE Office of Science (SC)
Grant/Contract Number:
SC0011297
OSTI ID:
2372731
Journal Information:
Molecules, Journal Name: Molecules Journal Issue: 10 Vol. 23; ISSN MOLEFW; ISSN 1420-3049
Publisher:
MDPI AGCopyright Statement
Country of Publication:
Switzerland
Language:
English

References (47)

Scalable molecular dynamics with NAMD journal January 2005
The implementation of a fast and accurate QM/MM potential method in Amber journal January 2008
CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields journal January 2009
Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table : Molcas 8 journal November 2015
An approach to computing electrostatic charges for molecules journal April 1984
New spherical-cutoff methods for long-range forces in macromolecular simulation journal July 1994
QM/MM: what have we learned, where are we, and where do we go from here? journal July 2006
Long-range electrostatic corrections in multipolar/polarizable QM/MM simulations journal June 2016
The continuous fast multipole method journal November 1994
Approximate electrostatic interaction operator for QM/MM calculations journal April 2002
QM/MM studies of enzymes journal April 2007
Ambient-Potential Composite Ewald Method for ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulation journal May 2016
Periodic Boundary Conditions in QM/MM Calculations: Implementation and Tests journal July 2016
Single-Reference ab Initio Methods for the Calculation of Excited States of Large Molecules journal November 2005
Classical Electrostatics for Biomolecular Simulations journal August 2013
An Efficient Linear-Scaling Ewald Method for Long-Range Electrostatic Interactions in Combined QM/MM Calculations journal October 2004
Including Charge Penetration Effects in Molecular Modeling journal October 2010
Long-Range Electrostatic Effects in QM/MM Studies of Enzymatic Reactions: Application of the Solvated Macromolecule Boundary Potential journal December 2010
Optimization of the Additive CHARMM All-Atom Protein Force Field Targeting Improved Sampling of the Backbone ϕ, ψ and Side-Chain χ 1 and χ 2 Dihedral Angles journal August 2012
Tuned and Balanced Redistributed Charge Scheme for Combined Quantum Mechanical and Molecular Mechanical (QM/MM) Methods and Fragment Methods: Tuning Based on the CM5 Charge Model journal December 2012
Isotropic Periodic Sum Treatment of Long-Range Electrostatic Interactions in Combined Quantum Mechanical and Molecular Mechanical Calculations journal November 2013
Fitting Molecular Electrostatic Potentials from Quantum Mechanical Calculations journal March 2007
A General Boundary Potential for Hybrid QM/MM Simulations of Solvated Biomolecular Systems journal October 2009
Toward the Accurate Modeling of DNA: The Importance of Long-Range Electrostatics journal May 1995
Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices journal August 1998
All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins journal April 1998
Structural basis for the spectral difference in luciferase bioluminescence journal March 2006
NAMD goes quantum: an integrative suite for hybrid simulations journal March 2018
Generalized solvent boundary potential for computer simulations journal February 2001
A Hamiltonian electrostatic coupling scheme for hybrid Car–Parrinello molecular dynamics simulations journal April 2002
The Fourier transform Coulomb method: Efficient and accurate calculation of the Coulomb operator in a Gaussian basis journal November 2002
Reliable treatment of electrostatics in combined QM/MM simulation of macromolecules journal July 2005
Ewald mesh method for quantum mechanical calculations journal March 2012
Comparison of simple potential functions for simulating liquid water journal July 1983
Particle mesh Ewald: An N ⋅log( N ) method for Ewald sums in large systems journal June 1993
The effect of long‐range electrostatic interactions in simulations of macromolecular crystals: A comparison of the Ewald and truncated list methods journal November 1993
A smooth particle mesh Ewald method journal November 1995
Periodic boundary conditions for QM/MM calculations: Ewald summation for extended Gaussian basis sets journal December 2013
Erratum: “Periodic boundary conditions for QM/MM calculations: Ewald summation for extended Gaussian basis sets” [J. Chem. Phys. 139 , 244108 (2013)] journal February 2015
A new smoothing function to introduce long-range electrostatic effects in QM/MM calculations journal July 2015
Treating electrostatics with Wolf summation in combined quantum mechanical and molecular mechanical simulations journal November 2015
TINKTEP: A fully self-consistent, mutually polarizable QM/MM approach based on the AMOEBA force field journal September 2016
Atomic-level accuracy in simulations of large protein crystals. journal August 1994
Damping functions in the effective fragment potential method journal April 2009
Advances in molecular quantum chemistry contained in the Q-Chem 4 program package journal September 2014
A priori evaluation of aqueous polarization effects through Monte Carlo QM-MM simulations journal October 1992
Progress in Time-Dependent Density-Functional Theory journal May 2012