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.

Correction: A Δ-machine learning approach for force fields, illustrated by a CCSD(T) 4-body correction to the MB-pol water potential

Learning the Quantum Centroid Force Correction in Molecular Systems: A Localized Approach

Systematic Improvement of Empirical Energy Functions in the Era of Machine Learning

Correction: How fast do defects migrate in halide perovskites: insights from on-the-fly machine-learned force fields

A Transferrable Range-Separated Force Field for Water: Combining the Power of Both Physically-Motivated Models and Machine Learning Techniques

A semilocal machine-learning correction to density functional approximations

Correction: Three-dimensional Numerical Simulation and Experimental Investigation of Boundary-Driven Streaming in Surface Acoustic Wave Microfluidics

Correction: Kinetics and Dynamics of the C(3P) + H2O Reaction on a Full-Dimensional Accurate Triplet State Potential Energy Surface

Correction: Enhanced 5F-Δ Bonding in [U(C7H7)2]-: C K-edge XAS, Magnetism, and Ab Initio Calculations.

Correction: Enhanced 5F-Δ Bonding in [U(C7H7)2]−: C K-edge XAS, Magnetism, and Ab Initio Calculations

Self-consistent density-based basis-set correction: How much do we lower total energies and improve dipole moments?

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

Machine learning for predictive condensed-phase simulation

Correction: Supersolidity of undercoordinated and hydrating water

Δ-Machine Learning to Elevate DFT-Based Potentials and a Force Field to the CCSD( T ) Level Illustrated for Ethanol

Linear regression correction to first principle theoretical calculations – Improved descriptors and enlarged training set

Correction: Corrigendum: on the Speed of Piezostrain-Mediated Voltage-Driven Perpendicular Magnetization Reversal: a Computational Elastodynamics-Micromagnetic Phase-Field Study

Corrigendum: Exploring the role of beyond mean-field interaction in the structure and dynamics of one-dimensional quantum droplets (2023 J. Phys. B: At. Mol. Opt. Phys. 56 165302)

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

Geometrical correction for the inter- and intramolecular basis set superposition error in periodic density functional theory calculations