 Electrostatic model of ionic binding

For two point electric charges in vacuum, Coulomb's law gives the interaction potential energywith zero at infinite separation. Here $r_{ij}$ is the distance between them and ${ϵ}_0$ is the vacuum permittivity. Opposite signs give negative potential energy and attraction; equal signs give positive potential energy and repulsion. For a finite collection of distinct point electric charges, sum over $i<j$ to count each pair once. A bulk ionic crystal requires the extended lattice sum, including both signs, rather than multiplying one attractive pair by the coordination number.^[ref]

At short distances, overlapping electron density makes the simple point-charge picture inadequate. Quantum mechanics, including the Pauli exclusion principle, matters to the steep repulsive contribution. Pure attraction proportional to $-1/r$ would collapse as $r\to 0$; it cannot set an equilibrium spacing on its own. A useful illustrative model isHere $r$ scales an assumed fixed arrangement, $A$ packages its attractive contribution, and $B/r^n$ represents short-range repulsion. The constants have units energy times length and energy times length to the power $n$, respectively. This is a model, not a fitted description of a named salt.^[ref]

Differentiating gives $r_0=(nB/A{)}^{1/(n-1)}$ and $U^{\prime \prime }(r_0)=(n-1)A/r_0^3>0$, so the stationary spacing is a stable minimum. For the diagram's stipulated $n=9$, put $x=r/r_0$ and $E_0=A(1-1/n)/r_0>0$. ThenThe dimensionless curves show how competing interactions select a spacing. Their energy scale and exponent are synthetic; neither is a measured lattice energy or a melting prediction.