Nephelauxetic Effect and Racah Parameters

Take a free Ni²⁺ ion, gaseous and alone, and shine light through it. You'll find a transition around 8500–15000 cm⁻¹. This transition reflects electron-electron repulsion energies fixed by pure atomic physics. Now bind that same Ni²⁺ into [Ni(H₂O)₆]²⁺, then into [NiCl₄]²⁻, then into [Ni(CN)₄]²⁻. Repeat the measurement each time. Every single transition energy shifts. This happens not just because the crystal field splits the d orbitals differently in each case. It also happens because the electrons themselves seem to repel each other less once the metal is surrounded by ligands. This shrinkage of interelectron repulsion is quantified through the Racah parameters B and C. It is called the nephelauxetic effect — literally "cloud-expanding." It is one of the cleanest experimental windows chemists have onto covalency in coordination compounds. The effect matters because Interpreting Spectra with Tanabe-Sugano Diagrams only works quantitatively if you know which value of B to plug in. That value is never the free-ion value. A crystal-field treatment using free-ion B systematically overestimates transition energies. It also underestimates the degree of metal-ligand orbital mixing. The nephelauxetic effect converts this apparent "failure" of pure ionic crystal-field theory into a quantitative covalency probe. The more B (and C) shrink relative to the free ion, the more the metal d-electron cloud has expanded onto ligand orbitals. This means the bond is more covalent. This bridges the gap between the electrostatic language of crystal-field theory and the orbital-mixing language of ligand-field/molecular-orbital theory. The big picture is that two independent ligand-ordering series exist side by side. The Spectrochemical Series and Δo ranks ligands by how strongly they split d orbitals. This is a σ/π-bonding, largely orbital-overlap argument. The nephelauxetic series ranks ligands by how much they shrink interelectron repulsion. This is a covalency/softness argument tied to how easily ligand electron density delocalizes onto the metal. The two series are not the same. Their partial independence is itself chemically informative. F⁻ is a strong nephelauxetic-series "hard" ligand producing little cloud expansion. Yet it sits in the middle of the spectrochemical series. Br⁻ and I⁻ sit low in Δo but high in nephelauxetic reduction. This is because soft, polarizable halides delocalize metal d-density efficiently. This article builds the Racah-parameter framework from atomic multiplet theory. It derives the nephelauxetic ratio β and shows how it decomposes into ligand and metal contributions. It works through three fully quantitative examples: a free-ion baseline, a real octahedral complex, and a competition-level multi-parameter fit. It also maps out where these ideas surface today. This ranges from interpreting d-d spectra of biological metal centers to calibrating covalency in actinide and lanthanide separations chemistry.

Key takeaways
  • Interelectron repulsion in a free ion is captured by three Racah parameters A, B, C; only B and C matter for relative term energies.
  • Coordinating a metal ion to ligands always shrinks B and C relative to the free ion — the nephelauxetic effect — because metal d-density delocalizes onto ligands (covalency), reducing d-d repulsion.
  • The nephelauxetic ratio β=Bcomplex/Bfree ion\beta = B_{complex}/B_{free\ ion} quantifies this shrinkage; smaller β means more covalent bonding.
  • Ligands and metals each have their own empirical ranking; 1−β≈hligand×kmetal1-\beta \approx h_{ligand} \times k_{metal} (Jørgensen's equation) separates the two contributions.
  • The nephelauxetic series (roughly F⁻ < H₂O < NH₃ < en < ox²⁻ < NCS⁻ < Cl⁻ < CN⁻ < Br⁻ < I⁻) tracks ligand softness/polarizability, not Δo.
  • Key idea: B and C shrink upon complexation because covalent delocalization of metal d-electrons onto ligand orbitals lowers electron-electron repulsion, and the size of that shrinkage (β) is a direct spectroscopic measure of covalency independent of crystal-field splitting.

Where this lives today

Nephelauxetic analysis is not a historical curiosity; it remains the standard way chemists extract a quantitative covalency number directly from an absorption spectrum, with no need for computational modeling. In lanthanide and actinide separations chemistry — critically important for nuclear fuel reprocessing and rare-earth recycling — the central practical question is whether a given ligand binds An³⁺ ions more covalently than Ln³⁺ ions of similar ionic radius, since that small covalency difference is exploited to separate chemically similar trivalent actinides from lanthanides using soft N- or S-donor extractants. Nephelauxetic ratios measured from f-f transitions (though small in magnitude for these contracted orbitals) provide direct spectroscopic evidence for exactly this differential covalency, complementing computational bonding analyses. In bioinorganic spectroscopy, d-d and f-f transitions of metal centers in metalloproteins (Cu, Fe, Co, Mn active sites; Ln³⁺ substituted into Ca²⁺ sites as luminescent probes) are routinely fit with Racah-parameter/Tanabe-Sugano analysis, and the extracted β reports on how covalent the protein's native donor set (carboxylate O, imidazole N, thiolate S) is relative to simple aqua or ammine model complexes — informative for understanding catalytic activation of substrates at these centers. In materials science, the nephelauxetic effect underlies the tuning of phosphor and laser-gain-medium emission colors: Cr³⁺-doped oxide and fluoride hosts (ruby, alexandrite, emerald) show R-line and broadband emission energies that shift systematically with host covalency exactly as predicted by β and the C/B-dependent formula for the ²E state used in Example 3, and this same physics underlies modern rare-earth-doped LED phosphor engineering, where nephelauxetic shifts of f-electron transitions are used to fine-tune emission color. Finally, in fundamental ligand-field theory pedagogy and modern computational benchmarking, β remains a simple, experimentally accessible number against which more sophisticated density-functional or ligand-field molecular-orbital calculations of covalency (spin densities, natural bond orbital delocalization indices, X-ray absorption pre-edge intensities) are cross-checked, precisely because it requires nothing more than an optical absorption spectrum and a Tanabe-Sugano diagram to obtain.

Ask the tutor: For a given metal, why does the nephelauxetic ratio β not correlate more strongly with the ligand's position in the spectrochemical series, and can you show me a second worked case where β and Δo trends actually run in opposite directions?

The physical picture

In a free transition-metal ion, the energies of terms arising from a dⁿ configuration (²D, ³F, ³P, ¹G, ¹S, etc., for d²) are fixed entirely by electron-electron repulsion within the d shell — there is no external field to break the spherical symmetry. Racah showed that all these repulsion energies can be written as linear combinations of just three parameters, A, B, and C, each with units of energy (typically cm⁻¹). A shifts every term by the same amount (an overall constant that never affects relative spacings, and is usually dropped from spectroscopic analysis). B and C carry all the physically meaningful information: B measures repulsion between electrons in different orbitals with parallel-spin correlation effects, and C is sensitive mostly to the pairing energy contribution from electrons of opposite spin in the same orbital. The ratio C/B is roughly constant across a row of the periodic table (empirically close to 4–4.5 for first-row transition metals), which is why many treatments track B alone as a proxy for the whole repulsion manifold. When the free ion is placed in a ligand field, the d orbitals split (this splitting is the Δo of the spectrochemical series and is an orbital energy effect, not a repulsion effect). But something else happens simultaneously: the metal d orbitals mix with ligand orbitals to form molecular orbitals that are partly ligand in character. An electron that used to live entirely on the metal now spends part of its time on the ligands. Two electrons that would have strongly repelled each other while both confined to the compact metal 3d shell now repel each other less, because their combined probability density has been spread — "expanded" — over a much larger volume that includes ligand atoms. This physical delocalization is the microscopic origin of the nephelauxetic effect: the "electron cloud" (Greek nephelē, cloud) genuinely expands (auxesis, growth) in space, lowering the effective repulsion parameters B and C measured spectroscopically. Two limiting descriptions of this delocalization are useful. In the central-field covalency picture, the metal-ligand bond is treated as still largely ionic, but the metal's own d orbitals become more diffuse (radially expanded) because the effective nuclear charge felt by the d electrons is screened differently once ligand electron density surrounds the ion — this is sometimes framed as a reduction in effective charge. In the symmetry-restricted (molecular orbital) covalency picture, true antibonding metal-ligand molecular orbitals form, with the metal d-based MOs carrying appreciable ligand orbital character (mixing coefficients from group-theoretically allowed σ and π overlaps). Both pictures predict the same experimental signature — reduced B and C — and in practice both mechanisms operate together; the nephelauxetic ratio does not distinguish them without additional data (e.g., EPR hyperfine coupling constants, which are more diagnostic of MO mixing coefficients).

Figure
Free ion vs. complex: d-electron cloud expansion and its spectroscopic signature
Interelectron repulsion shrinks upon complexation Free ion M³⁺(g) compact 3d cloud B = B(free ion) ligands bind Octahedral complex ML₆ B = β·B(free ion), β < 1 d-density delocalizes onto ligand orbitals → less d–d repulsion → smaller B, C

The measurable consequence is that every spin-allowed and spin-forbidden d-d band in a complex sits at a different energy than the naive free-ion Racah parameters would predict when combined with the observed Δo. Historically this was noticed as an anomaly — complexes' spectra "didn't fit" pure crystal-field predictions using free-ion B — and Jørgensen's systematic cataloguing of the shrinkage, ligand by ligand and metal by metal, turned the anomaly into the nephelauxetic series, a predictive covalency scale.

Figure
Nephelauxetic series and spectrochemical series compared — different orderings
Two independent ligand rankings Nephelauxetic series (h, covalency ↑) F⁻ < H₂O < NH₃ < en < ox²⁻ < NCS⁻ < Cl⁻ < CN⁻ < Br⁻ < I⁻ Spectrochemical series (Δo ↑) I⁻ < Br⁻ < Cl⁻ < F⁻ < H₂O < NH₃ < en < NO₂⁻ < CN⁻ F⁻ and I⁻ swap relative rank: soft, polarizable ligands expand the d-cloud strongly but split it only weakly while F⁻ resists covalent delocalization yet still gives a moderate Δo through σ-donation.

The core ideas

The Racah parameterization writes free-ion term energies for a dⁿ configuration as linear combinations of A, B, C. For example, in a d² configuration the ground and excited free-ion terms include

E(3F)=A−8B,E(3P)=A+7B,E(1D)=A−3B+2C,E(1G)=A+4B+2CE(^3F) = A - 8B, \qquad E(^3P) = A + 7B, \qquad E(^1D) = A -3B + 2C, \qquad E(^1G) = A + 4B + 2C

where:

  • AA is an overall (orbital-independent) repulsion constant, dropped in relative-energy comparisons,
  • BB is the Racah parameter controlling repulsion between electrons with parallel-spin correlation ("interorbital" repulsion, sensitive to covalency),
  • CC is the Racah parameter controlling pairing/exchange energy contributions important for spin-multiplicity changes.

The central defining relation of the nephelauxetic effect is the nephelauxetic ratio,

β=BcomplexBfree ion,(often also β′=CcomplexCfree ion)\beta = \frac{B_{complex}}{B_{free\ ion}}, \qquad \left(\text{often also } \beta' = \frac{C_{complex}}{C_{free\ ion}}\right)

where:

  • BcomplexB_{complex} is the Racah B extracted by fitting observed d-d transition energies to a Tanabe-Sugano diagram for the complex,
  • Bfree ionB_{free\ ion} is the well-tabulated Racah B for the gas-phase Mⁿ⁺ ion (from atomic emission spectroscopy),
  • β\beta (and analogously β′\beta') is always ≤1\le 1 for real complexes, with smaller values signaling greater covalency.

Reading: β = 1 would mean a perfectly ionic complex where the metal d electrons behave exactly as they do in vacuum — never observed. β approaching, say, 0.3–0.5 (seen for very covalent, soft-donor complexes) means the effective repulsion the d electrons feel has been cut by half or more, a huge delocalization effect that no purely electrostatic point-charge model can produce. Jørgensen's empirical separation of ligand and metal contributions to this shrinkage is:

1−β=hligand⋅kmetal1-\beta = h_{ligand} \cdot k_{metal}

where:

  • hligandh_{ligand} is a dimensionless empirical parameter assigned to each ligand (the nephelauxetic series is a ranking of hh values, with h(F−)h(F^-) small and h(I−),h(CN−)h(I^-), h(CN^-) large),
  • kmetalk_{metal} is a dimensionless empirical parameter assigned to each metal ion (larger for more polarizable, more covalently inclined metals — e.g., kk increases going down a triad and with formal oxidation state trends),
  • the product form assumes ligand and metal contributions are separable and additive across the different ligands in a mixed complex (an approximation, analogous to the additivity rules for Δo).

Reading: this equation lets you predict β for an unmeasured complex from tabulated h and k values measured on other complexes of the same metal or same ligand — exactly analogous to using Δo additivity rules to predict crystal-field splittings for mixed-ligand complexes.

Derivation AFrom atomic multiplet theory to the Racah B parameter

Goal: show why interelectron repulsion energies in a dⁿ free ion collapse to just three numbers (A, B, C), and why B in particular is the parameter most sensitive to radial electron distribution (hence to covalent delocalization).

Step 1 — the repulsion Hamiltonian. For an n-electron atom/ion, beyond the one-electron (kinetic + nuclear attraction) terms, the electron-electron repulsion operator is

Hee=∑i<je24πε0rijH_{ee} = \sum_{i<j} \frac{e^2}{4\pi\varepsilon_0 r_{ij}}

Within a single dⁿ configuration (all electrons confined to nl = 3d orbitals), matrix elements of HeeH_{ee} between Slater determinants can be expanded in terms of two-electron Coulomb integrals F(0),F(2),F(4)F^{(0)}, F^{(2)}, F^{(4)} and exchange-type integrals, using the Slater-Condon parameterization for equivalent electrons of the same l=2l=2.

Step 2 — reduction to three independent radial integrals. For d electrons, angular momentum coupling (Wigner 3j/6j algebra applied to l=2l=2) shows that only three independent radial Slater integrals survive: F0,F2,F4F_0, F_2, F_4 (all positive, all decreasing with increasing radial diffuseness of the d orbital, since larger orbitals have smaller average 1/r121/r_{12}). Racah found that specific linear combinations of these three integrals appear repeatedly and universally across all d² through d⁸ term energies, and defined

A=F0−49F4,B=F2−5F4,C=35F4A = F_0 - 49F_4, \qquad B = F_2 - 5F_4, \qquad C = 35F_4

so that every term energy for every dⁿ configuration can be written purely in A, B, C (never needing F₂, F₄ individually). This is a change of basis, not an approximation — it is exact within the dⁿ configuration-interaction problem.

Step 3 — why B (not A) tracks covalency. A is dominated by F0F_0, the spherically averaged repulsion, which shifts all terms by the same amount and therefore never appears in observed transition energies (spectroscopy only measures energy differences between terms). B and C, built from F2F_2 and F4F_4, measure the angular/radial shape dependence of repulsion — specifically, how repulsion differs between terms with different orbital occupation patterns (e.g., ³F vs. ³P in d²). Because F2F_2 and F4F_4 scale as

Fk∝∫ ⁣ ⁣∫R3d(r1)2 r< kr> k+1 R3d(r2)2 r12r22 dr1 dr2F_k \propto \int\!\!\int R_{3d}(r_1)^2 \, \frac{r_<^{\,k}}{r_>^{\,k+1}} \, R_{3d}(r_2)^2 \, r_1^2 r_2^2\, dr_1\, dr_2

they are extremely sensitive to how compact the radial wavefunction R3d(r)R_{3d}(r) is. Any physical process that expands R3d(r)R_{3d}(r) radially outward (screening changes, or true delocalization onto ligand orbitals) reduces both F2F_2 and F4F_4, hence reduces B and C. This is the direct microscopic link between "electron cloud expansion" and the shrinkage of the spectroscopically extracted B parameter.

Step 4 — connecting to the complex. In a ligand field, the true metal-based molecular orbitals are linear combinations ψ=cdϕd−cLϕL\psi = c_d\phi_d - c_L\phi_L (antibonding combination, normalized, with cd2+cL2−2cdcLS≈1c_d^2 + c_L^2 - 2c_dc_L S \approx 1 for small ligand-metal overlap SS). The effective d-orbital radial density used in the repulsion integrals is no longer pure R3d(r)2R_{3d}(r)^2 but cd2R3d(r)2c_d^2 R_{3d}(r)^2 plus ligand-centered density — spread over a much larger effective volume. Since cd2<1c_d^2 < 1 whenever there is any covalent mixing, the effective F2,F4F_2, F_4 (and hence B, C) computed from this delocalized density are strictly smaller than the free-ion values. This is why β < 1 is essentially guaranteed for any real bond with nonzero covalent mixing, and why β is often used as an empirical proxy for cd2c_d^2 (though the quantitative relationship requires care, since β also depends on how the ligand density itself is distributed).

Go deeper · Derivation B · 2 minThe Jørgensen additivity relation and extracting h, k values from data
Derivation BThe Jørgensen additivity relation and extracting h, k values from data

Goal: derive how a single measured β for one metal-ligand combination, combined with data across a series of complexes, yields separable ligand parameters h and metal parameters k, and show how to use these to predict an unmeasured β.

Step 1 — set up the multiplicative ansatz. Jørgensen proposed, empirically (not derived from first principles, but strongly supported by the data), that the fractional reduction in B factorizes as

1−βML=hL⋅kM1 - \beta_{ML} = h_L \cdot k_M

for a metal M in ligand environment L. This is directly analogous in spirit to assuming Δo factorizes into ligand and metal contributions in the spectrochemical additivity rule — both are mean-field-like separability assumptions that work well because the physical origin (σ/π covalent mixing) is roughly separable between "how good a covalent partner is this ligand" and "how covalently inclined is this metal."

Step 2 — anchor the scale. The h and k scales each need one reference point since only their product is measured. By convention, k(Co3+)≡1k(\text{Co}^{3+}) \equiv 1 (or in some tabulations Ni2+\text{Ni}^{2+}, conventions vary by source) is fixed as a reference metal, which then fixes the absolute scale of all hLh_L values from measurements of Co3+L6\text{Co}^{3+}L_6 complexes (well-studied, high-spin-free d⁶ low-spin complexes with clean, well-resolved spin-allowed bands). Once hLh_L values are pinned down this way, measuring β for any other metal with any of these same ligands gives kMk_M directly.

Step 3 — solve for an unknown k or h. Suppose βobs\beta_{obs} is measured for metal M with ligand L, and hLh_L is already tabulated (from Co³⁺ data). Then

kM=1−βobshLk_M = \frac{1-\beta_{obs}}{h_L}

Once kMk_M is known for that metal, it can be combined with any other tabulated hL′h_{L'} to predict

βpredicted=1−hL′⋅kM\beta_{predicted} = 1 - h_{L'}\cdot k_M

for a different ligand L′ bound to the same metal M — without ever measuring that new complex's spectrum.

Step 4 — worked numeric check of internal consistency. Take illustrative tabulated values h(F−)=0.8h(\mathrm{F^-}) = 0.8, h(H2O)=1.0h(\mathrm{H_2O}) = 1.0, h(NH3)=1.4h(\mathrm{NH_3}) = 1.4, h(Cl−)=2.0h(\mathrm{Cl^-}) = 2.0 (Jørgensen's original units, dimensionless, scaled so h(H2O)=1.0h(\mathrm{H_2O})=1.0 by definition) and k(Co3+)=1.0k(\mathrm{Co^{3+}}) = 1.0 (reference), k(Cr3+)≈0.21k(\mathrm{Cr^{3+}}) \approx 0.21, k(Ni2+)≈0.12k(\mathrm{Ni^{2+}}) \approx 0.12 (illustrative, order-of-magnitude values consistent with the literature range for these ions). For [Cr(H2O)6]3+[\mathrm{Cr(H_2O)_6}]^{3+}:

1−β=hH2O⋅kCr3+=1.0×0.21=0.21⇒β≈0.791-\beta = h_{H_2O}\cdot k_{Cr^{3+}} = 1.0 \times 0.21 = 0.21 \quad\Rightarrow\quad \beta \approx 0.79

This is consistent (order of magnitude and direction) with the experimentally derived β ≈ 0.8 typically quoted for [Cr(H2O)6]3+[\mathrm{Cr(H_2O)_6}]^{3+}, confirming that the multiplicative ansatz, even though empirical, captures the real chemistry: aqua ligands are weak nephelauxetic agents (small h), and Cr³⁺ is a relatively "hard," electrostatically-bonding metal center (small k) compared to, say, Co³⁺ or Ni²⁺ with softer character or more diffuse valence d orbitals.

Step 5 — limiting case. As the ligand field becomes purely ionic in the limiting sense (hL→0h_L \to 0, an idealized "non-covalent" ligand that donates no electron density into metal-centered antibonding orbitals), β→1\beta \to 1 for every metal — the nephelauxetic effect vanishes identically, exactly as required, since with zero orbital mixing the metal d-electron cloud is untouched by the ligand field and B must equal the free-ion value. No real ligand achieves h=0h=0, but F⁻ and H₂O approach this limit most closely among common ligands, consistent with their position at the "hard" end of the nephelauxetic series.

Worked examples

Example 1 (foundational) — Extracting β for a real octahedral Ni²⁺ complex from raw transition energies. The complex [Ni(H2O)6]2+[\mathrm{Ni(H_2O)_6}]^{2+} (high-spin d⁸, Oh symmetry) shows three spin-allowed d-d bands (using the Tanabe-Sugano diagram for d⁸, ground term ³A₂g) typically reported at:

ν1(3A2g→3T2g)≈8500 cm−1,ν2(3A2g→3T1g(F))≈13800 cm−1,ν3(3A2g→3T1g(P))≈25300 cm−1\nu_1(^3A_{2g}\to{}^3T_{2g}) \approx 8500\ \text{cm}^{-1}, \quad \nu_2(^3A_{2g}\to{}^3T_{1g}(F)) \approx 13800\ \text{cm}^{-1}, \quad \nu_3(^3A_{2g}\to{}^3T_{1g}(P)) \approx 25300\ \text{cm}^{-1}

For d⁸ in Oh, the standard Tanabe-Sugano energy expressions (in the weak-to-moderate field regime relevant here) give, using Δo\Delta_o directly from ν1\nu_1 since 3T2g^3T_{2g} is the first excited term with energy exactly Δo\Delta_o above the ³A₂g ground term in this configuration:

Δo=ν1=8500 cm−1\Delta_o = \nu_1 = 8500\ \text{cm}^{-1}

The ratio ν2/ν1\nu_2/\nu_1 and ν3/ν1\nu_3/\nu_1 can then be matched against tabulated Tanabe-Sugano ratios (or read off the diagram at Δo/B≈9\Delta_o/B \approx 9–1010, typical for d⁸ aqua complexes) to extract B. Using the standard d⁸ energy matrix result for the ³T₁g(F) term,

E(3T1g(F))=12(15B+8Dq±(9B−8Dq)2+4⋅6B⋅8Dq )(plus a constant),E(^3T_{1g}(F)) = \tfrac{1}{2}\left(15B+8Dq \pm\sqrt{(9B-8Dq)^2 + 4\cdot 6B\cdot 8Dq}\,\right)\quad(\text{plus a constant}),

fitting the three observed energies self-consistently (the standard graphical/tabular method) yields

Bcomplex≈940 cm−1B_{complex} \approx 940\ \text{cm}^{-1}

Comparing to the free-ion value for Ni²⁺, Bfree=1080 cm−1B_{free} = 1080\ \text{cm}^{-1} (tabulated from gas-phase atomic spectroscopy):

β=9401080=0.87\beta = \frac{940}{1080} = 0.87

Reading: β = 0.87 means water reduces Ni²⁺'s interelectron repulsion by only 13%, consistent with H₂O sitting low in the nephelauxetic series (a largely σ-donor, weakly covalent, "hard" ligand). This is the expected baseline against which more covalent ligands are compared.

Example 2 (applied) — Comparing β across a real nephelauxetic series for Ni²⁺ and identifying the covalency trend. Using the same Ni²⁺ free-ion reference Bfree=1080 cm−1B_{free} = 1080\ \text{cm}^{-1}, tabulated complex B values (extracted the same way as Example 1, from Tanabe-Sugano fits to experimental spectra) are approximately:

B[NiF6]4−≈1010 cm−1,B[Ni(H2O)6]2+≈940 cm−1,B[Ni(NH3)6]2+≈890 cm−1,B[Ni(CN)4]2− (Ni(0)-like covalency, illustrative)B[\mathrm{NiF_6}]^{4-} \approx 1010\ \text{cm}^{-1}, \quad B[\mathrm{Ni(H_2O)_6}]^{2+} \approx 940\ \text{cm}^{-1}, \quad B[\mathrm{Ni(NH_3)_6}]^{2+} \approx 890\ \text{cm}^{-1}, \quad B[\mathrm{Ni(CN)_4}]^{2-}\ (\text{Ni(0)-like covalency, illustrative})

Computing β for the first three (Ni²⁺, Oh, comparable coordination environments):

βF−=10101080=0.935,βH2O=9401080=0.870,βNH3=8901080=0.824\beta_{F^-} = \frac{1010}{1080} = 0.935, \qquad \beta_{H_2O} = \frac{940}{1080} = 0.870, \qquad \beta_{NH_3} = \frac{890}{1080} = 0.824

Reading: β decreases monotonically F⁻ > H₂O > NH₃ — exactly the order predicted by the nephelauxetic series (F⁻ least cloud-expanding, NH₃ more so among these three). Note this order does not match the spectrochemical series ranking of Δo for these same three ligands in terms of relative spacing — NH₃ gives a noticeably larger Δo jump than H₂O relative to F⁻ than the modest β changes would suggest — demonstrating explicitly that Δo (orbital splitting, dominated by σ-donor strength and orbital energy match) and B-shrinkage (interelectron repulsion reduction, dominated by covalent delocalization/softness) are physically distinct, only loosely correlated quantities. A ligand can be a strong σ-donor (raising Δo a lot) while being a comparatively poor nephelauxetic agent, or vice versa, because the two effects probe different aspects of the metal-ligand bond (orbital energy overlap vs. radial delocalization/polarizability).

Example 3 (competition-level) — Using the Jørgensen relation to predict β for a mixed-ligand complex and cross-check with C/B ratio consistency. Consider [Cr(NH3)4Cl2]+[\mathrm{Cr(NH_3)_4Cl_2}]^{+} (d³, Cr³⁺, trans or cis isomer — assume approximately Oh-like average field for this estimate), and suppose you are given (illustrative literature-consistent values):

h(NH3)=1.4,h(Cl−)=2.0,k(Cr3+)=0.21h(\mathrm{NH_3}) = 1.4, \quad h(\mathrm{Cl^-}) = 2.0, \quad k(\mathrm{Cr^{3+}}) = 0.21

Step A — additivity for a mixed-ligand field. Treat the effective ligand nephelauxetic parameter as a coordination-number-weighted average (an extension of the additivity idea, exactly parallel to Δo additivity for mixed-ligand octahedra):

heff=4 h(NH3)+2 h(Cl−)6=4(1.4)+2(2.0)6=5.6+4.06=9.66=1.60h_{eff} = \frac{4\,h(\mathrm{NH_3}) + 2\,h(\mathrm{Cl^-})}{6} = \frac{4(1.4) + 2(2.0)}{6} = \frac{5.6+4.0}{6} = \frac{9.6}{6} = 1.60

Step B — predict β.

1−β=heff⋅kCr3+=1.60×0.21=0.336⇒βpredicted≈0.6641-\beta = h_{eff}\cdot k_{Cr^{3+}} = 1.60 \times 0.21 = 0.336 \quad\Rightarrow\quad \beta_{predicted} \approx 0.664

Step C — predict B(complex) and cross-check against a plausible C/B ratio. Using Bfree(Cr3+)=918 cm−1B_{free}(\mathrm{Cr}^{3+}) = 918\ \text{cm}^{-1} (tabulated free-ion value):

Bcomplex,predicted=β×Bfree=0.664×918 cm−1=609.6 cm−1B_{complex,predicted} = \beta \times B_{free} = 0.664 \times 918\ \text{cm}^{-1} = 609.6\ \text{cm}^{-1}

For consistency, C/B for Cr³⁺ complexes typically runs close to C/B≈4.5C/B \approx 4.5 (slightly ligand-dependent but roughly constant, as expected since both B and C shrink together via the same delocalization mechanism, keeping their ratio comparatively stable — a well-known empirical regularity). This predicts

Ccomplex,predicted≈4.5×609.6 cm−1=2743 cm−1C_{complex,predicted} \approx 4.5 \times 609.6\ \text{cm}^{-1} = 2743\ \text{cm}^{-1}

which can then be used directly as an input to the d³ Tanabe-Sugano diagram to predict the position of the spin-forbidden 4A2g→2Eg^4A_{2g}\to{}^2E_g ("ruby line" analog) transition, since that transition's energy in d³ is given to a good approximation by E(2Eg)≈9B+3CE(^2E_g) \approx 9B + 3C (nearly field-independent, which is precisely why this transition is famously sharp and used for ruby R-line thermometry). Plugging in:

E(2Eg)≈9(609.6)+3(2743)=5486.4+8229=13715 cm−1E(^2E_g) \approx 9(609.6) + 3(2743) = 5486.4 + 8229 = 13715\ \text{cm}^{-1}

Reading: this predicted value (≈ 13700 cm⁻¹, roughly 730 nm, deep red) is in the right ballpark for Cr³⁺ spin-forbidden emission/absorption lines in mixed-ligand ammine/chloro complexes, illustrating how the nephelauxetic framework lets you predict an entire spin-forbidden transition energy from tabulated h, k values alone, without ever measuring the actual complex — a genuinely predictive, falsifiable use of the theory. (Real ruby, Cr³⁺ in Al₂O₃, is a different — oxide — ligand environment and gives the famous R-lines near 14400 cm⁻¹/694 nm; the calculation here is for the illustrative ammine/chloro complex, not ruby itself, and the moderate difference is consistent with the different, more ionic oxide ligand field in Al₂O₃ having a larger β than the softer mixed N/Cl donor set used here.)

Limiting case check: if instead all six ligands were NH₃ (pure hexaammine, heff=1.4h_{eff}=1.4), 1−β=1.4×0.21=0.2941-\beta = 1.4\times0.21 = 0.294, β=0.706\beta = 0.706 — larger than the mixed Cl/NH₃ value of 0.664, correctly reflecting that replacing two NH₃ with the more nephelauxetic Cl⁻ increases covalent cloud expansion (lowers β further) exactly as expected from Cl⁻'s higher position in the nephelauxetic series.

Limiting cases, scope & common traps

SituationWhat happens to βWhy
Purely ionic, point-charge limit (hypothetical)β→1\beta \to 1No orbital mixing ⇒ no d-cloud delocalization ⇒ B unchanged from free ion
Very soft, highly polarizable ligand (I⁻, CN⁻ via π-backbonding channels)β\beta small (strong shrinkage)Large covalent mixing coefficient cLc_L delocalizes d-density extensively
Free ion (g), no ligands at allβ≡1\beta \equiv 1 by definitionThis is the reference state that defines the ratio
Comparing Δo and β for the same ligand seriesOften uncorrelated rankingsΔo reflects orbital energy/overlap matching (σ, π*); β reflects radial delocalization/softness — different physics
C/B ratio across many complexes of one metalRoughly constant (~4–4.5, metal-dependent)Both B and C shrink via the same underlying delocalization mechanism
Lanthanide/actinide f-block ionsβ close to 1 (small nephelauxetic effect)Contracted, core-like 4f orbitals barely engage in covalent mixing; nephelauxetic effect is a much smaller, more subtle probe here, historically important for detecting any covalency in f-block bonding
Using free-ion B in a Tanabe-Sugano fit to a real complexSystematic disagreement/failure to fit all bands simultaneouslyThe correct procedure always requires an adjustable, reduced BcomplexB_{complex} (equivalently, fitting β) — this is precisely why nephelauxetic corrections exist

The chemist behind it

Origins · 1 minChristian Klixbüll Jørgensen and how this was discovered
CK
Christian Klixbüll Jørgensen
1931–2001 · Danish inorganic and theoretical chemist working at Geneva and CERN-adjacent institutions in the mid-20th century boom of ligand-field theory

Jørgensen was the central systematizer of ligand-field spectroscopy in the 1950s-60s, working alongside contemporaries like Carl Ballhausen. He compiled enormous tables of experimentally extracted Racah parameters across essentially the entire transition series and lanthanide/actinide series, and it was Jørgensen who coined the term "nephelauxetic effect" and formalized the empirical h/k separability rule that bears his name. His broader project was to push beyond pure electrostatic crystal-field theory toward a semi-empirical, spectroscopically calibrated picture of covalency — arguing, controversially at the time, that even "ionic" complexes show measurable, quantifiable covalent character, using the nephelauxetic ratio as the primary evidence. He also made foundational contributions to charge-transfer spectroscopy and the systematics of lanthanide/actinide electronic structure, where nephelauxetic effects are subtle but diagnostically crucial for detecting rare covalent participation of f orbitals.

History
- 1943 — Giulio Racah publishes the parameterization of atomic term energies for equivalent electrons in terms of A, B, C (Racah parameters), building on Slater-Condon theory.
- 1954-1955 — Ballhausen and Jørgensen begin systematic ligand-field analyses of transition-metal complex spectra, noting that free-ion B values fail to fit observed spectra.
- 1962 — Jørgensen publishes his comprehensive treatment establishing the nephelauxetic series and the multiplicative 1−β=h⋅k1-\beta = h\cdot k relation across dozens of metal-ligand combinations.
- 1960s-1970s — Nephelauxetic analysis extended systematically to lanthanide and actinide complexes, becoming a key early spectroscopic tool for detecting subtle covalency in f-block bonding.
- Present — β and Racah-parameter fitting remain a standard quantitative covalency diagnostic taught alongside Tanabe-Sugano analysis in every modern inorganic spectroscopy curriculum.
Why is the nephelauxetic effect called "cloud-expanding" specifically, rather than just "covalency reduction"?

The name emphasizes the physical, spatial mechanism rather than just labeling the outcome. When metal d electrons delocalize onto ligand orbitals through covalent mixing, the electron probability density genuinely occupies a larger volume of space than it would in the free ion — the "cloud" of electron density literally expands outward toward and onto the ligands. This spatial expansion is what reduces the average 1/r121/r_{12} repulsion between d electrons, which is measured as a smaller Racah B (and C). Jørgensen chose the Greek-rooted name (nephelē = cloud, auxēsis = growth/increase) specifically to keep the mechanistic, spatial picture in view, distinguishing it from other, more diffuse uses of the word "covalency" that don't specify a concrete physical origin.

Can the nephelauxetic effect ever make β greater than 1?

No, not for genuine metal-ligand bonding interactions — β ≤ 1 always, by the physics of the effect. Delocalizing d-electron density onto ligand orbitals can only spread the electron cloud out and reduce average repulsion; there is no mechanism by which coordinating ligands would compress the metal d-electron density and increase repulsion above the free-ion value. Reported β values occasionally appear very close to 1 (weakly covalent complexes with hard ligands like F⁻) or even nominally slightly above 1 in some tabulations, but such cases reflect experimental/fitting uncertainty in extracting B from a limited number of resolved bands, not genuine physics — they should be treated as β ≈ 1 (negligible nephelauxetic effect) rather than as evidence of "anti-nephelauxetic" behavior.

How is B actually extracted from an experimental spectrum in practice?

You start by identifying the observed spin-allowed (and ideally some spin-forbidden) d-d band positions and assigning each to a term predicted by the relevant Tanabe-Sugano diagram for that dⁿ configuration. Because Tanabe-Sugano diagrams plot E/B against Δo/B, you scan across trial values of Δo/B until the ratios between your observed transition energies match the ratios predicted by the diagram at some specific Δo/B value; that matching value of Δo/B, combined with the known observed Δo (typically read directly off the lowest-energy spin-allowed band), gives B directly. Modern practice often does this by least-squares fitting the full set of observed bands simultaneously against the exact Racah-parameter energy-matrix expressions (solving the appropriate secular determinant) rather than reading a printed diagram by eye, giving more precise, self-consistent B and C values along with quoted uncertainties.

Does the nephelauxetic effect apply to tetrahedral and other non-octahedral geometries too?

Yes — the nephelauxetic effect is a consequence of covalent metal-ligand orbital mixing generally, and is not tied to any particular coordination geometry. Tetrahedral, square-planar, and lower-symmetry complexes all show reduced B and C relative to the free ion, extracted using the Tanabe-Sugano diagram appropriate to that geometry and dⁿ configuration (noting that a tetrahedral dⁿ complex behaves spectroscopically like an octahedral d^(10-n) complex by the standard hole/electron formalism, up to the reduced Δt ≈ (4/9)Δo scaling). The same h/k tabulated values are often transferable across geometries to reasonable approximation, since h and k mainly encode ligand softness and metal covalency-preference rather than geometric details, though small geometry-dependent corrections do appear in high-precision work.

Why do soft, polarizable ligands like I⁻ and CN⁻ sit high in the nephelauxetic series despite occupying opposite ends of the spectrochemical series?

Both properties trace back to how easily a ligand's electron density mixes with metal d orbitals, but they weight different aspects of that mixing. I⁻ has large, diffuse, polarizable filled p orbitals that overlap substantially with metal d orbitals even though the resulting σ/π* splitting (Δo) stays modest because the ligand orbital energies are poorly matched to the metal d orbital energies. CN⁻, by contrast, is both a strong π-acceptor (large Δo via metal-to-ligand π-backbonding stabilizing the t2g set) and a good covalent partner overall (strong nephelauxetic effect), so it happens to rank high on both series simultaneously. The general lesson is that strong covalent mixing (driving nephelauxetic reduction) can occur through channels — like soft, polarizable σ-donation — that don't necessarily produce a large orbital energy splitting, which is why the two series need not track each other.

How does the nephelauxetic effect relate to the spin-only magnetic moment of a complex?

Indirectly, through its effect on spin-pairing energetics rather than through the spin-only formula itself. The spin-only magnetic moment depends only on the number of unpaired electrons, which is set by whether the complex is high-spin or low-spin — a competition between Δo and the mean pairing energy P, where P itself is built from Racah B and C (roughly P≈(combination of B,C)P \approx (\text{combination of } B,C) depending on the configuration). Because the nephelauxetic effect reduces B and C upon complexation, it lowers the pairing energy P needed to force electrons into the same orbital, which shifts the high-spin/low-spin crossover point compared to a naive calculation using free-ion B — a real, quantitatively important correction when predicting whether a borderline complex (Δo comparable to P) will be high-spin or low-spin.

Is there a single "true" free-ion B value, or does it vary with oxidation state?

Free-ion B values genuinely vary with oxidation state (and with which specific dⁿ configuration is being probed), because higher oxidation states contract the d orbitals (stronger effective nuclear charge, less electron-electron screening), which increases the Slater integrals F₂, F₄ and hence increases free-ion B. This is why tabulated free-ion Racah parameters are always listed per specific ion (Cr²⁺ vs. Cr³⁺ have different B_free, not the same "chromium" value), taken from gas-phase atomic emission spectroscopy of that exact charge state. Failing to use the correct oxidation-state-specific free-ion B is one of the most common numerical errors when computing β for a new complex.

Before this