Original Research ARTICLE
Cloning of Chimera States in a Large Short-term Coupled Multiplex Network of Relaxation Oscillators
- Nonlinear dynamics department, Institute of Applied Physics of the Russian Academy of Sciences, Nizhny Novgorod, Russia
A new phenomenon of the chimera states cloning in a large two-layer multiplex network with short-term couplings has been discovered and studied. For certain values of strength and time of multiplex interaction, in the initially disordered layer, a state of chimera is formed with the same characteristics (the same average frequency and amplitude distributions in coherent and incoherent parts, as well as an identical phase distribution in coherent part), as in the chimera which was set in the other layer. The mechanism of the chimera states cloning is examined. It is shown that the cloning is not related with synchronization, but arises from the competition of oscillations in pairs of oscillators from different layers.
Study of the formation of chimera states, i.e., peculiar types of hybrid states consisting of oscillators with coherent and incoherent behavior is one of the hot problems of the modern non-linear dynamics. To date, the chimera states have been discovered not only in a variety of theoretical papers [1–18], but also in experimental systems of various natures, for example, mechanical [19–22], optical [23, 24], chemical [25–29], and radiotechnical ones [30–33]. Similar states have also been registered in the neural activity of animal brain networks [34, 35]. At present, great attention is paid to studying of interaction of chimera states. The effects of generalized synchronization of chimera states , synchronization of chimera states in ensembles with asymmetrical connections , synchronization of chimera states in multiplex networks with delays , synchronization of chimera states in a two-layer multiplex network with adaptive connections in each layer , synchronization of chimera states in modular networks [40, 41], interaction of chimera states with fully coherent or fully incoherent states , etc. were explored. Note that in all these works the interaction of chimera states led to the formation of new chimera states (in some cases with synchronous coherent parts) that are different from the pre-existing chimera state. Recently we presented an example  of a two-layer multiplex network with seven oscillators in each layer where due to short-term interaction, one more chimera is emerged which is identical to the initial one (excluding the phase distribution of the incoherent part). We called this effect the chimera states cloning. In this article, we generalize the results of Dmitrichev et al.  for the case of a multiplex network with an arbitrary dimension of the layer and give a theoretical justification for the cloning effect based on the study of the fast-slow dynamics of the model using the methods of Geometric Singular Perturbation Theory (GPST) [44, 45].
2. Model of Multiplex Network
We consider a two-layer multiplex network with the topology illustrated in Figure 1A. Each layer of the network is a ring of locally and linearly coupled relaxational oscillators with phase portrait shown in Figure 1B. The dynamics of the network is described by the following system:
where f(u) = −u(u2 − a2)(u2 − b2)(u2 − c2); the parameters controlling the dynamics of the layers are for definiteness fixed as a = 0.32, b = 0.79, c = 1.166, ε = 0.001, and dr = 0.006; d > 0 and Tc > 0 are the parameters controlling the strength and the time of inter-layer (multiplex) interaction.
Figure 1. (A) Topology of multiplex network. Red and blue dots depict correspondingly oscillators of first and second layer. Solid and doted lines depict correspondingly intra-layer and (multiplex) inter-layer couplings. (B) Qulitative phase portrait of a single oscillator. There exist two stable limit cycles (bold black lines) separated by an unstable limit cycle (bold blue line) and an unstable equilibrium state (bold blue dot) located at the origin.
If the oscillators do not interact with each other, i.e., dr = 0; dm(t) ≡ 0, then the dynamics of each oscillator is described by second-order equation. Two stable limit cycles with “low” and “high” amplitudes exist on the (u, v) phase plane (see Figure 1B). The regions of attraction of these cycles are separated by an unstable limit cycle. An unstable equilibrium state is located at the coordinate origin (u = 0, v = 0). Thus, an oscillator can be either in the regime of low-amplitude oscillations with a dimensionless frequency of 0.0039 or, in the case of initial conditions “outside” the unstable cycle, in the regime of high-amplitude oscillations with a dimensionless frequency of 0.0021.
In our previous paper  we showed that various chimera states exist in a separate layer of system (1) at the chosen parameters and obtained chimera states at the experimental system consisting of seven bistable self-exciting oscillators with linear couplings. The example of chimera state is presented in Figure 2A by black crosses (distribution of instant phases φ, average amplitudes < A >, and average frequencies < ω >). Coherent part of the chimera state is formed by first 99 (j = 1 − 99) oscillators with high-amplitude oscillations, and incoherent part is formed by oscillators with j = 101 − 199 that demonstrate alternately the low- and high-amplitude oscillations. The oscillation frequencies and phases of oscillators were calculated as follows. For definiteness, we consider the jth oscillator in the ith layer. Let be a sequence of time instants at which the output voltage of the oscillator increases and intersects the straight line ; i.e.,
Then, the oscillation phase of the jth oscillator at the time t is given by the expression
and is the instantaneous oscillation frequency. This definition of the phase is meaningful only if are constant or quite close to each other. In the former case, the oscillator undergoes regular oscillations and the phase always varies at the same rate. In the latter case, oscillations can be, in particular, irregular and the phase is a piecewise linear function of the time. Furthermore, the phase introduced in Equation (2) describes the dynamics of only an individual oscillator decoupled from the remaining system. For this reason, to describe the dynamics of the system as a whole, it is more convenient to use the parameter describing the oscillation phase of the jth oscillator with respect to oscillations of the reference kth oscillator. If is independent of time, this means the phase matching of oscillations of the kth and jth oscillators. In the general case of interaction between oscillators, instantaneous frequencies and amplitudes are not constant. For this reason, here and below, frequencies and amplitudes are calculated with averaging over a quite long time series
Figure 2. Cloning of chimera state in system (1). Distribution of instant phases φ, average amplitudes < A > and average frequencies < ω > (from up to bottom) at the initial instant of time (A); after the interaction (B). Indexes “1” (black cross) and “2” (red circle) correspond respectively to the first and second layer states. Parameters values: a = 0.32, b = 0.79, c = 1.166, ε = 0.001, dr = 0.006, d = 0.06, Tc = 1, 000, and N = 200.
Notice that in addition to the coherent and incoherent parts, the chimera state also contains two isolated oscillators at j = 100 and j = 200. Such solitary states exist due to the bistability of the network oscillators. The explanation of this phenomenon was given in Nekorkin et al. . The bistability leads to formation of amplitude clusters with high- and low-amplitude oscillations and so to an amplitude gap (see the amplitude distribution in Figure 2A). The amplitude gap, in turn, causes a frequency gap. Since the magnitudes of gaps are quite large, diffusive coupling between the oscillators leads to emergence of solitary oscillators whose dynamics “smooths out” the dynamics of clusters (amplitude and frequency) having significantly different characteristics.
3. Cloning of Chimera States
Let the chimera state exist in the first layer at the initial instant of time when there is no interaction between the layers (black crosses in Figure 2A). And the initial conditions in the second layer correspond to the low-amplitude oscillations with random initial phases (red circles in Figure 2A). Now if we switch on the interaction between the layers, and then after some time switch the interaction off, then we can obtain a clone of the initial chimera state in the second layer. For certain values of strength and time of multiplex interaction, in the second layer, a chimera state is formed with the same average frequency and amplitude distributions in coherent and incoherent parts, as well as an identical phase distribution in coherent part, as in the chimera which was set in first layer. The example of such cloning for interaction strength d = 0.06 and interaction time Tc = 1, 000 is shown in Figure 2. Notice that, by definition, instantaneous phases of the incoherent parts in the chimera state should be random, and integral characteristics, e.g., average frequencies and amplitudes, are of fundamental importance for the incoherent part. For this reason, we believe that the coincidence of the phases of the incoherent parts is not necessary for the cloning of chimera states. Note that the cloning of chimera states occurs with certain initial conditions The role of initial conditions in more detail is discussed in section 4.
Next we show that the cloning effect is structurally stable. To do this, we introduce a reference layer with the index “0." Let us set in the reference layer a chimera state with the same characteristics as the original chimera that we set in the first layer. Then after the interaction of the first and the second layers we compare states formed in those layers with one in the reference layer using the following characteristics:
where < ωi,j > and < Ai,j > are averaged frequencies and amplitudes of the oscillators with the number j of the ith layer; accordingly, < ω0,j > and < A0,j > are those in the reference layer.
The results of such calculations for Tc = 1, 000 are presented in Figure 3. They indicate that there is an interval of multiplex coupling strength 0.37 ≤ d ≤ 0.96 where the states in all layers have the same averaged characteristics as the reference chimera state. Thus, the cloning occurs in the large enough interval of strength and so the effect is structurally stable.
Figure 3. Dependence of the maximum errors between the average frequencies (A) and amplitudes (B) in the reference chimera state and the states occurring in first (dashed line) and second (solid line) layer. Parameters values: a = 0.32, b = 0.79, c = 1.166, ε = 0.001, dr = 0.006, and Tc = 1, 000.
4. The Cloning Mechanism
We showed above that cloning of chimera states takes place when strength of coupling between the elements of the same layer is much smaller than that between the elements of different layers (dr < < d, see Figure 3). So in the first approximation we can assume that the key role in cloning is played by the dynamics of (multiplex) pairs of elements taken from different layers. Next we consider the dynamics of a pair in more detail. It is described by the following system of equations:
where , , , .
Notice also that to realize the cloning, initial conditions in non-interacting layers must be formed in a special way. In particular, a chimera state is set in the first layer with coherent part formed by the oscillators, demonstrating high-amplitude oscillations, and incoherent part formed by the oscillators demonstrating alternately low- and high-amplitude oscillations. In the second layer all oscillators are set in the regime of low-amplitude oscillations whose phases are randomly distributed. Moreover, after interaction, the elements of the second layer should switch to the regimes the corresponding elements of the first layer were in initial moment. Thus, we need to consider the evolution of a pair only for two types of initial conditions:
(I.C.)1 An oscillator of the first layer is in the regime of high-amplitude oscillations, while that of the second layer is in the regime of low-amplitude oscillations;
(I.C.)2 The oscillators of both layers are in the regime of low-amplitude oscillations.
4.1. Dynamics of a Pair of Constantly Coupled Oscillators
First, we study dynamics of a pair for the case when interaction between its oscillators is not limited by time. Then we obtain the following system of equations:
Since 0 < ε ≪ 1, the system (4) belongs to the class of fast-slow systems. Such systems are characterized by the presence of two timescales (or speeds), namely, fast and slow ones. In the result, the trajectories of the systems have epochs of a slow and a fast movements. In our system u1 and u2 are fast variables, while v1 and v2 are slow variables. Next to study the dynamics of the system we use GPST theory. According to the GPST, the partition of phase space ℝ4 of system (4) into trajectories can be established by studying two subsystems. As ε → 0, the trajectories of system (4) converge during fast epochs to the trajectories of the fast subsystem (or layer equations)
where t = ετ. During slow epochs the trajectories of (4) converge to the trajectories of the reduced system (or the slow flow)
The goal of GPST is to use the fast and slow subsystems (5) and (6) to understand the dynamics of the full system (4) for 0 < ε ≪ 1.
4.1.1. Dynamics of the Fast Subsystem
From system (5) one can see that v1 = const and v2 = const, so they play the role of additional parameters (denote them by and correspondingly). Thus, system (5) can be rewritten in the following gradient form
the trajectories of system (7) [and so system (5)], except for equilibrium states, relax to one of the stable equilibrium states. Moreover, since the system is gradient their trajectories relax to the equilibrium states by the fastest way. The number and type of equilibrium states depend on the parameters and may change due to saddle-node bifurcations. For example, for there are 49 equilibrium states, among which there are 16 stable and 9 unstable nodes and 24 saddles. The qualitative phase plane of system (7) in this case is shown in Figure 4.
Figure 4. Qualitative phase portrait of fast subsystem (7) for , a = 0.32, b = 0.79, and c = 1.166. The dots mark unstable nodes (blue), stable nodes (yellow), and saddles (black). The lines mark the separatrices of the saddle equilibrium states.
4.1.2. Dynamics of the Reduced Subsystem
The first two algebraic equations in system (6) define a critical manifold . Let us rewrite (6) in terms of the fast variables u1 and u2 to obtain the slow flow on S. For this, we differentiate algebraic Equation (6) with respect to t and combine the result with the equations for and :
System (8) is a system of linear inhomogeneous algebraic equations for derivatives and . Determinant of its coefficient matrix
If Δ ≠ 0, then system (8) has the only solution
Note that in the four-dimensional phase space ℝ4 of system (4), algebraic equation Δ = 0 and two algebraic equations of system (6) define set . Any its element, , defines a curve, where the fast and slow trajectories are stitched. According to GPST, each equilibrium state of fast subsystem (5) [and so (7)] in the phase space ℝ4 of system (4) corresponds to a submanifold, whose stability with respect to the trajectories of the fast subsystem coincides with the stability of the corresponding equilibrium state. Since the coordinates of the equilibrium states of system (5) depend on two parameters and , any such equilibrium state corresponds to a two-dimensional submanifold in the phase space of (4). The boundary of the submanifold is given by one of the curves . Hence the curves decompose the critical manifold S into some number of submanifolds of different stability
where denotes one of the stable submanifolds, are the unstable submanifolds, and are the saddle submanifolds of the slow flow. Note that Δ > 0 on submanifolds and , Δ < 0 on submanifolds and Δ = 0 on submanifolds . The number of these submanifolds depends on the parameter d. For example, Figure 5A for d = 0.06 on the phase plane (u1, u2) depicts submanifolds , corresponding to stable nodes of the fast subsystem and curves .
Figure 5. (A) Mutual location of the stable slow submanifolds () of system (5) on the phase plane (u1, u2). The boundaries of stable submanifolds () defined by the equations Δ = 0. (B) Trajectories behavior on slow submanifold . (C) Dependence of motion time over submanifold to its boundary on the initial conditions. Parameter values: a = 0.32, b = 0.79, c = 1.166, and d = 0.06.
Since the system (4) on submanifolds has no equilibrium states and limit cycles, any trajectory starting on the submanifolds eventually leave them. For example, Figure 5B shows the behavior of trajectories on starting on entering part of its boundary (blue line) and going till the exit part of the boundary (red line). Note that, the time the trajectories stay on varies depending on the initial conditions. The dependence for trajectories on is shown in Figure 5C. Note that the dependence is a monotonically increasing function asymptotically tending to the value T ≈ 220.
4.1.3. Dynamics of System (4) for Initial Conditions (I.C.)1
Let us study the dynamics of system (4) for initial conditions (I.C.)1. Note, the dynamics of system (4) is formed by the alternating dynamics of fast and slow epochs, which results in a “stitched” trajectory.
220.127.116.11. Slow epoch of motion
The initial conditions (I.C.)1 corresponds to one of the stable submanifolds of slow motions (see Figure 5A). First let the initial conditions belong to . The motions on are defined by system (9). Since Δ > 0 on , we can de-singularize the slow flow near by rescaling time with the factor Δ. This gives the following system:
where dt = Δdtn. System (6), and hence system (10), have the only equilibrium point at the origin. Therefore, the trajectories starting on leave the submanifold at some points of its boundary . Let be one of such exit points (see Figure 5B). Since the fast and slow trajectories of system (4) are glued together on , the point Osn is also the equilibrium state of the fast system (5). From this condition we find
By using Equation (11) fast system (5) can be rewritten in the form
The eigenvalues of the Jacobian matrix of (12) in the point Osn are given by
Because of Equation (13), the point Osn is a saddle-node with an unstable separatrix and a stable nodal branch. Further we show that there are such values of the parameter d that the separatrix Wu(Osn) tends to the stable node as t → +∞.
18.104.22.168. Fast epoch of motion
Consider the level curves of the function G(u1, u2) = C = const, satisfying the condition C ≤ Csn, where . Obviously, the curves of the level pass through the saddle-node Osn. In Figure 6A, the curve of this level has a maximum value, i.e., the curves of the higher levels are not indicated (white color), and the curves of the level corresponding to the values of the lower levels (C < Csn) are marked with different colors. Each color corresponds to the same value of C. Solid black lines in Figure 6 show the critical lines
At the points of these lines the following conditions are satisfied:
Consider the asymptotic behavior of the separatrix of the saddle-node Osn. Taking into account Equation (8), the location of the level curves of the function and lines (14), we establish that for the parameter value d = 0.0075 (Figure 6A) the separatrix Wu(Osn) asymptotically tends to the equilibrium state . Without changing the coordinates of the point Osn, we increase the value of the parameter d = 0.0115 (Figure 6B). For this value of the parameter d the lines merge at one point, a saddle-node bifurcation of equilibrium states occurs, and a new equilibrium state O* appear (see Figure 6B). With the further increase in the parameter, the equilibrium state O* disappears, and following the arrangement of lines (14) and the level curves of the function G, we find that in this case the separatrix Wu(Osn) tends to the node (Figure 6C). A further increase in the parameter d leads to the merging of the lines (d = 0.020909 = d*). This corresponds to the merging of the node and the saddle Osad (Figure 6C) and emerging of the saddle-node equilibrium state. For d > d*, this equilibrium state disappears, and the separatrix Wu(Osn) tends to the equilibrium state (Figure 6D). It is clear that d* depends on the coordinates of the point Osn.
Figure 6. A part of level map of the function G(u1, u2) taken at the saddle-node equilibrium state of the boundary [i.e., when and ] for (A) d = 0.0075, (); (B) d = 0.0115, (); (C) d = 0.0150, (); (D) d = 0.0225, (). Levels depicted are below the one of Osn (marked by red color). White color marks region with higher levels of G. Curves (respectively ) defined by Equation (14) are minimum and maximum of G on variable u2 (respectively variable u1). O*, Osad are additional saddle-node equilibrium states and , and are the node equilibrium states. Parameter values: a = 0.32, b = 0.79, and c = 1.166.
To describe such possible transitions, we introduced the distance R on the plane (u1, u2) of system (9) from the origin to . Note that the largest value of R corresponds to the points of the submanifolds , and .
Figure 7A depicts the results of analyzing the behavior of the separatrix Wu(Osn) for different values of belonging to the line of escape with the submanifold (see Figure 5A). For the values of from the region marked by red color, the separatrix Wu(Osn) of any of the saddle-nodes in the fast system (5) asymptotically tends to the stable node . Since in the phase space ℝ4 the equilibrium state corresponds to a stable manifold , high-amplitude oscillations are established in both elements in the system (4). Note that the trajectories of the submanifold of the slow system (10) have a similar behavior. Since the function is even, system (10) does not change when converting (u1, u2) → (−u1, −u2). Thus, transitions from to exist in ℝ4, and high-amplitude oscillations are established in system (4). The initial conditions found by us do not exhaust the entire set of initial conditions under which the oscillation amplitude changes from low to high in the second oscillator.
Figure 7. The partition of - parameter plane into areas corresponding to different transitions [along the trajectories of the fast system (7)] from the border of submanifold (A) (B) . R is the distance from the origin of slow subsystem phase plane (u1, u2) to a next stable submanifold. The red color in (A) marks the region where the next manifold is . Non-white colors in (B) mark the region where the next manifold is only one of , , , . Parameters values: a = 0.32, b = 0.79, c = 1.166, and d = 0.06.
4.1.4. Dynamics of System (4) for Initial Conditions (I.C.)2
Let us study the dynamics of system (4) for initial conditions (I.C.)2. These conditions correspond to one of the stable submanifolds of slow motions , and . Similar to the case of (I.C.)1 we have analyzed the behavior of the trajectories leaving the submanifolds. Figure 7B depicts the transitions of trajectories starting from the points at the boundary of submanifold obtained for different strength of coupling d. One can see that there are no transitions to submanifolds corresponding to high-amplitude oscillations. We have established that such behavior is also typical for other submanifolds, namely , and .
4.2. Dynamics of a Pair of Short-term Coupled Oscillators
So far, we have considered the dynamics of system (4) without any restrictions on the interaction time. However, in the initial model, the layers interact only during the time Tc. We numerically investigated the dynamics of system (3). For the initial conditions such as (I.C.)1, Figure 8A shows the behavior of 1, 000 pair of oscillators, or in other words, 1, 000 different initial conditions (I.C.)1 type, interacting during Tc = 300 and d = 0.06. For all initial conditions, after some transition process, high-amplitude oscillations are established in the interacting pairs. Figure 8B illustrates competition of oscillations in the pairs in the case of initial conditions such as (I.C.)2 type. Here, interaction of the pairs does not lead to high-amplitude oscillations, and the regime of low-amplitude oscillations persists.
Figure 8. Temporal snapshots of the pair of short-term coupled elements [system (3)] for 1,000 initial conditions such as (A) (I.C.)1; (B) (I.C.)2. Parameters values: a = 0.32, b = 0.79, c = 1.166, ε = 0.001, d = 0.06, and Tc = 300.
The occurrence of high-amplitude oscillations for the initial conditions (I.C.)1 depends on the values of the parameters Tc and d. We examine this dependence numerically and the results are presented in Figure 9. The color gradation on the plane (d, Tc) shows the dependence of the probability of establishing high-amplitude oscillations. The area highlighted in black corresponds to the establishment of high-amplitude oscillations from any initial conditions. The area highlighted in shades of gray corresponds to the establishment of high-amplitude oscillations from only some initial conditions. And finally the area marked in white corresponds to those values of the parameters for which high-amplitude oscillations are not established at all. Note that there are threshold values for both parameters d and Tc. The existence of threshold value for d has already been discussed above. The presence of threshold value for Tc is associated with the motion time T over stable submanifolds (see Figure 6C). Note that the value of the critical value Tc is determined by the dynamics of both oscillators and it is not related to the periods of high and low oscillations of the isolated oscillator.
Figure 9. Dependence of establishing probability of high-amplitude oscillations in the pair of short-term coupled elements [system (3)] for initial conditions such as (I.C.)1 on the parameters (d, Tc). Parameters values: a = 0.32, b = 0.79, c = 1.166, and ε = 0.001.
4.3. Dynamics of the Multiplex Network
Thus, it has been established that in the case of initial conditions (I.C.)1, there are values of the parameters d and Tc corresponding to the emergence of high-amplitude oscillations in the pairs of interacting oscillators belonging to the different layers. On the other hand, it has been shown that in the case of initial conditions (I.C.)2, when the oscillators belonging to the different layers do not change their initial regimes after the interaction and keep demonstrating low-amplitude oscillations.
Now let us consider the dynamics of multiplex network (1) based on the findings of the previous subsections. It can be divided into two main stages.
(a) In the time interval 0 < t < Tc, oscillators of different layers interact with each other through inter-layer couplings with strengths, dc, greatly exceeding those of the diffusive intra-layer couplings, dr. Therefore, in this stage the main contribution to the dynamics of the system due to the dynamics of the interacting pairs. We have established that as a result of this dynamics, the pairs of oscillators with high-amplitude oscillations are formed in (1) from the initial conditions (I.C.)1. On the other hand, the pairs of oscillators with low-amplitude oscillations are formed in (1) from the initial conditions (I.C.)2. This means that the average amplitude distribution in the first layer does not change, while that of the second layer becomes the same as in the first one.
(b) For t > Tc, there are no inter-layer couplings, and the oscillators interact only through diffusive intra-layer ones. Under the influence of these couplings, the neighboring oscillators with similar amplitudes become phase-locked with each other at some average frequency and form the coherent part of the chimera state. The neighboring oscillators with different amplitudes do not become phase-locked with other oscillators and form the incoherent part of the chimera state with distinguished bell-shaped distributions of average frequencies and amplitudes. Thus, the same chimera state is formed in the second layer as in the first one. Note that a finite interaction time is required to stop the competition of oscillations of pairs of oscillators. Otherwise, new complex states arise in the layers and they differ from the initial chimera.
In a large two-layer multiplex network with short-term couplings, a new phenomenon of the chimera states cloning, has been discovered and studied. Each layer of the system has a ring topology and consists of relaxation oscillators having two stable limit cycles on their phase planes. The oscillators inside the layers interact through diffusive couplings, while those of different layers interact by means of multiplex couplings. When the chimera state existing in one of the layers interacts for a while with oscillations of the other layer having a random distribution of phases, the same chimera state appears in the latter layer. Note that the time of occurrence of the chimera state in the second layer is less than the minimal partial oscillation period. We have found that the phenomenon is not related with synchronization of oscillations existing in the layers, but instead is determined by the competition of high- and low-amplitude oscillations. Using GPST, we showed that competition of oscillations in each (multiplex) pair of oscillators in the multiplex network is controlled by switching four-dimensional slow-fast dynamics. We have analytically established the initial conditions leading to the trajectories in phase space which start from stable “competitive” submanifolds of slow motions and then transit to stable “winner” submanifolds. The “competitive” submanifold corresponds to the case where oscillations in different layers have different (low and high) oscillation amplitudes. The “winner” submanifold corresponds to the case where oscillations in different layers have high amplitudes. Transitions between stable submanifolds occur along the trajectories of a two-dimensional fast subsystem. The given initial conditions belong to the basin of attraction of both the initial chimera state and the clone. We found that strength, as well as time of multiplex interaction, play a crucial role in the existence of the cloning effect of chimera states. A chimera clone is formed with 100% probability if the strength and time of multiplex interaction exceed certain threshold values. Below these threshold values a chimera clone occurs with a certain probability. Note that the effect of chimera state cloning does not depend on the choice of boundary conditions, since the dynamics of pairs of oscillators plays a crucial role in its existence. We hope also that the cloning effect is not specific to considered model and exists in other models, since the conditions necessary for it to take place are fairly general.
VN and AD conceived the original ideas. VN supervised the project. AD performed the numerical computations. VN, AD, and DS wrote the manuscript.
This work was supported by the Russian Foundation for Basic Research (Project Nos. 18-29-10040 and 18-02-00406).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
1. Kuramoto Y. Reduction methods applied to non-locally coupled oscillator systems. In: Hogan J, Champneys A, Krauskopf B, Bernardo M, Wilson E, Osinga H, et al., editors. Nonlinear Dynamics and Chaos: Where Do We Go From Here? Bristol; Philadelphia, PA: Institute of Physics Publishing (2003). p. 209–27.
2. Kuramoto Y, Battogtokh D. Coexistence of Coherence and Incoherence in Nonlocally Coupled Phase Oscillators. Nonlinear Phenomena Complex Syst. (2002) 5:380–5. Available online at: http://www.j-npcs.org/abstracts/vol2002/v5no4/v5no4p380.html
16. Semenova NI, Strelkova GI, Anishchenko VS, Zakharova A. Temporal intermittency and the lifetime of chimera states in ensembles of nonlocally coupled chaotic oscillators. Chaos (2017) 27:061102. doi: 10.1063/1.4985143
17. Shepelev IA, Vadivasova TE, Bukh AV, Strelkova GI, Anishchenko VS. New type of chimera structures in a ring of bistable FitzHugh–Nagumo oscillators with nonlocal interaction. Phys Lett A (2017) 381:1398–404. doi: 10.1016/j.physleta.2017.02.034
22. Dudkowski D, Grabski J, Wojewoda J, Perlikowski P, Maistrenko Y, Kapitaniak T. Experimental multistable states for small network of coupled pendula. Sci. Rep. (2016) 6:29833. doi: 10.1038/srep29833
26. Schmidt L, Schönleber K, Krischer K, García-Morales V. Coexistence of synchrony and incoherence in oscillatory media under nonlinear global coupling. Chaos (2014) 24:013102. doi: 10.1063/1.4858996
27. Wickramasinghe M, Kiss IZ. Spatially organized dynamical states in chemical oscillator networks: synchronization, dynamical differentiation, and chimera patterns. PLoS ONE 8:e80586. doi: 10.1371/journal.pone.0080586
28. Wickramasinghe M, Kiss IZ. Spatially organized partial synchronization through the chimera mechanism in a network of electrochemical reactions. Phys Chem Chem Phys. (2014) 16:18360–9. doi: 10.1039/C4CP02249A
31. Gambuzza LV, Buscarino A, Chessari S, Fortuna L, Meucci R, Frasca M. Experimental investigation of chimera states with quiescent and synchronous domains in coupled electronic oscillators. Phys Rev E (2014) 90:032905. doi: 10.1103/PhysRevE.90.032905
42. Maksimenko VA, Makarov VV, Bera BK, Ghosh D, Dana SK, Goremyko MV, et al. Excitation and suppression of chimera states by multiplexing. Phys Rev E (2016) 94:052205. doi: 10.1103/PhysRevE.94.052205
43. Dmitrichev AS, Shchapin DS, Nekorkin VI. Cloning of chimera states in a multiplex network of two-frequency oscillators with linear local couplings. JETP Lett. (2018) 108:543–47. doi: 10.1134/S0021364018200079
Keywords: chimera states, chimera states cloning, dynamical mechanism, bifurcations, relaxation dynamics, multiplex networks
Citation: Dmitrichev A, Shchapin D and Nekorkin V (2019) Cloning of Chimera States in a Large Short-term Coupled Multiplex Network of Relaxation Oscillators. Front. Appl. Math. Stat. 5:9. doi: 10.3389/fams.2019.00009
Received: 30 October 2018; Accepted: 25 January 2019;
Published: 18 February 2019.
Edited by:Anna Zakharova, Technische Universität Berlin, Germany
Reviewed by:Jakub Sawicki, Technische Universität Berlin, Germany
Johanne Hizanidis, University of Crete, Greece
Vadim S. Anishchenko, Saratov State University, Russia
Copyright © 2019 Dmitrichev, Shchapin and Nekorkin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Aleksei Dmitrichev, email@example.com