The Born-Mayer and Kapustinskii Equations

Take a compound that does not exist as a "normal" salt in your cupboard. Say hypothetical NaCl₂, or a wildly unstable oxide. Ask: how much lattice energy would it release if it did form? You cannot look up its crystal structure in a database. It has never been made. Yet chemists routinely predict such lattice energies to within a few percent. They use nothing but the ionic charges and a handful of tabulated "effective" ionic radii. That predictive power comes from two closely related refinements of the The Born-Lande Equation. The Born-Mayer equation replaces a crude power-law repulsion with a more physically honest exponential form. The Kapustinskii equation strips away the need for crystal-structure data entirely. Lattice energy is the energy released when gaseous ions condense into an ionic solid. It is the single number that ties together melting points, solubilities, thermal stabilities, and the very existence of compounds. The Born-Landé treatment already captures the main physics. Coulombic attraction between ions stabilizes the lattice. Short-range Pauli repulsion between overlapping electron clouds prevents collapse. But Born-Landé's $1/r^n$ repulsion term is a mathematical convenience with a shaky physical basis. The Born-Exponential (Born-Mayer) form fixes this by using an exponential repulsion. This mirrors how overlapping wavefunctions actually decay. Kapustinskii then made an astonishing observation. If you rescale any lattice sum to a fictitious "generalized" structure, the real structure almost doesn't matter. Madelung constants across many different lattice types turn out to be nearly proportional to the number of ions per formula unit. This lets you compute a lattice energy from tabulated thermochemical radii alone. You never need to know the crystal structure. The big picture is this: lattice energy calculations sit at the crossroads of electrostatics, quantum mechanical repulsion, and empirical bookkeeping. The Born-Landé equation gets the conceptual physics right. But it needs a compressibility-derived exponent $n$ and a real Madelung constant $A$ from a known structure. The Born-Mayer equation improves the repulsion physics with an exponential term. This term is calibrated by a single universal length parameter. The Kapustinskii equation trades some rigor for universality. It lets you estimate lattice energies for compounds whose structures are unknown, poorly characterized, or purely hypothetical. This is a tool indispensable for predicting the stability of exotic oxidation states. It also applies to noble-gas compounds and low-lattice-energy salts of large complex anions. This article builds both equations from the ground up. It covers the exponential repulsion model and its physical justification. It covers the algebra connecting Born-Landé to Born-Mayer via matched curvature/compressibility conditions. It covers the Kapustinskii generalization with its universal Madelung-constant-per-ion approximation. Three worked examples show how the numbers are actually built up term by term. These are a textbook check against MgO, an applied estimate for a real but awkward salt, and a competition-style prediction for a hypothetical superoxide. This machinery is also exactly the ingredient needed to close the The Born-Haber Cycle. There, predicted lattice energies are compared against experimentally derived ones. This tests whether a solid is "normal" ionic or has hidden covalent character.

Key takeaways

The Born-Landé $1/r^n$ repulsion is a convenient fit, not a physical law; the Born-Mayer equation replaces it with $Be^{-r/\rho}$ , matching how overlapping orbitals actually decay exponentially.
Both equations share the same Coulombic term $-\dfrac{N_A A z_+ z_- e^2}{4\pi\varepsilon_0 r_0}$ ; they differ only in how the repulsion term is written and normalized.
The Born-Mayer form uses a single, nearly universal softness parameter $\rho \approx 30$ – $35$ pm instead of a structure-dependent Born exponent $n$ .
Kapustinskii's insight: Madelung constants $A$ for many common structures are almost proportional to the number of ions per formula unit $\nu$ , so $A/r_0$ can be estimated without knowing the actual structure.
The Kapustinskii equation uses tabulated thermochemical radii (extended even to polyatomic ions like NO₃⁻ or SO₄²⁻) so lattice energies can be predicted for compounds that have never been crystallographically characterized.
Both equations feed directly into the The Born-Haber Cycle, where comparing calculated and experimental lattice energies reveals covalent contamination of "ionic" bonding.
Key idea: replacing an ad hoc power-law repulsion with a physically motivated exponential (Born-Mayer), and then replacing structure-specific Madelung constants with a universal ions-per-formula-unit scaling (Kapustinskii), turns lattice-energy calculation from a structure-dependent curve-fit into a nearly structure-free predictive tool.

Where this lives today

Kapustinskii-type estimates remain the standard first move whenever a chemist needs a lattice energy for a compound lacking a solved crystal structure — a routine situation for many transient, high-energy-density, or exotic-oxidation-state materials. In inorganic synthesis, predicted lattice energies via Kapustinskii's equation are used to rationalize why certain noble-gas compounds (like XeF⁺ salts) or unusual low oxidation states are only stabilizable with specific large, low-charge-density counterions — a large, soft anion lowers the lattice-energy "reward" needed to compensate destabilizing ionization steps, and Kapustinskii's formula quantifies exactly how that reward scales with $r_++r_-$ . In battery and energy-storage chemistry, lattice-energy estimates for candidate cathode and solid-electrolyte materials — many of them mixed-anion or defect-rich phases without fully resolved structures — use these same semi-empirical sums as a fast screening tool before committing to expensive first-principles density functional theory calculations. In geochemistry and mineralogy, Kapustinskii-style estimates help predict the relative thermodynamic stability of mineral polymorphs and evaporite salts containing complex oxyanions (sulfates, phosphates, silicates) where thermochemical radii substitute for structures that may only exist as fine-grained or poorly crystalline natural samples. Finally, these equations remain a teaching bridge into the Born-Haber Cycle: comparing a Born-Mayer or Kapustinskii-predicted lattice energy against the experimentally derived value (obtained by closing a thermochemical cycle of sublimation, ionization, dissociation, and electron-affinity steps) is the standard diagnostic for detecting covalent character hiding inside a nominally "ionic" bond, as famously seen for silver halides and many transition-metal chalcogenides.

Ask the tutor: How would you use a Kapustinskii-predicted lattice energy inside a full Born-Haber cycle to test whether AgCl's bonding is "normal" ionic or has significant covalent character?

The physical picture

In any ionic crystal, every ion sits in an electrostatic potential well built from all the other ions in the lattice — nearest-neighbor counter-ions pulling it in, next-nearest like-charged ions pushing it away, and so on out to infinity. This infinite alternating sum converges (slowly, and only when handled carefully) to the Madelung constant $A$ , a pure number depending only on the geometric arrangement of ions. The electrostatic part of the lattice energy is exact once $A$ and the shortest cation-anion distance $r_0$ are known. What is not exact in the elementary treatment is the repulsion. As two ions approach, their filled electron shells begin to overlap, and the Pauli exclusion principle forces electrons into higher kinetic-energy states — an effect that rises extremely steeply once wavefunction overlap becomes significant, and is negligible at large separation. The Born-Landé equation mimics this using $B/r^n$ , a repulsion chosen because its exponent $n$ can be extracted from bulk compressibility data. But real overlap-driven repulsion is not a power law — it is the tail of overlapping exponential-decay wavefunctions, so it should itself look exponential: $B\,e^{-r/\rho}$ . Max Born and Joseph Mayer built exactly this into their 1932 refinement, giving the Born-Mayer equation. The parameter $\rho$ (with units of length, typically 30–35 pm) sets how "hard" or "soft" the repulsive wall is, and turns out to be close to universal across many rock-salt-type compounds, removing the need for a compound-specific Born exponent. Anton Kapustinskii's 1956 contribution operates one level up in abstraction. He noticed empirically that if you tabulate Madelung constants $A$ for a wide variety of real crystal structures (rock salt, fluorite, rutile, cesium chloride...) and divide each by the number of ions in the formula unit $\nu$ (e.g., $\nu = 2$ for NaCl, $\nu = 3$ for CaF₂), the ratio $A/\nu$ is nearly the same constant across all of them — roughly $0.85$ – $0.88$ . This means you never actually have to know the crystal structure of a salt to estimate its lattice energy to a useful accuracy: you only need the formula unit stoichiometry and the ionic charges. Combined with thermochemical radii — effective ionic radii back-calculated from known lattice energies, and extended by consistency to polyatomic ions that have no well-defined crystallographic "radius" at all — Kapustinskii's equation becomes a tool that works even for compounds whose structure has never been solved, or that do not exist yet.

Figure

From power-law to exponential repulsion: how the repulsive wall softens with a physically-motivated form

Figure

Kapustinskii's rescaling: real Madelung constants A collapse onto a near-universal line when divided by ions per formula unit ν

The core ideas

The Born-Landé lattice energy for one mole of an ionic solid with charges $z_+$ , $z_-$ (as unitless multiples of $e$ ) and equilibrium separation $r_0$ is

U_{\text{BL}} = -\frac{N_A A z_+ z_- e^2}{4\pi\varepsilon_0 r_0}\left(1 - \frac{1}{n}\right)

where: $N_A$ is Avogadro's number, $A$ is the Madelung constant for the specific structure, $e$ is the elementary charge, $\varepsilon_0$ is vacuum permittivity, $r_0$ is the equilibrium (shortest) cation-anion distance, and $n$ is the Born exponent obtained from compressibility measurements. The Born-Mayer equation replaces the $(1 - 1/n)$ repulsion correction with a term built from the softness parameter $\rho$ :

U_{\text{BM}} = -\frac{N_A A z_+ z_- e^2}{4\pi\varepsilon_0 r_0}\left(1 - \frac{\rho}{r_0}\right)

where: $\rho$ is a repulsion range parameter (typically $\rho \approx 30$ – $35$ pm, nearly transferable between compounds), and all other symbols are as before. Reading this equation: the Coulombic prefactor is identical to Born-Landé's, but the correction factor $(1-\rho/r_0)$ now comes from an exponential repulsive term whose decay length $\rho$ is largely material-independent, so a single universal value can be used across many salts without separately measuring a compressibility-derived $n$ for each one. The Kapustinskii equation goes a step further, folding the (now-unnecessary) crystal-structure-specific Madelung constant into a universal proportionality against the number of ions per formula unit, $\nu$ :

U_{\text{K}} = \frac{1.2025\times10^{5}\;\nu\,z_+z_-}{r_+ + r_-}\left(1 - \frac{34.5\text{ pm}}{r_+ + r_-}\right)\ \text{kJ mol}^{-1}

where: $\nu$ is the number of ions per formula unit (e.g. 2 for NaCl, 3 for CaF₂ or Na₂O), $z_+, z_-$ are the ionic charge numbers (as positive magnitudes), $r_+ + r_-$ is the sum of thermochemical radii in picometers, and the constants $1.2025\times10^5$ kJ·pm·mol⁻¹ and $34.5$ pm absorb $N_Ae^2/4\pi\varepsilon_0$ together with the near-universal Madelung-constant-per-ion ratio and repulsion range. Reading this: everything a real crystal-structure calculation would need — the actual Madelung constant, the actual repulsion exponent — has been replaced by two universal constants and a simple stoichiometric count, at the cost of a few percent accuracy.

Derivation AFrom Born-Landé repulsion to Born-Mayer exponential repulsion

Goal: show how matching an exponential repulsion term to the same equilibrium condition used in Born-Landé (zero net force at $r=r_0$ ) produces the Born-Mayer correction factor $(1-\rho/r_0)$ , and relate $\rho$ to the Born exponent $n$ .

Step 1 — write the total lattice energy per mole with a general repulsion term. Let the pairwise-summed lattice energy be

U(r) = -\frac{N_A A z_+ z_- e^2}{4\pi\varepsilon_0 r} + N_A B f(r)

where $f(r)$ is a repulsion function common to all ion pairs (already summed with appropriate geometric weighting into the constant $B$ ), and $r$ is the nearest-neighbor cation-anion distance. In Born-Landé, $f(r) = r^{-n}$ ; in Born-Mayer, $f(r) = e^{-r/\rho}$ .

Step 2 — impose the equilibrium condition. At the observed equilibrium separation $r=r_0$ , the net force is zero, i.e. $dU/dr = 0$ :

\frac{dU}{dr}\bigg|_{r_0} = \frac{N_A A z_+ z_- e^2}{4\pi\varepsilon_0 r_0^{2}} + N_A B f'(r_0) = 0

This fixes $B$ in terms of $A$ , $z_+z_-$ , $r_0$ , and the derivative of $f$ .

Step 3 — solve for B with the exponential form. For $f(r) = e^{-r/\rho}$ , $f'(r) = -\tfrac{1}{\rho}e^{-r/\rho}$ , so

\frac{A z_+ z_- e^2}{4\pi\varepsilon_0 r_0^{2}} = \frac{B}{\rho}e^{-r_0/\rho}

giving

B e^{-r_0/\rho} = \frac{\rho\, A z_+z_- e^2}{4\pi\varepsilon_0 r_0^{2}}

Note that $Be^{-r_0/\rho}$ is exactly the repulsive energy term $N_A B f(r_0)$ evaluated at equilibrium (divided by $N_A$ ), so this equation says the repulsive energy at $r_0$ equals $(\rho/r_0)$ times the magnitude of the Coulombic energy at $r_0$ .

Step 4 — substitute back into the total energy.

U(r_0) = -\frac{N_A A z_+z_- e^2}{4\pi\varepsilon_0 r_0} + N_A\left[\frac{\rho\, A z_+z_- e^2}{4\pi\varepsilon_0 r_0^{2}}\right]

Factor out the Coulombic prefactor:

U(r_0) = -\frac{N_A A z_+z_- e^2}{4\pi\varepsilon_0 r_0}\left(1 - \frac{\rho}{r_0}\right)

This is exactly the Born-Mayer equation. The repulsion contributes a fractional correction $\rho/r_0$ , typically of order $0.1$ , because $\rho \ll r_0$ — the repulsive wall is very short-ranged compared to the equilibrium spacing.

Step 5 — connect $\rho$ back to the Born exponent $n$ (consistency check). Repeating Steps 2–4 with $f(r) = r^{-n}$ instead reproduces the Born-Landé correction $(1-1/n)$ . Equating the two correction factors at the same equilibrium distance requires

\frac{\rho}{r_0} = \frac{1}{n} \quad\Longrightarrow\quad \rho = \frac{r_0}{n}

For typical alkali halides, $n \approx 8$ – $10$ and $r_0 \approx 280$ pm, giving $\rho \approx 28$ – $35$ pm — consistent with the "universal" range quoted for Born-Mayer. This is why a single transferable $\rho$ can replace a compound-specific $n$ : physically, $\rho$ is a genuine decay length of overlapping orbitals, whereas $n$ is only a disguised, structure-dependent way of expressing the same decay through a power law's local slope.

Go deeper · Derivation B · 2 minKapustinskii's rescaling: from real Madelung constants to a universal ions-per-formula-unit proportionality

Derivation BKapustinskii's rescaling: from real Madelung constants to a universal ions-per-formula-unit proportionality

Goal: derive how replacing the structure-specific Madelung constant $A$ with the empirical proportionality $A \approx 0.874\,\nu$ (with $\nu$ the number of ions per formula unit) turns the Born-Mayer equation into the structure-free Kapustinskii equation, and pin down the numerical constants.

Step 1 — start from the Born-Mayer lattice energy in convenient units.

U_{\text{BM}} = -\frac{N_A A z_+z_-e^2}{4\pi\varepsilon_0 r_0}\left(1-\frac{\rho}{r_0}\right)

Evaluate the physical constants: $N_A e^2/4\pi\varepsilon_0 = 1.3893\times 10^{-4}$ J·m·mol⁻¹ $= 1.3893\times10^{5}$ kJ·pm·mol⁻¹ (using $r_0$ in picometers gives $U$ directly in kJ mol⁻¹ once the sign convention for $z_+z_-$ , both taken as positive magnitudes with charges of opposite sign understood, is fixed).

Step 2 — replace $A$ using the Kapustinskii empirical trend. Tabulating Madelung constants for a broad range of real ionic structures against $\nu$ (rock salt $A=1.748,\ \nu=2$ ; CsCl $A=1.763,\ \nu=2$ ; fluorite $A=2.519,\ \nu=3$ ; rutile $A=2.408,\ \nu=3$ ; corundum-type oxides, $\nu=5$ , etc.) shows that the ratio $A/\nu$ clusters closely around

\frac{A}{\nu} \approx 0.874

across structure types differing by coordination number, geometry, and stoichiometry. Substituting $A = 0.874\,\nu$ into the Coulombic prefactor:

\frac{N_A A z_+z_-e^2}{4\pi\varepsilon_0} = 0.874\,\nu\, z_+z_- \times 1.3893\times10^{5}\ \text{kJ pm mol}^{-1}

Step 3 — collapse the numerical prefactor.

0.874 \times 1.3893\times10^{5} = 1.2143\times10^{5}\ \text{kJ pm mol}^{-1}

Kapustinskii's originally tabulated constant (incorporating a slightly different empirical average and rounding of the repulsion-inclusive prefactor) is quoted as

1.2025\times10^{5}\ \text{kJ pm mol}^{-1}

the small difference from the value above reflects that the fitted average of $A/\nu$ and the fitted average repulsion correction were optimized jointly against a training set of known lattice energies rather than derived from a single idealized structure — this is an empirically calibrated equation, not a purely first-principles one.

Step 4 — replace $r_0$ with the sum of thermochemical radii. Since a hypothetical or poorly characterized compound has no measured $r_0$ , Kapustinskii uses the sum of tabulated thermochemical radii, $r_+ + r_-$ , back-derived so that when substituted into this same formula they reproduce known experimental lattice energies for well-characterized salts. This is a self-consistent bootstrapping: a handful of "anchor" compounds with reliable experimental lattice energies and known crystal radii calibrate the thermochemical radii of polyatomic or unusual ions, and those radii are then used predictively for everything else.

Step 5 — assemble the final form.

U_{\text{K}} = \frac{1.2025\times10^{5}\,\nu z_+z_-}{r_++r_-}\left(1-\frac{34.5\ \text{pm}}{r_++r_-}\right)\ \text{kJ mol}^{-1}

The repulsion range constant $34.5$ pm plays the role of $\rho$ from Born-Mayer, again taken as a universal average value rather than compound-specific. As a limiting-case sanity check, note that as $r_+ + r_- \to \infty$ (extremely large, diffuse ions), $U_{\text{K}} \to 0$ — infinitely large, soft ions produce vanishingly small lattice energy, exactly as physically expected, since Coulombic attraction weakens as $1/r$ and the repulsion correction term also vanishes.

Worked examples

Example 1 — Foundational: Born-Mayer lattice energy of MgO, checked against Born-Landé. MgO crystallizes in the rock-salt structure, $A = 1.7476$ , with $z_+ = 2$ , $z_- = 2$ (both taken as positive magnitudes in the formula, with the understanding that one charge is positive and one negative), and an experimental equilibrium interionic distance $r_0 = 210$ pm. Take the Born-Mayer softness parameter as $\rho = 30$ pm (a typical transferable value for oxide/dense-halide lattices). Coulombic prefactor:

\frac{N_A A z_+z_-e^2}{4\pi\varepsilon_0 r_0} = \frac{1.3893\times10^{5}\ \text{kJ pm mol}^{-1}\times 1.7476\times 4}{210\ \text{pm}}

= \frac{9.712\times10^{5}}{210} = 4625\ \text{kJ mol}^{-1}

Repulsion correction:

1-\frac{\rho}{r_0} = 1-\frac{30}{210} = 1-0.1429 = 0.8571

Lattice energy:

U_{\text{BM}} = -4625 \times 0.8571 = -3964\ \text{kJ mol}^{-1}

Compare with Born-Landé using $n=8$ (typical for O²⁻/Mg²⁺, both with neon-like or intermediate configurations, chosen for illustration):

1-\frac{1}{n} = 1-\frac{1}{8} = 0.875

U_{\text{BL}} = -4625\times 0.875 = -4047\ \text{kJ mol}^{-1}

The two methods agree to within about 2%, which is typical: both capture the same dominant Coulombic term, and the repulsion corrections (0.857 vs 0.875) are close numerically even though built from different physical assumptions. The experimental Born-Haber value for MgO lattice energy is close to $-3795$ kJ mol⁻¹, showing that both models slightly overestimate the magnitude here — a reminder that real MgO bonding has a modest covalent contribution beyond the pure ionic model, a discrepancy the The Born-Haber Cycle is specifically built to expose.

Example 2 — Applied: Kapustinskii estimate for calcium carbonate's lattice energy, using a polyatomic thermochemical radius. CaCO₃ has $z_+ = 2$ (Ca²⁺), $z_- = 2$ (CO₃²⁻), and $\nu = 2$ ions per formula unit (one Ca²⁺, one CO₃²⁻ — note $\nu$ counts ions, not atoms). The thermochemical radius of Ca²⁺ is $r_+ = 114$ pm (close to its ionic radius, since it is a simple monatomic cation), and the thermochemical radius of CO₃²⁻, derived from fitting known carbonate lattice energies, is approximately $r_- = 178$ pm. Sum of radii:

r_+ + r_- = 114 + 178 = 292\ \text{pm}

Numerator:

1.2025\times10^{5} \times \nu \times z_+z_- = 1.2025\times10^{5}\times 2 \times 4 = 9.620\times10^{5}

Coulombic-type term:

\frac{9.620\times10^{5}}{292} = 3295\ \text{kJ mol}^{-1}

Repulsion correction:

1-\frac{34.5}{292} = 1-0.1182 = 0.8818

Lattice energy:

U_{\text{K}} = -3295\times 0.8818 = -2906\ \text{kJ mol}^{-1}

This value (around $-2900$ kJ mol⁻¹) is close to the experimentally derived lattice energy of calcite, roughly $-2800$ kJ mol⁻¹, obtained independently via a full Born-Haber thermochemical cycle. This agreement is precisely why Kapustinskii's method is valuable: CO₃²⁻ has no meaningful crystallographic "ionic radius" (it is not spherical), yet a thermochemical radius fitted once from a reliable carbonate lattice energy transfers usefully to other carbonates, sulfates, nitrates, and similar polyatomic-anion salts.

Example 3 — Exam/competition level: predicting the lattice energy of a hypothetical compound, "NaO₂"-type superoxide vs. a genuinely unmade salt, "CaCl₃" (Ca³⁺, hypothetical). Consider hypothetically forming CaCl₃ with an implausible Ca³⁺ ion (this ion does not exist due to the enormous third ionization energy of calcium, but the lattice-energy side of a Born-Haber argument can still be built to show why it doesn't form). Take $z_+ = 3$ , $z_- = 1$ for three Cl⁻, so $\nu = 4$ ions per formula unit. Use thermochemical/ionic radii $r_+(\text{Ca}^{3+}) \approx 99$ pm (estimated by extrapolating the isoelectronic radius trend down from Ca²⁺'s 114 pm, since higher charge contracts the ion) and $r_-(\text{Cl}^-) = 181$ pm. Sum of radii:

r_+ + r_- = 99 + 181 = 280\ \text{pm}

Numerator:

1.2025\times10^{5}\times 4 \times 3 = 1.4430\times10^{6}

Coulombic-type term:

\frac{1.4430\times10^{6}}{280} = 5154\ \text{kJ mol}^{-1}

Repulsion correction:

1-\frac{34.5}{280} = 1-0.1232 = 0.8768

Estimated lattice energy:

U_{\text{K}} = -5154\times 0.8768 = -4519\ \text{kJ mol}^{-1}

This is a genuinely enormous lattice energy — nearly 40% larger in magnitude than real CaCl₂'s roughly $-2258$ kJ mol⁻¹ — because lattice energy scales strongly with $z_+z_-$ (here $z_+z_- = 3$ instead of $2$ ) and with $\nu$ . Yet CaCl₃ still does not exist, because the third ionization energy of calcium (to strip a core electron from an already noble-gas-configuration Ca²⁺, requiring roughly $4900$ kJ mol⁻¹) vastly exceeds any conceivable extra lattice-energy payoff. This is exactly the kind of comparison the full Born-Haber Cycle makes rigorous: a Kapustinskii-predicted lattice energy, however large, is only half the story — it must be weighed against the ionization, atomization, and electron-affinity terms before concluding whether a compound is thermodynamically favorable.

Limiting cases, scope & common traps

Regime / question	Born-Mayer behavior	Kapustinskii behavior
Structure known precisely	Most accurate; uses real $A$ , real $r_0$ , calibrated $\rho$	Slightly less accurate; uses averaged $A/\nu \approx 0.874$ instead of exact $A$
Structure unknown or hypothetical	Cannot be used (no real $A$ , $r_0$ )	Designed for exactly this case
Polyatomic ions (CO₃²⁻, NO₃⁻, SO₄²⁻)	Needs an effective radius anyway; awkward	Thrives — thermochemical radii built for this purpose
$\rho \to 0$ (infinitely hard repulsion)	Correction factor $\to 1$ ; reduces to pure Coulomb (unphysical, ions could collapse)	Same limiting behavior via $34.5\to0$
$r_++r_- \to \infty$	$U\to 0$ (weak, diffuse interaction)	$U\to 0$ ; same physical limit
High-charge, small-radius ions (e.g. Al³⁺, N³⁻)	Both formulas still qualitatively fine, but covalent character grows and pure ionic models increasingly overestimate $\|U\|$	Same caveat; even more pronounced since real structure corrections are also lost
Comparing to experiment (Born-Haber cycle)	Small discrepancies flag covalency	Larger, systematic discrepancies expected due to averaging step

The chemist behind it

Anton Fedorovich Kapustinskii (1906–1960) was a Soviet physical chemist working at the Karpov Institute in Moscow, specializing in thermochemistry and crystal energetics during a period when Soviet chemistry was rapidly building its own thermodynamic data tables somewhat independently of Western literature. His 1933–1956 work on lattice energies grew from a very practical problem: experimentalists needed lattice energies for compounds — including many double salts, hydrates, and complex anions relevant to Soviet industrial mineral processing and radiochemistry — for which no single crystal structure determination existed or was practical to obtain with 1930s–1950s X-ray techniques. Rather than wait for structural crystallography to catch up, Kapustinskii built a semi-empirical shortcut, systematically tabulating what came to be called thermochemical radii precisely so that any chemist with a periodic table and a slide rule could estimate a lattice energy for a salt no one had ever grown a crystal of. Max Born (1882–1970) and Joseph E. Mayer (1904–1983), working independently and earlier (1932), refined the repulsion term of the Born-Landé model using arguments grounded in the emerging quantum mechanics of overlapping atomic orbitals — Born having already co-developed the original ionic-crystal lattice-energy framework with Alfred Landé a decade earlier.

Origins · 1 minAnton Fedorovich Kapustinskii and how this was discovered

Anton Fedorovich Kapustinskii

1906–1960 · Soviet thermochemistry and lattice-energetics, Karpov Physico-Chemical Institute, Moscow

History

- 1918 — Born and Landé propose the original lattice-energy expression with power-law

1/r^n

repulsion.

- 1932 — Born and Mayer replace the power-law repulsion with an exponential term, giving the Born-Mayer equation.

- 1933 — Kapustinskii begins systematic thermochemical studies of lattice energetics across many ionic and complex-ion crystals.

- 1956 — Kapustinskii publishes the generalized equation using universal Madelung-constant-per-ion scaling and thermochemical radii.

- Mid-20th century onward — thermochemical radii extended to a wide range of polyatomic anions and cations, cementing the equation's use in predicting lattice energies of uncharacterized or unstable salts.

Why isn't the Born-Landé equation simply "good enough" — why introduce Born-Mayer at all?

Born-Landé's power-law repulsion $B/r^n$ works numerically because $n$ is fitted from experimental compressibility data for each specific compound, so it reproduces the equilibrium distance and curvature correctly by construction. But the physical justification is weak: real repulsion arises from overlapping filled-shell wavefunctions, which decay exponentially with distance, not as an inverse power. Born-Mayer's exponential form matches this physical picture and, crucially, uses a softness parameter $\rho$ that is nearly transferable between compounds (roughly 30–35 pm for many oxides and halides), so you don't need a separate compressibility measurement for every new material. It is a more physically grounded, and often more transferable, refinement rather than a strictly more accurate one for any single well-characterized compound.

Is the Kapustinskii equation less accurate than Born-Landé or Born-Mayer?

For compounds with well-known crystal structures, yes — typically Kapustinskii's estimate differs from a proper Born-Mayer calculation by a few percent, because it replaces the exact Madelung constant with an averaged $A/\nu \approx 0.874$ proportionality. But this is not really a fair comparison, because Kapustinskii's equation is designed for the opposite situation: compounds whose crystal structure is unknown, ill-defined, or purely hypothetical, where Born-Landé and Born-Mayer simply cannot be applied at all. In that regime, Kapustinskii's few-percent uncertainty is far preferable to having no estimate whatsoever.

What exactly is a "thermochemical radius" and how is it obtained?

A thermochemical radius is an effective ionic radius derived by working the Kapustinskii equation backward: take a compound with a reliably known experimental lattice energy (from a closed Born-Haber cycle) and a known radius for one of its ions, then solve for the radius of the other ion that makes the equation reproduce the correct lattice energy. Once fitted from one or more "anchor" compounds, that radius is then reused predictively for the same ion in other, less well-characterized salts. This is especially valuable for polyatomic ions like NO₃⁻, ClO₄⁻, or SO₄²⁻, which have no meaningful single crystallographic radius because they are not spherical.

Why does the Kapustinskii equation not require the crystal structure at all?

Because its key empirical finding is that the Madelung constant $A$ , when divided by the number of ions per formula unit $\nu$ , is nearly the same (about 0.87) across a wide range of very different structure types — rock salt, cesium chloride, fluorite, rutile, and others. This means the detailed geometric arrangement of ions barely affects the electrostatic energy once you normalize per ion, so a generic average value can stand in for the true, structure-specific Madelung constant without needing to know which structure actually applies.

How does lattice energy from these equations connect to the Born-Haber cycle?

The Born-Haber cycle is a Hess's-law thermochemical loop connecting the lattice energy to experimentally measurable quantities: sublimation energy, ionization energies, bond dissociation energy, and electron affinity, all summing to the enthalpy of formation of the solid. Closing that cycle gives an experimentally derived lattice energy that can be directly compared against a Born-Mayer or Kapustinskii prediction. Close agreement (within a few percent) supports a largely ionic bonding picture; a substantial mismatch, with the experimental value more exothermic than predicted, signals extra stabilization from covalent or polarization effects not captured by a purely electrostatic ionic model, as famously seen in silver halides.

Does the Born-Mayer softness parameter ρ really stay constant across compounds?

Approximately, yes, and that transferability is exactly its practical appeal. Values around 30–35 pm reproduce experimental lattice energies reasonably well across many alkali halides and simple oxides, because the repulsion physics — overlap of similarly-sized noble-gas-like closed shells — is broadly similar from one ionic compound to the next. It is not exactly universal: highly polarizable or unusually soft ions (large, heavy halides or chalcogenides) can shift the effective $\rho$ somewhat, which is one reason Born-Mayer, while an improvement, is still an approximation rather than an exact theory.

Can Kapustinskii's equation predict whether a hypothetical compound will actually form?

Not by itself. It only estimates the lattice energy — the energy released upon forming the solid from gaseous ions. Whether a compound is thermodynamically accessible also depends on the ionization energies needed to create the cation, the electron affinities (or lack thereof) needed to create the anion, and atomization/sublimation energies, all combined via a Born-Haber-type cycle. A compound like hypothetical Ca³⁺Cl₃ can have an enormous predicted lattice energy yet still never form, because the extra ionization energy required to strip a third, core electron from calcium vastly outweighs any lattice-energy gain.

The Born-Mayer and Kapustinskii Equations

Where this lives today

The physical picture

The core ideas

Worked examples

Limiting cases, scope & common traps

The chemist behind it

Before this

Where this leads