Электронная структура и магнитные свойства квазидвумерных систем MeO (Me = Mn, Fe, Co, Ni) в ферромагнитном упорядочении (стр. 1 )

Электронная структура и магнитные свойства квазидвумерных систем MeO (Me = Mn, Fe, Co, Ni) в ферромагнитном упорядочении

И.Г.Попова, И.В. Ершов,

Донской государственный технический университет, Ростов-на-Дону

Аннотация: Приводятся расчёты электронной структуры квазидвумерных систем 2D MeO (Me = Mn, Fe, Co, Ni) выполненные с использованием метода псевдопотенциала в рамках теории функционала плотности (DFT). Расчёты электронной плотности выполнены в приближении коллинеарного магнетизма, исследована электронная структура, выполнен прогноз электронных и магнитных свойств в квазидвумерных системах 2D MeO (Me = Mn, Fe, Co, Ni) с ферромагнитным типом упорядочения. Выполненные оценки магнитного момента на атоме переходного 3d-металла в кристаллах оксидов согласуются с экспериментом. Обнаружен эффект спонтанной спиновой поляризации Me 3d и 2p-зон атомов кислорода в квазидвумерных системах 2D MeO (Me = Mn, Fe, Co, Ni).

Ключевые слова: Теория функционала плотности, метод псевдопотенциала, ферромагнетизм, электронная структура, зонная структура, оксиды переходных металлов.

Введение

В последние годы повышенное внимание к себе стали привлекать пленочные структуры оксидов переходных металлов, классифицируемые как квазидвумерные структуры. Успехи в выращивании ультратонких слоёв MnO, CoO, NiO и др. методом импульсного лазерного осаждения [1] позволили изготавливать качественные магнитные гетероструктуры [2]. Изучение обменного взаимодействия в многослойных магнитных системах представляет повышенный интерес, что обусловлено, в частности, возможностью управления спином. Особенности роста и атомная структура ультратонких слоёв, например 2D MnO, на поверхности металлов и полупроводников, полученных методом импульсного лазерного осаждения, существенно зависят от толщины (1, 3, 7, 67 атомных слоёв [3]) выращенного слоя. Методами абсорбционной (XAS) и фотоэлектронной спектроскопии (XPS) показано [4], что электронные и кристаллические свойства в ультратонкой плёнке 2D MnO/Ag(001) и кристалле 3D MnO различаются. Экспериментальные рентгеновские спектры XAS и XPS могут также содержать небольшое количество артефактов и поэтому необходим их детальный теоретический анализ. Электронная структура и магнитные свойства 2D и 3D оксидов марганца в ферромагнитном упорядочении изучены достаточно полно с использованием зонных методов [5,6], результаты расчётов которых находятся в хорошем согласии с экспериментальными данными. Зонная структура кристаллов 3D MnO в ферромагнитном (FM) упорядочении соответствует металлическому типу проводимости. Изучение зонной структуры оксидов переходных металлов типа 2D MeO (Me = Mn, Fe, Co, Ni) в ферромагнитном упорядочении целесообразно проводить в сопоставлении с 3D MeO (Me = Mn, Fe, Co, Ni), что позволяют отметить их некоторые особенности. В частности, расчёты зонной структуры 2D систем MeO (Me = Fe, Co, Ni) в ферромагнитном упорядочении, выполненные методом псевдопотенциала с использованием теории функционала плотности (DFT) [7], показали на возрастание Me3d-O2p-гибридизации заполненных состояний в 2D системах относительно кристалла. Следует ожидать, что электронные и магнитные свойства квазидвумерных систем FeO, CoO и NiO в ферромагнитном упорядочении будут иметь отличие от свойств объёмных кристаллов одноименных оксидов. Таким образом, неполнота данных о природе свойств в магнитных ультратонких слоях типа 2D MeO (Me = Mn, Fe, Co, Ni), об особенностях спинового состояния в валентной зоне, на уровне Ферми и межатомного взаимодействия в двумерных структурах побудили нас к теоретическому изучению особенностей электронной структуры и свойств в ультратонких слоях 2D MeO (Me = Mn, Fe, Co, Ni).

Метод расчета

Расчёты электронной структуры квазидвумерных систем 2D MeO (Me = Mn, Fe, Co, Ni) выполнены с использованием программного пакета Quantum Espresso (QE) [8]. В основе (QE) которого лежит метод псевдопотенциала с использованием плоских волн, в рамках теории функционала плотности (DFT) Кона-Шэма. Вид локального эффективного потенциала приведен, например, в работе [9]. Энергия обрезания плоских волн для самосогласованного расчета составляла 30 Рид, что обеспечивало сходимость по полной энергии ячейки не хуже Рид/яч. Число -точек в обратном пространстве для самосогласованного расчёта выбиралось из таких же соображений и составило 28 -точек в неприводимой части зоны Бриллюэна (ЗБ), при этом набор плоских волн составил величину порядка 50000. Влияние остовных электронов учитывалось путём использования ультрамягких псевдопотенциалов, рассчитанных из первых принципов. Для корректного воспроизведения спинового состояния магнитных материалов нами использован нелокальный обменно-корреляционный функционал в параметризации Perdew-Burke-Ernzerhof (PBE) аналогично работе [10]. Расчёты электронной плотности выполнены в спин-поляризованном варианте. Проводилось раздельное интегрирование электронной плотности для двух проекций спина и вычисление их разности позволило найти интегральный спиновой магнитный момент на каждом атоме. Поверхность квазидвумерных систем 2D MeO (Me = Mn, Fe, Co, Ni) моделировали как полосу, состоящую из 3-х слоёв оксида металла и включающую 27 атомов в суперячейке, и каждая полоса была отделена от другой областью вакуума величиной около 9 Ǻ. При построении полосы мы использовали те же кристаллические структуры, что и в случае 3D кристаллов (типа NaCl), но с некоторыми изменениями: во-первых мы ограничили трансляционную симметрию кристаллической структуры по одной из осей, а именно по оси c; во-вторых, при рассмотрении элементарной ячейки каждой из систем придется вводить большее число базисных атомов, что связано с понижением симметрии кристаллической структуры с гранецентрированной кубической до тетрагональной. Зона Бриллюэна при переходе 3D MnO →2D MnO изменяется от неправильного додекаэдра к прямоугольному параллелепипеду [6]. Вводя такие изменения, можно рассматривать именно те атомы, которые непосредственно находятся на границе интерфейса. Ограничение трансляционной симметрии в нашем случае будет достигнуто путем увеличения параметра трансляции элементарной ячейки по оси c в 2 раза. Результирующие значения параметров модифицированных квазидвумерных кристаллических структур сведены в табл. 1.

Таблица №1.

Параметры кристаллических решеток исследуемых соединений в 2D-структурах.

Соединение	Параметр решетки, Å	Объем примитивной ячейки, Å3
а	с
MnO	4.44	8.88	175.05
FeO	4.43	8.86	173.88
CoO	4.26	8.52	154.62
NiO	4.18	8.36	146.07

Результаты и обсуждение

Электронная структура квазидвумерных оксидов 2D MeO (Me = Mn, Fe, Co, Ni) рассчитана с использованием теории функционала плотности и представлена на рис. 1, в сопоставлении с электронной структурой кристаллов 3D MeO (Me = Mn, Fe, Co, Ni). На рис. 1а, б, в, г представлены парциальные плотности электронных состояний квазидвумерных систем 2D MeO (Me = Mn, Fe, Co, Ni) для электронных подсистем со спином вверх и спином вниз. Картина распределения электронных плотностей состояний по энергии для 2D MnO (рис. 1а) принципиально не отличается от 3D MnO (рис. 1д). Немного изменяется относительное положение уровня Ферми, который сдвигается в область более высоких энергий. Анализ распределения парциальных плотностей состояний по энергии для 2D-оксидов, представленные для атомов интерфейса (верхнего слоя), позволяет отметить, что аналогично 3D-системам в ряду MnO - FeO - CoO - NiO наблюдается смещение 3d-состояний металла со спином вниз N↓ относительно уровня Ферми в область низких энергий. Отмеченное определенным образом влияет на намагниченность (выравнивание магнитных моментов) в двумерных системах оксидов переходных 3d-металлов.

Рис.1. - Парциальные плотности электронных состояний для 2D-оксидов: MnO (а), FeO (б), CoO (в),NiO (г) и 3D-оксидов: MnO (д), FeO (е), CoO (ж), NiO (з) в ферромагнитном упорядочении.

В частности, карты распределения электронной плотности в квазидвумерных системах 2D MeO (Me = Mn, Fe, Co, Ni), представленные на рис. 2, иллюстрируют характер распределения плотности заряда между атомами в системе. В центре каждой карты расположен атом кислорода, а ближайшие к нему – атомы соответствующего металла. Из приведенных карт видно, что преобладающим для всех соединений является ионный тип связи, о чем свидетельствует высокая концентрация электронной плотности вокруг атомов. По сравнению с 3D системами наблюдается некоторое нарушение симметрии в распределении плотности, что связано с перераспределением последней в связи с наличием поверхности интерфейса.

Рис.2. - Карты распределения электронной плотности для 2D-оксидов: MnO (а), FeO (б), CoO (в),NiO (г) в ферромагнитном упорядочении

Отмеченные выше эффекты могут быть связаны с заполнением d-оболочки металла при увеличении атомного номера. Следует также отметить расщепление 3d-состояний со спином вниз N↓ на две и даже на три полосы, как например, для CoO и NiO соответственно. Данные о парциальных зарядах валентных электронов на атомах 2D-оксидов и интегральный спиновой магнитный момент (ММ) на каждом атоме представлены в табл. 2. Анализ представленных в таблице данных позволяет отметить, с одной стороны, возрастание парциального заряда 3d-электронов металла со спином вниз N↓ в ряду MnO - FeO - CoO – NiO, а с другой – уменьшение магнитного момента на атомах металла. Следует отметить достаточно заметный магнитный момент на атомах кислорода в оксидах железа и никеля. Полученные в настоящей работе данные о парциальных зарядах валентных электронов на атомах в объёмных кристаллах оксидов переходных 3d-металлов и интегральный спиновой магнитный момент на каждом атоме представлены в табл. 3.