Analytical energy gradients of a self-consistent reaction-field solvation model based on CM2 atomic charges

Analytical energy gradients have been derived for an SM5-type solvation model based on Hartree–Fock self-consistent reaction-field theory and CM2 atomic charges. The method is combined with an analytic treatment of the first derivatives of nonelectrostatic first-solvation-shell contributions to the free energy and implemented in the General Atomic and Molecular Electronic Structure System (GAMESS). The resulting equations allow one to use accurate class IV charges to calculate equilibrium geometries of solutes in liquid-phase solutions. The algorithm is illustrated by calculations of optimized geometries and solvation free energies for water, methanol, dimethyl disulfide, and 9-methyladenine in water and 1-octanol.

[1]  Manuel F. Ruiz-López,et al.  Computation of hydration free energies using a parameterized continuum model: Study of equilibrium geometries and reactive processes in water solution , 1996, J. Comput. Chem..

[2]  David J. Giesen,et al.  Solvation Model for Chloroform Based on Class IV Atomic Charges , 1997 .

[3]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[4]  Warren J. Hehre,et al.  AB INITIO Molecular Orbital Theory , 1986 .

[5]  Donald G. Truhlar,et al.  New Class IV Charge Model for Extracting Accurate Partial Charges from Wave Functions , 1998 .

[6]  Donald G. Truhlar,et al.  The interface of electronic structure and dynamics for reactions in solution , 1998 .

[7]  D. Rinaldi,et al.  Fast geometry optimizationin self‐cosistent reaction field computations on solvated molecules , 1992 .

[8]  David J. Giesen,et al.  A universal model for the quantum mechanical calculation of free energies of solvation in non-aqueous solvents , 1997 .

[9]  Richard H. Bartels,et al.  Algorithm 432 [C2]: Solution of the matrix equation AX + XB = C [F4] , 1972, Commun. ACM.

[10]  István Mayer,et al.  Charge, bond order and valence in the AB initio SCF theory , 1983 .

[11]  David J. Giesen,et al.  General Semiempirical Quantum Mechanical Solvation Model for Nonpolar Solvation Free Energies. n-Hexadecane , 1995 .

[12]  Donald G. Truhlar,et al.  Universal reaction field model based on ab initio Hartree–Fock theory , 1998 .

[13]  C. Cramer,et al.  PM3‐SM3: A general parameterization for including aqueous solvation effects in the PM3 molecular orbital model , 1992 .

[14]  Jacopo Tomasi,et al.  Continuum solvation models: A new approach to the problem of solute’s charge distribution and cavity boundaries , 1997 .

[15]  Jiabo Li,et al.  MIDI! basis set for silicon, bromine, and iodine , 1998 .

[16]  C. Cramer,et al.  An SCF Solvation Model for the Hydrophobic Effect and Absolute Free Energies of Aqueous Solvation , 1992, Science.

[17]  Donald G. Truhlar,et al.  Density functional solvation model based on CM2 atomic charges , 1998 .

[18]  Gregory D. Hawkins,et al.  Parametrized Models of Aqueous Free Energies of Solvation Based on Pairwise Descreening of Solute Atomic Charges from a Dielectric Medium , 1996 .

[19]  C. Cramer.  Where is the unpaired electron in the phosphoranyl radicals H3PS− and H3PSH? , 1993 .

[20]  M. Frisch,et al.  Comparison of Local, Nonlocal, and Hybrid Density Functionals Using Vibrational Absorption and Circular Dichroism Spectroscopy , 1996 .

[21]  A. Klamt,et al.  First Principles Implementation of Solvent Effects Without Outlying Charge Error , 1997 .

[22]  Jean-Louis Rivail,et al.  Analytical energy derivatives for a realistic continuum model of solvation: Application to the analysis of solvent effects on reaction paths , 1996 .

[23]  Donald G. Truhlar,et al.  omnisol: Fast Prediction of Free Energies of Solvation and Partition Coefficients , 1998 .

[24]  Ian M. Mills,et al.  Force Constants and Dipole-Moment Derivatives of Molecules from Perturbed Hartree-Fock Calculations. I , 1968 .

[25]  Clemens C. J. Roothaan,et al.  New Developments in Molecular Orbital Theory , 1951 .

[26]  David J. Giesen,et al.  The MIDI! basis set for quantum mechanical calculations of molecular geometries and partial charges , 1996 .

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

[28]  J. Tomasi,et al.  Analytical derivatives for geometry optimization in solvation continuum models. II. Numerical applications , 1998 .

[29]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[30]  David J. Giesen,et al.  Class IV charge models: A new semiempirical approach in quantum chemistry , 1995, J. Comput. Aided Mol. Des..

[31]  C. Cramer,et al.  General parameterized SCF model for free energies of solvation in aqueous solution , 1991 .

[32]  David J. Giesen,et al.  A Universal Organic Solvation Model. , 1996, The Journal of organic chemistry.

[33]  P. Löwdin.  On the Non‐Orthogonality Problem Connected with the Use of Atomic Wave Functions in the Theory of Molecules and Crystals , 1950 .

[34]  T. Truong,et al.  Analytical first and second energy derivatives of the generalized conductorlike screening model for free energy of solvation , 1995 .

[35]  Donald G. Truhlar,et al.  Improved methods for semiempirical solvation models , 1995, J. Comput. Chem..

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

[37]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[38]  Peter Pulay,et al.  Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules , 1969 .

[39]  Donald G. Truhlar,et al.  MODEL FOR AQUEOUS SOLVATION BASED ON CLASS IV ATOMIC CHARGES AND FIRST SOLVATION SHELL EFFECTS , 1996 .

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

[41]  Donald G. Truhlar,et al.  Generalized born fragment charge model for solvation effects as a function of reaction coordinate , 1989 .

[42]  Donald G. Truhlar,et al.  Universal Quantum Mechanical Model for Solvation Free Energies Based on Gas-Phase Geometries , 1998 .

[43]  D. Truhlar,et al.  Effect of nonequilibrium solvation on chemical reaction rates. Variational transition-state-theory studies of the microsolvated reaction Cl-(H2O)n + CH3Cl , 1990 .