Abstract:Simulations of chemical reactivity in condensed phase systems represent an ongoing challenge in computational chemistry, where traditional quantum chemical approaches typically struggle with both the size of the system and the potential complexity of the reaction. Here, we introduce a workflow aimed at efficiently training neural network potentials (NNPs) to explore energy barriers in solution at the hybrid density functional theory level. The computational burden associated with training at the PBE0-D3(BJ) level is bypassed through the use of active and transfer learning techniques, whereas extensive sampling of the transition state region is accelerated by well-tempered metadynamics simulations using multiple time-step integration. These NNPs serve to explore a puzzling solute--solvent reactivity route involving the ring opening of N-enoxyphthalimide experimentally observed in methanol but not in 2,2,2-trifluoroethanol (TFE). This reaction represents a challenging example characterized by intricate hydrogen bonding networks and structurally ambiguous solvent-sensitive transition states. The methodology successfully delivers detailed free energy surfaces and relative energy barriers in quantitative agreement with experiment. These barriers are associated with an ensemble of transition states involving direct participation of up to five solvent molecules. While this picture contrasts with the single transition state structure assumed by current static models, no drastic qualitative difference is observed between the formed hydrogen bonding networks and the number of participating solvent molecules in methanol or TFE. The dichotomy between the two solvents thus essentially arises from an electronic effect (i.e., distinct nucleophilicity) and from the larger conformational entropy contributions in methanol. This example underscores the critical role dynamic simulations at the ab initio levels play in capturing the full complexity of solute-solvent interactions.

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Molecular Dynamics Simulation of Zinc Ion in Water with an Ab Initio Based Neural Network Potential

Many-body interactions and deep neural network potentials for water

A "short blanket" dilemma for a state-of-the-art neural network potential for water: Reproducing experimental properties or the physics of the underlying many-body interactions?

Polarizable Water Model with Ab Initio Neural Network Dynamic Charges and Spontaneous Charge Transfer

Enhanced Deep Potential Model for Fast and Accurate Molecular Dynamics; Application to the Hydrated Electron

Enhanced deep potential model for fast and accurate molecular dynamics: application to the hydrated electron

Scalable molecular simulation of electrolyte solutions with quantum chemical accuracy

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

Probing the self-ionization of liquid water with ab initio deep potential molecular dynamics

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

Homogeneous water vapor condensation with a deep neural network potential model

Calculation of solvation force in molecular dynamics simulation by deep-learning method

Modelling infrared spectra of the O-H stretches in liquid H2O based on a deep learning potential, the importance of nuclear quantum effects

Understanding the Dielectric Relaxation of Liquid Water Using Neural Network Potential and Classical Pairwise Potential

Modeling Liquid Water by Climbing Up Jacob's Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Bottom-up Derived Flexible Water Model with Dipole and Quadrupole Moments for Coarse-Grained Molecular Simulations.

Capturing Dichotomic Solvent Behavior in Solute--Solvent Reactions with Neural Network Potentials

Efficient Molecular Dynamics Simulations of Deep Eutectic Solvents with First-Principles Accuracy Using Machine Learning Interatomic Potentials