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

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

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

120

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

[ГЛ . ГѴ

В работах [24—27] теория квадрупольной релаксации для J = 3/2 была использована для оценок средних вре­ мен жизни ротонов т0 и х'а в кристаллах. Несложные расче­ ты показывают, что при достаточно высоких температурах основной вклад в квадрупольную релаксацию вносят двух­ квантовые переходы (Дт = + 2). При 0 )QTO ]> 1 темпера­ турное изменение Тг для двухквантовых переходов опре­ деляется множителем

(ех — I)2 L

2 (ch X

 

} + 2сЪх

-1

(4.39)

-

1)

1J

Если X 1, то

 

1

 

 

1

 

т

• ®

-,

Тг

(4.40)

|Д т |= ± 2

j i j

 

|Д т |= ± і

Г '

Поэтому из температурной зависимости

можно сделать

вывод о роли двухквантовых и одноквантовых переходов в релаксации. Тогда на основании экспериментальных дан­ ных по Тг можно определить та.

Оценки та получаются более надежными, если известен спектр комбинационного рассеяния малых частот и струк­ тура молекулы. В некоторых случаях в силу изотропной поляризуемости вращательные качания неактивны в спект­ ре комбинационного рассеяния. В этом случае для оценок V, можно воспользоваться температурной зависимостью частот ЯКР. Моменты же инерции молекул можно рассчи­ тывать на основе рентгеноструктурных данных. Вклад в ширину линий спектра комбинационного рассеяния малых частот за счет конечного времени жизни определим по фор­ муле

Дѵа = 1/агта,

(4.41)

где а 1 — некоторая постоянная, равная л для лоренцева контура линии.

Экспериментальные данные по средней ширине линий комбинационного рассеяния света малых частот взяты из работ [28, 29] и приведены в табл. 5.

Из формулы (4.37) получается два значения та. Вопрос о том, какое из этих значений больше соответствует дейст­ вительности, требует особого рассмотрения. Тепловое воз­ буждение кристалла можно описать в виде волн, распро­ страняющихся в кристалле со скоростью звука и обладаю-

§ і]

д в у х у р о в н е в а я

с и с т е м а

 

121

 

 

 

 

 

Т а б л и ц а 5

 

Т 1, сек

1

v(

 

 

Дѵ а ,

Дѵ,

С о ед ин ен ие

ѵ„ dr

г-см2

V сеК

 

3U0 °К

C M * 1

см*1

САІ“ 1

 

 

град-i

 

 

 

 

 

л-СвШВгг

0,0004

1,1

40

2208

10-и

0,3

6

п-СвШСЬ

0,02

0,9

47

1309

10-13

0,3

5

NaClOs

0,04

2,25

107

87

5.10-1’

4

10

КСЮз

0,03

2,35

102

89

5-10-1’

4

10

Ва (СЮа).

0,058

0,9

150

90

ю -i’

2

Са (С10з)а

0,049

0,9

150

„ 89

і о -i’

2

Sr (СЮзЬ

0.051

0,9

150

89

ю -i’

2

Mg (С10з)з

0,048

0,9

150

89

і о - ’2

2

щих способностью переносить энергию. В гармоническом приближении эти волны перемещаются свободно без зату­ хания и, следовательно, имеют неограниченную длину пробега. На самом же деле их пробег ограничен рядом при­ чин: а) энгармонизмом решетки, б) загрязнениями и нару­ шениями правильности структуры кристалла, в) конечны­ ми размерами кристалла.

Вопрос о среднем времени жизни квантов трансляцион­ ных колебаний (фононов) был рассмотрен Пайерлсом [30]. Им была получена формула для коэффициента теплопро­

водности твердого тела:

 

V = уСск,

(4.42)

где С — теплоемкость на единицу объема, с — скорость звука в кристалле, к — длина свободного пробега фонона, у — численный множитель, равный примерно Ѵ3. При по­ нижении температуры теплопроводность быстро возраста­ ет, в отличие от теплоемкости, которая падает. Свободный пробег фононов становится очень большим. Обычно тепло­ вое сопротивление определяется столкновением между фононами. Пайерлс показал, что вероятность таких столкновений падает экспоненциально при понижении температуры, что вызывает экспоненциальный рост тепло­ проводности.

Вопрос же о средних временах жизни ротонов более сложен. Однако в силу того, что в обоих случаях ангармо­

122 К В А Д Р У П О Л Ь Н А Я РЕЛ А К С А Ц И Я IV

низм обусловлен одними и теми же причинами, различие между та и средним временем пробега фонона не должно быть велико. Это позволяет сделать необходимый выбор из двух значений ха.

Если учесть в простейшем приближении нелинейные

силы в решетке, то можно написать

 

тк

— а$к

(4.43)

где mh — масса молекулы, S h — смещение молекулы из положения равновесия, а и ß — упругие постоянные. Вве­

дем безразмерное время т = У ahnt.

Тогда

S\ + S k = —L sl

(4.44)

Решая нелинейное уравнение (4.44) методом медленно ме­ няющихся амплитуд, получим

ха— л Y тк/а .

(4.45)

Оценка ха по формуле (4.45) для хлората калия дает вели­ чину порядка ІО-12 сек.

Из табл. 5 видно, что вклад в ширину линий комбина­ ционного рассеяния света составляет 5—40% от наблюдае­ мой на опыте величины. Интересно отметить, что ха в па­ радихлорбензоле и парадибромбензоле одинаковы, что указывает на сходное строение решеток этих соединений. Эти кристаллы, как известно, изоморфны. Причем оценка ха из данных по Тх для изотопов Br79, Вг81 и Cl35, Cl37 дает одинаковую величину. Для Вг79 значение Тхпри комнатной температуре в га-СвН4Вг2 равно 280 мксек. Такое сокра­ щение Тг по сравнению с Т\ для Вг81 связано с тем, что частота ЯКР Вг79 при 300 °К равна 268 Мгц, а не 224 Мгц, как дляВг81. Подобные оценки ха неизбежно приближены,

поскольку считается,

что

 

В работе [27] сделана детальная проверка теорий Байе­

ра и Уоснера — Гутовского

без использования предпо­

ложения, что

1. В

ряде комплексов на основе

AsCl3 и AsBr3 время ха экспоненциально зависит от темпе­ ратуры, причем W2 ]> Wv Теория Байера лишь в отдель­ ных случаях позволила объяснить экспериментальные ре­ зультаты. Вебер [23] для J = 3/.2 предложил способ оце­ нок ха из экспериментальных данных ЯКР, если измерено

§ 2]

М Н О Г О У Р О В Н Е В Ы Е СИСТЕМЫ

123

отношение у = W J W lt где W2 и W l — релаксационные вероятности. Однако метод Вебера требует применения мо­ нокристаллов, а также не обеспечивает достаточной точно­ сти, поскольку измерения проводятся в слабом магнитном поле для различных компонент мультиплета.

Анализ экспериментальных данных в ЯКР, особенно при низких температурах, мало что может сказать о кон­ кретных деталях модели решетки. Квадрупольная релак­ сация является интегральной характеристикой решетки, поэтому она мало чувствительна к различным моделям.

§ 2. Квадрупольная релаксация в многоуровневых системах

Если для двухуровневой системы можно ввести поня­ тие об одной спиновой температуре Ts, то для многоуров­ невой придется ввести несколько парциальных спиновых температур, характеризующих разности заселенности соот­ ветствующих пар уровней. Многоуровневую систему в ЯКР можно получить, наложив на монокристалл внешнее магнитное поле (для J = 3/2). Если J 3/2 (например, в/2, 7/2, 9/2), то за счет только чисто квадрупольного взаи­ модействия возникает система неэквидистантных уровней энергии. Неэквидистантная многоуровневая система будет и в случае J = 1, ц =j=0. Во всех этих случаях нельзя огра­ ничиться введением одной релаксационной константы для характеристики процесса спин-решеточной релаксации.

При исследовании спин-решеточной релаксации в мно­ гоуровневых системах можно рассмотреть случаи одно­ частотного и многочастотного возбуждения спин-системы. Полуцелые спины имеют в ЯКР особенно большое значе­ ние, так как основной экспериментальный материал полу­

чен для галогенов.

Гд) основное влия­

При высоких температурах

ние на температурную зависимость времени квадрупольной спин-решеточной релаксации в лабораторной системе координат, как уже отмечалось выше, оказывают враща­ тельные качания молекул. По мере увеличения амплитуд вращательных качаний молекул в решетке кристалла возрастает плотность спектра «молекулярного шума» на частоте резонанса, что приводит к интенсивному восста­ новлению теплового равновесия в системе, нарушенного в

124

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

[ГЛ. IV

результате поглощения квантов электромагнитного поля. Если спин ядра полуцелый и больше 3/2, то число квадрупольных уровней для J = 5/2; 7/2 и 9/2 равно соответствен­ но 3, 4 и 5, причем уровни дважды вырождены.

Обобщение теории Байера [15] на многоуровневые сис­ темы было дано в работах [31—34].

Если не учитывать связи между нижней и верхней па­ рой уровней, мера которой дается релаксационными ве­ роятностями, то для / = 5/2 можно получить [31]

т

, ,

. 64

 

х

 

 

 

X

2 (ch х — 1)

2ch X — 1

1

 

 

 

i)*

+

 

 

 

L 1 + (üjr* (e

'

(4.46)

 

 

я

(яѵ./.)2

 

 

 

 

 

 

 

T1 , м = 4 т а г г ( e * - l )2 X

 

| Д т | = 2

 

s

(n v Q> T a

 

 

 

 

 

X

2 (ch X — 1)

2ch X — 1

 

 

 

Ь W2Ta (e

—1)

+

J

 

 

1

1 + « K

Таким образом, спин-решеточная релаксация описывается двумя релаксационными константами, связь между кото­ рыми может быть установлена из решения кинетических

уравнений [33].

Т2 для каж­

Рассмотрим случай / = 7/2. Пусть Тх

дого из возможных переходов. Тогда уравнение (4.26) может быть приведено к виду [1] *)

d N

(4.47)

S T - = 2 (WnmNn - WmnNm),

П

 

где W nm, Wmn — вероятности прямого и обратного пере­ ходов между уровнями п и иг; N n и — населенности соответствующих уровней. Применимость этого кинети­ ческого уравнения для различных систем, обладающих диссипативными свойствами, исследовалась в общей тео­ рии необратимых процессов.

Будем считать, что диссипативная система (решетка) находится в равновесии, а динамическая (спины) выведена

*) Возможность использования этого уравнения, в частности, обусловлена тем, что в твердых телах обычно Тх > Тя.

§ 2]

М Н О Г О У Р О В Н Е В Ы Е СИСТЕМЫ

125

из равновесного состояния. Для сшш-сшшовой релакса­ ции в твердых телах такой подход, как правило, некоррек­ тен [35], поскольку система спинов не может находиться в равновесном состоянии. Как и раньше недиагональные элементы матрицы плотности учитывать не будем.

Для чисто квадрупольного резонанса нужно рассмат­ ривать две вероятности релаксационных переходов [1]:

W { m ± \ , m ) =

(2^ 1};^ 1(у т+1) W* (4-48)

W (т ± 2 , т ) =

(І+ т ) (JТ - ■-■/2)^ Т±п!;Ѣ (J ±m + 2) W t,

где m — магнитное квантовое число, / — спин ядра, Wt — вероятность одноквантового механизма релаксации, W2 — вероятность двухквантового механизма релаксации.

♦ b

л

 

 

~ 2

 

 

 

fa m

2/Wz

jj(Mâ)W2

jf{h5â)WB

7 rv, fymw,

'' 5

 

 

 

 

+ f -

 

 

 

2

 

 

 

—W

§Lw §m)w,

 

 

 

 

 

 

147п/

. 3 , Pt

Тд( МЩ

’l —Wv

Ж щ ^(/*3d)W2

,, 49WZ

■ - І

é : mm,

Ж,

 

 

M.W 1

,, /47

щ ( Ш г

и

w ’

♦ Г

 

 

2

Рис. 21. Вероятности прямых и обратных переходов между чисто квадрупольными уровнями в случае J = у 2.

Картина всех возможных релаксационных переходов представлена для J = Ѵ2 на рис. 21. Отсюда получим для

т} = О

 

 

dN.,

 

 

± Ѵ г

~ W 1(1 + ЗА) + J L W2(1 + 5A)

+

dt

+ — W1N±V2 + - ± - W ,N ±Vt,

126 К В А Д Р У П О Л Ь Н А Я РЕЛ А К С А Ц И Я [ГЛ. ГѴ

dN.± Ч г

32

Wl(l +

 

2A) + -g - W 2(l + 3A) +

dt

147

 

 

32

 

 

 

 

 

 

 

 

 

^ i] w ±./2+

 

WiN.

 

10 WtN±Vt +

 

 

 

 

147

 

 

±SA

49

 

 

 

 

 

 

 

 

 

 

 

+ 4 - ^

х(і +

за) іѵ±,/2,

d N ,

10

 

 

 

 

 

 

 

 

14740 W2(1+A) +

-d( _

І47. ^ г(1 + А)

 

21

W,

3 2 -Wl] N d

'

10

■WXN,

 

 

32

W 1(i+ 2 A )N ±Vt +

±•4

_

147

'± s/2

 

147

 

 

 

147

 

 

 

 

 

 

 

2

W 2{i

 

5Д) N±7/t +

W2N цз/г,

 

 

 

21

 

 

 

 

 

 

 

 

dN.±V.

10

Wt

10

 

 

40

(1 +

A)PF2 ^±V.+

dt

147

49 W2

147

+

10

+

A)

 

 

4 -

10

W2(l +

3A)iV±./2 +

147

 

 

49

 

 

 

 

 

 

 

 

40

- ^ 2(1 +

А)ІѴт а , (4.49).

 

 

 

 

 

 

 

 

ш

где А =

hvQ/кТ, VQ— частота ЯКР для перехода

Ѵ2 —>■3/2,

Т— абсолютная температура. Введем обозначения

п±і = N ±iTi/2— N ±і±у„ п+і = п-г, N+i = 7Ѵ_;.

(4.50)

Тогда из четырех уравнений для заселенностей (4.49) по­ лучаем три уравнения для разностей заселенностей уров­ ней:

*8 +

-Я - (6И7! +

w 2)

П

з - ^

(2W1+

W2) п2-

 

 

---- W W.«i =

-п г

( 4 7 ^

-

5 ^ 2) по,

** -

 

-

^ * )и» +

і |т - ( 16^ і

+

1 W 2)«2 -

 

- ^

( ^ i

+ W 2K

= ^

( - W

i

+

15W2)n0, (4.51)

% -

4 -

-

- щ - (2Wi -

Ws) n2

+

 

(2Я\ +

 

 

+

 

 

 

l W i + 2ЪШ2) щ,

где

«з =

AT./, — Ny„

п„ = Ny, Ny„

rax = ATV, — Ny,.

§ 2]

 

М Н О Г О У Р О В Н Е В Ы Е СИСТЕМЫ

127

п0 =

ДІѴ/4,

Ni = N14, N — полное

число резонирую­

щих

ядер.

 

 

W2/W 1и Wttl147 =

Если теперь ввести обозначения у =

— X (безразмерное время), то уравнения (4.51) перепишут­

ся в

виде

 

 

 

 

п3 +

14 (6 +

у) п3 — 16 (2 +

у) щ — 30 упх =

 

 

 

 

 

= 4 (47 — 5 у) п0,

?г2 — 14 (3 — у) п3 -j- 4 (16 +

11 у) п2 — 10 (1 +

у) щ =

 

 

 

= 8 (—1 +

15 у) п0,

п1 — 14 уп3— 16 (2 — у) n2 -f- 10 (2 +

11 у) /гх =

 

 

 

 

= 4 (-1 1 + 25 у) п0. (4.52)

Для решения этой системы применим операционный метод Лапласа. Тогда характеристическое уравнение по­ лучим в виде

14(6 +

у)

— 16 (2 + у)

— ЗОу

 

— 14(3 — у)

р + 4(16 + 11у)

— 10(1-)-у) =

0.

— 14у

 

— 16(2 — у)

 

р + 10 (2 + 11у)

 

 

 

 

 

(4.53)

Корни этого определителя рі =

147/T liW1могут быть най­

дены численно

для различных

значений параметра

у

Рис. 22. Зависимость отношения времен спин-решеточной

м

/

релаксации у /Т

 

от Y = ИѴѴГ,.

1

1

Отсюда T1i—i41lpiW l. Попарное отношение корней урав­ нения (4.53) приведено на рис. 22. Из рисунка видно, что

128 К В А Д р У П О Л ЬН А Я р е л а к с а ц и я [ГЛ . XV

процесс релаксации носит сложный характер. Изменение разности заселенностей не может быть определено лишь од­ ной экспонентой. Лишь при у — 1 и 7 ^ 4 - ь 8 имеет место вырождение, когда два корня уравнения совпадают друг с другом.

Рассмотрим теперь влияние различных начальных ус­ ловий (рис. 23). В случае 1 насыщающее радиочастотное

/

1

 

п 3

З П д

пз

4п0

п 3

 

пг

2,5п0

П г ' 10

Пг

 

л,

t 0

Л и г п а

п,Г

 

4

 

1

 

t

П3

Щ п д

Пд

' 0

п

' 0

п

Л о

п г

['0

пг

Ь п 0

П/ * 0

лг

Щл0

JhA о

‘ 0

ijfl

Л д

7

п3 ' о

пг ' 0

П; і о

Рис. 23. Различные случаи начальных условий для разностей населенностей уровней энергии для J = 7/г.

поле воздействует лишь на переход + Ѵг —»• + 3/2, а в случае 4 насыщаются в начальный момент времени пере­ ходы + Ѵ2 -* + 3/2 и 3/2 —> + 5/2. Затем радиочастот­ ные поля быстро выключаются, и система спинов приходит в равновесие с решеткой.

На рис. 24, 25, 26 представлены решения уравнений (4.52) для Y = 2, 0,1 и 5.

Из этих рисунков видно, что случаи 1 и 4 сильно отли­ чаются друг от друга при у — 0,1, а при у = 2 отличие менее заметно. В некоторых случаях процесс приближения к равновесию носит немонотонный характер. Таким об­ разом, налицо зависимость в поведении вектора ядерного намагничения от характера радиочастотного возму­ щения.

Очевидно, что изменение начальных условий не влияет на величину релаксационных констант, но существенно изменяет величину ісоэффициентов при экспонентах.

§ 2]

М Н О ГО У РО В Н Е В Ы Е СИСТЕМ Ы

129

5 В. С. Гречишкин

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