
книги из ГПНТБ / Гречишкин В.С. Ядерные квадрупольные взаимодействия в твердых телах
.pdf40 |
УРОВНИ ЭНЕРГИИ И ЧАСТОТЫ Я КР |
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) |
$ 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. |
|
||
Поскольку |
2я |
и оц <^; со, членом с |
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-го |
(спиновое |
эхо) |
импульсов, |