The Born-Lande Equation

The Born-Landé framework, despite being a century-old classical electrostatic model, remains the standard first-pass teaching and estimation tool wherever chemists need a quick, physically transparent lattice energy without running a full quantum-chemical calculation. In industrial inorganic chemistry, lattice energy estimates feed directly into predicting solubility trends — the balance between lattice energy (which must be overcome) and hydration or solvation energy (which compensates it) explains why, for example, all Group 1 halides are soluble in water while many Group 2 fluorides and carbonates are not, a trend rationalized entirely by the $z^+z^-$ and $r_0$ dependence built into this equation. In materials science, the same electrostatic reasoning underlies the design of solid electrolytes and battery cathode materials (e.g., layered oxides and phosphates for Li-ion and Na-ion batteries), where researchers estimate the lattice energy penalty for extracting or inserting an alkali-metal ion during charge and discharge cycles. In geochemistry and mineralogy, lattice energy calculations of this type — refined into more elaborate atomistic simulations that still use a Born-type short-range repulsion as their starting ansatz — are used to predict mineral stability, cation-ordering preferences in silicates, and the thermodynamics of high-pressure mantle phases. Explosives and propellant chemistry uses lattice energy trends to rationalize the sensitivity and stability of ionic energetic salts. In pharmaceutical science, the same electrostatic bookkeeping (adapted to organic ionic co-crystals and salts) helps predict which counter-ion will give a drug candidate the right solubility and crystallinity for formulation. And pedagogically, the Born-Landé equation remains the essential stepping stone that every student climbs before reaching the two natural refinements it foreshadows: the exponential-repulsion Born-Mayer Equation, which replaces the power-law $B/r^n$ with a more physically justified $Be^{-r/\rho}$ term, and the Kapustinskii Equation, which sidesteps the need to know the Madelung constant or even the exact crystal structure at all, making lattice-energy estimation possible for compounds whose structures have never been solved.

Ask the tutor: Why does the Born exponent get assigned per-ion from noble-gas-like electron configurations rather than calculated from first principles, and how does the Born-Mayer equation improve on this?

Picture a single Na⁺ ion sitting at the origin of an infinite crystal of rock-salt NaCl. It has six Cl⁻ nearest neighbors at distance $r_0$, each pulling it inward with Coulomb attraction. Beyond them lie twelve Na⁺ ions at distance $r_0\sqrt2$, pushing it away. Beyond those, eight Cl⁻ ions at $r_0\sqrt3$ pull again, and so on outward forever. The infinite alternating sum of attractions and repulsions converges — slowly, and only when the terms are grouped correctly — to a single dimensionless number: the Madelung constant $A = 1.7476$ for the rock-salt structure. This number is a purely geometric fact about the arrangement of point charges; it says nothing about how strongly the ions are actually held together; that job falls to the prefactor built from charges and distance. Multiplying $A$ by the elementary Coulomb energy of one nearest-neighbor pair, $z^+z^-e^2/(4\pi\varepsilon_0 r_0)$, gives the full electrostatic (Madelung) energy of the crystal per ion pair. This term is purely attractive (energy becomes more negative as $r_0$ shrinks) and, taken alone, predicts that the crystal should collapse to zero separation — clearly unphysical. What stops the collapse is Pauli repulsion: as the filled electron shells of adjacent ions begin to overlap, the energy rises steeply, over a range of only a few tenths of an Ångström. Born and Landé modeled this repulsion phenomenologically as $B/r^n$, a term with no direct microscopic derivation at this level but chosen so that it turns on sharply at short range and is negligible at long range — exactly the qualitative behavior of orbital overlap repulsion. The genius of the Born-Landé construction is recognizing that you never need to know $B$ explicitly. The crystal automatically finds its own equilibrium distance $r_0$ — the value at which the net force is zero, i.e., where attraction and repulsion balance. This is precisely the same $r_0$ measured by X-ray crystallography. So the equilibrium condition $dU/dr\big|_{r_0} = 0$ becomes an equation that expresses $B$ in terms of $A$, the charges, $r_0$, and $n$. Substituting this expression for $B$ back into $U(r_0)$ purges the unknown constant entirely, leaving a lattice energy that depends only on quantities you can either measure or look up: charges, structure type (via $A$), interionic distance, and the Born exponent $n$, which itself is read off from the ions' electron configurations using simple empirical rules (Table below). Everything downstream — numbers in kJ/mol, comparisons with thermochemical cycles, trends across the periodic table — flows from this one elimination trick.

The electrostatic energy of one ion pair separated by $r$, embedded in a lattice whose long-range sum is captured by the Madelung constant $A$, is:

where: $N_A$ is Avogadro's number (putting the energy on a per-mole basis), $A$ is the dimensionless Madelung constant for the specific structure type, $z^+, z^-$ are the (signed, so $z^-<0$) ionic charge numbers, $e$ is the elementary charge, $\varepsilon_0$ is the vacuum permittivity, and $r$ is the center-to-center ion separation. Reading it: this term is negative (stabilizing) whenever the charges have opposite sign, and it becomes more negative — more stabilizing — as the ions get closer, diverging unphysically as $r\to0$. The short-range Born repulsion, summed over the lattice in the same way, is modeled as:

where: $B$ is an empirical positive constant absorbing the lattice geometry and the strength of orbital overlap repulsion (never measured directly), and $n$ is the Born exponent, a small integer (roughly 5 to 12) that controls how steeply the repulsion rises as clouds interpenetrate. Reading it: unlike the Coulomb term, this one grows without bound as $r\to0$ and vanishes rapidly as $r$ increases — it acts only over a short range, exactly like real electron-cloud overlap. The total lattice potential energy is the sum:

where: every symbol is as defined above. Reading it: this single-variable function of $r$ has exactly one minimum for $r>0$ (attraction dominates at large $r$, repulsion dominates at small $r$), and that minimum is the equilibrium state the real crystal adopts.

Example 1 — Foundational: lattice energy of NaCl from first principles. Rock-salt NaCl has $A = 1.74756$, $z^+ = +1$, $z^- = -1$, and the experimental Na–Cl distance $r_0 = 2.820\times10^{-10}\text{ m}$ (from X-ray diffraction). Both Na⁺ and Cl⁻ have the electron configuration of neon and argon respectively — noble-gas cores — so the Born exponent for this pair is taken as the mean of the individual assignments; standard tables give $n(\text{Na}^+)=7$ (Ne-like, using the empirical rule below) and $n(\text{Cl}^-)=9$ (Ar-like), averaging to $n = 8$. Plugging into the boxed equation, using $e = 1.602\times10^{-19}\,\text{C}$, $\varepsilon_0 = 8.854\times10^{-12}\,\text{C}^2\text{N}^{-1}\text{m}^{-2}$, $N_A = 6.022\times10^{23}\,\text{mol}^{-1}$:

Divide by $r_0$:

Multiply by $A z^+z^- = (1.74756)(1)(-1) = -1.74756$:

Apply the Born correction factor $(1-1/n) = (1-1/8) = 0.875$:

The experimental (Born-Haber-cycle-derived) lattice energy of NaCl is about $-787$ kJ/mol for the lattice enthalpy of dissociation convention used in many UK textbooks, but note carefully: different sourcing conventions and slightly different $r_0$/$n$ choices in the literature give values from about $-753$ to $-790$ kJ/mol for NaCl depending on exactly which experimental inputs are used; the calculation above, using $r_0=2.82\,\text{Å}$ and $n=8$, is the standard textbook Born-Landé result and is typically quoted as approximately $-756$ to $-770$ kJ/mol in most sources — discrepancies of this size (a few percent) are normal and reflect the sensitivity of the formula to the precise $r_0$ and $n$ used, not an error in the derivation. (Careful recomputation with $r_0 = 2.814\times10^{-10}$ m as sometimes tabulated, and rounding conventions, reproduces the commonly cited $-756$ kJ/mol figure closely.) The key pedagogical point is the method, and that the predicted magnitude — hundreds of kJ/mol, comparable to strong covalent bond energies — correctly captures why ionic solids are so hard to melt or dissolve.

Example 2 — Applied: comparing MgO and NaCl to explain relative melting points. MgO adopts the same rock-salt structure as NaCl, so $A = 1.74756$ for both. But MgO has doubly charged ions: $z^+=+2$, $z^-=-2$, and a shorter equilibrium distance, $r_0(\text{MgO}) = 2.10\times10^{-10}\,\text{m}$ (versus $2.82\times10^{-10}\,\text{m}$ for NaCl). Both Mg²⁺ and O²⁻ are also Ne-like, so take $n=7$. The Coulomb prefactor for MgO:

Multiply by $A|z^+z^-| = 1.74756 \times 4 = 6.9902$:

Apply $(1-1/n) = (1-1/7) = 0.8571$:

Compare directly with the NaCl estimate above (roughly $-1250$ kJ/mol on the same convention). MgO's lattice energy is over five times larger in magnitude, driven almost entirely by the $z^+z^-$ factor of $4$ (from $(2)(2)$ versus $(1)(1)$) partly offset, but not nearly enough, by the smaller distance boosting it further. This is the electrostatic reason MgO melts at 2852 °C while NaCl melts at only 801 °C: doubling both ionic charges quadruples the leading Coulomb term, and the accompanying lattice contraction (smaller, more highly charged ions pack closer) only amplifies the effect rather than compensating it.

Example 3 — Exam/competition level: back-calculating the Born exponent from a measured lattice energy. Suppose a hypothetical 1:1 ionic solid MX crystallizes in the rock-salt structure ($A=1.74756$), has a measured (Born-Haber) lattice energy of $U_0 = -693.0$ kJ/mol, and an X-ray-determined interionic distance $r_0 = 3.15\times10^{-10}$ m. Find the Born exponent $n$ implied by these data, and decide whether it is chemically reasonable for two singly-charged, noble-gas-like ions. First compute the pure-Coulomb (uncorrected) term:

Multiply by $A z^+z^- = 1.74756\times(-1) = -1.74756$ (charges $+1,-1$):

This is the value the equation would give if $n\to\infty$ (no repulsion correction). The actual measured lattice energy is only $-693.0$ kJ/mol in magnitude, so the correction factor must be:

Solve for $n$:

This is a striking result: $n\approx2$ is far below the physically sensible range of 5–12 tabulated for real noble-gas-configured ions (see the table below). A value this low signals that the simple ionic point-charge model is failing badly for this hypothetical solid — the large discrepancy between the "pure Coulomb" estimate and the experimental value is too big to be explained by ordinary Born repulsion alone. In a real exam or research context, this is exactly the diagnostic used to flag significant covalent character (as happens for AgI, CuCl, and similar "soft-soft" ion pairs) or an incorrectly assumed structure type/coordination number, prompting a switch to a more complete treatment such as the Born-Mayer Equation with an explicitly exponential repulsion, or acceptance that the ionic model itself has broken down for this pair.