The Madelung Problem of Finite Crystals

TL;DR

This work tackles the Madelung problem for finite crystals by decomposing the pairwise electrostatic interaction into a bulk periodic term $\\nu_{pbc}$, a shape-dependent boundary term $\\nu_b$, and a leading finite-size correction $\\nu_{\\rm corr}$. It derives closed-form expressions for $\\nu_b$ and $\\nu_{\\rm corr}$ for orthogonal lattices, with the cubic case enabling $p^{-2}$ scaling for the finite-size term and $p^{-4}$ accuracy when both corrections are included. The authors present a rapidly convergent direct-summation scheme that remains renormalization-free and achieves nine-digit precision for modest crystal sizes (e.g., $p=60$), outperforming several traditional methods for common ionic crystals such as NaCl, ZnS, CaF$_2$, and CaTiO$_3$. The approach clarifies the distinct roles of boundary and finite-size effects, provides an additive framework for constructing Madelung constants via symmetry-inequivalent ion pairs, and connects to both dipole/quadrupole invariants and Clifford-based methods, enabling practical, high-accuracy calculations from small finite systems.

Abstract

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 ${\mathbf 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^4-40(x^4+y^4+z^4)]/[9\sqrt{3} (2p+1)^2]$ 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^3$ unit cells), enabling hands-on calculations of Madelung constants for a wide range of ionic crystals.

Figures (1)

• Figure 1: Two-dimensional cross-sectional view of two particles with unit charges ($\pm 1$) separated by vector ${\mathbf r}=(x,y,z)$ and their periodic images under standard periodic boundary conditions (left) and in the CS method (right). For the conventional direct sum, the potential at the reference particle (black dot) is obtained by summing the Coulomb interactions $1/\left| {\mathbf r}+{\mathbf n} \right|$ from all periodic images, as described by Eq. \ref{['eq:pw']}. In the CS method ($K^3$ unit cells with unit side length), the conventional Euclidean distance $\left| {\mathbf r} + {\mathbf n}\right|$ is replaced by the renormalized distance $d(\mathbf{r}, \mathbf{n})$ (Eqs. \ref{['eq:rd']} and \ref{['eq:cspw']}).