.RU

СОДЕРЖАНИЕ ДИССЕРТАЦИИ - Компьютерная и математическая модель ядерного спиновОго эха 01. 04. 02 теоретическая физика


^ СОДЕРЖАНИЕ ДИССЕРТАЦИИ
Во введении дана общая характеристика работы, приводится обоснование актуальности работы, сформулирована её цель и сущность методов исследования, практическая ценность и научная новизна, а также сформулированы результаты, выносимые автором на защиту.

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

^ Во второй главе описывается расчётная модель ядерного квадрупольного спинового эха при наличии слабого локального или внешнего постоянного магнитного поля и её непосредственное использование при моделировании “биений” огибающей амплитуды спинового эха во внешнем или локальном магнитном поле (ОСЭ) для определения ориентации градиента электрического поля на монокристалле CdSb, направления и величины локального магнитного поля в соединениях Bi2Ge3O9 и Bi3B5O12.

В ЯКР экспериментах исследуемый образец помещался в катушку радиочастотного поля и подвергался воздействию радиочастотных импульсов. Радиочастотное поле катушки действовало на образец в течение времени , затем оно выключалось на время  -. В интервале  < t < + действовал второй импульс радиочастотного поля длительностью (считалось, что  >>,). После этого через время  регистрировался сигнал спинового эха, интенсивность которого пропорциональна производной d/dt от среднего магнитного момента ядерной подсистемы, направленного вдоль оси радиочастотной катушки. Схема уровней ядерного находилась с помощью вычисления уровней энергии квадрупольного гамильтониана [1, 2]. В главных осях (x, y, z) симметричного тензора градиента электрического поля (ГЭП) квадрупольный гамильтониан имеет вид [1, 2]:

, (1)

где eQ – ядерный электрический квадрупольный момент,

qzz – максимальная компонента тензора ГЭП,

eQqz'z'/h – константа квадрупольной связи,

I – спин ядра (например, для 209Bi он равен 9/2),

 =(qx'x'-qy'y')/qz'z' – параметр асимметрии ГЭП,

I=Ix  iIy.

Предполагалось, что в системе может присутствовать постоянное внешнее He или локальное Hi магнитное поле, взаимодействующее с магнитным моментом ядер образца. Независящий от времени гамильтониан имеет вид:

, (2)

Взаимодействие с радиочастотным полем Hrf представлялось в виде:

. (3)

Полный гамильтониан является суммой гамильтонианов H0 и H1:

. (4)

В результате импульсного воздействия радиочастотного поля операторы спинов ядер будут эволюционировать во времени, моделирование ядерного квадрупольного спинового эха проводилось с помощью матрицы плотности  [3, 12]. Зависимость  от времени находилась с помощью квантового уравнения Лиувилля [12]

, (5)

где H – полный гамильтониан системы (4), .

Предполагалось, что в начальный момент времени матрица плотности (0) задаётся нормированным распределением Больцмана:

, (6)

где ^ N – число ядер в образце,

k – постоянная Больцмана,

T – абсолютная температура.

Используя высокотемпературное приближение (eQqz'z' << kT), выражение (6) можно упростить

. (7)

Решение уравнения (5) с помощью оператора эволюции S [12] записывалось следующим образом:

, (8)

Оператор эволюции ^ S удовлетворяет уравнению Шредингера с начальным условием – единичной матрицей E: , . Ранее в работах по теории квадрупольного спинового эха, в которых применялся метод матрицы плотности (см., например, [4]), считалось, что оператор эволюции S имеет вид произведения матриц преобразования R и D, относящихся, соответственно, к периодам действия радиочастотного поля и к промежуткам времени между радиочастотными импульсами:

S = DnRn …D2R2D1R1. (9)

В формуле (9) Rn (t,tn0) – решение уравнения Шредингера в момент времени t в период действия n-го импульса ():

, (10)

где tn0 – время начала действия n-го импульса и tn - его продолжительность. Начало действия 1-го импульса t10 полагалось равным нулю. В формуле (9) Dn(t,tn0+tn)=exp[-(i/ħ)H0(t-tn0-tn)] – решение уравнения Шредингера в промежутке между n и n+1 импульсами (H1 = 0). В работе [4] в качестве начального условия при решении уравнений Шредингера для всех матриц R и D в (9) брали единичную матрицу E.

Используя представление взаимодействия [12], уравнение (10) можно переписать следующим образом [4]:

, , (11)

где .

В нашей работе использовался как подход с перемножением матриц (9) с последующим вычислением матриц Rn(t,tn0) по формуле (10) с единичным начальным условием, так и прямой численный расчёт оператора эволюции S, когда в качестве начального условия при расчёте оператора эволюции S после действия n-го радиочастотного импульса следует брать не E, а Dn-1Rn-1. Было показано, что расчётную схему можно значительно упростить и свести уравнение (11) к дифференциальному уравнению с постоянными коэффициентами в случае малого радиочастотного поля Hrf ,1=Hrf <<,  - гиромагнитное отношение для ядра, - частота квадрупольного перехода при выполнении дополнительных условий на длительность радиочастотных импульсов >>1 и ограничения на величину внешнего или локального магнитного поля He,i t1,2<<1.

В качестве конкретного примера рассматривался ЯКР на ядрах висмута и сурьмы (на поликристаллическом образце Bi2Ge3O9, монокристаллах CdSb и Bi3B5O12). Спин ядра 209Bi равен I=9/2, спин ядра сурьмы равен I=5/2, /2=6.92310-4 МГц/Гс - гиромагнитное отношение для ядра Bi, /2=1.02610-3 МГц/Гс - гиромагнитное отношение для ядра Sb. Условие малости магнитных полей с большим запасом выполняется в ЯКР экспериментах на оксидных соединениях висмута. Так, например, при наблюдении ЯКР на ядрах 209Bi в соединении Bi2Ge3O9 значение константы квадрупольной связи eQqzz/h равно 605.9 МГц, а частота перехода между уровнями ядерного спина 1/2-3/2 составляет ~160 МГц. Нетрудно оценить величину магнитного поля H, при котором H будет порядка – это поля в области 3-4 Тл. Данные поля намного превышают величины локальных полей в оксидных соединениях висмута (30-250 Гс), обнаруженных в ЯКР экспериментах [5, 8]. В зеемановских ЯКР измерениях обычно используются поля He до 500 Э. Значения Hrf в импульсных ЯКР спектрометрах также, как правило, не превышают несколько сотен эрстед, поскольку заметный сигнал спинового эха наблюдается лишь при условии ~1, где длительность радиочастотных импульсов , как правило, порядка 1 мкс (10-6 с). Для достижения оптимальных условий наблюдения сигнала спинового эха значение частоты 1 должно быть порядка 1106 с-1, что дает для Hrf оценку ~ 200 Гс. Поэтому условие малости магнитных полей можно считать выполненным. Кроме того, с хорошей степенью точности выполняется условие >>1.

Согласно общим правилам метода матрицы плотности [4, 12], средний спин ядерной подсистемы вычислялся следующим образом:

. (12)

В высокотемпературном приближении матрица плотности в начальный момент времени определяется формулой (7). Значение матрицы плотности в любой последующий момент времени находилось с помощью выражения (8). В случае малых магнитных полей пользуясь представлением базиса собственных функций оператора (2) и вводя обозначение (где , Ei – собственные значения гамильтониана (2)), были получены выражения для вычисления среднего спина ядерной подсистемы после первого (формула (13)) и второго (формула (14)) радиочастотных импульсов

, (13) где , ;

, (14) где , .

Формулы (13), (14) использовались для вычисления среднего ядерного спина при наличии постоянного магнитного поля. Чтобы построить ОСЭ, при каждом значении  находился абсолютный максимум выражения (14) при t~2. Причём при вычислении амплитуды спинового эха в момент времени t~2 в выражении (14) оставлялись только вклады при которых отношение частот , а частота соответствует резонансной частоте или близкой к ней.

Чтобы смоделировать форму линии сигнала спинового эха, в силу неоднородности градиента электрического поля в образце проводилось дополнительное усреднение по всем частотам прецессии с форм фактором , где  – ширина линии гауссовой формы. При численном расчёте считалось, что константа квадрупольной связи в квадрупольном гамильтониане (1) распределена по Гауссу. После такой процедуры в выражении (14) выделялись вклады в сигнал индукции после первого и второго импульсов, а так же сигнал спинового эха при t~2.

Поскольку в гамильтониане (4) эффекты релаксации не учитывались, то простейшим образом моделирование релаксации сигнала спинового эха осуществлялось путем умножения амплитуды спинового эха на фактор . Значение времени релаксации T2 подбиралось из эксперимента.

Процедура моделирования сигналов квадрупольного спинового эха учитывала элементарные, но громоздкие геометрические факторы, описывающие взаимное расположение осей лабораторной () и кристаллографической (x, y, z) систем координат. Эту связь можно описать с помощью углов Эйлера ,  и . Оси лабораторной системы координат обычно определяются направлением радиочастотного поля (ось ) и внешнего постоянного магнитного поля (ось ). Связь между системой главных осей тензора ГЭП () и кристаллографической системой координат задается еще одним набором углов Эйлера ,  и . Ориентация вектора локального магнитного поля по отношению к главным осям тензора ГЭП определяется двумя углами i и i. Задание углов , , , , , , i, i полностью определяет направление магнитных полей , и относительно системы осей (), в которой написан гамильтониан (2). Все скалярные произведения , и , входящие в гамильтонианы (2) и (3), также следует записывать в той же системе координат. Для этого в расчётную схему включалась матрица вращений W, являющаяся функцией углов , , , ,  и , которая проектирует векторы и на главные оси тензора ГЭП. Поскольку в соединениях может быть несколько типов тензоров ГЭП, каждому из которых отвечает собственная система главных осей, расчёт среднего спина ядерной подсистемы производился для каждого направления главных осей ГЭП в отдельности (например, для двух направлений тензора ГЭП в соединении Bi3B5O12 и четырёх направлений тензора ГЭП в соединении CdSb). Результирующее значение среднего спина ядерной подсистемы являлось суммой вкладов от всех типов тензоров ГЭП. Для поликристаллического образца (порошка) вычисления производились тем же самым образом. Отличие заключалось в том, что рассчитанный средний спин усреднялся по всем ориентациям кристаллитов относительно осей лабораторной системы координат. Это достигалось усреднением по трем углам Эйлера ,  и .

Тем самым, расчётную схему для моделирования ОСЭ можно разбить на следующие этапы:

  1. Симметрийный анализ тензоров ГЭП. Нахождение матрицы вращений W. При необходимости расчёты проводятся для нескольких типов ГЭП.

  2. Вычисление собственных значений и собственных функций гамильтониана (2).

  3. Расчёт матричных элементов операторов в представлении собственных функций оператора (2).

  4. Решение уравнения Шредингера (11) для нахождения матрицы с единичным начальным условием.

  5. Вычисление коэффициентов и в формулах (13, 14).

  6. Нахождение абсолютного максимума выражения (14) при построении огибающей амплитуды спинового эха.

  7. Учет релаксации сигнала спинового эха путем умножения на фактор .

Во второй части второй главы изложено применение расчётной схемы для анализа ЯКР экспериментов на оксидных висмутовых соединениях и соединении сурьмы. Эксперименты на импульсном ЯКР спектрометре проводились в Институте Общей и Неорганической Химии РАН группой под руководством Э.А. Кравченко. На поликристаллическом образце Bi2Ge3O9 была получена отчетливая картина биений на огибающей амплитуды спинового эха в отсутствие внешнего магнитного поля. Данное наблюдение является с


Рис. 1. Экспериментальные (1) и

рассчитанные (2) огибающие спинового эха линий переходов 1 (a) и 2 (b) для порошка Bi2Ge3O9 в нулевом магнитном поле.
видетельством наличия в соединении Bi2Ge3O9 локальных магнитных полей. Результаты моделирования отражены на Рис. 1. Установлено, что модельная форма кривой ОСЭ для перехода 1 = 25.25 МГц, хорошо согласуется с экспериментальной в предположении, что ло­кальное магнитное поле в соединении Bi2Ge3O9 равно Hi = 65 ± 5 Гс и направлено под углом i = 83 ± 1° к оси qz'z' ГЭП. Полученная из этого сравнения оценка локального магнитного поля согласуется со значением магнитного поля 60 ± 10 Гс, полученного из анализа формы линий ЯКР в нулевом поле в работе [8].Что касается переходов выше, чем 1, то для них, как показано в работе [4], биения на кривых ОСЭ должны отсутствовать, если параметр асимметрии  = 0. Однако моделирование ОСЭ показало, что в предположении о ничтожно малом отклонении ГЭП от аксиальной симметрии ( ~ 0.005) моду­ляции ОСЭ линии перехода 2 = 50.5 МГц хорошо воспроизводятся при тех же характеристиках локального поля Hi, которые были найдены при моделировании ОСЭ линии перехода 1 (см. Рис. 1). Естест­венно, столь малое значение  лежит за пределами точности измерений ЯКР и указывает на высокую точность измерения свойств вещества методом огибающей спинового эха. Представляется интересным вопрос о воз­можности обнаружения соответствующего отклонения в рас­пределении электронной плотности от аксиальной сим­метрии с помощью рентгеновского эксперимента.

Также было выполнено моделирование формы ОСЭ линий переходов для соединения Bi3B5O12 с более сложной структурой (Рис. 2). Эксперимент выявил отчётливую модуляцию огибающей линии 2 в нулевом внешнем магнитном поле для атома Bi(1) (см. Табл. 1) как в монокристалле, так и в порошкообразном образце Bi3B5O12 (Рис. 3), что свидетельствует о наличии в соединении локального магнитного поля.



Расчёт формы ОСЭ с помощью разработанной программы позволил оценить значение локального магнитного поля (Hi = 2.25 Гс, см. Рис. 2) в соединении Bi3B5O12 и его ориентацию по отношению к главной оси qz'z' тензора ГЭП (i,= 22.5 2.5), а так же чувствительность метода квадрупольного спинового эха. Изменение локального поля на 0.5 Э или угла, ориентирующего это поле, на 2.5 приводило к существенному изменению сигнала амплитуды спинового эха.





На монокристалле CdSb в отсутствие внешнего магнитного поля биений амплитуды спинового эха не обнаружено. Однако, во внешнем поле, равном 16 Э и параллельном радиочастотному полю, на нижнем квадрупольном переходе 1/2-3/2 наблюдалась отчетливая картина биений (Рис. 4). Моделирование модуляции ОСЭ во внешнем поле равном 16 Э (параллельном радиочастотному полю) с использованием матрицы вращений W позволило определить ориентацию градиента кристаллического электрического поля в данном соединении. Углы =53.6 и =57 определялись из анализа спектров ЯКР 121Sb. Угол  невозможно было определить по спектрам ядерного квадрупольного резонанса. При моделировании ОСЭ оказалось, что наилучшее соответствие экспериментальной картинке происходит при угле =111 (см. Рис. 4).



Все расчёты проводились в соответствии с двухимпульсной методикой экспериментальной регистрации спинового эха. Предполагалось, что длительности радиочастотных импульсов равны мкс и мкс соответственно. Радиочастотное магнитное поле Hrf подбиралось из условия ~1, при котором сигнал спинового эха максимален. Например, при моделировании ОСЭ на ядрах 121Sb в монокристалле CdSb это соответствовало радиочастотному полю ~ 80 Гс. При сравнении экспериментальных и численно рассчитанных кривых ОСЭ важным критерием являлось количество максимумов и минимумов (например, см. Рис. 4), а так же совпадение основных особенностей на временной шкале расстояния между импульсами . При этом амплитуды рассчитанного и измеренного сигнала спинового эха могли отличаться, поскольку эффекты релаксации были учтены феноменологически, путем умножения рассчитываемой амплитуды спинового эха на экспоненциальный коэффициент, содержащий эффективное время спин-спиновой релаксации, определяемое по экспериментальному времени затухания сигнала в импульсном ЯКР спектрометре. Хорошее совпадение экспериментально измеренных и теоретически рассчитанных кривых ОСЭ на рисунках 1, 2 и 4 позволяет утверждать о состоятельности компьютерной модели квадрупольного спинового эха. Модель применима также и в случае ЯМР.

^ В третьей главе излагается аналитический метод описания явления спинового эха, основанный на теореме Флоке-Ляпунова [13, 14], применимой к решению уравнений с периодическими коэффициентами. Ранее эта теорема была использована Б.Н. Провоторовым и Э.Б. Фельдманом для теоретического описания динамики спиновой системы в многоимпульсных экспериментах [15]. Третья глава разделена на шесть частей. В первой части описана структура третьей главы, формулируются основные отличия резонансного и нерезонансного случая теории возмущений с применением теоремы Флоке-Ляпунова, перечисляются основные свойства оператора эволюции в случае ЯМР со спином I = 1/2 и ЯКР со спином ядра I =1. В первом разделе отмечается, что теорема Флоке-Ляпунова позволяет выделить периодическую и непериодическую части решения уравнения (10), разработать корректную теорию возмущений по малому параметру, в которой не содержалось бы так называемых вековых членов, и тем самым, справедливую на всём интервале эволюции системы. В соответствие с теоремой Флоке-Ляпунова решение системы (10) может быть представлено в виде:

. (15)

Постоянная матрица ^ K в выражении (15) определяется решением системы уравнений (10) за период системы T=2/, . Матрица F является периодической функцией времени , F(0)=Y(0)=E.

После проведения сравнительного анализа с существующими моделями можно заключить, что к началу нашего исследования аналитических моделей спинового эха (2004 год) стояла задача обобщить разработанный и подробнейшим образом изложенный в работах [10, 11] нерезонансный метод на основе теоремы Флоке-Ляпунова на резонансный случай, когда частота радиочастотного поля совпадает с одной или несколькими частотами переходов основного стационарного гамильтониана (2). Данную проблему удавалось решить в задачах о спиновом эхо в ЯКР экспериментах со спином ядра I=1 или ЯМР спиновом эхо со спином ядра I=1/2 путём перехода в представление взаимодействия (11). В этом случае уравнение (11) можно решить, комбинируя теорему Флоке-Ляпунова (15) со стандартной схемой решения, используемой в нерезонансном случае и описанной во втором разделе третьей главы, а также в работе [10]. После перехода в представление взаимодействия, введения безразмерного времени и параметра система (11) в случае ЯКР, I=1 представлялась в виде:

, . (16)

В этом же разделе было показано, что система дифференциальных уравнений для матрицы в случае ЯМР, I=1/2 может быть получена из выражения (16) заменой , и отбрасыванием тривиального нулевого второго столбца и второй строчки в периодической матрице A.

В случае малого параметра , использовался следующий алгоритм нахождения аналитических выражений для матриц F и K в формуле (15) [10, 13, 14]. Матрицы F и K искались в виде разложений в ряд по малому параметру  :

, , (17)

, , i=1, 2, 3,…; , j=0, 1, 2,… (18) Применение теоремы Флоке-Ляпунова (15) позволяет переписать уравнение с периодическими коэффициентами (16) в виде

(19)

После подстановки разложений (17) в систему уравнений (19) получаем рекуррентные соотношения:

(20) Решение уравнения нулевого порядка в (20) может быть записано в виде . В силу периодичности , следовательно, и . Для того, чтобы определить матрицы (n=1, 2,…), необходимо использовать условие периодичности матриц (см. (18)). Проинтегрировав систему дифференциальных уравнений (20) по периоду T=, получим систему алгебраических уравнений для нахождения матрицы по уже вычисленным на предыдущем шаге матрицам и . После того, как матрица вычислена, выражение для матрицы-функции получается из уравнений (20) простым интегрированием правой части за произвольный промежуток времени.

Используя описанный выше подход, во второй части третьей главы были получены выражения для вычисления среднего момента после действия второго импульса в случае ЯКР, I=1:



где коэффициенты A3 и B3, определяющие амплитуду спинового эха в момент времени t~2, с точностью до с точностью до членов первого приближения по параметру задавались выражениями:

(21)

(22)

Ф


Рис.5. Рассчитанный сигнал ядерной индукции и спинового эха при значении параметра =0.78. Длительности импульсов t1=210-6 с, t2=410-6 с, время начала действия второго импульса  =1.210-5 с, eQqzz/h=4 MHz, =0 (Случай ЯКР на ядрах азота 14N, спин ядра I=1).
ормулы (21, 22) отчётливо демонстрируют появление множественных сигналов индукции вблизи интервала времени t ~  и множественных сигналов спинового эха при t~2 (см. Рис. 5) в случае сильного радиочастотного поля. Этот результат может быть использован в качестве основы для создания нового типа спектрометров, в которых сигнал спинового эха измерялся бы не только в момент времени 2, но и на всём промежутке времени после воздействия второго импульса. В коэффициенты Ai, Bi входят вклады вида и , содержащие резонансную частоту . В силу неоднородности градиента электрического поля после свертки момента с распределением Гаусса , , получаем дополнительные сигналы, удалённые от основных сигналов индукции (t = 0, t = ) и спинового эха (t = 2) на времена кратные длительностям радиочастотных импульсов , , и так далее. Амплитуды таких сигналов будут пропорциональны соответствующим порядкам параметра (, , и так далее). Таким образом, явление спинового эха в двухимпульсной методике возникает в результате интерференции действия когерентных радиочастотных полей от двух радиочастотных импульсов на ядерный спин, находящийся в кристаллическом электрическом поле. В простейшем случае ЯКР, I=1 и линейной поляризации радиочастотного поля постоянная матрица K (15) имела чисто мнимые собственные значения , называемые показателями Флоке системы уравнений (11). Показатели Флоке отражают свойства физической системы, в частности мнимость этих показателей свидетельствует о том, что релаксационные эффекты не учитывались в общем гамильтониане (4). На Рис. 6, изображающего зависимость характеристического показателя Флоке от параметра , отчётливо видно существование двух режимов поведения ядерного спина в радиочастотном поле , и . В случае слабого радиочастотного поля (, ) установлена возможность разделения поведения ядерного спина на быстрые (с частотой, определяемой переходами между уровнями основного гамильтониана) и медленные (обусловленные поворотом спина вокруг направления радиочастотного поля на угол ) движения. В случае сильного радиочастотного поля такое разделение невозможно. Формулы моделей Блума-Хана-Герцога и Даса-Саха [3, 4] в нулевом приближении по при условии совпадают с более точными выражениями (21, 22).

Р


Рис.6. Зависимость характеристического показателя Флоке от параметра .
ассматривая линейную поляризацию как суперпозицию двух круговых поляризаций, в третьем разделе третьей главы предложена простая интерпретация формул вида (21, 22), полученных во втором разделе. В четвёртом разделе третьей главы излагается обобщение стандартного алгоритма использования теоремы Флоке-Ляпунова для решения задач спектроскопии в нерезонансном случае на случай произвольного стационарного гамильтониана (2) с неэквидистантными уровнями спина ядра. В этом случае ввиду появления значительного количества несоотносимых периодов в преобразованном гамильтониане (11), переход в представление взаимодействия неэффективен для использования теоремы Флоке-Ляпунова и Теорема Флоке-Ляпунова применялась непосредственно к уравнению (10). В случае, когда частота радиочастотного поля совпадает с одной или несколькими частотами стационарного гамильтониана (2), итоговые формулы, получаемые в этом разделе, а также в работах [10, 11], не определены. В пятом разделе третьей главы рассматривается модификация методов, изложенных в предыдущих разделах, на резонансный случай с произвольным спином ядра. Для того чтобы не менять стандартную схему вычислений изложенную в четвёртом разделе, в нашей работе в резонансном случае было предложено специальное преобразование системы (10). Например, если считать, что постоянный гамильтониан в формуле (10) в собственном представлении имеет диагональный вид и частота перехода совпадает с частотой радиочастотного поля , то необходимо произвести следующее преобразование матрицы в матрицу :

, .

Благодаря такой замене, полученные уравнения для матрицы будут по-прежнему представлять собой дифференциальные уравнения с периодическими коэффициентами с периодом T=2/, для решения которых допустимо использовать стандартный подход с использованием теоремы Флоке-Ляпунова, описанный в четвёртом разделе третьей главы. В результате были получены обобщённые формулы, позволяющие вычислять матрицы F и К с точностью до первого приближения по параметру  в резонансном случае.

В заключительном шестом разделе третьей главы разработанный в пятом разделе резонансный метод использовался для описания динамики двух спинов в модели Изинга при конечной температуре. Данная задача актуальна в связи с использованием ЯМР спектроскопии в теории квантовых вычислений [16]. В результате, был установлен вид оператора эволюции в нулевом приближении по параметру . Поскольку линейную поляризацию радиочастотного поля можно рассматривать как суперпозицию двух круговых поляризаций, в нулевом приближении, как и ожидалось, полученный результат сходился с выражением для оператора эволюции, получаемого в точно решаемой задаче о динамике двух спинов в модели Изинга с круговой поляризацией радиочастотного поля [16]. Разработанный нами метод позволяет получать поправки более высокого порядка по параметру для оператора эволюции и может быть использован для оценки реализуемости квантовых вычислений на ядерных спинах.

^ В четвёртой главе излагается заключительный этап работы — включение релаксационных явлений в разработанную компьютерную модель ядерного квадрупольного спинового эха. В качестве основного гамильтониана взаимодействия рассматривался гамильтониан вида:

(23)

Первый член в гамильтониане (23) - квадрупольный гамильтониан взаимодействия ядерного спина с градиентом кристаллического электрического поля (1), определяющий систему уровней ядерного спина и не зависит от времени. Последующие три члена в формуле (23) - гамильтонианы взаимодействия ядерного спина с флуктуирующим магнитным полем. Проекции флуктуирующего магнитного поля задавались тремя независимыми, нескоррелированными между собой стохастическими телеграфными процессами с амплитудами [17]. Случайные функции могут принимать значения равные 1 и -1. Основной характеристикой случайного процесса являлось время корреляции, зависящее от параметров системы - температуры, частот переходов квадрупольного гамильтониана. Времена корреляции определяли характер зависимости от времени среднего значения телеграфного процесса и его корреляционной функции : , , где .

Механизм релаксации, вызванный взаимодействием с флуктуирующим магнитным полем, мы называли магнитным и описывали с помощью магнитного времени корреляции . В качестве немагнитного механизма релаксации рассматривалось флуктуационное взаимодействие, определяемое орторомбическими искажениями кристаллической решётки и задаваемое гамильтонианом: , где . Такой немагнитный механизм релаксации мы называли квадрупольным механизмом релаксации, описываемым квадрупольным временем корреляции . Поведение случайной величины задавалось телеграфным процессом со средним значением и функцией корреляции, определяемыми выражениями: , .

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

При описании физической системы использовался формализм метода матрицы плотности . Зависимость  от времени находилась с помощью квантового уравнения Лиувилля (уравнения Неймана) [12]. В четвёртой главе методами теории Блоха — Вангснесса — Редфилда [18, 19, 20] были получены обобщённые стохастические уравнения Блоха с последующим их решением в случае ЯКР и спина ядра I = 1. Несмотря на то, что с помощью методов теории теории Блоха — Вангснесса — Редфилда давно уже получены уравнения, описывающие релаксацию наблюдаемых значений квантовых операторов спина ядра , и в случае ЯМР (см. [18]), случай ЯКР недостаточно изучен. Это связано с неэквидистантностью уровней спина ядра в ЯКР и, следовательно, значительным количеством времён релаксации. Некоторые общие формулы, описывающие релаксацию ядерного спина в ЯКР приведены в книге [2], однако детального рассмотрения процесса релаксации и получения выражений для времён поперечной и продольной релаксации аналогичных тем, которые хорошо известны в ЯМР (см. [18]) не осуществлено.

Работоспособность теории Блоха — Вангснесса — Редфилда была проверена численно путём моделирования зависимости z-проекции момента от времени в случае ЯМР со спином ядра I = 1/2 при наличии флуктуирующего магнитного поля, направленного вдоль оси x, с амплитудой x0, где - случайная величина, описываемая телеграфным процессом и принимающая значения 1 и  1 (см. Рис. 7 и Рис. 8).

Уравнение, позволяющее рассчитать вероятность в момент времени t случайной величине иметь значение 1 или -1, называется управляющим уравнением [17]. Генерирование амплитуд 1 и -1 случайной величины с вероятностью, рассчитанной из управляющего уравнения, с последующей подстановкой в динамическое уравнение на матрицу плотности (5), позволило получить одну из возможных реализаций матрицы плотности , а так же вероятность возникновения данной реализации. Пользуясь квантомеханическим средним можно получить одну из в


Рис.7. Сплошная кривая - одна из возможных реализаций зависимости z-проекции момента в случае ЯМР со спином ядра I=1/2. Расчет был произведён при следующих значениях параметров: частота Зеемана, определяемая постоянным магнитным полем, ωZ=2· 106 с-1; время корреляции телеграфного процесса tM=10-6 c, амплитуда флуктуирующего магнитного поля gx0= 105 с-1; начальное значение проекции момента (t=0) = 4.5, gx0<<ωZ. Пунктирная линия - это теоретически расчитанная кривая релаксации усреднённого по всем реализациям момента методами теории Блоха — Вангснесса — Редфилда c временем релаксации Т1=500 мкс.



Рис.8. Сплошная кривая - усредненная по 200 реализациям с помощью выборочного среднего зависимость z-проекции момента от времени в случае ЯМР со спином ядра I=1/2. Расчет кривой на рисунке 8 был произведён при тех же значениях параметров, что и при моделировании одной из возможных реализаций z-проекции момента на рисунке 7. Пунктирная линия - это теоретически рассчитанная кривая релаксации усреднённого по всем реализациям момента методами теории Блоха — Вангснесса — Редфилда c временем релаксации Т1=500 мкс.

озможных реализаций для наблюдаемой величины Iz (см. Рис. 7).

Усредняя по 200 возможным реализациям, получаем оценочное полное статистическое среднее спинового ансамбля наблюдаемой величины Iz (см. Рис. 8). Поскольку число возможных реализаций даже в случае использования простейшего телеграфного процесса для описания случайной величины экспоненциально увеличивается с ростом времени t, численный расчёт среднего значения наблюдаемой величины описанным выше способом требует значительных вычислительных ресурсов. Другим подходом к решению данной задачи является аналитический вывод уравнений, описывающих поведение среднего значения наблюдаемой величины — квантовых стохастических уравнений Блоха.

Для решения обобщённых уравнений Блоха в случае ЯКР, I=1 потребовалось разработать алгебру частичных спиновых операторов, позволяющую получить замкнутые уравнения, описывающие релаксацию наблюдаемых значений спиновых операторов , и . Основным свойством данных частичных операторов является их связь с обычными операторами спина : , , . Разности частичных операторов, умноженных на мнимую единицу и обладающие свойством эрмитововсти, были названы дополнительными спиновыми операторами: , , . Например, замкнутые уравнения, описывающие релаксацию наблюдаемой величины в случае ЯКР, I=1 примут вид:

(24)

где ,, , , - частоты переходов между уровнями квадрупольного гамильтониана (1) в случае спина ядра I = 1.

Путём диагонализации матрицы правой части уравнений (24) были получены два собственных значения, действительная часть которых определяет обратные времена релаксации . Времена релаксации имели наиболее простой вид в случае нулевого параметра асимметрии, когда квадрупольный гамильтониан (1) определяет одну ненулевую частоту перехода, ,:

, .

Аналогично были получены два поперечных времени релаксации, соответствующих релаксации x компоненты и два времени , соответствующих релаксации y компоненты. В простейшем случае, когда параметр асимметрии равен нулю и флуктуирующее поле одинаково вдоль осей x и y получалось два времени поперечной релаксации:

,

где .

Тем самым было установлено, что в случае ЯКР со спином ядра I=1, по сравнению с ЯМР, количество времён релаксации удваивается, а вклады в выражения для обратных времён релаксации, связанные с магнитными и квадрупольными механизмами релаксации, складываются.

^ В пятой главе проводится подробное обсуждение результатов работы.

В заключении коротко просуммированы основные результаты и приведены главные выводы работы.

Работа содержит три приложения: "Упрощение уравнения Шредингера в случае малости магнитных полей и большой константы квадрупольной связи. Явное выражение для элементов матрицы эволюции", "Стохастические процессы. Телеграфный процесс", "Свойства частичных спиновых операторов в случае ЯКР, I=1". В приложения входят не вошедшие в основной текст математические выкладки и определения.

^ Основные результаты и выводы

  1. Установлено, что результирующие формулы существующих моделей Блума-Хана-Герцога и Даса-Саха, положенные в основу работы большинства ЯКР и ЯМР спектрометров, справедливы в нулевом приближении по параметру , - частота Раби, - амплитуда радиочастотного поля,  - гиромагнитное отношение ядра,  - частота радиочастотного поля, при условии выполнения соотношения на длительность радиочастотного импульса tw: .

  2. Разработана компьютерная модель ядерного спинового эха, в которой уравнения для оператора эволюции решаются численно как в случае сильного, так и в случае слабого радиочастотного поля. Численная модель была протестирована при анализе экспериментально измеренных временных зависимостей огибающих амплитуды спинового эха на поликристаллическом образце Bi2Ge3O9, на монокристалле Bi3B5O12 и монокристалле CdSb. Сравнение экспериментальной кривой огибающей амплитуды спинового эха и рассчитанной численно на ЭВМ дало возможность определить с высокой точностью величину и направление локального магнитного поля в соединениях Bi2Ge3O9 и Bi3B5O12, ориентацию градиента кристаллического электрического поля в соединении CdSb.

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

  4. Аналитически и численно исследована амплитуда ядерного квадрупольного спинового эха на ядрах азота 14N с ядерным спином I=1 в случае сильного и слабого радиочастотного поля. Было показано, что сильное радиочастотное поле (~1kЭ) вызывает множественные сигналы спинового эха в радиочастотной катушке спектрометра.

  5. Удалось обобщить уравнения Блоха, описывающие релаксацию ядерного спина, на случай произвольного стационарного гамильтониана, определяющего систему уровней ядерного спина. Как следствие, была разработана алгебра частичных операторов и получены уравнения, описывающие поведение наблюдаемых значений квантовых операторов спина ядра , и в случае ЯКР со спином ядра I = 1. Показано, что в случае ЯКР со спином ядра I = 1, в отличие от ЯМР, количество времён релаксации удваивается. Установлено, что вклады в выражения для обратных времён релаксации, связанные с магнитными и квадрупольными механизмами релаксации, складываются.

^ Цитируемая литература

[1] Абрагам A. Ядерный магнетизм. - М.: Издательство иностранной

литературы, 1963. - 552 с.

[2] Гречишкин В.С. Ядерные квадрупольные взаимодействия в твёрдых телах. -

М.: Наука, 1973. - 264 с.

[3] Bloom M., Hahn E.L., Herzog B. Free Magnetic Induction in Nuclear Quadrupole

Resonance // Phys. Rev. - 1955. - Vol. 97. - P. 1699-1709.

[4] Das T.P., Saha A.K. Electric Quadrupole Interaction and Spin Echoes in Crystals

// Phys. Rev. - 1955. - Vol. 98. - P. 516-524.

[5] Kravchenko E.A., Orlov V.G. Local Magnetic Fields in Some Bismuth

Compounds. A Survey of Experimental Evidences // Z. Naturforsch. A: Phys. Sci.

- 1994. - Vol. 49. - P. 418-424.

[6] Ainbinder N.E., Volgina G.A., Kravchenko E.A., Osipenko A.N., Gippius A.A.,

Fam Suan Hai, Bush A.A. 209Bi NQR powder spectra influenced by local and

applied magnetic fields // Z. Naturforsch. A: Phys. Sci. - 1994. - Vol. 49. - P. 425-

432.

[7] Кравченко Э.А., Фам Суан Хай, Каргин Ю.Ф. Спектры ЯКР 209Bi

соединений Bi2M4O9 (M = Al, Ga), Bi2Ge3O9 // Неорганические материалы. -

1997. - Т. 33. - С. 1001-1003.

[8] Kravchenko E.A., Orlov V.G., Fam Suan Hai, Kargin Yu. F. 209Bi NQR and

Magnetic Properties of Bismuth Oxide-Based Compounds // Z. Naturforsch. A:

Phys. Sci. - 1998. - Vol. 53. - P. 504-513.

[9] Shirley J.H. Solution of the Schrödinger Equation with a Hamiltonian Periodic in

Time // Phys. Rev. - 1965. - Vol. 138. - P. B979- B987.

[10] Barone S.R., Narcowich M.A., Narcowich F.J. Floquet theory and applications //

Phys. Rev. A. - 1977. - Vol. 15. - P. 1109-1125.

[11] Maricq M. Application of average Hamiltonian theory to the NMR of solids //

Phys. Rev. B. - 1982. - Vol. 25. - P. 6622-6632.

[12] Блум К. Теория матрицы плотности и её приложения. - М.: Мир, 1983. -

248 с.

[13] Якубович В.А., Старжинский В.М. Линейные дифференциальные

уравнения с периодическими коэффициентами и их приложения. - М.:

Наука, 1972. - 720 с.

[14] Еругин Н.П. Линейные системы обыкновенных дифференциальных

уравнений с периодическими и квазипериодическими коэффициентами. -

Мн.: АН БССР, 1963. - 272 с.

[15] Провоторов Б.Н., Фельдман Э.Б. Теоретическое изучение динамики спиновой

системы в многоимпульсных экспериментах // Материалы VI всесоюзной школы

и симпозиума по магнитному резонансу. Радиоспектроскопия. - Пермь: ПГУ,

1981. - С. 34-53.

[16] Кокин А.А. Твердотельные квантовые компьютеры на ядерных спинах. -

М.- Ижевск: Институт компьютерных исследований, 2004. - 204 с.

[17] Гардинер К.В. Стохастические методы в естественных науках. - М.: Мир,

1986. - 528 с.

[18] Сликтер Ч. Основы теории магнитного резонанса. - М.: Мир, 1967. - 324 с.

[19] Redfield A.G. On the Theory of Relaxation Processes // IBM Journ. Res.

Develop. - 1957. - Vol. 1. - P. 19-31.

[20] Bloembergen N., Purcell E.M., Pound R.V. Relaxation Effects in Nuclear

Magnetic Resonance Absorption // Phys. Rev. - 1948. - Vol. 73. - P. 679-712.


Основное содержание диссертации изложено в следующих публикациях:

  1. Кравченко Э.A., Орлов В.Г., Шлыков М.П. Магнитные свойства кислородных соединений висмута (III) // Успехи химии. - 2006. - Т. 75. - С. 86-104.

  2. Kravchenko Е.А., Morgunov V.G., Kargin Yu. F., Egorysheva A.V., Orlov V.G., Shlikov M.P. NQR indications of unconventional magnetism in some bismuth-based diamagnets // Applied Magnetic Resonance. - 2004. - Vol. 27. - P. 65-75.

  3. Kravchenko E.A., Orlov V.G., Morgunov V.G., Kargin Yu.F., Egorysheva A.V., Shlikov M.P. Local magnetic fields in some bismuth-based diamagnets. A survey of NQR data // Hyperfine Interactions. - 2004. - Vol. 158. - P. 181-187.

  4. Orlov V.G., Shlikov M.P., Kravchenko E.A., Marenkin S.F., Varnavskii S.A. Magnetism-related properties of CdSb revealed by the Zeeman 121Sb NQR spectra // Hyperfine Interactions. - 2004. - Vol. 159. - P. 173-179.

  5. Kravchenko E.A., Orlov V.G., Morgunov V.G., Shlykov M.P. Zero-field splittings of NQR spectra for bismuth(III) oxy compounds revealed by quadrupole spin echo envelopes // Hyperfine Interactions. - 2007. - Vol. 180. - P. 7-10.

  6. Orlov V.G., Shlykov M.P. Nuclear spin echo model based on Floquet-Lyapunov theory // Hyperfine Interactions. - 2007. - Vol. 180. - P. 11-18.

  7. Шлыков М.П. Стохастическое описание релаксационных явлений в модели квадрупольного спинового эха // Труды 51-й научной конференции МФТИ "Современные проблемы фундаментальных и прикладных наук": Часть X. - М.: МФТИ, 2008. - С. 18-21.

soderzhanie-prakticheskih-zanyatij-rabochaya-programma-po-patologii-ochnoe-obuchenie-po-specialnosti-farmaciya.html
soderzhanie-pravoobladatel-programmi-togou-spo-priborostroietlnij-kolledzh-g-tambov-morshanskoe-shosse-17.html
soderzhanie-predislovie.html
soderzhanie-problemi-i-obosnovanie-ee-resheniya-programmnimi-metodami-doklad-o-sostoyanii-i-ispolzovanii-zemel-v.html
soderzhanie-profilakticheskogo-konsultirovaniya-pri-izbitochnoj-masse-tela-i-ozhirenii.html
soderzhanie-programma-professionalnoj-perepodgotovki-professionalnoj-podgotovki-perepodgotovki-povisheniya-kvalifikacii.html
  • essay.bystrickaya.ru/dosokratim-po-vesne-podskazhite-kakie-viplati-polozheni-materi-odinochke-ya-rabotayu-sejchas-zhdu-vtorogo-rebenka.html
  • spur.bystrickaya.ru/marvanova-z-m-prepodavatel-medicinskoj-mikrobiologii.html
  • pisat.bystrickaya.ru/syomka-podzemnih-kommunikacij.html
  • uchit.bystrickaya.ru/tnk-kak-osnovnoj-subekt-mezhdunarodnogo-marketinga.html
  • bukva.bystrickaya.ru/o-predmetnoj-sfere-marketingovih-issledovanij.html
  • pisat.bystrickaya.ru/tema-prielbruse-nacionalnij-park-kbr-pourochnoe-planirovanie-kursa-kultura-narodov-kbr.html
  • knowledge.bystrickaya.ru/obshestvennoe-mnenie-7.html
  • znanie.bystrickaya.ru/apparatno-programmnij-kompleks-izmereniya-parametrov-potokov-mnogofaznih-sred-n-v-belova-nauchnij-l-a-steshina-k-t-n-docent.html
  • doklad.bystrickaya.ru/uchebnoe-posobie-soderzhanie-k-145-r-145.html
  • textbook.bystrickaya.ru/kalendarno-tematicheskoe-planirovanie-urokov-geometrii-na-200-200-uchebnij-god.html
  • otsenki.bystrickaya.ru/spasibo-metodicheskoe-posobie-moskva-metodicheskoe-posobie-podgotovleno-v-ramkah-obrazovatelnoj-programmi-razrabotka.html
  • doklad.bystrickaya.ru/uchebno-metodicheskij-kompleks-disciplini-specialnost-050400-62-socialno-ekonomicheskoe-obrazovanie-profil-istoriya.html
  • paragraph.bystrickaya.ru/lekciya-xvi-i-xvii-vilyam-dzhems.html
  • ucheba.bystrickaya.ru/pozharov-sravnitelno-nebolshaya-tem-ne-menee-na-segodnya-eto-bolee-dvuh-tisyach-gektarov-na-tushenie-pozhar.html
  • teacher.bystrickaya.ru/glava-4-tajna-grafa-hvostova-mednij-vsadnik-eto-vam-ne-mednij-zmij.html
  • bukva.bystrickaya.ru/o-v-sharipov-ea-pirogov-itarnij-nauchnij-fond-trudi-vserossijskogo-filosofskogo-seminara-molodih-uchenih-im.html
  • write.bystrickaya.ru/glava-3-sovremennie-organizacionno-pravovie-osnovi-sozdaniya-i-deyatelnosti-pii-v-rossii.html
  • urok.bystrickaya.ru/praktikum-viznachiti-socotipi-5-6-chlenv-kolektivu-za-dopomogoyu-testu-1-dodatok-a-viznachiti-nayavnst-kvadr-u-kolektiv-opisht-osoblivost-nayavnih-kvadr-dodatok-g-zrobiti-visnovki-plan.html
  • tasks.bystrickaya.ru/-1-ponyatie-i-vidi-processualnih-srokov-uchebnik-pod-redakciej.html
  • desk.bystrickaya.ru/per-s-it-i-komment-f-m-dvin-287-s-16-l-il-21-sm-m-iskusstvo-1984-avtobiografiya-velikogo-italyanskogo-rezhissera-federiko-fellini-vknige-on-rasskaziv-stranica-12.html
  • notebook.bystrickaya.ru/inostrannij-yazik-pravila-dorozhnogo-dvizheniya-65-russkij-yazik-66.html
  • zadachi.bystrickaya.ru/osnovi-proektuvannya-avtomatizovanih-nformacjnih-sistem.html
  • pisat.bystrickaya.ru/struktura-rabochej-uchebnoj-programmi-disciplini-6.html
  • zanyatie.bystrickaya.ru/razrabotka-tehnologicheskogo-processa-uprochneniya-kulachka-glavnogo-vala-s-ispolzovaniem-lazernogo-izlucheniya-chast-5.html
  • znanie.bystrickaya.ru/62-cel-i-zadachi-duhovno-nravstvennogo-razvitiya-i-vospitaniya-obuchayushihsya-na-stupeni-nachalnogo-obshego-obrazovaniya.html
  • pisat.bystrickaya.ru/telekommunikacionnie-i-informacionnie-tehnologii-v-transportnoj-logistike.html
  • write.bystrickaya.ru/glava-7-rabota-v-centre-yanvar-1943-maj-1945-g-geroi-nevidimogo-fronta.html
  • abstract.bystrickaya.ru/1-specifika-filosofskogo-poznaniya-socialnoj-dejstvitelnosti.html
  • kontrolnaya.bystrickaya.ru/programma-proizvodstvennoj-praktiki-dlya-studentov-4-kursa-specialnosti-mirovaya-ekonomika.html
  • knowledge.bystrickaya.ru/metodicheskie-ukazaniya-i-kontrolnie-zadaniya-dlya-studentov-specialnosti-220301-avtomatizaciya-tehnologicheskih-processov-i-proizvodstv-v-neftegazovoj-otrasli-instituta-distancionnogo-obrazovaniya-tomsk-2008.html
  • thesis.bystrickaya.ru/poslanie-otdohnuvshim-vedomosti-maksim-glikin-vedomosti-30102008-206-str-a2.html
  • uchitel.bystrickaya.ru/proveryaemie-rezultati-korrekcii-effektivnost-ispolnitelnost-i-pribil-vashej-organizacii-mozhno-povisit-za-schet.html
  • prepodavatel.bystrickaya.ru/t-i-konnova-aktualnie-problemi-prepodavaniya-v-nachalnoj-shkole-kiryushkinskie-chteniya.html
  • universitet.bystrickaya.ru/tematichen-plan-osnovavane-na-prvite-evropejski-kolonii-v-severna-amerika-nova-franciya-loyalnite-britanskite-kolonii-i-amerikanskata-vojna-za-nezavisimost.html
  • education.bystrickaya.ru/3-klass-102-chasa-osnovnaya-obrazovatelnaya-programma-nachalnogo-obshego-obrazovaniya-mou-babarikinskaya-sosh.html
  • © bystrickaya.ru
    Мобильный рефератник - для мобильных людей.