Метод моментов в определении ширины линии магнитного резонанса
СОДЕРЖАНИЕ: Оглавление А.Введение ..2 §1.Локальное поле …2 §2.Общая теория магнитного поглощения …2Оглавление
А.Введение…………………………………………………………………………..2
§1.Локальное поле…………………………………………………………………2
§2.Общая теория магнитного поглощения………………………………………2
Б. Уширение, вызванное взаимодействием между одинаковыми спинами…….5
§3.Диполь-дипольное взаимодействие…………………………………………...5
§4.Определение моментов…………………………………………………………6
§5.Метод вычисления моментов………………………………………………….7
В. Кинетические свойства….………………………………………………………10
§6.Кинетическое уравнение………………………………………………………10
§7.Электропроводность…………………………………………………………...11
А. ВВЕДЕНИЕ
Линия магнитного резонансного поглощения системы спинов, находящихся в неоднородном магнитном поле, обладает некоторой шириной, обусловленной разбросом ларморовских частот. Аналогичное уширение может иметь место в неидеальных кристаллах благодаря взаимодействию ядерных квадрупольных моментов с малыми градиентами электрического поля, значения которых изменяются от одного узла решетки к другому случайным образом. В обоих случаях ширина линии обусловливается различием резонансных частот отдельных спинов, а не взаимодействиями между ними. Соответствующее уширение линии называется неоднородным уширением.
Положение существенно изменяется, если уширение линии обусловлено взаимодействием между соседними спинами. Эта задача и рассматривается в настоящей работе.
§ 1. ЛОКАЛЬНОЕ ПОЛЕ
Энергия взаимодействия между двумя ядерными спинами зависит от величины и ориентации их магнитных моментов, а также от длины и направления вектора, описывающего их относительное расположение. Влияние такого взаимодействия на ширину линии поглощения существенным образом зависит от того, зафиксирован ли этот вектор в пространстве или его положение быстро меняется со временем вследствие относительного движения ядер.
Последний случай, как правило, встречающийся в жидкостях и газах, будет рассмотрен позднее. В этой главе мы ограничимся случаем жесткой решетки, в которой ядра можно считать неподвижными. Такое приближение разумно для многих твердых тел при комнатной температуре, в частности для ионных кристаллов.
Энергия диполь-дипольного взаимодействия двух магнитных моментов m1 =g1 ћI1 и m2 =g2 ћI2 описывается хорошо известным выражением
(1)
которое можно переписать в виде
W12 = – m2 H12 = – g2 ћI2 H12 ,
где H12 — локальное поле, созданное первым спином в месте расположения второго спина. (Введение в рассмотрение понятия локального поля очень удобно.) Поскольку ядерные магнитные моменты имеют порядок 10-3 магнетона Бора, или 10-23 CGS, а между ядерные расстояния порядка нескольких ангстрем, то локальные поля в жесткой решетке в общем случае имеют порядок нескольких эрстед.
Взаимодействие двух одинаковых диполей в сильном поле Н0 может быть описано с классической точки зрения следующим образом. Первый диполь m1 прецессирует с ларморовской частотой вокруг поля Н0 и, следовательно, обладает постоянной составляющей вдоль этого поля и составляющей, которая вращается в плоскости, перпендикулярной полю. Постоянная составляющая m1 создает в месте расположения диполя m2 слабое постоянное поле, ориентация которого относительно Н0 зависит от взаимного расположения спинов. Если поле Н0 сильное, то на него заметно влияет только параллельная или антипараллельная ему составляющая слабого поля. Так как каждый спин в решетке имеет несколько соседей с различными относительными положениями и ориентациями, постоянная составляющая локального поля имеет разные значения в различных местах, что приводит к разбросу ларморовских частот и уширению линии.
Вращающаяся составляющая m1 создает в месте расположения m2 локальное магнитное поле, вращающееся с ларморовской частотой m1 , которая совпадает с ларморовской частотой для m2 . В свою очередь она имеет составляющую в плоскости, перпендикулярной Н0 и, следовательно, может заметно изменять ориентацию m2 благодаря явлению резонанса. Соответствующая ширина линии должна быть порядка величины вращающегося поля. В рассматриваемом случае оно того же порядка величины, что и локальное постоянное поле и, следовательно, вносит в уширение вклад сравнимой величины.
Необходимо отчетливо понимать, что механизмы, обусловливающие эти вклады в ширину линии, в действительности различны. Если два спина не являются одинаковыми, то вращающееся поле, созданное m1 , не является резонансным для m2 и оказывает на него пренебрежимо малое влияние, в то время как постоянное поле, созданное m1 , в месте расположения m2 является столь же эффективным, как и в случае одинаковых спинов. При прочих равных условиях одинаковые соседние спины оказывают более сильное влияние на уширение резонансной линии, чем неодинаковые.
§ 2. ОБЩАЯ ТЕОРИЯ МАГНИТНОГО ПОГЛОЩЕНИЯ
Для количественного описания формы линии, обусловленной дипольным уширением, необходимо развить формализм.
Когда все спины образца связаны друг с другом дипольным взаимодействием, представление об отдельных независимых спинах, находящихся в стационарных состояниях, становится неверным. Этот вывод следует хотя бы из того факта, что вращающееся локальное поле, созданное одним спином, приводит к переориентации его соседей. Поэтому образец приходится рассматривать как единую большую систему спинов, а переходы, вызванные радиочастотным полем, — как переходы между различными энергетическими уровнями этой системы. Соответственно изменяется и статистическое описание с использованием матрицы плотности. Вместо статистического ансамбля спинов, описываемых (2I +1) (2I +1) матрицей плотности, весь образец, содержащий N спинов, теперь становится одним элементом статистического ансамбля и описывается (2I +1)N (2I +1)N матрицей плотности. Такое видоизменение никоим образом не ограничивается ядерным магнетизмом, напротив, оно весьма часто встречается в статистической физике» а именно всякий раз, когда переходят от описания систем со слабыми взаимодействиями, например, таких, как молекулы газа при низком давлении, к описанию сильно взаимодействующих систем, таких, как атомы Кристалла. Первый подход соответствует методу Максвелла – Больцмана, а второй — методу Гиббса.
Стационарное состояние, следуя методу Гиббса, можно описать следующим образом. Если к системе спинов приложено линейно поляризованное вдоль оси Ох радиочастотное поле Н1 cos wt , то при стационарных условиях система приобретает намагниченность, составляющая которой вдоль этой же оси равна
Мх = H1 {c (w) cos wt +c (w) sin wt }. (la)
Условие линейности или отсутствия насыщения предполагает, что c и c не зависят от H0 . c и c можно измерить отдельно, а c пропорционально скорости поглощения радиочастотной энергии образцом.
Выведем общую формулу для c (w). Выше было показано, что в линейной теории резонанса между c (w) и c (w) существуют независимо от природы рассматриваемой системы общие соотношения (соотношения Крамерса – Кронига), позволяющие вычислить одну из этих величин, когда для всех значений частоты известна другая.
Ниже, чтобы избежать путаницы, мы будем обозначать через М макроскопическое значение намагниченности образца и через M — соответствующий квантовомеханический оператор. Между ними имеет место соотношение
М = M = Sp {rM }, (2)
где r – статистический оператор, или матрица плотности, описывающая систему спинов. Пусть H — полный гамильтониан системы в отсутствие внешнего радиочастотного поля. Если до приложения радиочастотного поля система находится в тепловом равновесии при температуре Т, то ее статистический оператор определяется выражением
(3)
которое просто означает, что статистическое поведение системы можно описать, если ее энергетическим уровням En приписать населенности, пропорциональные exp(— En / kT).
При наличии радиочастотного поля уравнение движения для r имеет вид
(4)
где V – объем образца. Чтобы решить (4) относительно r, сделаем подстановку
r* = ei H t r e – i H t , (5)
которая преобразует (4) в уравнение
. (6)
Предположим, что радиочастотное поле было включено в момент, когда образец находился в тепловом равновесии и
r (–) = r = r* (–).
В момент t решение (6) в линейном приближении относительно Н1 имеет вид
( 7)
Поэтому, возвращаясь к r [см. (5)], находим
(8)
Если предположить, что до включения радиочастотного доля намагниченность вдоль оси x была равна нулю, т. е.
Мх (–) = Sp {r0 M x } =0,
то
(9)
и, согласно определению (1 а),
(10)
Учтем, что температура обычно достаточно высока для того, чтобы для равновесной матрицы плотности (3) можно было использовать линейное разложение
где e – единичный оператор; тогда восприимчивость c(w) становится равной
(11)
откуда, интегрируя по частям, получаем
(12)
Выражение (12) можно преобразовать к более компактной форме двумя способами.
В первом способе, вводя в рассмотрение оператор Гейзенберга
M x (t) = e i H t M x e – i H t , (12a)
можно переписать (12) в виде
(13)
где
G(t) = Sp{M x (t) M x }, (13a)
Функцию G(t) назовем функцией корреляции, или функцией релаксации намагниченности системы.
Во втором способе выражение (12) можно переписать в виде
Отсюда после применения хорошо известной формулы для d-функции
получаем
(14)
где суммирование S производится только по тем энергетическим уровням, для которых | En — En | = w. Обычно, вводя в рассмотрение вероятности переходов, выражение (14) используют как отправную точку для вывода (13) с помощью интегрального представления d-функции. Из равенства (14) в общем виде следует, что функция формы f (w), определяющая форму линии, пропорциональна сумме S | п | M x | n ’ |2 . Точная зависимость этого выражения от co вытекает из условия, ограничивающего суммирование только по тем уровням, для которых | En — En | = w. Формулы (13) и (14) являются весьма общими и справедливы в случае, когда спектр магнитного поглощения системы содержит одну или несколько острых резонансных линий, т. е. в случае ядерного магнитного резонанса. Математически это условие может быть сформулировано следующим образом.
Гамильтониан H системы представляет собой сумму главной части H0 и малой возмущающей части, которую удобно записать в виде eH1 , где e — параметр малости возмущения. В отсутствие H1 спектр поглощения системы состоит из одной или нескольких бесконечно острых линий c частотами wa , a восприимчивость c(w) может быть записана в форме
c(w) = S Aa d(w-wa ); (15)
при этом функция релаксации G ( t ), пропорциональная фурье-преобразованию c(w), имеет вид
(15a)
Если существует возмущение eH1 , то функция релаксации принимает вид G(e, t) и может быть в принципе вычислена вплоть до любого порядка по e методом возмущений; восприимчивость c(w, e) получается как фурье-преобразование G(e, t).
Прежде чем производить детальный расчет, кратко рассмотрим соотношение между c(w) и поведением намагниченности после окончания действия радиочастотного импульса. Хорошо известно и достаточно очевидно, что для линейных систем стационарная реакция на возбуждение coswt представляется фурье-преобразованием нестационарной реакции на бесконечно острый импульс d(t). Однако на практике для аппроксимации такого импульса к системе спинов необходимо приложить кратковременно действующее магнитное поле, значительно большее постоянного поля Но .
Для системы взаимодействующих ядерных спинов в магнитном поле, характеризующейся острой резонансной линией на частоте w0 , действие бесконечно острого импульса постоянного поля можно аппроксимировать радиочастотным импульсом частоты w = w0 со значительно большей длительностью t и меньшей амплитудой H 1 . Поскольку в системе координат, вращающейся с частотой w, отлично от нуля только постоянное поле H 1 , то для аппроксимации бесконечно острого импульса конечной амплитуды достаточно того, чтобы H 1 было значительно больше локального поля; последнее представляет собой гораздо менее жесткое условие.
Б. УШИРЕНИЕ, ВЫЗВАННОЕ ВЗАИМОДЕЙСТВИЕМ
МЕЖДУ ОДИНАКОВЫМИ СПИНАМИ
§ 3. ДИПОЛЬ-ДИПОЛЬНОЕ ВЗАИМОДЕЙСТВИЕ
Полный гамильтониан системы одинаковых взаимодействующих спинов в сильном внешнем поле может быть записан в виде
H = (H 0 + H 1 ). (16)
Основной гамильтониан
H 0 = Sj Zj = – gH 0 Sj Ij z (16a)
описывает энергетические уровни, определяемые выражением E0 M = – gН0 M, где M — собственное значение оператора
Iz = Sj Ij z
Гамильтониан возмущения H1 , ответственный за уширение, имеет вид
(16б)
Прежде всего, рассмотрим несколько подробнее взаимодействие между двумя спинами, которые будем обозначать для краткости i и i’. Пусть q и j — полярные координаты вектора r , описывающего их взаимное положение, причем ось z направлена параллельно внешнему полю. Тогда Wii можно записать в виде
Wii = {i i — 3[iz cos q + sin q (ix cos j + iy sin j)]x[iz cos q + sin q (ix cos y + +iy sinj)]}g2 2 /r3 = {i i — 3[iz cos q + sin q (i+ e- i j + i- ei j )/2]x[iz cos q + sin q (i+ e- i j + + i- ei j )/2)]}g2 2 /r3 = (A+B+C+D+E+F)g2 2 /r3 , (17)
где
A = iz iz (l – 3cos2 q),
B = – (l – 3cos2 q) (i+ i– + i – i+ ) = (l – 3cos2 q)(iz iz – i i )/2,
C = – 3sinq cosq e- i j (iz i+ + i + iz )/2, (18)
D = С* = – 3sinq cosq e i j (iz i– + i – iz )/2,
E = – 3sin2 q e-2 i j i+ i+ /4,
F = E* = – 3sin2 q e-2 i j i – i– /4,.
Запись W в такой форме вызвана следующими причинами. Согласно формуле (14),
c(w) ~ S | п | M x | n ’ |2 .
Это приводит к необходимости определить изменение в положении энергетических уровней, отвечающих H0 , обусловленное наличием H1 . Операторы А, В, С, D, E, F дают качественно различным вклады в это изменение. Упомянутые операторы, действуя на состояние невозмущенного гамильтониана, характеризующееся значениями iz =т , i z =т, приводят к следующему изменению этого состояния:
(19)
Рассмотрим теперь энергетический уровень E0 M = – gH0 M, соответствующий гамильтониану (16a). Этот уровень сильно вырожден, так как существует много способов, которыми можно скомбинировать отдельные значения Ij z = mj , чтобы получить величину M = S mj . Таким образом, уровень E0 M соответствует вырожденному множеству состояний |М, причем вырождение снимается (по крайней мере частично) возмущением, описываемым гамильтонианом H1 , который расщепляет уровень E0 M на много подуровней. Согласно первому приближению теории возмущений, вклад первого порядка в расщепление уровня E0 M дают лишь те члены гамильтониана возмущения, которые обладают отличными от нуля матричными элементами внутри множества |М, т. е. те, которые, действуя на состояние |М, не вызывают изменения величины М. Обращаясь к формуле (19), мы видим, что только те части W, которые отвечают операторам А и В, удовлетворяют этому условию и должны быть сохранены для вычисления энергетических уровней H методом возмущений.
Член А имеет тот же вид, что и выражение для взаимодействия двух классических диполей и описывает упомянутое в разделе А взаимодействие одного диполя со статическим локальным полем, создаваемым другим диполем. Член В описывает взаимодействие, при котором возможно одновременное переворачивание двух соседних спинов в противоположных направлениях. Эта часть гамильтониана, названная «переворачивающей» частью, соответствует описанному в разделе А резонансному действию вращающегося локального поля. Влияние такого члена, как С, заключается в примешивании к состоянию |М с невозмущенной энергией E0 M = – gH0 M малой доли состояния |М— 1. Таким образом, точное собственное состояние H0 следует представить в виде
| М + a | М – 1 + …,
где a — малая величина. Взаимодействие системы спинов с радиочастотным полем, приложенным вдоль оси ох, пропорционально Ix = S Ij x и может индуцировать только переходы с DМ = ± 1. Слабые переходы знежду состоянием, скажем, |M – 2 + малая примесь, энергия которого приблизительно равна – gH0 (M —2), и состоянием | М + a | М – 1 + … становятся возможными с вероятностью порядка a2 . Разность энергии между этими состояниями приблизительно равна 2w0 . Следовательно, таким переходам на частоте 2w0 соответствует очень слабая линия, которую обычно трудно наблюдать экспериментально. Легко видеть, что линии сравнимых интенсивностей появляются на частотах 0 и 3w0 .
Доказательство справедливости сохранения в гамильтониане H1 только членов А и В, которые коммутируют с H0 обычно называются адиабатической или секулярной частью H1 и которые впредь будут обозначаться как H’0 , может быть также дано следующим способом. Так как c(w) пропорционально фурье-преобразованию G(t)=Sp{M x (t) M x }, то оно может быть вычислено, если известно M x (t) = е i H t M x е– i H t . В этом случае M x (t) удовлетворяет уравнению
(1/i) d M /dt = [H 0 +H 1 , M x (t) ]. (20)
§ 4. ОПРЕДЕЛЕНИЕ МОМЕНТОВ
Для резонансной кривой, описываемой нормированной функцией формы f(w) с максимумом на частоте w0 , n-й момент Mn относительно w0 определяется выражением
Мn = (w – w0 )n f(w)dw.
Если f(w) симметрична относительно w0 , то все нечетные моменты равны нулю. Знание моментов дает некоторую информацию о форме резонансной кривой и, в частности, о скорости, с которой она спадает до нуля на крыльях вдали от w0 .
Достоинство метода моментов состоит в том, что моменты могут быть вычислены на основании общих принципов без определения собственных состояний общего гамильтониана H . Прежде чем останавливаться на вычислении моментов, рассмотрим два примера резонансных кривых разном формы. Гауссова кривая описывается нормированной функцией
(24)
для которой легко найти
М2 = D2 , M4 =3D4 ,
М2n = 1, 3, 5, ..., (2n – 1) D2n ,
причем нечетные моменты равны нулю. Полуширина на половине высоты d определяемая соотношением f(w0 + d) = f(w0 )/2, или ехр( – d2 /2D2 ) = 1/2 оказывается равной
Отсюда видно, что значение второго момента M2 = D2 для гауссовой кривой обеспечивает удовлетворительное приближение для ширины линии d.
Другой формой линии, которая часто наблюдается в магнитном резонансе, является лоренцева форма, описываемая нормированной функцией
(25)
где d — полуширина на половине высоты.
В этом случае ни второй, ни более высокие моменты не могут быть определены, так как соответствующие интегралы расходятся. Однако иногда теория дает конечные значения для второго и четвертого моментов линий, которые в экспериментально наблюдаемой области имеют лоренцеву форму. В соответствии с конечными значениями M2 и М4 далеко на крыльях линии, где невозможно произвести достаточно точные измерения поглощения вследствие его малой величины, линия должна изменяться более быстро, чем это следует из лоренцевой формы.
Грубая, но удобная пробная модель состоит в описании кривой по формуле (25) внутри интервала |w – w0 |a, где ad и в предположении о том, что она равна нулю вне этого интервала. Тогда, пренебрегая членами порядка d/a, найдем
M2 = D2 = 2ad /p, M4 = 2a3 d /(3p), (IV.25a)
откуда, если известны M2 и M4 можно вычислить d и a. Поскольку
M4 /( M2 )2 = pa /6d,
упомянутая модель может быть использована лишь, когда теоретическое отношение M4 /( M2 )2 оказывается большим числом., В этом случае
(IV.25б)
Ширина на половине высоты значительно меньше, чем среднеквадратичная ширина. С другой стороны, предположение о гауссовой форме линии может быть разумным всякий раз, когда отношение M4 /( M2 )2 порядка 3.
§ 5. МЕТОД ВЫЧИСЛЕНИЯ МОМЕНТОВ
Основной недостаток метода моментов состоит в том, что важный вклад в значение момента (вклад тем существеннее, чем выше момент) дают крылья кривой, которые на практике не наблюдаются. Необходимо из вычисленных моментов линии магнитного резонанса с центром на ларморовской частоте w =w0 исключить вклады от сопутствующих линий на частотах w = 0, 2w0 , 3w0 о которых упоминалось ранее. Легко видеть, что, несмотря на их малую интенсивность (благодаря удаленности от центральной частоты w0 ) вклад во второй момент сравним с вкладом от главной линии и тем больше, чем выше порядок момента. Для исключения вкладов от них следует рассматривать в гамильтониане возмущения H1 ответственного за уширение, только его секулярную часть H 0 , которая коммутирует с H0 и, следовательно, не может отвечать перемешиванию состояний с различными полными М; такое смешивание является причиной появления побочных линий. Таким образом, сокращение дипольного гамильтониана до его секулярной части
не только упрощает вычисление моментов, но и делает его более точным.
Прежде чем начать расчет, отметим, что линия магнитного резонанса симметрична относительно центральной частоты w0 . Убедимся в правильности этого утверждения. Если | а и | b — два собственных состояния (H0 +H 1 ) с разностью энергии (Еа — Еb ) = w0 + dab , то два состояния | а~ и | b~ , полученные из | а и | b соответственно путем поворота всех спинов в обратном направлении, будут также собственными состояниями (H0 +H 1 ) с (Еb ~ – Еa ~) = w0 + dab . Таким образом, каждому переходу с частотой w0 + u соответствует переход равной интенсивности с частотой w0 – u. Если f(w) — функция формы, то h (u) = f(w0 + u)— четная функция u. Поскольку моменты кривой пропорциональны производным в начале координат от их фурье-преобразования, мы будем применять для их вычисления формулу (13). Вследствие узости линии ядерного магнитного резонанса можно пренебречь изменением величины w в пределах ширины линии и предположить, что форма линии описывается c(w)/w, так же как и c(w). Тогда, поскольку f(w) — нормированная функция формы, (13) может быть переписано в виде
f(w) = A G(t) cos wt dt, (IV.26)
где постоянная A определяется из условия нормировки f(w), а определенная ранее четная функция G (t) равна Sp{M x (t) M x }. Обратно
G(t) = 2/(pA) f(w) cos wt dw, (IV.27)
Согласно вышеизложенному, в выражении
M x (t) = е i H t M x е– i H t .
следует вместо H = H0 +H1 подставить H = H0 +H 1 что значительно упрощает вычисления. Поскольку H0 и H 1 коммутируют, можно записать
exp{i(H0 +H 1 )t} = exp(iH0 t) exp(iH 1 t).
Учитывая, что зеемановский гамильтониан H0 равен w0 Iz функцию G (t) можно переписать в виде
(IV.28)
Шпур произведения операторов инвариантен относительно циклической перестановки, поэтому
(IV.28a)
В этом выражении оператор exp(iw0 Iz t) определяет поворот на угол w0 t вокруг оси z, и, следовательно, можно записать
(29)
Легко видеть, что второй член в (29) равен нулю, так как поворот спинов на 180°, например вокруг оси ох, не изменяет H 1 и M x но преобразует M у в – M y .
Заменяя в (27) G (t) на G1 (t) cosw0 t, где
G1 (t)=Sp{е xp ( i H ‘1 t ) M x е(– i H ‘1 t )M x }
называется сокращенной функцией автокорреляции, и вводя обозначение
h (u) = f(w0 + u),
получаем
Заменяя нижний предел на – , что допустимо для узких линий, найдем
Поскольку h (и) является четной функцией, второй интеграл равен нулю и
G1 (t)=Sp{е xp ( i H ‘1 t ) M x е(– i H ‘1 t )M x }
(30)
Различные моменты кривой распределения h (и) относительно резонансной частоты w =w0 определяются выражением
Нечетные моменты равны нулю, а четные определяются формулой
(31)
Таким образом, для вычисления моментов резонансной кривой достаточно разложить G1 (t) в выражении (30) по степеням t. При этом коэффициенты разложения представляют собой шпуры от операторов, которые являются полиномами от H 1 и M x .
Сущность метода заключается в том, что значения упомянутых шпуров не зависят от выбора основных состояний и могут быть вычислены, например, в представлении, где значения mj = Ij z отдельных спинов (поэтому представление называется mj -представлением) являются хорошими квантовыми числами. Таким образом, нет необходимости решать проблему отыскания собственных состояний | n полного гамильтониана. Из определения (30) функции G1 (t) вытекает, что значение ее р-й производной в момент t = 0 определяется выражением
(IV.32)
Формула (32) просто находится из дифференциального уравнения
(33)
которому удовлетворяет зависящий от времени оператор
M x (t) = е(i H1 t) M x е(–i H1 t)t .
Решение этого уравнения может быть представлено в виде ряда
M x (t) = M x + M (1) x (t) + M (2) x (t) + …+ M (n) x (t) ,
отдельные члены, которого получаются методом индукции с помощью соотношения
из последнего сразу же следует (32). Из (31) и (32) для первых двух четных моментов находим
(34)
(34a)
B (34) M x заменено полным спином Ix , пропорциональным M x . Поскольку мы определили гамильтониан в виде H , следует помнить, что эти моменты соответствуют ширинам линии, измеренным в единицах w = 2pn.