Известия РАН. Серия физическая, 2022, T. 86, № 9, стр. 1247-1250

Изменение магнитосопротивления спин-туннельного элемента при неоднородном перемагничивании с образованием доменов

В. С. Шевцов 12*, В. В. Амеличев 3, Д. В. Васильев 3, Ю. В. Казаков 3, С. И. Касаткин 2, Д. В. Костюк 3, О. П. Поляков 1, П. А. Поляков 1

1 Федеральное государственное бюджетное образовательное учреждение высшего образования “Московский государственный университет имени М.В. Ломоносова”, физический факультет
Москва, Россия

2 Федеральное государственное бюджетное учреждение науки “Институт проблем управления имени В.А. Трапезникова Российской академии наук”
Москва, Россия

3 Федеральное государственное бюджетное научное учреждение “Научно-производственный комплекс “Технологический центр”
Москва, Россия

* E-mail: vs.shevtcov@physics.msu.ru

Поступила в редакцию 18.04.2022
После доработки 13.05.2022
Принята к публикации 23.05.2022

Полный текст (PDF)

Аннотация

Проведено сравнение типичных экспериментальных зависимостей гигантского магниторезистивного сопротивления спин-туннельного элемента с теоретически кривыми, полученными в рамках классической теории Стонера–Вольфарта когерентного перемагничивания, а также в случае некогерентного перемагничивания с образованием магнитных доменов. Установлено, что перемагничивание посредством образования доменной структуры объясняет различие экспериментальной гистерезисной кривой с теорией Стонера–Вольфарта.

ВВЕДЕНИЕ

В настоящее время широкое применение в спинтронике находят спин-туннельные гетероструктуры [13]. Типичные экспериментальные зависимости изменения сопротивления таких структур, обусловленного наличием гигантского магниторезистивного (ГМР) эффекта при их перемагничивании, имеют вид, изображенный на рис. 1 (пунктирные линии). Теоретический анализ процесса квазистационарного перемагничивания указанных структур требует решения вариационной задачи на нахождение такого распределения вектора намагниченности $\vec {M}\left( {x,y,z} \right)$, которое минимизирует функционал полной магнитной энергии

(1)
$W = \int\limits_V {w\left( {\vec {M}\left( {x,y,z} \right)} \right)~{\text{d}}V} ,$
где $w\left( {\vec {M}\left( {x,y,z} \right)} \right)$ – объемная плотность полной магнитной энергии в заданной точке магнетика с координатами x, y и z. В случае справедливости модели когерентного вращения вектора намагниченности свободного слоя спин-туннельного наноэлемента можно воспользоваться теорией Стонера–Вольфарта [4]. Тогда для плотности магнитной энергии получим следующее выражение [5]:

(2)
$\begin{gathered} w\left( {{\varphi }} \right) = - {{H}_{{0x}}}M\cos {{\varphi }} - {{H}_{{0y}}}M\sin {{\varphi }} + \\ + \,\,\frac{1}{2}{{H}_{{an}}}M{{\sin }^{2}}{{\varphi }} + 2{{\pi }}\left( {{{n}_{x}}{{M}^{2}}{{{\cos }}^{2}}{{\varphi }} + {{n}_{y}}{{M}^{2}}{{{\sin }}^{2}}{{\varphi }}} \right). \\ \end{gathered} $
Рис. 1.

Зависимости изменения ГМР сопротивления R СТМР элемента при H0y = 0 (а) и H0y = 10 Э (б) от проекции внешнего магнитного поля H0x. Пунктирными линиями показаны экспериментальные зависимости, сплошными линиями – теоретические, рассчитанные в рамках модели Стонера–Вольфарта.

Здесь φ – угол поворота вектора намагниченности $\vec {M}$ относительно оси легкого намагничивания (ОЛН), вдоль которой направлена ось X декартовой системы координат. Соответственно координатная ось Y направлена вдоль направления трудного намагничивания в плоскости свободного слоя. Величины H0x и H0y – проекции вектора напряженности ${{\vec {H}}_{0}}$ внешнего магнитного поля на координатные оси X и Y соответственно, Han – эффективное поле одноосной анизотропии, nx и ny – диагональные компоненты тензора размагничивающего фактора вдоль координатных осей X и Y соответственно.

Согласно выражению (1), полная магнитная энергия может иметь два равновесных состояния с наименьшим значением энергии: когда составляющая вектора намагниченности $\vec {M}$ вдоль оси X направлена по этой оси или против. При перемагничивании свободного слоя внешним магнитным полем, направленным вдоль оси X, при некотором критическом значении этого поля происходит изменение ориентации вектора намагниченности. Причем при прямом и обратном ходе этого поля критические значения имеют разные величины, что приводит к гистерезису кривой перемагничивания. Критические значения магнитного поля, согласно теории Стонера–Вольфарта [4], можно получить из условий [5]:

(3)
$\frac{{\partial w\left( {{\varphi }} \right)}}{{\partial {{\varphi }}}} = 0;\,\,\,\,\frac{{{{\partial }^{2}}w\left( {{\varphi }} \right)}}{{\partial {{{{\varphi }}}^{2}}}} = 0.$

На рис. 1 сплошными линиями показаны теоретические зависимости сопротивления R спин-туннельного элемента от проекции H0x внешнего магнитного поля, полученные в рамках этой теории. Качественно эта теория удовлетворительно описывает ход кривой и основные особенности гистерезиса ГМР сопротивления, но количественно экспериментальные графики изменяются более плавно. В этой работе сделана попытка объяснить это расхождение наличием некогерентного изменения вектора намагниченности в свободном слое.

ПАРАМЕТРЫ ИССЛЕДУЕМОГО ОБРАЗЦА

В работе исследуется ГМР сопротивление спин-туннельного магниторезистивного (СТМР) элемента следующей конфигурации: Ta(10 нм)/ CoFe(2 нм)/CoFeB(3 нм)/MgO(3.5 нм)/CoFeB(3 нм)/ Ru(0.85 нм)/CoFe(2 нм)/FeMn(10 нм)/Ta(10 нм). Данный СТМР элемент содержит так называемую синтетическую антиферромагнитную структуру, состоящую из двух противоположно намагниченных слоев (CoFe и CoFeB), связанных косвенным антиферромагнитным взаимодействием через слой рутения (Ru). Это позволяет исключить (или существенно уменьшить) влияние магнитостатического поля фиксирующего магнитного слоя на свободный слой СТМР элемента.

Нанослои СТМР элемента имеют форму эллипсов с полуосями a = 12 мкм и b = 6 мкм. Ось обменной (однонаправленной) анизотропии (ООА) между антиферромагнитным слоем (FeMn) и ферромагнитным слоем (CoFe) и ось легкого намагничивания (ОЛН) свободного слоя CoFeB направлены вдоль большей полуоси a эллипсов. Внешнее магнитное поле ${{\vec {H}}_{0}}$ можно разложить на две проекции: постоянная H0y, направленная вдоль оси трудного намагничивания (т.е. вдоль полуоси b эллипсов), а также H0x, направленная вдоль ОЛН. Изменение проекции H0x внешнего магнитного поля приводит к изменению ориентации вектора намагниченности $\vec {M}$ свободного слоя, в результате чего возникает неоднородное распределение намагниченности.

Физические параметры материала свободного слоя, использующиеся в дальнейших расчетах: намагниченность насыщения Ms = 931 Гс, эффективное поле одноосной анизотропии Han = 15 Э, константа обменного взаимодействия A = 1.3 ⋅ 10–6 эрг/см.

ТЕОРЕТИЧЕСКИЙ РАСЧЕТ

Неоднородное перемагничивание моделировалось в рамках микромагнитной двумерной теории с использованием метода динамического установления, основанного на решении динамического уравнении ЛандауЛифшицаГильберта:

(4)
$\frac{{\partial{ \vec {M}}}}{{\partial t}} = - {{\gamma }}\left[ {\vec {M} \times {{{\vec {H}}}_{{ef}}}} \right] + \frac{{{\alpha }}}{M}\left[ {\vec {M} \times \frac{{\partial{ \vec {M}}}}{{\partial t}}} \right],$
где t – время, γ – гиромагнитное отношение, α – параметр затухания, $\vec {M}$ – вектор намагниченности, ${{\vec {H}}_{{ef}}}$ – эффективное магнитное поле, равное вариационной производной от функционала полной магнитной энергии. Ферромагнитный образец разбивается на N одинаковых прямоугольных ячеек, при этом их размер подбирается таким образом, чтобы в пределах этих ячеек было практически однородное распределение намагниченности $\vec {M}$. Магнитный момент каждой такой ячейки равен ${{\vec {p}}_{{mi}}} = {{\vec {M}}_{i}} \cdot \Delta V,$ где ∆V – объем ячейки. Поместив эквивалентный диполь в центр ячейки можно свести континуальное уравнение Ландау–Лифшица–Гильберта (4) к системе N дискретных уравнений:

(5)
$\frac{{\partial {{{\vec {p}}}_{{mi}}}}}{{\partial t}} = - {{\gamma }}\left[ {{{{\vec {p}}}_{{mi}}} \cdot {{{\vec {H}}}_{{ef~i}}}} \right] + \frac{\alpha }{{{{p}_{{mi}}}}}\left[ {{{{\vec {p}}}_{{mi}}} \cdot \frac{{\partial {{{\vec {p}}}_{{mi}}}}}{{\partial t}}} \right].$

Эта система решается численно при задании начальных и граничных условий. С течением времени система стремится к некоторому стационарному распределению магнитных моментов ${{\vec {p}}_{{mi}}}$ в ферромагнетике. Это происходит за счет диссипативного члена, выраженного последним слагаемым в уравнении (5). Данное распределение является результирующим при заданных условиях.

Для реализации описанного выше численного метода использовался распространенный программный пакет OOMMF [6]. В результате расчетов действительно было зафиксировано возникновение неоднородных магнитных структур (доменов) в свободном слое. Одна из таких структур представлена на рис. 2.

Рис. 2.

Неоднородное распределение намагниченности в свободном слое, рассчитанное с использованием программного пакета OOMMF, при H0y = 0 и H0x = 16 Э.

Расчет зависимости ГМР сопротивления СТМР элемента R от компоненты H0x напряженности внешнего магнитного поля производился на основе полученных распределений намагниченности с использованием известной формулы:

(6)
$\frac{{\Delta R}}{{{{R}_{0}}}} = 0.5\frac{{\Delta {{\rho }}}}{{{\rho }}}\left( {1 - \cos {{\varphi }}} \right),$
где φ – угол между векторами намагниченности фиксированной и свободной соседних магнитных пленок, R0 – минимальное сопротивление СТМР элемента при параллельном расположении векторов намагниченности в соседних слоях (φ = 0), ΔR – изменение сопротивления СТМР элемента, Δρ/ρ – величина ГМР эффекта.

На рис. 3 представлены результаты теоретических расчетов изменения ГМР сопротивления при некогерентном перемагничивании свободного слоя спин-туннельного элемента, которая значительно лучше приближается к экспериментальной кривой.

Рис. 3.

Зависимости изменения ГМР сопротивления R СТМР элемента при H0y = 0 (а) и H0y = 10 Э (б) от величины внешнего магнитного поля H0x. Пунктирными линиями показаны экспериментальные зависимости, сплошными линиями – теоретические, рассчитанные в рамках модели некогерентного перемагничивания с использованием пакета OOMMF.

ЗАКЛЮЧЕНИЕ

Можно отметить, что график на рис. 3 имеет ступенчатый вид, что связано с особенностью неоднородного перемагничивания эллипсоида, а именно с возникновением доменов (рис. 2). Вследствие этого, при изменении внешнего магнитного поля происходит резкое перемагничивание определенных участков эллипсоида, в результате чего направление намагниченности в этих областях меняется практически на противоположное. Было установлена, что величина скачков зависит от выбора размеров ячейки при численном счете, однако дальнейшее уменьшение размеров ячейки в данной работе сопряжено с трудностями в виде кратного увеличения времени счета. Решение проблемы требует гораздо больших вычислительных ресурсов. Тем не менее показано, что расчет на основе модели некогерентного перемагничивания демонстрирует лучшее сходство с экспериментальными данными.

Работа поддержана Фондом развития теоретической физики и математики “БАЗИС”.

Список литературы

  1. Scheike T., Xiang Q., Wen Z. et al. // Appl. Phys. Lett. 2021. V. 118. No. 4. Art. No. 042411.

  2. Limeira V.P.C., Nagamine L.C.C.M., Geshev J. et al. // J. Phys. Cond. Matter. 2019. V. 31. No. 26. Art. No. 265802.

  3. Luong V.S., Nguyena A.T., Hoang Q.K. et al. // J. Sci. Adv. Mater. Devices. 2018. V. 3. No. 4. P. 399.

  4. Stoner E.C., Wohlfarth E.P. // Phil. Trans. Royal Soc. A. 1948. V. 240. No. 826. P. 599.

  5. Поляков О.П., Касаткин С.И., Амеличев В.В. и др. // Изв. РАН. Сер. физ. 2021. Т. 85. № 11. С. 1554; Polyakov O.P., Kasatkin S.I., Amelichev V.V. et al. // Bull. Russ. Acad. Sci. Phys. 2021. V. 85. No. 11. P. 1217.

  6. https://math.nist.gov/oommf.

Дополнительные материалы отсутствуют.