Конспект лекций Казань 2011

ФГАОУВПО
«КАЗАНСКИЙ (ПРИВОЛЖСКИЙ) ФЕДЕРАЛЬНЫЙ УНИВЕРСИТЕТ»
ИНСТИТУТ ФИЗИКИ

М. В. Еремин

Микроскопические модели в

конденсированных средах

Конспект лекций

Казань 2011

УДК 538.9
Печатается по решению Редакционно-издательского совета ФГАОУВПО

«Казанский (Приволжский) федеральный университет»
Учебно-методического совета Института физики

Протокол № 3 от 28 апреля 2011 г.
заседания кафедры Квантовой электроники и радиоспектроскопии

Протокол № 9 от 13 апреля 2011 г.

Рецензенты:

доктор физ.-мат. наук, проф. КФУ Б.З. Малкин

доктор физ.-мат. наук, зав. лаб. КФТИ КазНЦ РАН В. Ф. Тарасов

Еремин М. В.

Микроскопические модели в конденсированных средах: Учебное пособие / М.В. – Казань: Казанский (Приволжский) федеральный университет, 2011. – 113с.

Учебное пособие является конспектом лекций по курсу «Микроскопические модели в конденсированных средах», который читался автором для студентов-магистров первого года обучения в 2010-2011 годах. Данный курс входит в группу «Ф» профессионального цикла подготовки магистров по направлению 011200.68 - «Физика» и является необходимым для изучения в рамках магистерской программы «Современные проблемы физики».

Еремин М. В., 2011

Содержание

1. Метод канонических преобразований

в теории возмущений……………………………..5

2. Суперобменное взаимодействие

2.1. Кинетический антиферромагнитный обмен………..13

2.2. Кинетический ферромагнитный обмен……………..17

2.3. Косвенное суперобменное взаимодействие.

Квантовая интерференция …………………………….…22

2.4. Анизотропное симметричное обменное

взаимодействие…… ……………………………..…….….27

2.5. Антисимметричное обменное взаимодействие…..…30

2.6. Биквадратичное обменное взаимодействие……….. .36

2.7.Двойной обмен (Double exchange) ……………………38

2.8.Обменное взаимодействие ионов смешанной

валентности.………………………………………………...40

2.9. Обменное взаимодействие ионов в кристаллах

типа KMeF₃…………………..……………………….……43

2.10. Обменно-дипольные переходы ……………………..46

3. Взаимодействие электронов проводимости через поле фононов………………………………….………… .. 49

3.1 Гамильтониан модели Бардина,

Купера, Шриффера …………………………………….. …54

3.2. Уравнение для параметра порядка при переходе в

состояние с волнами зарядовой плотности ………………56

4. Модели перехода металл- диэлектрик

4.1. Модель Пайерлса…………………………………….. ..57

4.2. Переход Мотта…… …………………… ..…….…63

4.3. Модель Хаббарда…… …………… ………..….....65

4.4. Метод функций Грина… ………… ………………69

5. Сверхпроводники.

5.1. Суперобменный механизм спаривания

носителей тока………………………………………………74

5.2. Спиновая восприимчивость сверхпроводников……...81

5.3. Оператор сверхпроводящего тока ……………….....…87

5.4. Оператор диамагнитной компоненты тока………….. 90

5.5. Среднее значение парамагнитного тока……………... 91

5.6. Среднее значение диамагнитной компоненты тока… 94

6. Функция диэлектрической проницаемости

6.1. Продольная и поперечная проницаемости.

Оптическая проводимость…….….…………………… … 96

6.2. Плазмоны…………………………….………………...103

6.3. Перенормировка фононных мод. Экранирование

кулоновского взаимодействия…….………….…..………106

Литература ………..……………………………………..…. 112

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

( феноменологических) спиновых гамильтонианов.

Рис. 1 Схема уровней энергии.

Блочный вид матрицы гамильтониана

выглядит следующим образом:

. (1.1)
Предполагается, что .

Вспомним классическое определение: «переход к другой “системе координат” путем преобразования подобия для операторов и одновременной заменой волновых функций на называется каноническим преобразованием» [1.1] ».

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

. (1.4)

При этом, как правило, интересуются лишь видом блока . Условие унитарности оператора преобразования

(1.5)
сводится к условию . т. е. матрица S - антиэрмитова. Для обеспечения эрмитовости из оператора S иногда явно выделяют мнимую единицу, т. е. пишут S=iL. Ниже мы используем обозначения как в книге [1.2], в которой дается наиболее подробное изложение данного метода. Ниже, в учебных целях, дадим упрощенное изложение метода. В книге [1.3] данная тема фигурирует под названием “ Операторная форма теории возмущений”. Там терминология и обозначения несколько другие. Эффекты, связанные с возбужденными конфигурациями, “проектируются” на состояния основной конфигурации. Суть же метода та же самая.

Итак, подставляя (1.5) в (1.2), находим

(1. 6)
План дальнейших преобразований состоит в следующем. Матрица H разбивается на блочно-диагональную и недиагональную части . Матрица H⁽²⁾ не содержит диагональных матричных элементов по состояниям электронных конфигураций. Матрица содержит только блочно - диагональную часть возмущения.

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

(1.7)

Матрицу S будем искать методом последовательных приближений в виде . В (1.7) группируем члены одного порядка малости:

(1.8)

Диагональная часть оператора должна содержать коммутаторы и с четным числом операторов S и коммутаторы с нечетным числом операторов S. Недиагональная же часть правой части (1.8) (та, что должна обратиться в ноль) должна содержать коммутаторы и с нечетным числом операторов S и коммутаторы с четным числом операторов S . Учитывая это обстоятельство, условие обращения в нуль недиагонального блока записывается в виде цепочки равенств. Число уравнений соответствует рассматриваемому порядку теории возмущений.

(1.9)
(1.10)
(1.11)
Уравнение (1.11) с учетом (1.9) можно переписать в виде
. (1.12)
Учитывая, что , находим правило вычисления матричных элементов оператора S
(1.13)
Из (1.8) находим выражение для матричных элементов эффективного гамильтониана
(1.14)
Для упрощения (1.14) используем соотношения (1.9) и (1.10)

(1.15)
Нередко оказывается так, что энергетические разности велики по сравнению с расщеплениями состояний основной и возбужденной конфигураций. Тогда энергетические знаменатели можно считать не зависящими от индекса m , т.е. . В этом случае поправочные члены к эффективному гамильтониану записываются компактнее:

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

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

(1.17)
Отсюда, после использования формулы (1.13), во втором приближении теории возмущений имеем:
(1.18)

В дальнейшем нам потребуется вид эффективного гамильтониана с точностью до шестого порядка теории возмущений. Дополнительные к (1.13) уравнения для определения матричных элементов оператора S имеют вид:

(1.19)
Общая формула для эффективного гамильтониана такова:
(1.20)
Задача. Найти поправки к оператору (1.18) в третьем порядке теории возмущений.

Суперобменное взаимодействие

2.1. Кинетический антиферромагнитный обмен

Для пояснения основной идеи, следуя работе Anderson P.W. Phys. Rev, 115, 2-13 (1959), рассмотрим два магнитных иона в точках (a) и (b). На каждом из центров находится по одному спину. Каждый из спинов имеет два состояния – вверх и вниз. Так что основное состояние пары ионов четырёхкратно вырождено. Найдем вид эффективного спинового оператора, обусловленного перескоками электрона с одного центра на другой. Схема перескоков электрона показана на Рис. 2.1

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

Рис. 2.1 . Схема перескоков электрона с центра а на центр b и обратно.

. (2.1)
Используя (1.13), пишем

(2.2)
Подставляя их в эффективный оператор,
(2.3)
получаем

. (2.4)

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

Итак, в единицах матрица имеет вид:
(2.5)

Легко проверить, что все матричные элементы могут быть воспроизведены следующим спиновым оператором

(2.6)
Следовательно, эффективный оператор обменного взаимодействия имеет вид:
(2.7)
Словесное пояснение смысла квадратной скобки таково. Из-за принципа Паули, переброс спина с одного узла на другой возможен лишь при антипараллельной их ориентации. В этом случае скалярное произведение спинов отрицательно. Энергия системы спинов понижается. При параллельной ориентации спинов перескоки запрещены. Квадратная скобка обращается в ноль.

Обобщение теории Андерсона на случай нескольких электронов у взаимодействующих атомов рассматривалось в ряде работ. Хороший обзор полученных результатов в этом направлении приведен в книге [2.2]. Сжатое изложение многоэлектронной теории кинетического обменного взаимодействия в кристаллах с учетом, как спиновых, так и орбитальных моментов взаимодействующих ионов дано в [2.3]. Там же проведено сравнение расчётов с экспериментальными данными для пар ионов, находящихся в основных и возбужденных состоянияx. Четкая корреляция силы обменного взаимодействия с величиной энергетического знаменателя в убедительно свидетельствует, что механизм Андерсона является доминирующим механизмом обменного взаимодействия в соединениях с незаполненными 3d- оболочками.

2.2. Кинетический ферромагнитный обмен
Встречаются случаи, когда интегралы перескока между полузаполенными орбиталями магнитных ионов равны нулю. Пример пара ионов Ni²⁺ и V²⁺ в .

Рис.2.2. Схематическое распределение волновых функций для спинов, занимающих состояния и . Знаки соответствуют фазам волновых функций.

Рассмотрим процесс, связанный с перескоком электрона из заполненной оболочки иона а в незаполненную оболочку иона b.
(2.8)

Здесь индексы и относится к заполненной орбитали иона а.

(2.9)

Рис. 2.3. Схема перескоков электрона из заполненного состояния иона a в полузаполненное состояние иона b. Двойная пунктирная линия соответствует оператору кулоновского взаимодействия электронов.

(2.10)
Отметим, что первая поправка в данном случае обращается в ноль. Привлекаем к рассмотрению
(2.11)
(2.12)
Вычисляем коммутатор

(2.13)

Затем запишем произведение операторов, входящих в коммутатор

(2.14)
Здесь индекс относится к заполненной орбитали иона а. Оператор предназначен для расчетов в основном состоянии. Если оператор переставить налево, то соответствующее выражение должно обратиться в ноль, так как орбиталь заполнена. Там нет пустых мест. Это соображение позволяет

провести упрощения.

(2.15)
Так как оператор симметричен относительно перестановки индексов, можем провести смену номеров электронов

(2.16)
В первой сумме правила отбора по спиновым квантовым числам таковы

(2.17)

Этот оператор может быть воспроизведен в представлении спиновых функций оператором

(2.18)

Во второй сумме правила отбора иные

Поэтому эта сумма сводится к числу .

В итоге имеем

(2.19)

Здесь дополнительно учтено, что для эквивалентных электронов

. В (2.19) вошел обменный интеграл от кулоновского взаимодействия электронов. В теории свободного атома именно такого рода интегралы приводят к тому, что состояния атомов с максимальным спином имеют минимальную энергию. Учитывая это обстоятельство, механизм ферромагнитной связи можно интерпретировать следующим образом. Перескок электрона из заполненной орбитали осуществляется так, чтобы: 1) не войти в противоречие с принципом Паули на центре в , 2) при этом промежуточное состояние иона а должно соответствовать правилу Хунда, т. е. процесс перескока будет соответствовать минимальной энергии, если спины электронов на центре a будут параллельными.
Задача. Соединение EuO является ферромагнетиком. Температура Кюри составляет 69.4 K. Основная электронная конфигурация у иона Eu²⁺ - 4f⁷ . Нижайшая из возбужденных конфигураций -4f⁶5d. Интеграл перескока электрона из 4f- оболочки одного иона в незаполненную 5d оболочку другого много больше, чем в 4f –оболочку. В этой связи естественно ожидать, что процессы виртуального переноса электрона из 4f- оболочки одного иона в незаполненную 5d-оболочку другого являются доминирующими в механизме обменного взаимодействия между ионами Eu. Найти вид эффективного спинового гамильтониана, соответствующего такому механизму. Получить формулу для параметра обменного взаимодействия через интегралы перескока и обменной интеграл кулоновского взаимодействия 4f- и 5d- электронов, находящихся на одном и том же узле.
2.3. Косвенное суперобменное взаимодействие. Квантовая интерференция
Максимально упрощенная модель три центра – четыре электрона изображена на рис. 2.4

Рис. 2. 4. Схема каскадных перескоков электрона с иона a на ион b через диамагнитный ион с (механизм Андерсона). Греческие буквы обозначают состояния электронов на соответствующих ионах. Цифры и стрелки соответствуют последовательности и направлению перескоков электрона.

Рис. 2.5. Схема возбуждения, соответствующая

каскадным перескокам на рис. 2.4.

Матричные элементы эффективного оператора определяются следующим образом:

(2.20)
Здесь много слагаемых, но для описания суперобменного взаимодействия нужны только два

(2.21)

Здесь роль выполняет оператор

. (2.22)
Матричные элементы определяются оператором

(2.23)

Учитывая все это, для приведенной на рис.2.5 схемы возбуждения, имеем:

(2.24)

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

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

(2.27)

Оператор суперобменного взаимодействия приобретает вид:

, (2.28)

где

(2.29)
Обратим внимание, что по сравнению с формулой (2.8), здесь имеется новое качество. Для пояснения его рассмотрим случай, когда в роли промежуточного иона имеется не один, а два иона с₁ и с₂. Варианты возможных перескоков электронов иллюстрируются на рис.2.6 (ниже). В тех случаях, когда произведение интегралов перескока , каскадные процессы перескока электронов способствуют антиферромагнитной связи. Если же , то данный процесс способствует ферромагнитной связи.

Рис. 2.6. Схематическое изображение двухмостиковых механизмов суперобменного взаимодействия между магнитными ионами а и b через диамагнитные ионы с₁ и с₂ . Цифрами нумеруются последовательности перескоков электрона с одного центра на другой.
Здесь в механизме суперобмена особенно ярко проявляются процессы квантовой интерференции. Перебирая возможные варианты распределения электронной плотности основных состояний магнитных ионов, можно убедиться, что когда состояния магнитных ионов ортогональны друг другу, суммарный эффект изображенных виртуальных перескоков обращается в ноль. Антиферромагнитный характер суперобмена полностью погашен из-за квантовой интерференции. Кстати сказать, нетрудно понять, что аналогичное явление имеет место и в наноконтактах, соединенных двумя проводниками (мостиками).
Задача для самостоятельного рассмотрения. Составить схему процессов переноса электрона, соответствующую пятому порядку теории возмущений. Считать, что один из спинов находится в состоянии , а другой в состоянии . Углы связи a-с₁-b и a-с₂-b равны 90^о. Линия связи магнитных ионов совпадает с осью x. Схема распределения электронной плотности изображена на рис. 2.7

Рис.2.7. Два варианта распределения электронной плотности на магнитных ионах. Левый соответствует

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

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

2.4. Анизотропное обменное взаимодействие
Рассмотрим вначале упрощенную модель, состоящую из двух центров и двух электронов [ B. Bleaney and K. D. Bowers, Proc. R. Soc. A 214, 451 (1952)]

Рис. 2.9. Схема виртуальных процессов, соответствующая возникновению анизотропного обменного взаимодействия. Пунктирные стрелки соответствуют спин- орбитальному взаимодействию.

Пусть между основным и возбужденным состояниями иона a имеется матричный элемент от оператора спин-орбитального взаимодействия

. (2.30)

Проанализируем операторный вид поправок третьего порядка теории возмущений

(2.31)

Зависящая от операторов спинов часть имеет вид . В этой связи рассмотрим произведение операторов (2.32)

Используем свойства матриц Паули , и т. д. и получаем

(2.33)
Члены, содержащие скалярное произведение спинов и линейные по спину , можно опустить. Первые дают малую поправку к изотропному обменному взаимодействию , а последние сократятся, так как линейные члены по спину в гамильтониане в отсутствие магнитного поля фигурировать не могут. Это следствие инвариантности уравнения Шредингера относительно операции обращения знака времени.

Таким образом, спин – гамильтониан, соответствующий рассматриваемой поправке, может быть записан в виде:

(2.34)

Рис. 2.10. Схемы симметричного обменного взаимодействия между спинами в состояниях и через возбужденные состояния и Штриховые стрелки соответствуют переходам, обусловленным спин-орбитальным взаимодействием. Процесс номер (v), более детально проиллюстрирован внизу рисунка, с явным указанием соответствующих интегралов перескока и матричных элементов орбитального момента, входящего в спин-орбитальное взаимодействие.

2.5. Антисимметричное обменное взаимодействие.
Часто его называют еще как взаимодействие Дзялошинского-Мория. Эффективный оператор имеет вид
. (2.35)

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

Для вывода микроскопического выражения для вектора рассмотрим вначале произведение операторов
(2.36)
Первый член, будучи линейным по спину, не дает вклада в эффективный гамильтониан. Второй же по векторной структуре совпадает с оператором антисимметричного обмена. Это обстоятельство и было отмечено Мория. Согласно Мория, оператор должен быть пропорционален матричным элементам от оператора спин-орбитального взаимодействия.

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

Рассмотрим модель из трех центров и четырех электронов. На рис. буквой c обозначен диамагнитный ион, выступающий в качестве мостика между магнитными ионами a и b.

Рис. 2.11. Правило определения направления вектора Дзялошинского – Мория при взаимодействии магнитных ионов a и b через диамагнитный ион с, выведенное в работе [А.С. Москвин, И. Г. Бострем, ФТТ, Том. 19, стр. 1616 ( 1977) ]

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

Рис. 2.12. Схемы возникновения взаимодействия Дзялошинского- Мория между спинами в состояниях и на ионах a и b соответственно через возбужденные орбитальные состояния.

(2.37)

Далее воспользуемся выражением для эффективного оператора суперобменного взаимодействия

(2.38)
заменив в нем интегралы перескока по схеме
(2.39)

Члены, линейные по константе спин-орбитальной связи, имеют вид:

(2.40)
Сопоставляя матричные элементы этого оператора с элементами матрицы
(2.41)
найдем искомое выражение для компоненты вектора . Видно, что для определения компоненты вектора необходимо вычислить матричный элемент

(2.42)

После несложного расчета получаем

(2.43)
Матричные элементы оператора чисто мнимые (пример будет приведен ниже). С учетом этого в случае одинаковых магнитных ионов после объединения подобных членов в (2.43) находим

(2.44)
Рассмотрим пример пары ионов меди. Волновые функции основных состояний в локальных системах координат обычно имеют вид:

В качестве возбужденного состояния, например на центре а, выступает

Необходимый нам матричный элемент равен

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

2.6. Биквадратичное обменное взаимодействие
Наиболее яркое спектроскопическое проявление этого взаимодействия – нарушение правила интервалов.

(2.45)
в энергетическом спектре обменно-связанных пар.

Главная причина появления биквадратичного обменного взаимодействия связана с обменной стрикцией.

Для пояснения её рассмотрим оператор энергии отдельной пары

(2.46)
Из-за появления обменного взаимодействия равновесное расстояние между ионами может изменяться. Равновесное расстояние должно соответствовать минимуму полной энергии. Производная по расстоянию между ионами должна быть равна нулю

(2.47)
Отсюда находим новое равновесное расстояние между спинами
(2.48)
Каждому собственному значению оператора соответствует свое равновесное расстояние. Подставляя в исходный гамильтониан, находим

(2.49)
Отсюда видно, что параметр биквадратичного обмена - оценивается по формуле

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

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

2.7. Двойной обмен (Double exchange )
Понятие о двойном обмене ввел Зинер в 1951 г.[1], анализируя магнитные свойства .

Рис. 2.14. Фазовая диаграмма La_1-xSr_xMnO₃. PM – парамагнетик, FM- ферромагнетик, CA- неколлинеарный антиферромагнетик, M- металл, I- изолятор, R, O, O’, O’’ – фазы, соответствующие различным пространственным группам симметрии.
Оператор двойного обмена совпадает с оператором кинетической энергии перескоков от одного магнитного иона к другому.

(2.51)

Отметим, что оператор (2.51) не может быть представлен через спиновые операторы! Но для характеристики уровней энергии пары ионов можно использовать полный спин пары.

Рис. 2.15. Схема двойного обмена в паре Mn⁴⁺-Mn³⁺. Оба иона находятся в высокоспиновом состоянии в соответствии с правилом Хунда. При перескоках электрона между e_g орбиталями направление спина не меняется, поэтому происходит кинематическое подмагничивание спинов. Спины остовов ионов подстраиваются параллельно (ферромагнитно).

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

(2.52)

Собственные значения - модуль этих недиагональных элементов. 6-j символ, фигурирующий в правой части, может быть представлен простой аналитической формулой. В итоге расчетов (предлагается сделать самостоятельно в виде задачи) зависимость энергии от суммарного спина пары определяется формулой [2.8]

(2.53)

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

Задача. Проверить справедливость формулы (2.53) на примере модели два центра - три электрона. Для расчета матричных элементов оператора использовать волновые функции в представлении вторичного квантования.
2.8.Обменное взаимодействие ионов смешанной

валентности
Рассмотрим простую модель с мигрирующим электроном с одного центра на другой

(2.54)

Здесь q – колебательная координата, описывающая смещение локального окружения при локализации “лишнего” электрона около центра a или b. Или для простоты - это смещение промежуточного (мостикового) иона. V – параметр электронно–колебательной связи. Блочный вид матрицы потенциальной энергии имеет вид:
(2.55)

Отсюда находим собственные значения

(2.56)

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

(2.57)
Расстояния между минимумами

(2.58)
Каждая из ям соответствует локализации электрона либо на центре b, либо на центре a.

Предположим, что ямы достаточно глубокие. Уравнение Шредингера имеет вид:

(2.59)

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

(2.60)

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

(2.61)
Полагая

, и, вводя новую переменную , имеем

(2.62)

и, следовательно,

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

(2.63)
Множитель называется фактором поляронного подавления. Он меньше единицы, поэтому миграция электрона затруднена. При сильной связи с деформациями мигрирующему «лишнему» электрону выгоднее локализоваться у одного из центров и тогда вместо ферромагнитных пар будут образовываться антиферромагнитно связанные пары и в равных пропорциях. Вещество из металлического ферромагнетика перейдет в состояние антиферромагнитного изолятора. Это явление наблюдается в манганитах. См. на рис. 2.14 область малых концентраций x.

2.9. Обменное взаимодействие ионов в кристаллах типа KMeF₃.
В качестве простого примера применения теории суперобменного взаимодействия рассмотрим соединение KMgF₃, легированное ионами переходной группы железа. Фрагмент структуры представлен на Рис. 2.16. Прежде чем анализировать суперобменное взаимодействие надо определить основное состояние ионов в данном соединении. Используем теорию кристаллического поля. В октаэдрическом кристаллическом поле состояния 3d – электрона расщепляются на триплет () и дублет (). Нижним по энергии является триплет. Основным состоянием иона ванадия является спиновый триплет . Обменное взаимодействие описывается оператором вида . Согласно

(2.29 ) параметр определяется суммой антиферромагнитных вкладов по - связям через промежуточный ион фтора . У ионов никеля основное состояние -. - оболочка полностью заполнена. Суперобменное взаимодействие реализуется только между электронами -подоболочки. . У ионов марганца в основном состоянии полузаполнены оби подоболочки. Кроме антиферромагнитных вкладов в суперобменное взаимодействие есть ферромагнитный вклад между ортогональными состояниями и подоболочек различных ионов. Оператор взаимодействия имеет вид . Аналогично получаем выражения для параметров суперобменного взаимодействия ионов Co²⁺ и Fe²⁺

, (2. 64)
. (2. 65)
В последнем выражении . Здесь учтено, что основное состояние иона Co²⁺ является смесью состояний .

Рис. 2. 16. Фрагмент структуры кристаллов KMeF₃ и KZnF₃: Me, где Me=V, Mn, Co, Ni. Ионы Me²⁺ замещают позиции Zn²⁺ и находятся в центре октаэдров, ионы фтора - в вершинах. Позиции ионов K⁺ не показаны.

Приближенно можно считать, что параметры , и примерно равны для всех пар ионов в KMeF₃ и оценить их по экспериментальным данным, т. е. либо то температуре Неля, либо по спектроскопическим данным. Затем эти параметры можно использовать для предсказания силы суперобменного взаимодействия для других ионов. Пример такого расчета приведен в Таблице 2.1. Там же для сравнения приведены экспериментальные данные. Видно, что качественное согласие теории суперобменного взаимодействия и эксперимента имеется.
Таблица 2.1. Значения параметров суперобменного взаимодействия для некоторых пар ионов.

Основное состояние	Вещ-во	J_эксп. ( в см^-1)	J_расчет, ( , )
	КСuF₃
	KMg F₃:Ni²⁺		62
	КCoF₃		19.3
	КFeF₃		9.3
	KMg F₃:Mn²⁺		4.5
	KMg F₃:V²⁺		4.5

2.10. Обменно-дипольные переходы
При исследовании оптических спектров поглощения кристаллов типа Al₂O₃, легированных ионами Cr³⁺ было замечено, что по мере увеличения концентрации ионов хрома наряду с линиями поглощения от одиночных ионов хрома, в спектре поглощения появляются дополнительные линии, интенсивность которых пропорциональна квадрату концентрации ионов Cr³⁺. Пропорциональность интенсивностей квадрату концентраций можно было объяснить исходя из предположения, что дополнительные линии обусловлены переходами между уровнями энергии обменно-связанных пар ионов Cr³⁺- Cr³⁺. Число таких пар, образующихся при легировании кристалла Al₂O₃ ионами Cr³⁺, должно быть как раз пропорционально квадрату концентрации. Однако, это при этом возникает вопрос, а почему интенсивность оптических переходов в обменно-связанных парах сравнима с интенсивность переходов у одиночных ионов Cr³⁺ . В самом деле, при концентрации ионов хрома вероятность образования обменно-связанных пар пропорциональна и, следовательно, интенсивности линий, принадлежащих обменно-связанным парам ионов Cr³⁺- Cr³⁺, должны быть ничтожно малыми.

Объяснение этого, казавшегося загадочным противоречия, было найдено японскими физиками /Y. Tanabe, T. Moriya, S. Sugano, Phys. Rev. Lett., Vol. 15, p. 1023 (1965)/. К этому времени появились новые экспериментальные данные по исследованию оптических спектров кристаллов KZnF₃, активированных ионами Mn²⁺ и Ni²⁺ . Электрические дипольные переходы между состояниями dⁿ- конфигурации в этих кристаллах запрещены. У одиночных ионов Mn²⁺ и Ni²⁺ реализуются лишь слабые магнитно дипольные переходы. У обменно-связанных пар запрет по четности снимается из-за обменной связи. В этом случае имеет место новый механизм электро-дипольных переходов. Его называют обменно-дипольным.

Для пояснения основной идеи рассмотрим упрощенную модель двух центров и двух электронов.

Рис. 2. 17. Схема обменно-дипольного перехода

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

, (2.66)

где в качестве оператора выступает оператор

. (2.67)

Здесь - оператор дипольного момента. Электроны могут перескакивать с центра в на центр а и наоборот . В этом случае, как и в задаче о кинетическом обменном взаимодействии

. (2.68)

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

, (2.69)
Матричные элементы эффективного дипольного момента определяются по правилу

. (2.70)

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

Рис. 2.18. Пример обменно-дипольных переходов в паре ионов Mn²⁺-Ni²⁺ в кристалле KZnF₃. /J. Ferguson, H. Guggenheim, Y. Tanabe, Phys. Rev. Lett., 14. 737 (1965)/. Символ Mn^* обозначает возбужденное состояние иона Mn²⁺. Параметр J' - параметр суперобменного взаимодействия иона Ni²⁺ с возбужденным состоянием иона Mn²⁺ через ион фтора.

Взаимодействие электронов проводимости через поле фононов.

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

Рис. 3. 1. Диаграммы процесса взаимодействия квазичастиц через поле фононов. Волнистая линия соответствует фонону с волновым вектором q.

Исходный гамильтониан задачи имеет вид:

(3.1)

Последнее слагаемое описывает взаимодействие электронов проводимости с фононами, , () – операторы уничтожения (рождения) фононов.

Запишем эрмитово сопряженное выражение оператора взаимодействия

. (3.2)
Сравнивая (3.1) и (3.2) , находим, что условие самосопряженности соблюдается, если
. (3.3)

Составим уравнение

, (3.4)

где

. (3.5)

Оператор будем искать в виде

. (3.6)

Вычислим коммутатор

. (3.7)
Используя антикоммутационные соотношения операторов носителей тока

, (3.8)

находим

(3.9)

Используя коммутационные соотношения для фононных операторов

, (3.10)

вычисляем

(3.11)

С учетом (3.9) и (3.11) уравнение (3.4) переписывается в виде

(3.12)

Выражения в квадратных скобках, очевидно, должны быть равными нулю. Отсюда находим

(3.13)

Итак, оператор имеет вид

(3.14)

Теперь можем найти вид эффективного гамильтониана по формуле

(3.15)

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

(3.16)

Данный оператор предназначен для описания основного состояния системы. Т. е., какой фонон рождается (уничтожается ), такой и должен уничтожаться (рождаться). С учетом этого имеем

(3.17)

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

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

(3.18)

Отметим, что операторная структура в (3.18) точно такая же как и у суперобменного взаимодействия.

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

Вычислить коммутатор

Затем учесть, что средние значения выражаются через функции Ферми

Вид поправки в эффективном операторе энергии

означает, что частота фононов изменяется по правилу:

где

Отметим, что

это поляризуемость (или зарядовая восприимчивость) электронов проводимости. Если зона полностью заполнена (или пуста), то разность функций Ферми равна нулю. Отметим также, что если параметры электрон-фононной связи не зависят от квазиимпульсов электронов проводимости, то перенормировка частоты фононов пропорциональна полной зарядовой восприимчивости электронов в зоне проводимости на этой частоте. Из проведенного рассмотрения ясно, что для исследований вида зарядовой восприимчивости носителей тока (а затем и диэлектрической проницаемости) можно использовать измерения частот фононных мод, например, методами неупругого рассеяния нейтронов. Это важная идея лежит в основе современных исследований диэлектрической восприимчивости в физике конденсированного состояния.

3.1 Гамильтониан модели Бардина, Купера, Шриффера .
Для описания сверхпроводимости выделяются аномальные средние , т.е. из оператора полного взаимодействия выделяются лишь члены с и . Такой оператор лежит в основе теории Бардина, Купера, Шриффера (БКШ)

(3.19)
Величины и получаются из сравнения (3.19) и (3.18)

(3.20)

(3.21)

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

(3.21)

3.2 Уравнение для параметра порядка при переходе в состояние с волнами зарядовой плотности.
Предположим, что имеется модуляция в плотности распределения электронов проводимости с волновым вектором . В представлении вторичного квантования оператор плотности заряда имеет вид:
(3.22)
Отсюда видно, средние значения плотности будут отличны от нуля, если . Эти средние зависят от температуры и рассчитываются самосогласованно на основе оператора взаимодействия через поле фононов

(3.23)
Для компонент и получаем следующие формулы:

(3.24)
(3.25)
Частота фононов в кристалле ограничена некой предельной частотой (частота Дебая). Это обстоятельство учтено путем введения тета-функции . Кроме того учтено, что во взаимодействии могут участвовать лишь состояния вблизи . Последнее обстоятельство и привело к еще одной тета – функции .

Модели перехода металл- диэлектрик.

4.1. Модель Пайерлса.
Гамильтониан имеет вид:

(4.1)

Энергия отсчитывается от уровня Ферми. Спиновое квантовое число для краткости записи опущено (подразумевается). Составляем уравнения движения

(4.2)

Здесь предполагается, что

(4.3)

Решение уравнений (4.2) ищем в виде.

(4.4)

Для краткости записи полагаем . После подстановки (4.4) в (4.2) имеем

(4.5)
Система линейных уравнений (5.5) имеет нетривиальное решение, если её определитель равен нулю
(4.6)

Корни уравнения равны

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

Рис.4.1. Два возможных варианта выбора зоны. а) – первая зона Бриллюэна. Область изменения волнового вектора уменьшена в два раза по сравнению с областью в исходной зоне. б) - первая и вторая зоны Бриллюэна объединены. Удобство такого выбора заключается в том, что случаи и можно обсчитывать в рамках одной программы, не меняя области изменения волновых векторов.

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

Общее решение системы уравнений - это сумма частных решений

(4.8)

Они определяются из начальных условий. Положив в (4.8), имеем

(4.9)

Еще одна пара уравнений получается из уравнений для производных по времени при

(4.10)
Фигурирующие в правой части операторы (4.10) соответствуют времени . Из (4.10) с учетом (4.9) находим

(4.11)
Уравнения (4.10) относятся к классу однородных уравнений. Они решаются с точностью до постоянного (нормировочного множителя). По виду временной зависимости операторов в (4.8) видно, что

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

Найдем нормировочные коэффициенты и так, чтобы

(4.13)

Вычисления дают

(4.14)

Итак, имеем

(4.15)

Здесь удобно ввести обозначения

(4.16)

Заметив, что

(4.17)
формулы (4.15) удобно переписать в более компактном виде
(4.18)

Формула (4.18) аналогична той, что используется в теории сверхпроводимости (преобразования Боголюбова).

Формулы для обратного преобразования имеют вид:

(4.19)

Для дальнейшего полезно отметить следующие соотношения:

(4.20)

Для контроля и лучшего понимания результата полезно осуществить проверку справедливости формул (4.18) в произвольный момент времени. Непосредственной подстановкой (4.19) в (4.1) убедимся, что в представлении новых операторов гамильтониан (4.1) приобретает вид:

(4.21)

Перепишем оператор (4.1) в виде

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

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

4.2. Переход Мотта

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

(4.24)

где параметр экранирования примерно пропорционален концентрации носителей

. (4.25)

Если он будет достаточно мал, то кулоновское поле иона способно обратно захватить электрон и локализовать его на первой боровской орбите с радиусом

. (4.26)

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

. (4.27)

Это условие удобно переписать в виде

(4.28)

Заметим, что при увеличении температуры число свободных электронов увеличивается. Отсюда заключаем, что при низких температурах наиболее логично ожидать диэлектрическое состояние, а при высоких – металлическое. Этот сценарий перехода диэлектрик - металл был предложен Моттом. Создание количественной теории перехода диэлектрик - металл составляет одно их актуальных направлений современной физики конденсированного состояния.

4.3. Модель Хаббарда.
Рассмотрим упрощенную модель вещества с гамильтонианом вида:
(4.29)

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

Для частичной диагонализации гамильтониана (6) удобно ввести операторы проектирования (операторы Хаббарда)

. (4.31)
Удобство их использования связано с тем, что квадратичные члены, соответствующие одному узлу, сводятся к линейным:
(4.32)
Любой оператор может быть выражен через операторы Хаббарда по формуле

(4.33)

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

(4.34)

Можно, конечно, и наоборот, выразить операторы Хаббарда через обычные фермионные операторы. Так, например,

(4.35)

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

(4.36)

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

Найдем также

Учитывая все это, находим вид гамильтониана в представлении операторов:

(4.37)

Составим теперь уравнение движения для операторов рождения. Для пояснения основной идеи Хаббарда считаем, что . Тогда

(4.38)

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

Рис. 4.2. Схема расщепления хаббардовских подзон.

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

Переход диэлектрик – металл произойдет, когда “размытие” подзон из-за перескоков электрона (параметр t_ij) окажется сравнимым с величиной U, т.е. когда графики плотности состояний верхних и нижних подзон перекроются. Или иными словами диэлектрическая щель закроется. Подчеркнем , что в отличие от сценария Пайерлса переход происходит без изменения периода решетки. Сценарии Мотта и Хаббарда взаимно дополняют друг друга. В этой связи соединения с описанным механизмом перехода металл- диэлектрик выделяют в отдельный класс мотт-хаббардовских диэлектриков.

Картина перестройки плотности состояний в зоне по мере увеличения параметра U иллюстрируется на рис. 4.3

Рис.4.3 Эволюция плотности состояний в модели Хаббарда при изменении параметра перескока (- ширина зоны при ), рассчитанная в работе: Yu. Irkhin, A. V. Zarubin, Eur. Phys. J., B, 38, 563 (2004). Имеется в открытом доступе по адресу: arXive: cond-mat/0103054v3. Нижний рисунок соответствует металлическому состоянию с одним электроном в элементарной ячейке, верхний – диэлектрическому.
Для дальнейшего знакомства с теорией перехода металл-диэлектрик нужны знания по методу функций Грина.
4.4 Метод функций Грина.
Функции Грина введены в работах С.В. Тябликова и Н. Н. Боголюбова.

Запаздывающие функции Грина определяются следующим образом

(4.39)

Здесь - тета-функция. При t>0 она равна 1, а при t

(4.40)
Значение или -1 выбирается в зависимости от решаемой задачи. Обычно для бозонных операторов и для фермионных.

Составляя уравнения движения по обычному правилу

, (4.41)
можно получить уравнения для функций Грина. Решения уравнений ищут в виде

(4.42)

Фурье-образы функций Грина оказываются связанными цепочкой уравнений:

(4.43)

Если найдена функция Грина, то можно найти корреляционную функцию по формуле (спектральная теорема)

(4.44)

В правой части (20) берется предел . Полезно посмотреть, как этот формализм работает для электронов без учета их взаимодействия, т. е. на задаче, ответ которой мы уже знаем. Итак, пусть

. (4.45)

Переходим к новым переменным

. (4.46)

Оператор кинетической энергии принимает вид:

, (4.47)

здесь

(4.48)

имеет смысл энергии квазичастиц.

Составим уравнение для фурье-образа функции Грина . Используя (4.43), имеем

(4.49)

Отсюда находим

(4.50)
Теперь по формуле (4.40) находим корреляционную функцию, которая совпадает с функцией распределения Ферми.

(4.51)

Вернемся теперь к модели Хаббарда. Оператор перескоков в нижней подзоне имеет вид:

(4.52)

Вычислим коммутатор, используя правило

(4.53)

Здесь в качестве антикоммутаторов берутся их средние значения (это, так называемое, приближение - Хаббард I , использованное им в его первой статье на эту тему)

, (4.54)

а среднее значение .

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

, (4.55)

где

(4.56)

Вместо (4.48) имеем

(4.57)

Видно, что ширина зоны уменьшена (редуцирована) на множитель . Другой новый эффект в модели Хаббарда – необычное поведение спектрального веса или иначе сказать “емкости зоны”. Более точно, спектральный вес - это отношение числа носителей к максимально возможному их числу в зоне. В нижней хаббардовской подзоне число электронов в расчете на один узел не может быть больше единицы.

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

Условие полноты имеет вид

(4.58)

Оператор концентрации носителей тока определяется выражением

(4.59)

Здесь соответствует заполнению верхней подзоны, а появлению дырок в нижней. В парамагнитной фазе . С учетом этого из (4.33) и (4.34) находим

(4.60)

Теперь формула, аналогичная (4.51), запишется в виде:

(4.61)

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

(4.62)

В случае полностью заполненной зоны все ферми-функции, фигурирующие справа в (4.62), равны единице. При половинном заполнении левая часть (4.62) равна 1/2, что соответствует индексу заполнения (дырочному) .

Рис. 4.4. Графики спектральных весов. Левая сторона ( ) соответствует нижней подзоне Хаббарда. Правая ( ) – верхней. Половинное заполнение зон достигается при концентрации носителей .

Сверхпроводники
5.1. Суперобменный механизм спаривания носителей тока

Суперобменное взаимодействие связано с перескоками электронов. Эффективный оператор имеет вид:

( 5.1)
Здесь интеграл перескока электрона (или дырки) между узлами кристаллической решетки. Переходим к новым операторам

( 5.2)

( 5.3)

Подставляем их в (5.1) и анализируем суммы по номерам узлов решетки. В решетке с трансляционной симметрией сумма

(5.4)
не зависит от номера узла i. Это обстоятельство позволяет вычислить сумму

(5.5)

В итоге оператор суперобмена приобретает вид

(5.6)

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

Далее в (5.6) проводим частичное усреднение следующим образом
(5.7)

Из всех возможных средних оставляем лишь такие, которые соответствуют куперовским парам с нулевым суммарным импульсом

(5.8)

Полученный таким образом эффективный оператор имеет вид

,

(5.9)

где

Первое слагаемое учитывает кинетическую энергию электронов проводимости. Диагонализация модельного гамильтониана (5.9)

Бардина-Купера-Шриффера (БКШ) может быть проведена различными методами. Для учебных целей полезно использовать самый простой.

Составим уравнения движения для операторов рождения и уничтожения электронов (дырок)

(5.10)

Решение ищем в виде

(5.11)

Подставляем (5.10) в (5.9)

(5.12)

Система имеет решение, если её определитель равен нулю. Отсюда находим

(5.13)

Общее решение системы (5.9) это сумма частных решений

(5.14)

Здесь . Для нахождения полной системы для операторов C подставляем (5.14) в (5.9) и для простоты расчетов берем

(5.15)

Из первого и третьего уравнений находим

(5.16)
Из второго и четвертого уравнений получаем
(5.17)

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

(5.18)

Множитель подчиним условию

(5.19)

Составляем антикоммутатор

(5.20)

Следовательно, оператор уничтожения определяется выражением

(5.21)

Оператор рождения найдем через C₄.

Вычисляем коммутатор

(5.22)

Отсюда находим

(5.23)

Отметим, что оператор (5.9) при использовании новых квазичастичных операторов записывается в виде:

(5.24)

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

Далее уточним уравнение (5.8) для самосогласованного расчета параметра энергетической щели. Исходя из (5.23) и (5.21) можно записать

(5.25)

Здесь удобно ввести обозначения

. (5.26)

Полезно отметить, что

(5.27)

Используя (5.25), находим

(5.28)

Подставляя (5.28) в (5.8), получаем уравнение

(5.29)

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

Видно, что зависимость параметра сверхпроводящей щели от волнового вектора имеет вид:

(5.30)
Уравнение (5.29) распадается на систему четырех уравнений
(5.31)
(5.32)

Из (5.29) следует, что . В таком случае из (5.32) следует, что . Решение (5.30) должно соответствовать какому-то одному неприводимому представлению точечной группы симметрии зоны Бриллюэна. Из (5.30) видно, что возможны два варианта решений

(5.33)

Первый из них инвариантен относительно всех операций симметрии в зоне Бриллюэна. Говорят, что параметр щели имеет симметрию s- типа. Параметр меняет знак при повороте на 90 градусов. Он ведет себя как электронное состояние типа . В этой связи говорят, что имеет симметрию d- типа. Численные решения системы уравнений (5.31) показывают, что для медь-кислородных слоев в соединениях типа YBa₂Cu₃O₇ (T_c=92 K) реализуется решение d- типа со значением . Соответствие расчета экспериментальным данным как по зависимости от волнового вектора, так и по отношению свидетельствует в пользу доминирующей роли суперобменного механизма спаривания носителей тока в сверхпроводящих купратах. Исследования в этом направлении продолжаются, так как пока не ясно, какую роль при этом играют взаимодействие через оптические фононы и экранированное кулоновское взаимодействие между носителями тока.

Задача. При рассмотрении куперовских пар с ненулевым суммарным импульсом в операторе (5.7) следует сохранить аномальные средние вида

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

.

0>следующая страница >>