Mössbauer Data Interpretation

Introduction

A raw Mössbauer spectrum is simply transmitted gamma counts against source velocity. Turning it into chemistry means separating three hyperfine interactions that can be present at once: the electric monopole term (isomer shift), the electric quadrupole term and the magnetic dipole term. When several iron sites overlap, the task becomes a careful piece of detective work. This page sets out a systematic approach to reading line positions and judging whether a fit makes sense.

Core explanation

Fitting the lines. Each absorption line is well described by a Lorentzian function with a position, a width and a depth. Software fits a model of one or more components to the data by least squares. The chemist chooses the model; the computer only optimises it. A good fit has random residuals, sensible widths (usually 0.2–0.4 mm s⁻¹) and physically reasonable parameters.

Doublets. For a quadrupole doublet with lines at v₁ and v₂: - isomer shift δ = (v₁ + v₂)/2 - quadrupole splitting ΔE Q = v₂ − v₁

Sextets. Number the lines 1 to 6 from negative to positive velocity. A pure magnetic sextet is symmetric. When an electric-field gradient is also present and is weak compared with the magnetic interaction, lines 1 and 6 shift one way and lines 2 to 5 shift the other way by the same amount. This gives three useful formulae: - δ = average of all six line positions (or of the centroids of pairs 1,6 and 2,5) - quadrupole shift 2ε = [(v₆ − v₅) − (v₂ − v₁)]/2 - hyperfine field B hf is proportional to v₆ − v₁; with α-iron as the standard, B hf = 33.0 T × (v₆ − v₁)/10.62 mm s⁻¹

The quadrupole shift in a sextet is not the same as ΔE Q in a doublet; it depends on the angle between the EFG axis and the hyperfine field.

Overlap and ambiguity. Two overlapping doublets can be fitted in different ways: nested (one inside the other) or crossed. For example, lines at positions a < b < c < d could be (a,d)+(b,c) or (a,c)+(b,d). Choose the pairing that gives chemically reasonable δ values, and test it with spectra at another temperature, where ΔE Q of high-spin Fe(II) usually changes while that of Fe(III) barely moves.

Temperature effects. Lowering temperature raises δ slightly because of the second-order Doppler shift (typically by about 0.1 mm s⁻¹ from 300 K to 4 K). It can also reveal magnetic ordering or slow relaxation, turning doublets into sextets.

Areas. The fractional area of each component estimates the fraction of iron in that environment, provided the recoil-free fractions are similar and the absorber is thin.

Step-by-step reasoning

A reliable interpretation workflow:

1. Check the velocity calibration and the reference (α-Fe at room temperature). 2. Classify the pattern: singlet, doublet, sextet or a mixture. 3. Fit the simplest chemically plausible model first. 4. Extract δ, ΔE Q or 2ε, and B hf for each component. 5. Compare with known ranges and test alternatives with temperature or applied-field spectra.

Visual explanation

Sketch a sextet and mark the six line positions with vertical ticks. Draw a bracket across lines 1 and 6 labelled "span ∝ B hf". Mark the midpoint of lines 1 and 6 and the midpoint of lines 2 and 5: their average is δ, and their difference reveals the quadrupole shift.

Real-world analogy

Interpreting overlapped spectra is like working out who is singing in a recording of several voices. Each voice (component) has its own pitch pattern. You propose a set of singers, subtract their parts and check whether anything unexplained remains in the recording.

Real-world example

Corrosion products on steel often contain goethite, lepidocrocite and magnetite together. At room temperature some of these give doublets and others sextets, and small particles blur the pattern. Recording spectra down to liquid-helium temperature, where each phase orders magnetically with a characteristic field, lets conservators quantify each phase.

Why?

Why is the separation of lines 1 and 6 used for the field rather than lines 3 and 4? The outer lines are the most widely separated and most intense, so their positions are measured most precisely, and a first-order quadrupole shift moves both outer lines equally, leaving their separation unchanged.

Common misconception

"More components always give a better answer." Adding components always lowers the residual, but may fit noise. Every component must be justified by chemistry, by known structure or by changes with temperature.

Worked example

Question: A sextet has lines at −8.00, −4.33, … , +5.27 and +8.54 mm s⁻¹ (lines 1, 2, 5 and 6 given). Find δ, 2ε and B hf.

Reasoning: Centroid of lines 1 and 6 = (−8.00 + 8.54)/2 = 0.27; centroid of lines 2 and 5 = (−4.33 + 5.27)/2 = 0.47; δ = (0.27 + 0.47)/2 = 0.37 mm s⁻¹. 2ε = [(8.54 − 5.27) − (−4.33 + 8.00)]/2 = (3.27 − 3.67)/2 = −0.20 mm s⁻¹. B hf = 33.0 × (16.54/10.62) ≈ 51.4 T.

Answer: δ = 0.37 mm s⁻¹, 2ε = −0.20 mm s⁻¹, B hf ≈ 51 T, matching haematite at room temperature.

Quick check

1. A doublet has lines at +0.05 and +1.95 mm s⁻¹. What are δ and ΔE Q? Answer: δ = 1.00 mm s⁻¹ and ΔE Q = 1.90 mm s⁻¹, consistent with high-spin iron(II).

Exam focus

Practise extracting δ from centroids, ΔE Q from doublet separations, and B hf from the outer-line span scaled to α-iron. Explain why overlapping doublets can be paired in more than one way and how temperature helps choose.

Advanced insight

When the quadrupole and magnetic interactions are comparable, the simple first-order formulae fail and the full nuclear Hamiltonian must be diagonalised; line positions and intensities then depend on η and on the orientation of the field relative to the EFG. Paramagnetic samples at low temperature may show relaxation spectra that are neither doublets nor sextets and require dynamic lineshape models.

Summary

Mössbauer interpretation separates three contributions: δ from the centroid, ΔE Q from doublet splitting or 2ε from sextet asymmetry, and B hf from the outer-line span. Fitting uses Lorentzian components chosen for chemical sense. Overlaps can be ambiguous, so temperature dependence, applied fields and relative areas are used to confirm assignments.

Practice questions

1. Give the formulae for δ and ΔE Q in terms of doublet line positions v₁ and v₂. Answer: δ = (v₁ + v₂)/2 and ΔE Q = v₂ − v₁. 2. A sextet has v₁ = −5.33 and v₆ = +5.29 mm s⁻¹. Estimate B hf and identify a likely material. Answer: Span = 10.62 mm s⁻¹, so B hf ≈ 33.0 T; this is α-iron metal. 3. Four lines appear at 0.10, 0.30, 0.60 and 2.40 mm s⁻¹. Suggest a chemically reasonable nested pairing into two doublets. Answer: (0.10, 2.40) gives δ = 1.25, ΔE Q = 2.30 mm s⁻¹ (high-spin Fe(II)); (0.30, 0.60) gives δ = 0.45, ΔE Q = 0.30 mm s⁻¹ (high-spin Fe(III)). A crossed pairing is also possible, so a second temperature should be measured to confirm. 4. Why does δ increase slightly when a sample is cooled from 300 K to 4 K? Answer: The second-order Doppler shift, caused by the mean-square velocity of the vibrating nucleus, decreases as thermal motion is reduced, so the measured centre moves to slightly more positive velocity.