Abstract:Direct molecular dynamics (MD) simulation with ab initio quantum mechanical and molecular mechanical (QM/MM) methods is very powerful for studying the mechanism of chemical reactions in a complex environment but also very time-consuming. The computational cost of QM/MM calculations during MD simulations can be reduced significantly using semiempirical QM/MM methods with lower accuracy. To achieve higher accuracy at the ab initio QM/MM level, a correction on the existing semiempirical QM/MM model is an attractive idea. Recently, we reported a neural network (NN) method as QM/MM-NN to predict the potential energy difference between semiempirical and ab initio QM/MM approaches. The high-level results can be obtained using neural network based on semiempirical QM/MM MD simulations, but the lack of direct MD samplings at the ab initio QM/MM level is still a deficiency that limits the applications of QM/MM-NN. In the present paper, we developed a dynamic scheme of QM/MM-NN for direct MD simulations on the NN-predicted potential energy surface to approximate ab initio QM/MM MD. Since some configurations excluded from the database for NN training were encountered during simulations, which may cause some difficulties on MD samplings, an adaptive procedure inspired by the selection scheme reported by Behler [Behler Int. J. Quantum Chem. 2015, I IS, 1032; Behler Angew. Chem., Int. Ed. 2017, 56, 12828] was employed with some adaptions to update NN and carry out MD iteratively. We further applied the adaptive QM/MM-NN MD method to the free energy calculation and transition path optimization on chemical reactions in water. The results at the ab initio QM/MM level can be well reproduced using this method after 2-4 iteration cycles. The saving in computational cost is about 2 orders of magnitude. It demonstrates that the QM/MM-NN with direct MD simulations has great potentials not only for the calculation of thermodynamic properties but also for the characterization of reaction dynamics, which provides a useful tool to study chemical or biochemical systems in solution or enzymes.

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

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

Machine Learning for Parameter Auto-tuning in Molecular Dynamics Simulations: Efficient Dynamics of Ions near Polarizable Nanoparticles

OpenMM 8: Molecular Dynamics Simulation with Machine Learning Potentials

Molecular Dynamics with Neural-Network Potentials

Machine Learning for Performance Enhancement of Molecular Dynamics Simulations

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

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

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

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

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

Accurate Free Energy Calculation via Multiscale Simulations Driven by Hybrid Machine Learning and Molecular Mechanics Potentials

Machine Learning in QM/MM Molecular Dynamics Simulations of Condensed-Phase Systems

Molecular Dynamics Simulations with Quantum Mechanics/Molecular Mechanics and Adaptive Neural Networks.

Advancing Multiscale Molecular Modeling with Machine Learning-Derived Electrostatics

Accelerating ab initio QM/MM Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semi-empirical QM/MM Hamiltonian

Learning Interatomic Potentials at Multiple Scales

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

Accelerating Metadynamics-Based Free-Energy Calculations with Adaptive Machine Learning Potentials.

A Generator of Neural Network Potential for Molecular Dynamics: Constructing Robust and Accurate Potentials with Active Learning for Nanosecond-scale Simulations