Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Гречишкин В.С. Ядерные квадрупольные взаимодействия в твердых телах

.pdf
Скачиваний:
26
Добавлен:
24.10.2023
Размер:
10.56 Mб
Скачать

н о

КВАДРУІІОЛЬНАЯ РЕЛАКСАЦИЯ

 

[ГЛ. 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], где приведены ее основные результаты.

Соседние файлы в папке книги из ГПНТБ