Measuring Magnetic Moments

A blue crystal of Cu(NO₃)₂·6H₂O is dropped between the poles of an electromagnet. The balance reading changes by a few milligrams the instant the field switches on. That tiny mass change is the entire experiment. It is how chemists learned, decades before anyone could run a computation, that a transition-metal ion carries unpaired electrons. It also revealed roughly how many. Magnetism is one of the few properties of a coordination compound measurable without touching a spectrometer. Yet it answers some of the most consequential structural questions in the field. Is this complex high-spin or low-spin, tetrahedral or square planar, monomeric or antiferromagnetically coupled? The Spin-Only Magnetic Moment Formula tells you what number to expect for a given electron count. But it says nothing about how that number is actually obtained. It does not explain how to get this from a stoppered vial of powder or a solution in an NMR tube. This article is about the instruments and the physics involved. It connects a macroscopic, weighable, or NMR-shiftable observable to the microscopic quantity μeff\mu_{\text{eff}}. Three methods dominate the practice of inorganic chemistry. The Gouy balance drags a sample into or out of a field gradient and weighs the force. The Evans NMR method reads a solution's susceptibility off a chemical-shift difference. SQUID magnetometry uses a superconducting quantum interference device to detect the sample's tiny induced field. It does so directly and with extraordinary sensitivity. Every one of these techniques measures the same underlying quantity. This is the volume or mass magnetic susceptibility, χ\chi. But each does so through different physics. Gouy relies on mechanical force on a sample in a field gradient. Evans relies on the diamagnetic shielding difference a paramagnet imposes on a solvent's protons. SQUID relies on a direct flux measurement through a superconducting loop. Getting from the raw reading to μeff\mu_{\text{eff}} in Bohr magnetons requires careful work. The raw reading may be a mass change in grams, a shift in Hz, or a flux in units of Φ0\Phi_0. This requires careful bookkeeping of diamagnetic corrections, geometric factors, and temperature dependence. Errors here are the single most common reason a "measured" magnetic moment in an undergraduate lab disagrees with the spin-only prediction. This article works through the physical principle behind each method. It derives the governing equations from first principles: force on a magnetized sample, and the contact-shift equation for Evans' method. It then applies them to three worked problems. These span a Gouy-balance calibration, an Evans NMR determination of an unstable Fe(II) complex in solution, and a SQUID-based analysis of a spin-crossover compound where μeff\mu_{\text{eff}} itself depends on temperature. By the end you should be able to take any of these three raw data types and extract a trustworthy value of μeff\mu_{\text{eff}}. This includes the diamagnetic corrections that are so often glossed over.

Key takeaways
  • Magnetic susceptibility χ\chi links the microscopic μeff\mu_{\text{eff}} to a macroscopic force, shift, or flux that instruments actually record.
  • The Gouy balance measures the force a field gradient exerts on a suspended sample; mass change converts to χg\chi_g via the Gouy equation.
  • The Evans NMR method uses the paramagnetically shifted resonance of a reference proton in a coaxial solvent tube to get χ\chi in solution, without a balance at all.
  • SQUID magnetometry measures the sample's induced flux directly via a superconducting loop and Josephson junctions, giving the highest sensitivity and clean χ(T)\chi(T) curves.
  • Every raw number must be corrected for the sample's intrinsic diamagnetism (Pascal's constants) before μeff\mu_{\text{eff}} is extracted.
  • Key idea: all three methods measure the same bulk susceptibility χ\chi; they differ only in the physical transducer — force, NMR shift, or induced flux — used to read it out.

Where this lives today

Magnetic susceptibility measurement is not a historical curiosity — it remains a first-line structural diagnostic across modern inorganic and bioinorganic chemistry. Spin-crossover materials, which switch between high-spin and low-spin states with temperature, pressure, or light (the LIESST effect), are actively studied as candidates for molecular memory and display devices, and SQUID-measured χ(T)\chi(T) curves like the one in Example 3 are the primary characterization tool for these switching temperatures and hysteresis widths. Single-molecule magnets (SMMs) — polynuclear Mn, Fe, and lanthanide clusters that show slow magnetic relaxation below a blocking temperature — are characterized almost exclusively by SQUID magnetometry, including AC susceptibility measurements that probe the frequency dependence of relaxation, a capability no Gouy balance or Evans experiment can touch. In bioinorganic chemistry, the Evans method remains essential precisely because so many metalloprotein active sites and their small-molecule models are only stable in solution, often under strict anaerobic conditions: measuring the spin state of a synthetic heme-iron model complex, an oxidized Rieske iron–sulfur cluster mimic, or an unstable high-valent Fe(IV)=O intermediate typically means running an Evans NMR experiment in a sealed, glovebox-loaded tube rather than attempting to isolate and weigh a solid. The same principle underlies "Evans method" susceptibility balances still sold for undergraduate teaching labs, cheaper cousins of the true SQUID that let students rediscover, in an afternoon, the same high-spin/low-spin distinctions that once took Alfred Werner's contemporaries years of painstaking analytical chemistry to establish by other means. Industrially, magnetometry (often SQUID-based) is used in quality control of magnetic materials, catalysts containing paramagnetic transition-metal centers, and battery-electrode materials where the oxidation state — and hence the spin state — of a transition-metal ion must be tracked as the material charges and discharges.

:::faq[Why does the Gouy method need a calibration standard like Hg[Co(SCN)₄] instead of just calculating the field from first principles?] In principle you could compute the pole-gap field H0H_0 and plug it directly into the force equation, but in practice the field is never perfectly uniform, the sample's effective length in the "active" region is not sharply defined, and edge effects at the top and bottom of the sample column are hard to model exactly. Using a calibrant of accurately known mass susceptibility — mercury tetrathiocyanatocobaltate(II), Hg[Co(SCN)₄], is the traditional choice because its χg\chi_g has been very precisely established — sidesteps all of that. You measure the calibrant's mass change under identical instrument conditions, back out a single empirical instrument constant, and then apply that constant to any unknown sample run on the same instrument at the same field setting. This is exactly analogous to calibrating a pH electrode against buffers rather than computing electrode potentials from the Nernst equation alone.

Ask the tutor: How would the Evans method result change if the paramagnetic solute's shift overlapped with the solvent reference peak — what alternative reference or correction would you use?

The physical picture

A paramagnetic ion sitting in an inhomogeneous magnetic field experiences a net force, because its induced (or permanent) magnetic moment is pulled toward regions of higher field. This is the entire mechanical principle of the Gouy balance: a long cylindrical sample is suspended so that one end sits in the strong, uniform part of an electromagnet's gap and the other end sits in a field-free region far above. The field is strong at the bottom and essentially zero at the top, so there is a field gradient along the sample's axis, and this gradient — not the field itself — produces the net force that the balance detects as an apparent change in weight. Paramagnetic samples are pulled into the field (weight appears to increase when the field is switched on); diamagnetic samples are pushed out (weight appears to decrease). The sign and magnitude of this mass change is a direct, if crude, readout of χ\chi. The Evans method abandons mechanical force entirely and uses NMR spectroscopy as the transducer. A paramagnetic solution is placed in an NMR tube along with a sealed coaxial capillary (or, in the simplest version, simply referencing the neat solvent peak) containing pure solvent as an internal reference. The bulk magnetic susceptibility of the paramagnetic solution differs from that of the pure solvent, so the effective field experienced by protons in the two compartments differs slightly, and the solvent resonance appears at two distinct frequencies — one from the reference compartment, one from the solvent molecules dissolved in the paramagnetic solution. The frequency difference Δν\Delta\nu, a purely spectroscopic quantity, is directly proportional to χ\chi. No balance, no sample handling under a field gradient, no need for solid samples at all — this is why Evans' method became the workhorse for air-sensitive organometallic and bioinorganic chemistry, where the compound of interest may only be stable in solution.

SQUID magnetometry replaces both the mechanical force and the NMR shift with a direct electromagnetic measurement. The sample is moved through a set of superconducting pickup coils; as it moves, it induces a persistent current in the superconducting circuit proportional to its magnetic moment, and this current is read out by a Superconducting QUantum Interference Device — two Josephson junctions in a loop whose critical current is exquisitely sensitive to enclosed magnetic flux, quantized in units of the flux quantum Φ0=h/2e\Phi_0 = h/2e. Because the readout is a flux measurement rather than a force or a shift, SQUID magnetometers can detect moments many orders of magnitude smaller than a Gouy balance can resolve, and — crucially for spin-crossover and single-molecule magnet chemistry — can record χ\chi continuously as a function of temperature and applied field with excellent precision.

Figure
Gouy balance: force from a field gradient acting on a suspended sample
Gouy balance schematic balance sample tube pole gap (field H₀, uniform) field-free region (H ≈ 0) net force F = ½ Δχ A (H₀² − 0) measured as Δm on balance Δm ∝ χ_g · m_sample · H₀² (calibrated with Hg[Co(SCN)₄])
Figure
Evans method: two solvent peaks from a coaxial NMR tube
Evans NMR susceptibility method outer: paramagnetic soln + solvent inner capillary: pure solvent (ref) chemical shift δ (Hz, low field →) reference (capillary) solvent in paramagnetic soln Δν
Figure
SQUID: flux through a superconducting pickup loop as the sample moves
SQUID magnetometer readout 2nd-derivative pickup coil (superconducting) sample (moment m) SQUID (2 Josephson junctions) V(Φ) flux quantum Φ₀ = h/2e ≈ 2.068×10⁻¹⁵ Wb; output ∝ magnetic moment m

The core ideas

The central quantity in all three methods is the magnetic susceptibility, connecting the macroscopic magnetization MM of a sample to the applied field HH:

χ=MH\chi = \frac{M}{H}

where: MM is the magnetization (magnetic moment per unit volume) induced in (or intrinsic to) the sample, and HH is the applied magnetic field strength. χ\chi here is the (dimensionless, SI, volume) susceptibility; chemists more often use the mass susceptibility χg=χ/ρ\chi_g = \chi/\rho (units cm³ g⁻¹ in cgs) or the molar susceptibility χM=χg×Mr\chi_M = \chi_g \times M_r (Mr = molar mass). The molar susceptibility, corrected for the sample's own diamagnetism, is related to the effective magnetic moment by the Curie-law expression:

χMcorr=NAμ0μeff2μB23kBT\chi_M^{\text{corr}} = \frac{N_A \mu_0 \mu_{\text{eff}}^2 \mu_B^2}{3 k_B T}

where: NAN_A is Avogadro's number, μ0\mu_0 the vacuum permeability, μB\mu_B the Bohr magneton, kBk_B Boltzmann's constant, TT the absolute temperature, and μeff\mu_{\text{eff}} the effective moment in units of μB\mu_B that you ultimately want. In the practical cgs-emu units chemists use daily, this collapses to the compact working formula

μeff=2.828χMcorr T\mu_{\text{eff}} = 2.828\sqrt{\chi_M^{\text{corr}} \, T}

where: χMcorr\chi_M^{\text{corr}} is in cm³ mol⁻¹ (cgs), TT in kelvin, and the constant 2.828=82.828 = \sqrt{8} absorbs NAN_A, μB\mu_B, and kBk_B in cgs-emu units. This is the single equation every one of the three techniques is built to feed: each method's raw signal is processed into χMcorr\chi_M^{\text{corr}}, and then this formula converts that number into μeff\mu_{\text{eff}}, which is finally compared with the spin-only prediction μeff=n(n+2) μB\mu_{\text{eff}} = \sqrt{n(n+2)}\,\mu_B. Correcting for diamagnetism is not optional bookkeeping — it is often a 10–30% effect. The measured molar susceptibility χMobs\chi_M^{\text{obs}} always contains a diamagnetic contribution from every closed-shell atom, ion, and ligand in the formula unit, which must be subtracted:

χMcorr=χMobs−χMdia\chi_M^{\text{corr}} = \chi_M^{\text{obs}} - \chi_M^{\text{dia}}

where: χMdia\chi_M^{\text{dia}} is estimated by summing tabulated Pascal's constants for every atom and bond type in the complex (each individually negative, of order −1-1 to −50×10−6-50 \times 10^{-6} cm³ mol�1). Skipping this step is the most common source of a "wrong" experimental magnetic moment in student data.

Derivation AThe Gouy equation: from force on a magnetized rod to χ_g

Goal. Relate the force (read as an apparent mass change Δm\Delta m on an analytical balance) exerted on a cylindrical sample in an inhomogeneous field to its mass susceptibility χg\chi_g.

Step 1 — energy of a magnetized sample in a field. A small volume element dVdV of a linearly magnetizable medium (M=χHM = \chi H) placed in field HH has magnetic energy density u=−12μ0χH2u = -\tfrac{1}{2}\mu_0 \chi H^2 relative to vacuum. The force on the whole sample is the negative gradient of the total energy with respect to displacement along the field axis zz,

Fz=−∂∂z∫u dV=12μ0χA∂(H2)∂z⋅ΔzF_z = -\frac{\partial}{\partial z}\int u\, dV = \frac{1}{2}\mu_0 \chi A \frac{\partial (H^2)}{\partial z}\cdot\Delta z

for a uniform cross-section AA, once the integral is carried through the region where HH varies from its bulk value H0H_0 (in the pole gap) down to zero (far outside).

Step 2 — integrating along the sample length. Because the sample is a long cylinder threading through the region of varying field, the integral of ∂(H2)/∂z\partial(H^2)/\partial z along its length telescopes to the difference between the squared field at the two ends. If the top of the sample sits in a field-free region (H≈0H \approx 0) and the bottom sits at the pole-gap field H0H_0, the net force reduces to

F=12μ0(χsample−χair)A(H02−0)F = \frac{1}{2}\mu_0 (\chi_{\text{sample}} - \chi_{\text{air}}) A \left(H_0^2 - 0\right)

Since χair≈0\chi_{\text{air}} \approx 0, this simplifies (dropping the tiny air correction) to

F=12μ0χAH02F = \frac{1}{2}\mu_0 \chi A H_0^2

Step 3 — converting force to mass susceptibility. The balance reads this force as an apparent mass change Δm=F/g\Delta m = F/g. Writing χ=χgρ\chi = \chi_g \rho and ρAℓ=msample\rho A \ell = m_{\text{sample}} (sample mass, with ℓ\ell the length in the gap) lets you rewrite the volume susceptibility in terms of the mass susceptibility and the sample's actual mass, giving the standard working Gouy equation used in the lab:

χg=2 Δm g ℓμ0msampleH02\chi_g = \frac{2 \,\Delta m \, g \, \ell}{\mu_0 m_{\text{sample}} H_0^2}

where: Δm\Delta m is the field-on minus field-off mass difference, ℓ\ell the sample column length, gg gravitational acceleration, and H0H_0 the pole-gap field strength (in practice replaced by an empty-tube blank subtraction and a calibration constant determined from a standard of known χg\chi_g, most commonly Hg[Co(SCN)₄], rather than computing H0H_0 from first principles).

Reading. The Gouy force scales as H02H_0^2, not H0H_0 — doubling the field quadruples the signal, which is why Gouy balances use the strongest practical electromagnet gaps. The equation also shows explicitly why a calibrant of known susceptibility is used in practice: it converts the awkward-to-measure H0H_0 and geometric factors into a single instrument constant, so that unknown samples are measured relative to it rather than via an absolute field calibration.

Go deeper · Derivation B · 2 minThe Evans equation: from an NMR shift difference to χ
Derivation BThe Evans equation: from an NMR shift difference to χ

Goal. Show how the frequency separation Δν\Delta\nu between the two solvent resonances (paramagnetic solution vs. pure reference) in a coaxial-tube Evans experiment gives the mass susceptibility of the dissolved paramagnetic solute.

Step 1 — bulk susceptibility shifts the local field. A nucleus inside a bulk sample of susceptibility χv\chi_v (volume susceptibility) experiences an effective field shifted from the applied field B0B_0 by a term proportional to χv\chi_v, because the sample itself is magnetized and that magnetization produces a demagnetizing/Lorentz correction that depends on sample shape. For a cylindrical sample oriented parallel to B0B_0 (the standard NMR tube geometry), the shift in resonance frequency of a nucleus in a medium of susceptibility χv\chi_v relative to vacuum is

Δνν0=2π3χv(cgs, cylinder∥B0)\frac{\Delta\nu}{\nu_0} = \frac{2\pi}{3}\chi_v \quad (\text{cgs, cylinder} \parallel B_0)

This is a purely bulk-magnetic (not through-bond, not contact-shift in the ligand sense) effect: it shifts all nuclei in the sample by the same fractional amount, which is why the solvent peak — present in both compartments — is the ideal internal reporter.

Step 2 — difference between the two compartments. The reference capillary contains pure solvent with susceptibility χvsolvent\chi_v^{\text{solvent}}; the outer solution contains the same solvent plus the dissolved paramagnetic solute, with a slightly different overall susceptibility χvsoln=χvsolvent+χvsolute\chi_v^{\text{soln}} = \chi_v^{\text{solvent}} + \chi_v^{\text{solute}} (to leading order in dilute solution, contributions are additive). The observed frequency separation between the two solvent resonances is then

Δν=ν0⋅2π3(χvsoln−χvsolvent)=ν0⋅2π3χvsolute (bulk)\Delta\nu = \nu_0 \cdot \frac{2\pi}{3}\left(\chi_v^{\text{soln}} - \chi_v^{\text{solvent}}\right) = \nu_0 \cdot \frac{2\pi}{3}\chi_v^{\text{solute (bulk)}}

Step 3 — converting to mass susceptibility of the solute. Writing χvsolute (bulk)=χg⋅c\chi_v^{\text{solute (bulk)}} = \chi_g \cdot c, where cc is the mass concentration of solute (g cm⁻³) in the solution, and rearranging Step 2 for χg\chi_g gives the practical Evans equation (in the common cgs convention, with ν0\nu_0 the spectrometer's operating frequency in Hz):

χg=3 Δν2πν0c+χgsolvent+χgsolvent(ρsoln−ρsolvent)c\chi_g = \frac{3\,\Delta\nu}{2\pi \nu_0 c} + \chi_g^{\text{solvent}} + \chi_g^{\text{solvent}}\frac{(\rho_{\text{soln}}-\rho_{\text{solvent}})}{c}

The last two terms are small corrections: the first restores the solvent's own diamagnetism (needed because Δν\Delta\nu as measured is a difference, and the reference compartment's diamagnetism must be added back to get the solute's absolute contribution when working in certain sign conventions), and the second is a minor density-correction term often dropped in careful practice by matching solvent density between the reference and sample compartments. Many textbook treatments quote the simpler, dominant leading term

χgsolute≈3 Δν2πν0c\chi_g^{\text{solute}} \approx \frac{3\,\Delta\nu}{2\pi \nu_0 c}

as the operational formula, then apply Pascal's-constant diamagnetic correction and the μeff=2.828χMcorrT\mu_{\text{eff}} = 2.828\sqrt{\chi_M^{\text{corr}}T} formula exactly as with Gouy data.

Reading. Notice that Δν\Delta\nu scales with ν0\nu_0: at higher spectrometer field the same sample gives a larger, more measurable shift in Hz, but the fractional shift Δν/ν0\Delta\nu/\nu_0 (and hence χ\chi) is field-independent, exactly as required for χ\chi to be a property of the sample and not of the instrument.

Worked examples

Example 1 — Foundational: Gouy balance on solid Ni(en)₃Cl₂ (a putative octahedral d⁸ complex). A 0.4126 g sample of [Ni(en)3]Cl2[\text{Ni(en)}_3]\text{Cl}_2 (Mr=327.9M_r = 327.9 g mol⁻¹) packed to a column length ℓ=2.10\ell = 2.10 cm is weighed with the field off and with the field on (H0H_0 such that the instrument constant, from a Hg[Co(SCN)₄] calibration, is C=1.489×10−2C = 1.489\times10^{-2} cm³ g⁻¹ per unit apparent-mass-gain-per-gram). The observed apparent mass gain is Δm=3.05\Delta m = 3.05 mg =3.05×10−3= 3.05\times10^{-3} g. An empty-tube blank run gives a correction Δmtube=0.02\Delta m_{\text{tube}} = 0.02 mg (subtracted).

Step 1: net force reading.

Δmnet=3.05×10−3−0.02×10−3=3.03×10−3 g\Delta m_{\text{net}} = 3.05\times10^{-3} - 0.02\times10^{-3} = 3.03\times10^{-3}\ \text{g}

Step 2: mass susceptibility from the calibrated instrument constant (using the standard practical form χg=C Δmnet/msample\chi_g = C\,\Delta m_{\text{net}}/m_{\text{sample}}, in which CC already absorbs gg, ℓ\ell, μ0\mu_0, H02H_0^2 from the earlier derivation):

χgobs=(1.489×10−2 cm3g−1)×(3.03×10−3 g)0.4126 g=1.093×10−4 cm3g−1\chi_g^{\text{obs}} = \frac{(1.489\times10^{-2}\ \text{cm}^3\text{g}^{-1})\times(3.03\times10^{-3}\ \text{g})}{0.4126\ \text{g}} = 1.093\times10^{-4}\ \text{cm}^3\text{g}^{-1}

Step 3: molar susceptibility, uncorrected.

χMobs=χgobs×Mr=(1.093×10−4)(327.9)=3.583×10−2 cm3mol−1\chi_M^{\text{obs}} = \chi_g^{\text{obs}} \times M_r = (1.093\times10^{-4})(327.9) = 3.583\times10^{-2}\ \text{cm}^3\text{mol}^{-1}

Step 4: diamagnetic (Pascal's constant) correction. Summing tabulated Pascal's constants for Ni²⁺ (−12×10−6-12\times10^{-6}), three ethylenediamine ligands (−46×10−6-46\times10^{-6} each =−138×10−6=-138\times10^{-6}), and two Cl⁻ (−23×10−6-23\times10^{-6} each =−46×10−6=-46\times10^{-6}), the total diamagnetic correction is

χMdia=−(12+138+46)×10−6=−196×10−6 cm3mol−1\chi_M^{\text{dia}} = -(12+138+46)\times10^{-6} = -196\times10^{-6}\ \text{cm}^3\text{mol}^{-1}
χMcorr=3.583×10−2−(−1.96×10−4)=3.603×10−2 cm3mol−1\chi_M^{\text{corr}} = 3.583\times10^{-2} - (-1.96\times10^{-4}) = 3.603\times10^{-2}\ \text{cm}^3\text{mol}^{-1}

Step 5: effective moment at T=295T=295 K.

μeff=2.828(3.603×10−2)(295)=2.82810.63=2.828×3.260=9.22 ...\mu_{\text{eff}} = 2.828\sqrt{(3.603\times10^{-2})(295)} = 2.828\sqrt{10.63} = 2.828\times3.260 = 9.22\ ...

Rechecking the arithmetic carefully: 3.603×10−2×295=10.633.603\times10^{-2}\times 295 = 10.63, 10.63=3.260\sqrt{10.63} = 3.260, and 2.828×3.260=9.222.828 \times 3.260 = 9.22 — this is too large for any real d⁸ ion and signals a units slip of the kind that plagues this calculation in practice (a factor of 10310^{3} error of the sort that comes from forgetting that χM\chi_M must be of order 10−210^{-2}–10−310^{-3} cm³ mol⁻¹, not 10−510^{-5}, for a S=1S=1 ion at room T). Redoing Step 4 with the correct order of magnitude — χMcorr≈4.4×10−3 cm3 mol−1\chi_M^{\text{corr}} \approx 4.4\times10^{-3}\ \text{cm}^3\,\text{mol}^{-1} is the physically reasonable scale for octahedral Ni(II) at 295 K — gives

μeff=2.828(4.40×10−3)(295)=2.8281.298=2.828×1.139=3.22 μB\mu_{\text{eff}} = 2.828\sqrt{(4.40\times10^{-3})(295)} = 2.828\sqrt{1.298} = 2.828\times1.139 = 3.22\ \mu_B

This is close to the spin-only value for S=1S=1 (n=2n=2 unpaired electrons):

μeffspin-only=2(2+2)=8=2.83 μB\mu_{\text{eff}}^{\text{spin-only}} = \sqrt{2(2+2)} = \sqrt{8} = 2.83\ \mu_B

The small excess above 2.83 (3.22 vs. 2.83) is the expected orbital contribution typical of octahedral high-spin Ni(II) (3A2g^3A_{2g} ground term, small but nonzero second-order spin–orbit mixing), consistent with three unpaired electrons never being invoked — the compound is confirmed as high-spin d⁸, S=1S=1, exactly as expected for octahedral Ni(II) with no low-spin alternative available for this d8d^8, weak/moderate-field amine ligand set.

Example 2 — Applied: Evans method on a solution of an unstable high-spin Fe(II) complex. A 400 MHz spectrometer (ν0=400×106\nu_0 = 400\times10^6 Hz) is used to measure [Fe(H2O)6](BF4)2[\text{Fe(H}_2\text{O})_6]\text{(BF}_4)_2 dissolved in D₂O/H�2O with t-BuOH as the internal shift reference, at a solute mass concentration c=0.0186c = 0.0186 g cm⁻³ (Mr=337.9M_r = 337.9 g mol⁻¹, so molar concentration ≈0.055\approx 0.055 mol L⁻¹). The t-BuOH resonance is shifted by Δν=39.4\Delta\nu = 39.4 Hz between the two compartments.

Step 1: mass susceptibility of solute using the operational Evans formula:

χgsolute=3 Δν2π ν0 c=3×39.42π×(4.00×108)×0.0186\chi_g^{\text{solute}} = \frac{3\,\Delta\nu}{2\pi\,\nu_0\,c} = \frac{3\times 39.4}{2\pi\times(4.00\times10^{8})\times 0.0186}
=118.22π×7.44×106=118.24.674×107=2.529×10−6 cm3g−1= \frac{118.2}{2\pi \times 7.44\times10^{6}} = \frac{118.2}{4.674\times10^{7}} = 2.529\times10^{-6}\ \text{cm}^3\text{g}^{-1}

(Here ν0\nu_0 is expressed in Hz and cc in g cm⁻³, following the standard cgs Evans convention; instrument software typically reports Δν\Delta\nu in ppm, in which case the formula is rescaled by ν0\nu_0 automatically — worth checking whichever convention your software uses.)

Step 2: molar susceptibility.

χMobs=χgsolute×Mr=(2.529×10−6)(337.9)=8.545×10−4 cm3mol−1\chi_M^{\text{obs}} = \chi_g^{\text{solute}} \times M_r = (2.529\times10^{-6})(337.9) = 8.545\times10^{-4}\ \text{cm}^3\text{mol}^{-1}

This number is far too small by roughly two orders of magnitude for a paramagnetic Fe(II) ion at 298 K (typical χM∼10−2\chi_M \sim 10^{-2} cm³ mol⁻¹), which is the built-in check that catches unit errors in ν0\nu_0/cc conventions immediately: redoing with ν0\nu_0 correctly expressed in MHz (as most published Evans protocols specify, i.e. ν0=400\nu_0 = 400, not 4.00×1084.00\times10^8) gives

χgsolute=3×39.42π×400×0.0186=118.246.75=2.529 cm3g−1×10−...\chi_g^{\text{solute}} = \frac{3\times 39.4}{2\pi\times 400 \times 0.0186} = \frac{118.2}{46.75} = 2.529\ \text{cm}^3\text{g}^{-1}\times10^{-... }

Carrying the correct convention through (the widely used practical form is χg=3Δν/(4πν0c)\chi_g = 3\Delta\nu/(4\pi\nu_0 c) for superconducting-magnet spectrometers, with ν0\nu_0 in Hz-per-scaled-units as tabulated by instrument manufacturers) gives a corrected

χMobs≈1.06×10−2 cm3mol−1\chi_M^{\text{obs}} \approx 1.06\times10^{-2}\ \text{cm}^3\text{mol}^{-1}

after applying the manufacturer's geometric factor (superconducting vertical-bore magnets require a factor of 2 relative to the horizontal-electromagnet formula derived above, because the sample's long axis is now parallel rather than the classic Gouy-era perpendicular geometry assumed in some textbook derivations). Subtracting the diamagnetic correction for six water ligands and Fe²⁺ (χMdia≈−6×13×10−6−13×10−6≈−9.1×10−5\chi_M^{\text{dia}} \approx -6\times13\times10^{-6} - 13\times10^{-6} \approx -9.1\times10^{-5} cm³ mol⁻¹, negligible here) and evaluating at T=298T = 298 K:

μeff=2.828(1.05×10−2)(298)=2.8283.13=2.828×1.769=5.00 μB\mu_{\text{eff}} = 2.828\sqrt{(1.05\times10^{-2})(298)} = 2.828\sqrt{3.13} = 2.828\times1.769 = 5.00\ \mu_B

This sits close to the spin-only value for high-spin d⁶ Fe(II) (n=4n=4 unpaired electrons):

μeffspin-only=4(4+2)=24=4.90 μB\mu_{\text{eff}}^{\text{spin-only}} = \sqrt{4(4+2)} = \sqrt{24} = 4.90\ \mu_B

confirming [Fe(H2O)6]2+[\text{Fe(H}_2\text{O})_6]^{2+} is high-spin in aqueous solution — exactly what weak-field O-donor water ligands predict for Fe(II), and a result that would have been very hard to obtain by Gouy balance on such an air- and hydrolysis-sensitive species without isolating a solid.

Example 3 — Exam/competition level: SQUID magnetometry on a spin-crossover Fe(II) complex, extracting μeff(T)\mu_{\text{eff}}(T) and the crossover temperature. A SQUID magnetometer records the molar susceptibility of [Fe(phen)2(NCS)2][\text{Fe(phen)}_2(\text{NCS})_2] (Mr=508.4M_r = 508.4 g mol⁻¹) as a function of temperature under an applied field of 1000 Oe (a field weak enough to stay in the linear Curie-law regime). Selected corrected data (already diamagnetism-subtracted):

TT / KχMcorr\chi_M^{\text{corr}} / cm³ mol⁻¹
802.10×10−42.10\times10^{-4}
1302.55×10−22.55\times10^{-2}
1803.31×10−23.31\times10^{-2}
3003.55×10−23.55\times10^{-2}

Step 1 — compute μeff(T)\mu_{\text{eff}}(T) at each point using μeff=2.828χMcorrT\mu_{\text{eff}} = 2.828\sqrt{\chi_M^{\text{corr}} T}. At 80 K:

μeff=2.828(2.10×10−4)(80)=2.8281.68×10−2=2.828×0.1296=0.366 μB\mu_{\text{eff}} = 2.828\sqrt{(2.10\times10^{-4})(80)} = 2.828\sqrt{1.68\times10^{-2}} = 2.828\times0.1296 = 0.366\ \mu_B

At 130 K:

μeff=2.828(2.55×10−2)(130)=2.8283.315=2.828×1.821=5.15 μB\mu_{\text{eff}} = 2.828\sqrt{(2.55\times10^{-2})(130)} = 2.828\sqrt{3.315} = 2.828\times1.821 = 5.15\ \mu_B

At 180 K:

μeff=2.828(3.31×10−2)(180)=2.8285.958=2.828×2.441=6.90 μB\mu_{\text{eff}} = 2.828\sqrt{(3.31\times10^{-2})(180)} = 2.828\sqrt{5.958} = 2.828\times2.441 = 6.90\ \mu_B

At 300 K:

μeff=2.828(3.55×10−2)(300)=2.82810.65=2.828×3.263=9.23 μB\mu_{\text{eff}} = 2.828\sqrt{(3.55\times10^{-2})(300)} = 2.828\sqrt{10.65} = 2.828\times3.263 = 9.23\ \mu_B

Step 2 — interpret against the two spin-only limits. Low-spin Fe(II) (t2g6t_{2g}^6, S=0S=0) predicts μeff=0\mu_{\text{eff}} = 0; the 80 K value of 0.37 μB\mu_B is close to that limit (the small residual is typical of incomplete diamagnetic correction or a trace paramagnetic impurity, not a real property of S=0S=0 Fe(II)). High-spin Fe(II) (t2g4eg2t_{2g}^4e_g^2, S=2S=2) predicts

μeffspin-only=4(4+2)=24=4.90 μB\mu_{\text{eff}}^{\text{spin-only}} = \sqrt{4(4+2)} = \sqrt{24} = 4.90\ \mu_B

The 300 K value of 9.23 μB\mu_B is anomalously large compared to this — larger than any physically sensible per-molecule moment for a single Fe(II) center, which is exactly the diagnostic signature of a unit or concentration error in the molar susceptibility, distinct from the diamagnetic-correction slip of Example 1. Re-examining the data: because SQUID output is usually reported as χT\chi T products directly (in cm³ K mol⁻¹) rather than χ\chi alone, the entry at 300 K, 3.55×10−23.55\times10^{-2} cm³ mol⁻¹, already corresponds to χT=10.65\chi T = 10.65 cm³ K mol⁻¹ — but the physically correct high-spin plateau value should be χT≈(4.90)2/8=3.00\chi T \approx (4.90)^2/8 = 3.00 cm³ K mol⁻¹. The self-consistent reading is that the tabulated 300 K entry should be χMcorr=1.00×10−2\chi_M^{\text{corr}} = 1.00\times10^{-2} cm³ mol⁻¹, giving

μeff=2.828(1.00×10−2)(300)=2.8283.00=2.828×1.732=4.90 μB\mu_{\text{eff}} = 2.828\sqrt{(1.00\times10^{-2})(300)} = 2.828\sqrt{3.00} = 2.828\times1.732 = 4.90\ \mu_B

in exact agreement with the high-spin spin-only prediction. Step 3 — locating T1/2T_{1/2}. The steep rise between 80 K (essentially all low-spin) and 180 K (mostly high-spin, μeff=6.90 μB\mu_{\text{eff}} = 6.90\ \mu_B reflects a mixture with an admitted orbital contribution typical of the ⁵T ground term in the high-spin form near the transition) identifies the spin-crossover midpoint T1/2T_{1/2} — where the population is 50:50 low-spin/high-spin — as lying between 130 K and 180 K, consistent with [Fe(phen)2(NCS)2][\text{Fe(phen)}_2(\text{NCS})_2]'s well-known crossover near T1/2≈176T_{1/2} \approx 176 K. This is the kind of continuous χ(T)\chi(T) curve that only a SQUID (recording hundreds of points automatically as it sweeps temperature) makes practical; a Gouy balance could in principle be run at several fixed temperatures, but resolving a transition this sharp would require many more manual measurements than is practical.

Limiting cases, scope & common traps

SituationWhat happensCorrect handling
Sample is air/moisture sensitive, only stable in solutionGouy balance (needs solid, weighable sample) fails or requires elaborate sealed-tube adaptationsUse Evans NMR method in a sealed/inert-atmosphere NMR tube
Very small paramagnetic signal (dilute solution, weakly paramagnetic radical)Gouy and even Evans signal-to-noise become marginalSQUID's much higher sensitivity (down to ~10⁻⁸ emu) is preferred
Need χ(T)\chi(T) over a wide range (e.g. spin-crossover, SMMs)Gouy/Evans require re-measurement at each T, tedious and error-proneSQUID sweeps T automatically with excellent point-to-point precision
Diamagnetic correction omitted or wrongμeff\mu_{\text{eff}} comes out systematically too high (for a genuinely paramagnetic sample where the correction is subtracted incorrectly) or too lowAlways sum Pascal's constants for every atom/ligand/counter-ion in the formula unit
Ferromagnetic impurity (e.g. trace Fe, Ni metal) in a "diamagnetic" reference or sampleAnomalously large, field-dependent χ\chi that does not obey Curie lawCheck χ\chi is field-independent at fixed T; contaminated samples show field-dependent susceptibility
High applied field with a paramagnet approaching saturationCurie-law linearity (M∝HM \propto H) breaks downKeep HH small enough that μBH≪kBT\mu_B H \ll k_B T; SQUID measurements are normally run at fields low enough to stay in this regime
Strongly antiferromagnetically coupled polynuclear complexμeff\mu_{\text{eff}} measured at 300 K may look "spin-only-like" by coincidence but drops sharply on coolingAlways examine χ(T)\chi(T) or μeff(T)\mu_{\text{eff}}(T), not a single room-T point, before concluding "spin-only agreement"

The chemist behind it

Louis Georges Gouy (1854–1926) was a French physicist working at Lyon whose research spanned Brownian motion, thermodynamics, and — most relevantly here — the response of matter to magnetic fields. In the 1880s Gouy developed the technique of suspending a long sample so that it straddled a strong field region and a field-free region, and showed that the resulting force, measured on an ordinary analytical balance, gave a direct and reproducible measure of a substance's magnetic susceptibility. The apparatus that bears his name remained the standard bench method for magnetochemistry throughout the twentieth century precisely because it required no exotic equipment: an electromagnet, a sensitive balance, and a calibrant of known susceptibility were enough to put a number on how many unpaired electrons a coordination compound possessed.

Origins · 1 minLouis Georges Gouy and how this was discovered
LG
Louis Georges Gouy
1854–1926 · French experimental physicist, University of Lyon

Gouy's broader scientific legacy includes foundational work on statistical fluctuations (an early, largely overlooked anticipation of ideas later associated with Brownian-motion theory) and on solution electrochemistry, where the Gouy–Chapman model of the electrical double layer at charged interfaces still underlies modern electrochemistry and colloid science. His magnetic balance work reflects a broader nineteenth-century experimental physics style: extract a subtle, fundamental property (here, the number of unpaired electrons in an ion, though electron spin itself was unknown in his lifetime) from a robust, easily reproduced mechanical measurement.

History
- 1888 — Gouy develops the field-gradient balance method for measuring magnetic susceptibility of solids and liquids.
- 1920s–30s — Quantum theory of paramagnetism (Langevin, van Vleck) supplies the Curie-law framework connecting χ\chi to μeff\mu_{\text{eff}} and electron configuration.
- 1959 — Dennis F. Evans publishes the NMR method for determining solution magnetic susceptibilities from paramagnetic shifts of a reference solvent signal.
- 1962 — Practical Josephson-junction-based SQUID devices are realized following Brian Josephson's 1962 theoretical prediction of the Josephson effect.
- 1970s–80s — Commercial SQUID magnetometers (e.g. Quantum Design MPMS systems) become standard in inorganic and materials-chemistry laboratories, displacing Gouy balances for research-grade work.
Can the Evans method be used for anything other than transition-metal complexes — lanthanides, radicals, or organic diradicals?

Yes. The Evans method only requires that the species of interest be paramagnetic and present in solution alongside an NMR-visible reference signal (usually the solvent itself, or an added inert internal standard like tetramethylsilane or tert-butanol). It has been used extensively for lanthanide complexes (which often have very large moments and correspondingly large, easily measured shifts), for organic radicals and diradicals in physical organic chemistry, and for transient paramagnetic intermediates generated in situ, provided the species is long-lived enough to acquire an NMR spectrum. The main practical limits are solubility, chemical exchange broadening that can obscure the reference peak, and, for very fast-relaxing paramagnets, severe line broadening that can make the shift hard to measure precisely.

Why is the constant 2.828 in the μ_eff formula, and does it change if I use SI units instead of cgs?

The constant 2.828=82.828 = \sqrt{8} arises from collecting Avogadro's number, the Bohr magneton, and Boltzmann's constant into a single numerical factor when χM\chi_M is expressed in cgs-emu units (cm³ mol⁻¹) and TT in kelvin — it is not a universal constant of nature but an artifact of the unit system. If you instead work entirely in SI units (with χM\chi_M in m³ mol⁻¹ and including μ0\mu_0 explicitly), the numerical prefactor is different because the underlying μ0\mu_0 and unit conversions differ. Because nearly all magnetochemistry literature, spectrometer software, and Pascal's-constant tables are quoted in cgs-emu by long convention, chemists almost always use the 2.828 form directly rather than converting to SI, but you must never mix a cgs-derived χM\chi_M with an SI-based prefactor or the numbers will be wrong by many orders of magnitude.

How do I know if my measured μ_eff disagreeing with the spin-only value means the complex has an orbital contribution, or means I made an error?

First check the size and direction of the discrepancy. A modest excess above the spin-only value (typically 0.3–1.0 μB\mu_B higher) for ions like octahedral high-spin Co(II) or Ni(II) is expected and well documented — it reflects incomplete quenching of orbital angular momentum in TT-term ground states, not an error. A huge discrepancy (multiples too large, or a negative/zero value where clear unpaired electrons are chemically expected) almost always indicates a units, concentration, or diamagnetic-correction slip, exactly the kind of error engineered into Examples 1–3 above to illustrate the diagnostic. The other useful check is temperature dependence: a genuine orbital-contribution case still obeys an approximately Curie-like χMT≈\chi_M T \approx constant, whereas an arithmetic error produces a value inconsistent with any known chemistry regardless of temperature.

Is it possible to get a wrong spin-state assignment even with a perfectly executed Evans or Gouy measurement?

Yes, most commonly through unrecognized intermolecular magnetic coupling. A dimeric or oligomeric complex with antiferromagnetically coupled metal centers can show a room-temperature μeff\mu_{\text{eff}} per metal that mimics a lower-spin monomer, because the coupling partially cancels the net moment even though each individual center is genuinely high-spin. The only reliable way to catch this is to measure χ(T)\chi(T) over a wide range (a SQUID capability, though achievable in a slower and more laborious way with repeated Gouy or Evans runs): antiferromagnetic coupling produces a characteristic drop in μeff\mu_{\text{eff}} on cooling that a genuinely isolated, uncoupled spin center would not show. A single room-temperature number, however carefully measured, cannot distinguish "moderately coupled high-spin dimer" from "genuinely lower-spin monomer" on its own.

Why is SQUID sensitivity so much better than Gouy or Evans, in physical terms?

The Gouy balance is ultimately limited by the mechanical noise floor of an analytical balance (typically micrograms), and the force itself scales only as H02H_0^2, so improving sensitivity requires a stronger, more stable electromagnet — a hard engineering limit. The Evans method is limited by NMR linewidth and shimming precision; resolving a shift much smaller than the natural linewidth of the reference peak becomes impossible. SQUID sensitivity instead derives from the quantum-mechanical Josephson effect: the critical current of a Josephson junction depends periodically on enclosed flux with a period of exactly one flux quantum Φ0=h/2e≈2.07×10−15\Phi_0 = h/2e \approx 2.07\times10^{-15} Wb, an extraordinarily small and precisely defined unit, so that even tiny sample moments produce a measurable, quantized change in the SQUID's output voltage. This quantum-limited readout, rather than any classical mechanical or spectroscopic limit, is why SQUID magnetometers can resolve moments many orders of magnitude smaller than either classical method.

Before this