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

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

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

30 УРОВНИ ЭНЕРГИИ И ЧАСТОТЫ Я КР [ГЛ. I

Тогда

где

 

Фі (^i)

= Сіфі +

 

г2фа,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с,

. -

 

 

 

q

іал2 + г)2 ’

 

 

 

У 2

(Д2 + д2) + 2 л

 

 

_________ Д 4- /

Д2 +

If_______.

(1.61)

 

 

У 2 (Л2 +

ті2) + 2Л "КЛ2 +

т)2

 

 

 

где

 

Фг (Х2) =

Сіфі +

с2ф2,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

 

 

 

 

Л

 

 

 

 

 

 

1

 

У 2 (Л2 +

т)2) — 2Л / Л

2 +

г)2

 

'

 

_______Д — V Л2 +

д2______ .

(1.62)

 

 

У 2 (Л2

 

г)2) — 2Л /

Л2 +

г|2

 

 

 

 

 

 

 

Фз (*з)

= Фз-

 

 

 

 

 

Если ось радиочастотной катушки направлена вдоль X,

то Ѳі = я/2, фі = 0.

 

 

 

 

 

 

 

 

 

Рассмотрим переход с уровня Х3 на уровень А*:

 

 

 

 

/ 0

 

 

0

 

 

1 / / 2 \ / с л

 

Рзд = (0, 0, 1) I

0

 

0

 

 

І / У2

J U J . (1.63)

Отсюда

 

 

Ч ' / У 2

1/|Л2

0

/ Ѵо /

 

 

р?. -= -L

 

 

і

 

 

 

 

 

 

 

- г

 

 

г|2 '

 

 

 

 

 

3,1

2

 

~1~

2 "К Л2 +

 

 

 

При #х I X возбуждаются

переходы ѵ3)1 и ѵ3)2. Аналогич­

но при # і I Y ЯКР

возникает на частотах v3jl и ѵ3>2. При

Нх I Z возбуждается

«запрещенный»

переход ѵ2)1, вероят­

ность которого

равна т}2/ (R 2 + ц2).

 

 

 

 

Аналогичным

способом

можно

 

рассмотреть возбужде­

ние линий ЯКР с помощью радиочастотного поля при про­ извольных ядерном спине и ориентации вектора Н 0 в кри­ сталле [19]. Отметим, что «запрещенные» переходы можно наблюдать, как правило, в параллельных полях (Нх || Н 0). Если сделать Н 0 достаточно большим, то характер смеши­ вания волновых функций изменится, и отдельные переходы могут не наблюдаться.

§ 4]

ТЕМ ПЕРАТУРНАЯ ЗАВИСИМОСТЬ ЧАСТОТ ЯКР

31

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

§ 4. Температурная зависимость частот ЯКР

Зависимость частот ЯКР от температуры была обнару­ жена еще в первых работах Демельта и Крюгера [23]. По мере понижения температуры в большинстве соединений частота ЯКР возрастает. Температурные коэффициенты частоты для различных соединений неодинаковы; по по-

иядку величины —- ^ — ІО'4 К-1.

Демельт и Крюгер объяснили температурную зави­ симость частот ЯКР влиянием решеточных колебаний. Если кристалл построен из молекул, то естественно пред­ положить, что колебания молекул как целого будут про­ исходить на более низких частотах по сравнению с часто­ тами внутримолекулярных колебаний, поскольку силы сцепления между молекулами в решетке слабее, чем между атомами в самой молекуле. Все решеточные колебания при­ нято подразделять на вращательные качания и транс­ ляционные колебания. В свою очередь трансляционные колебания подразделяются на акустические и оптические ветви колебаний. При вращательных качаниях молекула испытывает в решетке малые колебания вокруг неподвиж­ ного центра тяжести. Частоты этих качаний в молекулярных кристаллах достигают ІО12 гц, что значительно превы­ шает частоты ЯКР (ІО7 ~ ІО8 гц). Поскольку с увеличе­ нием температуры амплитуда вращательных качаний воз­ растает, то среднее значение градиента электрического поля (ГЭП) должно уменьшаться. Именно вращательные качания молекул в решетке оказывают наибольшее влия­ ние на температурную зависимость частот ЯКР, что впер­ вые было показано Байером [24].

Поскольку в предыдущих параграфах вращательные качания молекул в решетке не учитывались, то необходимо найти вид qés при учете качаний; именно q'iS и определяют

32 УРОВНИ ЭНЕРГИИ И ЧАСТОТЫ Я К Р [ГЛ. I

теперь наблюдаемые частоты ЯКР. При повороте системы координат тензор ГЭП преобразуется по закону

з

Qas —

2

 

і.te=l

где г - > 1 , у 2, z 3,

а

a si, ask — направляющие ко­

синусы подвижной системы координат относительно не­ подвижной.

Если оси неподвижной системы координат совпадают с главными осями тензора ГЭП, то тензор ГЭП имеет диаго­ нальный вид. При переходе к штрихованной системе коор­

динат (рис. 5) тензор преобразуется по

закону

 

= Яхх<£’г + ѴѵуаЬѵ + Qzz^z'z,

(1.64)

или

 

 

Qz-z’= Qzzcos2 ѳ + Qxx sin2 Ѳ= qzz (l —

sin2 б).

(1.65)

В данном случае вращательные качания молекулы проис­ ходят вокруг одной оси, перпендикулярной к направлению химической связи (направлению

оси Z).

Байер при построении своей тео­ рии сделал следующие допущения:

1)Тензор qtj аксиально сим­ метричен, т] = 0.

2)В процессе колебаний мак­ симальное значение градиента электрического поля qzz не изме­ няется.

3)Учитывается лишь один вид колебаний.

 

Тогда

частота ЯКР

опреде­

Рис. 5. Колебания молеку­

ляется из

соотношения

 

лы в решетке вокруг одной

 

 

 

оси.

ѵ = ѵ0| і ---- !~sm2e j,

(1.66)

 

где ѵ0 — частота ЯКР для неподвижной молекулы. По­ скольку углы вращательных качаний малы, то

ѵ = ѵ „ ( і - 4 ѳ > ) .

( 1 .6 7 )

S 4]

ТЕМ ПЕРАТУРНАЯ ЗАВИСИМОСТЬ ЧАСТОТ ЯКР

33

Учитывая, что Ѳ = Ѳ^іп co^, где <Dj — частота враща­ тельных качаний молекулы, находим после усреднения

= - Г . ѵ = ѵ 0( і — J.0*).

(1.68)

Значение Ѳ? можно вычислить из энергетических сообра­ жений, приравнивая энергию вращательных качаний средней энергии гармонического осциллятора:

 

hvx 1

+ ,

(1.69)

 

 

exp (hvi/кТ) ■ 1

 

где / j — момент

инерции

молекулы относительно

осей

качаний. Тогда

 

 

 

V - Ѵ о [ 4

2 J i(02 (2 + e x p ( A v i y W ) - l ) ] -

^1-70)

Из формулы (1.70) видно, что частота ЯКР с повышением

/Тч

3

hv1

температуры уменьшается. Фактор — - - —; есть вклад

 

4

J лсо:2

 

 

іші

нулевых колебаний в частоту ядерного квадрупольного резонанса. При достаточно высоких температурах этот член мал и им можно пренебречь, а экспоненту exp (h vJkT)

можно разложить в ряд до линейных по h \J k T

членов.

Тогда

 

’ - ' ( ‘ “ Ь д г ) -

<1Л1>

Таким образом, в высокотемпературном приближении ча­ стота ЯКР линейно зависит от температуры.

Формула (1.71), конечно, не может претендовать на большую точность при сравнении с экспериментальными данными. Ее основное достоинство в том, что она качествен­ но правильно описывает температурную зависимость ча­ стот ЯКР. В реальных же кристаллах на температурную зависимость могут оказывать влияние и другие виды ко­ лебаний кристаллической решетки.

В 1955 г. Ф. И. Скрипов учел влияние трансляцион­ ных колебаний решетки на температурную зависимость частот ЯКР [25] *). Частотная плотность поперечных2*

•) Здесь воспроизводится расчет, проделанный Ф. И. Скриповым.

2 В. С, Гречишкин

34

УРОВНИ ЭНЕРГИИ

И ЧАСТОТЫ Я К Р

[ГЛ. X

акустических

колебаний для

изотропного твердого тела

определяется из соотношения

 

 

 

а dz,

V'

(1.72)

 

~WT\д =

 

 

 

 

где Z( — число нормальных колебаний в интервале частот от ѵ( до v t + dvt, ct — скорость распространения попе­ речных акустических волн, W — объем кристалла, ѵ( — частоты поперечных упругих волн. В упругой среде име­ ется две поперечные волны с одной и той же частотой v t. В изотропной упругой среде тепловые возмущения кри­ сталла образуют систему стоячих волн, распространяю­ щихся со скоростью звука. Уравнение такой волны за­ пишем в виде

у = А sin 2я [vtt ---- j ,

(1.73)

у — Л2яѵг cos 2я (vtt ---- 2-

J/max — A2nvt.

Вычислим плотность кинетической энергии в волне:

- = 2я*ѵ?і4*рУ;,

(1.74)

где р — плотность кристалла. Для нахождения функции А (Vj) приравняем плотность энергии в волне средней энергии квантовомеханического осциллятора:

2ЛМ *рГ = *ѵ1(

4 - + н д а ^

п ) .

(1.75)

Отсюда

 

 

 

А2(V,) = 2п*ѵ{рѴ

( іГ + exp (hvt/kT) — l) '

(4-76)

Вычислим теперь угол поворота элемента длины, ориенти­ рованного вдоль оси X, при распространении акустиче­ ской волны в кристалле:

ѳ=

<1Л7>

Тогда амплитудное значение Ѳ равно j- А.

§ 4]

ТЕМ ПЕРАТУРНАЯ ЗАВИСИМОСТЬ ЧАСТОТ Я КР

35

Для вычисления Ѳ2 воспользуемся методом Дебая. Бу­ дем считать, что спектральная функция распределения dzt имеет вид (1.72) во всей области частот, и во всей этой области суммирование заменим интегрированием, прово­ дя его до предельной (дебаевской) частоты vD, зависящей от числа частиц в кристалле. Вычисляя Q2dzt, получим

Ѵп

V о_..2 / _Т7/.

где Тв — hvDlk — дебаевская температура,

о

— функция Дебая. Подставляя значение Ѳ2 из (1.78) в вы­ ражение (1.67), получим формулу для температурной за­ висимости частот ЯКР

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

Ф. И. Скриповым [25] *) были получены также фор­ мулы температурной зависимости частот ЯКР для следую­ щих случаев: 1) (а) направление оси качаний произвольно,

(б) участие молекулы в нескольких колебательных движе­

*) Работа [25] является лишь кратким изложением доклада на конференции и не содержит всех соотношений, независимо опубликованных Кушидой [27|.

2*

36 УРОВНИ ЭНЕРГИИ И ЧАСТОТЫ Я К Р [ГЛ. I

ниях; 2) конечная величина параметра асимметрии гради­ ента электрического поля; 3) влияние ангармоничности вращательных качаний; 4) влияние внутримолекулярных деформационных колебаний.

Как показано в [25], продольные акустические волны влияют на температурную зависимость частот ЯКР в 10 — 20 раз слабее, чем поперечные. Численные оценки пока­ зали, что при комнатной температуре поперечные аку­ стические колебания могут обусловливать до г!ъ всего температурного изменения частот ЯКР, а при очень низких температурах этот фактор может стать даже преобладающим.

Рассмотрим теперь, как влияет на температурный ход частот ЯКР участие молекулы одновременно в нескольких

 

вращательных

движениях,

 

например в двух вращатель­

 

ных качаниях вокруг

взаим­

 

но перпендикулярных осей

 

[8]. Такая ситуация,

напри­

 

мер, имеет место для мо­

 

лекулы парадихлорбензола

 

(п-С„Н4С12).

 

 

 

Использование теории ли­

 

нейных преобразований дает

 

возможность получить табли­

 

цу направляющих косинусов.

 

Качания вокруг оси Z не да­

Рис. 6. Колебания молекулы в ре­

ют вклада в температурную

шетке вокруг двух осей.

зависимость

частот

ЯКР.

Вращения вокруг осей Х и Y можно рассмотреть по отдельности, а затем перемножить матрицы преобразований.

Пусть молекула поворачивается вокруг оси X на угол Ф, а вокруг оси Y на угол Ѳ(рис. 6). Рассмотрим поворот в две стадии — сначала поворот вокруг оси Y, а затем вокруг оси X. Составим таблицы направляющих косину­ сов для обоих случаев в отдельности (табл. 2).

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

§ 4]

ТЕМ ПЕРАТУРНАЯ ЗАВИСИМОСТЬ ЧАСТОТ ЯКР

37

 

 

 

 

 

 

 

 

Т а б л и ц а 2

Поворот вокруг оси Y

 

Поворот вокруг оси X

 

X

Y

 

Z

 

 

X

У

z

X'

cos Ѳ

0

sin Ѳ

 

X'

1

0

0

Y'

0

1

 

0

 

Y'

0

COS Ф

sin ф

Z'

sin Ѳ

0

cos Ѳ

 

Z'

0

sin ф

COS Ф

 

 

 

 

 

 

Т а б л и ц а 3

 

 

 

X

 

Y

 

z

 

 

X'

cos Ѳ

 

sin Ѳ

sin ф

sin Ѳ

cos ф

 

 

Y'

0

 

COS ф

 

s i n

Ф

 

 

Z'

sin6

—cos Ѳ

sin Ф

cos Ѳ

COS Ф

 

электрического поля

 

 

 

 

 

 

qz.,z. =

qzz cos2 Ѳcos2 ф -

-і- qzz cos2 Ѳ• sin2 ф ---- Y qzz sin2 Ѳ;

 

Q zz (1 -

ѳ2) (1 -

ф2) -

-

f

(21- - 4- 0ѳ2) Ф

 

 

 

 

 

Q zz

1 - - И3 Ѳ 2 + Ф2)

(1.80)

(Считая углы 0 и ф малыми, мы в (1.80) пренебрегли членами второго порядка малости вида Ѳ2ф2.) Отсюда

пн

■(Ѳ? + Ф?)

(1.81)

у = ѵ0 1

т. е. в этом приближении эффекты суммируются (см. (1.66), (1.67)).

38

 

 

 

УРОВНИ ЭНЕРГИИ И ЧАСТОТЫ Я КР

[ГЛ. I

Найдем теперь Ѳх и фх:

 

 

 

 

 

 

 

 

 

2

Ѳі — Лѵі(

2

+

exp (hvi/кТ) — 1

(1.82)

 

 

 

 

 

 

 

 

 

 

 

 

 

-Г -!-

= hv2 (“Г

+

exp (hVi/кГ) — i

 

где

Ox — частота вращательных

качаний вокруг

оси Y,

со2 — частота вращательных

качаний молекулы

вокруг

оси X ,

и / 2

— моменты инерции молекулы относительно

осей Y

и X

соответственно. Отсюда получим

 

 

 

 

3h

 

 

 

 

 

 

 

 

 

8я2

1

exp (hvilkT) — 1 j

 

 

 

 

 

 

(-L +

_____ *___ _

(1.83)

 

 

 

 

 

+ /гѴг

\ 2

'

exp (hvz/kT) — 1

Продифференцировав (1.83) по T, получим формулу для

температурного

коэффициента

частот ЯКР

 

1

іѵ

_

 

 

 

 

 

 

 

 

Ѵо

dT

~~

 

 

 

 

 

 

exp (hvilkT )

 

_

 

ЗА2

'

 

exp (hv/кТ)

 

 

 

(1.84)

 

8л?кТ2

 

J i [exp (hvt/кТ) — l]2

 

/ 2 [exp (hvilkT) — l]2

Значение v0 получают обычно экстраполяцией экспе­ риментальной кривой у (Г) в О °К. Соотношение (1.84) можно применить для оценок порядка величины средних частот вращательных качаний и средних моментов инер­ ции молекул. Если положить vx Ä V2 = ѵ0 и J x ~ / 2 = = / 0, то

1 dv

3h?

exp (hv0/kT )

(1.85)

v^~dT

/tn?ki"iJ u

[exp (hxB/k.T) — l]2

 

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

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

§ *1

ТЕМ ПЕРАТУРНАЯ ЗАВИСИМОСТЬ ЧАСТОТ Я КР

39

[25,

26]. Для этого удобно

ввести параметр а г:

 

 

Ѳ| =

а?ф?,

(1.86)

где Ѳг — угловое смещение оси Z градиента электриче­ ского поля, вызванное вращательным движением вокруг оси, перпендикулярной Z, фг — угловое смещение оси Z, обусловленное вращением вокруг произвольной оси. Тог­ да температурная зависимость константы квадрупольного взаимодействия определится из уравнения

eQqz7,-eQqlz

щ

у

«?________ 1

,,

eQqzz

8я2

&

/.Vj exp (hv^kT) — 1 ’ '

'

верхним индексом 0 обозначена константа квадруполь­ ного взаимодействия при О °К.

Рис. 7. Колебания молекулы вокруг произвольной оси.

Если вращение молекулы происходит вокруг произ­ вольной оси, составляющей угол а с направлением Z, то в случае малых углов углы Ѳ и ф связаны соотношением QR = фг (рис. 7), откуда

 

Ѳ= ф-^- = фа{= фвіпа.

Таким

образом, а,- = sina, где sin а — синус угла

между

осью вращения и направлением химической связи

В общем случае для вычисления Ѳ2 (t) нельзя пользо­ ваться уравнением (1.69), поскольку оно справедливо, лишь если Ѳ(t) является нормальной координатой реше­ точных колебаний. Поэтому Q(t) целесообразно разложить по нормальным координатам решеточных колебаний [27]:

6 = 2 а*£і> Чгг — 7о(і + 2 ßi£i + + •••) > (1.88)

І

'

І

І,j

'

4P

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