Wavelet algorithm for solving integral equations of molecular liquids. A test for the reference interaction site model

A new efficient method is developed for solving integral equations based on the reference interaction site model (RISM) of molecular liquids. The method proposes the expansion of site–site correlation functions into the wavelet series and further calculations of the approximating coefficients. To solve the integral equations we have applied the hybrid scheme in which the coarse part of the solution is calculated by wavelets with the use of the Newton–Raphson procedure, while the fine part is evaluated by the direct iterations. The Coifman 2 basis set is employed for the wavelet treatment of the coarse solution. This wavelet basis set provides compact and accurate approximation of site–site correlation functions so that the number of basis functions and the amplitude of the fine part of solution decrease sufficiently with respect to those obtained by the conventional scheme. The efficiency of the method is tested by calculations of SPC/E model of water. The results indicated that the total CPU time to obtain solution by the proposed procedure reduces to 20% of that required for the conventional hybrid method. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1369–1377, 2004

[1]  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 .

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

[3]  Reinhold Schneider,et al.  Wavelet approximation of correlated wave functions. I. Basics , 2002 .

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

[5]  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 .

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

[7]  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 .

[8]  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 .

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

[10]  T. Ichiye,et al.  Accurate integral equation theory for the central force model of liquid water and ionic solutions , 1988 .

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

[12]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

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

[14]  Stanislav Labík,et al.  A rapidly convergent method of solving the OZ equation , 1985 .

[15]  R. Coifman,et al.  Fast wavelet transforms and numerical algorithms I , 1991 .

[16]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

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

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

[19]  P. A. Monson.  Numerical solution of the RISM equations for the site-site 12-6 potential , 1982 .

[20]  W. M. Carey,et al.  Digital spectral analysis: with applications , 1986 .

[21]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. II. Solvation structure from the three-dimensional reference interaction site model approach, and comparison with simulations , 2000 .

[22]  T. Arias.  Multiresolution analysis of electronic structure: semicardinal and wavelet bases , 1998, cond-mat/9805262.

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

[24]  Johnson,et al.  Modified Broyden's method for accelerating convergence in self-consistent calculations. , 1988, Physical review. B, Condensed matter.

[25]  P. Pulay.  Improved SCF convergence acceleration , 1982 .

[26]  M. Fedorov,et al.  Wavelet treatment of radial distribution functions of solutes. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  P. Calmettes,et al.  Molecular solvent model for a cryptate solution in acetonitrile: A hypernetted chain study , 1994 .

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

[29]  K. Ng.  Hypernetted chain solutions for the classical one‐component plasma up to Γ=7000 , 1974 .

[30]  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 .

[31]  P. Rossky,et al.  Molecular conformational equilibria in liquids , 1986 .

[32]  Wolfgang Hackbusch,et al.  Wavelet approximation of correlated wave functions. II. Hyperbolic wavelets and adaptive approximation schemes , 2002 .

[33]  R. Fletcher.  Practical Methods of Optimization , 1988 .

[34]  S. Mallat.  A wavelet tour of signal processing , 1998 .

[35]  Yuko Okamoto,et al.  Calculation of hydration free energy for a solute with many atomic sites using the RISM theory: A robust and efficient algorithm , 1997 .

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

[37]  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 .

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

[39]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[40]  G. Zérah.  An efficient newton's method for the numerical solution of fluid integral equations , 1985 .

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

[42]  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 .

[43]  M. Gillan.  A new method of solving the liquid structure integral equations , 1979 .

[44]  B. Montgomery Pettitt,et al.  Integral equation predictions of liquid state structure for waterlike intermolecular potentials , 1982 .

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

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

[47]  Sujit K. Ghosh,et al.  Essential Wavelets for Statistical Applications and Data Analysis , 2001, Technometrics.

[48]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

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

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