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

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

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

130

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

[ГЛ. XV

Рис. 25. Изменение разности населенностей трех переходов в зависимости от величины W it при у = 0,1 для семи начальных условий, приведенных на рис. 23.

§ 2]

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

131

 

 

5*

132

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

[ГЛ . IV

Для произвольного полуцелого спина решение кине­ тических уравнений можно записать в виде

щ (t) — 2 « iÄ)exp(— Х^Ч) + Пі(оо) (к = 1, 2, 3 , ... , J —1/2).

(4.54)

Здесь коэффициенты а™ и щ (оо) зависят от начальных условий, а Х(к) являются собственными значениями опре­ делителя однородной части уравнения.

Экспериментальные данные по спин-решеточной релак­ сации в многоуровневых системах ЯКР были получены в

ряде работ [33,

34, 36—38]. Для измерения времени Тг

 

 

 

 

удобно

в

этом

случае

 

 

 

 

пользоваться

методом

 

 

 

 

стимулированного

эха.

 

 

 

 

Методика обработки эк­

 

 

 

 

спериментальных

дан­

 

 

 

 

ных ясна из рис. 27.

 

 

 

 

В табл. 6 представле­

 

 

 

 

ны результаты

измере­

 

 

 

 

ний при 77 °К для ряда

 

 

 

 

соединений

[34, 38]

(ре­

 

 

 

 

зонанс

ядра Sb123, J =

 

 

 

 

= 7/2). Из таблицы видно,

 

 

 

 

что

релаксационный

Рис. 27.

Изменение амплитуды стимули­

процесс

 

описывается

несколькими

релак­

рованного

эха в 2SbC]3(C«H5)2CH2 в

зави­

 

симости

от t.

 

сационными константа­

чаев эксперимент

позволяет

 

ми. Однако в ряде слу­

найти лишь одну константу.

Конечно, следует

иметь в виду, что для определения всех

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

Поскольку время квадрупольной релаксации выража­ ется через температуру Т, сог — частоту вращательных качаний молекулы, J t — момент инерции молекулы и т„, то в многоуровневых системах также возможно определе­ ние среднего времени жизни ротонов. Отношение y = W J W 1

§

2]

 

МНОГОУРОВНЕВЫЕ

СИСТЕМЫ

133

 

 

 

 

 

 

 

Т а б л и ц а 6

 

Соединение

Переход

Частота

Тш,

ТЬ

Т1С

 

Мгц

мксек

мсек

мсек мсек

 

2SbCls(CeH5)2CH2

V*

3h

35,706

280

30

1,9

 

2SfcCl3-CeH5OCH3

3h

—>6h

68,70

1200

13,7

4,7

 

4 t —*3h

38,521

260

0,8

 

 

 

3/2 —*bh

40,110

320

8,4

 

2SbCl8-CioH8

72,432

700

18

6,8

 

4 t

3h

36,058

250

10

 

SbCla

 

3h —»•It

72,151

2600

10,6

 

 

4 t —»3h

39,096

340

 

 

 

 

3h

&h

68,6

1600

13,7

 

 

SbBrä

bh

7/2

104,5

1520

12,4

 

 

4 t —>3h

31,1

130

3,4

 

 

2SbBr3CioH8

•It —j.*h

59,9

700

4,2

 

 

3h

6h

60,8

700

5,2

 

в этом случае можно

оценить из

отношения релаксацион­

ных констант.

 

 

 

 

 

 

 

Для /

= 7/2

релаксация описывается тремя релакса­

ционными

константами.

Действительно,

N+i/t + N+»/, +

+

^ ± " / 2 +

N±y, — Na, где N ±i — число спинов в данном

квантовом состоянии, N 0 — полное число резонирующих ядер. Таким образом, одно из кинетических уравнений яв­ ляется комбинацией остальных. Производя вычисление ре­

лаксационных

вероятностей для J = 7/2, получим, следуя

предположениям

Байера

[15],

 

 

W l = m

HW2IQ

(в* -

[(2 cha: — 1) + (ех

1 )1.

 

4

J A Xa

Б2

 

(4.55)

 

 

 

 

 

 

 

 

_ 441 * 4

^

1

 

 

 

 

8

J?G>f

(ex— 1)2 X

 

 

X

l+co

2 (ch X— 1)

 

2ch X— 1

 

 

T“ (ex — l)2

l + co;

 

 

 

m+2 ,w

 

 

m+2,mTa_

 

где ©tTa

1,

 

 

 

©Ql — частота ЯКР для пере­

хода + V2 -> + ®/2.

 

 

 

 

Таким образом, хотя W2 и является функцией частоты

перехода ©w+2,m, этой зависимостью при Т >

20 °К мож-

134

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

[ГЛ . IV

НО

пренебречь,

ПОСКОЛЬКУ (0TO:j-2,m ха ^

1 ПРИ высоких

температурах.

С другой стороны, при Т

20 °К в силу

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

Из температурной зависимости частот ЯКР двух пере­ ходов можно определить температурную зависимость кон­ стант квадрупольного взаимодействия eQqzz. В случае комплексов SbCl3 [38] в предположении слабой связи меж­ ду компонентами [39] можно определить среднюю частоту вращательных качаний молекулы SbGl3; момент инерции этой молекулы J t = 501-10-40 г-см2 (вычислено на основе

рентгеноструктурных

данных).

 

 

 

Действительно [401.

 

 

 

 

eQ?2Z—

167,4 ■ К Г 11

 

0,6Jtii’

(4.56)

e

— 1Г1

«Q4ZZ

V i

 

 

 

(здесь vt в см_1, J t в г-см?). Формула (4.56) получена в предположении, что пирамида SbCl3 в комплексе соверша­ ет вращательные качания вокруг двух главных осей тен­ зора инерции X и Y. Качания вокруг оси симметрии треть его порядка, проходящей через атом сурьмы, неактивны в спектре ЯКР. Поскольку частоты вращательных качаний vt — \ѵх + Ѵу)/2 при 77 °К не были измерены из спектров комбинационного рассеяния света, то приходится восполь­ зоваться оценками из данных ЯКР. Небольшие искажения формы пирамиды SbCl3 в комплексах не приводят к суще­ ственному изменению J t = (J х + /у)/2. Результаты рас­ четов при 77 °К приведены в табл. 7. В этой таблице gt — коэффициент в уравнении

со, = со? (1 - gt АТ),

(4.57)

где со°і — частота вращательных качаний молекулы при определенной температуре Т0, АТ — Т Т0.

Знание коэффициентов gt позволяет производить экст­ раполяцию частот vt на другие температуры. Сравнение с

§ 2}

М Н О ГО У РО В Н Е В Ы Е

СИСТЕМ Ы

135

 

 

 

Т а б л и ц а 7

 

С о е д и н е н и е

 

g {, г р ад *1

Y

 

 

 

 

2 SbCb(CeH6)2CH2

53

0,00137

4

 

2SbCls-CeH5OCH3

37

0,00058

2

 

2 SbCle-CioHe

32

5

данными работы [41] показывает, что полученные значения vt совпадают с результатами спектров комбинационного

рассеяния света

при

290 °К.

 

Из уравнений (4.36) для ravt )> 1 и VQ1T(1 <[ 1, получаем

г

=

2ех+ Г * - 3

(4.58)

Т°

 

Аѵ( 2ех+ 2е~х— 3

 

 

Для комплекса 2 SbCl3-CeH 5OCH3 время т„ = 0,64-10-11 сек при 77 °К; для 2 SbCl3*C10H8 та = 0,95'К Н 1 сек; и, наконец, для 2 SbCl3-(CeH 5)2CH2 ха = 0,77-10-11 сек.

Полученные значения ха примерно совпадают со сред­ ними временами жизни фононов, а также не противоречат данным о ширинах линий в спектрах комбинационного рас­ сеяния света малых частот. Время жизни молекулы на вращательно-колебательном уровне, таким образом, боль­ ше периода вращательных качаний, т. е. т0 )> 1/ѵ(.

И. Пригожин рассмотрел вопрос о приближении к рав­ новесному состоянию одиночного нормального колебания при условии, что остальные степени свободы уже достигли равновесия [6]. На основании решения своего основного кинетического уравнения для ангармонической модели твердого тела Пригожин показал, что время релаксации tr по порядку величины равно Ит2СкТ, где со — частота

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

136

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

[ГЛ . IV

Следует отметить, что сами релаксационные константы Тц в многоуровневой системе могут быть функцией радио­ частотной мощности [42], действующей на соседний пере­ ход. Если насыщается переход на частоте ѵ2 или ѵх, то в ки­ нетические уравнения необходимо ввести вероятность пере­ ходов Р под действием радиочастотного поля. Для J = 8/2 при насыщении перехода на частоте ѵ2 можно получить для релаксационных констант

у(1.2)

=

41 С1 + Ч) + Р ±

У -5(7 6 - 88Т +

61г 2) +

+ І

^

~ 7Г) + Р2У \

P - P I W u г

(4.59)

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

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

§ 3. Поперечная релаксация в ЯКР

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

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

Следует отметить, что при отсутствии внутренних дви­ жений в решетке, обычно поперечные компоненты намагни­ чения затухают не по экспоненте. Неэкспоненциальное за­ тухание естественно не позволяет строго ввести время Т2.

*)

Подробнее этот вопрос обсуждается в работе В. С. Г р е ч и ш ­

кин,

Е. М. Шишк и н, Изв. вузов, физика, № 3, 82 (1973).

§ з]

П О П Е Р Е Ч Н А Я РЕЛ А К С А Ц И Я В Я К Р

137

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

Рассмотрим две группы спинов А и В [43]. Спины А принадлежат к квадрупольной системе и дают сигнал квадрупольного спинового эха, а спины В создают в точке / расположения спина Aj (это будем отмечать индексом/ у А : Aj) флуктуирующее локальное поле. В системе спинов В, например, может развиваться какое-либо движение. Конечно, взаимодействие типа AjAft тоже существенно. Однако предположим, что спины В имеют большой маг­ нитный момент; тогда их вклад в общее диполь-дипольное взаимодействие станет преобладающим:

Ж = ГлГвОТ« 3 ( 4 ^ -

аГл1°;4А‘и{м1

I . (4.60)

1 = 1 І/ А І В І

r A i B j

J

где rAtBj — расстояние между і-м А-спином и /-м В-спи- ном; уд, YB — гиромагнитные отношения ядер. За счет этих взаимодействий линия ЯКР уширяется на величину

(АЯр)г =

(4.61)

Та

где б; (t) = ГаТв^ 2 msj(t) GU,

3=1

1— 3 cos20AiBj

 

Gij =

(4.62)

 

'AiBj

Ѳаіві — угол между гдів}и осью Ъ градиента электричес­ кого поля на ядрах А, тв} — магнитное квантовое число /-го спина В. Пусть суммарный магнитный момент спинов В подчиняется уравнению

dM,

S T - = - RM bp,

(4.63)

тогда Мър (t) = M 5 „ (0) e~Bt, где параметр

R характери­

зует скорость скачков локального поля. Отсюда среднее

значение б =

8пе~ш (здесь R — 1/Г2).

найти вероятность

Для расчета сигналов эхо нужно

Р (б, t, б0),

которая дает частотное

расщепление б на

138

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

[ГЛ . IV

спинах А в момент времени t, если при t = 0 оно имело значение б0. Функция Р (б, t, б0) должна удовлетворять условию

оо

5 P(80)P(ö,t,Ö0)db0 = P(ö),

(4.64)

—00

 

где Р (6) есть стационарное гауссово распределение рас­ щеплений. Далее

оо

б = 5 öP(8,t,ö0)dö = 80e~Rt.

(4.65)

—оо

 

Условие (4.65) устанавливает связь между функцией Р (б, t, б0) и поведением вектора Mgp. Наконец,

ОО

 

5 P(ö,t,ö0)dö=: 1.

(4.66)

—о

 

Можно показать, что всем этим условиям удовлетворяет функция

P(Ö, і,б0) = J_______ 1_______\

f

(5— 6ne~Rt)8 )

бА [2л (1 -

е-іШ)]‘,гJ еХР 1

26а (1 - e-*Rt)]'

 

 

 

(4.67)

К подобному распределению, например, приводит урав­

нение Ланжевена

 

 

 

Ж +

£ ( , =

АЩ,

(4.68)

где и — скорость молекул газа, т — масса молекул, %і — коэффициент трения, А (t) — флуктуирующая сила, воз­ никающая за счет столкновений.

Пусть 90°-импульс поворачивает вектор М 0 на 90°. Тогда каждая компонента момента М0Р(б„) будет в фазе при

t =

tw на частоте прецессии б0. Весь интервал времени

t

tw разобьем на к равных интервалов малой длитель­

ности At< ^ 1/7?, где kAt = t. За

время At функция б

изменяется медленно, а функция А

(t) успевает сильно из­

мениться. Функцию б найдем из

уравнения

 

 

(4.69)

§ 3] П О П Е Р Е Ч Н А Я РЕЛ А К С А Ц И Я В Я К Р 139

где F (t) — флуктуирующая функция, смысл которой ана­ логичен А (t).

В первый интервал

At

вероятность того, что спины А

образуют момент | F0

| =

Р (80) М 0 на новой частоте 6Х

дается функцией Р (б,

At, 80) из уравнения (4.64). В пе­

риод At средняя частота прецессии F0будет - (б0 + бг)/2, однако, в результате ряда быстрых изменений частоты за счет флуктуаций в конечном результате мы получим не­ большую ошибку, предполагая, что 60 — средняя частота. Часть каждого возможного момента F0дает вклад в новый момент I Fx I = Р (öj) М оф, имеющий частоту бх при t — At, где функция ер <С і включает затухание и фазовый мно­ житель, который дает Fx < F0 при том же самом положе­ нии спектра бх = б0. Это является следствием фазовых флуктуаций в системе ядер, возникающих из-за их пере­ мещений от области начальных частот 80 до конечных час­ тот бх. Поэтому

со

= іМ„ $ Р (б0) Р (б1; At, б0) a««*d60.

(4.70)

—оо

 

Вследующий момент времени At между t = At и t =

=2 At распределение моментов Fyдает вклад в новый мо­

мент Fi на частоте б2 и т. д. Тогда

О О

FH= $

JW > (6*. At,

(4.71)

—оо

 

Сигнал индукции

Рѵтя после 90°-импульсов будет пропор­

ционален сумме

прецессирующих

моментов Fk:

Т’инд - ] Fkdök =

М 0ехр { _ ( ^ - ) 2

(«Н» + R t - 1)}. (4.72)

Аналогичные рассуждения позволяют получить сигнал эха:

Е — М 0ехр ( “S')

1 +

е~ш +

R (t — т)

j .

 

 

 

 

 

 

(4.73)

При

t —

амплитуда эха

достигает максимума:

Ятах =

м йехр {

-

[Дт -

1 + (1 +

Rx) <г*дч} .

(4.74)

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