CFSE and Thermodynamic/Structural Trends

The double-hump analysis is not a historical curiosity; it underlies working design principles across inorganic and materials chemistry. In hydrometallurgy and separation chemistry, the anomalously large CFSE of Ni²⁺ and Co²⁺ (relative to Mn²⁺/Zn²⁺ neighbors) contributes to the different selectivities of solvent-extraction and ion-exchange resins used to separate these metals from ore leachates, since extraction equilibria are governed in part by relative hydration/complexation free energies that inherit the hump structure. In geochemistry and mineralogy, CFSE governs which transition-metal ions preferentially partition into octahedral versus other crystallographic sites in silicate melts and minerals — Cr³⁺ (d³, large octahedral CFSE) and Ni²⁺ (d⁸) show strong octahedral site preference energies that control trace-element partitioning between mantle minerals (olivine, pyroxene) and melt, a tool used in tracing planetary differentiation. In solid-state and materials chemistry, the same $t_{2g}$/$e_g$ occupancy logic (extended to extended-solid band pictures) explains anomalies in the lattice parameters and formation enthalpies of transition-metal oxide and spinel series (e.g., why NiO and Cr₂O₃ deviate from ionic-radius-based predictions of lattice constants more than MnO or ZnO do), directly informing the design of battery cathode materials (layered oxide and spinel Li-ion cathodes) where octahedral-site CFSE contributes to the thermodynamic stability of a given transition-metal oxidation state in the host lattice. In bioinorganic chemistry, the same reasoning explains why certain metalloenzyme active sites show a strong preference for particular metal ions in particular geometries — the CFSE of a $d^6$ low-spin Fe(II) center in a heme pocket, for instance, is part of why iron rather than, say, manganese is favored in oxygen-carrier and oxidase active sites despite otherwise similar ionic radii, tying the abstract double-hump curve to why biology chose the metals it chose. And pedagogically, this analysis is the standard bridge used to introduce the Spectrochemical Series: since the height of each hump scales with $\Delta_o$, and $\Delta_o$ depends sensitively on the ligand (H₂O vs. NH₃ vs. CN⁻), the same double-hump machinery immediately generalizes to compare hydration enthalpies (weak-to-moderate field) against ammine or cyanide complex formation constants (moderate-to-strong field), where much larger $\Delta_o$ values produce correspondingly larger — and, for strong-field ligands, low-spin — CFSE humps.

Ask the tutor: Why does changing from high-spin to low-spin (e.g., with a strong-field ligand like CN⁻) shift the position and height of the CFSE humps, and how does that connect to the ordering of ligands in the spectrochemical series?

Consider the free, gas-phase $M^{2+}$ ion with configuration $d^n$: in the absence of any ligand field, all five $d$ orbitals are degenerate, and the ion's energy, radius, and bonding depend only on nuclear charge and total electron count — a smooth function of atomic number. Now place six ligands around it in an octahedral arrangement. The negative point-charges (or dipole ends) of the ligands raise the energy of all five $d$ orbitals on average (electrostatic repulsion, the "barycenter" shift), but they raise the two orbitals that point directly at the ligands ($d_{z^2}, d_{x^2-y^2}$, called $e_g$) more than the three orbitals that point between the ligands ($d_{xy}, d_{yz}, d_{xz}$, called $t_{2g}$). The splitting energy is $\Delta_o$ (also written $10\,Dq$), with the $t_{2g}$ set lying $0.4\,\Delta_o$ below the barycenter and the $e_g$ set lying $0.6\,\Delta_o$ above it. If you now fill in $n$ electrons according to Hund's rules (subject to whether the field is strong enough to force low-spin pairing), the total electronic energy relative to the hypothetical spherical field barycenter is the CFSE. Because $t_{2g}$ is favorable and $e_g$ is unfavorable, CFSE is zero exactly when the configuration is "electronically neutral" with respect to the splitting — $d^0$ (no electrons to place), $d^{10}$ (both sets completely and symmetrically filled), and high-spin $d^5$ (one electron in each of the five orbitals, again symmetric). Between these three zero points, CFSE rises and falls — twice — because the $t_{2g}$ set fills first (stabilizing) and then the $e_g$ set must start filling once $t_{2g}$ is full or once high-spin pairing forces it (destabilizing relative to further $t_{2g}$ filling). This is the mechanical origin of the double hump: it is baked into the orbital-filling order, exactly analogous to how ionization energies of atoms show a sawtooth pattern from subshell filling order. Because a $t_{2g}$ electron is, on average, "hiding" from the ligands (its lobes point away from them), a metal ion with more $t_{2g}$ occupancy behaves as if it has a smaller effective ionic radius for the purposes of ligand approach — the ligands can get closer before repulsion balances attraction, shortening the $M$–$L$ bond. An $e_g$ electron sits in the ligands' path and pushes them away, effectively inflating the ionic radius. This is why the radius anomaly (structural) and the enthalpy anomaly (thermodynamic) are the same underlying phenomenon viewed through two different measurement lenses: shorter bonds mean stronger (more exothermic) hydration/lattice interactions, by the inverse-radius dependence built into Coulomb's law.

The total lattice (or hydration) enthalpy of a transition-metal ion is modeled as a smooth ionic term plus a CFSE correction:

where: $\Delta H_{\text{ionic baseline}}$ is the value predicted by a purely electrostatic (Born-type) model varying smoothly with effective nuclear charge and ionic radius across the period, and CFSE is the extra stabilization (always $\le 0$, i.e. exothermic contribution) from non-spherical $d$-orbital filling in the ligand field.

Reading: the baseline alone would give a smooth curve (e.g., roughly linear in $Z_{\rm eff}$, since radius contracts roughly linearly with atomic number across a transition series at fixed charge and spin state). The CFSE term is what bends that line into two humps. The crystal-field stabilization energy itself, for an octahedral high-spin or low-spin $d^n$ ion with $t_{2g}^{a}e_g^{b}$ occupancy ($a+b=n$):

where: $a$ = number of $t_{2g}$ electrons, $b$ = number of $e_g$ electrons, and the pairing term $P$ is added once per electron pair beyond the free-ion's baseline pairing (relevant chiefly for comparing high-spin vs. low-spin at the same $d^n$; it does not by itself create the hump pattern, which survives even ignoring $P$).

Reading: because $a$ maxes out at 6 ($t_{2g}^6$) and $b$ maxes out at 4 ($e_g^4$), and because for high-spin configurations $t_{2g}$ fills before $e_g$ starts doubly occupying, $|{\rm CFSE}|$ rises through $d^1$–$d^3$, drops sharply back to zero at $d^5$ (high spin), rises again through $d^6$–$d^8$, and drops back to zero at $d^{10}$ — the double hump, in closed form. The ionic radius correction tracks CFSE with a sign flip, since orbitals populated preferentially contract (t₂g) or expand (eg) the observed radius relative to the "spherical" Goldschmidt/Shannon trend:

where: $r_{\text{smooth}}$ is the radius the ion would have if $d$ electrons contributed no directional shielding effect (i.e., extrapolated smoothly from Ca²⁺ through Zn²⁺ using only nuclear charge), and $\delta r$ is positive (radius larger than smooth trend) when CFSE is small/zero and negative (radius contracted below trend) when $|{\rm CFSE}|$ is large — because strongly-stabilized $t_{2g}$-heavy ions pull ligands in closer.

Reading: since CFSE is most negative (largest magnitude) near $d^3$ and $d^8$ and zero at $d^0, d^5, d^{10}$, the radius correction $\delta r$ is most negative (radius most contracted) at $d^3/d^8$ and zero at $d^0, d^5, d^{10}$ — exactly the double-welled dip seen experimentally.

Example 1 (foundational) — CFSE and the sign/magnitude of the hydration-enthalpy anomaly for Cr²⁺. Cr²⁺ is high-spin $d^4$ in its hexaaqua complex ($t_{2g}^3 e_g^1$). From the CFSE table, CFSE $= -0.6\,\Delta_o$. Using a typical $\Delta_o(\text{Cr(H}_2\text{O)}_6^{2+}) \approx 13{,}900\ \text{cm}^{-1}$ (this ion is unusually high because $d^4$ is Jahn–Teller active and $\Delta_o$ estimates from spectra carry some ambiguity, but we use it for the exercise), convert to energy units:

Then:

This means the measured hydration enthalpy of Cr²⁺ should be roughly 100 kJ mol⁻¹ more exothermic than the smooth baseline value interpolated between Ca²⁺ ($d^0$, no CFSE) and Mn²⁺ ($d^5$, high-spin, no CFSE). Checking limiting cases: if $\Delta_o \to 0$ (a purely ionic, field-free limit — e.g., an extremely weak-field ligand or a gas-phase free ion), CFSE $\to 0$ regardless of $d^n$, and the anomaly vanishes identically, correctly recovering the "no ligand field" baseline for every ion, not just $d^0/d^5/d^{10}$. This is the expected limiting behavior: CFSE is a ligand-field effect and must disappear when the field does.

Example 2 (applied) — Comparing the CFSE-corrected lattice energy of NiO vs. the Born–Landé prediction (Kapustinskii-type sanity check). NiO crystallizes in the rock-salt structure with Ni²⁺ high-spin $d^8$ (in an octahedral oxide environment, $t_{2g}^6 e_g^2$), so CFSE $= -1.2\,\Delta_o$, the maximum-magnitude hump value. Spectroscopic $\Delta_o(\text{NiO}) \approx 8{,}500\ \text{cm}^{-1}$ (roughly consistent with the weaker field of oxide vs. water/ammonia at Ni²⁺, though NiO's band structure is more complex than a simple ionic $\Delta_o$; we proceed with the ionic-model approximation for the exercise). Converting:

The Born–Landé lattice energy of NiO computed from a purely ionic model (charges $+2,-2$; $r_0 \approx 2.084$ Å; Born exponent $n\approx 8$) gives roughly $-3850$ to $-3900$ kJ mol⁻¹ depending on the Madelung constant convention used; adding the CFSE correction brings the ionic-model-plus-CFSE estimate to roughly $-3970$ to $-4020$ kJ mol⁻¹, closer to the experimental Born–Haber-cycle lattice energy of NiO (around $-3900$ to $-4080$ kJ mol⁻¹ depending on the thermochemical cycle used for $\Delta H_f$, sublimation, and ionization data). The key qualitative point survives the numerical uncertainty: the CFSE correction is in the correct direction (more exothermic/more negative lattice energy than the naive ionic model alone) and of the correct rough magnitude (a few percent of the total), exactly as expected for a $d^8$ ion sitting at one of the two humps. Compare the limiting case of ZnO: Zn²⁺ is $d^{10}$, CFSE $=0$ identically, so the ionic Born–Landé model alone (no correction needed) already gives good agreement with experiment — consistent with Zn²⁺ never showing any CFSE-driven anomaly in any lattice or hydration series.

Example 3 (exam/competition level) — Predicting which of two isoelectronic-charge ions has the larger hydration enthalpy magnitude, using CFSE alone. Compare $\text{Fe}^{2+}$ (high-spin $d^6$) and $\text{Mn}^{2+}$ (high-spin $d^5$) hexaaqua ions. Both have charge $+2$ and adjacent atomic number ($Z=26$ vs. $25$), so their smooth ionic baselines differ only slightly (Mn²⁺ has one less proton, very slightly larger radius, very slightly less exothermic baseline hydration enthalpy). Using Shannon radii: Mn²⁺(hs) $=0.83$ Å, Fe²⁺(hs) $=0.78$ Å — note Fe²⁺ is already smaller, consistent with both a higher $Z_{\rm eff}$ and nonzero CFSE. From the CFSE table: Mn²⁺ ($d^5$) has CFSE $=0$; Fe²⁺ ($d^6$) has CFSE $=-0.4\,\Delta_o$. Using $\Delta_o(\text{Fe(H}_2\text{O)}_6^{2+}) \approx 10{,}400\ \text{cm}^{-1} \approx 124.5\ \text{kJ mol}^{-1}$:

So even though Mn²⁺ has one more proton's worth less nuclear charge than Fe²⁺ (which alone would make Mn²�sup>+ hydration slightly less exothermic — a small, expected baseline effect), Fe²⁺'s CFSE adds a further ≈50 kJ mol⁻¹ of exothermicity on top of that, so the total measured hydration enthalpy gap between the two ions ($\Delta H_{\rm hyd}(\text{Mn}^{2+}) \approx -1845$ kJ mol⁻¹ vs. $\Delta H_{\rm hyd}(\text{Fe}^{2+}) \approx -1920$ kJ mol⁻¹, a difference of about 75 kJ mol⁻¹) is larger than the baseline alone would predict from atomic number change alone (which typically accounts for only 15–25 kJ mol⁻¹ per unit $Z$ in this region). The residual (roughly 50–60 kJ mol⁻¹) is attributable to CFSE, matching the calculated $-49.8$ kJ mol⁻¹ within the precision expected from this level of ionic-model approximation. As a limiting-case check: if we instead compared Mn²⁺ ($d^5$, CFSE $=0$) to the hypothetical isoelectronic-charge ion one step further, Zn²⁺ ($d^{10}$, CFSE $=0$ also), both baseline points, the entire enthalpy difference between them (roughly $-2044$ vs. $-1845$ kJ mol⁻¹, a gap of about 200 kJ mol⁻¹ across five atomic numbers) should be attributable almost entirely to the smooth baseline (nuclear charge/radius contraction) with no CFSE contribution at either endpoint — and dividing by 5 gives roughly 40 kJ mol⁻¹ per unit $Z$, a sensible per-proton baseline slope consistent with the smaller Mn–Fe gap analysis above.