The Coulomb potential at an interior ion in a finite crystal of size p is given by a linear superposition of contributions from displacement vectors r= (x, y, z) to its neighbors. This additive structure underlies universal relationships among Madelung constants and applies to both standard periodic boundary conditions and alternative Clifford supercells. Each pairwise contribution decomposes into three physically distinct components: a periodic bulk term, a quadratic boundary term, and a finite-size correction whose leading order term is 24r⁴-40 (x⁴+y⁴+z⁴) /93 (2p+1) ² for cubic crystals with unit lattice constant. Combining this decomposition with linear superposition yields a rapidly convergent direct-summation scheme, accurate even at p=1 (3³ unit cells), enabling hands-on calculations of Madelung constants for a wide range of ionic crystals.
Zhao et al. (Thu,) studied this question.