| Abstract: | ![]()
We have developed an algorithm for solving the reference interaction site model (RISM) equations for water near a solute molecule with many atomic sites (interaction sites). It is a hybrid of the Newton–Raphson and Picard methods and is judiciously constructed. Various considerations are given so that the computer time can be saved as much as possible. The robustness and high efficiency of the algorithm has been demonstrated for calculating hydration free energies of Met-enkephalin (a peptide with 75 sites) with different conformations. The Jacobian matrix is treated as part of the input data, and it has been found that the same matrix can be used for a considerably large set of different conformations of the solute molecule. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1320–1326, 1997 |
