Hydration of ionic species studied by the reference interaction site model with a repulsive bridge correction

We have tested the reference interaction site model (RISM) for the case of the hypernetted chain (HNC) and the partially linearized hypernetted chain (PLHNC) closures improved by a repulsive bridge correction (RBC) for ionic hydrated species. We have analyzed the efficiency of the RISM/HNC+RBC and RISM/PLHNC+RBC techniques for decomposition of the electrostatic and the nonpolar hydration energies on the energetic and the enthalpic parts for polyatomic ions when the repulsive bridge correction is treated as a thermodynamic perturbation, and investigate the repulsive bridge effect on the electrostatic potential induced by solvent on solute atoms. For a number of univalent and bivalent atomic ions, molecular cations, and anions, the method provides hydration energies deviating only by several percents from the experimental data. In most cases, the enthalpic contributions to the free energies are also close to the experimental results. The above models are able to satisfactory predict the hydration energies as well as the electrostatic potential around the ionic species. For univalent atomic ions, they also provide qualitative estimates of the Samoilov activation energies. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2008

[1]  Fumio Hirata,et al.  A hybrid approach for the solvent effect on the electronic structure of a solute based on the RISM and Hartree-Fock equations , 1993 .

[2]  S. Rast,et al.  Structure and thermodynamics of liquid acetonitrile via Monte Carlo simulation and Ornstein-Zernike theories , 1997 .

[3]  C. Millot,et al.  A molecular Ornstein–Zernike study of popular models for water and methanol , 1999 .

[4]  Jason Crain,et al.  Improved estimates for hydration free energy obtained by the reference interaction site model , 2007 .

[5]  F. Hirata,et al.  Three-dimensional density profiles of water in contact with a solute of arbitrary shape: a RISM approach , 1998 .

[6]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[7]  D. Ives,et al.  Structure of aqueous electrolyte solutions and the hydration of ions , 1965 .

[8]  David Chandler,et al.  Optimized Cluster Expansions for Classical Fluids. II. Theory of Molecular Liquids , 1972 .

[9]  S. Rast,et al.  LIQUID ACETONE AND CHLOROFORM : A COMPARISON BETWEEN MONTE CARLO SIMULATION, MOLECULAR ORNSTEIN-ZERNIKE THEORY, AND SITE-SITE ORNSTEIN-ZERNIKE THEORY , 1998 .

[10]  A. Kovalenko,et al.  Self-consistent combination of the three-dimensional RISM theory of molecular solvation with analytical gradients and the Amsterdam density functional package. , 2006 .

[11]  Arieh Warshel,et al.  Langevin Dipoles Model for ab Initio Calculations of Chemical Processes in Solution: Parametrization and Application to Hydration Free Energies of Neutral and Ionic Solutes and Conformational Analysis in Aqueous Solution , 1997 .

[12]  Maxim V. Fedorov,et al.  Wavelet algorithm for solving integral equations of molecular liquids. A test for the reference interaction site model , 2004, J. Comput. Chem..

[13]  M. Fedorov,et al.  Wavelet treatment of structure and thermodynamics of simple liquids. , 2004, The Journal of chemical physics.

[14]  F. Hirata,et al.  Ion Hydration: Thermodynamic and Structural Analysis with an Integral Equation Theory of Liquids , 1997 .

[15]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[16]  M. B. Pinto,et al.  Optimized δ expansion for relativistic nuclear models , 1997, nucl-th/9709049.

[17]  Benoît Roux,et al.  An Integral Equation To Describe the Solvation of Polar Molecules in Liquid Water , 1997 .

[18]  Fumio Hirata,et al.  Chemical Processes in Solution Studied by an Integral Equation Theory of Molecular Liquids. , 1998 .

[19]  S. Ten-no,et al.  On the connection between the reference interaction site model integral equation theory and the partial wave expansion of the molecular Ornstein–Zernike equation , 1999 .

[20]  Benoît Roux,et al.  NUMERICAL SOLUTION OF THE HYPERNETTED CHAIN EQUATION FOR A SOLUTE OF ARBITRARY GEOMETRY IN THREE DIMENSIONS , 1995 .

[21]  H. Krienke,et al.  The solvation of ions in acetonitrile and acetone: A molecular Ornstein–Zernike study , 1998 .

[22]  Fumio Hirata,et al.  Self-consistent field, ab initio molecular orbital and three-dimensional reference interaction site model study for solvation effect on carbon monoxide in aqueous solution , 2000 .

[23]  Shigeyoshi Sakaki,et al.  Comparison of electronic structure theories for solvated molecules: RISM-SCF versus PCM , 2004 .

[24]  C. Brooks.  Computer simulation of liquids , 1989 .

[25]  A. Morita,et al.  The Charge Response Kernel with Modified Electrostatic Potential Charge Model , 2002 .

[26]  S. Sakaki,et al.  Alternative couplings of solute–solvent interaction in RISM–SCF method , 2007 .

[27]  William L. Jorgensen,et al.  Free Energies of Hydration from a Generalized Born Model and an All-Atom Force Field , 2004 .

[28]  Fumio Hirata,et al.  Potential of Mean Force between Two Molecular Ions in a Polar Molecular Solvent: A Study by the Three-Dimensional Reference Interaction Site Model , 1999 .

[29]  F. Hirata,et al.  Hydration structure and stability of Met-enkephalin studied by a three-dimensional reference interaction site model with a repulsive bridge correction and a thermodynamic perturbation method , 2000 .

[30]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[31]  S. Sakaki,et al.  New generation of the reference interaction site model self-consistent field method: introduction of spatial electron density distribution to the solvation theory. , 2007, The Journal of chemical physics.

[32]  David Chandler,et al.  Free energy functions in the extended RISM approximation , 1985 .

[33]  Fumio Hirata,et al.  Combination of molecular dynamics method and 3D‐RISM theory for conformational sampling of large flexible molecules in solution , 2008, J. Comput. Chem..

[34]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[35]  T. Truong,et al.  Thermochemistry of solvation: A self-consistent three-dimensional reference interaction site model approach , 2000 .

[36]  Maxim V Fedorov,et al.  Solvent effects and hydration of a tripeptide in sodium halide aqueous solutions: an in silico study. , 2007, Physical chemistry chemical physics : PCCP.

[37]  G. Morriss,et al.  Recent progress in the statistical mechanics of interaction site fluids , 1990 .

[38]  Gerhard Hummer,et al.  Multistate Gaussian Model for Electrostatic Solvation Free Energies , 1997 .

[39]  Gennady N. Chuev,et al.  A quasilinear RISM approach for the computation of solvation free energy of ionic species , 2006 .

[40]  Maxim V. Fedorov,et al.  Mathematik in den Naturwissenschaften Leipzig Low-Rank wavelet solver for the Ornstein-Zernike integral equation , 2005 .

[41]  Dennis R. Salahub,et al.  Optimization of Gaussian-type basis sets for local spin density functional calculations. Part I. Boron through neon, optimization technique and validation , 1992 .

[42]  F. Hirata,et al.  Application of the reference interaction site model theory to analysis on surface‐induced structure of water , 1996 .

[43]  P. Cummings,et al.  Exact asymptotic form of the site-site direct correlation function for rigid polar molecules , 1981 .

[44]  Fumio Hirata,et al.  An extended rism equation for molecular polar fluids , 1981 .

[45]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. I. Three-dimensional reference interaction site model approach , 2000 .

[46]  J. R. Pliego,et al.  Gibbs energy of solvation of organic ions in aqueous and dimethyl sulfoxide solutions , 2002 .

[47]  Ruth M. Lynden-Bell,et al.  From hydrophobic to hydrophilic behaviour: A simulation study of solvation entropy and free energy of simple solutes , 1997 .

[48]  G. Morriss,et al.  Recent Progress in the Statistical Mechanical Mechanics of Interaction Site Fluids , 2007 .

[49]  Fumio Hirata,et al.  Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model , 1999 .

[50]  D. Blankschtein,et al.  Liquid-state theory of hydrocarbon-water systems: application to methane, ethane, and propane , 1992 .

[51]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[52]  P. Rossky,et al.  A THREE-DIMENSIONAL REDUCTION OF THE ORNSTEIN-ZERNICKE EQUATION FOR MOLECULAR LIQUIDS , 1997 .

[53]  S. Sakaki,et al.  An integral equation theory for 3D solvation structure: A new procedure free from 3D Fourier transform , 2006 .

[54]  R. Gurney.  Ionic processes in solution , 1953 .

[55]  B. Roux,et al.  Solvation Free Energy of Polar and Nonpolar Molecules in Water: An Extended Interaction Site Integral Equation Theory in Three Dimensions , 2000 .

[56]  Fumio Hirata,et al.  Solution of three‐dimensional reference interaction site model and hypernetted chain equations for simple point charge water by modified method of direct inversion in iterative subspace , 1999 .

[57]  S. Sakaki,et al.  A new method to reconstruct three-dimensional spatial distribution function from radial distribution function in solvation structure. , 2005, The Journal of chemical physics.

[58]  O. Samoilov.  A new approach to the study of hydration of ions in aqueous solutions , 1957 .

[59]  Fumio Hirata,et al.  Analytical energy gradient for the reference interaction site model multiconfigurational self‐consistent‐field method: Application to 1,2‐difluoroethylene in aqueous solution , 1996 .

[60]  M. Born.  Volumen und Hydratationswärme der Ionen , 1920 .

[61]  B. Montgomery Pettitt,et al.  Application of an extended RISM equation to dipolar and quadrupolar fluids , 1982 .

[62]  Fumio Hirata,et al.  Hydration free energy of hydrophobic solutes studied by a reference interaction site model with a repulsive bridge correction and a thermodynamic perturbation method , 2000 .

[63]  P. Rossky,et al.  The coupling of long and short range correlations in ISM liquids , 1983 .

[64]  Jiali Gao,et al.  Hybrid Quantum and Molecular Mechanical Simulations: An Alternative Avenue to Solvent Effects in Organic Chemistry , 1996 .

[65]  Gerhard Hummer,et al.  Free Energy of Ionic Hydration , 1996 .

[66]  A. Kornyshev,et al.  Unravelling the solvent response to neutral and charged solutes , 2007 .

[67]  L. Blum,et al.  Invariant Expansion for Two‐Body Correlations: Thermodynamic Functions, Scattering, and the Ornstein—Zernike Equation , 1972 .

[68]  R. Friesner,et al.  Efficient recursive implementation of the modified Broyden method and the direct inversion in the iterative subspace method: Acceleration of self-consistent calculations , 1998 .

[69]  Fumio Hirata,et al.  Reference interaction site model self-consistent field study for solvation effect on carbonyl compounds in aqueous solution , 1994 .

[70]  Gennady N. Chuev,et al.  Comparative Study of Electrostatic Solvent Response by RISM and PCM Methods , 2007 .