Lattice energy ($U$) is the single most important thermodynamic quantity for characterising the stability of an ionic solid. It represents the energy released when one mole of gaseous cations and anions condense from infinite separation into a crystalline lattice, and its magnitude governs melting point, hardness, solubility, and thermal decomposition behaviour.
Because $U$ cannot be measured directly by calorimetry, chemists rely on two theoretical frameworks: the Born-Landé equation (1918), a first-principles electrostatic model requiring full knowledge of the crystal geometry, and the Kapustinskii equation (1956), a semi-empirical shortcut that works when the structure is unknown or when polyatomic ions are involved. This calculator implements both formulations in parallel, allowing rapid, reproducible estimation of $U$ without resorting to lengthy Born-Haber thermochemical cycles.
Required Input Parameters
To perform a Lattice Energy calculation, the following crystallographic and electronic data must be provided:
- Cation charge magnitude ($z^{+}$) — the absolute positive charge (e.g., $+1$ for Na⁺, $+3$ for Al³⁺).
- Anion charge magnitude ($z^{-}$) — the absolute negative charge (e.g., $1$ for Cl⁻, $2$ for O²⁻).
- Cation radius ($r^{+}$) in picometres — the effective ionic radius, typically taken from the Shannon-Prewitt tabulation.
- Anion radius ($r^{-}$) in picometres — the effective ionic radius of the negative species.
- Madelung constant ($M$) — required only for the Born-Landé model; a dimensionless geometric factor reflecting the crystal structure type.
- Born exponent ($n$) — an integer between 5 and 12 reflecting the compressibility of the ions' electron clouds (Born-Landé only).
- Number of ions per formula unit ($\nu$) — required only for the Kapustinskii equation (e.g., $\nu = 2$ for NaCl, $\nu = 5$ for Al₂O₃).
Theoretical Foundation & Formulas
The Coulombic Origin of Lattice Energy
The stability of an ionic crystal ultimately derives from the electrostatic attraction between oppositely charged ions, tempered by short-range Pauli repulsion between their closed electron shells. For a single ion pair, the electrostatic potential energy is given by Coulomb's Law:
$$E_{\text{pair}} = -\frac{z^{+} z^{-} e^{2}}{4\pi\varepsilon_{0} r}$$
In a three-dimensional lattice, however, every ion interacts with every other ion in the crystal. Summing these infinite attractive and repulsive Coulombic terms over a specific geometry produces a dimensionless number known as the Madelung constant ($M$).
The Born-Landé Equation
Max Born and Alfred Landé formalised the full crystal expression in 1918. The attractive Coulombic term is combined with a repulsive term of the form $B/r^{n}$, and minimisation with respect to $r$ at the equilibrium distance $r_{0}$ yields:
$$U = \frac{N_{A} \cdot M \cdot z^{+} \cdot z^{-} \cdot e^{2}}{4\pi\varepsilon_{0} \cdot r_{0}} \left( 1 - \frac{1}{n} \right)$$
For computational convenience, the cluster of fundamental constants is collapsed into a single coefficient. When $r_{0}$ is supplied in picometres and $U$ is returned in kJ/mol, the working form becomes:
$$U = \frac{1389.3 \cdot M \cdot z^{+} \cdot z^{-}}{r_{0}} \left( 1 - \frac{1}{n} \right)$$
The factor $\left(1 - \frac{1}{n}\right)$ is the Born correction. Because $n$ is typically between 5 and 12, this term sits near $0.85-0.92$, meaning short-range repulsion reduces the pure Coulombic attraction by roughly 10-15%.
The Kapustinskii Equation
A. F. Kapustinskii observed that the ratio $M/\nu$ (Madelung constant divided by the number of ions in the formula unit) is nearly constant across common structure types. Exploiting this regularity, he proposed a formula that eliminates the need to know the crystal structure:
$$U = \frac{1202 \cdot \nu \cdot z^{+} \cdot z^{-}}{r^{+} + r^{-}} \left( 1 - \frac{34.5}{r^{+} + r^{-}} \right)$$
Here $r^{+} + r^{-}$ replaces $r_{0}$, both expressed in picometres, and $34.5$ pm is the empirical repulsion parameter $d$ (originally $0.345$ Å). The Kapustinskii constant $1202$ kJ·pm·mol⁻¹ incorporates the averaged Madelung factor and fundamental constants. Published lattice energies obtained from this method typically agree with Born-Landé values to within 5%, which is more than sufficient for most thermochemical prediction work.
The Born Exponent Explained
The exponent $n$ reflects how steeply the Pauli repulsion increases as ion electron clouds begin to overlap. Its value is determined by the noble-gas electron configuration of the ions:
- Light ions (He-like, e.g., Li⁺): $n = 5$
- Ne-like (e.g., Na⁺, F⁻, Mg²⁺, O²⁻): $n = 7$
- Ar-like (e.g., K⁺, Cl⁻, Ca²⁺): $n = 9$
- Kr-like (e.g., Rb⁺, Br⁻): $n = 10$
- Xe-like (e.g., Cs⁺, I⁻): $n = 12$
When the cation and anion belong to different noble-gas configurations, the arithmetic mean is conventionally used (hence $n = 8$ for NaCl, which combines Ne-like Na⁺ and Ar-like Cl⁻).
Technical Specifications & Reference Data
Madelung Constants for Common Structure Types
The following values are built into the calculator and represent the consensus crystallographic literature:
| Structure Type | Example Compound | Coordination (M:X) | Madelung Constant ($M$) |
|---|---|---|---|
| Rock Salt | NaCl | 6:6 | 1.74756 |
| Cesium Chloride | CsCl | 8:8 | 1.76267 |
| Zinc Blende | ZnS (cubic) | 4:4 | 1.63805 |
| Wurtzite | ZnS (hexagonal) | 4:4 | 1.64132 |
| Fluorite | CaF₂ | 8:4 | 2.51939 |
| Rutile | TiO₂ | 6:3 | 2.40800 |
| Corundum | Al₂O₃ | 6:4 | 4.15900 |
Note that for structures containing unequal numbers of ion types (e.g., fluorite, rutile, corundum), the Madelung constant is referenced to a single formula unit, which is why its numerical value exceeds those of simple binary 1:1 salts.
Representative Shannon Effective Ionic Radii (6-Coordinate, in pm)
| Cation | Radius | Anion | Radius |
|---|---|---|---|
| Li⁺ | 76 | F⁻ | 133 |
| Na⁺ | 102 | Cl⁻ | 181 |
| K⁺ | 138 | Br⁻ | 196 |
| Cs⁺ | 167 | I⁻ | 220 |
| Mg²⁺ | 72 | O²⁻ | 140 |
| Ca²⁺ | 100 | S²⁻ | 184 |
| Al³⁺ | 53.5 | N³⁻ | 146 |
Typical Calculated Lattice Energies (kJ/mol)
| Compound | Born-Landé $U$ | Experimental $U$ | Deviation |
|---|---|---|---|
| LiF | 1030 | 1036 | < 1% |
| NaCl | 766 | 788 | ~2.8% |
| KCl | 686 | 717 | ~4.3% |
| MgO | 3795 | 3795 | ~ 0% |
| CaF₂ | 2609 | 2651 | ~1.6% |
Engineering Analysis & Real-World Application
How Each Variable Drives the Result
The governing equations reveal three dominant levers controlling lattice energy magnitude. Understanding their relative weight is essential for interpreting calculator output and for making rational predictions about new compounds.
1. The charge product ($z^{+} \cdot z^{-}$) dominates. Because charges multiply rather than add, doubling both charges (from NaCl-type $+1/-1$ to MgO-type $+2/-2$) quadruples the lattice energy. This is why MgO melts at 2852 °C while NaCl melts at 801 °C despite nearly identical interionic distances. When predicting whether a compound will form an ionic solid at all, the charge product is the first quantity to examine.
2. Interionic distance ($r_{0}$) enters inversely. A reduction in $r_{0}$ from 280 pm (NaCl) to 210 pm (LiF) increases $U$ by roughly $280/210 \approx 33%$. Down a group in the periodic table, ionic radii increase monotonically, which explains the well-documented trend of decreasing lattice enthalpy from LiF → NaF → KF → RbF → CsF.
3. Madelung constant ($M$) is a secondary but non-trivial correction. Differences between rock-salt (1.748) and CsCl (1.763) structures are marginal — under 1%. However, for higher-coordination structures such as fluorite (2.519) or corundum (4.159), the geometric multiplier dramatically inflates $U$.
When to Choose Born-Landé Over Kapustinskii
The two methods are complementary rather than competing. The correct choice depends on what information is available and the required precision:
- Use Born-Landé when the compound adopts a known, well-characterised binary structure (rock salt, fluorite, etc.). This method provides the most rigorous electrostatic treatment and is the default in crystallographic reference work.
- Use Kapustinskii when the structure is unknown, when polyatomic ions such as SO₄²⁻, NO₃⁻, or PO₄³⁻ are present (no clean Madelung constant exists), or when making rapid comparative estimates across a large series of hypothetical compounds.
The Kapustinskii approach is the standard tool for inorganic chemists estimating the stability of novel or unisolated species via Born-Haber thermochemistry, a methodology pioneered by Kapustinskii himself and later refined by H. D. B. Jenkins and others.
Interpreting the Energy Well Curve
The visualisation of lattice potential as a function of interionic distance — often called the Lennard-Jones-style energy well — communicates three facts simultaneously. First, the depth of the well at $r_{0}$ equals $-U$: a deeper well means a more stable crystal. Second, the steepness of the wall on the repulsion (left) side is governed by $n$: larger $n$ produces a stiffer, less compressible solid. Third, the asymmetry of the curve explains thermal expansion: as temperature rises, vibrational amplitude increases preferentially on the shallower attractive side, shifting the mean $r_{0}$ outward.
Limitations of Purely Ionic Models
Both equations implicitly assume 100% ionic bonding, which is rarely true in practice. Significant covalent character — present in compounds such as AgI, ZnS, or any halide of a highly polarising cation — causes the model to underestimate the true bond strength. The lattice energy of AgCl, for instance, exceeds the Born-Landé prediction by approximately 50 kJ/mol due to covalent stabilisation.
Fajans' rules provide the qualitative framework: small, highly charged cations (e.g., Al³⁺, Be²⁺) paired with large, polarisable anions (e.g., I⁻, S²⁻) promote covalent character and, consequently, deviations from ionic-model predictions. When a calculation departs by more than 10% from experimental data, covalency is the most likely culprit.
Frequently Asked Questions
This reflects a longstanding sign convention ambiguity in thermochemistry. The original Born-Landé derivation computes the potential energy of the crystal relative to separated gaseous ions, which is inherently negative (the system is stabilised).
Modern inorganic textbooks, following the IUPAC recommendation, define lattice energy as the energy required to dissociate the crystal into gaseous ions — an endothermic, positive quantity. This calculator reports $U$ as a positive value consistent with this modern convention, where a larger number indicates a more stable crystal. The Coulombic term is shown separately as negative to preserve the physical intuition that attraction stabilises the system.
Three common sources account for most discrepancies. Incorrect ionic radii are the most frequent culprit — the Shannon-Prewitt values depend on coordination number, and using the 6-coordinate radius for a genuinely 8-coordinate structure can introduce errors of 5-10 pm. Always match the radius to the actual coordination environment.
A second cause is unaccounted covalent character, as discussed above. If your compound contains transition metals, $d^{10}$ cations (Cu⁺, Ag⁺, Zn²⁺), or polarisable heavy anions (I⁻, S²⁻, Se²⁻), expect the purely ionic model to underestimate $U$ by anywhere from 2% to 20%.
Finally, check the Born exponent. A mixed-configuration compound such as LiI (He-like Li⁺, $n=5$; Xe-like I⁻, $n=12$) requires the mean value $n=8.5$, not the value of either ion alone. Small changes in $n$ shift the Born correction by several percent.
Yes, but only via the Kapustinskii pathway, and with an important caveat. Polyatomic ions such as SO₄²⁻ or NH₄⁺ have no well-defined Madelung constant, so the Born-Landé route is inapplicable.
For Kapustinskii calculations, substitute the thermochemical radius — an effective spherical radius derived by treating the polyatomic ion as a sphere of equivalent electrostatic behaviour. Common values include SO₄²⁻ ≈ 230 pm, NO₃⁻ ≈ 189 pm, CO₃²⁻ ≈ 185 pm, and PO₄³⁻ ≈ 238 pm. These values were originally tabulated by Kapustinskii and later refined by Jenkins and Thakur. Using them, one can successfully predict the lattice energies of a vast range of salts for which no crystallographic radius data exists.
Professional Conclusion
Lattice energy is the thermodynamic keystone of solid-state inorganic chemistry, and its accurate estimation underpins any rational prediction of compound stability, solubility, and reactivity. The Born-Landé and Kapustinskii equations represent two complementary approaches: the former rigorous but structure-dependent, the latter empirical but universal.
Manual execution of these equations is straightforward in principle but tedious in practice, with the conversion of SI fundamental constants into working units (kJ/mol from pm) a notorious source of arithmetic error. Automated calculation eliminates this friction, delivers reproducible energy decompositions into attractive and repulsive terms, and enables rapid comparative studies across structure types and charge states.
For professional inorganic chemists, materials scientists, and advanced students preparing Born-Haber analyses, this tool replaces an hour of careful arithmetic with a reproducible numerical estimate — one that, when interpreted with awareness of covalent contributions and radius conventions, captures the physics of ionic bonding to the precision required by modern research.