CC BY
58
Электронная структура и магнитные свойства квазидвумерных систем MeO (Me = Mn, Fe, Co, Ni) в ферромагнитном упорядочении
И.Г.Попова, И.В. Ершов, В.В. Илясов Донской государственный технический университет, Ростов-на-Дону
Аннотация: Приводятся расчёты электронной структуры квазидвумерных систем 2D MeO (Me = Mn, Fe, Co, Ni) выполненные с использованием метода псевдопотенциала в рамках теории функционала плотности (DFT). Расчёты электронной плотности выполнены в приближении коллинеарного магнетизма, исследована электронная структура, выполнен прогноз электронных и магнитных свойств в квазидвумерных системах 2D MeO (Me = Mn, Fe, Co, Ni) с ферромагнитным типом упорядочения. Выполненные оценки магнитного момента на атоме переходного 3d-металла в кристаллах оксидов согласуются с экспериментом. Обнаружен эффект спонтанной спиновой поляризации Me 3d и 2р-зон атомов кислорода в квазидвумерных системах 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-02р-гибридизации заполненных состояний в 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) которого лежит метод псевдопотенциала с
использованием плоских волн, в рамках теории функционала плотности (ББТ) Кона-Шэма. Вид локального эффективного потенциала приведен, например, в работе [9]. Энергия обрезания плоских волн для самосогласованного расчета составляла 30 Рид, что обеспечивало сходимость по полной энергии ячейки не хуже 10-4 Рид/яч. Число к -точек в обратном пространстве для самосогласованного расчёта выбиралось из таких же соображений и составило 28 к -точек в неприводимой части зоны Бриллюэна (ЗБ), при этом набор плоских волн составил величину порядка 50000. Влияние остовных электронов учитывалось путём использования ультрамягких псевдопотенциалов, рассчитанных из первых принципов. Для корректного воспроизведения спинового состояния магнитных материалов нами использован нелокальный обменно-корреляционный функционал в параметризации Ре^е,№-Вигке-Ег^еАоГ (РВЕ) аналогично работе [10]. Расчёты электронной плотности р(г) выполнены в спин-поляризованном варианте. Проводилось раздельное интегрирование электронной плотности р(г) для двух проекций спина и вычисление их разности позволило найти интегральный спиновой магнитный момент на каждом атоме. Поверхность квазидвумерных систем 2Б МеО (Ме = Мп, Бе, Со, N1) моделировали как полосу, состоящую из 3-х слоёв оксида металла и включающую 27 атомов в суперячейке, и каждая полоса была отделена от другой областью вакуума величиной около 9 А. При построении полосы мы использовали те же кристаллические структуры, что и в случае 3Б кристаллов (типа №С1), но с некоторыми изменениями: во-первых мы ограничили трансляционную симметрию кристаллической структуры по одной из осей, а именно по оси с; во-вторых, при рассмотрении элементарной ячейки каждой из систем придется вводить большее число базисных атомов, что связано с понижением симметрии кристаллической структуры с гранецентрированной кубической до тетрагональной. Зона Бриллюэна при переходе 3Б МпО МпО
изменяется от неправильного додекаэдра к прямоугольному параллелепипеду [6]. Вводя такие изменения, можно рассматривать именно те атомы, которые непосредственно находятся на границе интерфейса. Ограничение трансляционной симметрии в нашем случае будет достигнуто путем увеличения параметра трансляции элементарной ячейки по оси с в 2 раза. Результирующие значения параметров модифицированных квазидвумерных кристаллических структур сведены в табл. 1.
Таблица №1.
Параметры кристаллических решеток исследуемых соединений в 2Б-структурах.
Соединение Параметр решетки, А Объем
примитивной ячейки, А3
а с
МпО 4.44 8.88 175.05
БеО 4.43 8.86 173.88
СоО 4.26 8.52 154.62
№О 4.18 8.36 146.07
Результаты и обсуждение
Электронная структура квазидвумерных оксидов 2Б МеО (Ме = Мп, Бе, Со, N1) рассчитана с использованием теории функционала плотности и представлена на рис. 1, в сопоставлении с электронной структурой кристаллов 3Б МеО (Ме = Мп, Бе, Со, N1). На рис. 1а, б, в, г представлены парциальные плотности электронных состояний квазидвумерных систем 2Б МеО (Ме = Мп, Бе, Со, N1) для электронных подсистем со спином вверх и спином вниз. Картина распределения электронных плотностей состояний по энергии для 2Б МпО (рис. 1а) принципиально не отличается от 3Б МпО (рис. 1д). Немного изменяется относительное положение уровня Ферми, который сдвигается в область более высоких энергий. Анализ распределения парциальных плотностей состояний по энергии для 2Б-оксидов,
представленные для атомов интерфейса (верхнего слоя), позволяет отметить, что аналогично 3Б-системам в ряду МпО - БеО - СоО - №О наблюдается смещение 3d-состояний металла со спином вниз N1 относительно уровня Ферми в область низких энергий. Отмеченное определенным образом влияет на намагниченность (выравнивание магнитных моментов) в двумерных системах оксидов переходных 3d-металлов.
Рис.1. - Парциальные плотности электронных состояний для 2Б-оксидов: МпО (а), БеО (б), СоО (в),№О (г) и 3Б-оксидов: МпО (д), БеО (е), СоО (ж), №О (з) в ферромагнитном упорядочении.
В частности, карты распределения электронной плотности в квазидвумерных системах 2Б МеО (Ме = Мп, Бе, Со, №), представленные на рис. 2, иллюстрируют характер распределения плотности заряда между атомами в системе. В центре каждой карты расположен атом кислорода, а ближайшие к нему - атомы соответствующего металла. Из приведенных карт видно, что преобладающим для всех соединений является ионный тип связи, о чем свидетельствует высокая концентрация электронной плотности вокруг атомов. По сравнению с 3Б системами наблюдается некоторое нарушение симметрии в распределении плотности, что связано с перераспределением последней в связи с наличием поверхности интерфейса.
Рис.2. - Карты распределения электронной плотности для 2Б-оксидов: МпО (а), БеО (б), СоО (в),№О (г) в ферромагнитном упорядочении
Отмеченные выше эффекты могут быть связаны с заполнением d-оболочки металла при увеличении атомного номера. Следует также отметить расщепление 3d-состояний со спином вниз N1 на две и даже на три полосы, как например, для СоО и №О соответственно. Данные о парциальных зарядах валентных электронов на атомах 2Б-оксидов и интегральный спиновой магнитный момент (ММ) на каждом атоме представлены в табл. 2. Анализ представленных в таблице данных позволяет отметить, с одной стороны, возрастание парциального заряда 3d-электронов металла со спином вниз N4 в ряду МпО - БеО - СоО - №О, а с другой - уменьшение магнитного момента на атомах металла. Следует отметить достаточно заметный магнитный момент на атомах кислорода в оксидах железа и никеля. Полученные в настоящей работе данные о парциальных зарядах валентных электронов на атомах в объёмных кристаллах оксидов переходных 3d-металлов и интегральный спиновой магнитный момент на каждом атоме представлены в табл. 3.
Таблица №2.
Парциальные заряды валентных электронов и интегральный спиновой магнитный момент (ММ) на атомах 2Б-оксидов для спиновых подсистем ^(Е) и Щ(Е) в ферромагнитном упорядочении.
MnO slab FeO slab CoO slab NiO slab
Mn O Fe O Co O Ni O
si 1.1775 0.8604 1.1839 0.864 9 0.2445 0.9366 0.2462 0.9418
sj 1.1256 0.8334 1.1463 0.831 3 0.2102 0.9014 0.2412 0.9126
pi 3.4103 2.4019 3.4023 2.567 1 2.5361 2.7065
Pi 3.3704 2.3818 3.3926 2.277 - 2.4910 - 2.4731
8
dT 4.8153 - 4.9688 - 4.8717 - 4.9352 -
d| 0.4569 - 1.3215 - 2.6809 - 3.4203 -
ММ, 4.4502 0.0472 3.6945 0.322 9 2.2252 0.0803 1.5200 0.2627
Таблица №3.
Парциальные заряды валентных электронов и спиновой магнитный момент на атомах 3Б-оксидов переходных 3d-металлов для спиновых подсистем ^(Е) и N |(Б) в ферромагнитном упорядочении.
МпО БеО СоО №О
Мп О Бе О Со О № О
вТ 1.1553 0.860 3 1.158 7 0.8568 0.2329 0.940 9 0.236 8 0.9383
1.1110 0.807 7 1.127 4 0.8134 0.2094 0.884 3 0.236 3 0.8982
рТ 3.4352 2.539 0 3.425 7 2.5748 0.0000 2.594 0 0.000 0 2.7648
р! 3.3967 2.331 3 3.416 8 2.2852 0.0000 2.460 6 0.000 0 2.4597
dТ 4.9156 4.976 5 4.8125 4.972 2
d| 0.4219 1.332 7 2.7704 3.364 6
ММ, № 4.5764 0.260 2 3.683 9 0.3329 2.0656 0.190 0 1.608 1 0.3451
Из табл. 2 можно видеть, что, как и следовало ожидать, наибольший вклад в намагниченность соединений вносят 3d-состояния металла. Также наблюдается перетекание заряда с атомов металла на атомы кислорода, например, в Мп (14.432), Бе (15.434), Со (16.023), N1 (16.808). Наибольшее перетекание заряда наблюдается у №О, что указывает на преобладание ионного типа связи.
Из сопоставления данных, приведенных в табл. 2 и 3, можно заключить, что колебания намагниченностей в двумерных системах рассмотренных оксидов по сравнению с 3Б-системами минимальны. Однако наблюдается сильные колебания спиновой поляризации в слое предшествующем интерфейсу для Со (0.98 ^В/се11) и Мп (3.6 ^В/се11). Эти колебания влияют также на магнитный момент на кислороде.
Таким образом, расчёты из первых принципов показали, что при синтезе 2Б МеО (Ме = Мп, Бе, Со, N1) наблюдается перестройка зонной структуры, ответственная за электромагнитные свойства данной системы. Рассчитанные значения магнитных моментов на атомах металла в рассмотренных ультратонких оксидных плёнках представлены в табл. 4 в сопоставлении с кристаллическим состоянием (3Б), известными экспериментальными и теоретическим данными. на основе разных аппроксимаций (ЬББА, ООА). Анализ табл. 4 показывает, что использование нелокального обменно-корреляционного функционала в параметризации РВЕ оказалось оправданным.
Таблица №4.
Магнетизм ультратонких оксидных пленок (2Б) и кристаллов (3Б) типа МеО (Ме = Мп, Бе, Со, N1) в ферромагнитном упорядочении.
Фаза Магнитный момент на атоме металла,
РВЕ ЬББА [8] ООА [8] Эксперимент [11]
МпО 2Б 4.45 - - -
3Б 4.58 3.72 4.15 4.58; 4.892
БеО 2Б 3.70 - - -
3Б 3.68 3.33 3.38 3.32
СоО 2Б 2.23 - - -
3Б 2.07 2.36 2.42 3.35; 3.8
№О 2Б 1.52 - - -
3Б 1.61 1.21 1.38 1.64; 1.9
Заключение
С помощью спин-поляризованного метода псевдопотенциала в рамках теории функционала плотности (БЕТ) исследована электронная структура, выполнен прогноз электронных и магнитных свойств в квазидвумерных системах 2Б МеО (Ме = Мп, Бе, Со, №) с ферромагнитным типом упорядочения.
Установлен анизотропный характер электронной структуры в направлении перпендикулярном полосе 2Б МеО (Ме = Мп, Бе, Со, №), что может обуславливать различие в магнитных состояниях атомов металла и кислорода. Особенности строения электронной энергетической структуры прифермиевских зон в двумерных системах должны оказывать влияние на электронные и магнитные свойства данных материалов. Выполненные оценки магнитного момента на атоме переходного 3d-металла в криталлах оксидов согласуются с экспериментом. Обнаружен эффект спонтанной спиновой поляризации Ме 3й и 2р-зон атомов кислорода в квазидвумерных системах 2Б МеО (Ме = Мп, Бе, Со, №), играющий, на наш взгляд, доминирующую роль в образовании магнитных моментов на атомах переходного 3d-металла и кислорода, что согласуется с известными физическими представлениями.
Литература
1. Müller F., R. de Masi, Reinicke D., Steiner P., Hüfner S., Stöwe K. Epitaxial Growth of Mn0/Ag(001) Film // Surface Science. 2002. Vol.520: pp. 158-172.
2. Csiszar S.I., Haverkort M.W., Hu Z., Tanaka A., Hsieh H.H., Lin H.-J., Chen C.T., Hibma T., Tjeng L. H. Controlling orbital moment and spin orientation in CoO layers by strain // Physical Review Letters. 2005. Vol.95. : pp. 187205-1-4.
3. Chassé A., Langheirch Ch., Müller F., Hüfner F. Growth and structure of thin MnO films on Ag(001) in dependence on film thickness // Surface Science. 2008. Vol.602. : pp. 597 - 606.
4. Nagel M., Biswas I., Nagel P. Pellegrin E., Schuppler S., Peisert H, Chassé T. Ultrathin transition-metal oxide films: Thickness dependence of the electronic structure and local geometry in MnO // Physical Review B -2007. Vol. 75: pp. 195426-1 - 6.
5. Ilyasov V., Meshi B., Ryzhkin A., Ershov I., Nikiforov I., Ilyasov A. Materials for spintronics: magnetic and transport properties of ultrathin (monolayer graphene)/Mn0(001) and Mn0(001) films // Journal of Modern Physics 2011. Vol. 2. : pp. 1120-1135.
6. Илясов В.В., Великохатцкий Д.А., Ершов И.В., Никифоров И.Я., Жданова Т.П. Первопринципные расчеты параметров химической связи и зонной структуры двумерной системы графен/Mn0(001) // Журнал структурной химии. 2011. №5. С. 879-890.
7. Ершов И.В., Попова И.Г., Илясов В.В. Зонная структура 2D оксидов переходных металлов группы железа // Труды 2 международного молодежного симпозиума «Физика бессвинцовых пьезоактивных и родственных материалов LFPM-2013».- С. 192-194.
8. Dufek P., Blaha P., Sliwko V., Schwarz K. Generalized-gradient-approximation description of band splittings in transition-metal oxides and fluorides. //Phys. Rev. B. 1994-I. Vol. 49. : pp. 10170-10175.
9. Каспржицкий А.С., Лазоренко Г.И., Явна В.А. Моделирование ab initio электронной структуры слоистых алюмосиликатов. // Инженерный вестник Дона, 2013, №3 URL: ivdon.ru/magazine/archive/n3y2013/1861.
10. Положенцев О.Е., Сухарина Г.Б., Гуда А.А., Шаповалов В.В.,. Подковырина Ю.С, Чайников А.П., Бугаев А.Л., Поль А. Наноразмерная атомная и электронная структура наноструктурированного конденсированного материала для возобновляемых источников тока на основе нанокомпозита V2O5/Fe/LiF в цикле зарядка-разрядка. // Инженерный вестник Дона, 2013, №4 URL: ivdon.ru/magazine/archive/n4y2013/2016.
11. Fang Z., Solovyev I.V., Sawada H., Terakura K. First-principles study on electronic structures and phase stability of MnO and FeO under high pressure. // Phys. Rev. B. - 1999-II. Vol. 59. : pp. 762-774.
References
1. Müller F., R. de Masi, Reinicke D., Steiner P., Hüfner S., Stöwe K. Surface Science. 2002. Vol.520. : pp. 158-172.
2. Csiszar S.I., Haverkort M.W., Hu Z., Tanaka A., Hsieh H.H., Lin H.-J., Chen C.T., Hibma T., Tjeng L. H. Physical Review Letters 2005. Vol. 95: pp. 187205-1-4.
3. Chassé A., Langheirch Ch., Müller F., Hüfner F. Surface Science - 2008. Vol.602: pp. 597 - 606.
4. Nagel M., Biswas I., Nagel P. Pellegrin E., Schuppler S., Peisert H, Chassé T. Physical Review B - 2007. Vol. 75: pp. 195426-1 - 6.
5. Ilyasov V., Meshi B., Ryzhkin A., Ershov I., Nikiforov I., Ilyasov A. Materials for spintronics: magnetic and transport properties of ultrathin (monolayer graphene)/MnO(001) and Mn0(001) films // Journal of Modern Physics - 2011. Vol. 2: pp. 1120-1135.
6. Ilyasov V.V., Velikohazkiy D.A., Ershov I.V., Nikiforov I.Ya., Zhdanova T.P. Jurnal structurnoy himii. 2011. №5. pp. 879-890.
7. Ershov I.V., Popova I.G., Ilyasov V.V. 2 mezhdunarodnyi molodezhnyi simpozium «Lead-free piezoceramic materials LFPM-2013».- pp. 192-194.
8. Dufek P., Blaha P., Sliwko V., Schwarz K. Phys. Rev. B. - 1994-I. - Vol. 49. : pp. 10170-10175.
9. Kasprjitzkyi A.A., Lasorenko G.I., Yavna V.A. Inzenernyj vestnik Dona (Rus), 2013, №3 URL: ivdon.ru/magazine/archive/n3y2013/1861.
10. Polozhentsev O.E., Suharina G.B., Guda A.A., Shapovalov V.V., Podkovyrina Yu.S, Chainikov A.P., Bugaev A.L., Pol' A. Inzenernyj vestnik Dona (Rus), 2013, №4 URL: ivdon.ru/magazine/archive/n4y2013/2016.
11. Fang Z., Solovyev I.V., Sawada H., Terakura K. Phys. Rev. B. - 1999-II. -Vol. 59. : pp. 762-774.
