Abstract:A powerful tool to study the mechanism of reactions in solutions or enzymes is to perform the ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations. However, the computational cost is too high due to the explicit electronic structure calculations at every time step of the simulation. A neural network (NN) method can accelerate the QM/MM-MD simulations, but it has long been a problem to accurately describe the QM/MM electrostatic coupling by NN in the electrostatic embedding (EE) scheme. In this work, we developed a new method to accelerate QM/MM calculations in the mechanic embedding (ME) scheme. The potentials and partial point charges of QM atoms are first learned in vacuo by the embedded atom neural networks (EANN) approach. MD simulations are then performed on this EANN/MM potential energy surface (PES) to obtain free energy (FE) profiles for reactions, in which the QM/MM electrostatic coupling is treated in the mechanic embedding (ME) scheme. Finally, a weighted thermodynamic perturbation (wTP) corrects the FE profiles in the ME scheme to the EE scheme. For two reactions in water and one in methanol, our simulations reproduced the B3LYP/MM free energy profiles within 0.5 kcal/mol with a speed-up of 30-60-fold. The results show that the strategy of combining EANN potential in the ME scheme with the wTP correction is efficient and reliable for chemical reaction simulations in liquid. Another advantage of our method is that the QM PES is independent of the MM subsystem, so it can be applied to various MM environments as demonstrated by an SN2 reaction studied in water and methanol individually, which used the same EANN PES. The free energy profiles are in excellent accordance with the results obtained from B3LYP/MM-MD simulations. In future, this method will be applied to the reactions of enzymes and their variants.

Accurate and efficient molecular dynamics based on machine learning and non von Neumann architecture

A Heterogeneous Parallel Non-von Neumann Architecture System for Accurate and Efficient Machine Learning Molecular Dynamics

High-speed and low-power molecular dynamics processing unit (MDPU) with ab initio accuracy

Accurate interatomic force field for molecular dynamics simulation by hybridizing classical and machine learning potentials

Scaling Molecular Dynamics with ab initio Accuracy to 149 Nanoseconds per Day

Ab-Initio Molecular Dynamics Acceleration Scheme with an Adaptive Machine Learning Framework

Pushing the limit of molecular dynamics with ab initio accuracy to 100 million atoms with machine learning

NNP/MM: Accelerating molecular dynamics simulations with machine learning potentials and molecular mechanic

Accurate and scalable graph neural network force field and molecular dynamics with direct force architecture

NNP/MM: Accelerating Molecular Dynamics Simulations with Machine Learning Potentials and Molecular Mechanics

Accurate and scalable multi-element graph neural network force field and molecular dynamics with direct force architecture

Extending the limit of molecular dynamics with ab initio accuracy to 10 billion atoms

TensorMD: Scalable Tensor-Diagram based Machine Learning Interatomic Potential on Heterogeneous Many-Core Processors

Efficient Machine Learning Force Field for Large-Scale Molecular Simulations of Organic Systems

Advancing Nonadiabatic Molecular Dynamics Simulations for Solids: Achieving Supreme Accuracy and Efficiency with Machine Learning

Deep Potential Molecular Dynamics: a scalable model with the accuracy of quantum mechanics

Constructing accurate and efficient general-purpose atomistic machine learning model with transferable accuracy for quantum chemistry

Accelerated Quantum Mechanics/Molecular Mechanics Simulations Via Neural Networks Incorporated with Mechanical Embedding Scheme.

Simulate Time-integrated Coarse-grained Molecular Dynamics with Geometric Machine Learning

Ab initio Accuracy Neural Network Potential for Drug-like Molecules

Molecular Dynamics with Neural-Network Potentials