Mixed Direct-Iterative Methods for Boundary Integral Formulations of Dielectric Solvation Models
S. A. Corcelli, J. D. Kress, L. R. Pratt, G. J. Tawa
Abstract
This paper describes a mixed direct-iterative method for boundary integral formulations of dielectric solvation models. We give an example for which a direct solution at thermal accuracy is nontrivial and for which Gauss-Seidel iteration diverges in rare but reproducible cases. This difficulty is analyzed by obtaining the eigenvalues and the spectral radius of the iteration matrix. This establishes that the nonconvergence is due to inaccuracies of the asymptotic approximations for the matrix elements for accidentally close boundary element pairs on different spheres. This difficulty is cured by checking for boundary element pairs closer than the typical spatial extent of the boundary elements and for those pairs performing an `in-line' Monte Carlo integration to evaluate the required matrix elements. This difficulty are not expected and have not been observed when only a direct solution is sought. Finally, we give an example application of these methods to deprotonation of monosilicic acid in water.
Create a lesson
Related papers
A Computational Study of Thirteen-atom Ar-Kr Cluster Heat Capacities
Don D. Frantz
A Variational Approach to the Structure and Thermodynamics of Linear Polyelectrolytes with Coulomb and Screened Coulomb Interactions
B. Jönsson, C. Peterson, B. Söderberg
Applications of the Optimized Effective Potential Method of Density Functional Theory to Atomic and Molecular Systems
T. Grabo, E. K. U. Gross
Relation between a Screened Polyelectrolyte and a Field Theory
Bo Soderberg
ΔI=4 and ΔI=8 bifurcations in rotational bands of diatomic molecules
Dennis Bonatsos, C. Daskaloyannis, G. A. Lalazissis et al.
A Dynamical Theory of Electron Transfer: Crossover from Weak to Strong Electronic Coupling
Juergen T. Stockburger, C. H. Mak