Abstract:Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations have been developed to simulate molecular systems, where an explicit description of changes in the electronic structure is necessary. However, QM/MM MD simulations are computationally expensive compared to fully classical simulations as all valence electrons are treated explicitly and a self-consistent field (SCF) procedure is required. Recently, approaches have been proposed to replace the QM description with machine learned (ML) models. However, condensed-phase systems pose a challenge for these approaches due to long-range interactions. Here, we establish a workflow, which incorporates the MM environment as an element type in a high-dimensional neural network potential (HDNNP). The fitted HDNNP describes the potential-energy surface of the QM particles with an electrostatic embedding scheme. Thus, the MM particles feel a force from the polarized QM particles. To achieve chemical accuracy, we find that even simple systems require models with a strong gradient regularization, a large number of data points, and a substantial number of parameters. To address this issue, we extend our approach to a delta-learning scheme, where the ML model learns the difference between a reference method (DFT) and a cheaper semi-empirical method (DFTB). We show that such a scheme reaches the accuracy of the DFT reference method, while requiring significantly less parameters. Furthermore, the delta-learning scheme is capable of correctly incorporating long-range interactions within a cutoff of 1.4 nm. It is validated by performing MD simulations of retinoic acid in water and the interaction between S-adenoslymethioniat with cytosine in water. The presented results indicate that delta-learning is a promising approach for (QM)ML/MM MD simulations of condensed-phase systems.

Learning emergent models from $\textit{ab initio}$ many-body calculations

Learning Effective Molecular Models from Experimental Observables.

Downfolding from Ab Initio to Interacting Model Hamiltonians: Comprehensive Analysis and Benchmarking of the DFT+cRPA Approach

Ab-Initio Approach to Many-Body Quantum Spin Dynamics

Ab Initio Validation on the Connection Between Atomistic and Hydrodynamic Description to Unravel the Ion Dynamics of Warm Dense Matter

Renormalization of States and Quasiparticles in Many-body Downfolding

Learning from the Density to Correct Total Energy and Forces in First Principle Simulations

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

Requirements for building effective Hamiltonians using quantum-enhanced density matrix downfolding

Embedding parameters in ab initio theory to develop approximations based on molecular similarity

Density Matrix Renormalization Group Approach Based on the Coupled-Cluster Downfolded Hamiltonians

Reduced density matrix functional theory from an ab initio seniority-zero wave function: Exact and approximate formulations along adiabatic connection paths

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

From Many-Body Ab Initio to Effective Excitonic Models: A Versatile Mapping Approach Including Environmental Embedding Effects

Semiempirical Hamiltonians learned from data can have accuracy comparable to Density Functional Theory

Atom-wise formulation of the many-body dispersion problem for linear-scaling van der Waals corrections

Ab initio electron-lattice downfolding: potential energy landscapes, anharmonicity, and molecular dynamics in charge density wave materials

Hierarchical Mapping for Efficient Simulation of Strong System-Environment Interactions