All issues

В статье рассматривается квартичное семейство планарных векторных полей, соответствующее рациональной системе Холлинга, которая моделирует динамику популяций типа «хищник–жертва» в данной экологической или биомедицинской системе и которая обобщает классическую систему Лотки–Вольтерры. В простейших математических моделях изменение концентрации жертв в единицу времени в расчете на одного хищника, которое характеризуется так называемой функцией отклика, прямо пропорционально концентрации жертв, т. е. функция отклика в этих моделях линейная. Это означает, что в системе нет насыщения хищников, когда количество жертв достаточно велико. Однако было бы более реалистично рассматривать нелинейные и ограниченные функции отклика, и в литературе действительно используются различные виды таких функций для моделирования отклика хищников. После алгебраических преобразований рациональную систему Холлинга можно записать в виде квартичной динамической системы. Для исследования характера и расположения особых точек в фазовой плоскости этой системы используется разработанный нами метод, смысл которого состоит в том, чтобы получить простейшую (хорошо известную) систему путем обращения в нуль некоторых параметров (обычно параметров, поворачивающих поле) исходной системы, а затем последовательно вводить эти параметры, изучая динамику особых точек (как конечных, так и бесконечно удаленных) в фазовой плоскости. Используя полученную информацию об особых точках и применяя наш геометрический подход к качественному анализу, мы изучаем бифуркации предельных циклов квартичной системы. Чтобы контролировать все бифуркации предельных циклов, особенно бифуркации кратных предельных циклов, необходимо знать свойства и комбинировать действия всех параметров, поворачивающих векторное поле системы. Это может быть сделано с помощью принципа окончания Уинтнера–Перко, согласно которому максимальное однопараметрическое семейство кратных предельных циклов заканчивается либо в особой точке, которая, как правило, имеет ту же кратность (цикличность), либо на сепаратрисном цикле, который также, как правило, имеет ту же кратность (цикличность). Применяя этот принцип, мы доказываем, что квадричная система (и соответствующая рациональная система Холлинга) может иметь не более двух предельных циклов, окружающих одну особую точку.

In this paper, we consider a quartic family of planar vector fields corresponding to a rational Holling system which models the dynamics of the populations of predators and their prey in a given ecological or biomedical system and which is a variation on the classical Lotka–Volterra system. For the latter system, the change of the prey density per unit of time per predator called the response function is proportional to the prey density. This means that there is no saturation of the predator when the amount of available prey is large. However, it is more realistic to consider a nonlinear and bounded response function, and in fact different response functions have been used in the literature to model the predator response. After algebraic transformations, the rational Holling system can be written in the form of a quartic dynamical system. To investigate the character and distribution of the singular points in the phase plane of the quartic system, we use our method the sense of which is to obtain the simplest (well-known) system by vanishing some parameters (usually field rotation parameters) of the original system and then to input these parameters successively one by one studying the dynamics of the singular points (both finite and infinite) in the phase plane. Using the obtained information on singular points and applying our geometric approach to the qualitative analysis, we study the limit cycle bifurcations of the quartic system. To control all of the limit cycle bifurcations, especially, bifurcations of multiple limit cycles, it is necessary to know the properties and combine the effects of all of the rotation parameters. It can be done by means of the Wintner–Perko termination principle stating that the maximal one-parameter family of multiple limit cycles terminates either at a singular point which is typically of the same multiplicity (cyclicity) or on a separatrix cycle which is also typically of the same multiplicity (cyclicity). Applying this principle, we prove that the quartic system (and the corresponding rational Holling system) can have at most two limit cycles surrounding one singular point.

Мы проводим глобальный бифуркационный анализ квартичной модели типа «хищник–жертва». В частности, исследуя глобальные бифуркации особых точек и предельных циклов, мы доказываем, что соответствующая динамическая система имеет не более двух предельных циклов.

We complete the global bifurcation analysis of a quartic predator–prey model. In particular, studying global bifurcations of singular points and limit cycles, we prove that the corresponding dynamical system has at most two limit cycles.

В приближении однородной намагниченности построена математическая модель ячейки памяти MRAM c осью анизотропии, расположенной в плоскости запоминающего ферромагнитного слоя ячейки и ориентированной параллельно ее краю (продольная анизотропия). Модель базируется на уравнении Ландау–Лифшица–Гильберта с токовым членом в форме Слончевского–Берже. Выведена система обыкновенных дифференциальных уравнений в нормальном виде, описывающая динамику намагниченности в трехслойной вентильной структуре Co/Cu/Co в зависимости от величины тока инжекции и внешнего магнитного поля, параллельного оси анизотропии магнитных слоев. Показано, что при любых токах и полях система имеет два основных состояния равновесия, расположенных на оси, совпадающей с осью анизотропии. Проведен анализ устойчивости этих состояний равновесия. Выписаны уравнения для определения дополнительных состояний равновесия. Показано, что в зависимости от величины внешнего магнитного поля и тока инжекции система может иметь всего два, четыре и шесть симметричных относительно оси анизотропии положений равновесия. Построены бифуркационные диаграммы, характеризующие основные типы динамики вектора намагниченности свободного слоя. Проведена классификация фазовых портретов на единичной сфере в зависимости от управляющих параметров (тока и поля). Изучены особенности динамики вектора намагниченности в каждой из характерных областей бифуркационной диаграммы и численно построены траектории переключения. Для построения траекторий использовался метод Рунге–Кутты. Найдены параметры, при которых существуют неустойчивые и устойчивые предельные циклы. Установлено, что неустойчивые предельные циклы существуют вокруг основного устойчивого равновесия на оси, совпадающей с осью анизотропии, а устойчивые циклы — вокруг неустойчивых дополнительных равновесий. Граница области существования устойчивых предельных циклов рассчитана численно. Обнаружены новые типы динамики под влиянием внешнего магнитного поля и спин-поляризованного тока инжекции: случайное и неполное переключение намагниченности. Аналитически определены значения пороговых токов переключения в зависимости от внешнего магнитного поля. Численно выполнены оценки времени переключения в зависимости от величин управляющих параметров.

The mathematical model of the magnetic memory cell MRAM with the in-plane anisotropy axis parallel to the edge of a free ferromagnetic layer (longitudinal anisotropy) has been constructed using approximation of uniform magnetization. The model is based on the Landau–Lifshits–Gilbert equation with the injection-current term in the Sloncžewski–Berger form. The set of ordinary differential equations for magnetization dynamics in a three-layered Co/Cu/Cu valve under the control of external magnetic field and spin-polarized current has been derived in the normal coordinate form. It was shown that the set of equations has two main stationary points on the anisotropy axis at any values of field and current. The stationary analysis of them has been performed. The algebraic equations for determination of additional stationary points have been derived. It has been shown that, depending on the field and current magnitude, the set of equations can have altogether two, four, or six stationary points symmetric in pairs relatively the anisotropy axis. The bifurcation diagrams for all the points have been constructed. The classification of the corresponding phase portraits has been performed. The typical trajectories were calculated numerically using Runge–Kutta method. The regions, where stable and unstable limit cycles exist, have been determined. It was found that the unstable limit cycles exist around the main stable equilibrium point on the axis that coincides with the anisotropy one, whereas the stable cycles surround the unstable additional points of equilibrium. The area of their existence was determined numerically. The new types of dynamics, such as accidental switching and non-complete switching, have been found. The threshold values of switching current and field have been obtained analytically. The estimations of switching times have been performed numerically.

В данной статье, используя наш бифуркационно-геометрический подход, мы изучаем глобальную динамику и решаем проблему о максимальном числе и распределении предельных циклов (автоколебательных режимов, соответствующих состояниям динамического равновесия) в планарной полиномиальной механической системе типа Эйлера–Лагранжа–Льенара. Такие системы используются также для моделирования электротехнических, экологических, биомедицинских и других систем, что значительно облегчает исследование соответствующих реальных процессов и систем со сложной внутренней динамикой. Они используется, в частности, в механических системах с демпфированием и жесткостью. Существует ряд примеров технических систем, которые описываются с помощью квадратичного демпфирования в динамических моделях второго порядка. В робототехнике, например, квадратичное демпфирование появляется при управлении с прямой связью и в нелинейных устройствах, таких как приводы с переменным импедансом (сопротивлением). Приводы с переменным сопротивлением представляют особый интерес для совместной робототехники. Для исследования характера и расположения особых точек в фазовой плоскости полиномиальной системы Эйлера–Лагранжа–Льенара используется разработанный нами метод, смысл которого состоит в том, чтобы получить простейшую (хорошо известную) систему путем обращения в нуль некоторых параметров (обычно параметров, поворачивающих поле) исходной системы, а затем последовательно вводить эти параметры, изучая динамику особых точек в фазовой плоскости. Для исследования особых точек системы мы используем классические теоремы Пуанкаре об индексе, а также наш оригинальный геометрический подход, основанный на применении метода двух изоклин Еругина, что особенно эффективно при исследовании бесконечно удаленных особых точек. Используя полученную информацию об особых точках и применяя канонические системы с параметрами, поворачивающими векторное поле, а также используя геометрические свойства спиралей, заполняющих внутренние и внешние области предельных циклов, и применяя наш геометрический подход к качественному анализу, мы изучаем бифуркации предельных циклов рассматриваемой системы.

In this paper, using our bifurcation-geometric approach, we study global dynamics and solve the problem of the maximum number and distribution of limit cycles (self-oscillating regimes corresponding to states of dynamical equilibrium) in a planar polynomial mechanical system of the Euler–Lagrange–Liйnard type. Such systems are also used to model electrical, ecological, biomedical and other systems, which greatly facilitates the study of the corresponding real processes and systems with complex internal dynamics. They are used, in particular, in mechanical systems with damping and stiffness. There are a number of examples of technical systems that are described using quadratic damping in second-order dynamical models. In robotics, for example, quadratic damping appears in direct-coupled control and in nonlinear devices, such as variable impedance (resistance) actuators. Variable impedance actuators are of particular interest to collaborative robotics. To study the character and location of singular points in the phase plane of the Euler–Lagrange–Liйnard polynomial system, we use our method the meaning of which is to obtain the simplest (well-known) system by vanishing some parameters (usually, field rotation parameters) of the original system and then to enter sequentially these parameters studying the dynamics of singular points in the phase plane. To study the singular points of the system, we use the classical Poincarй index theorems, as well as our original geometric approach based on the application of the Erugin twoisocline method which is especially effective in the study of infinite singularities. Using the obtained information on the singular points and applying canonical systems with field rotation parameters, as well as using the geometric properties of the spirals filling the internal and external regions of the limit cycles and applying our geometric approach to qualitative analysis, we study limit cycle bifurcations of the system under consideration.

В приближении однородной намагниченности построена математическая модель трехслойной ячейки памяти MRAM c осью анизотропии, расположенной перпендикулярно запоминающему ферромагнитному слою ячейки (перпендикулярная анизотропия). Предполагается, что первоначально намагниченность свободного слоя ячейки ориентирована вдоль оси анизотропии и соответствует состоянию «нуль». Одновременное мгновенное включение спин-поляризованного тока и магнитного поля воздействует на намагниченность свободного слоя и может перевести ее в противоположное положение, соответствующее состоянию «единица». Математическое описание эффекта основано на классическом векторном уравнении Ландау–Лифшица с диссипативным членом в форме Гильберта. В нашей модели учтены взаимодействия намагниченности с внешним магнитным полем и эффективными полями анизотропии и размагничивания, а также с током инжекции в форме Слончевского–Берже. Выведена система обыкновенных дифференциальных уравнений, описывающая динамику намагниченности в трехслойной вентильной структуре Co/Cu/Co в зависимости от управляющих параметров: величины тока инжекции и внешнего магнитного поля, параллельного оси анизотропии магнитных слоев. Показано, что при любых токах и полях система имеет два основных состояния равновесия, расположенных на оси, совпадающей с осью анизотропии. Установлено, что в данной системе, в отличие от системы с продольной анизотропией, дополнительные состояния равновесия отсутствуют. Проведен анализ устойчивости основных состояний равновесия по первому приближению. Построены бифуркационные диаграммы, характеризующие типы динамики вектора намагниченности свободного слоя. Проведена классификация фазовых портретов на единичной сфере в зависимости от управляющих параметров (тока и поля). Изучены особенности динамики вектора намагниченности в каждой из характерных областей бифуркационной диаграммы и численно, методом Рунге–Кутты, построены траектории переключения. Найдены комбинации управляющих параметров, при которых переключение невозможно. Найдены области существования устойчивых и неустойчивых предельных циклов системы. Аналитически определены значения пороговых токов переключения в зависимости от внешнего магнитного поля. Проведено сравнение значений порогового тока в моделях с продольной и перпендикулярной анизотропией при нулевом магнитном поле и показано, что в модели с перпендикулярной анизотропией ток переключения почти на порядок ниже, чем в модели с продольной анизотропией.

The mathematical model of a three-layered Co/Cu/Co nanopillar for MRAM cell with one fixed and one free layer was investigated in the approximation of uniformly distributed magnetization. The anisotropy axis is perpendicular to the layers (so-called perpendicular anisotropy). Initially the magnetization of the free layer is oriented along the anisotropy axis in the position accepted to be “zero”. Simultaneous magnetic field and spinpolarized current engaging can reorient the magnetization to another position which in this context can be accepted as “one”. The mathematical description of the effect is based on the classical vector Landau–Lifshits equation with the dissipative term in the Gilbert form. In our model we took into account the interactions of the magnetization with an external magnetic field and such effective magnetic fields as an anisotropy and demagnetization ones. The influence of the spin-polarized injection current is taken into account in the form of Sloczewski–Berger term. The model was reduced to the set of three ordinary differential equations with the first integral. It was shown that at any current and field the dynamical system has two main equilibrium states on the axis coincident with anisotropy axis. It was ascertained that in contrast with the longitudinal-anisotropy model, in the model with perpendicular anisotropy there are no other equilibrium states. The stability analysis of the main equilibrium states was performed. The bifurcation diagrams characterizing the magnetization dynamics at different values of the control parameters were built. The classification of the phase portraits on the unit sphere was performed. The features of the dynamics at different values of the parameters were studied and the conditions of the magnetization reorientation were determined. The trajectories of magnetization switching were calculated numerically using the Runge–Kutta method. The parameter values at which limit cycles exist were determined. The threshold values for the switching current were found analytically. The threshold values for the structures with longitudinal and perpendicular anisotropy were compared. It was established that in the structure with the perpendicular anisotropy at zero field the switching current is an order lower than in the structure with the longitudinal one.

Рассматривается известное эволюционное уравнение математической физики, которое в современной математической литературе принято называть уравнением Курамото–Сивашинского. В данной работе это уравнение изучается в первоначальной редакции авторов работ, где оно было предложено, вместе с однородными краевыми условиями Неймана. Изучен вопрос о существовании и устойчивости локальных аттракторов, сформированных пространственно-неоднородными решениями изучаемой краевой задачи. Данный вопрос стал особенно актуален в последнее время в связи с моделированием процесса формирования наноструктур на поверхности полупроводников под воздействием потока ионов или лазерного излучения.

Изучен вопрос о существовании и устойчивости состояний равновесия второго рода двумя различными способами. В первом из них использован метод Галёркина. Второй подход основан на использовании строго обоснованных методов теории динамических систем с бесконечномерным фазовым пространством: метод интегральных многообразий, теория нормальных форм, асимптотические методы.

В работе в целом повторен подход из известной работы Д. Армбрустера, Д. Гукенхеймера, Ф.Холмса, где использован подход, основанный на применении метода Галёркина. Результаты такого анализа расширены и развиты. Использование возможностей современных компьютеров помогло существенно дополнить анализ этой задачи. В частности, найти все решения в четырех- и пятичленных аппроксимациях Галёркина, которые для изучаемой краевой задачи следует интерпретировать как состояния равновесия второго рода. Также дан анализ их устойчивости в смысле определения А. М. Ляпунова.

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

The well-known evolutionary equation of mathematical physics, which in modern mathematical literature is called the Kuramoto – Sivashinsky equation, is considered. In this paper, this equation is studied in the original edition of the authors, where it was proposed, together with the homogeneous Neumann boundary conditions.

The question of the existence and stability of local attractors formed by spatially inhomogeneous solutions of the boundary value problem under study has been studied. This issue has become particularly relevant recently in connection with the simulation of the formation of nanostructures on the surface of semiconductors under the influence of an ion flux or laser radiation. The question of the existence and stability of second-order equilibrium states has been studied in two different ways. In the first of these, the Galerkin method was used. The second approach is based on using strictly grounded methods of the theory of dynamic systems with infinite-dimensional phase space: the method of integral manifolds, the theory of normal forms, asymptotic methods.

In the work, in general, the approach from the well-known work of D.Armbruster, D.Guckenheimer, F.Holmes is repeated, where the approach based on the application of the Galerkin method is used. The results of this analysis are substantially supplemented and developed. Using the capabilities of modern computers has helped significantly complement the analysis of this task. In particular, to find all the solutions in the fourand five-term Galerkin approximations, which for the studied boundary-value problem should be interpreted as equilibrium states of the second kind. An analysis of their stability in the sense of A. M. Lyapunov’s definition is also given.

In this paper, we compare the results obtained using the Galerkin method with the results of a bifurcation analysis of a boundary value problem based on the use of qualitative analysis methods for infinite-dimensional dynamic systems. Comparison of two variants of results showed some limited possibilities of using the Galerkin method.

В работе изучена динамика двух искусственных генетических осцилляторов — репрессиляторов, — связанных диффузией аутоиндуктора. Выбрана модель генетической сети, в которой производство, диффузия и ген-мишень для аутоиндуктора обеспечивают расталкивающее взаимодействие между фазовыми точками. Исследовано появление периодических режимов, устойчивых неоднородных стационарных состояний в зависимости от главных бифуркационных параметров: силы связи и скорости синтеза мРНК. Показано, что добавление в генетическую схему аутоиндуктора приводит к исчезновению предельного цикла через бифуркацию бесконечного периода в изолированном осцилляторе, если скорость синтеза мРНК велика. Найден гистерезис между предельным циклом и стационарным состоянием, размер которого зависит от соотношения времен жизни мРНК и белков. Взаимодействие двух осцилляторов приводит к появлению устойчивого противофазного предельного цикла, который может переходить в хаотический режим через «тор-хаос» или путем каскада Фейгенбаума.

We have investigated dynamics of synthetic genetic oscillators — repressilators — coupled through autoinducer diffusion. The model of the system with phase-repulsive coupling structure is under consideration. We have examined emergence of periodic regimes, stable inhomogeneous steady states depending on the main systems’ parameters: coupling strength and maximal transcription rate. It has been shown that autoinducer production module added to the isolated repressilator cause the limit cycle to disappear through infinite period bifurcation for sufficiently large transcription rate. We have found hysteresis of limit cycle and stable steady state the size of which is determined by ratio between mRNA and protein lifetimes. Two coupled oscillators system demonstrates stable anti-phase oscillations which can become a chaotic regime through invariant torus emergence or via Feigenbaum scenario.