Метод вторичного квантования с неортогональным базисом и его приложение к теории локальных магнитных полей на ядрах диамагнитных ионов в кристаллах с незаполненными 3d- и 4f-оболочками тема автореферата и диссертации по физике, 01.04.02 ВАК РФ

Аникеенок, Олег Алексеевич АВТОР
доктора физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Казань МЕСТО ЗАЩИТЫ
2015 ГОД ЗАЩИТЫ
   
01.04.02 КОД ВАК РФ
Автореферат по физике на тему «Метод вторичного квантования с неортогональным базисом и его приложение к теории локальных магнитных полей на ядрах диамагнитных ионов в кристаллах с незаполненными 3d- и 4f-оболочками»
 
Автореферат диссертации на тему "Метод вторичного квантования с неортогональным базисом и его приложение к теории локальных магнитных полей на ядрах диамагнитных ионов в кристаллах с незаполненными 3d- и 4f-оболочками"

На правах рукописи №

Аникеенок Олег Алексеевич

МЕТОД ВТОРИЧНОГО КВАНТОВАНИЯ С НЕОРТОГОНАЛЬНЫМ БАЗИСОМ И ЕГО ПРИЛОЖЕНИЕ К ТЕОРИЙ ЛОКАЛЬНЫХ МАГНИТНЫХ ПОЛЕЙ НА ЯДРАХ ДИАМАГНИТНЫХ ИОНОВ В КРИСТАЛЛАХ С НЕЗАПОЛНЕННЫМИ 3(1- и 4Г-ОБОЛОЧКАМИ

01.04.02. - теоретическая физика

1 9 ФЕВ 2015

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

КАЗАНЬ-2015

005559322

005559322

Работа выполнена на кафедре квантовой электроники и радиоспектроскопии ФГАОУ ВПО «Казанский (Приволжский) федеральный университет»

Официальные оппопепты: Джепаров Фридрих Саламоиович

доктор физико-математических наук, профессор, ГНЦ РФ ИТЭФ, г. Москва, заведующий лабораторией нейтронной физики

Вальков Валерий Владимирович

доктор физико-математических наук, профессор, ИФ СО РАН, г. Красноярск, заведующий лабораторией теоретической физики

Пухов Константин Константинович

доктор физико-математических наук, Институт общей физики им. A.M. Прохорова, г. Москва, ведущий научный сотрудник

Ведущая организация: ФГАОУ ВПО «Уральский федеральный

университет», г.Екатеринбург

Научный консультант: Еремин Михаил Васильевич

доктор физико-математических наук, профессор, Казанский (Приволжский) федеральный университет

Защита состоится 21 мая 2015 г. в 14 час. 30 мин. на заседают диссертационного совета Д 212.081.15 при ФГАОУ ВПО «Казанский (Приволжский) федеральный университет» по адресу: 420008, г. Казань, ул.Кремлевская, д. 16а, ауд. 110.

С диссертацией можно ознакомиться в библиотеке имени Н.И. Лобачевского Казанского (Приволжского) федерального университета (г.Казань, ул.Кремлевская, д. 35). Электронная версия размещена на официальных сайтах ВАК при Министерстве образования и науки РФ (vak2.ed.gov.ru) и Казанского (Приволжского) федерального университета kpfu.ru.

Автореферат разослан 14 февраля 2015 г.

Ученый секретарь у

диссертационного совета 7 sfa^ М.В. Еремин

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. При нахождении спектра многочастичных систем практически всегда используют одночастичное приближение, т.е. предполагается, что состояние каждого электрона может быть описано некоторой волновой функцией. В этом случае спектр и состояние системы могут быть найдены как экстремумы функционала среднего значения гамильтониана такой системы. При этом волновую функцию системы естественно представить как некоторую суперпозицию слетеровских детерминантов, составленных из этих искомых одночастичных функций. Как правило, на эти одиочастнчные функции накладывается условие ортонормированности. В то же время очевидно, что функционалы с условием ортонормированности одночастичных функций и с отсутствием условий ортогональности одночастичных функций являются разными функционалами и в общем случае должны иметь рапные экстремумы. В работе В.А. Фока [1] отмечалось, что при учете возбужденных конфигураций следует потребовать только частичную ортогональность орбиталей. В ряде работ, например, [2,3] данный вопрос изучался численным решением уравнений Хартри - Фока для простейших конфигураций при разных условиях относительно ортогональности радиальных орбиталей. Был сделан вывод, что чем больше условий ортогональности, тем результата хуже. Дальнейшее развитие теории и практическое применение метода самосогласованного поля Хартри - Фока с использованием неортогональных орбиталей задерживается из-за математических трудностей. Это относится как к вычислению матричных элементов операторов на слетеровских детерминантах, построенных с использованием частично неортогональных орбиталей, так и к вычислениям, при решении самих уравнений Хар1:ри - срока на подобных детерминантах в координатном представлении.

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

вторичного квантования также как и в случае ортонормированного одночастичного базиса. В работах [4,5] была предпринята попытка развития метода вторичного квантования с частично неортогональным одночастичным базисом для фермионов. В этих работах «бра»-вектора дуального базиса не являются эрмитово-сопряженными к «кет»-векторам и как следствие одночастичные и двухчастичные операторы не являются эрмитовыми. Антикоммутатор операторов рождения и уничтожения орбиталей в отличие от обычного символа Кронекера, получался равным соответствующему интегралу перекрывания. В силу этого возникают трудности с нормировкой многочастичных состояний и т.д. Каких-либо приложений предлагаемого метода вторичного квантования к реальным системам не последовало и вопрос о так называемой «катастрофе неортогональности» остался открытым. Таким образом, проблема нахождения точного вида эрмитовых операторов в представлении вторичного квантования с частично неортогональным одночастичным базисом осталась не решенной.

Все вышеизложенное указывает на актуальность темы диссертации.

Цель работы. Целью настоящей работы является построение метода вторичного квантования с линейно независимым, неортогональным одночастичным базисом. На настоящем этапе рассматривается приложение этого метода к системам, в которых возможно применение теории возмущений. К таким системам могут быть отнесены ионные кристаллы. В качестве основных объектов исследования выбраны примесные центры с незаполненными Зс1- и 4^оболочками, а именно, лигандные сверхтонкие взаимодействия (ЛСТВ) в них.

Проблемы в интерпретации экспериментальных данных по ЛСТВ возникают уже для ионов группы железа. Применение стандартного варианта метода молекулярных орбиталей в ряде случаев затруднено. Среди примесных центров с незаполненной 3<1-оболочкой к ним относятся: парамагнитные центры с отсутствующими а-связями, комплексы с орбитальным вырождением и большой класс низкосимметичных систем. В

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

Более того, выполненные для ионов С<13+ расчеты параметров спиновой плотности, перенесенной на ядра лигандов, с использованием замены волновых функций 4Р-злектроноп на соответствующие молекулярные орбитали, даже по знаку не согласуются с экспериментальными данными. В этой связи возникает задача об учете пространственного распределения еще и внешних электронных оболочек. Однако интегралы перекрывания 6й- , 6р-и 5с]-оболочек с лигандами достаточно большие. При этом возникает проблема сходимости рядов с интегралами неортогональности. Данная проблема впервые возникла в теории химической связи и получила название «катастрофы неортогональности». Она известна ещё по работам Слетера, но так и осталась неразрешенной до настоящего времени.

Положения, выносимые на защиту.

1. Получены выражения для одночастичного и двухчастичного операторов в методе вторичного квантования с линейно независимым, частично неортогональным одночастичным базисом. Операторы имеют вид ряда по п — ым степеням коммутатора, введенного нами оператора <3 с оператором Н, с коэффициентами сп = Е2П/[2"п(2п)!], где Е2п являются числами Эйлера. Матричные элементы оператора Н являются линейной комбинацией матричных элементов гамильтониана Н системы. Коэффициенты этой линейной комбинации выражаются через матричные элементы матрицы (1+8)"', где I - единичная матрица, а Б - матрица интегралов перекрывания выбранного одночастичного базиса. Матричные элементы оператора (} также выражаются через матричные элементы матрицы (1+8)"'. Для сходимости полученного ряда требуется только ограниченность матричных элементов множества п — ых степеней коммутаторов.

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

3. Дано обобщение выражений для амплитуд перехода электрона металл-лиганд (аналога параметра ков&гентности в методе молекулярных орбиталей) позволяющих рассчитать их значения, не предполагая малости соответствующих интегралов перекрывания. Предложенные выражения для амплитуд перехода объясняют близость экспериментальных значений по лигандной сверхтонкой структуре для примесных центров УЬ3+:СзСаР3 и УЬ3+:С82№УРб, несмогря на то, что ионы УЪ3+ внедряются в разные кристаллы и в одном случае замещение изовалентно, а в другом неизовалентно.

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

5. Продемонстрировано, что развитая теория и предложенные механизмы перенесенных магнитных полей на ядра диамагнитных ионов позволяет объяснить основные особенности формирования локальных полей на ядрах фтора в ряде фторидов: ¥Ъ3+:СвСаРз и УЬ3+:С82МаУРб.

6. Рассчитаны значения локальных магнитных полей в ЬаМпОз, обогащенном изотопом пО, при различных вариантах кооперативного упорядочения орбиталей ионов Мп3+. Наилучшее согласие с экспериментами достигается в предположении, что при кооперативном упорядочении электроны е„ -подоболочки находятся в состояниях с волновой функцией вида с^Ъг1-г^ + с^х1-у1) (определена в локальных осях октаэдрических

фрагментов Мп06). Так, при Т=298К коэффициенты равны: С1 = 0.995, с2 = - 0.10. Отмечена важность эффектов квантовой интерференции пропорциональных произведению С1С2 Таким образом, опираясь на экспериментальные данные ЯМР в совместной с экспериментаторами работе, однозначно определены волновые функции основных состояний ионов Мп, и тем самым установлена структура кооперативного орбитального упорядочения в этом соед1шении в парамагнитной фазе.

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

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

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

Апробация работы. Основные результаты работы докладывались и обсуждались на конференциях «XIX International Seminar on Modern Magnetic Resonances, RAMIS» (Poznan-Bedlewo, Poland, 2001), «XI-th Feofilov symposium on spectroscopy of crystals activated by rare earth and transition metal ions» (Kazan, Russia, 2001), «European Conference Physics of Magnetism» (Poznan, Poland, 2002), «XII Всероссийская конференция. Оптика и спектроскопия конденсированных сред» (Краснодар, Россия, 2006), «International conference. Modern development of magnetic resonance» (Kazan, Russia, 2007), «International conference. Modern development of magnetic resonanc» (Kazan, Russia, 2007), «XV всероссийская конференция. Оптика и спектроскопия конденсированных сред» (Краснодар, Россия, 2009), «XIV International Feofilov Symposium. Spectroscopy of crystals doped with rare earth and transition metal ions» (St.-Petersburg, Russia, 2010), «VII Международная конференция. Фундаментальные проблемы оптики» (Санкт-Петербург, Россия, 2012), «Actual problems of magnetic resonance and its application. XV International Youth Scientific School» (Kazan, Russia, 2012). «XV International Feofilov Symposium. Spectroscopy of crystals doped with rare earth and transition metal ions» (Kazan, Russia, 2013).

Публикации. Основные результаты диссертации опубликованы в 17 статьях в российских и международных журналах, включенных в базу SCOPUS (15 журналов входят в перечень ВАК), а также в материалах и тезисах вышеперечисленных конференций.

Личный вклад. Все теоретические результаты, представленные в диссертации, получены лично автором.

Структура и объем диссертации. Диссертация состоит из введения, 6 глав и заключения, изложена на 191 страницах машинописного текста, содержит 3 рисунка, 2 диаграммы, 18 таблиц и 3 приложения. Список используемой литература содержит 100 наименований.

Содержание диссертации.

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

В первой главе обсуждаются результаты [А1-А7], [6-11], которые получены методом «приближенного вторичного квантования», в котором операторы в представлении вторичного квантования получены с точностью до квадратов интегралов перекрывания. Подчеркивается, что, даже рассматривая параметры ковалентности как подгоночные, без привлечения новых механизмов формирования локальных магнитных полей на ядрах диамагнитных ионов, невозможно объяснять имеющиеся экспериментальные данные по JICTB примесных центров. В особенности это четко видно в случае редкоземельных соединений.

Однако при выполнении этой программы возникает ряд принципиальных трудностей, которые фактически формулируют постановку задачи. В качестве иллюстрации приведем некоторые качественные оценки. Запишем гамильтониан лигандного сверхтонкого взаимодействия Нлсгв, т.е. оператор взаимодействия эффективного сшша примесного иона с ядерным спином выделенного лиганда, через изотропную и анизотропную константы As и Ар соответственно.

Нлств = 4 (SI) + Ар (2SZIZ-SXIX - SyIy), где

Ар = А'р + А,,

4,=

ШпЕп

Р УЬ3+

— спиновые плотности, п1, х - квантовые числа орбиталей ионов.

Знак с^ может быть как больше, так и меньше нуля в зависимости от механизма возникновения спиновой плотности.

Рис 1. Модель примесного центра СБСаРз: УЪ3+.

Приведенный оператор взаимодействия относится к одной из выделенных пар УЪ3+ — Р" и записан в локальной системе координат с осью г вдоль оси данной пары. Наблюдаемая величина Ар заметно отличается от

магнитного диполь-дипольного взаимодействия Аа примесного иона с ядром лиганда, что явно указывает на проявление эффектов ковалентности

А'р = + сМ + сШ,УР , ар =

На диаграмме 1.1 приведен процесс соответствующий образованию ковалентной связи 4Г-электрона с 2р (2э) — электронами фтора 0).

4\

2э,2р

Диаграмма 1.1. Формирование перенесенного сверхтонкого поля на ядре фтора, в результате образования ковалентной связи.

В случае примесного центра КМ^3: Тт2+, УЪ3~, Б = 1/2 она является определяющей [8]. Цифрами обозначена последовательность операторов в поправке второго порядка теории возмущений.

Видно, что магнитное поле, создаваемое оставшимся электроном на ядре лиганда, имеет тот же знак, что и поле, создаваемое электроном центрального иона, т.е. изотропная константа А1 > 0, что и наблюдается на эксперименте.

На диаграмме 1.2 приведен процесс с участием вышележащих незаполненных оболочек (с'''< 0).

Диаграмма 1.2. Вклад в перенесенное сверхтонкое поле от процессов переноса в пустые оболочки 5с1- или бе-оболочки для конфигурации .

Из диаграммы 1.2 видно, что магнитное поле, создаваемое оставшимся электроном, имеет на ядре лпганда знак противоположный полю, создаваемому ковалентным вкладом. Энергия перехода в 5(1-, бе-оболочки больше энергии перехода в валентную оболочку. Однако следует учесть, что перекрывание 5(1-, бе-оболочек с орбиталями лиганда заметно больше перекрывания валентной 41-оболочьи. Кроме того, в отличие от ковалептного вклада, где согласно правилам отбора, возможен только переход 2з - 4ГО, в процессах с пустыми оболочками участвуют все семь орбиталей конфигурации ЛС. Таким образом, процессы диаграммы 1.2 могут объяснить

отрицательное значение Д< О, которое наблюдается на эксперименте в случае примесных центров СаР2: Еи2+, Ос13+, Б = 7/2 [12]. Легко видеть, что процесс, начинающийся переходом электрона с направлением спина противоположным направлению, приведенному на диаграмме 1.2, запрещен правилами отбора.

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

• интегралы перекрывания (бр\в),(бз\в),(5<1\в) и т.д., где |б>) -

орбитали лигандов не малы и ряды неортогональности становятся расходящимися;

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

Исходя из этого, цель работы можно сформулировать следующим образом:

а) развить метод, позволяющий проводить расчеты при произвольных интегралах неортогональности;

в) избежать катастрофы неортогоналъности, т.е. определить критерий возможности использования того или иного одночастичного базиса;

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

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

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

V»; + = ^. я л + я А' = а>1 + аХ■ = 0. а 1°> = 0 (5)

и вводятся функции

({£}Н0Кай~абма&> 1М>=«-,••« 1°)- (6)

Показано, что если _ слетеровские детерминанты с

орбиталями неортогонального базиса и квантовыми числами {<?}, {7}, то матричный элемент

К>Нфм) = «,-«I0)' (7)

где А — одночастичный оператор, а оператор /г имеет вид:

/1 = ^

ехр

и аналогично g — двухчастичный оператор, а оператор

g = N

ехр

Уz

(8)

(9)

где N — знак нормального произведения.

В разделе 2.3, используя результаты предыдущего пункта, рассматривается интеграл перекрывания между слетеровскими

детерминантами и |Ф{7))» кот°рый может быть вычислен по формуле:

ехр

П< 1°)- (10)

Показано, что оператор, входящий в (10), может быть представлен в

виде:

¡ф = N

ехр ^а>Лв\&)

= ехр(2),

(И)

где оператор б - определяется следующим выражением:

в=I«;«, И") - (?>(/, (12)

п=..

= - матричные элементы матрицы перекрывания 5,

составленной в базисе одноэлектронных орбиталей. Дано определение матрицы 1п(/ + 5) не основанное на обычно подразумеваемом операторном ряде типа ряда Тейлора [А11]. Мы представляем 1п(/ + £) в виде интеграла, который затем записывается через сумму:

1 _ Л-

1п (/ + 5) = 51 (I + аБ)'1 йа ~ (/ + «А) , (13)

О -1

где (/ + «,5) матрица обратная к (1 + а^) и « = (/ N. Выражение в правой части (13) является интефальной суммой этого оператора. Оно существует всегда, когда существует обратная матрица (/ + 5) . В то же время определитель матрицы (/+5), в случае выбора базиса из линейно независимых функций, положителен и не равен нулю, так как является определителем Грама [13]. Легко показать прямым вычислением, что если ряд (12) сходится, то уже при значении N равном 106 и Да = Ю-6 матричные элементы обоих представлений 1п(/ + 5') совпадают с точностью до 8-10 знаков после запятой. Таким образом, 1п(/ + 5) всегда определен и не зависит от величины интегралов перекрывания, если выбранный базис является системой линейно независимых функций.

Простой пример, иллюстрирующий изложенный выше способ устранения катастрофы неортогоналыюсти, приведен в Приложении I.

В разделе 2.4 найден вид произвольного оператора в базисе ортонормированных многочастичных функций. Систему

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

<Р = Фхехр —£> ), (14)

V 2 у

где Ч'и Ф матрицы-строки. Метод построения системы функций (14) называется симметричной ортогонилизацией [14], но сейчас проводится в операторном виде. Показано, что матричный элемент одночастичного и двухчастичного операторов на волновых функциях (14) может быть вычислен по формулам:

= (15)

Ят = ехр(^хЯ хехр|-Ае). (16)

Здесь

н=щ?)+{т\

и + - матричный элемент матрицы, обратной к матрице

(/ + 5).

В разделе 2.5 найдено представление оператора в виде ряда. В выражение (16) входит неэрмитов оператор Н и, кроме того, оператор О стоит в показателе экспоненты. Работать непосредственно с таким выражением сложно. Поэтому желательно записать его в виде достаточно быстро сходящегося ряда, все члены которого представлялись бы эршттовыми операторами. Показано, что оператор Н^ может быть представлен в виде:

где числа Е2 = -1, Е4 = 5, Е6 = - 61, Е% = 1385 , Ею = -50521, Еп = 2702765, Еи = —199360981 и т.д. являются числами Эйлера.

Для сходимости ряда (17) достаточно чтобы множество матричных элементов операторов [б, Н] было ограниченным. Из выражения для коэффициентов сп видна его быстрая сходимость. Например, уже коэффициент с7 к 1.39576хКГ7.

В разделе 2.6 найден вид оператора числа частиц в развиваемом формализме. Показано, что оператор числа частиц имеет свой обычный вид:

(18)

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

а<а*. —а*.ас а.а„ —а,.ас = —а+,.а* =0,

ь £ £ ь ' в 4 ъ £ £ ь '

а все вычисления проводились на так называемых перманентах, т.е. симметризованных относительно парных перестановок произведениях

одночастичных орбиталей. Тогда в обозначениях настоящей работы функции | {7}) должны иметь вид:

: >]"■:

где п, — числа заполнения орбиталей.

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

В разделе 3.1 рассматривается теория возмущений для случая, когда основной уровень нулевого приближения является вырожденным или квазивырожденным [16] в приложении к задачам, решаемым в диссертации. Если в качестве нулевого одночастичного базиса используются орбитали Хартри - Фока, то оператор возмущения записывают в виде:

1 Л' 1 .V

(19)

где оператор и (г,) выбирается из условия исключения параметров Рака /■""'•' для одноцентрового взаимодействия и матричных элементов прямого кулоновского взаимодействия для двухцентрового взаимодействия. В то же время, обычно используемое выражение и(г:) (см. например [15]) решает эту задачу лишь частично. В данном разделе предлагается оператор в представлении вторичного квантования, который удовлетворяет искомым условиям.

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

-{(йФог+ГхГа^Ше), (20)

где + +

(7^. ={£,\С\^) - амплитуда перехода электрона с лиганда на центральный ион, |Д{{}{4.-}| - энергия перехода системы из основного в возбужденное состояние.

В разделе 3.3 рассматривается следующий процесс третьего порядка. Электрон с лиганда переходит в валентную оболочку центрального иона, затем электростатическим полем возникающей дырки на лиганде, т.е. электрон-дырочным взаимодействием переводится в одну из вышележащих орбиталей \ф) и возвращается назад. Эффективный оператор У2, описывающий эти процессы, имеет вид: (ь\Кь\<Р)

1+

(21)

где - квантовые числа орбиталей валентной оболочки, <р — квантовые числа орбиталей пустых оболочек, О, С — квантовые числа орбиталей лиганда, Ьек — оператор электронно-дырочного взаимодействия. В диссертации отмечается также возможность построения операторов и диаграммным методом.

В разделе 3.4 рассмотрен следующий процесс возникновения спиновой плотности на лиганде. Электрон с лиганда переходит в вышележащие пустые оболочки |(р), \<р') центрального иона. Из-за обменного взаимодействия с

электронами валентной оболочки амплитуда перехода электрона со спином вверх и со спином вниз будет разная, что приводит к возникновению не равного нулю магнитного поля на ядре лиганда. Оператор Уъ, соответствующий этому процессу, имеет вид:

(£\у\0) + 1г.с. (22)

В разделе 3.5 рассмотрен следующий процесс третьего порядка теории возмущений. Соответствующий ему оператор обозначим как . Электрон с лиганда переходит в валентную оболочку центрального иона. Затем дырка на лиганде электронно-дырочным взаимодействием переводится на другую орбиталь и электрон возвращается назад. Этот процесс может быть сопоставлен с поляризацией лиганда:

Ы..........(23)

где — квантовые числа орбитален валентной оболочки, -

квантовые числа орбиталей лиганда, - оператор электронно-дырочного взаимодействия. Видно, что оператор У4 является поправкой следующего порядка к оператору ^. Форма оператора также может быть определена диаграммным методом.

В разделе 3.6 рассмотрен вклад в ЛСТВ, обусловленный процессами поляризации остова на центральном ионе. Влияние поляризации остова на сверхтонкую структуру примесного иона обсуждается уже достаточно давно [17,18]. Роль поляризации остова центрального иона в формировании лиганднои сверхтонкой структуры методом «приближенного вторичного

квантования» исследовалось в работе [19]. В работе [All] влияние поляризации остова центрального иона hei лигандную сверхтонкую структуру нами оценивалось следуя работе [19]. Затем в рамках метода развиваемого в данной диссертации, нами получен следующий оператор, учитывающий процессы поляризации остова [А12]:

2IVi

±{a\q\C){e\q\p)

{C\v\e)+h.c„ (24)

где — амплитуда перехода электрона с орбитали лиганда \С) на

освободившуюся орбиталь |«), принадлежащую основной конфигурации,

при условии, что на бывшей пустой орбитале | ф) находится электрон.

Видно, что полученные выражения существенно отличаются от тех, что получены ранее в расчетах с точностью до квадратов интегралов перекрывания. Главное отличие — перенормировка произведений различных двухцентровых интегралов и параметров ковалентности. Это обстоятельство особенно важно учитывать при определении параметров ковалентной связи с лигандами по экспериментальным данным об ЛСТВ.

В четвертой главе рассматриваются амплитуды перехода электронов между ионами, а также даны выражения для спиновых плотностей.

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

h = h> -\\\-Ъ\Г1Г\- Е ' (25)

lr-Ra| * lr-Ril .'«¿Ir-R/I

где hk - оператор кинетической энергии, Za,Zb - заряды ядер центрального иона и лиганда, соответственно, <7, - заряды ионов в беспримесном кристалле. Запишем заряд ядра иона как: Z,- =qi+ni + mii где я,- - число электронов i - ого иона для рассматриваемое! конфигурации, число -назовем зарядовым дефектом и оно будет отлично от нуля, если заряд внедряемого иона не совпадает с зарядом <7;. Тогда mi > 0 соответствует недостатку, а Щ < 0 избытку электронов. Для амплитуды перехода нами получено следующее выражение [All]:

lr_Ki.| 4i*S

"е+те

Ir-Rj

Г) + РШ)Ш\Пе)

г — К

(26)

где суммирование по орбиталям центрального нона отмечается индексом е, а по орбиталям лиганда индексом Ь, £\г и ^ энергии Хартри - Фока электрона на орбитали центрального иона и орбитали лиганда, соответственно. Величина:

где he - энергия кулоновского взаимодействия электрона, на орбитале |<?) центрального иона, с бесконечной кристаллической решеткой, взятой в ионном приближении, h^ - энергия на лиганде, пе+те,пь + ть - числа электронов и зарядовые дефекты в основной конфигурации на центральном ионе и на лиганде соответственно.

В разделе 4.2 получены выражения для амплитуд перехода электрона с лиганда в незаполненные и заполненные в основной конфигурации оболочки центрального иона. Выводы этих выражений аналогичны выводам амплитуды в валентную оболочку. В них (<p\G\£) = G(/%- - амплитуда перехода электрона из орбитали лиганда в орбиталь \<р), незаполненной в основной конфигурации оболочки, {cc\G9\^) = G^ — амплитуда перехода электрона в освободившуюся орбиталь заполненной оболочки.

В разделе 4.3 выражения для вкладов от рассмотренных процессов в константы суперсверхтонкого взаимодействия записаны через спиновые плотности [All], подобно тому, как это делается, например, в монографии [20]. Обозначим q\^) = q£( и в) = р.в. Тогда спиновые плотности могут быть записаны следующим образом.

Например, для ковалентного вклада, т.е. оператора V[, имеем: ДО 1 2 - -2

Jfi -^Л-РадУац+Гал, i = s,a,n^ (28а)

Л0 -1 1 _ 1_ _ _

Jfcm - 4 fall fx ~2Р*/°-/4fx + YafoY*fir ■ (286)

Обозначения s,ct,k соответствуют обозначениям, принятым в теории ковалентной связи [20]. Спиновые плотности, соответствующие процессам возбуждения электрона на центральном ионе, вызванные электрон —

дырочным взаимодействием и описываемые оператором К,, запишутся следующим образом:

И:

•'л I» I I = 5,<Х,7Г.

7(:> =

(1<77Т I .

1

^ РыРаА 2 РмУ*'

1 _

1+

УмУьр

4/.2р |

I _1 - I- _- -

. Р^ЛаР^/я ~ РьйггУЛ/ж Г, УыаР^/х УыяУА/Я

(29а)

1 + !

-ЧлМ

Л*) = 7(2) , 7(2)

J йети У ¿етт У </тгсг '

(296)

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

Вклад в спиновые плотности, обусловленный поляризацией лиганда, т.е. соответствующий оператору , записывается как:

7(2)

М _

Jfas ~ ]

1

1

^ Р4/аР4~ ^РА^УЛ " У4/<тР4 /з + У А /аУл j .

(30)

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

Таким образом, вклады от рассматриваемых процессов для всех ионов ряда редкоземельных элементов или ряда переходных металлов будут выражаться через введенные в настоящем разделе спиновые плотности.

В пятой главе получены выражения для вычисления матричных элементов операторов, входящих в амплитуды перехода электронов между ионами.

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

орбитали (г)) возьмем в виде разложения по гауссовскому типу

орбиталей (ОТО):

К/ = ®ф(-«,г2).

(31)

Такое разложение, как известно, является чрезвычайно удобным для расчета двухцентровых интегралов [21]. Однако, в отличие от метода расчета, описанного в работе [21], мы не используем предварительного преобразования оператора кулоновского взаимодействия к его фурье-образу с последующим интегральным представлением. В предлагаемом нами методе расчета формулы для двухцентровых матричных элементов получаются более компактными.

Пусть первый ион находится в узле с радиусом-вектором Л<)+Гу,И0=0. Второй ион - в узле с радиусом вектором гй. А в узле с

радиусом-вектором где - вектор «-ой элементарной ячейки

кристалла находится заряд др и Тр - радиус-вектор иона в элементарной ячейке. Тогда, как показано в работе [А9], матричный элемент вида:

где [и,/2] — целая часть от числа в скобках; хт = х0, х02 = у0, х03 = г0 — обозначают координата вектора г0 = гл — г;; а:к - « +ак, с. - /ал -обозначают компоненты вектора С; Я, -Я2= Ку, А, = - обозначают

может быть выражен через функции:

(32)

координаты вектора Н = Н„-(г,--гр). Например, матричный элемент, вычисленный на функциях |рг (г)), будет иметь вид:

Вычисления по формулам

(32) не требуют специального программирования и легко вычисляются в режиме пользователя, например в среде «Математика».

В разделе 5.2 получены выражения для вычисления одноцентровых и двухцентровых матричных элементов вида (27) и входящих в выражения для амплитуд перехода [А13]. Определяются матричные элементы оператора

11 ,.п

- Г/ (г = («■ - Г/)| - ^(г СЗЗ)

где сумма по п - есть сумма по элементарным ячейкам бесконечной кристаллической решетки, а сумма по р - по ионам в элементарной ячейке. В отличие от метода Эвальда [22] в конечных выражениях суммирование проводится только по векторам обратной решетки.

Показано: что матричные элементы вида (33) могут быть выражены через функции:

^л-2 ( 1 («1«2»3 ) = - "-«1!'Ь!"3!Х Я А —

ус V ал

3

^ ( и ^

2 I аА -.2

ехр| -

а,к у

g V У 8

где ёз) являются полиномами п - порядка от компонент gs вектора

обратной решетки g. Например,

где ZJ = Ркх0, / а& - igs / 2ал. Вегагчина (ё) определяется как структурный фактор, зависящий только от радиуса вектора Г. и вектора обратной решетки.

Таким образом, могут быть вычислены все матричные элементы оператора II1Я для е-, р-, (1- и Г-орбиталей ионов.

В разделе 5.3 получены выражения для одноцентровых матричных элементов оператора Я„, для е-, р-, с1-электронов выделенного иона [А14]. Матричные элементы такого типа возникают при вычислении амплитуд перехода в 5р- и 5с1-оболочки редкоземельных ионов, а также при переходах электронов анион-катион в ЬаМп03. В случае в-электронов матричный элемент является энергией взаимодействия я-электрона и определяется формулой:

1

V«*

е,(8) '

^ 8*0 g

V 4«*

+ 2д,

(35)

где vc - объем элементарной ячейки, а,Ь,с - постоянные решетки, <7,^) -структурный фактор, — заряд замещаемого иона. Например, в случае р-электронов недиагональный матричный элемент на функциях имеет

вид:

V СЛ8) '

а

4 ая

(36)

В разделе 5.4 получены выражения для матричных элементов оператора 11 и,, входящего в амплитуду перехода электрона с аниона в 4£-оболочку редкоземельного иона [А 16]. Например, матричный элемент:

(/,2|ЯМ |/,2) = (402) +К (042)+ 2^. (222)] + Зд,5>,

16

1

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

В Приложении III получены формулы для вычисления интегралов перекрывания: прямых, гибридных и обменных кулоновских интегралов, используя результаты Приложения II. Все вычисления проводятся в режиме пользователя, например в среде «Математика».

В шестой главе вычислен тензор ЛСТВ примесных центров УЬ3+: СэСаБз, С82№УР6 и пО: ЬаМпОэ.

В разделе 6.1 общие выражения для амплитуды перехода электрона в валентную оболочку конкретизируются для примесного центра УЪ3+:С5СаР3. Приведены численные значения матричных элементов, необходимых для вычисления этих амплитуд. Для самих амплитуд получены следующие значения:

(4/0[0|25) = 0.01465, (4/0|С|2/?0) = 0.01738, (4/1|С|2р1) = -0.007943

В разделе 6.2 общие выражения для амплитуд перехода электрона в 5(3-и бэ-оболочки конкретизируются для примесного центра УЪ :СзСаР3. Приведены численные значения матричных элементов необходимых для вычисления этих амплитуд. Для самих амплитуд получены следующие значения:

<5£/0|С?[2^) = - 0.2293, (5с?0|С|2р0) = - 0.1133, (5с1\\0\2р\) = 0.06164, (6* | С12^} = - 0.4642, (б* | в 12р0) = - 0.1293.

В разделе 6.3 общие выражения для амплитуд перехода электрона в 5э-и 5р-оболочки конкретизируются для примесного центра УЬ3,:СяСаР3. Приведены численные значения матричных элементов, необходимых для вычисления этих амплитуд. Вычисления для амплитуд (7^., приводят к

значениям: я «0.003, как минимум на порядок меньше остальных амплитуд перехода, и для амплитуд были получены следующие значения: =0, С& =-0.0638, = -0.0581.

В полуэмпирических расчетах практически всегда включают 6р-оболочку. Однако, для включения этих орбиталей в схему теории возмущений, возможно недостаточно приближений, которые используются выше для вычисления амплитуд перехода. Перекрывание орбиталей 6р-оболочки с орбиталями лиганда достаточно большое, поэтому необходимо оценить слагаемые, отброшенные при выводе выражений для амплитуд перехода. Для оценки порядка величины и знака вкладов от бр-оболочки с базисом настоящей работы и в обозначениях для спиновых плотностей, принято: д(1р; = = у,ч>< = т.е. значения равные интегралам перекрывания, взятые с соответствующим знаком. Проводя вычисления, аналогичные вычислениям для амплитуд С", для амплитуды перехода (7^ получим:

«0.0662.3, С^«-0.026487.

В разделе 6.4 вычислены параметры ковалентности для примесных центров УЪ^СвСаРз, СвгКаУБб. В случае неизовалентного замещения, т.е. для примесных центров УЬ3+:СзСаРз. для параметров ковалентности валентной оболочки получены следующие значения:

П/з = - 0.0134, уЧа = - 0.0653, = 0.0298.

Для примесного центра УЬ'^СБгКаУРс,, т.е. изовалентного замещения, зарядовые дефекты тс,ть в основной конфигурации равны нулю, в отличие от то же время в возбужденной конфигурации возникает

взаимодействие электрон-дырка [23]. Там же отмечается, что оно заметно зависит от реального распределения электронной плотности. Приведем для сравнения несколько матричных элементов от этих взаимодействий:

(2511 /|г — Ке|12.у) = 0.24051, 2%|2я) =0.24053 ,

к

(4/0|1/|г-Н.||2*) = -0.00372, <£,■4/0|2|&,2*) = -0.00373,

где = нижний крамерсов дублет. Для остальных матричных

элементов выполняются аналогичные соотношения. Видно, что замена левого столбца на правый, при переходе от

УЬ .СвСаРз к УЪ3 :С82№УР6, практически не изменяет значение амплитуд перехода. Незначительное изменение происходит так же из-за сумм матричных элементов оператора Н,к. Таким образом, показано, что параметры ковалентности в рамках сделанных приближений, могут бьпъ взятыми одинаковыми для УЪ3+:С5СаР3 и УЬ3 :С52НаУР6. Для остальных: параметров ковалентности получены следующие значения:

7$^ = 0.152, уиа = 0.167, 0.0906, г6„==0.290, ^ = 0.165.

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

В разделе 6.5, используя формализм спинового гамильтониана, проведено сравнение с экспериментом для примесных центров УЬ3+:С8СаР3, СэгНаУРб. В случае кубических примесных центров, в которых примесным ионом является ион УЪ3+, гамильтониан лигандного сверхтонкого взаимодействия запишется следующим образом:

Я« = 7-5./, + 71 (V, + ЗД- (37)

Переход от операторов V,, в представлении вторичного квантования к спиновому гамильтониану, проводится стандартным образом [8,20]. Волновая функция дублета |ГЙ) имеет следующий вид:

= |4/0,+) + ^|4/1,-Н-] |4/ — 3,—). (38)

Приведем выражения для вкладов в компоненты тензора ЛСТВ от некоторых процессов. Вклад от эффектов неортогональности и ковалентности 4{'-оболочки будет определяться формулами [А 11]:

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

[All]:

^=f /14 if}

rW=_2 .(2) (2 {2) 3 и'4) 13 /з (2)

2Jds s I 3 ■>da ryJd* \ P g y 2 J*fa* Рш

В экспериментальных работах обычно приводится изотропная часть тензора ЛСТВ Л =(^ + 2rL)/3 и анизотропная часть Ар = (^-Г,)/3. В

табл.4 даны вклады от всех рассмотренных процессов, а также экспериментальные данные.

Таблица 4

Теорепгческие величины ЛСТВ параметров As и Ар (в MHz) первой координационной сферы F" для Yb3+ в CsCaF3 и Cs2NaYF6. Экспериментальные данные взяты из работы [All]

Кристалл Ad V, V, v3 v4 V, Сумма Экспер.

CsCaF3 As 0 9.88 4.34 0.68 0 5.18 20.08 20.577

Cs2NaYF6 0.9 20.98 22.111

CsCaF3 AP -9.0 2.89 1.38 -0.95 0 0.62 -5.06 - 5.963

Cs2NaYF6 -0.8 -5.86 - 6.208

В разделе 6.5 для кристалла ЬаМп03 вычислена изотропная компонента тензора ЛСТВ иона пО для комнатной температуры. Проведено сравнение с экспериментально наблюдаемыми сдвигами линии ЯМР. Получено достаточно хорошее согласие с экспериментом [А15,А17]. Таким образом, по данным ЯМР была определена картина кооперативного упрядочения в этом соединении.

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

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

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

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

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

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

позволяют объяснить основные особенности формирования локальных полей

3+ 3+

на ядрах фтора в примесных центрах: Yb в CsCaF3 и Yb в Cs^NaYF^.

17

6. Рассчитаны величины локальных магнитных полей на ядрах О в

3+

LaMnO^, при различном упорядочении орбиталей ионов Ми . Путем

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

Список публикаций по теме диссертации

[AI] Anikeenok, O.A. ENDOR and transferred spin densities of the 4f" ions in fluorides / O.A.Anikeenok, M.V.Eremin, M.L. Falin,A.L. et al. // J. Phys. C: Solid State Phys. - 1984. - 17. - N.15. - P.2S13-2823.

[A2] Еремин, M.B. Косвенное взаимодействие 4^электронов с лигандами через заполненные 5р-оболочки./М.В.Еремин, A.A. Каминский, O.A. Аникеенок//ФТТ. - 1985. - Т.27. - №2. - С.455-458.

[A3] Аникеенок, O.A. Операторная техника в теории взаимодействия редкоземельных ионов с ядрами лигандов / О.А.Аникеенок, М.В.Еремин,О.Г.Хуцишвили // ФТТ. - 1986. - Т.28. - №6. - С.1690-1697.

[A4] Аникеенок, O.A. ДЭЯР и спиновые плотности на лигандах 4f°-ионов во фторидах /О.А.Аникеенок, И.Р.Ибрагимов, В.А.Уланов, М.Л.Фалин // ФТТ. - 1986. - Т.28. - №3. - С.821-826.

[А5] Anikeenok, O.A. Derealization of Cu2+ unpaired electronon the next nearest ligands in LaSrGa0.995Cul).oo504 single crystal / O.A.Anikeenok, M.A.Augustyniak-Jablokow, T.A.Ivanova, P.Reiche, R.Uecker, Yu.V.Yablokov //Phys. stat. sol. (b). - 2001. - V.226/ - №1. - R1-R3.

[A6] Anikeenok, O.A. Supertransferred hyperfine interactions in layerLaSrGao.995Cuo.oo504 / O.A.Anikeenok, M.A.Augustyniak-Jablokow, T.A.Ivanova, P.Reiche, R.Uecker, Yu.V.Yablokov // Physica B. - 2003 - V.325. -P.246-255.

[А7] Уланов, В.А. Электронная структура нецентрального комплекса двухвалентной меди в кристалле SrF2 по данным ЭПР и ДЭЯР / В.А.Уланов, О.А.Аникеенок, М.М.Зарипов, И.И.Фазлижанов // ФТТ. - 2003. - Т.45. -№10. - С.1814-1817.

[А8] Аникеенок, О.А. Вычисление из первых принципов сверхтонких полей на лигандах во фторидах / О.А.Аникеенок // ФТТ. — 2003. - Т.45. — №5. — С.812-816.

[А9] Аникеенок, О.А. Кристаллическое поле на примесных центрах в ионных кристаллах / О.А.Аникеенок // ФТТ. - 2005. - Т.47. - №6. - С. 10651070.

[А 10] Аникеенок, О.А. Вычисление из первых принципов амплитуд перехода электрона с лиганда в 5с1-оболочку Yb3+:KZnF3 / О.А.Аникеенок // ФТТ. - 2006 - Т.48. - №10. - С. 1771-1776.

[All] Falin, M.L. Transferred hyperfme interactions forYb3+ ions in CsCaF3 andCs2NaYF6 single crystals: Experimental and ab initio study / M.L.Falin, O.A.Anikeenok, V.A.Latypov, N.M.Kliaidukov, F.Callens, H.Vrielinck, A.Hoefstaetter // Phys. Rev. B. - 2009. - V.80. - №17. - P.l74110 (1-11).

[A12] Аникеенок, О.А. Лигандное сверхтонкое взаимодействие Yb3+ в кристаллах CsCaF2 и Cs2NaYF6 / О.А.Аникеенок // ФТТ. - 2011. - Т.53, -№11. - С.2209-2215.

[А13] Anikeenok, О.А. Approach to calculation of long-range Coulomb interactionmatrix elements in ion crystals / O.A.Anikeenok // Magn. Resonance in Solids.EJ. - 2011. - V. 13. - №2. - P.27-35.

[A14] Аникеенок, О.А. Далънодействующее кулоновское взаимодействие в ионных кристаллах / О.А. Аникеенок // ФТТ. — 2012. — Т.54. -№9. - С.1733-1738.

[А 15] Trokiner, A. S. Melting of the orbital order in ЬаМпОЗ probed by NMR / A.Trokiner, S.Verkhovskii, A.Gerashenko, Z.VoIkova, O.Anikeenok, K. Mikhalev, M. Eremin, L. Pinsard-Gaudard. // Phys. Rev. B. - 2013. - V.87. -№12. — P.125142(1-6).

[A16] Аникеенок, О.А. Длльнодействующее кулоновское взаимодействие электронов 4f-op6irraneft в примесных центрах Yb3+:KZnF3, CsCaF3 и Sm3+:CaF2 / О.А. Аникеенок // ФТТ. - 2013. - Т.55. - №11. - С.2190-2195.

[А17] Anikeenok, O.A. Transferred hyperfine interactions for 017:LaMn03/ OA.Anikeenok // Magn. Resonance in Solids. EJ. - 2014. - V.16. - №1. -P.14101(1-7).

Список цитированной литературы

[1] Фок, В.А. Приближенный способ решения квантовой задачи многих тел /В.А.Фок // УФН. - 1967. - Т.93. - №.10. - С.342-361.

[2] Froese, Ch. The orthogonality assumption in the Hartree-Fock approximation / Ch.Froese // Canad. J. Phys. - 1967. - V.45, no.l, P.7-12.

[3] Wilson, R.G. The third order Hartree-Fock energy for the ls2s states of helium like ions / R.G.Wilson, I.Eliezer // J. Phys. B. - 1971. - V.4. - N10. -P.1255-1256.

[4] Moshinsky, M. Group theory and second quantization for nonorthogonal orbitals / M.Moshinsky, T.H.Seligman // Annals of Phys. - 1971. - V.66. - N.l. -P.311-334.

[5] Artacho, E. Nonorthogonal basis sets in quantum mechanics: representation and second quantization / E.Artacho, L.M.Bosch // Phys. Rev. -1991. -V.43. - N.l 1. -P.5770-5777.

[6] Еремин, M. В. Техника вторичного квантования в обобщенном методе Гайтлера-Лондона / М. В. Еремин, А. М. Леушин // ФТТ. - 1974. -Т.16. №7. -С. 1917-1923

[7] Eremin, M.V. The Superposition Model in Crystal Field Theory/ M.V. Eremin, A.A. Kornienko //phys stat sol (b) . - 1977. - V. 79. -N.2. - P.775-785.

[8] Anikeenok, O.A. ENDOR and transferred spin densities of the 4f13 ions in fluorides / O.A.Anikeenok, M.V.Eremin, M.L.Falin, V.P.Meiklyar // J. Phys. C: Solid State Phys. - 1982. - V.15. -N.7. - P. 1557-1567.

[9] Eremin, M.V. Charge transfer process contribution to the zero-Field splitting of the S-state transition ions / M.V.Eremin, I.I.Antonova // J. Phys.: Condens. Matter . - 1998. - V.l 0. - N.25. - P.5567-5575.

[10] Аникеенок, O.A. Теория электронно-ядерных взаимодействий парамагнитных ионов с лигандами при отсутствии а- или я-связей / О.А.Аникеенок, М.В.Еремин // ФТТ. - 1981. - Т.23. - №3. - С.706-713.

[11] Anikeenok, О.А. The peculiarities of the transferred hyperfine interaction in CaF2:Ti2+, V2+ and SrF2: Ti2+ / O.A.Anikeenok, M.V.Eremin, M.L.Falin, V.P.Meiklyar // J. Phys. C: Solid State Phys. - 1982. - V.15. - N.5. -L105-L107.

[12] Baker, J.M. Fluorine electron nuclear double resonance of divalent thulium in calcium fluoride / J.M. Baker, J.P.Hurrel // Proc. Phys. Soc. - 1963. -V. 82.-N.3.-P. 742-756.

[13] Шилов, Г.Е. Конечномерные линейные пространства / Г.Е.Шилов. - М.: Наука, 1969. - 432 с.

[14] Lowdin, Р.-О. On the nonorthogonality problem / P.-O.Lowdin // Adv. Quantum. Chemistry. - 1970. - V.5. - N. 1. - P. 185-199.

[15] Judd, B.R. Second Quantization and Atomic Spectroscopy. - Baltimore: The Johns Hopkins Press, 1967. - 210 p.

[16] Бир, Г.Л. Симметрия и деформационные эффекты в полупроводниках / Г.Л. Бир, Г.Е. Пикус. - М.: Наука, 1972. - 584 с.

[17] Watson, R.E. Hyperfine Interactions / R.E.Watson, A.J.Freeman. -New York, 1967,- P.53.

[18] Moser, C.M. Hyperfine Interactions / C.M.Moser. - New York, 1967. -

P.95.

[19] Falin, M.L. ENDOR and Transferred Hyperfine Interaction of Impurity Rare-Earth Ions with Nearest Diamagnetic Ions in Crystals / M.L.Falin, M.V.Eremin, H.Bill, D.Lovy // Appl. Magn. Reson. - 1995. - V.9. - N.3. - P.329-354.

[20] Абрагам, А. Электронный парамагнитный резонанс переходных ионов / А. А брагам, Б.Блини. - М.: Мир, 1973. - Т.1. - 651 е.; Т. 2. - 349 с.

[21] Taketa, Н. Gaussian-expansion methods for molecular integrals / H.Taketa, S.Huzinaga, K.Ohata / J. Phys. Soc. Japan. - 1966. - V.21. - N.ll. -P.2313 - 2324.

[22] Ewald, P.P. Die berechnung optischer und elektrostatischer gitterpotentiale / P.P. Ewald // Ann. der Physik - 1921. - V.64. - N.3. - P.253-287.

[23] Маделунг, О. Теория твердого тела / О.Маделунг. - М.: Наука, 1980.-416 с.

Подписано в печать 06.01.2015. Bjniara офсетная. Печать цифровая. Формат 60x84 1/16. Гарнитура «Times New Roman». Тираж 150 экз. Заказ 42/2.

Отпечатано с готового оригинал-макета в типографии Издательства Казанского университета

420008, г. Казань, ул. Профессора Нужина, 1/37 тел. (843) 233-73-59, 233-73-28