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

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

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

40

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

t r n .' I

где і г =

liSin (югг + А,). Тогда в записи через нормаль­

ные координаты получим

 

 

ѵ = ѵ0(' і - 4 - 2І^ Г )' .

(1-89)

где А і = а г2— 2/3бгг. Среднюю энергию каждого нормаль­ ного колебания можно приравнять средней энергии кван­ товомеханического осциллятора:

4я2ѵ^Г

«... (1

+

_____1_____ )

 

(1.90)

2

ІМі ^ 2

~

exp (hv^kT) — 1/ ’

V — v0|l —

 

i vi

L

exp (hv./k'T) — 1

 

 

 

 

 

 

 

 

 

 

Если h \i/k T <^ 1,

то

 

 

 

 

V = v0

 

 

 

32яh?Ѣ Т 24

(1.91)

В общем случае А С1 имеет размерность и порядок величи­ ны момента инерции молекулярного колебания. Если т] ф 0, то

V ~ ѵ° [і —

Ѳ2 (*)]»

4 = (l - - f ) a f - А й « . (1.92)

 

Кушида, Бенедек и Бломберген [28] показали, что ча~ стота ЯКР является функцией давления. Поэтому фор­ мула (1.92) описывает температурную зависимость частоты ЯКР при постоянном объеме. В связи с этим можно на­ писать следующие термодинамические соотношения:

 

^

f

dqo

dv_ di\

,

-y dv_ /

 

\ j

/ dV \

~

 

V

[9g0

dV'

ÖT)

 

d V 'j

d£?

\ d V /Т)

V \ дР / т ’(1.

1

/

дѵ

\

V

[ dv

dqo

 

dv

df\

 

 

 

 

V

\

дТ

V

\ dq0

d V

dr\

dV'

 

 

 

 

 

 

 

 

/ ^ i \ 1 1

/

 

d V \

1

у дѵ / Щ \

 

 

 

+

^

дЦ {д Ѵ '/т }

V'

\

дТ

j p '

V

 

\

дТ) у .ш

Таким образом, в общем случае необходимо знать как температурную, так и объемную зависимость частоты

§ 4]

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

41

ЯКР, поскольку V = V (Р, W, Т). Наличие сильных объ­ емных эффектов в некоторых случаях даже обусловли­ вает положительный температурный коэффициент. С уменьшением объема (увеличением внешнего давления) амплитуда решеточных колебаний уменьшается, что и вызывает увеличение частоты ЯКР. Однако при уменьше­ нии объема увеличивается степень межмолекулярного взаимодействия, что приводит к уменьшению частоты ЯКР. Поэтому на кривых зависимости частоты ЯКР от давления в молекулярных кристаллах можно наблюдать максимум, положение которого зависит от температуры. В ионных кристаллах такой максимум не наблюдается.

Благодаря присутствию в членах суммы (1.91) фактора 1/ѵг2 вклад высокочастотных колебаний решетки значи­ тельно меньше, чем низкочастотных. Бломберген, Кушида и Бенедек используют уравнение (1.91) в виде

ѵ = а ( і + Ь Т + -%г ),

(1.94)

где

а = ѵ„, Ъ=

Константы а, Ъи с являются функциями объема. Суммиро­ вание в уравнении (1.94) проведено вплоть до частот, удовлетворяющих условию кшіІкТ = 1, где Т — низшая температура, использованная в эксперименте.

Чтобы определить а, Ь, с и ѵ через ^-изотермы (экспе­ риментально можно снять лишь зависимость ѵ (Р)), необ­ ходимо знать уравнение состояния твердого тела. Правда, в молекулярных кристаллах в силу слабого взаимодей­ ствия между молекулами в элементарной ячейке зави­ симость частоты ЯКР от объема невелика. Эта зависи­ мость может быть учтена с использованием линейной температурной зависимости решеточных частот [29, 30]

(йі = со°г (1 — qit),

(1.95)

где Qi — константа, зависящая от типа решеточных коле­ баний. Значения qt обычно заключены между 0,0050 и 0,0015 “К “1.

Существенное влияние на температурную зависимость частот ЯКР могут оказывать и решеточные движения с

42

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

[ГЛ. I

узким частотным спектром, например возникновение реориентациочной подвижности молекул. Такие движения мо­ гут значительно увеличивать температурный коэффициент и приводить к пропаданию линий квадрупольного резо­ нанса. Этот вопрос будет рассмотрен особо. Если известна кристаллическая структура вещества, то можно учесть влияние симметрии кристалла на температурную зави­ симость частот ЯКР [31]. Учет взаимодействия между ро­ таторами в элементарной ячейке приводит к следующему выражению для среднеквадратичного отклонения оси Z от равновесного положения:

П

<Ѳг>п = 2 ^<Ѳ?>п0,

(1.96)

J—1

 

где п — число молекул в элементарной ячейке, <Ѳ*)П — амплитудные отклонения в простой теории Байера, Bj — фактор, учитывающий межмолекулярные взаимодействия; он может быть протабулирован.

Иногда необходимо также учитывать и изотопический эффект в температурной зависимости частот ЯКР [8]. Так,

eQf„(Cl*)

например, отношение

в одном и том же соеди­

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

Г Jj Д В А ІІ

КВАДРУПОЛЬНОЕ СПИНОВОЕ ЭХО

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

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

рот вектора ядерного намагничения Ж = 2 Мч (где |аг —

г

магнитный момент ядра) вокруг направления магнитной компоненты электромагнитного поля. Импульс, который вызывает поворот вектораЖ на 90°, называется 90°-импуль- сом; тогда импульс в два раза большей длительности вызы­ вает поворот вектора Ж на 180° (180°-импульс), поскольку угловая скорость поворота Ж есть уНі, где — ампли­ туда радиочастотного поля. После воздействия 90°-им- пульса вектор Ж поворачивается на 90° и оказывается в плоскости X Y , прецессируя с частотой ѵ0 вокруг оси Z (магнитное поле направлено вдоль оси Z). Прецессирую­ щая намагниченность наводит в катушку, окружающую образец, сигнал индукции, который затухает во времени в силу различия скоростей прецессии отдельных слагаю­ щих вектора Ж за счет внутренних неоднородностей ло­ кальных полей в кристалле. Если через время т после дей­ ствия 90°-импульса, подействовать 180°-импульсом, то все векторы повернутся на 180°, так что процесс начнется

44

КВАДРУПОЛЬНОЕ СПИНОВОЕ ЭХО

[ГЛ. п

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

§ 1. Квадрупольное спиновое эхо

ваксиально-симметричном

иаксиально-несимметричном градиентах электрического ноля

При импульсном возбуждении ЯКР спиновые опера­ торы будут эволюционировать во времени, поэтому удобно провести рассмотрение эффекта с помощью статистического оператора матрицы плотности [1] с элементами

Pmm' = IР I т').

(2.1)

где индекс т нумерует состояния. Оператор р содержит всю необходимую информацию для описания статисти­ ческого спинового ансамбля. Изменение р во времени дается квантовым уравнением Лиувилля (уравнение Ней­ мана)

»»-§■ = [Ж, р] = Жр - рЖ,

(2.2)

где Ж — гамильтониан системы, і = Y — 1- Таким об­ разом, расчет спинового эха в ЯКР упирается в аккурат­ ное рассмотрение отклика системы, описываемой квадрупольным гамильтонианом MQ при том или ином выборе гамильтониана взаимодействия с радиочастотным полем Жі (решение уравнения Неймана или Шредингера).

Решение уравнения (2.2) можно записать в виде

Р (0 =

и {t) р (0)

и - 1 (t),

(2.3)

где U (t) — унитарный

оператор,

р (0) — начальное

зна­

чение матрицы плотности.

 

 

Если Ж не зависит от времени, то

 

U =

ехр | — і —

fJ .

(2.4)

Знание р позволяет вычислить математическое ожидание любой наблюдаемой величины Q для данной системы:

<<?> = Sp {рС}.

(2.5)

Рис. 8. Векторная модель квадрупольного спинового эха.

$ 1]

с и м м е т р и ч н ы й и н е с и м м е т р и ч н ы й г э п

45

До воздействия радиочастотных импульсов спиновая система находилась в равновесии, поэтому

Р (°) = W T T e^ { - ^ r ) ’

<2 *6>

где N — число частиц в системе, J —спин ядра. Вычислим полный магнитный момент системы <М2):

{Mzy = Sp {р (0) Mz) ж - j2j N^ kt Sp ( З Д ) - 0. (2.7)

Здесь использовано высокотемпературное приближение и выписан лишь первый член разложения. Таким образом,

вслучае произвольных спина

ипараметра асимметрии пол­ ный магнитный момент спино­ вой системы в ЯКР равен нулю. В этом и состоит главное отли­ чие ЯКР от ЯМР, где (Ж*) ф 0. С другой стороны, это разли­ чие не столь уж существенно, поскольку в ЯКР магнитные моменты подсистем отличны от нуля. Так, например, для J =

=3/2 для подсистем (-І-Ѵа, + 3/2)

и( —1/2, —3/2) магнитные мо­

менты равны

соответственно

+ М+т и — М_т.

Векторная

модель квадру-

польного спинового эха созда­ на Блюмом, Ханом и Эрцогом [2] (рис. 8). Поскольку линейно поляризованное радиочастотное поле Н1х может быть представле­ но в виде суперпозиции двух по­ лей, поляризованных по кругу с противоположным направле­

нием поляризации, то импульс радиочастотного поля вы­ зывает поворот на 90° магнитных моментов подсистем М+т и М_т. В отличие от ЯМР в ЯКР существенны оба направ­ ления поляризации радиочастотного поля. В дальнейшем

46

КВАДРУПОЛЬНОЕ СПИНОВОЕ ЭХО

ІГЛ. II

процесс образования спинового эха полностью совпадает с ЯМР.

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

и решеткой называют временем продольной релаксации Та)

втвердом теле возможен лишь при наличии сильных дви­ жений (случай «сильного сужения» линии).

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

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

рощенного варианта теории спинового эха.

присут­

Если внешнее магнитное поле Н 0 = 0, то в

ствий радиочастотного поля H x (t) — 2.fl\cos at

гамиль­

тониан запишем в виде [3]

 

Ж = Жя - ч М - Н х (1) = Ж а+Ж ъ

(2.8)

т. е. ось радиочастотной катушки параллельна оси X тен­ зора градиента электрического поля.

§ 13 СИММЕТРИЧНЫЙ И НЕСИММЕТРИЧНЫЙ ГЭП 47

Временная картина действующих импульсов приве­ дена на рис. 9. Радиочастотное поле действует на образец при 0 < t < tw, а для t > tw поле Н х = 0.

Рис. 9.

Сигнал спинового эха после двух импульсов.

Пусть J = 3/2, ц = 0.

Найдем решение

уравнения

Шредингера с

гамильтонианом (2.8) для 0 <it < tw:

 

ih'Vp

= Ж^>.

(2.9)

Решение уравнения (2.9) ищем в виде

 

 

Ф = 2

ст (t) <рте~ш™1,

(2.10)

 

т = —Ч ,

 

 

где сот = с0 — EJH,

Е т — собственные

значения

квадрупольного гамильтониана.

После подстановки (2.8) и (2.10) в (2.9) получим, ис­ пользуя ортогональность функций tpm,

ст = іоі \(т 11+I т — 1) ег<Ш,п

°>nb'l)<cm_1 +

+ I /_ I т +

1) e,(“m' “m+l) cm+1] cos соt, (2,11)

где а ! = у # ,.

48

КВАДРУПОЛЬНОЕ СПИНОВОЕ ЭХО

[ГЛ. II

Для J — 3/2 при точном резонансе

 

 

I ® -

(ш*л — «о»/,) I < т 4 г

 

получим

 

w 1

 

 

 

 

 

У з .

 

С;Н,2 — --- 2---

1(01С+ 1:2>

 

'

у з

 

( 2. 12)

С-у/, ——2 1 С±*/, + 2tCOiCzf1 /2cos сot.

 

Поскольку

и оц <^; со, членом с

cosco£ можно

пренебречь. Тогда решения уравнения (2.12) даются вы­ ражениями

с±ч, it) =

с+,/2 (0) COS

С О +

іс±ѵ, (0) sin

 

CO^ ,

c+4t it) =

с±*/2(0) cos

(O jt) +

fc+./2(0) Sin

 

(Olt} ,

 

 

 

 

 

 

 

 

(2.13)

 

 

 

2 K ( * ) ia =

i.

 

 

 

 

При t = 0

 

 

 

 

 

 

 

 

*/.

= y~2 ’

C~'k ^ ^ V

l ' C±Vs ^

=

° ’

(2 -14)

где ß — случайный фазовый фактор.

=

tw)

имеем

В конце действия 1-го импульса it

 

 

 

 

 

 

 

 

(2.15)

с-*/, (U

= еі(3 с+у, itw), с_у, (*„) =

eißс+ѵ, (О-

 

Для области t tw поле Н ± = 0 и решение уравнения Шредингера ищется в виде

•/.

 

Ф = 2 Фm*«i(Ue"*“m('~ 4

(2.16)

m = —3/2

где ст itw) определяются соотношениями (2.15). Теперь мы

§ 1]

СИММЕТРИЧНЫЙ И НЕСИМ М ЕТРИЧНЫ Й ГЭП

49

можем вычислить </ху после 1-го импульса *)

 

<Іху =

(Г I/+ ІФ) + (Ч>* I Н>) =

 

X

2

 

 

 

=

sin (/З м ^ ш) sin «о (t — Q ,

(2.17)

где ©о =

©у, «у, =

eQqJ2.

 

Всилу неоднородности градиента электрического поля

вобразце необходимо провести дополнительное усредне­ ние по всем частотам прецессии:

 

оо

\-Гх/“ср == § S (ö>0 ©о) (.^хУ Лщ =

 

—о

=

sin ( /З о н у sin ©о (г — Q exp [ — -у- {t — *ш)2] , (2.18)

где б — ширина линии гауссовой формы, g (щ0— ю0) — форм-фактор линии поглощения.

Аналогично получим решение уравнения Шредингера после 2-го импульса. Тогда для спинового эха имеем

<ІХУср = — - ^ s in (/3©i*w) sin2

X

 

X sin ©о (t — 2r) exp

(t — 2r)2- y j .

(2.19)

Это выражение получено в предположении, что

tw<^ т.

В дальнейшем, чтобы не нагромождать формул, мы

будем

везде считать, что длительность импульсов мала по срав­ нению с интервалом времени между ними.

В выражениях (2.18) и (2.19) также предполагается, что линия поглощения имеет гауссову форму [4]. Такой закон распределения частот поглощения очень часто при­ нимается для твердых тел. Однако если в твердом теле при­ сутствуют какие-либо быстрые внутренние движения, то

форма

линии

остается

лоренцевой.

В

этом случае

*)

<МЖ> = тгй</ж>,

поэтому сигнал индукции

или эха на за-

жимах

катушки

будет

d (.М

п — число витков катушки,

Ап — ^ — ,

А — площадь поперечного

сечения.

Поэтому

достаточно знать

</*> после 1-го

(индукция)

и 2-го

(спиновое

эхо)

импульсов,

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