Square-Planar Crystal Field Splitting

Dissolve a nickel(II) salt in water and add cyanide. The pale green octahedral aqua complex collapses into a new form. It becomes an intensely yellow, diamagnetic, perfectly flat ion, [Ni(CN)₄]²⁻. Nothing about nickel's electron count changed — it is still d⁸. Yet the magnetism flipped from paramagnetic to diamagnetic. The geometry flattened from a sphere-like octahedron to a square. The culprit is a change in coordination geometry. The consequence is a completely different splitting of the five d orbitals. This is the subject of square-planar crystal field theory. It examines what happens to the five-fold-degenerate d orbitals when six ligands become four. These four ligands arrange in a plane. It also explains why this particular geometry is so strongly and specifically favored by d⁸ metal ions. Square-planar coordination is not a geometric curiosity. It is the geometry of some of the most important molecules in chemistry. Vaska's complex, Wilkinson's catalyst's relatives, and Zeise's salt all depend on it. The anticancer drug cisplatin does too. So do the heme-like macrocycles of vitamin B₁₂ and chlorophyll analogs. All rely on square-planar or near-square-planar metal centers. Square-planar platinum(II) and palladium(II) complexes are the workhorses of cross-coupling catalysis. Square-planar nickel(II), palladium(II), and gold(III) complexes are textbook diamagnetic d⁸ systems. They are used to teach ligand field theory. This is precisely because their splitting pattern is derivable from first principles. One takes an octahedron and pulls two trans ligands off to infinity. The big picture is this: square-planar splitting is not an independent pattern to be memorized from scratch. It is what an Octahedral Crystal Field Splitting diagram becomes under a specific distortion. This is an extreme tetragonal (Jahn–Teller-like) distortion. It involves elongation along one axis carried to its logical limit. At that limit, the two axial ligands are removed entirely. This connects smoothly to octahedral theory. Square-planar geometry is also the endpoint of the exact distortion mode. That mode drives Jahn-Teller Distortion in octahedral d⁹ and high-spin d⁴ complexes. So square-planar splitting sits at a conceptual crossroads. It explains why certain octahedral complexes distort. It also explains why certain d⁸ ions abandon octahedral coordination altogether. This article builds the square-planar dd-orbital splitting diagram from the octahedral one. It does so by correlating orbital energies as axial ligands are withdrawn. It derives the ordering and approximate spacing of all four energy levels (sometimes drawn as five). It explains why strong-field d⁸ configurations achieve exceptional stabilization in this geometry. It also works through Ni(II), Pd(II)/Pt(II), and a diagnostic magnetism problem in full quantitative detail.

Key takeaways
  • Square-planar splitting is derived by removing the two axial ligands from an octahedron and following each dd orbital's energy as zz-ligands recede.
  • The dx2−y2d_{x^2-y^2} orbital points directly at all four ligands and is by far the highest in energy; dz2d_{z^2} drops sharply as axial ligands leave.
  • The ordering (low → high) is typically dxz,dyz<dz2<dxy<dx2−y2d_{xz},d_{yz} < d_{z^2} < d_{xy} < d_{x^2-y^2}, giving four resolvable levels from five orbitals.
  • d⁸ ions (Ni²⁺, Pd²⁺, Pt²⁺, Au³⁺, Rh⁺, Ir⁺) with strong-field ligands place 8 electrons in the lower four orbitals, leaving dx2−y2d_{x^2-y^2} empty — a large HOMO–LUMO gap favors diamagnetic square-planar geometry over tetrahedral or octahedral alternatives.
  • Crystal field stabilization energy (CFSE) analysis shows square-planar d⁸ can exceed octahedral CFSE for the same ligand set once electron-pairing costs are included, explaining the geometric preference.
  • Key idea: Square-planar splitting is the octahedral pattern carried to the limit of removing two trans ligands, and it is uniquely stabilizing for low-spin d⁸ ions because it opens a large gap above exactly eight electrons' worth of orbitals.

Where this lives today

Square-planar Pt(II) and Pd(II) complexes are the backbone of modern homogeneous catalysis: Pd(0)/Pd(II) cycles in Suzuki, Heck, and Negishi cross-coupling reactions pass through square-planar Pd(II) intermediates whose ligand substitution kinetics (associative, via a five-coordinate trigonal-bipyramidal transition state) are a direct consequence of the accessible, empty, high-energy dx2−y2d_{x^2-y^2} orbital pointing exactly at the site of incoming nucleophile attack. Wilkinson-type and related Rh(I)/Ir(I) square-planar d⁸ catalysts (themselves isoelectronic with Pt(II)/Pd(II)) drive industrial hydrogenation and hydroformylation, their oxidative-addition chemistry likewise gated by the same empty high-lying orbital. In medicine, cisplatin — cis-diamminedichloroplatinum(II) — is a square-planar Pt(II) d⁸ complex whose geometry is mechanistically essential: only the cis isomer can chelate adjacent guanine bases on DNA to form the crosslink responsible for its anticancer activity, while the trans isomer (same square-planar splitting, different substitution pattern) is pharmacologically inactive. The rigid, substitution-inert-but-tunable square-planar Pt(II) center, governed by exactly the splitting diagram derived above, is the structural reason cisplatin and its descendants (carboplatin, oxaliplatin) function as targeted DNA-damaging agents rather than indiscriminately reactive species. In materials science, square-planar d⁸ complexes of Pt(II), Pd(II), and Au(III) with extended aromatic or cyanide ligands are prized for their photoluminescence and metal–metal (often called metallophilic, dz2d_{z^2}–dz2d_{z^2}) interactions in the solid state, exploited in OLED phosphors and molecular electronics; the low-lying, partially populated dz2d_{z^2} orbital derived directly in this splitting diagram is precisely the orbital responsible for these intermolecular Pt···Pt and Au···Au stacking interactions. Bioinorganic chemistry relies on the same picture for vitamin B₁₂'s cobalt corrin (formally octahedral but with a strongly tetragonally distorted, near-square-planar equatorial field) and for synthetic metalloporphyrin mimics of heme and chlorophyll, where the square-planar-like equatorial ligand field sets the metal's spin state and reactivity toward axial ligand binding.

:::quiz[Why is octahedral d⁸ (e.g. [Ni(H2O)6]2+) always paramagnetic regardless of ligand field strength?]

  • Because dz² and dx²−y² are never degenerate in Oh
  • Because the eg set has only two orbitals, forcing one electron into each by Hund's rule for the last two electrons ✓
  • Because water is always a weak-field ligand
  • Because Ni2+ cannot achieve a d8 configuration in octahedral fields

Ask the tutor: How would the square-planar splitting diagram and its energies change if the four ligands were strong π-acceptors like CO instead of pure σ-donors?

The physical picture

Begin with an octahedral complex, ligands along ±x,±y,±z\pm x, \pm y, \pm z, giving the familiar splitting into t2gt_{2g} (three orbitals: dxy,dxz,dyzd_{xy}, d_{xz}, d_{yz}, lower) and ege_g (two orbitals: dz2,dx2−y2d_{z^2}, d_{x^2-y^2}, higher), separated by Δo\Delta_o. Now slide the two ligands on the zz-axis away from the metal, symmetrically, while keeping the four equatorial ligands fixed. This is a tetragonal elongation. Any orbital with electron density along zz feels reduced repulsion from the departing ligands and drops in energy; any orbital with density concentrated in the xyxy-plane feels increased relative repulsion from the four ligands that are now proportionally "closer" (or at least unchanged while their axial competitors have weakened) and rises. Four orbitals split apart from the original two levels. From ege_g: dz2d_{z^2} (strong axial lobes, big ±z\pm z density) drops steeply; dx2−y2d_{x^2-y^2} (density concentrated in-plane, pointing directly at the four remaining ligands) rises further still, ending up the highest orbital of all. From t2gt_{2g}: dxyd_{xy}, whose lobes lie in the xyxy-plane between the ligand directions, rises somewhat because equatorial ligand repulsion is now felt more acutely without competing axial destabilization elsewhere in the diagram; dxzd_{xz} and dyzd_{yz}, which have equal density split between the zz-axis and the xyxy-plane, drop because the axial component of their density loses its repulsive partner and remains a degenerate pair (by the surviving C4C_4 symmetry, xzxz and yzyz stay equivalent under reflection through the plane). Carried to the limit of complete axial ligand removal (idealized square-planar, D4hD_{4h} symmetry), the four-level pattern is

dxz,dyz (eg)  <  dz2 (a1g)  <  dxy (b2g)  <  dx2−y2 (b1g)d_{xz},d_{yz}\ (e_g) \;<\; d_{z^2}\ (a_{1g}) \;<\; d_{xy}\ (b_{2g}) \;<\; d_{x^2-y^2}\ (b_{1g})

read low to high energy. Crucially, note that dz2d_{z^2} drops so far that it can fall below dxyd_{xy} — a nontrivial, ligand-dependent detail that some textbooks omit; the ordering of dz2d_{z^2} relative to dxyd_{xy} can invert depending on the specific ligand's σ\sigma- and π\pi-donor/acceptor character, but the identity of the extremes — dxz/dyzd_{xz}/d_{yz} lowest, dx2−y2d_{x^2-y^2} highest by a wide margin — is essentially universal for real square-planar d⁸ complexes.

Figure
Correlation diagram: octahedral → square-planar splitting as axial ligands recede
Octahedral (Oh) Square planar (D4h) axial ligands receding → eg: dz², dx²−y² t2g: dxy, dxz, dyz dx²−y² (b1g) — highest dz² (a1g) dxy (b2g) dxz, dyz (eg) — lowest, degenerate energy ↑

The second key physical fact is electron count. In a strong-field square-planar d⁸ ion, the eight d electrons fill the lower four orbitals — dxz,dyz,dz2,dxyd_{xz}, d_{yz}, d_{z^2}, d_{xy} — completely, each doubly occupied, and leave dx2−y2d_{x^2-y^2} empty. Because dx2−y2d_{x^2-y^2} sits far above the rest, this is a large HOMO–LUMO gap, all electrons are paired, and the complex is diamagnetic. Compare an octahedral d⁸ ion: it must place two electrons in ege_g (three in each t2gt_{2g} orbital, then one in each ege_g orbital by Hund's rule for a modest Δo\Delta_o, or paired in one t2gt_{2g}... actually for d⁸ in OhO_h, six electrons fill t2gt_{2g} completely and two electrons fill ege_g singly, giving t2g6eg2t_{2g}^6 e_g^2, invariably paramagnetic with two unpaired electrons regardless of field strength, since ege_g has only two orbitals to distribute two electrons into with parallel spin preferred). Square-planar geometry is the only common geometry in which a d⁸ ion can achieve full electron pairing in exactly four low-lying orbitals with a genuinely large gap to the fifth — this is the physical root of d⁸ ions' famous preference for square-planar coordination.

Figure
Electron filling: octahedral d⁸ (paramagnetic) vs square-planar d⁸ (diamagnetic)
Octahedral d⁸: t2g⁶eg² Square planar d⁸ ↑↑ ↑ eg (2 unpaired) ↑↓ ↑↓ ↑↓ t2g (filled) dx²−y² empty (LUMO) ↑↓ dxy filled ↑↓ dz² filled ↑↓ ↑↓ dxz, dyz filled All 8 electrons paired → diamagnetic

The core ideas

The correlation between octahedral and square-planar splitting parameters is normally expressed relative to Δo\Delta_o, the octahedral splitting parameter, using empirical/semi-quantitative relations first systematized for tetragonally distorted complexes:

E(dx2−y2)≈+1.28 Δo,E(dz2)≈−0.43 ΔoE(d_{x^2-y^2}) \approx +1.28\,\Delta_o, \qquad E(d_{z^2}) \approx -0.43\,\Delta_o
E(dxy)≈+0.23 Δo,E(dxz)=E(dyz)≈−0.54 ΔoE(d_{xy}) \approx +0.23\,\Delta_o, \qquad E(d_{xz})=E(d_{yz}) \approx -0.54\,\Delta_o

where:

  • Δo\Delta_o is the parent octahedral splitting for the same metal–ligand pair (a hypothetical or actual six-coordinate reference),
  • energies are measured relative to the spherical (barycenter) average, so the weighted mean over all five orbitals is zero,
  • these particular coefficients come from angular-overlap/point-charge treatments of the limiting square-planar geometry (axial ligands fully removed) and vary modestly with the specific ligand and metal.

Reading this: the barycenter rule from Octahedral Crystal Field Splitting still applies — the sum of all five orbital energies, weighted by degeneracy, is zero, because crystal field splitting redistributes energy without changing the total (it is not a real bond-energy input, just an electrostatic reshuffling of the dd-orbital manifold). Check the weighted sum with the coefficients above:

1(+1.28)+1(−0.43)+1(+0.23)+2(−0.54)=1.28−0.43+0.23−1.08=0.001(+1.28) + 1(-0.43) + 1(+0.23) + 2(-0.54) = 1.28 - 0.43 + 0.23 - 1.08 = 0.00

which confirms internal consistency (the barycenter condition holds exactly with these coefficients). The square-planar splitting parameter analogous to Δo\Delta_o or Δt\Delta_t is usually defined as the HOMO–LUMO-type gap between dxyd_{xy} (or, in some conventions, the whole lower set) and dx2−y2d_{x^2-y^2}:

Δsp=E(dx2−y2)−E(dxy)≈(1.28−0.23) Δo=1.05 Δo\Delta_{sp} = E(d_{x^2-y^2}) - E(d_{xy}) \approx (1.28 - 0.23)\,\Delta_o = 1.05\,\Delta_o

where:

  • Δsp\Delta_{sp} is loosely called Δ1\Delta_1 in older spectroscopic literature on square-planar d⁸ complexes,
  • the relation shows that even though square-planar has fewer ligands than octahedral, the top orbital dx2−y2d_{x^2-y^2} is pushed up almost as far above the mean as Δo\Delta_o itself, because all four ligand point-charges/dipoles are now concentrated in-plane with no axial ligands to "share" the repulsive burden.

The overall crystal field stabilization energy (CFSE) for a configuration in square-planar geometry is computed exactly as in octahedral or tetrahedral cases: sum each electron's orbital energy (in units of Δo\Delta_o, using the coefficients above) using the observed spin state, and compare to zero (no CFSE, spherical field).

Derivation ACFSE of low-spin d⁸ in square-planar geometry, and comparison to octahedral d⁸

Goal. Compute the crystal field stabilization energy of a d⁸ ion in square-planar geometry using the correlation coefficients above, and compare it quantitatively to the same d⁸ ion in a (hypothetical, same-ligand) octahedral field, to show why square-planar wins for strong-field d⁸ systems once pairing energy is included.

Step 1 — assign electrons to square-planar orbitals. For a strong-field d⁸ ion, all four lower orbitals are filled: dxz2dyz2dz22dxy2d_{xz}^2 d_{yz}^2 d_{z^2}^2 d_{xy}^2, and dx2−y2d_{x^2-y^2} is empty. Using the energies above (each in units of Δo\Delta_o): dxz=dyz=−0.54d_{xz}=d_{yz}=-0.54, dz2=−0.43d_{z^2}=-0.43, dxy=+0.23d_{xy}=+0.23.

Step 2 — sum orbital energies weighted by occupation. Each of the four filled orbitals holds 2 electrons:

CFSEsp=2(−0.54)+2(−0.54)+2(−0.43)+2(+0.23)\text{CFSE}_{sp} = 2(-0.54) + 2(-0.54) + 2(-0.43) + 2(+0.23)
CFSEsp=−1.08−1.08−0.86+0.46=−2.56 Δo\text{CFSE}_{sp} = -1.08 - 1.08 - 0.86 + 0.46 = -2.56\,\Delta_o

This is the electrostatic stabilization before accounting for the extra electron-pairing energy that the strong field forces relative to a free-ion, all-parallel-spin reference; that pairing cost, PP, is the same physical quantity used in low-spin/high-spin octahedral analysis. In d⁸ square-planar, four electron pairs beyond the single pair already required by any d⁸ configuration must be counted; conventionally one compares apples to apples against octahedral d⁸.

Step 3 — compute octahedral d⁸ CFSE for the same ligand (same Δo\Delta_o). Octahedral d⁸ is always t2g6eg2t_{2g}^6 e_g^2 (this configuration is forced regardless of field strength, because there are only two ege_g orbitals for the last two electrons and Hund's rule keeps them unpaired unless Δo\Delta_o vastly exceeds pairing energy, which does not happen for real octahedral d⁸ systems — hence octahedral Ni(II) is always paramagnetic, high-spin in the sense that mattes here). Using standard octahedral energies eg=+0.6Δoe_g = +0.6\Delta_o, t2g=−0.4Δot_{2g} = -0.4\Delta_o:

CFSEoh=6(−0.4)+2(+0.6)=−2.4+1.2=−1.2 Δo\text{CFSE}_{oh} = 6(-0.4) + 2(+0.6) = -2.4 + 1.2 = -1.2\,\Delta_o

Step 4 — compare electronic (non-pairing) stabilization. The square-planar arrangement already gives more electrostatic stabilization from the dd-orbital pattern alone:

Δ(CFSE)=CFSEsp−CFSEoh=−2.56−(−1.20)=−1.36 Δo\Delta(\text{CFSE}) = \text{CFSE}_{sp} - \text{CFSE}_{oh} = -2.56 - (-1.20) = -1.36\,\Delta_o

Step 5 — include the extra pairing energy cost. Square-planar d⁸ has four filled (paired) orbitals: that is 4 electron pairs total in dd-orbitals (dxz2,dyz2,dz22,dxy2d_{xz}^2,d_{yz}^2,d_{z^2}^2,d_{xy}^2). Octahedral high-spin d⁸ has t2g6eg2t_{2g}^6 e_g^2: three pairs in t2gt_{2g} and the two ege_g electrons unpaired (one pair total among t2gt_{2g} orbitals times three, i.e. 3 pairs, plus zero pairs in ege_g) — so octahedral d⁸ has 3 electron pairs, square-planar d⁸ has 4 electron pairs. The extra pairing energy cost of forming square-planar from octahedral is +1P+1P (one additional forced pairing), where PP is the mean pairing energy (typically 15,000–20,000 cm⁻¹ for 3d metals, smaller relative to Δ\Delta for 4d/5d metals). Net stabilization of square-planar over octahedral becomes

ΔEsp−oh=−1.36 Δo+P\Delta E_{sp-oh} = -1.36\,\Delta_o + P

For a strong-field ligand where Δo\Delta_o is large (cyanide, strong-field phosphines, or intrinsically for 4d/5d metals where Δo\Delta_o is much larger than 3d), −1.36Δo-1.36\Delta_o dominates over +P+P, and square-planar geometry is strongly favored electronically — consistent with Ni(CN)₄²⁻, PdCl₄²⁻, PtCl₄²⁻, and virtually all Pd(II)/Pt(II) complexes being square-planar. For weak-field ligands on Ni²⁺ (small Δo\Delta_o), PP can dominate and octahedral (or tetrahedral) geometry is competitive — consistent with NiCl₄²⁻ being tetrahedral and [Ni(H₂O)₆]²⁺ being octahedral, both paramagnetic.

Limiting-case check. As Δo→0\Delta_o \to 0 (a purely ionic, very weak field, e.g., very large soft counter-ions with negligible field), ΔEsp−oh→P>0\Delta E_{sp-oh}\to P>0, meaning octahedral (or a geometry avoiding forced pairing) is favored — correctly predicting that weak-field d⁸ ions like high-spin Ni²⁺ with weak ligands remain octahedral or tetrahedral, never square-planar, matching experiment (Ni²⁺ square-planar complexes are essentially always found with strong-field ligands: CN⁻, low-coordination-number phosphines, dimethylglyoximate chelates).

Go deeper · Derivation B · 2 minWhy dx2−y2d_{x^2-y^2} is the unique highest orbital: an angular-overlap argument
Derivation BWhy dx2−y2d_{x^2-y^2} is the unique highest orbital: an angular-overlap argument

Goal. Show quantitatively, using a simple point-charge/angular-overlap argument, why dx2−y2d_{x^2-y^2} rises so far above the other four orbitals when axial ligands are removed, rather than relying purely on the correlation-diagram picture.

Step 1 — set up the angular-overlap model (AOM) framework. In AOM, each metal–ligand σ\sigma-bond contributes an energy raise eσe_\sigma to whichever dd orbital points most directly at that ligand, weighted by the square of the angular overlap factor for that orbital in that ligand's direction. For an octahedral field with six ligands (each contributing eσe_\sigma), the standard AOM results reproduce eg=3eσe_g = 3e_\sigma and t2g=0t_{2g}=0 (pure σ\sigma-only octahedral field), consistent with Δo=3eσ\Delta_o = 3e_\sigma.

Step 2 — remove the two axial ligands. In square-planar (four equatorial ligands along ±x,±y\pm x,\pm y, no axial ligands), each equatorial σ\sigma-donor contributes eσe_\sigma to dx2−y2d_{x^2-y^2} with full angular weight (since dx2−y2d_{x^2-y^2} points exactly along xx and yy), giving

E(dx2−y2)=4×34eσ×(normalization)E(d_{x^2-y^2}) = 4\times \tfrac{3}{4} e_\sigma \times (\text{normalization})

Working through the standard AOM angular factors for dx2−y2d_{x^2-y^2} with respect to ligands along x,yx,y (angular overlap factor =1=1 for each of 4 ligands in-plane, using the convention that a single on-axis σ\sigma donor contributes eσe_\sigma fully) gives the standard textbook result

E(dx2−y2)=3eσE(d_{x^2-y^2}) = 3e_\sigma

for four in-plane σ\sigma-only donors, essentially unchanged in magnitude from its octahedral ege_g contribution (ege_g in OhO_h also receives 3eσ3e_\sigma total, shared with dz2d_{z^2}) — but now dx2−y2d_{x^2-y^2} does not have to share destabilization with an axial partner, because dz2d_{z^2}'s axial lobes have lost their ligands.

Step 3 — track dz2d_{z^2}. In octahedral geometry, dz2d_{z^2} receives contributions from both axial ligands (large lobes) and equatorial ligands (small "doughnut" lobes), summing (by the standard AOM decomposition) to the same ege_g energy as dx2−y2d_{x^2-y^2}, 3eσ3e_\sigma, when averaged over the octahedral field by symmetry. Removing the axial ligands removes the dominant contribution to dz2d_{z^2}'s destabilization (its axial lobes were the primary overlap direction), leaving only the smaller equatorial "doughnut" overlap. The AOM equatorial contribution to dz2d_{z^2} from four in-plane donors is eσe_\sigma (roughly one-third the in-plane contribution to dx2−y2d_{x^2-y^2}, because the equatorial angular factor for dz2d_{z^2}, (cos⁡2θ−12sin⁡2θ)2(\cos^2\theta - \tfrac12\sin^2\theta)^2 evaluated at θ=90°\theta=90°, is 1/41/4 compared to dx2−y2d_{x^2-y^2}'s equatorial factor of 3/43/4):

E(dz2)≈eσ(much less than E(dx2−y2)=3eσ)E(d_{z^2}) \approx e_\sigma \quad (\text{much less than } E(d_{x^2-y^2})=3e_\sigma)

Step 4 — the qualitative payoff. Even though both dz2d_{z^2} and dx2−y2d_{x^2-y^2} derive from the octahedral ege_g set and would be degenerate in a spherical or octahedral field, losing the axial ligands strips dz2d_{z^2} of most of its destabilizing overlap while leaving dx2−y2d_{x^2-y^2}'s in-plane overlap untouched — hence the large observed gap between them in square-planar complexes, consistent with the correlation-diagram coefficients (+1.28Δo+1.28\Delta_o vs −0.43Δo-0.43\Delta_o) derived independently above.

Limiting-case check. If instead of removing the axial ligands we only lengthened them slightly (partial tetragonal elongation, not full removal), AOM predicts a smooth, continuous drop in E(dz2)E(d_{z^2}) and rise in E(dx2−y2)E(d_{x^2-y^2}) as the axial σ⊥\sigma_\perp term shrinks — recovering the octahedral limit exactly when axial and equatorial bond lengths become equal, i.e. eσaxial=eσequatoriale_\sigma^{axial}=e_\sigma^{equatorial} collapses back to t2g/egt_{2g}/e_g degeneracy. This continuity is exactly the physical basis for treating square-planar splitting as the endpoint of a tetragonal Jahn-Teller Distortion pathway rather than a wholly separate phenomenon.

Worked examples

Example 1 (foundational): construct the CFSE of low-spin square-planar Ni(CN)₄²⁻ and confirm diamagnetism. Ni²⁺ is d⁸ ([Ar]3d8[Ar]3d^8). Cyanide is an exceptionally strong-field ligand (very high in the spectrochemical series), driving Ni²⁺ into a square-planar, low-spin, diamagnetic configuration rather than the tetrahedral or octahedral alternatives typical of weaker-field Ni²⁺ complexes. Fill the four square-planar orbitals from lowest energy up, obeying the Pauli principle, with 8 electrons:

(dxz)2(dyz)2(dz2)2(dxy)2(dx2−y2)0(d_{xz})^2(d_{yz})^2(d_{z^2})^2(d_{xy})^2(d_{x^2-y^2})^0

All eight electrons occupy four completely filled, spin-paired orbitals. There are zero unpaired electrons — the complex must be diamagnetic, exactly as observed experimentally for [Ni(CN)₄]²⁻ (measured magnetic moment μ≈0\mu \approx 0, versus μ≈2.83 μB\mu \approx 2.83\ \mu_B expected for two unpaired electrons in a hypothetical paramagnetic alternative). Using the correlation coefficients, the CFSE is (from the derivation above) −2.56 Δo-2.56\,\Delta_o, a large stabilization consistent with cyanide's very large Δo\Delta_o for Ni²⁺ (on the order of 30,000–34,000 cm⁻¹ range for related strong-field Ni(II)–CN systems, though the exact experimental Δo\Delta_o for a hypothetical six-coordinate Ni(CN)₆⁴⁻ reference is not directly measurable since that species does not form; the comparison is conceptual/thermodynamic rather than a literal spectroscopic subtraction).

Limiting case check: if cyanide were instead a very weak-field ligand (hypothetically, keeping only the geometry fixed), Δo→\Delta_o \to small, and the pairing-energy penalty computed in the derivation would dominate, making square-planar geometry electronically unfavorable relative to tetrahedral — correctly predicting that weak-field Ni²⁺ complexes (e.g., [NiCl₄]²⁻, [Ni(H�2O)₆]²⁺) are not square-planar, but tetrahedral or octahedral respectively, both paramagnetic with two unpaired electrons.

Example 2 (applied): why PtCl₄²⁻ is square-planar and diamagnetic while the isoelectronic-in-charge NiCl₄²⁻ is tetrahedral and paramagnetic. Both ions are d⁸ metal dichlorides with four chloride ligands, yet they adopt opposite geometries: [PtCl₄]²⁻ is square-planar and diamagnetic; [NiCl₄]²⁻ is tetrahedral and paramagnetic (two unpaired electrons, μ≈2.9\mu\approx 2.9–3.2 μB\mu_B typically observed for tetrahedral Ni(II)). The geometric divergence, despite an identical ligand and identical formal d-electron count, is a direct test of the CFSE-versus-pairing-energy competition derived above. Two factors both push in the same direction. First, Δo\Delta_o (and by extension Δsp\Delta_{sp}, roughly 1.05Δo1.05\Delta_o by the relation derived above) scales up sharply down a group: for the same ligand set, Δo(Pt2+)≫Δo(Ni2+)\Delta_o(\text{Pt}^{2+}) \gg \Delta_o(\text{Ni}^{2+}), typically by a factor of roughly 1.5–2× per row descended (3d → 4d → 5d), because the more diffuse 5d orbitals overlap ligand orbitals more strongly and 5d orbitals are less contracted by electron–electron repulsion. Second, the mean pairing energy PP decreases down a group (heavier d-block metals have smaller interelectronic repulsion per orbital due to more diffuse orbitals), so the pairing-energy penalty in the derivation's Step 5, ΔEsp−oh=−1.36Δo+P\Delta E_{sp-oh} = -1.36\Delta_o + P, becomes overwhelmingly negative (square-planar favored) for Pt²⁺: Δo\Delta_o is large and PP is comparatively small. For Ni²⁺ (3d), Δo\Delta_o from chloride (a weak/moderate-field ligand low in the spectrochemical series) is small, while PP for the compact 3d orbitals is comparatively large, so −1.36Δo+P>0-1.36\Delta_o + P > 0: pairing costs dominate, and tetrahedral geometry — which has an even smaller CFSE magnitude but avoids the extra forced pairing altogether, and also minimizes ligand–ligand steric repulsion for only four large chloride ions — wins out. This is a clean illustration of how the same qualitative formula, applied with correctly scaled parameters across a group, flips the predicted ground-state geometry.

Example 3 (exam/competition level): estimate Δsp\Delta_{sp} for [PtCl₄]²⁻ from its d–d absorption spectrum and check internal consistency with the correlation formula. Square-planar Pt(II) complexes typically show three spin-allowed d–d transitions (from filled dxz,yzd_{xz,yz}, dz2d_{z^2}, dxyd_{xy} up to the empty dx2−y2d_{x^2-y^2} level, in D4hD_{4h} symmetry with appropriate orbital-symmetry labels 1A1g→1A2g,1B1g,1Eg^1A_{1g}\to {}^1A_{2g}, {}^1B_{1g}, {}^1E_g manifold, though d–d bands in Pt(II) are often obscured by intense charge-transfer absorption). Suppose (for a representative strong-field square-planar Pt(II) chloride-type complex) the highest-energy, most clearly assigned spin-allowed d–d band — corresponding to the dxy→dx2−y2d_{xy}\to d_{x^2-y^2} transition, i.e. Δsp\Delta_{sp} itself — is observed at ν~=20,300 cm−1\tilde\nu = 20{,}300\ \text{cm}^{-1}. Using the correlation-diagram relation derived above, Δsp≈1.05 Δo\Delta_{sp}\approx 1.05\,\Delta_o, back out the equivalent octahedral splitting parameter that would characterize a hypothetical six-coordinate PtCl₆⁴⁻-type reference with the same metal–ligand bond character:

Δo≈Δsp1.05=20,300 cm−11.05=19,333 cm−1\Delta_o \approx \frac{\Delta_{sp}}{1.05} = \frac{20{,}300\ \text{cm}^{-1}}{1.05} = 19{,}333\ \text{cm}^{-1}

Compare this to the experimentally known Δo\Delta_o for the genuinely octahedral hexachloroplatinate(IV) ion, [PtCl₆]²⁻ (Pt in the +4 oxidation state, admittedly a different oxidation state and thus not a strictly rigorous comparison, but illustrative of the right order of magnitude): Δo\Delta_o for Pt(IV)–Cl is reported in the range of roughly 20,000–23,000 cm⁻¹ for octahedral hexahalide platinum(IV) complexes. The back-calculated value of 19,300 cm⁻¹ sits close to, if somewhat below, this range — a reasonable outcome given that Pt(II) has one less positive charge than Pt(IV) and would be expected to produce a somewhat smaller Δo\Delta_o for the same ligand, everything else equal (lower oxidation state → weaker ligand field, consistent with periodic trends established for Octahedral Crystal Field Splitting).

Consistency/limiting check: if Δsp\Delta_{sp} were instead measured as exactly 1.05Δo1.05\Delta_o with Δo\Delta_o taken as a round number, say Δo=20,000 cm−1\Delta_o = 20{,}000\ \text{cm}^{-1}, the predicted Δsp\Delta_{sp} would be 21,000 cm−121{,}000\ \text{cm}^{-1} — demonstrating that the correlation coefficient of 1.051.05 is a genuinely small perturbation on Δo\Delta_o (square-planar splitting is comparable in magnitude to octahedral splitting for the same ligands, not dramatically larger or smaller), which is why d–d spectra of square-planar and octahedral complexes of the same metal/ligand family often fall in overlapping spectral windows.

Limiting cases, scope & common traps

Regime / limitBehaviorPhysical reason
Axial ligands only partially removed (moderate tetragonal elongation)Smooth intermediate splitting between OhO_h and D4hD_{4h} patternsContinuous correlation diagram; recovers this article's endpoint only when axial ligands are fully removed
Axial ligands fully removed, ligand is weak-fieldSquare-planar geometry usually not favored; tetrahedral or octahedral competesPairing-energy penalty PP exceeds electronic CFSE gain, per the derivation
Axial ligands fully removed, ligand is strong-field, metal is 4d/5d d⁸Square-planar strongly favored, diamagneticLarge Δo\Delta_o, small PP (diffuse 4d/5d orbitals) — dominant regime for Pd(II), Pt(II), Au(III)
d⁴ (high-spin, e.g. Cr²⁺, Mn³⁺) or d⁹ (Cu²⁺) in octahedral fieldDistorts toward tetragonal elongation but usually stops short of true square-planarJahn-Teller Distortion driven by degenerate ege_g occupancy, but electronic driving force is smaller than for d⁸; genuinely four-coordinate square-planar d⁹ and d⁴ complexes exist (many Cu²⁺ complexes) but retain weak axial interactions (elongated octahedral/"4+2" coordination) rather than rigorously D4hD_{4h} square-planar
d⁸ configuration generallyNearly always square-planar with strong-field ligands, tetrahedral or octahedral with weak-field ligandsOnly geometry giving full pairing in four low orbitals with a genuinely large HOMO–LUMO gap
n=3n=3 (first row) d⁸, weak-field ligand (e.g. NiCl₄²⁻, NiBr₄²⁻)Tetrahedral, paramagneticSmall Δo\Delta_o for weak first-row ligand field; pairing-energy cost of square-planar not repaid
n=4,5n=4,5 (Pd, Pt) d⁸, virtually any ligandSquare-planar essentially without exceptionBoth larger Δo\Delta_o and smaller PP favor square-planar strongly and robustly

The chemist behind it

Origins · 1 minHans Bethe and how this was discovered
HB
Hans Bethe
1906–2005 · Crystal field theory origin, 1929, extended through the 1930s–1950s by Van Vleck, Ilse & Hartmann, and Orgel for transition-metal spectroscopy

Hans Bethe's 1929 paper on the splitting of atomic terms in crystalline electric fields was written for a solid-state physics audience — describing how an ion's electronic terms split under the point-group symmetry of a crystalline lattice — and contained no chemistry at all. It was not until the 1950s, through the work of John Van Vleck (who clarified the balance between purely electrostatic crystal field theory and covalent/molecular-orbital contributions), and inorganic chemists such as Leslie Orgel, that Bethe's group-theoretical machinery was adapted wholesale into ligand field theory for transition-metal complexes, including the square-planar and tetrahedral limiting cases built directly from his original correlation diagrams.

Square-planar splitting itself is not attributed to a single discoverer in the way that, say, a named reaction is. It emerged as a natural corollary once Bethe's crystal-field group theory was combined with the specific geometric case of a D4hD_{4h} point group, and once inorganic chemists in the 1950s–60s (notably including Orgel's widely used qualitative treatments of transition-metal stereochemistry) explained the observed magnetism and color of square-planar Ni(II), Pd(II), Pt(II), and Au(III) complexes with exactly this orbital-splitting picture. The experimental anchor — diamagnetic square-planar Ni(CN)₄²⁻ versus paramagnetic tetrahedral NiCl₄²⁻, both d⁸ — became one of the standard textbook demonstrations that geometry, not merely electron count, determines magnetism in the d-block, cementing square-planar crystal field theory's place in the coordination chemistry canon.

Origins · 1 minhow this was discovered
History
- 1893 — Alfred Werner proposes coordination theory, including square-planar geometry as a distinct stereochemical possibility for four-coordinate complexes.
- 1929 — Hans Bethe publishes the group-theoretical treatment of ion term splitting in crystalline fields, the mathematical ancestor of all crystal field diagrams.
- 1930s–1940s — John Van Vleck develops the electrostatic/ligand-field synthesis and magnetic-susceptibility theory needed to connect splitting diagrams to measurable magnetism.
- 1950s–1960s — Leslie Orgel and coworkers popularize qualitative crystal/ligand field diagrams (including square-planar and tetrahedral limiting cases) for teaching and rationalizing transition-metal spectra and magnetism.
- 1965 — Barnett Rosenberg's serendipitous discovery of cisplatin's antitumor activity draws square-planar Pt(II) coordination chemistry into pharmacology, cementing its practical importance beyond pure inorganic theory.
Why does removing axial ligands from an octahedron split t2g and eg into four levels instead of keeping two levels?

Removing the two axial ligands breaks the full octahedral (Oh) symmetry down to tetragonal/square-planar (D4h) symmetry. Orbitals that were degenerate purely because of Oh symmetry are no longer required to stay degenerate once that symmetry is lowered. Group theory shows that in D4h, dz2d_{z^2} and dx2−y2d_{x^2-y^2} (formerly degenerate eg) transform as different irreducible representations (a1ga_{1g} and b1gb_{1g}) and split apart, while dxyd_{xy} (formerly part of t2g, now b2gb_{2g}) also separates from dxz/dyzd_{xz}/d_{yz} (which remain degenerate as ege_g in D4h). The physical driver is simply that different orbitals feel different amounts of remaining ligand repulsion once the axial ligands are gone.

Is square-planar splitting always exactly the same shape, or does it depend on the ligand?

The qualitative ordering — dxz,dyzd_{xz},d_{yz} lowest, dx2−y2d_{x^2-y^2} highest by a wide margin — is essentially universal for real square-planar d⁸ complexes. The relative position of dz2d_{z^2} versus dxyd_{xy}, however, is genuinely ligand-dependent: strong π-acceptor ligands (like CO or CN⁻) and strong π-donor ligands shift these two intermediate levels differently because π-bonding effects (absent from the pure σ-only correlation diagram) add extra stabilization or destabilization specifically to dxyd_{xy}, dxzd_{xz}, and dyzd_{yz}, which have the correct symmetry to engage in π-interactions with ligand p or π* orbitals, while dx2−y2d_{x^2-y^2} and dz2d_{z^2} are purely σ-antibonding in character.

Why don't all d8 ions become square-planar — why are tetrahedral d8 complexes like NiCl4²⁻ still common?

Square-planar geometry only pays off electronically when the ligand field is strong enough that the CFSE gained by opening the large dx²−y² gap outweighs the extra electron-pairing energy cost of forcing all eight electrons into four paired orbitals. For weak-field ligands (like chloride on first-row Ni2+), the pairing-energy penalty dominates, and the ion instead adopts tetrahedral geometry, which needs no forced pairing at all and also minimizes ligand-ligand steric repulsion. Square-planar therefore correlates strongly with strong-field ligands (CN−, phosphines, dimethylglyoximate) and with heavier 4d/5d d8 metals (Pd2+, Pt2+, Au3+), where Δo is intrinsically larger and pairing energy P intrinsically smaller.

How is square-planar splitting related to Jahn-Teller distortion?

Square-planar geometry is the extreme endpoint of a tetragonal elongation distortion mode — the same distortion mode responsible for Jahn-Teller distortion in octahedral complexes with unevenly occupied eg orbitals (classically Cu2+, d9, and high-spin Mn3+/Cr2+, d4). In those cases the distortion typically stops at a moderate elongation (a "4+2" coordination environment) because the electronic driving force, while real, is not as strong as it is for d8. For d8 ions with strong-field ligands, the driving force is strong enough that the distortion runs all the way to complete axial ligand loss, producing a genuinely four-coordinate square-planar complex rather than a merely distorted octahedron.

Do square-planar complexes ever have five-coordinate or six-coordinate character in reality?

Yes, frequently. Many "square-planar" d8 complexes retain weak, long-range axial interactions with solvent molecules, counter-ions, or adjacent molecules in the solid state (sometimes called "square-pyramidal" or "elongated octahedral/4+2" coordination). This is common for Ni(II) and less common but not absent for Pd(II) and Pt(II). These weak axial contacts are the chemical manifestation of the small residual stabilization the metal can gain by re-admitting some axial ligand density, and they are mechanistically important: axial approach of a fifth ligand is exactly the first step of associative ligand substitution at square-planar centers, which proceeds via a five-coordinate trigonal-bipyramidal transition state.

Why is Δsp comparable in magnitude to Δo for the same ligand, rather than much smaller (since there are fewer ligands)?

Naively, fewer ligands might suggest a weaker overall field and a smaller splitting. But the correlation-diagram analysis shows the opposite: removing the axial ligands does not weaken the in-plane field at dx²−y² at all — the four equatorial ligands are unchanged, and they alone are enough to push dx²−y² up to roughly the same distance above the barycenter that the full eg set reached in the octahedral case. What actually changes most dramatically is dz², which drops steeply because it loses its axial ligands, not dx²−y² which was always primarily an in-plane orbital. So Δsp (defined via dxy to dx²−y²) ends up close to, even slightly larger than, Δo for the same ligand set.

What experimental evidence distinguishes square-planar from tetrahedral d8 complexes besides magnetism?

Beyond the sharp diamagnetic-vs-paramagnetic distinction (unambiguous for d8: zero unpaired electrons only in strong-field square-planar), UV-Vis spectra differ characteristically: square-planar d8 complexes typically show three spin-allowed d-d bands (from the three distinct filled-orbital-to-dx²−y² transitions) often with much higher molar absorptivities and more structured splitting patterns than the single broad, more Laporte-forbidden-relaxed band typical of tetrahedral complexes. X-ray crystallography directly confirms the coordination geometry, and bond angle measurements (90° at the metal for square-planar versus 109.5° for tetrahedral) are unambiguous once a structure is solved.

Before this

Where this leads