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.

Solvent Effects on the Symmetric and Asymmetric SN2 Reactions in the Acetonitrile Solution: A Reaction Density Functional Theory Study.

Dynamical Effects of SN2 Reactivity Suppression by Microsolvation: Dynamics Simulations of the F-(H2O) + CH3I Reaction on a 21-Dimensional Potential Energy Surface.

Dynamical density functional theory for solvation dynamics in polar solvent: Heterogeneous effect of solvent orientation

Solvent effects on Diels-Alder reaction in ionic liquids: A reaction density functional study

A Molecular Electron Density Theory Study of the Bimolecular Nucleophilic Substitution Reactions on Monosubstituted Methyl Compounds

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

DFT Study of Inner-Sphere Electron Transfer in a Gas-Phase SN2 Reaction at the Saturated Carbon

Advances in Solvent Effect of Quaternarization

Solvation Free Energies in Subsystem Density Functional Theory

Solvent Effect on Concertedness of the Transition State in the Hydrolysis of P-Nitrophenyl Acetate.

Modification of the Onsager Reaction Field and Its Application on Spectral Parameters.

Correlation between the energy and electron density representations of reactivity: mPW1K study of the asymmetric SN2 reactions at the saturated nitrogen

Tackling solvent effect by coupling electronic and molecular Density Functional Theory

Molecular Insights on Solvent Effects in Organic Reactions as Obtained through Computational Chemistry Tools

Investigating the role of halogen-bonded complexes in microsolvated Y-(H2O)n + CH3I SN2 reactions

Solvent reaction coordinate for an S$_N$2 reaction

Solvent Effect on Excited-State Intramolecular Proton Transfer Process Based on 'naked' Diazaborepins

Effect of solvent on the initiation mechanism of living anionic polymerization of styrene: A computational study

Reactivity and stability of ion pairs, dimers and tetramers versus solvent polarity: SNAr fluorination of 2-bromobenzonitrile with tetramethylammonium fluoride

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Solvation Effects on Quantum Tunneling Reactions