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.

Studies of Arithmetic Coupling Between Molecular Dynamics and Neural Network

Large Scale and Quantum Accurate Molecular Dynamics Simulation:Liquid Iron under Extreme Condition

The Development and Comparison of Molecular Dynamics Simulation and Monte Carlo Simulation

A Machine Learning Approach to Calculate Electronic Couplings Between Quasi-Diabatic Molecular Orbitals: the Case of DNA

Towards Full Ab Initio Simulations of Mechanical Response of Ti-N Binary System under Extremely Operational Conditions Via Neural Network Potential

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

Numerical analysis of the influence of initial and external conditions on the association of artificial atoms

A quantitative study on the approximation error and speed-up of the multi-scale MCMC (Monte Carlo Markov chain) method for molecular dynamics

Analyzing Atomic Interactions in Molecules as Learned by Neural Networks

First-principles Nonadiabatic Dynamics of Molecules at Metal Surfaces with Vibrationally Coupled Electron Transfer

Simulation of chemical reaction dynamics based on quantum computing

Ab initio Accuracy Neural Network Potential for Drug-like Molecules

Predicting scalar coupling constants by graph angle-attention neural network

Parallel Smoothed Molecular Dynamics Method and Coupling with Molecular Dynamics

Molecular Dynamics with Neural-Network Potentials

Neural network representation of electronic structure from $ab$ $initio$ molecular dynamics

Highly Accurate Description of Long-Range Interactions through the Combination of Neural Networks and Physical Models

Coupling the molecular motion and collision processes in numerical simulations

Simulate Time-integrated Coarse-grained Molecular Dynamics with Multi-Scale Graph Networks

Recent Development In Quantum Mechanics/Molecular Mechanics Modeling For Materials