Abstract:Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype SN2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

What problem does this paper attempt to address?

### What problems does this paper attempt to solve? This paper aims to solve the difficult problem of understanding the mechanism of chemical reactions in complex solutions, especially the interactions between the atomic scale and the mesoscale. Specifically, the authors have developed a new hybrid method that combines ab initio molecular dynamics (AIMD) with classical density functional theory (cDFT) to overcome the limitations of existing solvation models in terms of time, length, and system complexity. #### Main problems: 1. **Challenges of multi - scale simulation**: - Existing molecular simulation techniques have limitations in time and space scales and are difficult to describe the reaction mechanisms in complex solutions. - Currently available solvation models can rarely describe the interactions between the molecular scale and the mesoscale. 2. **Understanding of solvent effects**: - The influence of solvents on chemical reactions is complex, especially in solutions, where solvent molecules constantly interact with reactants, affecting the reaction paths and rates. - A more accurate method is needed to describe how solvents affect the reaction coordinates, especially for prototype reactions such as SN2 reactions. 3. **Calculation of solvation energy**: - Existing continuum solvation models (such as SMD) are efficient but not accurate enough, especially when dealing with complex molecules. - A new method that can provide reasonable solvation energy is needed to better match experimental results. #### Solutions: To meet these challenges, the authors proposed a hybrid AIMD/cDFT method with the following main features: - **Dynamic coupling**: By using AIMD to simulate the quantum - mechanical behavior of the region of interest and cDFT to describe the mesoscopic properties of the surrounding environment, a self - consistent ground state is achieved. - **High precision and efficiency**: cDFT can provide an accurate thermodynamic description of electrolyte solutions in a wide concentration range (up to 2M). - **Application verification**: The effectiveness of this method was verified by calculating the solvation energies of various molecules and ions and comparing them with the widely used SMD model and experimental data. - **Reaction path analysis**: The nucleophilic attack of chloride ions on chloromethane (SN2 reaction) was studied, revealing the specific influence of solvent effects on the reaction path energy barrier. Through this method, researchers can gain a deeper understanding of the role of solvents in complex chemical reactions and provide a powerful tool for future research.

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

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance

Coupling Molecular Density Functional Theory with Converged Selected Configuration Interaction Methods to Study Excited states in Aqueous Solution

Hybrid Functional and Plane Waves based Ab Initio Molecular Dynamics Study of the Aqueous Fe$^{2+}$/Fe$^{3+}$ Redox Reaction

Scalable molecular simulation of electrolyte solutions with quantum chemical accuracy

Hybrid Functional and Plane Waves based Ab Initio Molecular Dynamics Study of the Aqueous Fe2+ /Fe3+ Redox Reaction

Multi-scale theory in the molecular simulation of electrolyte solutions

Quantum mechanical and molecular mechanical simulation approaches bridging length and time scales for simulation of interface reactions in realistic environments

Speeding-up Hybrid Functional based Ab Initio Molecular Dynamics using Multiple Time-stepping and Resonance Free Thermostat

Tackling solvent effect by coupling electronic and molecular Density Functional Theory

Achieving an Order of Magnitude Speed-up in Hybrid Functional and Plane Wave based Ab Initio Molecular Dynamics: Applications to Proton Transfer Reactions in Enzymes and in Solution

Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions

Density-Functionalized QM/MM Delivers Chemical Accuracy For Solvated Systems

Modeling solvation effects in real-space and real-time within Density Functional Approaches

Exploring A Multi-Scale Method for Molecular Simulations in Continuum Solvent Model: Explicit Simulation of Continuum Solvent As An Incompressible Fluid

Doubly Hybrid Functionals Close to Chemical Accuracy for Both Finite and Extended Systems: Implementation and Test of XYG3 and XYGJ-OS.

Hybrid Combination of Quantum Mechanics with Quantum-Based Polarizable Reactive Force Field for Large Scale Full Solvent Simulations of Electrocatalysis

A Seamless Grid-Based Interface for Mean-Field QM/MM Coupled with Efficient Solvation Free Energy Calculations

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

Generating High-Precision Force Fields for Molecular Dynamics Simulations to Study Chemical Reaction Mechanisms using Molecular Configuration Transformer

Density Functionals with Quantum Chemical Accuracy: From Machine Learning to Molecular Dynamics