3D RISM theory with fast reciprocal-space electrostatics
The calculation of electrostatic solute-solvent interactions in 3D RISM (“three-dimensional reference interaction site model”) integral equation theory is recast in a form that allows for a computational treatment analogous to the “particle-mesh Ewald” formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems.
- OSTI ID:
- 22415524
- Journal Information:
- Journal of Chemical Physics, Vol. 142, Issue 11; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
Efficient implementation of the three-dimensional reference interaction site model method in the fragment molecular orbital method
Calculating the binding free energies of charged species based on explicit-solvent simulations employing lattice-sum methods: An accurate correction scheme for electrostatic finite-size effects