WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

Казанский (Приволжский) федеральный университет

Физический факультет

Кафедра квантовой электроники и радиоспектроскопии

Использование программного модуля

EasySpin в анализе спектров

магнитного резонанса.

Часть 1. Стационарный ЭПР

Учебно-методическое пособие

для студентов и аспирантов физического факультета

Казань 2010

УДК 537.635, 577.334

ББК 22, 28

Печатается по решению методической комиссии физического факультета Протокол №8 от 10 июня 2010 года заседания кафедры квантовой электроники и радиоспектроскопии Протокол №13 от 4 июня 2010 года Авторы:

канд. физ. мат. наук, доцент Г.В. Мамин канд. физ. мат. наук, доцент С.Б. Орлинский докт. физ. мат. наук, профессор Н.И. Силкин асп. И.Н. Субачева канд. физ. мат. наук, доцент Р.В. Юсупов Рецензент:

докт. физ. мат. наук, профессор КИББ РАН А. В. Анисимов Содержание Введение

1.

Установка модуля EasySpin

2.

Возможности модуля Easyspin

3.

Работа с модулем EasySpin

4.

Определение спиновой системы (переменная Sys)

4.1.

Углы Эйлера

4.2.

Задание параметров уширения

4.4.

Задание условий эксперимента (структура Exp)

4.5.

Структура Opt

4.6.

Пример использования модуля EasySpin: ион Yb3+ в кристалле BaF2...... 5.

Список литературы

6.

1. Введение Спектроскопические методы магнитного резонанса на сегодняшний день используются в очень широкой области научных исследований. Это и физика, и химия, и биология, и стремительно развивающаяся сфера нанотехнологий.

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

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

распространяемый модуль «EasySpin», представляющий собой набор подпрограмм, работающих в пакете Matlab. Данный модуль написан и развивается д-ром Стефаном Столлом (Dr. Stefan Stoll).

Настоящее методическое пособие предназначено для студентов старших курсов и аспирантов, осваивающих и выполняющих исследования методом стационарного электронного парамагнитного резонанса (ЭПР). В рамках пособия будут рассмотрены основные моменты и некоторые важные особенности расчета спектров ЭПР в модуле Easyspin. Выбор данного продукта обусловлен следующим рядом причин: общедоступность данного программного продукта, его мультиплатформенность (EasySpin работает на платформах Windows, Linux и Mac), наличие у Казанского университета официальной лицензии на пользование пакетом Matlab. Важным фактором явилось и то, что идеология разработчиков модуля EasySpin позволяет пользователю легко создавать собственные программы, используя подпрограммы модуля вкупе с широкими возможностями пакета Matlab.

Считаем необходимым отметить, что данное пособие предполагает вполне определенный уровень подготовки читателя. Он должен знать, что такое электронный парамагнитный резонанс, быть знаком с физическими основами ЭПР спектроскопии, особенностями используемых в экспериментах конфигураций спектрометров. Кроме того, пользователь должен осознавать ограничения, накладываемые приближением спинового гамильтониана, реализованного в модуле EasySpin, дабы не требовать от продукта «невозможного».

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

Столла с соавторами [1, 2].

Описание работы с модулем EasySpin начнем с рекомендаций по его установке. Предполагается, что работа ведется на платформе Microsoft Windows, программный пакет Matlab установлен и работоспособен.

Модуль EasySpin доступен для скачивания в сети Интернет по адресу http://www.easyspin.org. Модуль распространяется в виде ZIP-архива, после распаковки соответствующая файловая структура содержит папки documentation (документация), easyspin (папка подпрограмм) и examples (примеры). Для получения доступа к подпрограммам модуля EasySpin в основном окне пакета Matlab необходимо выбрать File Set Path..., и в открывшемся окне выбрать местонахождение подпапки easyspin из трех распакованных папок и нажать кнопку Add Folder, затем Save и Close.

Далее для компиляции подпрограмм в командном окне Matlab необходимо набрать easyspin (для версии easyspin-3.1.1) или easyspincompile (для версии 3.0.0) и нажать Enter. При предложении выбрать компилятор С для Windows 32 бит следует выбирать Lcc (входит в Matlab). С некоторыми сложностями можно столкнуться при работе на платформе Windows 64 бит.

Соответствующего компилятора в пакете Matlab нет, его необходимо устанавливать и настраивать дополнительно. Однако, пользователь может легко найти соответствующую информацию как на сайте компании MathWorks – производителя пакета Matlab, так и на иных страницах сети Интернет.

После компиляции подпрограмм модуля EasySpin файлы помощи будут доступны также через общую помощь пакета Matlab.

предназначенных для расчета спектров, соответствующих следующим основным случаям:

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

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

Функция pepper: стационарный ЭПР в твердом теле Позволяет производить расчеты спектров для порошков и кристаллических образцов с учетом пространственных групп симметрии с любым числом электронных и ядерных спинов. В расчетах учитываются все взаимодействия, в том числе операторы высокого порядка и ядерное квадрупольное взаимодействие. Есть возможность выбора метода вычислений: путем диагонализации матрицы либо с использованием теории возмущений второго порядка. Включены различные механизмы уширения линий: за счет разброса значений g, D и A, а также за счет неразрешенных сверхтонких расщеплений.

Можно задавать системы с неравновесной заселенностью уровней, учитывать параллельную и перпендикулярную ориентацию компонент магнитного поля.

Функция salt: двойной электронно-ядерный резонанс (ДЭЯР, или ENDOR) в твердом теле кристаллических образцов при произвольном числе ядер. Учитывает сверхтонкие взаимодействия. Имеет встроенную систему выбора ориентаций.

Позволяет выбирать метод вычисления: путем диагонализация или по теории возмущений.

Функция saffron: импульсный ЭПР/ДЭЯР в твердом теле Рассчитывает спектры систем с произвольным числом ядер для двух- и трех импульсных последовательностей. Есть возможность задавать последовательности самостоятельно. Система может обладать большим спином. Включен учет квадрупольных взаимодействий.

Функция esfit: аппроксимация методом наименьших квадратов.

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

Подпрограммы модуля EasySpin предоставляют достаточно широкие возможности и, в частности, позволяют учитывать такие параметры «эксперимента» при расчете спектров, как температура, направление магнитного поля, симметрия образца, конфигурация используемого резонатора, различные механизмы уширения линий. Есть возможность рассчитывать угловые зависимости, выделяя нужные переходы. При расчетах можно выбирать способ определения уровней энергии: путем диагонализации матрицы или же во втором порядке теории возмущений, что значительно экономит время в случае большого объема вычислений. Кроме того, у модуля EasySpin есть ряд дополнительных возможностей. К ним относятся:

краткий справочник по физике ЭПР, включающий в себя информацию о спиновых операторах, операторах Стивенса, спиновом гамильтониане и угловых моментах;

различные утилиты, например, для расчета ширины линий, углов Эйлера, спектрометров, база данных ядерных изотопов, утилита для пересчета магнитное поле/частота/g-фактор;

анализ данных: интегрирование RC цепочкой, модуляция поля, Фурьепреобразование, коррекция базовой линии;

расчет диаграмм энергетических уровней.

Работу с модулем EasySpin мы рассмотрим на примере функции pepper, работа с остальными функциями во многом аналогична. Функция pepper позволяет рассчитывать спектры стационарного ЭПР монокристаллических и порошковых образцов. Общий вид использования функции:

pepper(Sys,Exp);

pepper(Sys,Exp,Opt);

spec = pepper(...);

[B,spec] = pepper(...);

[B,spec,trans] = pepper(...);

В первых двух случаях выходные параметры отсутствуют, и функция pepper симулирует спектр ЭПР с заданными условиями. Выходная переменная spec является массивом ординат спектра (интенсивность сигнала). Переменная B представляет собой массив значений магнитного поля в миллитесла, в которых вычисляются величины сигнала ЭПР; trans содержит список пар номеров уровней, переходы между которыми представлены в спектре. Номера соответствуют уровням энергии, полученным при диагонализации матрицы, в возрастающем порядке, то есть уровень 1 будет обладать наименьшей энергией.

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

Расчет спектра при вызове любой из вышеперечисленных функций требует определения как самого парамагнитного центра (спиновой системы), так и «условий эксперимента». Рассмотрим эти два важнейших шага.

4.1. Определение спиновой системы (переменная Sys) Для этой цели служит входной параметр Sys. Он содержит параметры спинового гамильтониана и характеристику уширения спектральных линий. Sys представляет собой тип данных Matlab под названием structure. Такая структура содержит в себе набор полей, которые могут принимать различные значения: числовые, символьные или типа массив. Например, для структуры, определяющей спиновую систему, полями переменной Sys будут являться спин, g-тензор, тензор сверхтонкого взаимодействия с собственным ядром, типы ядер ближайшего окружения, их характеристики и т.д. Существует два способа определения данных типа structure в Matlab. В первом случае вся структура описывается одной строкой Sys = struct('S',1/2,'g',[1.9 2.0 2.1]).

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

Sys.g = [1.9 2.0 1];

С помощью такой записи задается парамагнитная система со спином и орторомбическим g-фактором.

Рассмотрим поля переменной Sys и способы их определения подробнее.

Здесь необходимо отметить несколько общих моментов, а именно:

1. Энергетические параметры в спиновом гамильтониане необходимо задавать в МГц. Все углы определяются в радианах, значения магнитного поля в мТ.

2. Если значение поля не задано, то в программе используется значение по умолчанию.

3. В пакете Matlab квадратная матрица вида может быть записана следующим способом:

A = [a11 a21 a31; a12 a22 a32; a13 a23 a33];

То есть, сначала записывается первый ряд, потом второй и так далее. При задании спиновой системы из нескольких спинов или ядер они записываются как члены матрицы-столбца, т. е. ряда.

Поле S – электронный спин системы. Значение по умолчанию 1/2.

Запись вида определяет систему, состоящую из одного парамагнитного центра со спином 3/2; запись задает систему из трех центров (ионов, радикалов), два из которых имеют спин 1 и один – спин.

Поле Nucs – задает ядра, присутствующие в системе. По умолчанию взаимодействующих с электронным спином ядер нет.

Запись вида Sys.Nucs = `Cu`;

используется для системы, содержащей природную смесь изотопов 63Cu и 65Cu (природная распространенность изотопов и спины ядер при расчетах амплитуд линий спектра учитываются автоматически);

Примеры:

Sys.Nucs = '59Co'; система содержит только один изотоп 59Co;

Sys.Nucs = '59Co,14N,14N'; система содержит три ядра (один Поле Nucs задает программе спины соответствующих ядер. Иногда удобнее использовать функцию nucspinadd (см. ниже).

Поле g – задает g-тензор системы. Значение этого поля по умолчанию равно 2.0023, что соответствует изотропному g-фактору свободного электрона.

В общем случае для каждого спина можно задавать матрицу 3х3, содержащую 9 компонент g-тензора. Однако чаще используется запись, содержащая только главные значения тензора. Общий вид такой записи:

Sys.g = [x1 y1 z1; x2 y2 z2;…];

здесь каждый ряд из трех значений соответствует одному спину в последовательности, заданной в поле S; запись вида Sys.g = [2 2.05 2.3];

соответствует орторомбическиму g-фактору для системы из одного спина;

Sys.g = [2 2.1 2.3; 1.9 1.95 2.01];

- системе из двух спинов с разными g-факторами.

Есть возможность не выписывать все значения g-фактора в случае относительно высокой симметрии, например, запись Sys.g = [2.25 2.03];

задает один спин с анизотропным g-фактором при наличии оси симметрии, первое значение соответствует направлению, перпендикулярному оси симметрии, второе – параллельному оси [x z];

Sys.g = 2.005; – изотропный g-фактор для одного спина.

Поле gpa – каждый ряд этого массива содержит три угла Эйлера в радианах, описывающие переход от собственной системы координат gтензора, в которой он имеет диагональный вид, к молекулярной системе координат; по умолчанию предполагается, что главные направления осей g-тензора и оси молекулярной системы координат совпадают.

Пример описания Sys.gpa = [0 10 0]*pi/180; – для одного спина;

Далее для справки схематично показаны углы Эйлера.

Три угла Эйлера (,, ) описывают переход от одной декартовой системы координат xyz к другой XYZ, при этом начало отсчета остается на месте.

Переход осуществляется в три этапа:

1) Поворот системы xyz вокруг оси z на угол, в результате получаем систему x'y'z';

2) Поворот системы x'y'z' вокруг оси y' на угол, что приводит к системе 3) Поворот системы x''y''z'' вокруг оси z'' на угол, получаем искомую В программном модуле EasySpin углы Эйлера записываются именно в такой последовательности [ ]. Для всех полей типа *pa (где звездочка обозначает соответствующий тензор) система xyz будет отвечать собственной системе координат тензора, где он имеет диагональный вид, а система координат XYZ – молекулярной системе координат.

Поле A – определяет тензор сверхтонкого взаимодействия, вводится значение A/h в МГц, где h – постоянная Планка; по умолчанию все компоненты этого тензора равны нулю.

Может быть задано по 9 значений тензора для каждого ядра из поля Nucs, обычно указываются только главные значения. Общий вид записи здесь первый ряд из трех главных значений соответствует первому ядру, записанному в поле Nucs, второй – второму ядру;

по аналогии с g-фактором при наличии аксиальной симметрии или в изотропном случае запись может быть сокращена, например, Sys.A = 34; что соответствует [34 34 34].

Поле Apa – каждый ряд матрицы содержит три угла Эйлера, определяющие ориентацию главных осей тензора константы сверхтонкой структуры относительно осей молекулярной сиcтемы координат, записывается аналогично полю gpa.

Поле Q – массив, содержащий главные значения тензора квадрупольного взаимодействия в виде Q/h в МГц, где h – постоянная Планка; по умолчанию полагается нулевым;

в случае аксиальной симметрии вводится одно значения для каждого ядра, например, Sys.Q = [0.7, 1.2]; – для двух ядер значения e2qQ/h;

в случае орторомбического тензора Q для каждого ядра задается два значения e2qQ/h и :

Sys.Q = [1.2 0.29; 0.1 0]; – для двух разных ядер;

в любом случае можно задавать три главных значения Q-тензора для каждого ядра.

Поле Qpa – здесь задаются углы Эйлера, определяющие относительную ориентацию главных осей тензора квадрупольного взаимодействия и молекулярной системы координат, записывается аналогично полю gpa для каждого ядра, например, Sys.Qpa = pi/180*[0 30 0; 0 -30 0]; – для двух различных Поле D – определяет главные значения тензора параметров тонкого расщепления в МГц; по умолчанию тензор предполагается нулевым.

Как и для других полей (A, g), каждый ряд этого массива соответствует одному из спинов, указанных в поле S. Любой ряд состоит из трех, двух или одного значения. Если указано два числа, то считается, что первое – это значение параметра D, второе – параметра Е. Если указывается лишь одно число, то считается, что Е равно нулю. В случае задания трех значений они представляют собой главные значения тензора D. Например, запись Sys.D = [120 50]; соответствует D = 120 МГц и Е = 50 Мгц.

Поле Dpa – данный массив содержит углы Эйлера, определяющие ориентацию главных осей тензора D относительно осей молекулярной сиcтемы координат, записывается аналогично полю gpa.

Поле aF – задает параметры расщепления в кристаллическом поле кубической симметрии в МГц, значения по умолчанию – нули.

Здесь первое значение соответствует параметру a, второе – параметру F; эти параметры считаются определенными в молекулярной системе координат, оси которой совпадают с осями симметрии четвертого порядка в случае кубической симметрии.

Поле ee – представляет собой массив, который содержит главные значения матрицы электрон-электронного взаимодействия, каждый ряд соответствует диагонали матрицы (в ее диагональном виде); вводятся значения, приведенные к постоянной Планка; по умолчанию массив предполагается нулевым.

Если n – число спинов, то число рядов равно R = n(n-1)/2; то есть для 4 спинов нужно задать 6 значений в следующей последовательности: 1-2, 1-3, 1-4, 2-3, 2где цифры – номера спинов, как они заданы в поле S; например, при учете только изотропного взаимодействия возможна запись Sys.S = [1/2 1/2 1/2 1/2]; –четыре электрона, Sys.ee = [10 0 0 5 -10 20]; – шесть значений для взаимодействий;

кроме того, можно записывать матрицу целиком для каждой пары электронов или один ее соответствующий ряд.

Поле eepa – каждый ряд этой матрицы соответствует ряду матрицы в поле e и содержит три угла Эйлера, описывающие переход от собственной к молекулярной системе координат, записывается аналогично подобным полям для других параметров (gpa, Apa) Задавать параметры ядер в системе можно как заполняя поле Nucs в структуре Sys, так и используя встроенную функцию nucspinadd.

Общий вид записи функции NewSys = nucspinadd(Sys,Nuc,A) NewSys = nucspinadd(Sys,Nuc,A,Apa) NewSys = nucspinadd(Sys,Nuc,A,Apa,Q) NewSys = nucspinadd(Sys,Nuc,A,Apa,Q,Qpa) Функция, заданная таким образом, к спиновой системе, описанной в структуре Sys, добавляет ядро Nuc. A определяет параметры сверхтонкой структуры, Q – тензор квадрупольного взаимодействия, параметры Apa и Qpa содержат углы Эйлера, определяющие переход от главных осей соответствующих тензоров к молекулярной системе координат.

Sys = struct('S',1/2,'g',[2 2 2.2]);

Sys = nucspinadd(Sys,'63Cu',[50 50 520]);

Sys = nucspinadd(Sys,'14N',[3 3 9],[],[-1 -1 2]) Такая запись соответствует системе из одного электронного спина и двух ядер 63Cu и 14N. Иногда использование функции nucspinadd является более удобным, поскольку для добавления ядра в систему или удаления его из системы достаточно изменить или удалить всего одну строку из программы.

В EasySpin существует два способа описания уширения линий ЭПР: для изотропного и анизотропного случая.

В случае изотропной ширины каждому вычисленному значению резонансного поля ставится в соответствие линия в спектре с формой Гаусса или Лоренца и с заданной пользователем шириной. Такой метод определения является простейшим, он не предполагает никаких физических оснований уширения и носит название «феноменологический». Параметры уширения задаются в структуре Sys с помощью заполнения соответствующих полей.

Поле lwpp – определяет ширину линии peak-to-peak в миллитесла.

Представляет собой характеристику ширины линии, равную горизонтальному расстоянию между максимумом и минимумом первой производной линии поглощения; если указывается одно значение, то форма линии считается гауссовой, в общем случае ряд из двух значений [Gauss Lorentz]. Например, Sys.lwpp = 10; - гауссова форма, в мТ Sys.lwpp = [0 12]; - лоренцева форма, в мТ Sys.lwpp = [10 12]; - форма линии типа Фойхт, в мТ Поле lw – полная ширина на половине высоты линии в миллитесла (иногда обозначается как FWHM – full width at half maximum).

Следует отметить, что значения по умолчанию для полей, ответственных за ширину линии ЭПР в симулируемом спектре не заданы, в любом случае одно из полей должно быть заполнено.

Часто бывает так, что уширение связано с неразрешенной сверхтонкой структурой или деформациями, присутствующими в образце, что выражается в некотором распределении параметров спинового гамильтониана между парамагнитными центрами. Такие механизмы могут в результате приводить к анизотропии ширины линии. Для их описания в EasySpin служат специальные поля структуры Sys.

Для задания параметров уширения, обусловленного неразрешенной сверхтонкой структурой, используется поле HStrain. Если механизм уширения связан с распределением значений параметров гамильтониана по объему, необходимо задать поля gStrain, AStrain или DStrain. Следует отметить, что использование этих трех последних полей наиболее корректно при описании системы, состоящей из одного электронного спина, в случае системы из нескольких спинов в расчетах предполагается, что эти поля соответствуют только первому заданному спину.

Поле HStrain – определяет полную ширину на половине высоты линии в МГц (FWHM), обусловленную неразрешенной сверхтонкой структурой;

представляет собой массив из трех значений, каждое из которых соответствует ширине гауссовой линии в МГц для трех главных направлений молекулярной системы координат x, y и z.

Пример записи:

Sys.HStrain = [10 10 50]; - описывает гауссово уширение в 10 МГц для направлений x и y и большее уширение в 50 МГц для выделенного направления z.

Результирующая ширина линии для произвольного направления постоянного магнитного поля, заданного направляющими косинусами (nx, ny, nz), определяются в поле HStrain.

Поле gStrain – определяет полную ширину на полувысоте для гауссовых распределений каждого из трех главных значений g-тензора первого электронного спина. Распределения полагаются полностью некоррелированными. Общий вид записи:

Sys.gStrain = [FWHM_gx FWHM_gy FWHM_gz];

Поле AStrain – представляет собой вектор, три компоненты которого соответствуют ширинам гауссовых распределений главных значений тензора A(x, y, z) в МГц для первого заданного ядра. Предполагается, что эти распределения независимы. Вид записи аналогичен полю gStrain.

Поле DStrain – определяет ширины гауссовых распределений параметров тонкого расщепления D и E в МГц, если значение для параметра Е не задано, оно предполагается равным нулю; эти распределения также независимы. Общий вид записи Sys.DStrain = FWHM_D; в случае E = 0, или Sys.DStrain = [FWHM_D FWHM_E];

4.5. Задание условий эксперимента (структура Exp) «экспериментальных условий» в структуре Exp. Записывается она аналогично структуре Sys. Стоит отметить, что последовательность записи на работу программы не влияет, т.е. можно сначала определить экспериментальные условия, а потом спиновую систему, или даже в процессе записи одной структуры определить часть параметров другой (в случае записи построчно).

Рассмотрим основные поля структуры, задающей условия эксперимента.

Поле mwFreq – задает частоту спектрометра в ГГц, является обязательным для заполнения. Пример записи:

Exp.mwFreq = 9.5.

Диапазон магнитных полей, в котором рассчитывается спектр ЭПР, в модуле EasySpin может быть задан двумя способами: через CenterSweep или через Range, если заданы оба, то программой используется CenterSweep.

Поле CenterSweep – cодержит два значения, первое – центральное значение магнитного поля в мТ, второе – ширина развертки также в мТ, например, магнитного поля от 300-200/2 = 200 мТ до 300+200/2 = 400 мТ.

Поле Range – содержит границы диапазона значений магнитного поля в мТ, первое значение – нижняя граница, второе – наибольшее значение магнитного поля; запись Exp.Range = [200 400]; – задает тот же самый диапазон, что и предыдущий пример.

магнитного поля B0, может принимать два значения: `perpendicular` и `parallel`; по умолчания предполагается, что эти компоненты перпендикулярны.

Exp.Detection = `perpendicular`; - B1 B0;

Exp.Detection = `parallel`; - B1 B0.

Поле nPoints – задает количество точек по оси абсцисс, на которые разбивается диапазон магнитных полей и в которых вычисляются значения спектра. Значение по умолчанию – 1024.

Поле Harmonic – определяет, в каком виде будет представлен искомый спектр. Может принимать только три значения: 0, 1 и 2. По умолчанию задано значение 1. Пример Exp.Harmonic = 0; – соответствует симуляции спектра поглощения;

1 и 2 отвечают первой и второй производной спектра поглощения соответственно.

Поле Ordering – определяет распределение ориентаций парамагнитных молекул в образце; если значение не задано или равно 0, то распределение считается равномерным, т. е. все направления ориентации имеют равную вероятность; если задано значение, то распределение вычисляется по формуле P( ) exp( (3cos 2 1) / 2), где – угол между осью z и направлением внешнего магнитного поля, а задается в поле Ordering.

Поле Temperature – это поле определяет населенность уровней энергии через значение температуры или напрямую; если не задано определенное значение, то все уровни считаются заселенными одинаково.

В случае теплового равновесия - если задан скаляр, то он определяет температуру системы в кельвинах, населенность рассчитывается согласно закону Больцмана; для неравновесного состояния системы, содержащей N уровней энергии, должно быть определено N значений поля Temperature, каждое из которых соответствует населенности одного уровня, например, Exp.Temperature = [0.85 0.95 1.2]; – для системы со спином 1, обладающей тремя уровнями энергии, первое значение отвечает уровню с наименьшей энергией.

Поле Crystal Symmetry – определяет симметрию кристалла, может задаваться номером пространственной группы (от 1 до 230) или символом пространственной группы, также можно указывать символ точечной группы симметрии. Значение по умолчанию 1 (Р1).

Поле Orientations – определяет направление внешнего магнитного Представляет собой массив, каждый ряд которого состоит из трех [ ] (или двух [ ]) углов Эйлера, описывающих переход от молекулярной системы координат к системе координат внешнего магнитного поля. Здесь первый угол – это угол между осью x и проекцией на плоскость xy вектора постоянного магнитного поля; второй угол – это угол между осью z спиновой системы координат и направлением внешнего магнитного поля; третий угол необязателен для определения, его значение не влияет на величину резонансного поля, но определяет интенсивность линий. В системе координат магнитного поля считается, что поле B0 направлено вдоль оси z, а поле B1 – вдоль оси x (когда они перпендикулярны). Примеры:

Exp.Orientations = [0;0;0]; – ось кристалла z совпадает с направлением поля B0;

Exp.Orientations = [0;pi/2;0]; - ось z перпендикулярна B0;

Если значение поля Orientations не задано или задано пустым [], то спектр рассчитывается не для кристалла, а для порошка.

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

Поле Verbosity – определяет количество информации, выводимой на экран помимо рассчитанного спектра. По умолчанию задано значение 0.

0 – только спектр;

1 – выводит основную информацию;

2 – более детальную.

Поле Output – значение этого поля определяет, в каком виде выводится спектр и выходная переменная spec, содержащая массив рассчитанных значений спектра.

Opt.Output = `separate`; – рисует отдельный спектр для каждого из переходов, переменная spec представляет собой матрицу, каждый ряд которой соответствует одному переходу;

Opt.Output = `summed`; – весь спектр записывается как один ряд, это значение принимается по умолчанию.

Поле nKnots - задает число ориентаций, для которых вычисляется спектр порошка.

Задается в виде [N1] или [N1 N2]. N1 – определяет число ориентаций между 0° и 90°, типичные значения для N1 лежат между 10 и 91, чем больше анизотропия и уже линии в спектре, тем большее количество ориентаций нужно задавать для получения гладкого спектра; N2 – фактор интерполяции; если N2 = 4, то между каждой парой точек добавляется по три точки, значения в которых находятся интерполированием, например, Opt.nKnots = 91; – промежуток в 1°, без интерполяции;

Opt.nKnots = [31 6]; – промежуток в 3°, 6-кратная интерполяция (добавляется по 5 точек) Поле Symmetry – определяет симметрию, используемую при симуляции спектра порошка. Может принимать значения `auto` – симметрия находится автоматически для данной спиновой системы (стоит по умолчанию);

`Dinfh` – соответствует линии от = 0° до 90° ( = 0);

`D2h` – один октант;

`C2h` – два октанта;

`Ci` – соответствует полусфере (4 октанта) Поле Trasitions – в этом поле задаются пары уровней энергии, переходы между которыми должны быть представлены в вычисляемом спектре.

Представляет собой матрицу m x 2, где m – число переходов, например, присутствовать линии, соответствующие двум переходам, между первым и вторым и третьим и четвертым уровнями (номером 1 всегда обозначается уровень с наименьшей энергией).

Поле Threshold – определяет граничное значение для амплитуды линии (по умолчанию стоит значение 10-4). Любой переход, чья относительная средняя амплитуда меньше этого значения, в вычислениях не участвует. Относительная средняя амплитуда наиболее интенсивного перехода равна 1.

Поле Intensity – управляет механизмом расчета интенсивности линий.

Со значением `on` (по умолчанию) интенсивности линий вычисляются в соответствии с правилами отбора, разрешенные переходы будут более интенсивны, нежели квази-запрещенные; значение `off` полагает все вероятности переходов равными 1, т. е. запрещенные и разрешенные переходы будут иметь одинаковую интенсивность.

Поле Method – определяет метод, которым функция pepper находит резонансные поля.

Значение `matrix` соответствует диагонализации матрицы, учитываются все взаимодействия, включая квадрупольное, подходит для любых систем;

`perturb` – соответствует теории возмущений второго порядка, подходит только для систем с одним электронным спином, не включает в себя ядерное зеемановское взаимодействие, вычисляются только разрешенные переходы;

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

используются диагонализация матрицы и теория возмущений, по умолчанию этот режим отключен (0).

Если задано значение Opt.Perturb = 1; – то все ядра учитываются в первом порядке теории возмущений, для них не вычисляются квадрупольное и ядерное зеемановское взаимодействия, учитывается только сверхтонкое взаимодействие; есть возможность вручную определять те ядра, для которых расчет следует производить методом теории возмущений, например, Sys.Nucs = '63Cu,14N,1H,1H';

рассмотрения методом теории возмущений.

5. Пример использования модуля EasySpin: ион Yb3+ в В качестве примера рассмотрим расчет спектра ЭПР тригонального центра иона Yb3+ в кристаллическом поле кубической симметрии (BaF2) с помощью модуля EasySpin. Известно, что при легировании матрицы BaF2 ионами иттербия образуется два типа центров: кубический с изотропным g-фактором (g Основное состояния иона Yb3+ 2F7/2. J = 7/2, и в кристаллическом поле уровни энергии распадаются на 4 крамерсовых дублета. Так как расстояние от основного до первого возбужденного дублета велико по сравнению с энергией тепловых колебаний при температуре жидкого гелия 4.2 К, то при рассмотрении можно ограничиться нижним дублетом и использовать в расчетах эффективный спин. Кроме зеемановского, будем учитывать также сверхтонкое взаимодействие. В природе существует всего 7 изотопов иттербия.

распространенность 14,4%) и Yb (I = 5/2, 16.2%). Известны величины констант сверхтонкого взаимодействия В рамках приближения, учитывающего только матричные элементы между сверхтонким и зеемановским взаимодействиями. Главные оси тензоров g и A совпадают, и выполняется соотношение Для иона Yb3+ gJ = 8/7. Таким образом, зная главные значения тензора g, мы легко находим главные значения тензора константы сверхтонкой структуры.

Спиновый гамильтониан для тригонального центра будет иметь вид (для одного типа изотопа):

Для симуляции спектра ЭПР с использованием модуля EasySpin необходимо параметры спинового гамильтониана представить в понятной программе форме.

Итак, спин равен, g-фактор обладает аксиальной симметрией, существует три типа изотопов иттербия с различными ядерными спинами и значениями констант сверхтонкой структуры, также имеющих аксиальную симметрию. Определим численные значения величин, которые затем в качестве переменных будем заносить в поля структуры Sys. Компоненты тензора сверхтонкой структуры задаем в МГц.

gz = 2.763; g1 = 3.768; gJ = 8/7;

A171 = 888.62; A173 = -241.85;

Далее рассчитаем главные значения тензора константы сверхтонкой структуры в соответствие с соотношением (1).

Atrz1 = gz*A171/gJ;

Atr11 = g1*A171/gJ;

Atrz2 = gz*A173/gJ;

Atr12 = g1*A173/gJ;

Ядра будем задавать с помощью функции nucspinadd. Структуру Sys описывать в виде строки. Будем считать, что линии в спектре имеют лоренцеву форму с шириной peak-to-peak в 3 мТ.

Sys = struct('S',1/2,'g',[g1 g1 gz],'lwpp',[0 3]);

Вычисления будем производить для каждого изотопа по отдельности, поэтому далее записываем Sys1 = nucspinadd(Sys,'171Yb',[Atr11 Atr11 Atrz1]);

Sys2 = nucspinadd(Sys,'173Yb',[Atr12 Atr12 Atrz2]);

Таким образом, с описанием спиновой системы закончено, и мы перейдем к заданию экспериментальных параметров. Предполагаем, что используемый нами ЭПР спектрометр работает в X-диапазоне (9.4 ГГц), измерения проводятся при температуре 5 К. Желаемый диапазон магнитных полей от 100 до 300 мТ.

Рассчитанный спектр будет представлен в виде первой производной спектра поглощения. Итак, Exp.mwFreq = 9.4;

Exp.Range = [100 300];

Exp.Temperature = 5;

Exp.Harmonic = 1;

Здесь мы используем задание структуры построчно, что иногда является более удобным. Осталось определить необходимые ориентации.

Оси тригонального центра направлены вдоль пространственных диагоналей куба элементарной ячейки кристалла. Следовательно, существует типа тригональных центров. Ось z молекулярной системы координат направим вдоль одной из диагоналей, тогда другие 3 оси третьего порядка будут составлять с осью z одинаковые углы по 60°. Расчеты проведем для одного направления магнитного поля, когда внешнее постоянное магнитное поле B направлено вдоль одной из осей C3, то есть для случая B0 z. В этой ситуации три тригональных центра с различными направлениями, оси которых составляют с осью z углы по 60°, являются, вообще говоря, эквивалентными. То есть достаточно вычислить значения резонансного поля для двух случаев: 1) ось магнитного поля совпадает с осью центра; 2) направление магнитного поля составляет с осью центра угол в 60°. Таким образом, вклад в спектр тригонального центра иона Yb3+ в кристаллическом поле кубической симметрии дают 6 групп линий, по две возможные ориентации для каждого из трех изотопов. Результирующий спектр будет просто их суммой (следует отметить, что в этом примере мы не учитывали реальную амплитуду линий).

Поэтому запишем Exp.Orientations = [0 0 0];

[Field,Spec0] = pepper(Sys,Exp);

[Field,Spec1] = pepper(Sys1,Exp);

[Field,Spec2] = pepper(Sys2,Exp);

Exp.Orientations = [0 acos(1/3) 0];

[Field,Spec03] = pepper(Sys,Exp);

[Field,Spec13] = pepper(Sys1,Exp);

[Field,Spec23] = pepper(Sys2,Exp);

Spec = Spec0+Spec1+Spec2+Spec03+Spec13+Spec23;

Наконец, полный текст программы clc;

clear all;

gz = 2.763; g1 = 3.768; gJ = 8/7;

A171 = 887.2; A173 = -243.3;

Atrz1 = gz*A171/gJ;

Atr11 = g1*A171/gJ;

Atrz2 = gz*A173/gJ;

Atr12 = g1*A173/gJ;

Sys = struct('S',1/2,'g',[g1 g1 gz],'lwpp',[0 3]);

Sys1 = nucspinadd(Sys,'171Yb',[Atr11 Atr11 Atrz1]);

Sys2 = nucspinadd(Sys,'173Yb',[Atr12 Atr12 Atrz2]);

Exp.mwFreq = 9.4;

Exp.Range = [100 300];

Exp.Temperature = 5;

Exp.Harmonic = 1;

Exp.Orientations = [0 0 0];

[Field,Spec0] = pepper(Sys,Exp);

[Field,Spec1] = pepper(Sys1,Exp);

[Field,Spec2] = pepper(Sys2,Exp);

Exp.Orientations = [0 acos(1/3) 0];

[Field,Spec03] = pepper(Sys,Exp);

[Field,Spec13] = pepper(Sys1,Exp);

[Field,Spec23] = pepper(Sys2,Exp);

Spec = Spec0+Spec1+Spec2+Spec03+Spec13+Spec23;

plot(Field, Spec);

Результатом работы такой программы будет представление спектра в графическом виде, как показано на рисунке 2.

Рис.2. Спектр ЭПР тригональных центров ионов Yb3+ в кристалле BaF2, рассчитанный с помощью программного модуля EasySpin.

1. Stoll S., Schweiger A. EasySpin, a comprehensive software package for spectral simulation and analysis in EPR // Journal of Magnetic Resonance. 2006. V.78. P.

42-55.

2. Stoll S., Schweiger A. EasySpin: Simulating cw ESR spectra // Biol. Magn.

Reson. 2007. V. 27. P. 299-321.



Похожие работы:

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ПЕЧАТИ Е. В. Никульчев ИДЕНТИФИКАЦИЯ ДИНАМИЧЕСКИХ СИСТЕМ НА ОСНОВЕ СИММЕТРИЙ РЕКОНСТРУИРОВАННЫХ АТТРАКТОРОВ Учебное пособие Москва 2010 1 УДК 519.711.3 ББК 22.18 Н 65 Рецензенты: А. И. Дивеев, д. т. н., профессор (Вычислительный центр РАН); Е. Е. Ковшов, д. т. н., профессор (МГТУ СТАНКИН). Н 65 Никульчев Е. В. Идентификация динамических систем на основе симметрий реконструированных аттракторов : учеб....»

«Министерство образования и науки Российской Федерации ТОМСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ УТВЕРЖДАЮ Декан геолого-географического факультета Г.М. Татьянин марта 2008 г. ВЫПУСКНАЯ РАБОТА БАКАЛАВРА ГЕОЛОГИИ Направление : 020300 - Геология Учебно-методическое пособие Томск 2008 Выпускная работа бакалавра геологии: учебно-методич. пособие / составители: А.И. Чернышов, Н.И. Савина: Том. гос. ун.-т. – 3-е изд., доп. и перераб. – Томск, 2008. – 33 с. Учебно-методическое пособие составлено на основе...»

«Министерство образования Российской Федерации ТОМСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ СИСТЕМ УПРАВЛЕНИЯ И РАДИОЭЛЕКТРОНИКИ (ТУСУР) Кафедра автоматизированных систем управления (АСУ) Е.Н. Сафьянова ДИСКРЕТНАЯ МАТЕМАТИКА Часть 1 Учебное пособие 2000 Сафьянова Е.Н. Дискретная математика. Часть 1: Учебное пособие. Томск: Томский межвузовский центр дистанционного образования, 2000. 106 с. Учебное пособие рассмотрено и рекомендовано к изданию методическим семинаром кафедры автоматизированных систем...»

«МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ИМЕНИ М.В.ЛОМОНОСОВА ФАКУЛЬТЕТ ВЫЧИСЛИТЕЛЬНОЙ МАТЕМАТИКИ И КИБЕРНЕТИКИ Анализ данных, обучение по прецедентам, логические игры, системы WEKA, RapidMiner и MatLab (ПРАКТИКУМ НА ЭВМ КАФЕДРЫ МАТЕМАТИЧЕСКИХ МЕТОДОВ ПРОГНОЗИРОВАНИЯ) УЧЕБНОЕ ПОСОБИЕ Дьяконов А.Г. Москва, 2010 г. Печатается по решению редакционно-издательского совета факультета вычислительной математики и кибернетики МГУ имени М.В. Ломоносова Рецензенты: Ю.И. Журавлёв, д.ф.-м.н., профессор,...»

«УДК 364.4(075.8) ББК 65.272я73 МИНОБРНАУКИ РОССИИ У 91 ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ ПОВОЛЖСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ СЕРВИСА (ФГБОУ ВПО ПВГУС) Рецензент Кафедра Социальные технологии к.ф.н., доц. Рузова Л. А. УЧЕБНО-МЕТОДИЧЕСКИЙ КОМПЛЕКС по дисциплине Социальная реабилитация для студентов направления 040100.62 Социальная работа Учебно-методический комплекс по дисциплине Социальная У 91 реабилитация / сост. Л. И....»

«Artrit_A5 copy-new.qxd 30/08/2007 16:14 Page 1 ВРАЧ – ПАЦИЕНТУ ИНФОРМАЦИЯ О РЕВМАТИЧЕСКИХ ЗАБОЛЕВАНИЯХ РЕВМАТОИДНЫЙ АРТРИТ (методическое пособие по материалам научно-практической конференции (школы) Ревматоидный артрит. Современные методы лечения 15 марта 2007 г.) Межрегиональная общественная организация инвалидов Ревматологическая Ассоциация Надежда www.revmo-nadegda.ru e-mail: revmo-nadegda@mail.ru Artrit_A5 copy-new.qxd 30/08/2007 16:14 Page Artrit_A5 copy-new.qxd 30/08/2007 16:14 Page...»

«Федеральное агентство по образованию Государственное образовательное учреждение высшего профессионального образования Казанский государственный технологический университет Институт технологии легкой промышленности, моды и дизайна Кафедра Дизайн ПРОГРАММА ПРОИЗВОДСТВЕННОЙ ПРАКТИКИ Методические указания для студентов IV курса Казань КГТУ 2007 УДК 665.6:033.28 Составители: проф. В.В. Хамматова доцент Е.В. Кумпан зав. худ. маст. А.И. Вильданова Программа производственной практики: метод. указания...»

«УДК 51(075.8) ББК 22.1я73 МИНОБРНАУКИ РОССИИ У 91 ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ ПОВОЛЖСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ СЕРВИСА (ФГБОУ ВПО ПВГУС) Кафедра Высшая математика Рецензент к.п.н., доц. Киричек Г. А. УЧЕБНО-МЕТОДИЧЕСКОЕ ПОСОБИЕ по дисциплине Математика. Элементы высшей математики для студентов всех специальностей СПО. Раздел Элементы теории рядов и дифференциальных уравнений Учебно-методическое пособие по...»

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ ИНСТИТУТ УПРАВЛЕНИЯ, ИНФОРМАЦИИ И БИЗНЕСА КАФЕДРА АИС ДИСКРЕТНАЯ МАТЕМАТИКА Методические указания к выполнению контрольной работы студентами специальности АИС Часть I Ухта 2003 УДК 60.5 Л41 Лихачев В.Н., Баскакова Ю.Л. Дискретная математика. Методические указания к выполнению контрольной работы студентами специальности АИС. – Ухта: Институт управления, информации и бизнеса, 2003. – 15 с. Методические указания предназначены для студентов...»

«1 Московский государственный университет геодезии и картографии МИИГАиК Кафедра высшей геодезии Шануров Геннадий Анатольевич Матрицы в геодезии. Применение матриц в обработке и оценке точности результатов геодезических измерений и определений Учебное пособие по курсам Высшая геодезия и Геотроника для студентов и аспирантов геодезических специальностей Москва 2013 год 2 Содержание Введение 1. Измеряемые величины и определяемые величины 1.1. Линейные и угловые измеряемые величины 1.2. Связь...»

«1 СОДЕРЖАНИЕ КУРСА УЧЕБНО-ЛЕТНОЙ ПОДГОТОВКИ СПОРТИВНЫХ АВИАЦИОННЫХ ОРГАНИЗАЦИЙ ДОСААФ СССР НА САМОЛЕТАХ (КУЛПа-САО-С-86)*. Курс учебно-летной подготовки спортивных авиационных организаций ДОСААФ СССР на самолетах является основным руководящим документом, определяющим содержание, объем, порядок и последовательность обучения постоянного** и переменного летного состава. Курс состоит из двух частей: часть первая - теоретическая подготовка; часть вторая - летная подготовка. Часть первая содержит...»

«Р. Н. Абалуев, Н. Е. Астафьева, Н. И. Баскакова, Е. Ю. Бойко, О. В. Вязовова, Н. А. Кулешова, Л. Н. Уметский, Г. А. Шешерина ИНТЕРНЕТ-ТЕХНОЛОГИИ В ОБРАЗОВАНИИ ИЗДАТЕЛЬСТВО ТГТУ Министерство образования Российской Федерации Тамбовский региональный центр Федерации Интернет-Образования Тамбовский областной институт повышения квалификации работников образования Р. Н. Абалуев, Н. Е. Астафьева, Н. И. Баскакова, Е. Ю. Бойко, О. В. Вязовова, Н. А. Кулешова, Л. Н. Уметский, Г. А. Шешерина...»

«Министерство образования Российской Федерации Нижегородский государственный университет им. Н.И. Лобачевского ЗАДАЧИ И МЕТОДЫ КОНЕЧНОМЕРНОЙ ОПТИМИЗАЦИИ Учебное пособие Часть 3. Д.И. Коган Динамическое программирование и дискретная многокритериальная оптимизация Издательство Нижегородского университета Нижний Новгород 2004 УДК 519.6 Коган Д.И. Динамическое программирование и дискретная многокритериальная оптимизация: учебное пособие. Нижний Новгород: Изд-во Нижегородского ун-та, 2004. 150 с....»

«МАТЕМАТИКА К УЧЕБНИКУ М.И. Моро и др. (М.: Просвещение) 2 -е издан ие, пер ер а б о тан н о е 4 класс МОСКВА • ВАКО УДК 372.851 ББК 74.262.21 К64 Контрольно измерительные материалы. Математика: К64 4 класс / Сост. Т.Н. Ситникова. – 2-е изд., перераб. – М.: ВАКО, 2011. – 96 с. – (Контрольно-измерительные материалы). ISBN 978-5-408-00454-6 В пособии представлены контрольно-измерительные материалы по математике для 4 класса. Все задания соответствуют программе общеобразовательных учреждений и...»

«Уважаемые выпускники! В перечисленных ниже изданиях содержатся методические рекомендации, которые помогут должным образом подготовить, оформить и успешно защитить выпускную квалификационную работу. Рыжков, И. Б. Основы научных исследований и изобретательства [Электронный ресурс] : [учебное пособие для студентов вузов, обучающихся по направлению подготовки (специальностям) 280400 — Природообустройство, 280300 — Водные ресурсы и водопользование] / И. Б. Рыжков.— СанктПетербург [и др.] : Лань,...»

«Федеральная таможенная служба Государственное казенное образовательное учреждение высшего профессионального образования Российская таможенная Академия RUSSIAN CUSTOMS ACADEMY ПРИКA3 27 12 2012г №1385 г. Люберцы О методических указаниях по разработке учебных планов (рабочих учебных планов) направлений подготовки (специальностей). В целях обеспечения соответствия содержания подготовки обучающихся требованиям Федеральных государственных образовательных стандартов высшего профессионального...»

«Kim Fleischer Michaelsen Кормление и питание грудных детей и детей раннего возраста Lawrence Weaver Francesco Branca Aileen Robertson Кормление и питание грудных детей и детей unicef раннего возраста Методические рекомендации для Европейского региона ВОЗ с особым акцентом на республики бывшего Советского Союза Региональные публикации ВОЗ, Европейская серия, 87 № Всемирная организация здравоохранения была создана в 1948 г. в качестве специализированного учреждения Организации Объединенных Наций,...»

«И.П. Заикин, А.В. Тоцкий, С.К. Абрамов, В.В. Лукин ПРОЕКТИРОВАНИЕ АНТЕННЫХ УСТРОЙСТВ СВЧ 2005 МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ УКРАИНЫ Национальный аэрокосмический университет им. Н.Е. Жуковского Харьковский авиационный институт И.П. Заикин, А.В. Тоцкий, С.К. Абрамов, В.В. Лукин ПРОЕКТИРОВАНИЕ АНТЕННЫХ УСТРОЙСТВ СВЧ Учебное пособие Харьков ХАИ 2005 УДК 621.396.67 Проектирование антенных устройств СВЧ / И.П. Заикин, А.В. Тоцкий, С.К. Абрамов, В.В. Лукин. – Учеб. пособие. – Харьков: Нац....»

«УДК 656.7(075.8) ББК 39.511я73 А44 ISBN 978-601-7086-59-6 АВИАЦИОННАЯ КЛИМАТОЛОГИЯ (глобальная, региональная и маршрутная) АКЫЛБАЕВА К.И. АЛМАТЫ 2011 АКЫЛБАЕВА К.И. АЛМАТЫ 2011 Акылбаева К.И. Авиационная Климатология (глобальная, региональная и маршрутная). Алматы, 2011, 399стр., 156рис., 35 табл., 1график. Работа сделана по результатам анализа зарубежных учебников. В работе кратко представлена глобальная, региональная и маршрутная Авиационная Климатология Земного шара, дана краткая...»

«ЦКП Материаловедение и диагностика в передовых технологиях при ФТИ им. А.Ф. Иоффе МЕТОДИКА ИЗМЕРЕНИЯ ПАРАМЕТРОВ ГЛУБОКИХ УРОВНЕЙ В ПОЛУПРОВОДНИКОВЫХ СТРУКТУРАХ Определение концентрации, энергии ионизации и сечения захвата дефектов с глубокими уровнями методом нестационарной спектроскопии глубоких уровней. Методические указания к лабораторным работам по диагностике материалов Санкт-Петербург 2010 Оглавление: Метод измерения 3 Установка для проведения измерений 10 Подготовка образцов к проведению...»










 
© 2013 www.diss.seluk.ru - «Бесплатная электронная библиотека - Авторефераты, Диссертации, Монографии, Методички, учебные программы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.