Laboratory of Mathematical Modeling of Nonlinear Processes 

Technical Center of the NAS of Ukraine


Technical Center

National Academy of Sciences of Ukraine

Laboratory of Mathematical Modeling of Nonlinear Processes

Pokrovsʹka St., 13, Kyiv, 04070, Ukraine


Recent Publications

Asymmetry induces critical desynchronization of power grids.    

Patrycja Jaros, Roman Levchenko, Tomasz Kapitaniak, Jurgen Kurths, and Yuri Maistrenko

Chaos 33, 011104 (2023);

Dynamical stability of the synchronous regime remains a challenging problem for secure functioning of power grids. Based on the symmetric alternative model [Nature Communication 11, 592 (2020)], we demonstrate that the grid stability can be destroyed by elementary violations (motifs) of the network architecture such as cutting a connection between any two nodes or removing a generator or a consumer. We describe the mechanism for the cascading failure in each of the damaging case and show that the desynchronization starts with the frequency deviation of the neighbouring grid elements followed by the cascading splitting of the others, distant elements and ending eventually in the bi-modal or a partially desynchronized state. Our findings reveal that symmetric topology underlines stability of the power grids, while local damaging can cause a fatal blackout 

Mixed-Mode Chimera States in Pendula Networks

P. Ebrahimzadeh, M. Schiek, and Y. Maistrenko

Chaos 32, 103118 (2022),

We report the emergence of peculiar chimera states in networks of identical pendula with global phase-lagged coupling. The states reported include both rotating and quiescent modes, i.e. with non-zero and zero average frequencies. This kind mixed-mode chimeras may be interpreted as images of bump states known in neuroscience in the context of modelling the working memory. We illustrate this striking phenomenon for a network of N = 100 coupled pendula, followed by a detailed description of the minimal non-trivial case of N = 3. Parameter regions for five characteristic types of the system behavior are identified consisting: two mixed-mode chimeras with one and two rotating pendula, classical weak chimera with all three pendula rotating, synchronous rotation and quiescent state. The network dynamics is multistable: up to four of the states can coexist in the system phase state as demonstrated through the basins of attraction. The analysis suggests that the robust mixed-mode chimera states can generically describe the complex dynamics of diverse pendula-like systems widespread in nature.

Chimera states for directed networks.    

Patrycja Jaros, Roman Levchenko, Tomasz Kapitaniak, and Yuri Maistrenko

Chaos 31, 103111 (2021);

We demonstrate that chimera behavior can be observed in ensembles of phase oscillators with unidirectional coupling. For a small network consisting of only three identical oscillators (cyclic triple), tiny chimera islands arise in the parameter space. They are surrounded by developed chaotic switching behavior caused by a collision of rotating waves propagating in opposite directions. For larger networks, as we show for a hundred oscillators (cyclic century), the islands merge into a single chimera continent, which incorporates the world of chimeras of different configurations. The phenomenon inherits from networks with intermediate ranges of the unidirectional coupling and it diminishes as the coupling range decreases.

Generalized splay states in phase oscillator networks.

Rico Berner, Serhiy Yanchuk, Yuri Maistrenko, and Eckehard Schöll

Chaos 31, 073128 (2021);


Networks of coupled phase oscillators play an important role in the analysis of emergent collective phenomena. In this article, we introduce generalized 𝑚-splay states constituting a special subclass of phase-locked states with vanishing 𝑚th order parameter. Such states typically manifest incoherent dynamics, and they often create high-dimensional families of solutions (splay manifolds). For a general class of phase oscillator networks, we provide explicit linear stability conditions for splay states and exemplify our results with the well-known Kuramoto–Sakaguchi model. Importantly, our stability conditions are expressed in terms of just a few observables such as the order parameter or the trace of the Jacobian. As a result, these conditions are simple and applicable to networks of arbitrary size. We generalize our findings to phase oscillators with inertia and adaptively coupled phase oscillator models.

Scroll ring chimera states in oscillatory networks.

Volodymyr Maistrenko, Oleksandr Sudakov, Ievgen Sliusar

We report the appearance of a scroll ring and scroll toroid chimera states from the proposed initial conditions for the Kuramoto model of coupled phase oscillators in the 3D grid topology with inertia. The proposed initial conditions provide an opportunity to obtain as single as well as multiple scroll ring and toroid chimeras with different major and minor diameters. We analyze their properties and demonstrate, in particular, the patterns of coherent, partially coherent, and incoherent scroll ring chimera states with different structures of filaments and chaotic oscillators. Those patterns can coexist with solitary states and solitary patterns in the oscillatory networks.
Illustrative video (mp4): VideoFig3a   VideoFig3b   VideoFig7 

Chimera complexity.

Serhiy Brezetsky, Patrycja Jaros, Roman Levchenko, Tomasz Kapitaniak, and Yuri Maistrenko

We show an amazing complexity of the chimeras in small networks of coupled phase oscillators with inertia. The network behavior is characterized by heteroclinic switching between multiple saddle chimera states and riddling basins of attractions, causing an extreme sensitivity to initial conditions and parameters. Additional uncertainty is induced by the presumable coexistence of stable phase-locked states or other stable chimeras as the switching trajectories can eventually tend to them. The system dynamics becomes hardly predictable, while its complexity represents a challenge in the network sciences. 

Solitary states in the mean-field limit.

N. Kruk, Y. Maistrenko, and H. Koeppl

Chaos 30, 111104 (2020);

We study active matter systems where the orientational dynamics of underlying self-propelled particles obey second order equations. By primarily concentrating on a spatially homogeneous setup for particle distribution, our analysis combines theories of active matter and oscillatory networks. For such systems, we analyze the appearance of solitary states via a homoclinic bifurcation as a mechanism of the frequency clustering. By introducing noise, we establish a stochastic version of solitary states and derive the mean-field limit described by a partial differential equation for a one-particle probability density function, which one might call the continuum Kuramoto model with inertia and noise.By studying this limit, we establish second order phase transitions between polar order and disorder. The combination of both analytical and numerical approaches in our study demonstrates an excellent qualitative agreement between mean-field and finite size models.

Minimal Chimera States in Phase-Lag CoupledMechanical Oscillators.

P. Ebrahimzadeh, M. Schiek, P. Jaros, T. Kapitaniak, S. van Waasen and Y. Maistrenko

Eur. Phys. J. Special Topics 229, 2205 (2020).

We obtain experimental chimera states in the minimal net-work of three identical mechanical oscillators (metronomes), by intro-ducing phase-lagged all-to-all coupling. For this, we have developed areal-time model-in-the-loop coupling mechanism that allows for flexibleand online change of coupling topology, strength and phase-lag. Thechimera states manifest themselves as a mismatch of average frequencybetween two synchronous and one desynchronized oscillator. We findthis kind of striking chimeric behavior is robust in a wide parameterregion. At other parameters, however, chimera state can lose stabilityand the system behavior manifests itself as a heteroclinic switching be-tween three saddle-type chimeras. Our experimental observations arein a qualitative agreement with the model simulation.

Spiral wave chimeras for coupled oscillators with inertia.

Volodymyr Maistrenko, Oleksandr Sudakov, and Yuri Maistrenko

Eur. Phys. J. Spec. Top. 229, 2327–2340 (2020).

We report the appearance and the metamorphoses of spiral wave chimera states in coupled phase oscillators with inertia. First, when the coupling strength is small enough, the system behavior resembles classical two-dimensional (2D) Kuramoto-Shima spiral chimeras with bell-shape frequency characteristic of the incoherent cores. As the coupling increases, the cores acquire concentric regions of constant time-averaged frequencies, the chimera becomes quasiperiodic. Eventually, with a subsequent increase in the coupling strength, only one such region is left, i.e., the whole core becomes frequency-coherent. An essential modification of the system behavior occurs, when the parameter point enters the so-called 'solitary' region. Then, isolated oscillators are normally present on the spiral core background of the chimera states. These solitary oscillators do not participate in the common spiraling around the cores; instead, they start to oscillate with different time-averaged frequencies (Poincar\'e winding numbers). The number and the disposition of solitary oscillators can be any, given by the initial conditions. At a further increase in the coupling, the spiraling disappears, and the system behavior passes to a sort of spatiotemporal chaos.

Chimera and solitary states in 3D oscillator networks with inertia.

Volodymyr Maistrenko, Oleksandr Sudakov, and Oleksiy Osiv

Chaos 30, 063113 (2020); 

We report the diversity of scroll wave chimeras in the three-dimensional (3D) Kuramoto model with inertia for N3 identical phase oscillators placed in a unit 3D cube with periodic boundary conditions. In the considered model with inertia, we have found patterns which do not exist in a pure system without inertia. In particular, a scroll ring chimera is obtained from random initial conditions. In contrast to this system without inertia, where all chimera states have incoherent inner parts, these states can have partially coherent or fully coherent inner parts as exemplified by a scroll ring chimera. Solitary states exist in the considered model as separate states or can coexist with scroll wave chimeras in the oscillatory space. We also propose a method of construction of 3D images using solitary states as solutions of the 3D Kuramoto model with inertia.

Network-induced multistability through lossy coupling and exotic solitary states.

Frank Hellmann, Paul Schultz, Patrycja Jaros, Roman Levchenko, Tomasz Kapitaniak, Jürgen Kurths and Yuri Maistrenko

Nature Communications 11:592 (2020);

The stability of synchronised networked systems is a multi-faceted challenge for manynatural and technologicalfields, from cardiac and neuronal tissue pacemakers to power grids.For these, the ongoing transition to distributed renewable energy sources leads to a pro-liferation of dynamical actors. The desynchronisation of a few or even one of those wouldlikely result in a substantial blackout. Thus the dynamical stability of the synchronous statehas become a leading topic in power grid research. Here we uncover that, when taking intoaccount physical losses in the network, the back-reaction of the network induces new exoticsolitary states in the individual actors and the stability characteristics of the synchronousstate are dramatically altered. These effects will have to be explicitly taken into account in thedesign of future power grids. We expect the results presented here to transfer to othersystems of coupled heterogeneous Newtonian oscillators.