Nonlinear resonance and synchronization in the ring of unidirectionally coupled Toda oscillators

Abstract
In the ring of unidirectionally coupled Toda oscillators the nonlinear resonance and the synchronization are investigated. It is shown how the nonlinear resonance affects the structure of the main synchronization region. As a result of nonlinear resonance we observe the coexistence of two stable limit cycles near the resonant frequency, which leads to coexistence of periodic and quasi-periodic regimes within the synchronization region.



Introduction
Oscillator chains can be divided into two groups: self-excited, which demonstrate stable limit cycle without external driving [1][2][3][4][5] and strictly dissipative chains, which have stable steady state without forcing [6][7][8][9]. For the ring of unidirectionally coupled Duffing oscillators [2] and for the ring of unidirectionally coupled Toda oscillators with nonlinear coupling function [3], it has been shown that with the increase of coupling coefficient the stable fixed point undergoes Andronov-Hopf bifurcation where the stable limit cycle may appear [3,4]. Hence, the system can reach a stable fixed point (when the coupling coefficient of the system is small) or stable limit cycle (when the coupling coefficient is above the threshold of Andronov-Hopf bifurcation). The systems of coupled Toda and Duffing oscillators demonstrate the effect of nonlinear resonance [6][7][8][9][10][11]. Simulationousely as for all self-oscillatory coupled systems we expect the appearence of synchronization [12][13][14][15]. It has been shown, that in the case of strong nonlinearity the main synchronization region of self-oscillatory systems can have some peculiarities, such as asymmetry and stretched form of synchronization tongue [16,17]. Neimark-Sacker curves for such systems do not meet saddle-node curves of limit cycles at its cusp points. It is interesting to see, how the nonlinear resonance affects the synchronization in the system and the peculiarities of the structure of the main synchronization region.
In this paper we consider unidirectionally coupled Toda oscillators with ring topology of coupling. It has been shown, that increasing the coupling strength leads to selfoscillation birth [2,3]. Hence, the system which is investigated in our paper, represents a special type of generator under external force. The most interesting effect which takes place in the forced generators is synchronization. We show that the main synchronization region has classical structure only for small amplitude of driving force. For larger values of amplitude, we observe ranges where different solutions coexist, hence it is non-classical behavior. To explain this non-typical structure, we study the resonance in the system under the threshold of self-oscillation birth.

The model
We consider the ring of three unidirectionally coupled Toda oscillators with the harmonic force acting on one of them: where x 1,2,3 , y 1,2,3 , are the dynamical variables of the system, α is the damping coefficient, γ is the coupling coefficient (γ > 0) and H(x) = exp(x) − 1 is the coupling function.
The investigations of the system have been conducted using the XppAut package [18]. Limit cycle branches and bifurcations of limit cycles have been obtained by Auto07p [19] software, time series and projections of phase portraits have been obtained by numerical integration of differential equations using Runge-Kutta method.
The system without external force (F = 0) with α = 0.1 has the stable steady state for the coupling parameter γ < 0.1189, Andronov-Hopf bifurcation occurs at γ = 0.1189 and the stable limit cycle is born. Bifurcation diagram is shown in Fig. 1(a). Further increase of coupling coefficient leads to Neimark-Sacker bifurcation of the limit cycle and the appearance of complex dynamics. In all figures we mark the equilibrium as EQ and limit cycles with letter C using numerical subscripts, if we observe more than one solution.
The periodic regime is characterized by phase shifts between partial oscillators equal to T /3, where T is period of oscillations (see Fig. 1(b)).

Nonlinear resonance in the system
To observe resonance in the system (1), we choose two values of the coupling coefficient γ, which correspond to stable steady state in the system without external forcing. The first value fixed as γ = 0 demonstrates the resonance properties of partial oscillator. Bifurcation structure of the driven Toda oscillator on the parameter space of amplitude and frequency of driving force is studied [20]. Nonlinear resonance leads to the appearance of «resonance horns» on the parameter space -areas, which are formed by lines of saddle-node bifurcations. Inside these areas two resonance limit cycles coexist. Figure 2(a) shows the resonance curve of the system with γ = 0, where the resonance 1:1 is clearly seen whereas Fig. 2(c) shows the bottom part of resonance horn 1:1. Figure 2(b) shows the resonance curve when γ = 0.11. Here, we can see two resonances, which are signed by 1 and 2. The corresponding two-parameter diagram is shown in Fig. 2(b). Here we can see two resonance horns for these two resonances labelled by 1 and 2.
Time series of x 1 (t), x 2 (t) and x 3 (t) for coexisting resonance limit cycles for the first and the second resonances are shown in Figs. 3(a, b) and 3(c, d) respectively. Let us call first resonance as partial oscillator resonance and the second resonance as ring resonance. We can see, that phase shifts between oscillations for the first resonance are near zero, and for the second resonance oscillations are similar to the self-oscillations in the autonomous system ( Fig. 1): the amplitudes of partial oscillators are different, but phase shifts are T /3, where T is the period of the driving force.
We can see, that the coupling in the system leads to the appearance of the new resonance, determined by the ring structure of the coupled oscillators.

The structure of the main synchronization region
In order to observe synchronization in the system we fix the internal parameter values as α = 0.1, γ = 0.12 (which corresponds to the stable limit cycle in an autonomous system) and investigate the main synchronization region on the plane (ω, A) of the external force parameters.
The results of the detailed investigations of dynamical regimes and transitions between them are shown in Fig. 4. In region A, there are three limit cycles: stable C 1 , saddle C 2 and C 3 (Fig. 5(a)). On the lines l 1 SN and l 1 SN there are saddle-node bifurcations of cycles C 1 and C 2 on torus ( Fig. 6(a) shows the bifurcation). On the line l 3 SN there is saddle-saddle bifurcation of cycles C 2 and C 3 (Fig. 6(b)), so in the region B only one cycle C 1 is left (Fig. 5(b)). In region C between the lines l 1 SN and l T D there are hysteretic transitions between the limit cycle and the torus, in this region the stable torus and the limit cycle coexist (Fig. 5(c)). On line l T D the coexisting torus disappears. On line l 3 NS the torus converges to cycle C 3 , which becomes stable and in region D two stable cycles C 1 and C 3 coexist (Fig. 5(d)). On line l 3 SN there is saddlenode bifurcation of cycles C 2 and C 3 , so in region B still only one limit cycle C 1 is left. Lines l 3 SN , l 3 SN , l 3 NS and l T D are connected in the point p 3 . The bifurcation diagram   for the case when we are under this point is shown in Fig. 6(c), when we are above this point -in Fig. 6(d). On line l 1 NS there is Neimark-Sacker bifurcation of cycle C 1 . Summarizing, in regions A and B we see only one stable limit cycle in the phase space of the system. In region C we see the coexisting limit cycle and the torus. In region D we observe two coexisting limit cycles. The regions where only the torus is stable are marked by T.
Such a structure of the synchronization area is closely connected with the appearance of the new resonance. Figure 7 shows the bifurcation lines of the system under the threshold of self-oscillation birth (a) and over it (b). Under the threshold of self-oscillation birth we can see saddle-node bifurcation lines of resonance limit cycles, which form the area of coexisting regimes. Over the threshold of self-oscillation birth in the bottom part of the parameter plane (by small values of external force amplitude) we see classical Arnold tongue formed by saddle-node bifurcation lines on a torus. In the upper part of the parameter plane we see the coexisting resonance cycles, as it was under the threshold of self-oscillation birth but now each of them undergoes Neimark-Sacker bifurcation.
Hence, in the bottom part of the parameter plane we observe the classical case of synchronization. In the upper part we see the effects of the nonlinear resonance in the ring.
In order to see the processes in the system which take place by the transition through the threshold of self-oscillation birth under external force, we fix the driving force amplitude as F = 0.015 and increase the coupling coefficient. We observe the dynamical regimes of the system on the parameter plane (ω, γ) (Fig. 8). For the weak coupling only one limit cycle is observed ( Fig. 9(a), Fig. 8 region B). With the increase of the coupling coefficient the softening characteristic of stiffness is visible, hence we  observe coexistence of the stable limit cycles ( Fig. 9(b)). In Fig. 8 between the lines l 1 SN and l 3 SN (region D) two stable limit cycles C 1 and C 3 and saddle cycle C 2 are observed. Point p 4 here is the cusp point of saddle-node bifurcation lines l 1 SN and l 3 SN . Over the threshold of self-oscillation birth (which occurs in the system without driving at γ = 0.1189) each of the stable limit cycles undergo an Neimark-Sacker bifurcation ( Fig. 9(c)) on lines l 1 NS and l 3 NS (for cycles C 1 and C 3 respectively). Over line l 3 NS in region C, there are coexisting limit cycle C 1 and torus (it disappears on line l T D ), so over this line in region A only one stable limit cycle C 1 is left. Point C 2 and C 3 . On the upper side of the parameter plane we see another saddle-node bifurcation line l 1 SN for cycles C 1 and C 2 . Bifurcation diagram for this case can be seen in the Fig. 9(d). Here, we can observe saddle-node bifurcations of cycles C 1 and C 2 on lines l 1 SN while cycle C 3 is unstable. Point p 2 is the point where line l 1 NS meets line l 1 SN . On line l 3 SN there is saddle-saddle bifurcation of limit cycles C 1 and C 2 . Hence, in region B we can only see only one stable limit cycle in the phase space of the system, which corresponds to the forced oscillations of the system. In region D we see two coexisting limit cycles, while in region C we observe the coexisting limit cycle and the torus. In region A there is one stable limit cycle. The regions where only the torus is stable are marked by T.
We can see that under the threshold of self-oscillation birth two coexisting stable limit cycles can be observed in the phase space of the system. Over the threshold of self-oscillation birth, the coexisting limit cycles undergo Neimark-Sacker bifurcations, which lead to the coexistence of periodic and quasi-periodic solutions.

Conclusion
The main conclusion from our paper is that the peculiarities in the structure of synchronization region such as coexisting periodic and quasi-periodic solutions are caused by nonlinear resonance. When the coupling parameter of the system is close to zero, we see only one resonance near the frequency of partial oscillators. Under the threshold of self-oscillation birth, the increase of the coupling strength leads to the appearance of the new resonance, which we call a ring resonance, which is caused by the ring structure of the system (three coupled systems). Nonlinear resonances lead to the coexistence of two limit cycles. Over the threshold of self-oscillation birth we see the interaction between external oscillations and the inner oscillatory mode of the system, which leads to the appearance of the synchronization region in the parameter plane. The bottom part of this synchronization region has the classic form of Arnold tongue, while the upper part of it is the result of nonlinear resonance (coexistence of different attractors).