
книги из ГПНТБ / Гречишкин В.С. Ядерные квадрупольные взаимодействия в твердых телах
.pdf120 |
К В А Д Р У П О Л Ь Н А Я РЕЛ А К С А Ц И Я |
[ГЛ . ГѴ |
В работах [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 |
|
х |
1у |
|
|
|
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 В. С. Гречишкин