QM/MM Simulations with the Gaussian Electrostatic Model: A Density-based Polarizable Potential.

Area of Science:

• Computational Chemistry
• Molecular Modeling
• Quantum Mechanics

Background:

• Advanced polarizable potentials enhance QM/MM simulation accuracy.
• The Gaussian Electrostatic Model (GEM) provides accurate QM/MM environments.
• Existing QM/MM methods have limitations in calculating intermolecular interactions.

Purpose of the Study:

• Introduce a new QM/GEM implementation for QM/MM simulations.
• Include all major intermolecular interaction energy components (except charge transfer).
• Assess the accuracy of the new QM/GEM implementation.

Main Methods:

• Developed a density-based potential, the Gaussian Electrostatic Model (GEM).
• Implemented QM/GEM to calculate Coulomb, exchange-repulsion, polarization, and dispersion energies.
• Tested the method using water dimer interactions and compared with AMOEBA potential and SAPT2+3/aug-cc-pVTZ calculations.

Main Results:

• The QM/GEM implementation accurately calculates intermolecular interaction energies.
• Achieved an average root-mean-square error <0.15 kcal/mol for all interaction energy components.
• Demonstrated comparable or improved accuracy over the classical AMOEBA potential.

Conclusions:

• The new QM/GEM implementation provides a highly accurate MM environment for QM/MM calculations.
• This method significantly improves the calculation of total intermolecular interaction energy.
• GEM represents a promising advancement for accurate molecular simulations.