Hybrid QM/MM Molecular Dynamics with AMOEBA Polarizable Embedding

Area of Science:

• Computational Chemistry
• Molecular Dynamics
• Quantum Mechanics/Molecular Mechanics (QM/MM)

Background:

• Accurate molecular dynamics (MD) simulations require sophisticated methods to model complex chemical systems.
• Hybrid quantum mechanics/molecular mechanics (QM/MM) approaches are essential for studying systems where electronic effects are crucial.
• Polarizable force fields, like AMOEBA, offer improved accuracy but pose computational challenges in QM/MM MD.

Purpose of the Study:

• To implement and evaluate a Born-Oppenheimer (BO) hybrid QM/MM MD strategy.
• To couple density functional theory (DFT) with the polarizable AMOEBA force field for enhanced simulation accuracy.
• To assess the stability and energy conservation of different BO formalisms in QM/MM MD.

Main Methods:

• Implementation of a BO hybrid QM/MM MD strategy linking Gaussian (for DFT) and Tinker (for AMOEBA) programs.
• Utilized a variational formalism for self-consistent relaxation of AMOEBA induced dipoles and DFT electron density.
• Compared the time-reversible BO (TR-BO) and extended BO Lagrangian (XL-BO) approaches for MD propagation.

Main Results:

• The extended BO Lagrangian (XL-BO) approach demonstrated stable, energy-conserving trajectories in QM/MM MD simulations.
• The coupling scheme allowed for full self-consistent relaxation of both electronic density and polarizable force field parameters.
• Computational efficiency was primarily limited by the DFT self-consistent field (SCF) cycles.

Conclusions:

• The XL-BO approach provides a robust and stable method for QM/MM MD simulations employing polarizable force fields.
• This implementation opens new avenues for accurate simulations of complex chemical and biological systems.
• The developed strategy enhances the capabilities of hybrid QM/MM methods in computational molecular science.