Solvent effects on the S(N)2 reaction: Application of the density functional theory-based effective fragment potential method.

The performance of the density functional theory (DFT)-based effective fragment potential (EFP) method is assessed using the S(N)2 reaction: Cl- + nH2O + CH3Br = CH3Cl + Br- + nH2O. The effect of the systematic addition of water molecules on the structures and relative energies of all species involved in the reaction has been studied. The EFP1 method is compared with second-order perturbation theory (MP2) and DFT results for n = 1, 2, and 3, and EFP1 results are also presented for four water molecules. The incremental hydration effects on the barrier height are the same for all methods. However, only full MP2 or MP2 with EFP1 solvent molecules are able to provide an accurate treatment of the transition state (TS) and hence the central barriers. Full DFT and DFT with EFP1 solvent molecules both predict central barriers that are too small. The results illustrate that the EFP1-based DFT method gives reliable results when combined with an accurate quantum mechanical (QM) method, so it may be used as an efficient alternative to fully QM methods in the treatment of larger microsolvated systems.

[1]  Anthony J. Stone,et al.  Distributed multipole analysis, or how to describe a molecular charge distribution , 1981 .

[2]  H. Schaefer.  Methods of Electronic Structure Theory , 1977 .

[3]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[4]  M. Krauss,et al.  Solvation and the excited states of formamide , 1997 .

[5]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[6]  B. Grigorenko,et al.  Hydrogen bonding at the diatomics-in-molecules level: Water clusters , 2000 .

[7]  Mark S. Gordon,et al.  A combined discrete/continuum solvation model: Application to glycine , 2000 .

[8]  Shridhar R. Gadre,et al.  Structure and Stability of Water Clusters (H2O)n, n ) 8-20: An Ab Initio Investigation , 2001 .

[9]  R. Levy,et al.  Molecular dynamics simulation of time-resolved fluorescence and nonequilibrium solvation of formaldehyde in water , 1990 .

[10]  Mark S. Gordon,et al.  Evaluation of Charge Penetration Between Distributed Multipolar Expansions , 2000 .

[11]  Mark S. Gordon,et al.  The Effective Fragment Potential Method: A QM-Based MM Approach to Modeling Environmental Effects in Chemistry , 2001 .

[12]  P. N. Day,et al.  A study of water clusters using the effective fragment potential and Monte Carlo simulated annealing , 2000 .

[13]  Jan H. Jensen,et al.  Accurate Intraprotein Electrostatics Derived from First Principles: An Effective Fragment Potential Method Study of the Proton Affinities of Lysine 55 and Tyrosine 20 in Turkey Ovomucoid Third Domain , 2001 .

[14]  H. L. Hartley,et al.  Manuscript Preparation , 2022 .

[15]  Andreas Klamt,et al.  Incorporation of solvent effects into density functional calculations of molecular energies and geometries , 1995 .

[16]  Jacopo Tomasi,et al.  An Integrated Effective Fragment—Polarizable Continuum Approach to Solvation: Theory and Application to Glycine , 2002 .

[17]  E. Clementi,et al.  Small Clusters of Water Molecules Using Density Functional Theory , 1996 .

[18]  Wolfram Koch,et al.  A Chemist's Guide to Density Functional Theory , 2000 .

[19]  Mark S. Gordon,et al.  An Approximate Formula for the Intermolecular Pauli Repulsion Between Closed Shell Molecules. II. Application to the Effective Fragment Potential Method , 1998 .

[20]  Orlando Tapia,et al.  Solvent effects and chemical reactivity , 2002 .

[21]  Han Myoung Lee,et al.  Structures, energies, vibrational spectra, and electronic properties of water monomer to decamer , 2000 .

[22]  J. J. Dannenberg,et al.  Effect of Basis Set Superposition Error on the Water Dimer Surface Calculated at Hartree-Fock, Møller-Plesset, and Density Functional Theory Levels , 1999 .

[23]  M. Gordon,et al.  On the Number of Water Molecules Necessary To Stabilize the Glycine Zwitterion , 1995 .

[24]  H. Rzepa,et al.  An AM1 and PM3 molecular orbital and self-consistent reaction-field study of the aqueous solvation of glycine, alanine and proline in their neutral and zwitterionic forms , 1991 .

[25]  Mark S. Gordon,et al.  The Effective Fragment Model for Solvation: Internal Rotation in Formamide , 1996 .

[26]  Jan H. Jensen,et al.  The Prediction of Protein pKa's Using QM/MM: The pKa of Lysine 55 in Turkey Ovomucoid Third Domain , 2002 .

[27]  Jacopo Tomasi,et al.  Molecular Interactions in Solution: An Overview of Methods Based on Continuous Distributions of the Solvent , 1994 .

[28]  Jan H. Jensen,et al.  Modeling intermolecular exchange integrals between nonorthogonal molecular orbitals , 1996 .

[29]  J. Langlet,et al.  Theoretical study of solvent effect on intramolecular proton transfer of glycine , 2000 .