Can dielectric spheres accurately model atomic-scale interactions?

We calculate the polarization portion of electrostatic interactions at the atomic scale using quantum-mechanical methods such as density functional theories (DFT) and the coupled cluster approach, and using classical methods such as a surface charge method and a polarizable force field. The agreement among various methods is investigated. Using the coupled clusters method CCSD(T) with large basis sets as the reference, we find that for systems comprising two to six atoms and ions in S-states the classical surface charge method performs much better than commonly used DFT methods with moderate basis sets such as B3LYP/6-31G(d,p). The remarkable performance of the classical approach comes as a surprise. The present results indicate that the use of a rigorous formalism of classical electrostatics can be better justified for determining molecular interactions at intermediate distances than some of the widely used methods of quantum chemistry.