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.

Density Functional Calculation and Monte Carlo Simulation of Solvating Effect of Tautomerism and Proton Transfer of 2-Amino-2-imidazoline

DFT study and Monte Carlo simulation on proton transfers of 2-amino-2-oxazoline, 2-amino-2-thiazoline, and 2-amino-2-imidazoline in the gas phase and in water.

Theoretical Study of the Effect of Water in the Process of Proton Transfer of Glycinamide

Tackling solvent effect by coupling electronic and molecular Density Functional Theory

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

Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model.

A molecular density functional theory to study solvation in water

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

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Theoretical Study on Tautomerism and Proton Transfer of 2,4-Dithiothymine

A Density Functional Theory Study On The Intramolecular Proton Transfer In The Enzyme Carbonic Anhydrase

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

Molecular Theories and Simulation of Ions and Polar Molecules in Water

Comparison of DFT, MP2/CBS, and CCSD(T)/CBS methods for a dual-level QM/MM Monte Carlo simulation approach calculating the free energy of activation of reactions in solution and “on water”: a case study

Density Functional Theory Study of Proton Transfer in Carbonic Anhydrase

Systematic Study of the Tautomerism of Uracil Induced by Proton Transfer. Exploration of Water Stabilization and Mutagenicity

Density Functional Studies on Formamide′s Water-influenced Proton Transfer in Solution

A theoretical study of solvent effects on tautomerism and electronic absorption spectra of 3-hydroxy-2-mercaptopyridine and 2,3-dihydroxypyridine.

Intramolecular and water mediated tautomerism of solvated glycine

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

Theoretical Calculations of the Tautomeric Equilibrium in 3-X-2(1H)-Pyridones