Bloembergen-Spin Relaxation Processes in A Two-Proton System
Bloembergen-Spin Relaxation Processes in A Two-Proton System
The general theory of nuclear spin relaxation, based on the Boltzmann transport equation for the density
I
matrix, is applied to the very simple, nontrivial system of two identical spins = ~. A proton pair undergoing
hindered rotation around one axis shows a resolvable doublet. The relaxation processes in this doublet are
investigated. Explicit formulas for the longitudinal and transverse relaxation times are derived, which
display a dependence on the angles of the axis of rotation with the external magnetic Geld and the radius
vector. Experimental evidence for this dependence is discussed. General expressions for the Overhauser
and saturation e6'ects with two applied radio-frequency Gelds are given, which may also be applied to a
I
magnetic ion or nucleus with = 1 and crystalline Geld splitting.
in the primed coordinate system. The coefficients in this Here 0~~) is the diagonal spin density matrix corre-
transformation are the irreducible representations sponding to thermal equilibrium at the lattice tem-
n&2&(o, fp, o) of the group of three-dimensional rota- perature T.
"
tions, described by Wigner. With the introduction of
the raising and lowering operators I+=I,&I„, one can o &r&=exp( —E /kT)/P exp( —E /kT). (13)
write
Equation (12) gives a complete dynamical description
xd v +f k r12 (I1 I2 3I1 I2 )I ' p (0) (6a) for any magnitude of the external field Ho, including
zero. It contains all interference and nonsecular per-
7 k r12 (Il I2 +Il I2 )++1 (0 4') (6b) turbation eBects between the five independent elements
—1y h, r12 2I1+I2+Py21 (8 $) ~
of the density matrix, whose trace is normalized to
(6C)
unity. The five independent coupled differential equa-
with
tions, represented by the operator relation Eq. (12),
2
z„'(e,y)= p n„„. would still be difficult to solve, but the problem is con-
(O, p„o)z . (s', y'). (7)
mf 2 siderably simpler than the complete motion of two
nonidentical spins, which is described by nine coupled
The P&') are the unnormalized spherical harmonics
equations. "
Pp&2& =-'(3 cos28 —1), Considerable simplification results if nonsecular per-
turbations are neglected and the time-independent
P~~&" = ~ sin8 cosine+'& Hamiltonian is diagonal in the representation chosen.
P &') =3 sin'Oe+"& The latter situation occurs when the external field is
large, 3CO&&K~;„and only first-order perturbation
In this form all matrix elements in the (+0 —) repre- theory to the time-independent part of Kp'p is applied.
sentation can be written down at once, and the time- One obtains the result of Gutowski and Pake':
independent part is clearly separated. It consists of
the terms 222'=0 which do not depend. on p'. The time-
—4ypkpr 12—'(3 cos2fp —1),
Ey = yhHp
dependent matrix elements occur in the form Ep +2y'k'r12 —
'(3 cos'Pp —1),
.
ac. .(t)=P, K. .' H(pt), E = —ykHp —41y2kpr12 '(3 cospp —1).
which is used in the theories of Bloch and Redfield. Redfield' shows that application of second-order per-
The E& are Hermitian spin operators which do not turbation theory with the random time-dependent part
contain the time explicitly and the H'(t) depend only leads to the following equation for the density matrix:
on the time-dependent lattice coordinate. The Hp(t) can
be chosen as the four real functions cosp', sin&', 2 cosp'
and sin2tt '. With random time variations of p', introduce
= —2k(E- —E- )o- +Z R- pp opp
Bt pp'
the generalized correlation functions and spectral densi-
x(~, ~', p, p'=+, 0, —). (15)
ties of the functions Hp(t) by
Only secular perturbations will be retained. The secular
k„(~)=-', ~ p'-(H (t)H*"(t —r))dr. elements of the relaxation matrix, which satisfy the
condition
i
E Ep —
R p— E+Ep
p (16)
i
and the nonvanishing spectral densities are wp —(pP) Py4N'iq 'r {2 simp cosQp sin'8' cos'0'
+8 simp sin'8'), (29)
k„(&p) = r/(1+to'r'). (21)
(b) There are a very large number of equilibrium ~P'5 rip '(1 cosQp) —
(
positions and a stochastic diGusion process describes
wi sin'8'~
&1+ '&o — p'r'&
the motion in azimuthal angle. The probability to find —3 cosQp+1)
the pair at an angle tt)p'+t)It' at time t, when it was at tttp' + ,'py4htrip '—(4cosQp
p(y', t) =
2(~tr ')'-
The correlation functions are in this case
exp( @'r/4t)dy'. —(22)
= same
9
expression with co+0
E1+~ r
replaced by oro,
)
(30)
&(sin't)' cos'0'
32
cosgp)
( 2r
(31)
Again the spectral densities vanish for q &q'. ~
(1+&p+ 'r') ~.
kq, ; =~i r/(1+to' r'), — (24a) These formulas hold for the case of stochastic di6usion
&q, q =+q=-'r/(1+ 's~'r') — (24b) of the azimuthal angle. They can immediately be
changed to the case of random jumps between three
The nonvanishing spin matrix elements E p& are equilibrium positions by using Eq. (21) rather than
obtained from Eqs. (6), (7), (8) and (9). (24). The small difference between (o~p and top in the
SP I N RELAXATION P ROCESSES
expression for m~ has been ignored. The 5)&') functions This corresponds to a partial application of the general
have been taken from Wigner. The normalized constants operator formulation, Eq. (12), which contains all
have been checked by direct trigonometric transfor- interference and nonsecular effects. Equations (35)
mation of the functions (8), which led to the results describe the frequency pulling and damping by cross-
about as quickly as the use of Wigner's general formula relaxation of off-diagonal elements which have nearly
in this case. the same frequency. The coupled system can readily
The set of Eqs. (26) describe the longitudinal relaxa- be solved for the two normal modes. This situation,
tion. There are two characteristic times, as could be which illustrates the principles and limitations involved
expected for a problem with effective spin 1. The in the distinction between secular and nonsecular per-
"magnetic moment relaxation"" is characterized by turbations is rather academic. The overlap of resonance
~( ++ — ) — 1
(~++ —o —)=0++Ti (32)
lines is the experimental counterpart of this theory of
cross-relaxation. In the case of a proton pair in a solid.
2fop+ tffi lattice the e6ect described by Eqs. (35) will be obscured
completely by the interaction from neighboring moments
The "quadrupole moment relaxation" is described by outside the pair.
ft(o+++o —2o pp) 1
(o'+++ 0' —20'pp) GI. SATURATION AND OVERHAUSER EFFECT
3'Ny Consider the case that two radio-frequency magnetic
o'+~+0 20'pp
fields are applied. Hi(v+p ) has a frequency close to the
(33)
resonant frequency v~p, H, (vp ') is close to the resonance
frequency vp of the other line. Assume that the am-
The intensity of the magnetic resonance absorption
—
lines is proportional to 0-++ C.pp ol Opp — and will
in general approach an equilibrium according to a
0, plitudes of both fields satisfy the relation
Ep- op/Tp'
to secular character takes place. The term can be
''
=
op-
carried along in a rigorous manner provided the cor- ffop o'p ZMp o'p
relation time r, is short compared to k(E+ Ep — -opp)
+ ', iyHf(vp ')V2(o (38b)— .
+E ) '. A set of two coupled equations for the off-
diagonal elements 0+p and o=p should be solved These equations replace Eqs. (27). There are corre-
sponding complex conjugate equations for 0-p+ and 0. p.
a~+ p/at=icoypo'yp+~~po'~p+~pp The steady-state solution for the diagonal com-
(35)
80'p /Bt=zÃp op +Ep p op+Re +po+p. ponents is obtained by putting the left-hand side of
'P F. Lurcat, Compt. rend. 240, 2402, 2517 (l.955). '4 F. Bloch, reference 8, Sec. III.
1546 N. BLOEM BERGEN
Eqs. (26) equal to zero and adding a term "
state Maser. In general, transitions between the +
— and —level 'could also occur. They can readily be
ss7H)(v+p')v2(o~ ~+p)
incorporated in the analysis.
to the right-hand side of (26a), The treatment presented here has to be modified if
one of the radio-frequency fields becomes larger than
,'i—
y—Hr (vp+')%2(o p a—
p) the intrinsic line width, yII~&T2 . A rigorous dis-
to (26c) and subtract the sum of these terms from the cussion has been given for the case than one radio-
right-hand side of (26b). Solution of the set of steady- frequency Geld is very large, while the other is relatively
state equations leads, by elimination of the off-diagonal weak. Transform to a rotating coordinate system to
elements, to eliminate the explicit time dependence of the strong
field. Red6eld" has discussed the transformation of the
dipolar interaction and relaxation terms in the rotating
system. The magnetic resonance response to the weak
field is split into components. Their spacing and relative
amplitudes depend on the value of the effective field
ACKNOWLEDGMENTS
' R. E. Richards (private communication). The author is
indebted to Dr, Richards for discussions on the experimental The author is indeted to Professor E. M. Purcell and
situation, which led to this work. Dr. A. G. Redfield for reading the manuscript and
20
J. G. Powles and H. S. Gutowski, J. Chem. Phys. 21, 1695
(1953). making valuable suggestions for improvement.