An improved iterative solution to solve the electrostatic problem in the polarizable continuum model

Abstract. We present a discrete iterative interpolation scheme (DIIS) to improve the convergence rate of electrostatic calculations in the polarizable continuum model (PCM) to describe solvent effects on molecular solutes. The electrostatic calculations may easily become the bottleneck of the calculation when the solute size is large. For large molecules iterative procedures turn out to be computationally more convenient than matrix inversion or closure methods. The DIIS scheme is compared here to another iterative procedure (DAMP) and to the biconjugate gradient (BCG) method. The comparisons show that DIIS leads to a sizeable saving of computational time for the C-PCM and IEF-PCM methods (average 40%) compared to DAMP, and more than 50% with respect to the BCG method.