Известия РАН. Серия физическая, 2019, T. 83, № 7, стр. 987-989

Особенности динамики магнитной и упругой подсистем в тонкой ферритовой пленке при магнитострикционном преобразовании частот

Д. А. Плешев 12*, Ф. Ф. Асадуллин 12, Н. А. Оганезова 2, В. С. Власов 2, Л. Н. Котов 2, В. И. Щеглов 3

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

2 Федеральное государственное образовательное учреждение высшего профессионального образования “Сыктывкарский государственный университет имени Питирима Сорокина”
Сыктывкар, Россия

3 Федеральное государственное бюджетное учреждение науки Институт радиотехники и электроники имени В.А. Котельникова Российской академии наук
Москва, Россия

* E-mail: dpleshev@gmail.com

Поступила в редакцию 07.09.2018
После доработки 31.01.2019
Принята к публикации 27.03.2019

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

Аннотация

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

ВВЕДЕНИЕ

Проблема преобразования частот и генерации мощных упругих и магнитных колебаний на сверхвысоких частотах до 100 ГГц представляют важную и актуальную задачу. В настоящее время область сверхбыстрой магнитоакустики является одной из бурно развивающихся областей современной физики [1, 2]. На сверхкоротких временах можно получить мощные упругие импульсы [2]. Однако проблема преобразования частот мощных магнитоупругих колебаний еще не до конца решена, в том числе из-за параметрического распада прецессии на обменные спиновые волны при больших мощностях переменного поля.

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

ГЕОМЕТРИЯ ЗАДАЧИ

Рассмотрим бесконечную магнитную пленку толщины d, которая помещена в постоянное магнитное поле H0, ориентированное перпендикулярно плоскости пленки. Значение поля превышает значение поля размагничивания H0>M0, где M0 – намагниченность насыщения пленки. Задача решается в декартовой системе координат. Плоскость Oxy лежит в плоскости пленки, а центр системы координат O находится в центре пленки. Оси системы координат Ox, Oy и Oz параллельны ребрам кристаллографической ячейки. При решении задачи рассматривались только поперечные упругие колебания.

ОСНОВНЫЕ УРАВНЕНИЯ

Полная энергия U пленки магнетика в постоянном магнитном поле $\vec {H} = \{ 0;0;{{\vec {H}}_{0}}\} $ равна сумме плотностей магнитной, упругой и магнитоупругой энергий:

(1)
$\begin{gathered} U = - {{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}}} + } \right. \\ \left. { + \,\,{{m}_{y}}{{m}_{z}}{{u}_{{yz}}} + {{m}_{z}}{{m}_{x}}{{u}_{{zx}}}} \right), \\ \end{gathered} $
где $\vec {m}$ = ${{\vec {M}} \mathord{\left/ {\vphantom {{\vec {M}} {{{M}_{0}}}}} \right. \kern-0em} {{{M}_{0}}}}$ – нормализованный вектор намагниченности, ${{M}_{0}}$ – намагниченность насыщения пленки, ${{u}_{{ij}}}$ – компоненты тензора деформаций, ${{c}_{{44}}}$ – упругая константа, ${{B}_{2}}$ – магнитоупругая константа.

Колебания намагниченности и упругого смещения определяются уравнением Ландау–Лифшица с диссипативным членом в форме Гильберта и уравнением для компонент вектора механических смещений:

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

Эффективные поля определяются следующим образом:

(4)
${{\vec {H}}_{{eff}}} = - \frac{1}{{{{M}_{0}}}}\frac{{\partial U}}{{\partial{ \vec {m}}}}.$

Граничные условия:

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

Система уравнений (2)–(3) с граничным условием (5) решалась численно методом Рунге–Кутты 7–8 порядка с контролем точности на каждом шаге интегрирования. При расчете использовались параметры материала, типичные для монокристалла ЖИГ: M0 = 138 Гс, H0 = 2750 Э, c44 = = 7.64 ⋅ 1011 эрг ∙ см–3, B2 = 6.96 ⋅ 106 эрг ∙ см–3, d = = 68.65 ∙ 10–6см. Параметры затухания магнитной и упругой подсистем были выбраны следующими: α = 0.001, β = 106 с–1.

ОСОБЕННОСТИ ПРЕОБРАЗОВАНИЯ ЧАСТОТ

Рассмотрим случай генерации колебаний компонент вектора намагниченности и упругого смещения, когда частота акустического резонанса пленки находится в определенном соотношении с частотой переменного магнитного поля, которое можно описать соотношением n = ω0ar, где n – кратность, ωar – частота акустического резонанса, а ω0 – частота переменного поля. При этом, когда частота акустического резонанса ниже частоты переменного магнитного поля, происходит деление частоты, в обратном случае умножение частоты.

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

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

Как видно из рис. 1, максимальная амплитуда колебаний компонент вектора намагниченности на частоте переменного поля при круговой поляризации (сплошная линия) меньше амплитуды колебаний компонент вектора намагниченности при линейной поляризации (пунктирная линия) за исключением деления с кратностью менее 0.5. При этом различие амплитуд колебаний компонент вектора намагниченности на частоте акустического резонанса минимально. Это может быть обусловлено тем, что процессы преобразования частот при линейной и круговой поляризации происходят различным образом. В случае круговой поляризации (сплошная линия) важное значение имеет взаимное влияние процессов релаксации магнитных и упругих колебаний, а именно оптимальный момент наблюдения преобразования частоты находится между временами релаксации магнитной и упругой подсистем. В случае линейной поляризации (пунктирная линия) переменного поля, подобная зависимость отсутствует.

Рис. 1.

АЧХ колебаний компоненты mx вектора намагниченности на частоте акустического резонанса (I) и переменного поля (II) при линейной (пунктирная линия) и круговой (сплошная линия) поляризации переменного магнитного поля.

Различие в амплитудах колебаний компонент вектора намагниченности на частотах переменного поля (линии II) и частоте акустического резонанса (линии I) пленки составляет два порядка.

На амплитудно-частотной характеристике колебаний упругого смещения, представленной на рис. 2, видно, что, как и в случае с колебаниями вектора намагниченности, амплитуда колебаний упругого смещения больше в случае возбуждения линейно поляризованным магнитным полем (пунктирная линия). Однако, в случае кратности деления менее 0.5 и кратности умножения более 2.5, амплитуды колебаний упругого смещения при круговой поляризации переменного поля больше или сопоставимы с амплитудами колебаний упругого смещения при линейной поляризации. Отметим также, что при делении частоты максимумы амплитуд колебаний вектора намагниченности и упругого смещения на частоте переменного поля отмечаются при соотношении частот кратном 0.5. При промежуточных значениях отмечаются минимумы амплитуды, а на частоте акустического резонанса все происходит с точностью до наоборот, но менее выражено. Изменение амплитуд колебаний упругого смещения от соотношения частот переменного поля и акустического резонанса связано с распределением моды упругой волны в образце.

Рис. 2.

АЧХ колебаний компоненты ux упругого смещения на частоте акустического резонанса (I) и переменного поля (II) при линейной (пунктирная линия) и круговой (сплошная линия) поляризации переменного магнитного поля.

ЗАКЛЮЧЕНИЕ

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

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

Работа выполнена при поддержке Российского фонда фундаментальных исследований (гранты №№ 17-57-150001 и 17-02-01138).

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

  1. Chang C.L., Lomonosov A., Janusonis J. et al. // Phys. Rev. B. 2017. V. 95. Art. № 060409.

  2. Temnov V.V., Razdolski I., Pezeril T. et al. // J. Opt. 2016. V. 18. Art. № 093002.

  3. Asadullin F.F., Poleshikov S.M., Pleshev D.A. et al. // J. Siberian Fed. Univ. Math. Phys. 2017. № 10(1). P. 36.

  4. Pleshev D., Vlasov V., Asadullin F. et al. // EPJ Web Conf. 2018. V. 185. Art. № 02004.

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