Scientific journal
International Journal of Applied and fundamental research
ISSN 1996-3955
ИФ РИНЦ = 0,593

SPIRAL MAGNETIC ORDERING AND PHASE SEPARATION IN METALLIC SYSTEMS WITH THE HYBRIDIZATION OF THE ELECTRONIC SPECTRUM

Igoshev P.A. 1 Igoshev P.A. 2 Antipin T.V. 1 Antipin T.V. 2 Irkhin V.Y. 1
1 Institute of Metal Physics
2 Ural Federal University
1182 KB
The ground state magnetic phase diagram of the Anderson lattice model is calculated within a generalized Hartree-Fock approximation with account of spiral magnetic ordering. The consistent treatment of phase transitions between different magnetic phases allows to describe the magnetic phase separation which is a general consequence of density and magnetization jump for first-order transitions. The physical picture of spiral magnetism in the cases of integer and intermediate valence appears to be substantially different. In the former case the perturbation theory in the inverse Coulomb interaction and d-level energy yields the effective s-d exchange model. The latter gives the picture of magnetic order, which is qualitative close to that obtained earlier within the Hubbard model.
spiral magnetic order
itinerant magnetism
phase separation
Anderson lattice model

В последние годы классическая проблема связи магнитных свойств сильно коррелированных соединений переходных металлов с геометрией решетки и электронной структурой продолжает привлекать внимание исследователей. Чтобы получить физическую картину магнитного порядка, необходимо учитывать, наряду с влиянием электронных корреляций и спиральным (неколлинеарным) магнитным порядком, магнитное фазовое расслоение, возникающее при переходах первого рода (первое свидетельство о существовании этого явлении было получено в работе [9]), и влияние многоорбитальных эффектов. Кроме того, спиральный магнитный порядок может играть важную роль в механизме перехода металл-изолятор [1]. Для исследования формирования магнитного порядка используются различные многоэлектронные модели, такие как модель Хаббарда, s-d обменная модель и периодическая модель Андерсона. Эти модели широко применяются для описания свойств соединений переходных и редкоземельных металлов, высокотемпературных сверхпроводящих купратов, систем с тяжелыми фермионами [2,5]. Формирование магнитных спиралей было исследовано в рамках невырожденной модели Хаббарда [1,7], однако влияние гибридизации электронных состояний на конфигурацию магнитного порядка требует дополнительных исследований. Целью настоящей работы является изучение магнитной фазовой диаграммы периодической модели Андерсона для квадратной решетки с учетом магнитного фазового расслоения, а также неколлинеарных магнитных состояний и их взаимосвязи. Поставленная задача позволит пролить свет на роль многоорбитальных эффектов в формировании спирального магнитного порядка.

Формализм

Гамильтониан системы имеет вид

igosh1.wmf (1)

где igosh2.wmf – операторы рождения (уничтожения) s-электрона в зоне проводимости и локализованного d-электрона на узле i, обладающих проекцией спина igosh3.wmf, соответственно, igosh4.wmf – оператор числа локализованных электронов на узле i, обладающих проекцией спина igosh5.wmf, igosh6.wmf – интеграл переноса между узлами i и j, igosh7.wmf – энергия локализованного электронного состояния, V – sd-гибридизация на одном узле, последний член отвечает за кулоновское взаимодействие d-электронов на одном узле. Полная концентрация электронов в системе

igosh8.wmf,

где igosh9.wmf – числа заполнения для зонного и локализованного состояний соответственно (⟨...⟩ – статистическое среднее).

Предположим, что имеет место спиральный магнитный порядок в плоскости xy, модулированный волновым вектором Q

igosh11.wmf (2)

где igosh12.wmf и igosh13.wmf – единичные векторы координатных осей, igosh14.wmf – оператор намагниченности на узле, igosh15.wmf – вектор матриц Паули. Выполняя локальный поворот в спиновом пространстве на узле i, что позволяет совместить направление намагниченности на разных узлах,

igosh16.wmf (3)

и используя стандартное приближение Хартри-Фока для представления кулоновского взаимодействия в новой системе координат

igosh17.wmf (4)

мы получим гамильтониан спирально упорядоченной системы в обобщенном приближении Хартри-Фока.

igosh18.wmf

igosh19.wmf (5)

где igosh20.wmf igosh21.wmf igosh22.wmf

При решении уравнений обобщенного приближения Хартри-Фока выполнялась численная диагонализация гамильтониана (5) на сетке 100×100 в зоне Бриллюэна. Для определения энергетически выгодной спиральной фазы (в число которых включаются ферромагнетизм с Q=0 и антиферромагнетизм с igosh23.wmf) вычислялся и минимизировался термодинамический потенциал igosh24.wmf для каждого волнового вектора Q спирали, выбранного на достаточно большой сетке в зоне Бриллюэна 20×20 [1,6], где µ – химпотенциал. Варьируя химпотенциал, можно определить зависимость магнитной конфигурации кристалла от концентрации носителей, в том числе и с учетом магнитных переходов первого рода, возникающих при условии igosh25.wmf. В этом случае при igosh26.wmf параметры магнитного упорядочения испытывают скачок от igosh27.wmfи система при igosh28.wmf неустойчива как однородная (явление «магнитного» фазового расслоения).

Если d-уровень лежит значительно ниже энергии Ферми и igosh29.wmf гамильтониан (1) может быть сведен с помощью преобразования Шриффера-Вольфа [8] у s-d модели с параметром обмена

igosh30.wmf (6)

где igosh31.wmf – уровень Ферми.

Результаты. Результаты вычислений, описанных в предыдущем разделе, представлены на рис. 1–4. На фазовых диаграммах спиральные фазы идентифицированы видом их волнового вектора. Закрашенная область ‘’Ф1+Ф2» означает фазовое расслоение между фазами Ф1 и Ф2. Толстые линии означают переходы второго рода, тонкие линии – границы области фазового расслоения.

На рис. 1 представлен пример магнитной фазовой диаграммы основного состояния в переменных igosh33.wmf для случая igosh34.wmf, когда имеет место частично-дырочная симметрия.

igos1.tiff

Рис. 1. Магнитная фазовая диаграмма основного состояния решеточной модели Андерсона в переменных n–U для случая igosh36.wmf «AFM» – антиферромагнитный порядок с igosh37.wmf «FM» – ферромагнетизм

При igosh38.wmf когда перенос тока осуществляется d-состояниями и наша модель каноническим преобразованием сводится к эффективной модели Хаббарда, для достаточно больших U доминируют ферро- и антиферромагнитное упорядочение (между ними имеется фазовое расслоение), при меньших U – спиральный магнитный порядок. При igosh39.wmf для всего диапазона значений U, за исключением достаточно малых, доминирует спиральный магнитный порядок. Вблизи половинного заполнения igosh40.wmf имеется антиферромагнитная область.

На рис. 2 представлен пример фазовой диаграммы основного состояния в переменных igosh41.wmf для случая умеренных значений igosh42.wmf и U, что обеспечивает нецелое заполнение локализованного состояния (см. рис. 4, случай промежуточной валентности). Видно, что распределение фаз на фазовой диаграмме для случаев igosh43.wmf существенно отличаются.

igos2.tiff

Рис. 2. Фазовая диаграмма основного состояния решеточной модели Андерсона с U=12tigosh45.wmf Слева: igosh46.wmf справа: igosh47.wmf Обозначения фаз те же, что и на рис. 1

При n=1 и n=2 имеет место шахматный антиферромагнитный порядок, при легировании в области igosh49.wmf возникает область магнитного фазового расслоения: при малых igosh50.wmf – с участием «параллельной» спиральной фазы, при большем V и igosh51.wmf – с участием ферромагнитной фазы, причем ее ширина в обоих случаях существенно превышает соответствующую область для невырожденной модели Хаббарда [1,6,7]. С другой стороны, при igosh52.wmf и значительном легировании фазовое расслоение практически исчезает. Таким образом, картина магнитного упорядочения является существенно ассиметричной по отношению к знаку носителей заряда.

Совершенно иная картина наблюдается для случая целой валентности, что реализуется для igosh53.wmf Пример такой фазовой диаграммы представлен на рис. 3. Видно, что фазовая диаграмма практически симметрична по отношению к знаку носителей заряда, а распределение фаз подобно случаю промежуточной валентности при igosh54.wmf. Физически реализуется случай решетки локализованных магнитных моментов, взаимодействующих друг с другом через электроны проводимости (подобно взаимодействию РККИ). Ясно, что в этом предельном случае реализуется спиральная фаза с волновым вектором, обеспечивающим максимальность функции Линдхарда

igosh55.wmf (7)

где igosh56.wmf – функция Ферми. Выражение для эффективного обмена дается формулой (6).

igos3.tif

Рис. 3. Фазовая диаграмма основного состояния решеточной модели Андерсона с U=200t и igosh57.wmf. Слева: в терминах V, справа: в терминах эффективного s-d обменного параметра I

Наши результаты находятся в согласии с результатами s-d модели [10] при малом обменном параметре I (отметим однако, что спиральная фаза (Q,Q) полученная в это работе, оказывается нестабильной по отношению к фазовому расслоению между неелевской антиферромагнитной фазой и спиральной фазой igosh58.wmf). Кроме того, как можно видеть, последовательность магнитных фаз на рис. 3 качественно согласуются с результатами для невырожденной модели Хаббарда [6] при промежуточных U (концентрация носителей равна igosh59.wmf для модели Хаббарда и igosh60.wmf для модели Андерсона). Однако, в отличие от модели Хаббарда, где при малых U реализуется парамагнитная фаза, в модели Андерсона (или s-d обменной модели) магнитные моменты не исчезают с уменьшением s-d обменного взаимодействия и парамагнитная фаза никогда не реализуется.

Чтобы проиллюстрировать поведение чисел заполнения для зонных и локализованных электронных состояний, на рис. 4 представлена их зависимость от концентрации электронов. Заметим, что направление намагниченности для s- и d-состояний противоположны. В случае igosh61.wmf соответствующем целой валентности igosh62.wmf фазовая диаграмма почти симметрична. С другой стороны, когда igosh63.wmf мы имеем случай промежуточной валентности, а для igosh64.wmf – случай почти целой валентности. В этом случае фазовая диаграмма становится существенно ассиметричной. Для различных спиральных состояний степень насыщенности магнетизма (отношение намагниченности к числу носителей практически не зависит от концентрации как для s-, так и для d-состояний).

igos4.tif

Рис. 4. Зависимость чисел заполнения (сплошная линия) и намагниченности на одном узле m (пунктирная линия) от концентрации электронов. Разрыв линий происходит при фазовом переходе первого рода, что связано со скачком волнового вектора Q спирального состояния. Тонкие линии соответствуют областям фазового расслоения. Обозначены типы спирального упорядочения. Слева: igosh65.wmf Справа: igosh66.wmf

Заключение

Мы исследовали формирование спирального магнитного порядка для периодической модели Андерсона в различных режимах, включая случай целой и промежуточной валентности, и установили связь с другими многоэлектронными моделями. Мы показали, что картина формирования спиральных состояний и магнитного фазового расслоения, важного для кинетических явлений [3,4], существенно меняется по сравнению со случаем невырожденной модели Хаббарда. Однако и в том, и в другом случае реализуется богатое многообразие спиральных фаз. Эти результаты могут быть применены для объяснения электронных свойств соединений переходных металлов.

Работа поддержана проектом УрО РАН №15-8-2-9.