Радиотехника и электроника, 2023, T. 68, № 9, стр. 893-896
Влияние нелинейности на особую точку в системе связанных осцилляторов Дуффинга
О. С. Темная a, *, А. Р. Сафин a, b, О. В. Кравченко a, c, С. А. Никитов a, d, e
a Институт радиотехники и электроники им. В.А. Котельникова РАН
125009 Москва, ул. Моховая, 11, стр. 7, Российская Федерация
b Национальный исследовательский университет “МЭИ”
111250 Москва, ул. Красноказарменная, 14, стр. 1, Российская Федерация
c Вычислительный центр им. А.А. Дородницына
Федерального исследовательского центра “Информатика и управление” РАН
119333 Москва, ул. Вавилова, 40, Российская Федерация
d Московский физико-технический институт (национальный исследовательский университет)
141700 Долгопрудный, Московской области, Институтский пер., 9, Российская Федерация
e Cаратовский национальный исследовательский государственный университет им. Н.Г. Чернышевского
410012 Саратов, ул. Астраханская, 83, Российская Федерация
* E-mail: ostemnaya@gmail.com
Поступила в редакцию 10.06.2023
После доработки 15.06.2023
Принята к публикации 25.06.2023
- EDN: SJESLH
- DOI: 10.31857/S0033849423090231
Аннотация
Исследовано влияние нелинейности на смещение особой точки в системе двух связанных осцилляторов Дуффинга при изменении коэффициентов связи и вносимых потерь. Показано, что смещение особой точки при изменении коэффициента нелинейности сопровождается уменьшением амплитуды возбуждаемых колебаний и сдвигом резонансной частоты. Численно найдены пороговые значения коэффициентов нелинейности, связи и вносимых потерь, при которых возникает особая точка. Показано, что увеличение коэффициента нелинейности приводит к уменьшению порогового значения вносимых потерь, необходимых для образования особой точки.
ВВЕДЕНИЕ
Особыми точками (ОТ) называют такие точки пространства параметров, в которых собственные векторы и собственные значения совпадают [1]. Такое вырождение характерно для неэрмитовых квантовых систем с симметрией “четность–время” (PT-симметрией). Вещественные спектры таких систем при переходе через особую точку становятся комплексными, что сопровождается спонтанным нарушением симметрии. В РТ-симметричных системах добиться появления ОТ можно за счет изменения величины вносимых потерь и силы связи между частями системы [2]. В ОТ могут возникать различные интересные с практической точки зрения эффекты, такие как резкое увеличение чувствительности, увеличенный диапазон генерации частот, лучшее шумоподавление [2, 3]. Таким образом, важной задачей является поиск условий возникновения ОТ и способов управления их положением в пространстве параметров системы.
Особые точки были обнаружены в различных физических средах и системах, включая оптику [2, 4, 5], акустику [6, 7], электронику [3], магнонику и спинтронику [8–11]. В большинстве приведенных выше примеров для наблюдения ОТ используются не квантовые, а классические системы, а именно колебательные и волновые. С точки зрения теории колебаний и волн две нормальные моды двух взаимосвязанных систем за счет вносимых положительных и отрицательных потерь вырождаются в одну собственную в ОТ. Для этого системы должны быть идентичными, вносимые потери в одной из них положительными (затухание усиливается), а в другой отрицательными (затухание компенсируется). При таком внесении потерь резонансные свойства колебательных систем в ОТ не ухудшаются и определяются собственным затуханием. Величина же вносимых потерь, необходимая для появления ОТ, прямо пропорционально зависит от величины связи между системами. В большинстве работ, посвященных ОТ в колебательных системах, анализ проводится в линейном случае, тогда как нелинейность предполагается малой (см., например, [3, 8]). Вместе с тем существенная нелинейность может быть дополнительной степенью свободы для возможности управления положением особой точки.
Цель данной работы – исследовать влияние консервативной кубической нелинейности на ОТ в системе связанных осцилляторов Дуффинга. Выбор данного типа осцилляторов обусловлен их широким применением для описания различных нелинейных колебательных систем – от электрических генераторов до спинтронных наноосцилляторов [12–14].
1. АНАЛИЗ МАТЕМАТИЧЕСКОЙ МОДЕЛИ
Рассмотрим систему неавтономных дифференциальных уравнений относительно координат x(t), y(t) двух взаимосвязанных нелинейных осцилляторов с коэффициентом связи μ и коэффициентом нелинейности β следующего вида:
(1)
$\left\{ {\begin{array}{*{20}{c}} {\ddot {x} + \left( {{{\Delta }_{0}} - \Delta } \right)\dot {x} + \omega _{0}^{2}\left( {1 + \beta {{x}^{2}}} \right)x + \mu y = f\sin \omega t,} \\ {\ddot {y} + \left( {{{\Delta }_{0}} + \Delta } \right)\dot {y} + \omega _{0}^{2}\left( {1 + \beta {{y}^{2}}} \right)y + \mu x = 0,} \end{array}} \right.$Проведем численное интегрирование системы уравнений (1) стандартными методами анализа
дифференциальных уравнений (метод Рунге–Кутты четвертого порядка) при следующих параметрах:
${{\omega }_{0}} = 1,$ ${{\Delta }_{0}} = 0.1,$ $f = 0.1,$ $\mu = 0.7,$ $\Delta ,\beta ,\omega = \operatorname{var} .$ При в системе двух связанных осцилляторов наблюдаются две нормальные частоты ${{\omega }_{{1,2}}}$ на амплитудно-частотной характеристике (рис. 1а), которые при увеличении коэффициента нелинейности $\beta $ > 0 также увеличиваются (при $\Delta = 0$), причем помимо смещения нормальных частот при некотором значении $\beta $ возникает гистерезисный эффект с образованием неустойчивых ветвей [12, 13]. При этом для фиксированной частоты внешнего воздействия амплитуда нижней моды системы
уменьшается при увеличении частоты вынужденных колебаний.
Рис. 1.
Семейство амплитудно-частотных характеристик двух связанных осцилляторов Дуффинга (а) и зависимости действительных частей нормальных частот от величины вносимых потерь (б) для трех значений коэффициента нелинейности: β = 0 (1), 1 (2) и 3 (3); штрихпунктирной линией обозначены неустойчивые области резонансной кривой.

Расстояние между резонансными пиками (рис. 1а) напрямую зависит от константы связи между осцилляторами. Как было сказано выше, ОТ возникает из-за вырождения нормальных частот ${{\omega }_{{1,2}}}$ в собственную ${{\omega }_{0}}$ за счет увеличения коэффициента вносимых потерь $\Delta $. В линейной системе для любой величины константы связи μ существует единственное значение вносимых потерь ${{\Delta }_{{{\text{ОТ}}}}}$, при котором возникает ОТ (рис. 1б). При этом имеется такое пороговое значение ${{\beta }_{{\text{п}}}}$ (для приведенных выше параметров ${{\beta }_{{\text{п}}}} \approx 0.79$), при котором смещение особой точки резко меняется, а при $\beta < {{\beta }_{{\text{п}}}}$ частота появления ОТ практически не меняется при увеличении $\Delta $. Критическое значение ${{\Delta }_{{{\text{ОТ}}}}}$ также остается неизменным, однако при $\beta > {{\beta }_{{\text{п}}}}$ частота ОТ ${{{{\omega }}}_{{{\text{ОТ}}}}}$ увеличивается, а ${{\Delta }_{{{\text{ОТ}}}}}$ уменьшается (см. рис. 1б). Таким образом, наличие кубической нелинейности приводит к появлению ОТ при меньших значениях вносимых потерь. Увеличение параметра нелинейности $\beta $ приводит к резкому уменьшению ${{\Delta }_{{{\text{ОТ}}}}},$ при котором появляется ОТ (рис. 2а). При увеличении порогового значения параметра нелинейности резонансные частоты смещаются, вследствие чего требуется меньшее вносимое затухание (рис. 2б).
Рис. 2.
Зависимости вносимых потерь ${{\Delta }_{{{\text{ОТ}}}}}$ (а) и частоты ${{\omega }_{{\,{\text{ОТ}}}}}$ (б), при которых образуется особая точка, от величины коэффициента нелинейности β; штрихпунктиром обозначена линия критического параметра нелинейности ${{\beta }_{{\text{п}}}} \approx 0.79,$ при котором происходит смещение вносимых потерь и резонансной частоты.

В заключение проведем численное моделирование влияния зависимости вносимых потерь $\Delta $ от константы связи $\mu $ и параметра нелинейности $\beta ,$ при которых появляется особая точка. В линейном случае при $\beta = 0$ зависимость ОТ от $\mu $ линейна, а именно $\mu = \Delta .$ Как указывалось ранее, увеличение коэффициента нелинейности β приводит к появлению ОТ при меньших значениях Δ, что также проявляется при изменении параметра связи μ. При этом зависимость μ(Δ) имеет следующий вид:
где $A > 1$ – крутизна этой зависимости (рис. 3). Следовательно, при большей связи между осцилляторами наличие нелинейности приводит к уменьшению вносимых потерь, при котором появляется ОТ.
ЗАКЛЮЧЕНИЕ
Таким образом, исследовано влияние нелинейности на положение особой точки на плоскости параметров в системе двух консервативно связанных осцилляторов Дуффинга при изменении коэффициентов связи и вносимых потерь. Показано, что смещение положения особой точки при изменении коэффициента нелинейности сопровождается уменьшением амплитуды возбуждаемых колебаний и сдвигом резонансной частоты. Численно найдены пороговые значения коэффициентов нелинейности, связи и вносимых потерь, при которых возникает особая точка (см. рис. 3). Показано, что увеличение коэффициента нелинейности приводит к уменьшению порогового значения вносимых потерь, необходимых для образования особой точки. Предлагаемый анализ может быть проведен для осцилляторов различной физической природы, разных типов связи между ними, а также большего числа осцилляторов с появлением особых точек высших порядков. Соответствующие задачи являются предметом отдельного детального исследования и в данной работе не рассматриваются.
Список литературы
Kato T. A Short Introduction to Perturbation Theory for Linear Operators. N.Y.: Springer., 2011. https://doi.org/10.1007/978-1-4612-5700-4
Wiersig J. // Photon. Res. 2020. V. 8. № 9. P. 1457. https://doi.org/10.1364/PRJ.396115
Weidong C., Wang C., Chen W. et al. // Nat. Nanotech. 2022. V. 17. Article No. 262268. https://doi.org/10.1038/s41565-021-01038-4
Зябловский А.А., Виноградов А.П., Пухов А.А. и др. // Успехи физ. наук. 2014. Т. 184. № 11. С. 1177. https://doi.org/10.3367/UFNr.0184.201411b.1177
Rüter C., Makris K., El-Ganainy R. // Nat. Phys. 2010. V. 6. Article No. 192195. https://doi.org/10.1038/nphys1515
Вилков Е.A., Бышевский-Конопко О.А., Темная О.С. и др. // Письма в ЖТФ. 2022. Т. 48. № 24. С. 38. https://doi.org/10.21883/PJTF.2022.24.54023.19291
Zhu X., Ramezani H., Shi C. et al. // Phys. Rev. X 2014. V. 4. Article No. 031042. https://doi.org/10.1103/PhysRevX.4.031042
Wittrock S., Perna S., Lebrun R. et al. // arXiv: 2108.04804.
Liu H., Sun D., Zhang C. et al. // Sci. Adv. 2019. V. 5. № 11. Article No. aax9144. https://doi.org/10.1126/sciadv.aax9144
Temnaya O.S., Safin A.R., Kalyabin D.V., Nikitov S.A. // Phys. Rev. Appl. 2022. V. 18. Article No. 014003. https://doi.org/10.1103/PhysRevApplied.18.014003
Sadovnikov A.V., Zyablovsky A.A., Dorofeenko A.V., Nikitov S.A. // Phys. Rev. Appl. 2022. V. 18. Article No. 024073. https://doi.org/10.1103/PhysRevApplied.18.024073
Rajasekar S., Sanjuan M. Nonlinear Resonances. Cham: Springer, 2015.
Рабинович И.М., Трубецков Д.И. Введение в теорию колебаний и волн. Ижевск: НИЦ РХД, 2000.
Moon K.-W., Chun B.S., Kim W. et al. // Sci. Reports. 2014. V. 4. Article No. 6170.
Дополнительные материалы отсутствуют.
Инструменты
Радиотехника и электроника