
книги из ГПНТБ / Гречишкин В.С. Ядерные квадрупольные взаимодействия в твердых телах
.pdfн о |
КВАДРУІІОЛЬНАЯ РЕЛАКСАЦИЯ |
|
[ГЛ. XV |
||
Ж — полный гамильтониан, Ѳ(t) = |
f 1 |
при |
t ]> О, |
||
{ n |
при |
J |
. _ При та- |
||
|
' |
[О |
t |
< 0 . |
ком подходе не учитывается насыщение.
Однако в ЯКР имеется одно существенное отличие от
ЯМР, где Жа— — ^Н0%2 I zi- В магнитном резонансе пере-
І
ход к представлению взаимодействия есть преобразование вращения. В квадрупольный гамильтониан входят квад
раты операторов 1%. Поэтому здесь необходимо так из менить схему расчетов, чтобы были применимы обычные правила унитарных преобразований. Эту трудность мож но обойти с помощью выбора подходящего представле ния [10].
Рассмотрим четырехранговую матрицу, каждый квад рант которой является в свою очередь двухранговой мат рицей:
' i f i У <4 -в >
Пусть матрицы а, и, I и Ъ могут быть представлены через матрицы Паули:
0 1 |
|
|
|
|
0 - |
і |
|
1 |
0 |
|
|
|
|
1 |
0 |
|
|
|
|
|
|
1 о |
(4.7) |
4 L |
|
|
|
|
|
|
|
|
|
|
|
|
0 1 |
|
|
Определим матрицу оах так: |
|
|
|
|
|
||
|
|
'0 |
1 |
0 |
(Г |
|
|
<5„ѵ = |
1 0 |
0 |
0 |
, |
(4.8) |
||
0 0 |
0 0 |
||||||
|
|
,0 |
0 |
0 |
0 |
|
Тогда для J — 3/2 матричное представление оператора І г принимает вид
72 |
0 |
0 |
0 |
|
|
0 |
7 2 |
0 |
0 |
~2~(aaz + 2еа -f- abz — 2еь). |
(4.8') |
0 |
0 |
- 7 2 |
I = |
||
0 |
|
|
|||
0 |
0 |
0 |
- 7 2 |
|
|
§ О |
Д В У Х У Р О В Н Е В А Я СИСТЕМА |
|
111 |
|||
Отсюда |
гамильтониан |
квадрупольного |
взаимодействия |
|||
можно |
записать |
в виде |
|
|
|
|
|
|
hQ = |
3К (оа2 — сЬ2), |
|
(4.9) |
|
где К = |
eQqJ\2. |
|
р (t) |
возьмем в |
представлении |
|
Матрицу плотности |
||||||
взаимодействия? |
|
|
|
|
|
|
|
р* (t ) = |
Тр (*) Г " 1 = |
e^Q'p (t ) e - i x Q \ |
(4.9') |
Тогда для р* (t) основное уравнение (2.2) запишется в ви де
|
|
|
г - С ^ = |
[ЗГ,Р*(*)], |
|
|
(4.10) |
||
где |
%= 1, |
Ж’ = |
е ^ ^ Ж е '1^ |
1- |
Жр. |
|
|
|
|
Если в гамильтониане (4.1) |
опустить члены Ж2и Жad, |
||||||||
то в |
выбранном |
представлении |
|
|
|
||||
Ж = |
з я 2 |
( 4 - |
4 ) - |
2т-Я, 2 |
(öL + |
4 ) + |
|
||
|
і |
|
|
|
г |
|
|
|
|
|
|
|
|
|
+ |
4" (a«- + |
4 |
) cos |
4 11) |
где |
a+ = |
ax + |
iay, |
Ö)0 — резонансная |
частота |
системы, |
— 2 уНх cos Wo t — линейно поляризованное радиочастот ное поле.
Из соотношений |
(4.1), |
(4.8') |
и (4.9') |
находим |
||
Т |
= |
ехр |
Ш |
2 ( 4 |
— 4 ) } |
= ПТ ь |
|
|
|
^ |
г |
J |
і |
|
|
|
|
|
|
(4.12) |
T i |
= |
exp j- |- m |
t (6eZ— 6bz)| . |
|
Если Ql есть некоторый оператор, действующий только на і-спин, то TQ'T~l = TiQ'Tf1. Поэтому, если матрица Q{имеет только квадранты а и Ъ(см. (4.6)), с нулями в квад рантах и и I, то
icofa* |
Qa ехр\ |
|
|
|
|
TQ'T-1 = ехр 2 |
|
|
|
||
+ |
ехр ( — |
fatal |
г®г4 |
. (4.13) |
|
2 |
2 |
||||
|
|
|
112 |
К В А Д Р У П О Л Ь Н А Я РЕЛ А К С А Ц И Я |
[ГЛ . IV |
Первый член в уравнении (4.13) дает вращение Ql вокруг
направления оси Z, а второй — вращение Ql в противоноложном направлении. Итак, линейно поляризованное по ле в (4.11) представляется суммой двух циркулярно поля ризованных полей с противоположным направлением вра щения. Каждое из этих полей взаимодействует лишь с теми спинами, прецессия которых совпадает по направлению с вращением соответствующего поля. Используя (4.11), по лучим
ЗГ = |
[ък - |
4 - ю) |
2 (*« -°ьг)-- |
Т1 3ТЯх 2 |
(бах + |
<&) + |
|
|
' |
' |
г |
г |
|
|
|
|
+ Жм + члены вида еш + |
члены вида |
е 2 |
Ш , |
(4.14) |
||
где |
Ж*м — секулярная часть диполь-дипольного |
взаимо |
действия, не зависящая от времени. При точном резонансе 3 К — Ѵ2 со = 0. Если мы опустим в (4.14) члены, завися щие от времени, т. е. несекулярную часть гамильтониана диполь-дипольного взаимодействия, то это означает, что мы интересуемся лишь линией поглощения на резонансной частоте [11]. Несекулярная часть диполь-дипольного взаи модействия дает вклад в ширину линий ЯКР на частотах 2ю0, Зю0 и т. д. *).
Диполь-дипольное взаимодействие между |
парой |
спи |
||||
нов можно записать в виде |
|
|
|
|
||
hik — Y*7iJt lAihlzilzh + |
B ih (/+j/-h + |
I-tl+h) |
+ |
|
||
+ |
I+ilzh) + Dik (lzil-h + |
I-ilzk) |
+ |
(4.15) |
||
где |
+ |
E ihI +tI , k + |
FihI |
^ |
h], |
|
|
|
|
|
|
|
|
Aik — 1 |
— 3cos2 Ѳ, |
Bik ------- (1 — 3cos2 Ѳ), |
|
|||
Cik = Dik — ---- sin Ѳ• cos Ѳ• e~iv, |
|
|
(4.16) |
|||
|
E ik = Eu = |
---- sin2 Ѳ• er2'*, |
|
|
|
Ѳ и ф — сферические координаты радиуса-вектора, соеди
*) Экспериментально поглощение на кратных частотах в ЯКР не наблюдалось. Исследование кратных резонансов позволяет су
дить о несекулярной части гамильтониана диполь-дипольного взаи модействия.
§ и |
Д В У Х У Р О В Н Е В А Я СИСТЕМА |
113 |
няющего спины і и к. Переходя к представлению взаимодей ствия, получим для секулярной части
к ш — |
У 1"ік ^ y - ^ i k h i h k + |
|
|
|
+ |
B i к |
+ |
öa_ — 5a+ -f- 5b+<3b- -f- Oj,_Öb+) -f- Ou-3; +~j- |
|
+ |
sl+ •öu- |
C ik { h 0 |
l - 4- a u - h k ) + C ik |
(l z0 f + + ° l + h k ) |
|
+ E. |
(Ga+Cb+ + öb+öa+) + |
öu-6«-^ + |
+ Я « [— (б«Д5ь- + Gb_Oa_) -j- Oi+Oi+J , (4.17)
где
M d â — 2 h k - i > k
Знание секулярной части диполь-дипольного взаимодейст вия позволяет вычислить второй момент линии ЯКР [11]:
S , |
Sp{[*&, У 2} |
(4.18) |
|
SP(^}
Такие расчеты по аналогии с ЯМР были проделаны для
различных |
спинов [12, |
13]. |
|
Если / |
= 3/2, то для простой кубической решетки полу |
||
чим |
S2 = |
60 ^ ß 4d-6, . |
(4.19) |
|
где g — гиромагнитный фактор ядра, ß — ядерный магне тон, d — длина элементарной ячейки кристалла.
Пусть на спин-систему действует радиочастотный им
пульс. Предположим, что член Жх Жаа, поэтому пре небрегаем диполь-дипольными взаимодействиями в период действия импульсов. При t = 0 начальное значение мат
рицы плотности |
имеет |
вид |
|
|
|
е-№ о |
( |
3ßL (0) |
К |
(4.20) |
|
Р(0) = |
Spp |
S p l |
—®bz)> |
||
|
|
где ßb (0) = 1ікТъ — обратная температура решетки. По скольку при р(0) = р*(0) спиновая температура постоянна,
114 К В А Д Р У П О Л Ь Н А Я РЕЛ А К С А Ц И Я [ГЛ . IV
то |
из (4.9) |
и (4.13) |
получим |
|
|
|
|
- |
3ßL (0) к |
exp iant |
aaz exp ! |
|
|
р* (0 = |
S p l |
Oax\ — |
||||
|
|
|
|
|
|
|
|
|
— exp (■—— öbxj öbz exp |
icojf Зъх) |
, (4.21) |
||
где |
<»! = / |
3yHv |
хотя в случае ЯКР |
при / = |
3/2 мы |
|
|
Таким |
образом, |
имеем две пары дважды вырожденных уровней, поведение
спин-системы аналогично случаю / |
= Ѵг в ЯМР. В конце |
||
действия 90°-импульса (tw = пІ2а>^ |
|
|
|
Р* (U = ( - -\ Т |
) S |
(a“v - 4 ) . |
(4.22) |
|
i |
|
|
Спиновая температура становится функцией времени і при t > tw.
Гамильтониан спин-решеточного взаимодействия запи
шем |
в виде |
[1] |
|
|
|
|
|
ж2 (t) = 2 h2i (t) = 2 |
F™ (t) AW, |
(4.23) |
|
|
|
I |
g |
|
|
где |
А2і (£) = |
Fi~g) (t) A\e). |
Будем |
считать, что вклад маг |
нитных диполь-дипольных взаимодействий за счет флук туаций решеточных координат мал. Функция F(~8> (t) — случайная функция решеточных колебаний, включающая флуктуирующие градиенты электрического поля. Флук туации градиента электрического поля могут возникать за счет вращательных качаний молекул, смещений ионов, из менений степени ковалентности химических связей, затор моженных вращений *). AW представляет операторы квадрупольного взаимодействия.
Если все спины в образце занимают эквивалентные по ложения, то гамильтониан (4.23) можно записать в виде
+ 2 |
|
h2(t)= 2 F(-e)(t)AW, |
(4.24) |
g=-2 |
|
*) Как будет показано ниже, здесь следует различать случай вращений молекул как целого и вращения фрагментов молекул.
§ И
где
Д в у х у р о в н е в а й с и с Дё м а |
115 |
Аі0) = 3l l - J ( J + |
1), |
А{±1) = 4 - / 6 - (I ZI ± + |
I ±I Z), |
||
л ( ± в = 4 - / б ‘4 , |
^ (о)= 4 - ^ > |
|
|||
• ^ ( ± 1 ) |
( £ ) = |
- |
р — =• |
(?3C Z r t * 7 l / z ) > |
( 4 . 2 5 ) |
Я±2) ^ = |
2 V I |
' ^ xx ~ Qvv — 2iqxv)' |
|
Компоненты тензора градиента электрического поля явля ются случайными функциями времени.
При решении релаксационной задачи воспользуемся ос новным кинетическим уравнением, полученным из урав нения Лиувиля (2.2) во втором порядке теории возмуще ний для диагональных элементов матрицы плотности *) [1, 14]:
- 4 г б * ( 0 = ------ |
Y s J A<*P) [Bpg\ [ Bf, в ( O l ] , ( 4 .2 6 ) |
|
g,p |
а* (t) = p* (t) — p* (o c ).
Это уравнение справедливо, если в системе можно выделить малый интервал времени тс. Впервые на важность такого условия в теории необратимых процессов обратил внима ние Боголюбов [7], введший понятие естественной времен ной шкалы системы. Роль тс в нашем случае выполняет время корреляции, в течение которого случайная функция не успевает существенно измениться. Основное кинетиче ское уравнение в форме (4.26) справедливо при тс Тг (для так называемой гидродинамической области релакса ции). Поскольку мы не учитываем недиагональных элемен тов матрицы плотности, то это также предполагает условие Тх ^> Т2, что довольно часто выполняется в твердых телах. Наше приближение, следовательно, эквивалентно прене брежению членами высоких порядков при разложении по степеням безразмерной величины тJ T X.
*) Метод матрицы плотности тесно связан с теорией возмуще ний. Аналогичное рассмотрение может быть сделано й для релакса ции </*> и </„>.
і і б К В А Д Р У П О Л Ь Н А Я РЕЛ А К С А Ц И Я £ГЛ. IV
В |
состоянии теплового равновесия |
ßs = ßt» поэтому |
||
матрица плотности определяется р (оо). |
||||
В уравнении |
(4.26) В (р8) так |
связан |
с А іе) (t): |
|
|
AW (t) = |
ехЖ^ А ^ е ' гЖ^ |
= 2 |
(4.27) |
|
|
|
V |
|
где |
Bp) — 2 B(s) (О- Суммирование здесь проводится по |
|||
|
І |
|
|
|
собственным частотам спин-системы; w(pg) — частоты реше точного спектра. Введем понятие о функции корреляции
ё8(т) = FW(t)FW{t + X). |
(4.28) |
|
Это — вещественная и |
четная функция т. Тогда |
фурье- |
образ этой функции |
будет |
|
|
ОО |
|
Jg(сорг)) = |
^ gg(т)ехр( ш р \) dt. |
(4-29) |
|
—оо |
|
Поскольку в рамках сделанных предположений спиновая
температура |
является медленной функцией времени, то |
|
|
P ( 0 - ( - ^ T ^ ) | (4 - ö Ü , |
(4.30) |
что эквивалентно уравнению (4.22). Отсюда |
|
|
о* (t) = |
а(о = ( - ^ - ) [ßs (о - ßL] 2 (4 - |
4 ) . (4.31) |
Это уравнение указывает на наличие релаксационного процесса в системе.
Запишем гамильтониан спин-решеточного взаимодейст
вия в |
представлении |
взаимодействия: |
|
||
h’2(t)= Th2{t) |
= |
|
|
|
|
= pm (о 3 (aaz - |
abz) + |
pi-v (t) |
|
- ab+e-W) + |
|
|
+ Я+« (0 |
(ба-еч“»‘ — <5ъ-еіы) + |
|||
+ |
F{-2)(t) ( - y r ~ ) Ke« + 6uz) |
+ (eu - |
auz) er**] + |
||
+ |
/A+2) (0 ( ' Y f 7') К8/ — *iz) еіш* + |
(e, + |
elz) er**]. (4.32) |
§ Ü Д В У Х У Р О В Н Е В А Я СИСТЕМА , Щ
Анализируя уравнение (4.27), можно определить опе
раторы В ^ на основе уравнения (4.32). Используя (4.31) и (4.32), вычисляем коммутаторы в уравнении (4.26); они дают либо нуль, либо пропорциональны (оаг — аЬг). Эти множители пропорциональности, будучи просуммирован
ными, дают ЗК/ТХ. Поскольку |
коммутирует с о* (t), |
||||
то решеточные колебания, связанные с F(0\ не вносят вкла |
|||||
да в |
релаксацию. |
|
|
|
|
Окончательно |
получаем |
|
|
|
|
|
|
Т? = 2 ( W 1 + |
W,), |
(4.33) |
|
где |
Wx = 18 J x (to0) для Am = + |
1 (m — zb V2 <-» m = |
|||
— + |
3/2), W2 = 18 / 2 (co0) для переходов c Am— + 2 (m— |
||||
= + |
Ѵг <-* m — zb 3/2). Спектральные плотности J опре |
||||
деляются уравнением (4.29). |
|
релаксации |
|||
Конкретные |
механизмы |
квадрупольной |
|||
были рассмотрены в ряде |
работ |
[15—18]. |
|
Байер [15] рассмотрел влияние вращательных качаний молекул на время Тх в молекулярных кристаллах. Функ
ции корреляции g (т) он взял |
|
в виде |
|
|||||
g. (t) = (Dl - |
^ W |
^ e x p |
(----Щ-) + |
|
||||
|
+ |
2 |
(Dl - |
^ W |
|
^ e x p ( - |
i l l ) , (4.34) |
|
|
|
r=l |
|
|
|
|
' |
r ' |
n r |
T„ |
|
1 |
|
h v t |
, v< частота вращатель |
||
где — = e TX, |
— = - |
_ 1, x--= |
|
|
ных качаний молекулы в решетке, Ѳ2 (£) = Ѵ2 Z)2 — угло_ вое смещение оси градиента Z в процессе вращательных ка.
чании, и г = _L_1 / |
-—=г-----— момент |
инерции молеку- |
|
2л Г |
J Гt |
|
— число моле |
ли, п — полное число молекул в 1 ел8, |
|||
кул в r-колебательном состоянии; |
|
|
|
gi (г) = -у- l^o cos (2яѵ,т) -у- exp J— І І І -J + |
|
||
Oo |
|
|
|
+ [ 2 |
Dl - f r exP ( — ^ |
) c°s (2яѵ,т)]}. (4.35) |
Если положить время жизни в возбужденном состоянии x r — х а = const, то можно получить после несложных
118 |
К В А Д Р У П О Л Ь Н А Я |
РЕЛ А К С А Ц И Я |
[ГЛ . ІѴ |
вычислении |
|
|
|
|
2 (ch X— 1) |
2ch X — 1 |
1 -X |
1 + 4я2ѵ2т2
(4.36)
-1
где та — среднее время жизни кванта вращательных кача ний, VQ — частота ЯКР. Таким образом, при вычислении времени спин-решеточной релаксации Байер считает, что все возбужденные состояния имеют одно и то же время жизни, и переход из одного состояния в другое всегда про исходит с предварительным возвращением в основное сос тояние. Естественно, эти предположения не всегда будут корректными.
Уосснер и Гутовский [16] несколько обобщили схему Байера, решив уравнения для заселенностей вращатель ных уровней энергии. Если молекула совершает качания в решетке на угол Ѳ(t) = Dr cos (ott, то при вычислении функций корреляции (4.29) необходимо на основе какойлибо модели гармонического осциллятора вычислить сред нее время жизни ротонов. В теории Байера ха не зависит от температуры, тогда как у Уосснера и Гутовского вводится
га = |
Ѵ2ег'та, |
где |
X = |
Ѣ&іІкТ. |
Для случая ѵ?т2 |
1 и |
||
ѵ2Та<^1, результаты |
Байера |
и |
Уосснера — Гутовского |
|||||
можно записать в компактном виде: |
|
|
||||||
|
|
|
А Ы — ха/2Ті + |
— О, |
|
(4.37) |
||
где |
в случае |
Байера |
|
|
|
|
|
|
|
|
|
|
|
|
3 - е _х |
|
|
|
|
|
|
|
|
X |
л |
|
|
|
|
|
|
|
е |
— 1 |
|
|
|
|
3 |
(/іѵо)2 |
Г 4ch X— 3 |
|
|
|
|
|
2 |
32 |
(яJt\\f |
|
(еж — I)2 J |
’ |
|
§ 1] Д В У Х У Р О В Н Е В А Я СИСТЕМА 119
а в случае |
Уосснера — Гутовского |
|
|
|
||||
Аі |
3 |
hvl |
* |
л |
3 |
|
(hvо)* |
1 |
8 |
яV V3 |
2 |
’ |
2 ~ 64 |
(2я.Гѵ()2 |
s h 4 - |
||
|
|
л J r t |
|
|
|
v |
t V |
В уравнение (4.37) в случае Уосснера — Гутовского вмес то та следует подставить т'а. Эти соотношения позволяют определить среднее время жизни квантов вращательных качаний по экспериментальным значениям Тг.
Теория Уосснера — Гутовского во многих случаях луч ше согласуется с экспериментом, чем теория Байера, пред сказывающая менее крутой ход кривых температурной за висимости времен спин-решеточной релаксации. Однако обе эти теории учитывают лишь вращательные качания мо лекул в решетке, поэтому в некоторых температурных ин тервалах уравнение (4.37) может иметь комплексные или чисто мнимые корни.
Ван-Кранендонк и Чанг [17,18] *) учли влияние транс ляционных колебаний решетки на Тх. Если использовать распределение Дебая с характеристической температурой Тв для решеточных колебаний, то температурная зависи мость времени квадрупольной спин-решеточной релакса ции будет следующей:
-^r---- Т2 ( а -----для J, > - i - 7 ,D,
(4.38)
—для 7’<0,027,d,
где а и Ъ — константы. Теория Ван-Кранендонка и Чанга может быть применена для ионных кристаллов.
Иосида и Мория [19] учли влияние решеточных колеба ний на изменение ковалентности, Кондо и Ямашита [20] — влияние перекрывания орбит ионов в ионных кристаллах. Кочелаев рассмотрел роль оптической ветви колебаний решетки в квадрупольной релаксации [21]. Эксперимен тальная проверка различных теорий квадрупольной ре лаксации в ЯКР была выполнена в ряде работ [16, 18, 22-27].
*) Работа Чанга недоступна русскому читателю, поэтому она цитируется по [18], где приведены ее основные результаты.