Green's function method for elimination of the spurious multipole interaction in the surface/interface slab model

When modeling a surface/interface using a slab geometry that imposes the periodic boundary condition, there are spurious dipole and higher-order multipole interactions with the image slabs. Here, we show that the effective screening medium (ESM) method [M. Otani and O. Sugino, Phys. Rev. B 73, 115407 (2006)] can deal with such electrostatic issues rigorously in a slab calculation. After benchmark calculations for a hypothetical water layer, we show that there is an error in the total energy caused by the multipole interaction, but this can be treated efficiently using the ESM method. The accurate treatment of the electrostatic interaction is essential for performing large-scale surface/interface calculations.