Journal of Physical Chemistry, Vol.100, No.5, 1515-1523, 1996
Solvation Free-Energy Calculations Using a Continuum Dielectric Model for the Solvent and Gradient-Corrected Density-Functional Theory for the Solute
Electrostatic solvation free energies are calculated using a self consistent-reaction field(SCRF) procedure that combines a continuum dielectric model of the solvent with both Hartree-Fock (HF) and density functional theory (DFT) for the solute. Several molecules are studied in aqueous solution. They comprise three groups : nonpolar neutral, polar neutral, and ionic. The calculated values of Delta G(el) are sensitive to the atomic radii used to define the solute molecular surface, particularly to the value of the hydrogen radius. However, the values of Delta G(el) exhibit reasonable correlation with experiment when a previously determined, physically motivated set of atomic radii were used to define the van der Waals surface of the solute. The standard deviation between theory and experiment is 2.51 kcal/mol. for HF and 2.21 kcal/mol for DFT for the 14 molecules examined. The errors with HF or DFT are similar. The relative difference between the calculated values.of Delta G(el) and experiment is largest for nonpolar neutral molecules, intermediate for polar neutral molecules, and smallest for ions. This is consistent with the expected relative importance of nonelectrostatic contributions to the free energy that are omitted in the model.
Keywords:REACTION-FIELD METHOD;ELECTRON CORRELATION;AQUEOUS-SOLUTION;VIBRATIONAL FREQUENCIES;WATER MOLECULE;FORCE-FIELD;ABINITIO;HYDRATION;CAVITY;POTENTIALS