NMR, ESR, Mössbauer (µSR)
(for solid- state physics magnetism)
- P. Mendels
- Lab. Physique des solides
- Univ. Paris-Sud Orsay
Frustrated magnetism Spin liquids Metals Correlated electrons Superconductivity Spin chains
NMR, ESR, Mssbauer (SR) (for solid- state physics magnetism) P. - - PowerPoint PPT Presentation
NMR, ESR, Mssbauer (SR) (for solid- state physics magnetism) P. Mendels Metals Lab. Physique des solides Correlated electrons Frustrated magnetism Univ. Paris-Sud Orsay Superconductivity Spin liquids Spin chains NMR, ESR,
(for solid- state physics magnetism)
Frustrated magnetism Spin liquids Metals Correlated electrons Superconductivity Spin chains
(for solid- state physics magnetism)
formalism
(ii) the time window, the field range (iii) sensitivity and pulsed versus continuum
Zeeman, Nobel Physics 1902 Rabi, Nobel Physics 1944
Nuclear spin Electronic spin
The probe Hamiltonian is a weak perturbation of the electronic system; acts like a spy. (µSR also)
More involved treatment
Sweep the frequency at a fixed external field
Sweep the field at a constant frequency
NMR, Mössbauer and ESR
Note: Highest similarity is between NMR and µSR, see D. Andreica Many thanks to:
. Bonville (Mössbauer, CEA Saclay);
(Some slides also borrowed from Carretta, Murad, Takigawa)
J
Bloch & Purcell, Nobel Physique 1952 Ernst, Nobel Chemistry 1991 Wuthrich, Nobel Chemistry 2002 Lauterbur & Mansfeld, Nobel Medecine 2003
7 Tesla 23 Tesla
1/2 3/2 5/2
Need for r.f. field h1(n) h1H0 H = - m.(H0 + h1 ) = - g ħ (H0 Iz + h1 Ix) nNMR = g/2p H0
Nuclear spin I in a magnetic field H0 Zeeman effect : H = - m.H0 = - g ħ H0 Iz Energy levels E = - m g ħ H0 , m=-I, -I+1 ... I-1, I
. . . . . .
1/2 3/2 5/2
Spatially resolved magnetometer Nuclear spin I in a magnetic field H0 r.f. field h1 Zeeman effect : H = - m.H0 = - g ħ H0 Iz Energy levels E = - m g ħ H0 , m=-I, -I+1 ... I-1, I
. . . . . .
nNMR = g/2p (H0 + Hlocal)
Screening of H0 by electrons modification of orbitals by the applied field
Indirect interaction between nuclear moments (electrons)
A very involved Hamiltonian…coupling to electronic moments A very involved Hamiltonian…coupling to electronic moments and surrounding charges
n e n e n e hf
2 5 3 2 3 2
Orbital effect Dipolar effect from An unpaired spin Contact contribution from an unpaired spin on a s orbital No!
Very strong - Isotropic weak - anisotropic weak - anisotropic
n e n e n e hf
2 5 3 2 3 2
, ,
polarizati core contact i dip i
z y x i
Orbital effect
Gyromagnetic ratio: depends on the nucleus
i
Spin-dipolar effect from an unpaired spin s Contact contribution from an unpaired spin on a s orbital
hf
dip hf
contact hf
Orbital shift Spin shift
(hyperfine interaction)
n e n e n e hf
2 5 3 2 3 2
, ,
polarizati core contact i dip i
z y x i
Orbital effect
Gyromagnetic ratio: depends on the nucleus
i
Spin-dipolar effect from an unpaired spin s Contact contribution from an unpaired spin on a s orbital
hf
dip hf
contact hf
Orbital shift Spin shift
(hyperfine interaction)
Orbital shift
Main features
Information
3 2
nucleus e
I s
electron e n hf spin
2 1
hf
The spin shift yields the local susceptibility near the nucleus: « atomic » resolved susceptibility susceptibility
Spin shift
Main features
Information
hf
hf
Linewidth DH : spatially inhomogeneous susceptibility (dilution) Line shift K : susceptibility frustr
, ,
polarizati core contact i dip i
z y x i
Gyromagnetic ratio Orbital or chemical shift Magnetic (« Knight ») shifts »
Nb of nuclei n reference n0=gH0/2p
Nb of nuclei Local field or frequency or shift H0 YBa2Cu3O7 CuO2
Y Ba Cu O
Ni Ni Ni Ni Ni Ni Ni Zn Ni Ni Ni Ni Ni
One impurity in a Haldane chain, YBa2NiO5 (S=1) with Zn impurities on Ni site
14580 14600 14620 14640 14660 0.0 0.2 0.4 0.6 0.8 1.0
n (kHz)
Y2 Ba Ni98% Zn2% O5 T = 200 K Intensité (unités arbitraires)
pur Zn 2%
Kspin yields a histogram of values, not a sum
Tedoldi et al., PRL 99; Das et al.PRB 04
Nb of nuclei local field
Zn 2%
Zn Ni Ni Ni Ni Ni Ni Ni
0,000 0,002 0,004 0,006
doping
Antiferro T Anti Ferro Supra
0.000 0.002 0.004 0.006
0,000 0,002 0,004
I=1/2
I > 1/2 If I>1/2, nuclear spin I is sensitive to any Electric Field Gradient from the lattice (non-sphericity of the nucleus)
I=1/2
I > 1/2
... 2 1 ) ( ) (
2 _ 3 _ 3
r j i directions i j i r i directions i i
x x V x x x V x V r V eQ z r d R
1 2 3 2
2 3
( )r
Quadrupolar moment of the nucleus
0 since center
charge coincides cst Quadrupole term We express it in principal axes where V is diagonal :
... 2 1 ) ( ) (
2 _ 3 _ 3
r j i directions i j i r i directions i i
x x V x x x V x V r V ) 3 ( ) ( ) 1 2 ( 4 ) ( ) (
2 2 2 2 2 2 2 2 2 2
I I z V I I y V x V I I eQ H r d r V r H
z y x Q n Q
r eQ z r d R
1 2 3 2
2 3
( )r
Quadrupolar moment of the nucleus
Wigner-Eckart theorem I=1/2
2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2
, ) ( 2 1 ) 1 ( 3 ) 1 2 ( 4 ) 3 ( ) ( ) 1 2 ( 4
, 2 2
x V y V z V z V y V x V y x z Q z y x Q
z V eq
I I I I I I I qQ e H I I z V I I y V x V I I eQ H
I=1/2
Quadrupolar nuclei: lifting the multiplicity of transitions on single crystals nNMR = g/2p H0 ~ nQ f(q) nNMR = g/2p H0 nNMR = g/2p H0
2I+1 levels Degeneracy of the transitions lifted by quadrupolar effects
nNMR nQ
EFG Ppal axis
nNMR = g/2p H0
2I+1 levels
Degeneracy of the transitions lifted by quadrupolar effects
nNMR = g/2p H0
nNMR nQ Quadrupolar nuclei: distribution of angles powder average
The EFG and hyperfine tensors may not have the same principal axis! One can manage, playing with isotopes, field …
nNMR = nQ
2 levels
Degeneracy of the transitions lifted by quadrupolar effects
m= 1/2 m= 3/2
nQ 0
18 20 22 24 26 28 30 32 34 0.0 0.2 0.4 0.6 0.8 1.0
x=1,0
×àñòî òà (Ì Ãö)
n (MHz)
Cuchain Cuplane
If I>1/2, nuclear spin I is sensitive to any Electric Field Gradient from the lattice
6.60 6.65 6.70 6.4 6.6 6.8 0.0 0.5 1.0
150 K 200 K 250 K 300 K
NMR Intensity (normalized)
H (Tesla)
300 K
Cu Zn/Cu Cl OH M.I.T., 2005 Herbertsmithite ZnCu3(OH)6Cl2 Cu2+, S=1/2
100 200 300
1 2
10 20
17O lineshift (%)
T (K)
SQUID (10
3/mol Cu)
17O NMR
electronic spins; associated with the existence of a magnetization
so that <S> ≠ 0; in an ordered phase H ~ Ahf <S>, <S> ≠ 0
applied field:
n0Ahf S
n0Ahf S
Zero Field H0=0 Null !
n0n0(1+K) n0+AhfS
n0+AhfS n0-AhfS
NMR under an applied fieldH0
Fe Kitagawa, Takigawa, JPSJ 08
pnictides BaFe2As2
20 40 60 80 100 100 200
(G) T (K)
TN TMIT 0.0 0.4 78500 78750 79000 79250
n (kHz)
Intensity (arbitrary units)
T=80 K
Antiferromagnetic ordering eg: Na0.5CoO2 Splitting of the lines Spin density wave eg: Cr: distributed field
Cu/Co ferromagnetic multilayers
Co surrounded by Co Co surrounded by Co and Cu NArath, Phys Rev 1965 Co entouré de Co Cr in the domain wall Cr inside the domain
Panissod et al., PRB 1992
Work from Panissod Co/Cu mutilayers (1992)
Co-Cu multilayers: Co resonance depends
Cu
Work from Panissod Co/Cu mutilayers (1992)
Grenoble High Magnetic Field Kodama, Science (2002)
T=25mK H=27 Tesla
Static
Dynamics <h+
loc(t) h- loc(0)>
Techniques
For I>1/2, quadrupolar effects, much better with single crystals
~ single crystals
Radiofrequency pulse ~ few msec H0 M H0
H0 U t
Dt
) ( )) ( sin( ) ( n n p n n p n D D t t f f Nb of nuclei
frequency or local field A pulse has a spectral width in Fourier space. Fourier transform yields the response of the sample in the frequency domain of the pulse.
Field range: 1T – 45 T T-range: 10 mK – 1000 K Sensitivity: 1 mMole… depends on sensitivity Misc: pressure (few GPa), in-situ rotation
transverse relaxation : T2
Enegy is conserved H0
Z Z m equilibriu Z
H M T M M dt dM g
1
Y X Y X Y X
H M T M dt dM
, 2 , ,
g
Longitudinal relaxation : T1
Energy exchange with the lattice
z
dt t i B t B T
n L L
exp ) ( ) ( ~ 1
1
Local magnetic fluctuations at n (Fermi golden rule)
) , ( . ) ( ) ( t r S I r A t B
i r nuclei coupled i hf
i
dt t i q s t q s q A T
n hf q
exp ) , ( ) , ( ) ( ~ 1
2 1
( ) ( )
n t n B n
q dt t i q S t q S T k , exp ) , ( ) , ( ) exp 1 ( 2 1
"
T kB
n
n n t q B B
q q A g T k T m ) , ( ) ( 1 1
" 2 2 2 1
Fourier transform Fluctuation
Cu O
A(q) form factor and favours some q.
( )
n n t q B B
q q A g T k T m ) , ( ) ( 1 1
" 2 2 2 1
i
r i i
r q i e r A q A ) . ( ) ( ) (
( )
) cos( ) cos( 2 1 1 2 ~ ) (
2
b q a q q A
y x
g
( )
2 2
) cos( ) cos( 2 ~ ) ( b q a q q A
y x
g
favours q=0, ferromagnetic fluctuations between Cu favours q=p,p, antiferromagnetic fluctuations
Takigawa et al., PRB (1991) Underdoped cuprate
n
Slowing down of fluctuations
In a weak metallic antiferromagnet
Kyogaku et al., JPSJ (1993)
AgVP2S6
T
D
1
Haldane gap
Shimizu et al., PRB (1995)
1 10 100 0.0 0.2 0.4 0.6
1/T1 (ms
T(K)
Frustrated 2DHAF S=1/2
M.-H. Julien, Phys. Rev. Lett. 83, 227 (1999)
1958: Discovery and interpretation by Rudolf Mössbauer 1961: Nobel Prize
→ transmission geometry for bulk samples
→ Surface studies: Conversion Electron Spectroscopy
Born 1929
Conservation of energy and momentum hn = E0 + P2/2Mn hn/c = P Recoil energy ER = P2/2Mn
2 n 2
Mnc2 ~ 100 GeV ; E0 ~ 100 keV; ER ~ 1 meV
Energy conservation:
if ER, no overlap
Optics X rays g rays
At temperature T: moving emitter/abs Doppler broadening D 2 kBT ER for g rays: << ER
For g rays, << ER and D ~ ER weak overlap On cooling, not on heating (D as T ), R.M. observes an increase of the resonant absorption in 191Ir … Interpretation: nucleus bound in a solid ER
2 2
2 1 c M E
solid
= ER
free /NA <<
photon emitted without recoil of the nucleus (for rigid atomic bonds)
Recoilless is for source and absorber!
finite probability f(T) of nuclear resonant absorption
Mössbauer spectroscopy of hyperfine (electro-nuclear) interactions (~10-6 eV) if << hf Lamb-Mössbauer factor: f(T) = |J0(kx0)|2 1 ½ k2x0
2
f(T) 1 k2<x2>T exp( )
Vary the speed of the source ~ sweep the frequency (Doppler) 1 mm/s ~30 MHz
46 elements, 89 isotopes, 104 Mössbauer transitions
1/2
1/2 3/2
1/2 ±3/2 ±1/2
1/2
1/2 3/2
1/2
1/2 3/2
Symmetric charge No magnetic field Asymmetric charge No magnetic field Symmetric or asymmetric charge Magnetic field (internal or external) Δ Bhf δ Isomer shift Quadrupole splitting Magnetic hyperfine field
2 4
Velocity (mm/s) Relative Transmission
D
2 4
Velocity (mm/s)
1.0 1.5 1
0.5 Isomer shift (mm/s) 2 3 4 0.0
[6]Fe(II) [6]Fe(III) [6]Fe3+ [4]Fe3+ [6]Fe2+ [4]Fe2+ [sq]Fe2+ [8]Fe2+ [5]Fe3+ [5]Fe2+
Isomer shifts and quadrupole splittings of Fe-bearing phases vary systematically as a function of Fe oxidation, Fe spin states, and Fe coordination. Knowledge
the Mössbauer parameters can therefore be used to “fingerprint” an unknown phase.
(Murad et al.)
intermediate spin state
first indication in favour of a commensurate Spin Density Wave
(Klauss, Luetkens et al.)
(Klauss, Luetkens et al.)
Note: effects are not the same on all lines: outer are more « protected »
Slow relaxation
Intermediate relaxation
S=5 AF interactions 3 ineq sites, EFG 0
155Gd in Gd3Ga5O12
No magnetic order
geometrical frustration
P.Bonville et al PRL 92 (2004) 167202
170Yb
2
Which sample?
Needs a source Many…needs time
Fluctuation rate
Few 10 GHz… MHz 100 MHz – fraction of Hz
Location/coupling
0.1 T – 10 T/ µB
0.1 T – 10 T/ µB
Observables
Magnetic transitions Magnetic susceptibilities
Temperature range
10 mK – … K 10 mK – 1000 K
Field range
0 – (few T) 1 – 45 T
Intrinsic drawback
Need a source r.f. field needed, field needed Tuning of the probe
B e e
e e
2 24 Am
e B
E H ms=+1/2 ms=-1/2 (1/2)geH (-1/2)geH hn = geH Dms= 1
ESR and NMR comparison!
electron proton ratio
Rest mass me =9.1094*10-28 g mp =1.6726*10-24 g 5.446*10-4 Magnetic dipole moment mS=-ge meS ge= 2.002322 me=eh/4pmec = 9.274*10-21 erg/G mS=-gN mNS gN= 5.5856 mN=eh/4pmNc = 5.0504*10-24 erg/G 1836.12
Frequency: Factor 1000 larger in EPR ! (GHz instead of MHz) Dipolar coupling: Factor 1 000 000 larger in EPR ! (MHz instead of Hz) Relaxation Times: Factor 1000 000 smaller in EPR ! (ns instead of ms)
= much higher techniqual requirements, but unique sensitivity to molecular motion
Sensitivity : Factor 1 000 000 better than in NMR !! (1nM instead of 1mM )
An ideal case, though
– X band; ~9-9.5 GHz, in most widespread use (l~3 cm). – K band; ~ 24 GHz (l~1 cm) – Q band; ~ 35 GHz (l~0.8 cm) – W band; ~ 95 GHz
– At g=2, about which most spectra are centered, X-band setups have resonances at 3,000-3,500 Gauss.
Use of a cavity (except at high frequencies)→ sweep the field
HeZ = Electron Zeeman interaction: g tensor HZFS = Zero-field splitting interaction: anisotropy, dipolar Hee = Interactions between electron moments: exchange Hen = Electron - nucleus interaction HNZ = Nuclear zeeman interaction Hilbert space of coupled electrons and nuclear spins has a dimension
n n m m H
4f: use gJ (free atom) instead of g (takes into account the spin orbit term)
B e B e LS eZ
B e B e LS eZ
2
n n j n n i ij e
ge=2.0023
hexagonal): gxx= gyy = g and gzz= g||
zz yy xx
g
hn = gxxmBB, B//x, hn = gyymBB, B//y hn = gzzmBB, B//z.
B eZ
99
E
B eZ
B // z (gzz = g//) B z (gxx = g)
S=1/2, I=0, gx=gy<>gz Axially symmetric g-factor
2 / 1 2 2 2 2
] sin cos [
q q m n m n g g h g h B
II B B eff r
q is the angle between a z-principal axis and the magnetic field direction The given solid angle W is defined to be the ratio of the surface area S to the total surface area on the sphere: W S/4pr2: dW/W2pr2sinqdq/4pr2 sinqdq/2
q qd dB B f sin ) ( q cos / 1 ) ( d dB B f
q q q n m cos ) ( ) sin cos ( ) (
2 2 2 / 3 2 2 2 2
g g g g h B f
II II B
q n m cos ) ( 1 ) (
2 2 3
g g B h B f
II r B
Absorption B g// g
gxx=2.0507; gyy=2.080; gzz=2.230.
Absorption Observed
B ZFS
E=0: 3 states, singlet(1) and triplet (2) All these terms resume in a quadratic form of Sx, Sy, Sz
2 2 2 2 2 2
y x z z z y y x x ZFS
) 1 3 3 )( 1 ( 5 1 6
2 4 4 4 '
S S S S S S S a H
z y x ZFS ZFS
(S > 2)
On a powder, broad line due to anisotropy
e
g
e
g
E = gmBBoMS + aMSmI
n e
Isotropic case MS=-1/2 MS=1/2 mI=1/2 mI=-1/2 mI=-1/2 mI=1/2 Zeeman only Hyperfine interaction
1𝑇 2
S=0 S=1 H
gµBS Only anisotropic part contributes to M2- e.g. dipolar, DM
under isotropic coupling.
+1/4 J
i i i i i
1 B
Zeeman energy isotropic exchange additional interactions e.g. crystal field anisotropic exchange dipole-dipole interaction hyperfine interaction
strong isotropic coupling averages local fields similar to fast movements of the spins “exchange narrowing“ of the ESR signal
local inhomogeneous fields local, static resonance shift inhomogenous broadening of the ESR signal
ij
j i ij
( )
5 2 2 2
3 2 1 ) (
ij ij ij ij B ij
r r r r g dip K
m
B pp 2 4 3 2
anisotropies completely contained in the second Moment M2: remaining task: calculate the second moment for the different contributions to the spin Hamiltonian, find the dominating line-broadening mechanism (and check for the anisotropy) Uncoupled spins
Special features:
(dipolar contribution suppressed by strong exchange coupling).
Similar to NMR with S = total spin of ferromagnet. Magnetic selection rule: Δ mS = 1.
Consider an ellipsoid sample of cubic ferromagnetic insulator with principal axes aligned with the Cartesian axes. Bi = internal field . B0 = external field. N = demagnetization tensor
i
B B M N
i j i j i
N N
Lorenz field = (4 π / 3)M. Exchange field = λ M.
( )
i
d dt g g M B M B M M N
Bloch equations: →
0 ˆ
B B z
( )
x y z y
dM B N N M M dt g
( )
y x z x
dM B N N M M dt g
z
M M
i t k k
M e M
→
( )
( )
y z x y x z
i B N N M M M B N N M i g g
ˆ ˆ ˆ
x y x x y y z
M M M N M N M B N M x y z
( don’t contribute to torque)
( )
( )
2 2 y z x z
B N N M B N N M g
FMR frequency: uniform mode
( )
( )
2 2 y z x z
B N N M B N N M g
For a spherical sample,
x y z
N N N
→
B g
For a plate B0 ,
0, 4
x y z
N N N p
→
( )
4 B M g p
For a plate // B0 ,
0, 4
x z y
N N N p
→
4 B M g p
Polished sphere of YIG at 3.33GHz & 300K for B0 // [111]
Shape-effect experiments determine γ & hence g. Fe Co Ni g 2.10 2.18 2.21
(H. Hurdequint)
ESR and NMR comparison!
electron proton ratio
Rest mass me =9.1094*10-28 g mp =1.6726*10-24 g 5.446*10-4 Magnetic dipole moment mS=-ge meS ge= 2.002322 me=eh/4pmec = 9.274*10-21 erg/G mS=-gN mNS gN= 5.5856 mN=eh/4pmNc = 5.0504*10-24 erg/G 1836.12
Frequency: Factor 1000 larger in EPR ! (GHz instead of MHz) Dipolar coupling: Factor 1 000 000 larger in EPR ! (MHz instead of Hz) Relaxation Times: Factor 1000 000 smaller in EPR ! (ns instead of ms)
= much higher techniqual requirements, but unique sensitivity to molecular motion
Sensitivity : Factor 1 000 000 better than in NMR !! (1nM instead of 1mM )
An ideal case, though
Very powerful, quite involved treatment, if detected: single crystals needed other information modeling
Spin waves of odd number of half- wavelenths can be excited in thin film by uniform Brf
( )
2
4 B M Dk g p
Condition for long wavelength SWR:
D = exchange constant
( )
2
4 n B M D L p g p
For wave of n half-lengths:
Permalloy (80Ni20Fe) at 9GHz
Consider a uniaxial antiferromagnet with spins on 2 sublattices 1 & 2. Le t
1
ˆ
A
B M z
BA = anistropy field derived from
2 1
sin
K
U K q
θ1 = angle between M1 & z- axis.
2
A
K B M →
1 2
M B B
2
ˆ
A
B M z
Exchange fields:
( )
1 2
ex l B M
( )
2 1
ex l B M l
For
a
B
1 2
ˆ
A
B l B M z
2 1
ˆ
A
B l B M z
1 2 z z
M M M
With the linearized Bloch equations become:
( )
( )
1 1 2 x y y A
dM M M B M M dt g l l
( ) ( )(
)
2 2 1 x y y A
dM M M B M M dt g l l
( )
( )
1 2 1 y x x A
dM M M M M B dt g l l
( )(
)
( )
2 1 2 y x x A
dM M M M M B dt g l l
x y i t j j j
M M i M e
→
( ) ( )
1 1 2 A
i M i M M B M M g l l
( ) ( )
2 2 1 A
i M i M M B M M g l l
( ) ( )
1 2 A E E E A E
B B B M B B B M g g g g
E
B M l
exchange field
( )
2 2
2
A A E
B B B g
AFMR frequency
., Principles of Magnetic Resonance, Springer Verlag, 1978
Press, New York, 1967) Chap. 7
. Carretta, A Lascialfari NMR-MRI, muSR, and Mössbauer spectroscopies in molecular magnets, Springer, 2007
England: Oxford University Press, 1970.