Физика
УДК 538.955 DOI: 10.14529/mmph240106
ЭЛЕКТРОННЫЕ И МАГНИТНЫЕ СВОЙСТВА СПЛАВА ГЕЙСЛЕРА Mn2CrSn: ПЕРВОПРИНЦИПНЫЕ ИССЛЕДОВАНИЯ
М.А. Обамби, М.А. Загребин, В.Д. Бучельников
Челябинский государственный университет, г. Челябинск, Российская Федерация e-mail: [email protected]
Аннотация. С помощью теории функционала плотности исследованы электронные и магнитные свойства стехиометрического сплава Гейслера Mn2CrSn в зависимости от приложенного давления. Показано, что в сплаве могут существовать два магнитных состояния с низким (LMS) и высоким (HMS) магнитным моментом при меньшем и большем объеме элементарной ячейки, разделенные энергетическим барьером. При приложении всестороннего давления 1 ГПа данные состояния вырождаются по энергии. Фаза LMS характеризуется полуметаллическим состоянием, тогда как фаза HMS обладает металлическими свойствами. Константы магнитного обменного взаимодействия, а также температура Кюри в фазе LMS принимают меньшие значения, чем в фазе HMS.
Ключевые слова: сплавы Гейслера; спинтроника; полуметаллические ферромагнетики; обменные взаимодействия; температура Кюри.
Введение
В последние годы ферримагнитные сплавы Гейслера на основе Mn2 привлекают особое внимание благодаря их полуметаллическим свойствам и потенциальному применению в области спинтроники в качестве спиновых инжекторов для магнитной оперативной памяти, а также устройств переноса спинового магнитного момента [1-4]. За последние несколько десятилетий было проведено большое количество теоретических и экспериментальных исследований металлических сплавов Mn2YZ, содержащих 3d, 4d и 5d переходные металлы Y, с целью найти и предсказать новые функциональные полуметаллические соединения. Большинство этих исследований посвящено изучению электронной структуры, магнитных, механических и полуметаллических свойств сплавов Гейслера в рамках расчетов ab initio. Многие из них теоретически предсказаны как полуметаллические ферро- или ферримагнетики с псевдощелью в одном из двух спиновых каналов, что делает их пригодными для технических приложений. Целью данной работы является исследование электронных и магнитных свойств сплава Гейслера Mn2CrSn, в зависимости от приложенного внешнего всестороннего давления.
Детали вычислений
Расчеты свойств сплава Mn2CrSn выполнены с использованием присоединенных проекционных плоских волн (Projector augmented wave - PAW) и гриновских функций Корринга-Кона-Ростокера, реализованных в программных пакетах VASP (Vienna Ab initio Simulation Package) [5, 6] и SPR-KKR (a Spin Polarized Relativistic Korringa-Kohn-Rostoker code) [7] соответственно. В качестве обменно-корреляционного функционала использовалось приближение обобщенного градиента электронной плотности в параметризации Пердью-Бурке-Эрнцерхофа [8].
Исследуемый сплав моделировался с помощью 16-атомной кубической элементарной ячейки регулярной структуры Гейслера L2i (прототип Cu2MnAl), которая состоит из четырех взаимопроникающих ГЦК-подрешеток: две из которых в равной степени заняты атомами Mn в позиции Вайкоффа 8с (1/4, 1/4, 1/4) и (3/4, 3/4, 3/4), а атомы Sn и Cr занимают позиции 4a (0, 0, 0) и 4b (1/2, 1/2, 1/2) соответственно. Следует отметить, что, как было показано в работе [9], инверсная структура ХА (прототип Hg2CuTi) не является выгодной для данного сплава.
В основном магнитном состоянии для соединения Мп2Сгёп энергетически выгодной является ферримагнитная (ФиМ) конфигурация [10]. Упорядочение ФиМ отвечает спиновой конфигурации с обратным магнитным моментом атомов Сг по отношению к Мп.
Оптимизированные с помощью программного пакета УЛ8Р параметры кристаллической решетки в дальнейшем использовались для расчета констант магнитного обменного взаимодействия (.Ту), плотностей электронных состояний (ПЭС), магнитных моментов и температуры Кюри с использованием пакета 8РЯ-ККЯ.
30
20
10
Результаты
В результате проведения геометрической оптимизации кристаллической структуры сплава Mn2CrSn получено, что в отсутствие давления на зависимостях энтальпии от параметра решетки наблюдается наличие глобального и локального минимумов, соответствующих магнитным состояниям с малым целым и высоким дробным магнитным моментом при меньшем и большем объеме кристаллической структуры L2i (рис. 1). Состояние с низким магнитным моментом (low magnetic state - LMS) характеризуется полуметаллическими свойствами, тогда как состояние с высоким магнитным моментом (high magnetic state -HMS) - металлическими свойствами (рис. 2).
В отсутствие давления глобальный минимум энтальпии соответствует фазе HMS c большим равновесным параметром решетки (a0 ~ 6,31 А), в то время как локальный минимум отвечает фазе LMS с меньшим a0 (6,02 А). При приложении всестороннего давления разница энтальпии между фазами LMS и HMS уменьшается и при давлении 1 ГПа фазы становятся вырожденными. Дальнейшее увеличение давления приводит к тому, что в качестве глобального минимума становится фаза LMS (см. рис. 1). Аналогичные результаты были получены в работе [10] для других сплавов Гейслера. Таким образом, внешнее давление позволяет управлять разницей энтальпии между фазами HMS и LMS и достичь вырожденных состояний.
Парциальные и полные ПЭС для фаз LMS и HMS в отсутствии давления представлены на рис. 2, а при наличии давления - на рис. 3.
-10
Mn2CrSn t j f
—0 ГПа \ / —1.0 ГПа 4 J » 1 5 ГПа___________________________________
5.9
6.0
6.1
6.2
6.3
6.4
6.5
Параметр решетки (А)
Рис. 1. Изменение энтальпии АН как функции параметра
решетки для сплава Гейслера Мп^^п относительно минимума энтальпии при меньшем объеме
Энергия (эВ) Энергия (эВ)
Рис. 2. Плотности электронных состояний сплава Mn2CrSn при нормальном давлении в фазе (а) LMS и (б) HMS
Рассмотрим случай состояния ЬМ8 более подробно (рис. 2, а). Состояния между -4 и -2,2 эВ обусловлены в основном ^-электронами атомов Мп и Сг, которые гибридизируются друг с дру-
гом. Полные ПЭС от «-2,2 эВ до Ер обусловлены в основном ^-состояниями атомов Мп. а?-состояния Сг также вносят вклад в суммарные ПЭС, но не так сильно, как атомы Мп. По мере приближения к уровню Ферми в сплаве Mn2CrSn появляется энергетическая псевдощель в состоянии с проекцией спина «вверх» с наименьшей заселенностью носителей вблизи уровня Ферми Ер, а состояния с проекцией спина «вниз» полностью заняты. Таким образом, можно предположить, что состояние LMS Mn2CrSn характеризуется почти полуметаллическим поведением. Что касается ПЭС-состояния HMS, то видно, что ПЭС для этих состояний имеет вид, характерный для металлов (рис. 2, б).
Рассмотрим далее, какое влияние оказывает внешнее давление на ПЭС как в LMS- и HMS-состояниях. На рис. 3 представлены ПЭС, вычисленные для давлений 0 и 1,5 ГПа. Из рис. 3 видно, что давление не оказывает существенного влияние на профиль ПЭС как в состоянии LMS, так и в HMS. Данная тенденция хорошо видна на зависимости спиновой поляризации от давления.
(а1 сс
- МП2Сг8П (ЬМв)
/* //
// у * // /
^--\
N V / -0 ГПа р V ----1.5 ГПа V \
V I \ — I
(б) сс
-2
-4
МП2СГ8П(НМ8) -
\ у* ** \ту' ч\ -0 ГПа Х^ц. / У
----1.5 ГПа V .1.1 / ■
-0.50
-0.25
0.00
0.25
0.50
-0.50
-0.25
0.00
0.25
0.50
Энергия (эВ)
Энергия (эВ)
Рис. 3. Плотности электронных состояний сплава Мп2СгЭп в фазах (а) ЬМЭ и (б) ИМЭ при нормальном давлении и давлении 1,5 ГПа вблизи уровня Ферми
В фазе LMS спиновая поляризация при нормальном давлении составляет 93 % и слабо убывает с приложением внешнего давления. Это говорит о том, что сплав сохраняет свое полуметаллическое состояние. Полученные ПЭС для состояния HMS свидетельствуют о заполнении уровня Ферми, характерного для металла. Спиновая поляризация оценена с использованием следующего выражения:
N Ер )- N Ер )
Р = -
■х 100%.
(1)
N Ер ) + N Ер )
Здесь Щ(Ер), ^|(ЕР) - значения ПЭС на уровне Ферми для электронов с проекцией спина «вверх» и спина «вниз», соответственно. Значения спиновой поляризации представлены в табл. 1. Видно, что давление слабо влияет на степень спиновой поляризации в обеих фазах. Такое поведение обусловлено малыми значениями давления.
Таблица 1
Давление 0 ГПа 1,0 ГПа 1,5 ГПа
LMS 93,7 92,9 92,4
HMS 52,6 47,7 46,9
На рис. 4 представлены рассчитанные параметры магнитного обменного взаимодействия Л] для фаз LMS и HMS в зависимости от приложенного давления. Видно, что все взаимодействия имеют осциллирующий характер и также слабо зависят от приложенного давления. В фазе LMS обменные взаимодействия между ближайшими соседними атомами Мп малы и составляют « 1,5 мэВ в третьей координационной сфере. В фазе HMS обменные взаимодействия Мп-Мп в первой координационной сфере значительно больше по сравнению с фазой LMS и составляют около 30 мэВ. С увеличением расстояния, начиная с 8-й координационной сферы, значения обменных взаимодействий становятся близкими к нулевым. Взаимодействие между атомами Мп и & в первой координационной сфере составляет « 6 и 9 мэВ в фазах LMS и HMS соответственно.
Что касается взаимодействия между атомами Cr, то здесь можно отметить, что во второй координационной сфере в фазе HMS обменное взаимодействие больше в 3 раза по сравнению с таковым в фазе LMS. Увеличение обменного взаимодействия в фазе HMS по сравнению с LMS связано c большим значением магнитного момента. Отметим, что аналогичное поведение обменного взаимодействия в фазах LMS и HMS было показано для серии сплавов Mn2 YSn (Y = Sc, Ti, V) в работе [10].
Рис. 4. Параметры магнитного обменного взаимодействия в зависимости от расстояния между парами атомов сплава Mn2CrSn в фазах (а) LMS и (б) HMS при нормальном давлении и давлении 1,5 ГПа
Рассчитанные обменные интегралы использовались для оценки температуры Кюри TC в модели Гейзенберга в рамках приближения среднего поля (mean field approximation) [11]. Отметим, что значения, полученные с помощью приближения молекулярного поля, обычно завышены по сравнению с экспериментальными значениям, но качественно повторяют наблюдаемые экспериментально зависимости от внешнего давления [12]. В табл. 2 приведены рассчитанные значения полного магнитного момента и температуры Кюри при различной величине внешнего давления. Видно, что фаза HMS характеризуется более высокой TC (более чем в 3 раза) по сравнению с фазой LMS. Большие значения TC в фазе HMS объясняются большими суммарными магнитными моментами и параметрами магнитного обменного взаимодействия Jij. Приложение давления приводит к уменьшению температуры Кюри.
Значения полного магнитного момента
Таблица 2
. (/(в) и температуры Кюри TC (K) в фазах LMS и HMS сплава Mn2CrSn
Фаза LMS HMS
Давление Аголн. (/B) Tc (K) Аголн. (Mb) Tc (K)
0 ГПа -0,585 292,4 2,076 689,0
1,0 ГПа -0,595 274,2 3,181 651,9
1,5 ГПа -0,603 258,3 3,122 625,9
Заключение
В данной работе с помощью теории функционала плотности и метода функций Грина теоретически исследованы электронные и магнитные свойства сплава Гейслера Mn2CrSn. Показано, что в сплаве могут иметь место две магнитные фазы с разной величиной магнитного момента и электронной структурой при разных объемах кристаллической решетки. Фаза с меньшим объемом решетки обладает малым магнитным моментом (фаза LMS), тогда как фаза при большем объеме решетки - большим магнитным моментом (фаза HMS). При нормальном давлении фаза HMS является стабильной (имеет глобальный минимум на зависимости энергии от параметра решетки), а фаза LMS - метастабильной (локальный минимум). Минимумы разнесены по энергии на ~ 10 мэВ. Анализ плотностей электронных состояний вблизи уровня Ферми показывает, что фаза LMS имеет ферримагнитное полуметаллическое поведение с практически 100 % степенью спиновой поляризации, в то время как HMS-состояние характеризуется металлическим поведением. Приложение внешнего гидростатического давления приводит к уменьшению разности
энтальпии между фазами LMS и HMS и при давлении 1 ГПа наблюдается вырождение фаз. Дальнейшее увеличение давления делает LMS-фазу выгодной, а HMS-фазу - метастабильной. Расчеты параметров обменного взаимодействия показали, что в металлической фазе HMS они имеют более высокие значения, чем в LMS-фазе. В связи с этим температура Кюри в HMS-фазе выше, чем в LMS-состоянии. Внешнее гидростатическое давление слабо влияет на параметры магнитного обменного взаимодействия и температуру Кюри.
Полученные результаты показывают, что индуцированный давлением механизм переключения между полуметаллическим и металлическим состояниями с различной намагниченностью и температурой Кюри может быть перспективным для разработки устройств спинтроники, таких как спиновые фильтры, датчики и переключатели.
Работа выполнена при финансовой поддержке гранта Российского научного фонда (Проект № 22-12-20032, https://rscf.ru/project/22-12-20032/, Челябинский государственный университет, Челябинская область).
Литература
1. Doping of Mn2VAl and Mn2VSi Heusler Alloys as a Route to Half-Metallic Anti-Ferromagneism / I. Galanakis, K. Özdogan, E. §a§ioglu, B. Akta§ // Physical Review B. - 2007. - Vol. 75, Iss. 9. -P.092407.
2. Design of Compensated Ferrimagnetic Heusler Alloys for Giant Tunable Exchange Bias / A.K. Nayak, M. Nicklas, S. Chadov et al. // Nature Materials. - 2015. - Vol. 14. - P. 679-684.
3. Basics and Prospective of Magnetic Heusler Compounds / C. Felser, L. Wollmann, S. Chadov et al. // APL Materials. - 2015. - Vol. 3. - P. 041518.
4. Marchenkov, V.V. Unusual Kinetic Properties of Usual Heusler Alloys / V.V. Marchenkov, V.Yu. Irkhin, A.A. Semiannikova // J. Supercond Nov. Magn. - 2022. - Vol. 35. - P. 2153-2168.
5. Kresse, G. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations using a Plane-Wave Basis Set / G. Kresse, J. Furthmüller // Physical Review B. - 1996. - Vol. 54, Iss. 16. - P. 11169.
6. Kresse, G. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method / G. Kresse, D. Joubert // Physical Review B. - 1999. - Vol. 59. - P. 1758.
7. Ebert, H. Calculating Condensed Matter Properties using the KKR-Green's Function Method -Recent Developments and Applications / H. Ebert, D. Ködderitzsch, J. Minar // Reports on Progress in Physics. - 2011. - Vol.74, no. 9. - P. 096501.
8. Perdew, J.P. Generalized Gradient Approximation Made Simple / J.P. Perdew, K. Burke, M. Ernzerhof // Physical Review Letters. - 1996. - Vol.77. - P. 3865-3868. (Erratum Phys. Rev. Lett. 78, 1396 (1997)).
9. Prediction of a Heusler Alloy with Switchable Metal-to-Half-Metal Behavior / V.D. Buchelnikov, V.V. Sokolovskiy, O.N. Miroshkina et al. // Physical Review B. - 2021. - Vol. 103, Iss. 5. - P.054414.
10. Sokolovskiy, V. Ab Initio Prediction of Coexistence of Two Magnetic States in Mn27Sn (7=Sc, Ti, and V) Heusler Alloys Under Applied Pressure / V. Sokolovskiy, M. Zagrebin, D. Baigutlin, V. Buchelnikov // Computational Materials Science. - 2023. - Vol. 228. - P. 112365.
11. Anderson, P.W. Theory of Magnetic Exchange Interactions: Exchange in Insulators and Semiconductors / P.W. Anderson // Solid State Physics. - 1963. - Vol. 14. - P. 99-214.
12. Ab Initio Studies of Phase Transformations in Fe100-xSix / M.A. Zagrebin, M.V. Matyunina, A.B. Koshkin et al. // Physics of the Solid State. - 2020. - Vol. 62, no. 5. - P. 739-743.
Поступила в редакцию 24 января 2024 г.
Сведения об авторах
Обамби Мишэ Антипас - студент, Челябинский государственный университет, г. Челябинск, Российская Федерация.
Загребин Михаил Александрович - доктор физико-математических наук, доцент, профессор, кафедра радиофизики и электроники, Челябинский государственный университет, г. Челябинск, Российская Федерация, e-mail: [email protected].
Бучельников Василий Дмитриевич - доктор физико-математических наук, профессор, заведующий кафедрой физики конденсированного состояния, Челябинский государственный университет, г. Челябинск, Российская Федерация, e-mail: [email protected].
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2024, vol. 16, no. 1, pp. 49-55
DOI: 10.14529/mmph240106
ELECTRONIC AND MAGNETIC PROPERTIES OF THE M^CrSn HEUSLER ALLOY: FIRST-PRINCIPLES INVESTIGATIONS
M.A. Obambi, M.A. Zagrebin, V.D. Buchelnikov
Chelyabinsk State University, Chelyabinsk, Russian Federation e-mail: [email protected]
Abstract. This work studies the electronic and magnetic properties of the stoichiometric Heusler alloy Mn2CrSn depending on the applied pressure, using density functional theory. It is shown that two magnetic states with a low (LMS) and high (HMS) magnetic moment can exist in the alloy at a smaller and larger unit cell volume, separated by an energy barrier. When a confining pressure of 1 GPa is applied, these states become degenerate in energy. The LMS phase is characterized by pseudo semi-metallic and metallic behavior, while the HMS phase has metallic properties. The magnetic exchange interaction constants, and the Curie temperature in the LMS phase take lower values than in the HMS phase.
Keywords: Heusler alloys; spintronics; half-metallic ferromagnets; exchange interactions; Curie temperature.
References
1. Galanakis I., Özdogan K., §a§ioglu E., Akta§ B. Doping of Mn2VAl and Mn2VSi Heusler Alloys as a Route to Half-Metallic Anti-Ferromagneism. Phys. Rev. B, 2007, Vol. 75, Iss. 9, p. 092407. DOI: 10.1103/PhysRevB.75.092407
2. Nayak A.K., Nicklas M., Chadov S., Khuntia P., Shekhar C., Kalache A., Baenitz M., Skourski Y., Guduru V.K., Puri A., Zeitler U., Coey J.M.D., Felser C. Design of Compensated Ferrimagnetic Heusler Alloys for Giant Tunable Exchange Bias. Nature Mater., 2015, Vol. 14, pp. 679-684. DOI: 10.1038/nmat4248
3. Felser C., Wollmann L., Chadov S., Fecher G.H., Parkin S.S.P. Basics and Prospective of Magnetic Heusler Compounds. APL Mater, 2015, Vol. 3, p. 041518. DOI: 10.1063/1.4917387
4. Marchenkov V.V., Irkhin V.Yu., Semiannikova A.A. Unusual Kinetic Properties of Usual Heusler Alloys. J. Supercond Nov. Magn., 2022, Vol. 35, pp. 2153-2168. DOI: 10.1007/s10948-022-06342-1
5. Kresse G., Furthmüller J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations using a Plane-Wave Basis Set. Phys. Rev. B., 1996, Vol. 54, Iss. 16, p. 11169. DOI: 10.1103/PhysRevB.54.11169
6. Kresse G., Joubert D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. RevB., 1999, Vol. 59, Iss. 3, p. 1758. DOI: 10.1103/PhysRevB.59.1758
7. Ebert H., Ködderitzsch D., Minar J. Calculating Condensed Matter Properties using the KKR-Green's Function Method - Recent Developments and Applications. Reports on Progress in Physics, 2011, Vol. 74, no. 9, p. 096501. DOI: 10.1088/0034-4885/74/9/096501
8. Perdew J.P., Burke K., Ernzerhof M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett., 1996, Vol. 77, Iss. 18, pp. 3865-3868. DOI: 10.1103/PhysRevLett.77.3865. (Erratum Phys. Rev. Lett. 78, 1396 (1997)).
9. Buchelnikov V.D., Sokolovskiy V.V., Miroshkina O.N., Baigutlin D.R., Zagrebin M.A., Barbiellini B., Lähderanta E. Prediction of a Heusler Alloy with Switchable Metal-to-Half-Metal Behavior. Phys. Rev. B., 2021, Vol. 103, Iss. 5, p. 054414. DOI: 10.1103/PhysRevB.103.054414
10. Sokolovskiy V., Zagrebin M., Baigutlin D., Buchelnikov V. Ab initio Prediction of Coexistence of Two Magnetic States in Mn27Sn (7=Sc, Ti, and V) Heusler Alloys under Applied Pressure. Comp. Mater. Sci., 2023, Vol. 228, p. 112365. DOI: 10.1016/j.commatsci.2023.112365.
11. Anderson P.W. Theory of Magnetic Exchange Interactions: Exchange in Insulators and Semiconductors. Solid State Phys., 1963, Vol. 14, pp. 99-214. DOI: 10.1016/S0081-1947(08)60260-X
12. Zagrebin M.A., Matyunina M.V., Koshkin A.B., Buchelnikov V.D., Sokolovskiy V.V. Ab Initio Studies of Phase Transformations in Fei00-xSix. Phys. Solid State, 2020, Vol. 62, pp. 739-743. DOI: 10.1134/S1063783420050327
Received January 24, 2024
Information about the authors
Obambi Miche is Student, Chelyabinsk State University, Chelyabinsk, Russian Federation. Zagrebin Mikhail Aleksandrovich is Dr. Sc. (Physics and Mathematics), Associate Professor, Professor, Department of Radiophysics and Electronics, Chelyabinsk State University, Chelyabinsk, Russian Federation, e-mail: [email protected].
Buchelnikov Vasily Dmitrievich is Dr. Sc. (Physics and Mathematics), Professor, Head of the Department of Condensed Matter Physics, Chelyabinsk State University, Chelyabinsk, Russian Federation, e-mail: [email protected].