Известия РАН. Серия физическая, 2022, T. 86, № 2, стр. 158-164
Моделирование поведения и магнитных свойств спин-вентильных наноструктур
Е. В. Дроворуб 1, *, В. В. Прудников 1, П. В. Прудников 1
1 Федеральное государственное бюджетное образовательное учреждение высшего образования
Омский государственный университет имени Ф.М. Достоевского
Омск, Россия
* E-mail: drovorub.egor@gmail.com
Поступила в редакцию 20.09.2021
После доработки 11.10.2021
Принята к публикации 22.10.2021
- EDN: FDMPTT
- DOI: 10.31857/S0367676522020089
Аннотация
Реализовано Монте-Карло описание магнитных свойств спин-вентильных структур от температуры и внешнего магнитного поля. Исследованы зависимости гистерезисных явлений от типа анизотропии в структурах, толщины магнитных пленок, величин внутрислоевого и межслоевого обменного взаимодействия.
ВВЕДЕНИЕ
Тенденцией современных исследований в области магнетизма стало смещение фундаментальных интересов от изучения свойств объемных образцов к изучению наноразмерных структур из ультратонких магнитных пленок [1, 2]. В данных низкоразмерных системах возрастает роль поверхностных эффектов, а также процессов, происходящих на границе раздела фаз, которые определяют такие свойства наноразмерных устройств, как гигантское магнитное сопротивление (ГМС) [3–5] и туннельное магнитосопротивление (ТМС) [6–8]. Магнитные структуры с эффектами ГМС и ТМС нашли широчайшее применение в сенсорах магнитного поля, считывающих головках жестких дисков и энергонезависимой магниторезистивной памяти [9, 10].
Синтетические антиферромагнитные наноструктуры с эффектами ГМС состоят из ферромагнитных слоев, разделенных немагнитной металлической прослойкой. Выбор толщины такой прослойки обеспечивает за счет дальнодействующего и осциллирующего межслоевого обменного взаимодействия эффективную антиферромагнитную связь между соседними магнитными слоями. Недостатком таких структур является большое поле насыщения (Hs ≈ 20 кЭ), необходимое для параллельной ориентации намагниченностей слоев и осуществления ГМС.
Заметным прогрессом в деле совершенствования структур с ГМС можно рассматривать создание спин-вентильных систем, которые, как и мультислойные структуры состоят из магнитных слоев, разделенных немагнитной металлической прослойкой. Однако направление намагниченности одного из ферромагнитных слоев в структуре закрепляется взаимодействием с введенным дополнительным слоем антиферромагнетика. Для ослабления величины обменной связи между магнитными слоями толщину немагнитной проводящей прослойки из благородных металлов (Cu, Ag или Au) подбирают достаточно большой (1.5–5 нм), чтобы ориентация намагниченности в магнитном слое, незакрепленном связью с антиферромагнетиком, могла свободно изменяться под действием слабого внешнего магнитного поля. Преимуществом спин-вентильных структур являются низкие поля насыщения (Hs = 5–50 Э), значительно меньшие, чем в мультислойных магнитных структурах. Это свойство спиновых вентилей делает их техническое использование наиболее выгодным. Введенные в 1991 г. [11], спин-вентильные структуры находят широкое применение, например, в конструкциях считывающих головок для жестких дисков с плотностью записи более 100 Гбайт/кв. дюйм. В спиновых вентилях магнитные слои изготавливают из мягкого ферромагнетика, например, сплава Ni80Fe20 (пермаллой) толщиной 1.5–5 нм. В качестве материала для антиферромагнитной пленки используют сплавы FeMn, NiO или IrMn [11–13]. Структура спинового вентиля, например, Ni80Fe20/Cu/Ni80Fe20/FeMn характеризуется тем, что намагниченность изолированного медью слоя Ni80Fe20 может свободно вращаться под действием внешнего магнитного поля относительно намагниченности слоя Ni80Fe20, чье направление закреплено обменной связью с антиферромагнетиком FeMn.
В настоящее время разработаны различные виды усложненных спин-вентильных структур [14], которые по сравнению со спиновым вентилем, введенным в [11], характеризуются наличием дополнительной трехслойной антиферромагнитной системы Co/Ru/Co, усиливающей действие антиферромагнетика FeMn на свободно перемагничивающийся слой ферромагнетика. В данной работе ставится целью осуществить моделирование магнитных свойств такой усложненной спин-вентильной структуры с выделением влияния магнитной анизотропии, величин внутрислоевого и межслоевого обменного взаимодействия на гистерезисные эффекты в спиновом вентиле при изменении толщины наноразмерных ферромагнитных пленок.
МОДЕЛИ И МЕТОДЫ ОПИСАНИЯ УСЛОЖНЕННОЙ СПИН-ВЕНТИЛЬНОЙ СТРУКТУРЫ
На рис. 1а изображена схематическая модель спин-вентильной структуры [11], в которой обычно применяются комбинированные ферромагнитные слои из сплава Ni80Fe20 и Co для повышения в них температуры магнитного упорядочения и увеличения магнитосопротивления. Отметим, что в слое ферромагнетика данной структуры, спаренном со слоем антиферромагнетика и закрепляющем ориентацию намагниченности слоя, под влиянием обменного взаимодействия со стороны антиферромагнетика происходит сдвиг петли гистерезиса на величину Hex, получившего название поля обменного сдвига [11]. На рис. 1б представлено изображение усложненной спин-вентильной структуры, использующей влияние обменного взаимодействия через пленку, например, рутения для образования дополнительной трехслойной антиферромагнитной системы Co/Ru/Co, которая совместно со слоем антиферромагнетика более эффективно закрепляет ориентацию намагниченности в слое кобальта Co, ближайшем к слою меди Cu.
В ультратонких ферромагнитных слоях Co и Ni80Fe20, контактирующих со слоями меди или рутения, намагниченность, как правило, ориентируется в плоскости пленок. Однако в настоящее время в технологических применениях ультратонких магнитных пленок особое внимание начинает уделяться структурам с намагниченностями, перпендикулярными к плоскости пленок. Это обусловлено тем, что целый ряд активных элементов спинтроники с эффектами ГМС, в частности, считывающие магнитные головки с жестких дисков с перпендикулярными к плоскости пленок магнитными конфигурациями демонстрируют значительно лучшие характеристики [15, 16] по сравнению со структурами с намагниченностью в плоскости пленок. Реальными примерами таких магнитных систем являются структуры Pt/Co/Ir/Co/Pt и Pt/Co/Cu/Co/Pt с ультратонкими пленками платины, нанесенными на пленки кобальта [17]. В данной работе мы будем изучать магнитные свойства наноструктур с намагниченностью, ориентированной как в плоскости пленок, так и перпендикулярной плоскости пленок.
В работах [18–20] нами был реализован теоретический подход, который показал, что магнитные свойства ультратонких пленок переходных металлов на основе Fe, Co и Ni, контактирующих с немагнитной металлической подложкой, адекватно описываются анизотропной моделью Гейзенберга. При этом, обменное взаимодействие ферромагнитных пленок в спин-вентильных структурах с намагниченностью, перпендикулярной плоскости ферромагнитных пленок, удобно характеризовать гамильтонианом вида [20]:
(1)
$\begin{gathered} H = - \sum\limits_{\left\langle {i,j} \right\rangle } {{{J}_{{ij}}}\{ (1 - {{\Delta }_{1}}(N))(S_{i}^{x}S_{j}^{x}} + S_{i}^{y}S_{j}^{y}) + \\ + \,\,\,S_{i}^{z}S_{j}^{z}\} - h\sum\limits_i {S_{i}^{z},} \\ \end{gathered} $а в структурах с намагниченностью в плоскости пленок гамильтонианом
(2)
$\begin{gathered} H = - \sum\limits_{\left\langle {i,j} \right\rangle } {{{J}_{{ij}}}\{ S_{i}^{x}S_{j}^{x}} + 0.8S_{i}^{y}S_{j}^{x} + \\ + \,\,(1 - {{\Delta }_{2}}(N))S_{i}^{z}S_{j}^{z}\} - h\sum\limits_i {S_{i}^{x}.} \\ \end{gathered} $В данных выражениях ${{\vec {S}}_{l}}$ = $\left( {S_{i}^{x},S_{i}^{y},S_{i}^{z}} \right)$ – классический трехмерный единичный вектор спина, зафиксированный в i-ом узле ГЦК-решетки ферромагнитной пленки (кобальта), Δ1,2(N) – параметры, учитывающие эффективное влияние анизотропии, создаваемой кристаллическим полем подложки, на магнитные свойства пленки с толщиной N в единицах монослоев. Параметром Δ1(N) мы характеризуем анизотропию типа “легкой оси”, перпендикулярной плоскости пленки, а параметром Δ2(N) – анизотропию типа “легкой плоскости”. Зависимости для параметров Δ1,2(N) брались из работ [19, 20]. Параметр h = gμBH характеризует влияние внешнего магнитного поля, ориентированного вдоль “легкой оси” (оси z) в случае (1) или в плоскости пленки с введенной слабой анизотропией вдоль оси x (случай (2)) для снятия вырождения.
Методика моделирования критических свойств ультратонких гейзенберговских ферромагнитных пленок и мультислойных структур на их основе была отработана нами ранее и представлена в работах [18–22]. В работе [23] нами осуществлялось моделирование поведения простой спин-вентильной структуры с акцентом на разработку методики по расчету ее магнитосопротивления. В данной работе, ориентированной на рассмотрение магнитных свойств усложненной спин-вентильной структуры (рис. 1б), моделирование проводилось для пленок с размерами L × L × N1, 2, 3 и наложенными периодическими граничными условиями в плоскости пленки. Значения обменных интегралов внутрислоевого взаимодействия задавались как J1/kBT = 1, J01/J1 = 0.4, J02/J1 = 0.75, а межслоевого взаимодействия пленок (1 и 2 сверху) J2/J1 = 0.01. Температура T системы измеряется при этом в единицах обменного интеграла J1/kB. Положительный знак обменного интеграла J2 указывает на реализацию ферромагнитной (параллельной) ориентации намагниченностей ферромагнитных пленок относительно друг друга в отсутствие внешнего магнитного поля, что достигается соответствующим подбором толщины немагнитной прослойки. Для моделирования свойств антиферромагнетика вводится отрицательный по знаку обменный интеграл внутрислоевого межподрешеточного взаимодействия J5/J1 = –2.0 с J1 = 1 для описания взаимодействия внутри ферромагнитных подрешеток. Данный слой антиферромагнетика фиксирует ориентацию намагниченности спаренного с ним слоя ферромагнетика за счет обменного взаимодействия с константой J4/J1 = –2.0. В модели усложненного спинового вентиля (рис. 1б) вводится дополнительная трехслойная антиферромагнитная система, образованная двумя ферромагнитными пленками (2 и 3 сверху), разделенными пленкой немагнитного металла, и характеризуемая межслоевым взаимодействием с J3/J1 = –3.0.
РЕЗУЛЬТАТЫ РАСЧЕТОВ И ИХ АНАЛИЗ
В данной работе нами осуществлялся расчет температурной зависимости приведенных намагниченностей пленок (на спин)
(3)
${{M}_{{n,\alpha }}} = \left\langle {\frac{1}{{{{N}_{n}}{{L}^{2}}}}\sum\limits_{i = 1}^{{{N}_{n}}{{L}^{2}}} {S_{i}^{\alpha }} } \right\rangle \,\,\,\,(\alpha = z,x)$для структур с различными толщинами пленок Nn при L = 32 в отсутствие внешнего магнитного поля с проявлением в них различного типа магнитной анизотропии. Это позволяло определить температуры спонтанного магнитного упорядочения для различных пленок в рассматриваемых структурах. В качестве примера на рис. 2а, 2б представлены температурные зависимости намагниченности пленок для структур с анизотропией “легкая плоскость” с толщинами N1 = N2 = N3 = 3 МС (рис. 2а) и анизотропией “легкая ось” с N1 = N2 = = N3 = 5 МС (рис. 2б). Из результатов расчета температурной зависимости спонтанной намагниченности пленок видно, что температуры магнитного упорядочения в 1-ой, 2-ой и 3-ей пленках Tcn данных структур отличаются друг от друга прежде всего за счет влияния внутрислоевого и межслоевого взаимодействий, характеризуемых различными значениями обменных интегралов, а также зависят от толщины этих пленок. В проведенных исследованиях толщину слоя антиферромагнетика мы брали фиксированной с NAF =15 МС.
С целью исследования гистерезисных явлений в спин-вентильных структурах нами было осуществлено моделирование их магнитных свойств при циклическом изменении внешнего магнитного поля, характеризуемого параметром h, в интервале его изменения от h = –2 до h = 2 и обратно. Для изучения выбирались температуры близкие, но ниже температуры Tc1 реализации ферромагнитного упорядочения в 1-ой пленке для осуществления перемагничивания в ней при слабых магнитных полях с минимальными гистерезисными эффектами. Полученные графики петель гистерезиса для структур с анизотропией “легкая плоскость” представлены на рис. 3 как для отдельных ферромагнитных пленок с толщинами N1 = N2 = N3 = 3 МС (рис. 3а) при температуре T = 0.7 J1/kB и с толщинами N1 = N2 = N3 = 7 МС (рис. 3в) при температуре T = 1.2 J1/kB, так и для всей структуры в целом на рис. 3б и 3г. Для данных структур магнитное поле применялось вдоль оси х в плоскости пленок в соответствии с применяемым для моделирования гамильтонианом (2) и введенной в нем слабой анизотропией вдоль оси x для снятия вырождения. Поэтому на рис. 3 представлены изменения компоненты намагниченности Mx от циклического изменения поля h как для отдельных пленок Mxn, так и намагниченности Mx всей структуры в целом.
Анализ графиков, представленных на рис. 3а и 3в, показывает, что в интервале изменения внешнего магнитного поля, определяемого |h| ≤ 2, ориентация намагниченности Mx3 в 3-ей пленке не меняется за счет ее закрепления межслоевым обменным взаимодействием J4/J1 = –2.0 со слоем антиферромагнетика. При выбранных температурах T = 0.7 J1/kB для структуры с N1 = N2 = N3 = = 3 МС и T = 1.2 J1/kB для структуры с N1 = N2 = = N3 = 7 МС, а также слабого межслоевого взаимодействия пленок 1-ой и 2-ой ферромагнитных пленок с J2/J1 = 0.01, достигаемого за счет достаточно толстого (≈25 А) немагнитного слоя меди, намагниченность 1-ой пленки легко переориентируется уже в слабых магнитных полях. Наиболее интересное поведение демонстрирует намагниченность 2-ой пленки, петля гистерезиса которой демонстрирует сильное обменное смещение в область полей, характеризуемых h ≈ 1.5, для структуры с N1 = N2 = N3 = 3 МС и более слабое обменное смещение с h ≈ 0.5 для структуры с N1 = = N2 = N3 = 7 МС при тех же параметрах межслоевого и внутрислоевого обменного взаимодействия. Данные эффекты обменного смещения обусловлены прежде всего введенной в структуру спинового вентиля трехслойной антиферромагнитной системой из 2-ой и 3-ей ферромагнитных пленок с прослойкой немагнитного металла (как правило, рутения с толщиной примерно 7 Å), и характеризуемой межслоевым взаимодействием с J3/J1 = –3.0. На рис. 3б и 3г представлено проявление гистерезисных эффектов в данных структурах в целом. Наблюдаемое обменное смещение петли гистерезиса для 2-ой пленки при узкой петле гистерезиса для 1-ой пленки обеспечивает лучшие условия для проявления эффектов ГМС в данных структурах со стабильными характеристиками.
Нами было также проведено исследование проявления гистерезисных эффектов в структурах с анизотропией “легкая ось”, обеспечивающей ориентацию намагниченности перпендикулярную плоскости пленок. Для моделирования применялся гамильтониан в форме (1). Для сопоставления влияния различных типов анизотропии на явления гистерезиса использовался одинаковый набор параметров, задающих внутрислоевое и межслоевое обменное взаимодействие в структурах. В качестве примера результатов исследования для спиновых вентилей с толщинами пленок N1 = N2 = N3 = 5 МС на рис. 4а–4в приведены графики полевой зависимости компонент намагниченности Mzn (n = 1, 2, 3) для отдельных пленок (рис. 4а) и Mz для всей структуры в целом (рис. 4б) при намагниченности, перпендикулярной плоскости пленок, и графики компонент намагниченности Mxn для отдельных пленок (рис. 4в) и Mx для всей структуры в целом (рис. 4г) при намагниченности, лежащей в плоскости пленок, при температуре T = 1.0 J1/kB. Сопоставление графиков для данных структур показывает их общую схожесть с проявлением заметного обменного смещения петли гистерезиса для 2-ой пленки при узкой петле гистерезиса для 1-ой пленки. Отличия связаны с более широкой петлей гистерезиса у 2-ой пленки для случая с намагниченностью, перпендикулярной плоскости пленок, чем для случая с намагниченностью в плоскости пленок. Это связано с большей величиной намагниченности 2-ой пленки при температуре T = 1.0 J1/kB для случая с намагниченностью, перпендикулярной плоскости пленок. Уменьшение величины интеграла J02 внутрислоевого обменного взаимодействия во 2-ой пленке приводит к уменьшению ее ширины петли гистерезиса, как это показано на рис. 4а и 4б для J02 = 0.5 штриховыми линиями.
ЗАКЛЮЧЕНИЕ
C применением метода Монте-Карло осуществлено численное моделирование усложненной спин-вентильной структуры, используемой в качестве активного элемента приборов спинтроники и характеризуемой эффектами гигантского магнитосопротивления. Для описания магнитных свойств ультратонких пленок, образующих данные структуры, была применена анизотропная модель Гейзенберга как с анизотропией “легкая плоскость”, позволяющая реализовать ориентацию намагниченности, лежащую в плоскости пленок, так и с анизотропией “легкая ось” с реализацией намагниченности, перпендикулярной плоскости пленок. Получена зависимость магнитных характеристик пленок в этих структурах от температуры и внешнего магнитного поля. Исследовано влияние на проявление гистерезисных эффектов толщины магнитных пленок, температуры и величины внутрислоевого и межслоевого обменного взаимодействия. Выделены факторы, позволяющие влиять на условия лучшего проявления в структурах эффектов гигантского магнитосопротивления со стабильными свойствами. Полученные зависимости магнитных характеристик пленок и структур в целом хорошо согласуются с результатами экспериментальных исследований спин-вентильных структур [11, 14, 24].
Исследования выполнены при финансовой поддержке РФФИ (проект № 20-32-70189), Минобрнауки РФ (соглашение № 0741-2020-0002) и Совета по грантам Президента РФ (проект МД-2229.2020.2).
Список литературы
Vaz C.A.F., Bland J.A.C., Lauhoff G. // Rep. Prog. Phys. 2008. V. 71. Art. No. 056501.
Lavrijsen R., Lee J-H., Fernández-Pacheco A. et al. // Nature. 2013. V. 493. P. 647.
Baibich M.N., Broto J.M., Fert A. et al. // Phys. Rev. Lett. 1988. V. 61. No. 21. P. 2472.
Binash G., Grunberg P., Saurenbach F., Zinn W. // Phys. Rev. B. 1989. V. 39. No. 7. P. 4828.
Ферт А. // УФН. 2008. Т. 178. № 12. С. 1336; Fert A. // Phys. Usp. 2008. V. 51. No. 12. P. 1336.
Julliere M. // Phys. Lett. A. 1975. V. 54. No. 3. P. 225.
Miyazaki T., Tezuka N. // J. Magn. Magn. Mater. 1995. V. 139. Art. No. L231.
Sousa R.C., Sun J.J., Soares V. et al. // Appl. Phys. Lett. 1998. V. 73. No. 22. P. 3288.
Prinz G.A. // J. Magn. Magn. Mater. 1999. V. 200. P. 57.
Chappert C., Fert A., Van Dau F. // Nature Mater. 2007. V. 6. No. 11. P. 813.
Dieny B., Speriosu V.S., Parkin S.S.P. et al. // Phys. Rev. B. 1991. V. 43. No. 1. P. 1297.
Anthony T.C., Brug J.A., Zhang S. // IEEE Trans. Mag. 1994. V. 30. No. 11. P. 3819.
Iwasaki H., Saito A.T., Tsutai A., Sahashi M. // IEEE Trans. Mag. 1997. V. 33. No. 5. P. 2875.
Marrows C.H., Stanley F.E., Hickey B.J. // J. Appl. Phys. 2000. V. 87. No. 9. Art. No. 5058.
Richter H.J. // J. Phys. D. 2007. V. 40. No. 9. Art. No. R149.
Suess D., Schrefl T., Fahler S. et al. // Appl. Phys. Lett. 2005. V. 87. Art. No. 12504.
Morgunov R., Hamadeh A., Fache T. et al. // Superlatt. Microstruct. 2017. V. 104. P. 509.
Прудников П.В., Прудников В.В., Медведева М.А. // Письма в ЖЭТФ. 2014. Т. 100. № 7. С. 501; Prudnikov P.V., Prudnikov V.V., Medvedeva M.A. // JETP Lett. 2014. V. 100. No. 7. P. 446.
Prudnikov P.V., Prudnikov V.V., Menshikova M.A., Piskunova N.I. // J. Magn. Magn. Mater. 2015. V. 387. P. 77.
Prudnikov P.V., Prudnikov V.V., Mamonova M.A., Piskunova N.I. // J. Magn. Magn. Matter. 2019. V. 482. P. 201.
Прудников В.В., Прудников П.В., Романовский Д.Е. // Вестн. Омск. ун-та. 2015. № 4. С. 27.
Прудников В.В., Прудников П.В., Борзилов В.О., Сайфутдинов И.К. // Вестн. Омск. ун-та. 2019. Т. 24. № 2. С. 25.
Прудников В.В., Прудников П.В., Самошилова А.А., Христовский К.А. // Вестн. Омск. ун-та. 2020. Т. 25. № 1. С. 22.
Dieny B. // J. Magn. Magn. Mater. 1994. V. 136. P. 335.
Дополнительные материалы отсутствуют.
Инструменты
Известия РАН. Серия физическая