Multigrid solver for the reference interaction site model of molecular liquids theory

In this article, we propose a new multigrid‐based algorithm for solving integral equations of the reference interactions site model (RISM). We also investigate the relationship between the parameters of the algorithm and the numerical accuracy of the hydration free energy calculations by RISM. For this purpose, we analyzed the performance of the method for several numerical tests with polar and nonpolar compounds. The results of this analysis provide some guidelines for choosing an optimal set of parameters to minimize computational expenses. We compared the performance of the proposed multigrid‐based method with the one‐grid Picard iteration and nested Picard iteration methods. We show that the proposed method is over 30 times faster than the one‐grid iteration method, and in the high accuracy regime, it is almost seven times faster than the nested Picard iteration method. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011

[1]  Maxim V Fedorov,et al.  An accurate prediction of hydration free energies by combination of molecular integral equations theory with structural descriptors. , 2010, The journal of physical chemistry. B.

[2]  Roberto Cammi,et al.  Continuum Solvation Models in Chemical Physics , 2007 .

[3]  M. Kinoshita.  Structure of aqueous electrolyte solutions near a hydrophobic surface , 2007 .

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

[5]  K. Dill,et al.  Orientation-dependent integral equation theory for a two-dimensional model of water , 2003 .

[6]  H. Krienke,et al.  Improvements of DRISM calculations: symmetry reduction and hybrid algorithms. , 2008, Physical chemistry chemical physics : PCCP.

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

[8]  M. Fedorov,et al.  Wavelet method for solving integral equations of simple liquids , 2005 .

[9]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

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

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

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

[13]  J. Goodman,et al.  To switch or not to switch: the effects of potassium and sodium ions on alpha-poly-L-glutamate conformations in aqueous solutions. , 2009, Journal of the American Chemical Society.

[14]  L. Scales,et al.  Efficient solution of liquid state integral equations using the Newton-GMRES algorithm , 1999 .

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

[16]  Wolfgang Hackbusch,et al.  Multi-grid methods and applications , 1985, Springer series in computational mathematics.

[17]  L. Gendre,et al.  Classical density functional theory of solvation in molecular solvents: Angular grid implementation , 2009 .

[18]  W. Hackbusch,et al.  F ¨ Ur Mathematik in Den Naturwissenschaften Leipzig a Multigrid Solver for the Integral Equations of the Theory of Liquids , 2022 .

[19]  Sauer,et al.  Multicanonical multigrid Monte Carlo method. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[20]  Galli,et al.  Real-space adaptive-coordinate electronic-structure calculations. , 1995, Physical review. B, Condensed matter.

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

[22]  Maxim V Fedorov,et al.  The effect of sodium chloride on poly-l-glutamate conformation. , 2009, Chemical communications.

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

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

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

[26]  Shigeyoshi Sakaki,et al.  A highly parallelizable integral equation theory for three dimensional solvent distribution function: application to biomolecules. , 2009, The Journal of chemical physics.

[27]  J. Guthrie,et al.  A blind challenge for computational solvation free energies: introduction and overview. , 2009, The journal of physical chemistry. B.

[28]  Maxim V Fedorov,et al.  Accurate calculations of the hydration free energies of druglike molecules using the reference interaction site model. , 2010, The Journal of chemical physics.

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

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

[31]  C. Kelley,et al.  A fast solver for the Ornstein-Zernike equations , 2004 .

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

[33]  Seiichiro Ten-no,et al.  Free energy of solvation for the reference interaction site model: Critical comparison of expressions , 2001 .

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

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

[36]  T. Torsti,et al.  Multigrid method for electronic structure calculations , 2001, cond-mat/0101233.

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

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

[39]  R. Gebauer,et al.  Density functional theory of solvation in a polar solvent: extracting the functional from homogeneous solvent simulations. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[40]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

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

[42]  You-nian Wang,et al.  Implicit and electrostatic particle-in-cell/Monte Carlo model in two-dimensional and axisymmetric geometry: I. Analysis of numerical techniques , 2010, 1001.2969.

[43]  M. Kinoshita,et al.  Entropic insertion of a big sphere into a cylindrical vessel , 2010 .

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

[45]  Shigeyoshi Sakaki,et al.  An analysis of 3D solvation structure in biomolecules: application to coiled coil serine and bacteriorhodopsin. , 2010, The journal of physical chemistry. B.

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

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

[48]  Herbert H. H. Homeier,et al.  Iterative solution of the Ornstein-Zernike equation with various closures using vector extrapolation , 1995, chem-ph/9509001.

[49]  S. Sakaki,et al.  Analytical energy gradient for reference interaction site model self-consistent field explicitly including spatial electron density distribution. , 2009, The Journal of chemical physics.

[50]  平田 文男.  Molecular theory of solvation , 2003 .

[51]  H. Fenniri,et al.  Structural water drives self-assembly of organic rosette nanotubes and holds host atoms in the channel. , 2010, Chemphyschem : a European journal of chemical physics and physical chemistry.

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

[53]  E. Meeron.  Erratum : Series Expansion of Distribution Functions in Multicomponent Fluid Systems , 1957 .

[54]  J. Lebowitz,et al.  Mean Spherical Model for Lattice Gases with Extended Hard Cores and Continuum Fluids , 1966 .