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.

Solubility Prediction of Organic Molecules with Molecular Dynamics Simulations

Solubility of Organic Salts in Solvent–Antisolvent Mixtures: A Combined Experimental and Molecular Dynamics Simulations Approach

A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals

Polymorph-Specific Solubility Prediction using Constant Chemical Potential Molecular Dynamics

Organic Solubility Prediction at the Limit of Aleatoric Uncertainty

Molecular Mechanism of Gas Solubility in Liquid: Constant Chemical Potential Molecular Dynamics Simulations

Predicting Solubility Limits of Organic Solutes for a Wide Range of Solvents and Temperatures

Polymorph-Specific Solubility Prediction of Urea Using Constant Chemical Potential Molecular Dynamics Simulations

Simulation of Crystal Growth by an Innovative Hybrid Density Functional Theory Continuum Solvation Approach: Kink Site Formation on Barite (001)

Studies on the Conformational Interconversion of Sulfonylurea Molecules in Solution by Molecular Dynamics Simulation Method

Solvent hydrodynamics speed up crystal nucleation in suspensions of hard spheres

Enhancing Predictions of Drug Solubility through Multidimensional Structural Characterization Exploitation

A general graph neural network based implicit solvation model for organic molecules in water

General Graph Neural Network Based Implicit Solvation Model for Organic Molecules in Water

Sampling of Organic Solutes in Aqueous and Heterogeneous Environments Using Oscillating Excess Chemical Potentials in Grand Canonical-like Monte Carlo-Molecular Dynamics Simulations

Solution Thermodynamics of L-Glutamic Acid Polymorphs from Finite-Sized Molecular Dynamics Simulations

Nucleation and relaxation processes in weak solutions: molecular dynamics simulation

Unravelling guest dynamics in crystalline molecular organics using solid-state NMR and molecular dynamics simulation

Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways

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

Modeling the solubility of pharmaceuticals in pure solvents and solvent mixtures for drug process design