Физика металлов и металловедение, 2022, T. 123, № 3, стр. 299-304

Характер, структура колебаний компонент намагниченности и упругого смещения пленки при амплитудно-модулированном возбуждении

Д. А. Плешев a*, В. С. Власов b, Ф. Ф. Асадуллин a, Л. Н. Котов b, И. А. Чупров b, В. И. Щеглов c

a СПбГЛТУ им. С.М. Кирова
192021 Санкт-Петербург, Институтский переулок, 5 Литер У, Россия

b СГУ им. Питирима Сорокина
167001 Сыктывкар, Октябрьский проспект, 55, Россия

c ИРЭ им В.А. Котельникова РАН
125009 Москва, ул. Моховая 11, корп. 7, Россия

* E-mail: dpleshev@gmail.com

Поступила в редакцию 08.07.2021
После доработки 07.11.2021
Принята к публикации 12.11.2021

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

Аннотация

Рассмотрены колебания магнитоупругой подсистемы ферритовой пленки, возбужденные амплитудно-модулированным переменным магнитным полем. Рассмотрена возможность детектирования колебаний на частоте модуляции за счет магнитострикционных свойств магнитной пленки. Определена зависимость возбуждения магнитной и упругой динамики от материальных параметров и параметров возбуждения, в том числе типа поляризации переменного магнитного поля.

Ключевые слова: детектирование колебаний, магнитные пленки, магнитоакустика, композитные пленки, магнитострикция, нелинейные колебания

ВВЕДЕНИЕ

Применение тонких магнитных периодических структур, композитных и однокомпонентных пленок с минимальными массогабаритными параметрами находят применение в широком спектре радиотехнических устройств [13]. Одним из перспективных материалов для изучения возможностей подобных устройств является железоиттриевый гранат (ЖИГ).

Значительное внимание исследователей привлекает возбуждение спиновых волн за счет магнитоупругости при воздействии на магнитные пленки мощным импульсом света от фемтосекундного лазера в схеме “накачки-зонда” [47]. Достаточно полный обзор представлен в работах [7, 8]. Однако открытым остается вопрос приема таких волн и выделения из СВЧ-сигнала материалов сообщения.

Для приема и регистрации высокочастотных сигналов в традиционной радиотехнике используется метод детектирования на нелинейном элементе, но в СВЧ-диапазоне его применение имеет низкую эффективность вследствие высокой междуэлектродной емкости. Использование магнитоупругих детекторов сигналов позволяет обойти данное ограничение, а использование геометрии тонкой нормально намагниченной пленки позволяет избежать влияния обменных спиновых волн при их параметрическом возбуждении [3].

Общая схема эффекта детектирования в рамках описанной задачи подробно описана в работе [9]. Однако зависимость динамики магнитной и упругой подсистемы от типа поляризации возмущающего переменного магнитного поля, характер параметрических портретов намагниченности и упругого смещения описаны лишь частично.

Настоящая работа посвящена рассмотрению характера и структуры возбуждения колебаний компонент намагниченности и упругого смещения переменным амплитудно-модулированным магнитным полем, а также развитию во времени и анализу параметрических портретов колебаний намагниченности и упругого смещения в различных режимах возбуждения.

ПОСТАНОВКА ЗАДАЧИ

Задача решается в декартовой системе координат, центр которой совпадает с центром пленки. Магнитная пленка толщиной d с магнитными, упругими и магнитоупругими свойствами обладает кубической кристаллографической симметрией (рис. 1). Пленка помещена в постоянное магнитное поле H0, ориентированное вдоль положительного направления оси Oz. Магнитная и упругая подсистемы возбуждаются приложенным в плоскости пленки переменным магнитным полем h0.

Рис. 1.

Геометрия задачи.

Пусть полная плотность энергии пластины U в поле H = {hx; hy; H0} равна сумме магнитной, магнитоупругой и упругой энергий:

(1)
$\begin{gathered} U = - {{M}_{0}}{{h}_{x}}{{m}_{x}} - {{M}_{0}}{{h}_{y}}{{m}_{y}} - \\ - \,\,{{M}_{0}}{{H}_{0}}{{m}_{z}} + 2\pi M_{0}^{2}m_{z}^{2} + \\ + \,\,2{{c}_{{44}}}\left( {u_{{xy}}^{2} + u_{{yz}}^{2} + u_{{zx}}^{2}} \right) + \\ + \,\,2{{B}_{2}}\left( {{{m}_{x}}{{m}_{y}}{{u}_{{xy}}} + {{m}_{y}}{{m}_{z}}{{u}_{{yz}}} + } \right.\left. {{{m}_{z}}{{m}_{x}}{{u}_{{zx}}}} \right), \\ \\ \end{gathered} $
где m = M/M0 – нормированный вектор намагниченности, M0 – намагниченность насыщения; c44 – константа (модуль) упругости; B2 – константа магнитоупругого взаимодействия.

Компоненты тензора деформаций равны:

(2)
${{u}_{{ik}}} = \frac{1}{2}\left( {\frac{{\partial {{u}_{i}}}}{{\partial {{x}_{k}}}} + \frac{{\partial {{u}_{k}}}}{{\partial {{x}_{i}}}}} \right),$
где ui – компоненты упругого смещения при i = x, y, z, xi – координаты x, y, z.

Для снижения сложности решения задачи введем два упрощения:

− прецессия намагниченности однородна в объеме пленки и не зависит от координат;

− продольные колебания упругого смещения вдоль оси Oz не связаны с колебаниями намагниченности, т.е. uz = 0.

При принятых предположениях система уравнений для намагниченности и упругих смещений имеет вид (3)–(4) и решается численно методом Рунге–Кутта 7–8 порядка с контролем точности на каждом шаге [10, 11]:

(3)
$\frac{{\partial {\mathbf{M}}}}{{\partial t}} = - \left| \gamma \right|\left[ {{\mathbf{M}} \times {{{\mathbf{H}}}_{{{\text{eff}}}}}} \right] + \frac{\alpha }{M}\left[ {{\mathbf{M}} \times \frac{{\partial {\mathbf{M}}}}{{\partial t}}} \right],$
(4)
$\frac{{{{\partial }^{2}}{{u}_{{x,y}}}}}{{\partial {{t}^{2}}}} = - 2\beta \frac{{\partial {{u}_{{x,y}}}}}{{\partial t}} + \frac{{{{c}_{{44}}}}}{\rho }\frac{{{{\partial }^{2}}{{u}_{{x,y}}}}}{{\partial {{z}^{2}}}},$
где γ – гиромагнитная постоянная, α – параметр магнитного затухания, ρ – плотность материала пластины, β – параметр упругого затухания.

Граничные условия выглядят следующим образом:

(5)
${{c}_{{44}}}{{\left. {\frac{{\partial {{u}_{{x{\text{,}}y}}}}}{{\partial z}}} \right|}_{{z = \pm {d \mathord{\left/ {\vphantom {d 2}} \right. \kern-0em} 2}}}} = - {{B}_{2}}{{m}_{{x,y}}}{{m}_{z}}.$

Амплитудная модуляция для возмущающего переменного магнитного поля вводится в классическом виде:

(6)
${{h}_{x}} = {{h}_{{0x}}}\left[ {1 + {{A}_{{\text{m}}}}\sin \left( {2\pi Ft} \right)} \right]\sin \left( {2\pi ft} \right),$
где Am – коэффициент (глубина) модуляции. Для неискаженной передачи сигнала Am ≤ 1.

Материальные параметры пленки выбраны характерными для ЖИГ: M0 = 139 Гс, H0 = 2750 Э, параметр диссипации магнитной и упругой подсистемы равны соответственно α = 0.02 и β = 2 × × 108 c–1, константа магнитоупругого взаимодействия B2 = 6.96 × 106 эрг см–3, константа (модуль) упругости с44 = 7.64 × 1011 эрг см–3, напряженность переменного магнитного поля составила h0x = h0y = (0–2) × 103 Э. Частота ФМР равна несущей частоте амплитудно-модулированного поля и составляет 2800 МГц, а частота модуляции равна частоте акустического резонанса при толщине пленки 6.865 мкм, что соответствует 280 МГц. Глубина модуляции Am = 1.

РЕЗУЛЬТАТЫ

Рассмотрим развитие колебаний вектора намагниченности в зависимости от поляризации и амплитуды переменного амплитудно-модулированного магнитного поля (рис. 2).

Рис. 2.

Параметрические портреты и развитие колебаний намагниченности во времени при круговой поляризации переменного магнитного поля (а–з) и при линейной поляризации (и–р). Напряженность переменного магнитного поля h0 = 10 Э – a, б, и, к; h0 = 500 Э – ж, з, п, р; h0 = 1000 Э – в, г, л, м; h0 = 2000 Э – д, е, н, о.

При амплитудах магнитного поля до h0 = 20 Э параметрические портреты колебаний вектора намагниченности представляют собой аттрактор, состоящий из внешней круговой окружности и нескольких вложенных (рис. 2a, 2б, 2и, 2к).

Дальнейшее увеличение амплитуды переменного поля в случае круговой поляризации не вызывает значительных изменений в структуре параметрических портретов и характере колебаний вектора намагниченности, амплитуда которых монотонно возрастает с полем (рис. 2в, 2д, 2ж).

При линейной поляризации ситуация кардинально меняется. При амплитуде переменного поля h0 > 20 Э параметрический портрет колебаний вектора намагниченности начинает испытывать существенные искажения, вызванные модуляцией и сильной нелинейностью при больших амплитудах переменного магнитного поля. Колебания вектора намагниченности приобретают форму эллипса, ориентированного вдоль оси Ox, т.е. вдоль направления поляризации поля.

При амплитуде поля более 300 Э растягивание эллипса вдоль оси поляризации приводит к образованию двух эллиптических траекторий движения вектора намагниченности, соединенных перемычкой, формирование которой завершается к h0 > 800 Э. При этом режим колебаний становится сильно нелинейным, а раскрытие конуса прецессии колебаний намагниченности приближается к 90 град., т.е. прецессия намагниченности происходит в плоскости пленки (рис. 2л, 2м, 2н, 2о).

Дальнейшее увеличение амплитуды переменного поля приводит к усложнению траектории колебаний вектора намагниченности, которые становятся неустойчивыми и переходят в квазигармонический режим (рис. 2п, 2р).

Анализ спектров колебаний намагниченности, представленных на рис. 3, показал, что в случае круговой поляризации переменного поля колебания намагниченности возбуждаются на частоте возбуждения и боковых частотах модуляции, а при линейной поляризации поля уже при амплитудах переменного магнитного поля h0 > 20 Э отмечается появление большого количества комбинационных частот.

Рис. 3.

Спектры колебаний намагниченности при круговой поляризации переменного магнитного поля – a, б, в и при линейной поляризации – г, д, е. Напряженность переменного магнитного поля h0 = 10 Э – a, г, h0 = 500 Э – б, д, h0 = = 1000 Э – в, е.

Отметим, что согласно данным, представленным на рис. 2 и 3, можно сделать вывод о том, что амплитуда колебаний компонент вектора намагниченности при возбуждении магнитной подсистемы пленки линейно поляризованным магнитным полем больше (более 40%), чем в случае круговой поляризации.

Рассмотрим динамику колебаний упругого смещения магнитной пленки при возбуждении амплитудно-модулированным переменным магнитным полем. При этом ограничимся случаем линейной поляризации поля, как наиболее эффективным для накачки магнитной подсистемы.

В первую очередь необходимо обратить внимание на то, что колебания компоненты вектора намагниченности mz имеют четко выраженную составляющую на частоте модуляции (рис. 4а), что с учетом (5) позволяет предположить наличие существенного вклада продольных колебаний намагниченности в поперечные упругие смещения. Таким образом, продольная компонента намагниченности, осциллируя на частоте модуляции, возбуждает поперечные колебания компоненты упругого смещения.

Рис. 4.

Колебания намагниченности mz (а) и упругого смещения ux (б) при амплитуде переменного поля h0 = 500 Э.

Представленные на рис. 5 графики демонстрируют развитие колебаний упругого смещения с увеличением амплитуды переменного поля.

Рис. 5.

Развитие колебаний упругого смещения при линейной поляризации переменного магнитного поля (а–г) и спектры этих колебаний (д–з). Напряженность переменного магнитного поля h0 = 10 Э (a, д), h0 = 500 Э (б, е), h0 = = 1000 Э (в, ж), h0 = 2000 Э (г, з).

При значениях амплитуды переменного поля менее 20 Э, колебания упругого смещения следуют за колебаниями намагниченности в установившемся режиме, а их спектры соответствуют рис. 3г.

Дальнейшее увеличение амплитуды колебаний магнитного поля приводит к возникновению нелинейных режимов колебаний, и при значениях поля h0 > 200 Э увеличивается амплитуда колебаний упругого смещения в области неустойчивости на частоте модуляции, соответствующей частоте акустического резонанса пленки (рис. 5е).

При значениях переменного магнитного поля в диапазоне h0 = 700–1000 Э в спектре колебаний компонент вектора намагниченности происходит увеличение количества комбинационных составляющих в области верхней боковой частоты (рис. 2е), что приводит к перераспределению энергии в спектре колебаний упругого смещения: уменьшению амплитуды колебаний на частоте акустического резонанса, равной частоте модуляции, и росту амплитуды на несущей частоте возбуждения (рис. 5в, 5ж).

После насыщения магнитной подсистемы динамика колебаний упругой подсистемы возвращается к увеличению амплитуды колебаний на частоте модуляции. Однако при амплитуде поля более h0 = 1600 Э начинают возникать квазигармонические колебания упругого смещения, обусловленные сильной нелинейностью системы, что приводит к уменьшению амплитуды колебаний упругого смещения на частотах, близких к несущей частоте переменного магнитного поля.

ЗАКЛЮЧЕНИЕ

В работе рассмотрены колебания магнитной и упругой подсистем магнитной пленки, возбужденные амплитудно-модулированным переменным магнитным полем, а также показана возможность выделения на упругой подсистеме частоты модуляции за счет резонансных магнитоакустических свойств пленки.

Выявлена роль поляризации возмущающего магнитного поля при детектировании. Определено, что эффективное детектирование сигнала на частоте модуляции, совпадающей с частотой акустического резонанса, возможно в нелинейных режимах работы магнитоупругой системы при амплитудах переменного магнитного поля h0 > 200 Э, за исключением диапазона амплитуд магнитного поля, соответствующих насыщению магнитной подсистемы и раскрытию конуса прецессии на угол 90 град.

Показано, что значительное влияние на процесс детектирования оказывает эффективное возбуждение продольной компоненты намагниченности mz, частота биений которой совпадает с частотой модуляции и частотой акустического резонанса пленки.

Исследование выполнено при поддержке Правительства Республики Коми и Российского фонда фундаментальных исследований (грант РФФИ и Правительства Республики Коми, № 20-42-110004, р_а).

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

  1. Kikuchi E. The ultra-sound converters. M.: Mir, 1972.

  2. Comstock R.L., LeCraw R.C. Generation of microwave elastic vibrations in a disk by ferromagnetic resonance // J. Appl. Phys. 1963. V. 34. № 10. P. 3022–3027.

  3. Sementsov D.I., Shuty A.M. Nonlinear regular and stochastic dynamics of magnetization in thin–film structures // Physics Uspekhi. 2007. V. 50. № 8. P. 793.

  4. Dreher L., Weiler M., Pernpeintner M., Huebl H., Gross R., Brandt M.S., Goennen-wein S.T.B. Surface acoustic wave driven ferromagnetic resonance in nickel thin films: theory and experiment // Phys. Rev. B. 2012. V. 86. №13. P. 134415(13).

  5. Bigot J.V., Vomir M. Ultrafast magnetization dynamics of nanostructures // Ann. Phys. (Berlin). 2013. V. 525. № 1–2. P. 2–30.

  6. Walowski J., Münzenberg M. Perspective: Ultrafast magnetism and THz spintronics // J. Appl. Phys. 2016. V. 120. № 14. P. 140901(16).

  7. Wei-Gang Yang, Holger Schmidt. Acoustic control of magnetism toward energy-efficient applications // Appl. Phys. Rev. 2021. V. 8. P. 021304.

  8. Kim S.-K. Micromagnetic computer simulations of spin waves in nanometer-scale patterned magnetic elements // J. Phys. D: Appl. Phys. 2010. V. 43. P. 264004(25).

  9. Власов В.С., Плешев Д.А., Шавров В.Г., Щеглов В.И. Нелинейное детектирование магнитоупругих колебаний в режиме амплитудной модуляции // Журн. радиоэлектроники. 2019. № 3. С. 1–68.

  10. Pleshev D.A., Asadullin F.F., Poleshikov S.M., Vlasov V.S., Kotov L.N., Shavrov V.G., Shcheglov V.I. Investigation of nonlinear dynamics of magnetoelastic oscillations in normal magnetized ferrite plate // Solid State Phenomena. 2015. V. 233–234. P. 471–475.

  11. Pleshev D., Asadullin F., Poleshchikov S., Vlasov V., Kotov L., Shavrov V., Shcheglov V. Conditions for division and multiplication of frequencies in ferrite // EPJ Web of Conferences. Moscow International Symposium on Magnetism, MISM. 2017, 2018. P. 02004.

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