Toward Hamiltonian Adaptive QM/MM: Accurate Solvent Structures Using Many-Body Potentials.

Area of Science:

• Computational Chemistry
• Multiscale Modeling
• Physical Chemistry

Background:

• Adaptive quantum mechanical (QM)/molecular mechanical (MM) methods are crucial for simulating chemistry in solution.
• Current adaptive QM/MM methods face challenges with structural distortions at the QM/MM boundary.
• Accurate modeling of solvent effects is essential for predicting molecular properties.

Purpose of the Study:

• To develop a novel Hamiltonian adaptive multiscale scheme for QM/MM simulations.
• To address and resolve the long-standing issue of structural distortions at the QM/MM interface.
• To enable accurate microcanonical adaptive QM/MM simulations for vibrational spectra and dynamical properties.

Main Methods:

• Implementation of a Hamiltonian adaptive multiscale scheme using many-body potentials.
• General adaptive expressions compatible with standard QM/MM embedding schemes.
• Application to a molecular system with two different molecular mechanics (MM) potentials and particle mesh Ewald for long-range interactions.

Main Results:

• The proposed Hamiltonian approach successfully eliminates structural distortions at the QM/MM boundary.
• Demonstrated energy conservation across the entire simulated system.
• Achieved accurate representation of solvent structure, even with different MM potentials.

Conclusions:

• The novel Hamiltonian adaptive multiscale scheme provides a robust solution for QM/MM simulations in solution.
• This method overcomes previous limitations, enabling reliable calculations of dynamical and spectral properties.
• The approach ensures both energy conservation and correct structural representation, advancing multiscale modeling capabilities.