Скачать 388.33 Kb.
|
На правах рукописи Шлыков Максим Павлович Компьютерная и математическая модель ядерного спиновОго эха 01.04.02 – теоретическая физика Автореферат диссертации на соискание учёной степени кандидата физико-математических наук Москва 2009 Работа выполнена в Российском научном центре «Курчатовский Институт» Научный руководитель: кандидат физико-математических наук, Орлов Валерий Георгиевич Официальные оппоненты: доктор физико-математических наук, Доброхотов Сергей Юрьевич доктор физико-математических наук, Якубовский Андрей Юрьевич Ведущая организация: Московский физико-технический институт (государственный университет) Защита состоится: 24 июня 2009 г. в 17 ч. 00 мин. на заседании диссертационного совета Д 520.009.03 при РНЦ «Курчатовский институт», 123182 г. Москва, пл. Академика Курчатова, 1. С диссертацией можно ознакомиться в библиотеке РНЦ «Курчатовский институт». Автореферат разослан «_____» мая 2009 г. Ученый секретарь диссертационного совета доктор физико-математических наук А.Л. Барабанов Общая характеристика работы Актуальность работы. Импульсное возбуждение ядерного магнитного резонанса (ЯМР) и ядерного квадрупольного резонанса (ЯКР) широко используется в физике как для изучения внутренних электрических и магнитных полей в кристаллах, так и при исследовании неравновесных состояний ядерных спинов (релаксационных процессов) [1, 2]. Явление спинового эха позволяет наиболее удобно и точно производить измерения времён релаксации, исследовать структуру кристаллов, на явлении спинового эха основана работа современных томографов, ЯКР и ЯМР спектрографов. Огромную роль явление спинового эха играет в физических моделях твёрдотельных квантовых компьютеров. Особое значение при изучении локальных магнитных и электрических полей в кристаллах имеет анализ экспериментальной картины квадрупольного спинового эха в магнитном поле. Если изменять интервал времени между радиочастотными импульсами, то амплитуда спинового эха при наличии постоянного внешнего или локального магнитного поля испытывает “медленные биения”, частота которых определяется зееман-расщеплением в стационарных экспериментах. Основные модели расчёта огибающей амплитуды спинового эха середины 50-x годов, разработанные Блумом, Ханом и Герцогом, а также Дасом и Сахой [3, 4] и положенные в основу работы большинства ЯКР и ЯМР спектрометров обладают рядом приближений. В частности, было выявлено, что результирующие формулы для интенсивности спинового эха, приведенные в этих моделях, справедливы при выполнении неравенства 1tw<<1, где 1=Hrf, - гиромагнитное отношение, Hrf - амплитуда радиочастотного поля, tw – длительность радиочастотного импульса. В то время как в ЯКР спектрометрах заметный сигнал наблюдается лишь при условии 1tw~1. Тем самым, вплоть до настоящего времени не был выполнен детальный анализ области применимости этих формул и не сформулирована математически корректная процедура моделирования процесса квадрупольного спинового эха, которая позволяла бы извлекать информацию о физических свойствах вещества – величину и ориентацию локального магнитного поля – на основании сравнения экспериментальных данных с результатами модельных расчетов. Актуальность разработки численной модели продиктована необходимостью не только качественно, но и количественно интерпретировать экспериментальные данные по квадрупольному спиновому эху, разработать универсальную схему, применимую при вычислении амплитуды спинового эха в случае целого и полуцелого спина ядра, сильного и слабого радиочастотного поля и учитывающую эффекты релаксации. Например, при анализе спектров ЯКР на ядрах 209Bi установлена уникально высокая чувствительность электронных характеристик к воздействию слабых внешних магнитных полей [5, 6, 7, 8], указывающая на сильную взаимосвязь магнитной и электронной подсистем соединений. Это говорит о принципиальной возможности изменять функциональные свойства соответствующих материалов, воздействуя на них слабыми магнитными полями. Сравнение экспериментальной кривой огибающей амплитуды спинового эха и рассчитанной численно на ЭВМ дало бы возможность определить величину и направление локального магнитного поля, ориентацию градиента электрического поля в этих соединениях. Кроме того, поскольку у некоторых ядер маленький квадрупольный момент (например, у ядра азота по сравнению с ядром висмута), частоты квадрупольного резонанса и интенсивности сигналов спинового эха малы и особую актуальность приобретает задача разработки ЯКР спектрометров, более эффективных по сравнению с существующими. Актуальность задачи измерения и анализа огибающей амплитуды ядерного квадрупольного спинового эха в соединениях азота обусловлена необходимостью разработки эффективных средств контроля и безопасности на транспорте. Решение задачи о расчёте амплитуды спинового эха в случае малых частот основного стационарного гамильтониана, определяющего систему уровней ядерного спина, и сильного радиочастотного поля требует усовершенствование существующих аналитических моделей для описания квантовых систем с периодически зависящим от времени гамильтонианом [9, 10, 11]. Большинство таких моделей обладают значительной громоздкостью и ограниченной применимостью результирующих формул в области близкой к резонансу. Цель работы. Первой целью работы было создание компьютерной модели спинового эха, применимой для моделирования амплитуды спинового эха и извлечения информации о ранее неизвестных физических свойствах вещества – величины и ориентации локального магнитного поля, направления градиента кристаллического электрического поля в образце. Второй целью являлся анализ условий формирования спинового эха, усовершенствование аналитических моделей, применяемых для описания динамики ядерного спина. Для выполнения поставленных целей необходимо было решить следующие основные задачи:
Научная новизна. Впервые разработан комплекс программ, позволяющий не только качественно, но и количественно интерпретировать экспериментальные данные по квадрупольному спиновому эху, моделировать временную зависимость амплитуды спинового эха. В рамках метода матрицы плотности, используемого для моделирования ядерного квадрупольного спинового эха, был предложен подход, основанный на теореме Флоке-Ляпунова, к решению систем дифференциальных уравнений, описывающих временную зависимость оператора эволюции. Отличительной особенностью разработанного метода по сравнению с ранее разработанными схемами [9, 10, 11], использующими теорему Флоке-Ляпунова, является его простота использования и применимость в области близкой к резонансу, в случае, когда уровни спина ядра, определяемые основным стационарным гамильтонианом, не являются эквидистантными. С помощью полученного метода было объяснено явление множественного спинового эха: моделирование амплитуды спинового эха для двухимпульсной методики измерения ЯКР в соединениях азота 14N с ядерным спином I = 1 предсказывает возможность обнаружения множественных сигналов индукции вблизи “классических” сигналов индукции после первого и второго импульсов и сигнала спинового эха в случае сильного радиочастотного поля. В результате произведённых аналитических расчётов было показано существование двух режимов поведения ядерного спина в радиочастотном поле. В случае слабого радиочастотного поля (частота 1=Hrf мала по сравнению с характерной частотой перехода между уровнями основного стационарного гамильтониана) установлена возможность разделения поведения ядерного спина на быстрые (с частотой, определяемой переходами между уровнями основного гамильтониана) и медленные (обусловленные поворотом спина вокруг направления радиочастотного поля на угол 1tw) движения. В случае сильного радиочастотного поля такое разделение невозможно, амплитуда спинового эха становится сложной функцией резонансной частоты, определяемой уровнями основного гамильтониана. Было показано, что формулы Даса и Саха, применяемые при расчётах сигнала спинового эха в современных спектрометрах, являются справедливыми в нулевом приближении по радиочастотному полю при условии выполнения соотношения на длительность радиочастотного импульса tw: . Установлено, что формулы Даса и Саха могут быть получены предельным переходом из общих формул, основанных на теореме Флоке. Таким образом, удалось определить область применимости результирующих формул существующих моделей Блума-Хана-Герцога и Даса-Саха и получить обобщённые формулы для расчёта амплитуды спинового эха. Разработав алгебру частичных спиновых операторов, удалось решить квантомеханические релаксационные уравнения Блоха в случае ЯКР со спином ядра I = 1. По аналогии со случаем ЯМР, в случае ЯКР со спином ядра I = 1 были получены выражения для времён спин-спиновой и спин-решёточной релаксации, в которые входят основные параметры задачи – времена корреляции, частоты переходов основного гамильтониана. Было показано, что в случае ЯКР со спином ядра I = 1, в отличие от ЯМР, количество времён релаксации удваивается. Установлено, что вклады в выражения для обратных времён релаксации, связанные с магнитными и квадрупольными механизмами релаксации, складываются. Практическая ценность работы. Возможности расчётной схемы были непосредственно использованы для моделирования “биений” огибающей амплитуды спинового эха во внешнем или локальном магнитном поле для определения ориентации градиента электрического поля на монокристалле CdSb, направления и величины локального магнитного поля в соединениях Bi2Ge3O9 и Bi3B5O12. В результате проведенных теоретических и экспериментальных исследований разработана компьютерная и аналитическая модель спинового эха, включающая в себя исчерпывающий набор составляющих: численный расчёт, анализ динамики ядерного спина и учёт эффектов релаксации при помощи стохастических моделей описания релаксационных явлений. Основные положения, выносимые на защиту.
Публикации и апробация работы. Результаты данной работы докладывались на научном семинаре лаборатории теоретической физики и лаборатории многочастичных систем Института общей и ядерной физики (ИОЯФ) РНЦ "Курчатовский институт" под руководством проф. В.Г. Вакса, на научном семинаре Института радиотехники и электроники им. В.А. Котельникова (ИРЭ) РАН под руководством проф. В.А. Ацаркина 17 октября 2006 г. Результаты данной работы докладывались на 13-й Международной конференции по сверхтонким взаимодействиям и 17-м Международном симпозиуме по ядерным квадрупольным взаимодействиям (Бонн, Германия, 2004 г.), 14-й Международной конференции по сверхтонким взаимодействиям и 18-м Международном симпозиуме по ядерным квадрупольным взаимодействиям (Бразилия, 2007 г.), а так же представлены в качестве докладов “Квадрупольное спиновое эхо в магнитном поле” на 1-й Курчатовской молодёжной школе в 2003 г., “Ядерное квадрупольное спиновое эхо на ядрах азота” на 3-й Курчатовской молодёжной школе в 2005 г., "Стохастическое описание релаксационных явлений в модели квадрупольного спинового эха" на 6 й Курчатовской молодёжной школе в 2008 г. (доклад отмечен дипломом лучшей работы в секции фундаментальных исследований), "Стохастическое описание релаксационных явлений в модели квадрупольного спинового эха" на 51 й Научной конференции МФТИ в 2008 г. По результатам диссертации опубликовано 7 статей, из них 6 статей опубликованы в рецензируемых российских и зарубежных журналах. Личный вклад автора в работу состоит в непосредственном и активном участии в постановке всех задач, проведении исследований и интерпретации результатов. Лично автором разработаны программы моделирования огибающей амплитуды квадрупольного спинового эха в магнитном поле и проведены соответствующие численные расчеты. Автором проведен анализ параметров существующих расчётных и аналитических моделей спинового эха. Автор приложил значительные усилия для улучшения аналитических методов расчёта амплитуды спинового эха, использующих теорему Флоке. При учёте релаксации ядерного спина в моделях спинового эха, автор самостоятельно разработал алгебру частичных спиновых операторов, необходимую для решения квантовых релаксационных уравнений Блоха в случае ЯКР и спина ядра I = 1, и получил выражения для времён спин-спиновой и спин-решёточной релаксации. Структура и объем работы. Диссертационная работа состоит из введения, пяти глав, заключения, списка литературы и трёх приложений. Объём диссертации составляет 143 страницы, в том числе 12 рисунков, 2 таблицы и библиографический список из 51 наименования. |
Программа по формированию навыков безопасного поведения на дорогах... Знать понятие о задаче в общем виде, форматизация задачи, математическая и компьютерная модель | Программа по формированию навыков безопасного поведения на дорогах... Формирование понятий: описательная информационная модель, формализованная модель, компьютерная модель, компьютерный эксперимент,... | ||
Урок по теме "Математическая культура". 7-й класс Вступительное слово учителя (сопровождение компьютерная презентация см. Приложение) | Реферат по курсу: Товарная политика предприятия на тему: Товарный... Товарный ассортимент. Товарная номенклатура. Управление ассортиментом. Жизненный цикл товара. Математическая модель планирования... | ||
Новая компьютерная модель для диагностики сотрясений мозга среди спортсменов Присутствовали: 11 человек из 15 членов Совета. От управы Бутырского района: Беляев А. А., Серёгин И. В., Старкова О. Ф., Волгин... | Программа по формированию навыков безопасного поведения на дорогах... ... | ||
Рефераты №11-12 2013 г Разработана многопараметрическая математическая модель износа ручья. Рассмотрен способ повышения долговечности канатоведущего шкива... | Математическое исследование баллистического движения Аннотация. В докладе проводится изучение баллистического движения методами математики, с использованием программного пакета «maple».... | ||
V. V. Poroshin, A. V. Onanko, D. G. Bogomolov трехмерная математическая... Порядку оказания медицинской помощи по профилю «хирургия (трансплантация органов и (или) тканей человека)», утвержденному приказом... | Я дрегля Нина студентка s-13 прочитала 150 страниц из книги «Семейная... Ют сегодня виды консультативной психологической помощи семье чрезвычайно разнообразны. В соответствии с ориентированностьюи характером... | ||
Доклад Секретариата: «Перспективы взаимодействия Израиля и других... Большим успехом стало создание зон, свободных от ядерного оружия («безъядерных зон») в Латинской Америке и Карибском бассейне, южной... | Математическая модель термоупругого деформирования слоистых композитных оболочек и пластин Закона Камчатского края от 18. 09. 2008 №122 «О дополнительных гарантиях и дополнительных видах социальной поддержки детей-сирот... | ||
Модель оценки альтернатив управления слабоструктурированными динамическими ситуациями 1 Рассмотрена интегрированная нечеткая система поддержки принятия решений в слабоструктурированных динамических ситуациях, включающая... | Программа по формированию навыков безопасного поведения на дорогах... Математическая экономика – это прикладная математическая дисциплина, в которой изучаются конкретные количественные отношения экономических... | ||
Математическая модель кристаллизации переохлажденных капель водных... «Речевой этикет» разработана на основе авторской программы по русскому языку для общеобразовательных учреждений. 5-11 классы: (автор-составитель... | Методические рекомендации к выполнению проекта по технологии 1 Проект это самостоятельная творческая завершенная работа. Она обычно состоит из двух частей: теоретической и практической. В качестве... |