## Finite Kuramoto System with Shear and Symmetry

**Arindam Chakraborty**^{1}, **Anirban Ray**^{2}, **A. Roy Chowdhury**^{2,}

^{1}Department of Physics, Swami Vivekananda Institute of Science and Technology, Sonarpur, Kolkata, India

^{2}High Energy Physics Division, Department of Physics, Jadavpur University, Kolkata, India

### Abstract

In this paper we study the locally coupled finite Kuramoto Oscillator system with shear under periodic boundary condition. We also show how analytical solutions can be obtained from symmetry conditions. Existence of attractors and bifurcation patterns are revealed in an elegant way through these solutions. The synchronized regions are identified in the parameter space and critical situations are discussed. an important outcome of present analysis is the derivation of analytical form of the Poincare map.

### At a glance: Figures

**Keywords:** discreet Kuramoto oscillators, bifurcation, Poincare section

*International Journal of Physics*, 2013 1 (4),
pp 94-100.

DOI: 10.12691/ijp-1-4-3

Received July 23, 2013; Revised August 02, 2013; Accepted August 05, 2013

**Copyright**© 2013 Science and Education Publishing. All Rights Reserved.

### Cite this article:

- Chakraborty, Arindam, Anirban Ray, and A. Roy Chowdhury. "Finite Kuramoto System with Shear and Symmetry."
*International Journal of Physics*1.4 (2013): 94-100.

- Chakraborty, A. , Ray, A. , & Chowdhury, A. R. (2013). Finite Kuramoto System with Shear and Symmetry.
*International Journal of Physics*,*1*(4), 94-100.

- Chakraborty, Arindam, Anirban Ray, and A. Roy Chowdhury. "Finite Kuramoto System with Shear and Symmetry."
*International Journal of Physics*1, no. 4 (2013): 94-100.

Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |

### 1. Introduction

Diverse biological and physical phenomena from the division of a cell to the cycle of an engine piston, can all be modeled by an ensemble of oscillators due to their basic periodic properties ^{[1]}. In nature we observe many instances of systems coupled with each other in various manners. Such a collection is best represented in the elegant formulation of Kuramoto oscillators system ^{[2]}. There has been a growing interest on globally coupled limit cycle oscillators [3-8]^{[3]}. Of course majority of the analysis related to Kuramoto model deals with the situation of infinite number of such coupled oscillators ^{[9, 10]}. The natural frequencies of the elements of the population are drawn from predefined prescribed distribution. It has been shown analytically that the system with a mean-field coupling exhibits a phase transition in the thermodynamic limit, if the coupling exceeds a certain

On the other hand, as shown by Strogatz and Mirollo ^{[11]}, an infinite 1D system does not synchronize. However, in real-life application only finite systems appear. It is, therefore, interesting to see whether finite 1D Kuramoto systems may synchronize. Recently some work has been done in this direction by ^{[12, 13, 14, 15, 16]}. Another interesting configuration is Kuramoto systems with shear. This was first described by Y. Kuramoto himself ^{[17]}. Recently, Diego Paz´o et. al. studied the same system under distributed shear ^{[18, 19]}. With this point in mind we have studied a finite system of phase scillators(actually *N *= 3 and latter we have extended it to the case of *N *= 4) with the effect of shear, and where their configuration is endowed with a property of symmetry. It is quite interesting to observe in many of the situations described in the following the existence of aforesaid symmetry helps to segregate the analytical situations in easier way. Since we started with case of three oscillators, in the later stage the number was increased to *N *= 4 to have an idea about the increase of the size of the population. We then analyzed the bifurcation pattern as the parameters were varied. In this connection we may mention the work of E.C.Paulo e.t al.^{[20]} who also studied a finite system of such oscillators but without shearing force. Here we have tried to visualize this type of system which is under the influence of such forces.

**F**

**igure**

**1.**Variation of real and imaginary part of solution of Eq. (15)

### 2. Formulation

We consider a set of Kuramoto Oscillators under the influence of shearing force, which can be written as

(1) |

where is the frequency of the -th oscillator. stands for the coupling between the -th and the -th oscillator. is the corresponding number of such oscillators. is the shear of the -th oscillator which actually express the nonsynchronicity or the dependence of the frequency on the amplitude. Here, stands for

(2) |

We start by choosing that we have only three such oscillators coupled with each other and consider the fixed points obtained by the condition . The fixed point condition can be written as;

(3) |

where, stands for

(4) |

with and

whence Eq. (3) can be written as

(5) |

so from Eq. (4)

(6) |

with On the other hand Eq. (6) can also be written as

(7) |

which can be seen to lead to

(8) |

Now let us specialize to the case of and impose periodic boundary condition. For this case we have chosen , ; , . and . With this restrictions we can easily observe that Eq. (5) implies

(9) |

or

(10) |

If we set ,then we get the second condition;

(11) |

At this becomes;

(12) |

(13) |

We can analyze the function regarding their extremum. For example in case of ,if we set , then

(14) |

leads to the equation

(15) |

Solution of Eq. (14) are given by

(16) |

Then, four roots can be numerically studied as shown in Figure 1. Both the real and imaginary part of the solutions are plotted separately. It shows that for values of ,two roots become purely real. Thus, has two minimum for in the real plane. Similarly, has a single minima at These three minima give birth to a pair of solutions for . There values will help to ascertain the nature of the eigenvalues for the study of the stability. Also for the four roots (16) reduce to two. In the first case, eigenvalues are

(17) |

(18) |

where as in the second case they turns out to be

(19) |

(20) |

For this system with few oscillators it is not necessary to impose symmetry in frequency distribution to obtain the result numerically. But it is done to make the analytical solution easier. We show the above calculation with asymmetric frequency distribution in the following example. If natural frequencies are randomly chosen, it is always possible to assign positive values to and and set (for . With this configuration fixed point solutions are described by the following set of equations,

(21) |

(22) |

We find solutions of above equations in the same fashion. Solutions are given below

(23) |

(24) |

**Fig**

**ure**

**2.**(a) Plot of functions where stable solutions are shown in solid black lines and the unstable solutions in dashed red lines and (thin blue lines) for and . (b) Plot of functions of where stable solutions are shown in solid black lines and the unstable solutions in dashed red lines and (thin blue lines) for and . (c)Variation of functions on vs. plane. (d) Stable region(black) and unstable region (gray) of functions on vs. plane. (e) Plot of functions where stable solutions are shown in solid black lines and the unstable solutions in dashed red lines and (thin blue lines) for and . (f) Plot of functions of where stable solutions are shown in solid black lines and the unstable solutions in dashed red lines and (thin blue lines) for and . (g)Variation of functions on vs. plane. (h) Stable region(black) and unstable region (gray) of functions on vs. plane.

In Figure 2 explicit form of for the first set of solution. For example, in Figure we see the explicit structure of as functions of for some value . The region of stability is shown in the surface plot given in . The similar situation for is given in Figure and Figure . On the other hand the second set solution leads to the situations for and depicted respectively in Figure . To make the result more transparent we have given the two dimensional plots of plane for the and second kind of solution in Figure 3 and Figure 4. In this figure the black region gives the stable zones, whereas grey one indicate the instability and the white zones represents the negative values. One may add that the first maximum at represents a saddle(stable)node bifurcation. Atthis point two solutions are born, one stable with and one unstable solution with in the limit . The solutions generated at the minima and are unstable as the real parts of the eigenvalues are positive at this points. One should note that is a saddle node bifurcation and is a node node bifurcation. Unstable bifurcation points are indicated with the rightmost branch of . At this point a pitch-fork bifurcation takes place. To have more insight in the bifurcation pattern we consider the two equations

(25) |

(26) |

where , under the assumption and . Also we have set . To proceed further we substitute

(27) |

(28) |

and obtain the following equations

(29) |

(30) |

In equation, we have used the following nomenclature,

(31) |

The Jacobian pertaining to the equation is

(32) |

It is now convenient to define the trace and determinant of as, the help of grey shedding, while stable one by black. It is now possible to observe that solutions are starting from

(33) |

From which after some algebraic manipulation we can deduce that the following relation holds when the synchronization condition is satisfied

(34) |

which is an ellipse in the parameter plane . Since the determinant of the Jacobian is less than zero at this fixed point, the point is a saddle which changes from stable to unstable direction when the boundary of the ellipse is crossed. On the otherhand when we get

(35) |

along with

(36) |

From which one gets

(37) |

**Fig**

**ure**

**3**

**.**(a) Stable region(black) and unstable region (gray) of functions on vs. plane. (b) Stable region(black) and unstable region (gray) of functions on vs. plane

This means that for all synchronization states loose stability due to Andronov-Hopf bifurcation. On the otherhand saddle node bifurcation lines are determined implicitly by

(38) |

and

(39) |

So implies

(40) |

and

(41) |

**Fig**

**ure**

**4.**Variation of Eq. (33) in parameter plane

So may imply either or . Similar conclusion can be drawn for . So for any value of q.

### 3. Poincare Map

Further we calculate the Poincare section analytically of coupled Kuramoto equations of dimension . We transfer the system in to co-ordinate using relation .

(42) |

(43) |

Where

(44) |

(45) |

(46) |

(47) |

(48) |

where and . Hence the Poincare map becomes

(49) |

### 4. N=4 with Specular Symmetry

Based on the same assumption as the previous case, we analyzed a ring of oscillators where the configuration of frequency follow the prescribed symmetry. We assumed the situation when oscillator oscillator present specular symmetry on the natural frequency with and both and positives. Now, we have a set of equations

(50) |

(51) |

We replace , and , . Thus Eqns. (??) transform into a new set of coupled equations in and .

(52) |

(53) |

For

(54) |

(55) |

The eigenvalues of the Jacobian is given by

(56) |

For

(57) |

(58) |

The eigenvalues of the Jacobian is given by

(59) |

### 5. Two Natural Frequencies

Here

(60) |

(61) |

(62) |

(63) |

(64) |

and

(65) |

### 6. Conclusion

In the present work, we studied a exhaustive description of locally coupled Kuramoto model with shear under presence of symmetries and explored to a large extent the analyticity limits of the system. In the following we summarize the salient features. Starting with the analysis of small system, we gave a full description of synchronized region, along with determination of critical coupling for some case. We showed how fixed point may be obtained from critical synchronization coupling for some case, and we showed how the fixed points may be obtained from general configuration of natural frequency independent of system size. We also deduced the bifurcation scenario along with Poincare section analytically. But it becomes difficult to pursue the analytical result for higher dimensional system when shear is present in the system.

The author ARC is thankful to UGC(University Grant Commission) BSR(Basic Science Research) grant of Govt. of India and AR is thankful to CSIR(Council of Scientific and Industrial Research) for SRF(Senior Research Fellowship).

### References

[1] | Arthur T. Winfree. The Geometry of Biological Time, volume 8 of Biomathematics. Springer-Verlag, Berlin, Germany, 1980. | ||

In article | |||

[2] | Yoshiki Kuramoto. Chemical oscillations, waves, and turbulence. Chemistry Series. Dover Publications, 2003. originally published: Springer Berlin, New York, Heidel- berg, 1984. | ||

In article | |||

[3] | G.B. Ermentrout. Oscillator death in populations of all to all coupled nonlinear oscillators. Physica D: Nonlinear Phenomena, 41(2):219-231, 1990. | ||

In article | CrossRef | ||

[4] | Hiroaki Daido. Order function theory of macroscopic phase-locking in globally and weakly coupled limit-cycle oscillators. International Journal of Bifurcation and Chaos, 07(04):807-829, 1997. | ||

In article | CrossRef | ||

[5] | Liwei Ren and Bard Ermentrout. Phase locking in chains of multiple-coupled oscillators. Physica D: Nonlinear Phenomena, 143(14):56-73, 2000. | ||

In article | |||

[6] | Joel T. Ariaratnam and Steven H. Strogatz. Phase dia- gram for the winfree model of coupled nonlinear oscilla- tors. Phys. Rev. Lett., 86:4278-4281, May 2001. | ||

In article | CrossRef | ||

[7] | J. A. Acebr´on, A. Perales, and R. Spigler. Bifurcations and global stability of synchronized stationary states in the kuramoto model for oscillator populations. Phys. Rev. E, 64:016218, Jun 2001. | ||

In article | CrossRef | ||

[8] | Edward Ott, Paul So, Ernest Barreto, and Thomas Antonsen. The onset of synchronization in systems of glob- ally coupled chaotic and periodic oscillators. Physica D: Nonlinear Phenomena, 173(12):29-51, 2002. | ||

In article | CrossRef | ||

[9] | Juan A. Acebr´on, L. L. Bonilla, Conrad J. P´erez Vicente, F´elix Ritort, and Renato Spigler. The kuramoto model: A simple paradigm for synchronization phenomena. Rev. Mod. Phys., 77:137-185, Apr 2005. | ||

In article | CrossRef | ||

[10] | Steven H. Strogatz. From kuramoto to crawford: exploring the onset of synchronization in populations of coupled oscillators. Physica D: Nonlinear Phenomena, 143(14):1-20, 2000. | ||

In article | |||

[11] | S. H. Strogatz zand R. E. Mirollo. Collective synchronization in lattices of nonlinear oscillators with random- ness. Journal of Physics A: Mathematical and General, 21(13):L699-L705, 1988. | ||

In article | CrossRef | ||

[12] | Yu. Maistrenko, O. Popovych, O. Burylko, and P. A. Tass. Mechanism of desynchronization in the ﬁnite- dimensional kuramoto model. Phys. Rev. Lett., 93:084102, Aug 2004. | ||

In article | CrossRef PubMed | ||

[13] | V. Maistrenko, A. Vasylenko, Yu. Maistrenko, and E. Mosekilde. Phase chaos and multistability in the discrete kuramoto model. Nonlinear Oscillations,11(2):229-241, 2008. | ||

In article | CrossRef | ||

[14] | O. E. Omel’chenko, M. Wolfrum, S. Yanchuk, Y. L.Maistrenko, and O. Sudakov. Stationary patterns of coherence and incoherence in two-dimensional arrays of non-locally-coupled phase oscillators. Phys. Rev. E, 85:036210, Mar 2012. | ||

In article | CrossRef PubMed | ||

[15] | Volodymyr Maistrenko, Anna Vasylenko, Yuri Maistrenko, and Erik Mosekilde. Phase chaos in the discrete kuramoto model. International Journal of Bifurcation and Chaos, 20(06):1811-1823, 2010. | ||

In article | CrossRef | ||

[16] | BenjaminI. Triplett, DanielJ. Klein, and KristiA. Mor- gansen. Discrete time kuramoto models with delay. In PanosJ. Antsaklis and Paulo Tabuada, editors, Networked Embedded Sensing and Control, volume 331 of Lecture Notes in Control and Information Science, pages 9-23. Springer Berlin Heidelberg, 2006. | ||

In article | |||

[17] | Hidetsugu Sakaguchi and Yoshiki Kuramoto. A soluble active rotater model showing phase transitions via mutual entertainment. Progress of Theoretical Physics,76(3):576-581, 1986. | ||

In article | CrossRef | ||

[18] | Ernest Montbri´o and Diego Paz´o. Shear diversity prevents collective synchronization. Phys. Rev. Lett., 106:254101, Jun 2011. | ||

In article | CrossRef PubMed | ||

[19] | Diego Paz´o and Ernest Montbri´o. The kuramoto model with distributed shear. EPL (Europhysics Letters), 95(6):60007, 2011. | ||

In article | CrossRef | ||

[20] | Paulo F.C. Tilles, Hilda A. Cerdeira, and Fernando F. Ferreira. Local attractors, degeneracy and analyticity: Symmetry effects on the locally coupled kuramoto model. Chaos, Solitons And Fractals, 49(0):32-46, 2013. | ||

In article | |||

[21] | A. Nama junas, K. Pyragas, and A. Tamaevicius. An electronic analog of the mackey-glass system. Physics Letters A, 201(1):42-46, 1995. | ||

In article | CrossRef | ||

[22] | Satoshi Sano, Atsushi Uchida, Shigeru Yoshimori, And Ra jarshi Roy. Dual synchronization of chaos in mackey- glass electronic circuits with time-delayed feedback. Phys. Rev. E, 75:016207, Jan 2007. | ||

In article | CrossRef PubMed | ||

[23] | S. Banerjee, D. Ghosh, A. Ray, and A. Roy Chowdhury. Synchronization between two different time-delayed sys- tems and image encryption. EPL (Europhysics Letters), 81(2):20006, 2008. | ||

In article | CrossRef | ||

[24] | S. Madruga, S. Boccaletti, and M. A. Matas. Effect of a variable delay in delayed dynamical systems. Inter- national Journal of Bifurcation and Chaos, 11(11):2875-2880, 2001. | ||

In article | |||

[25] | E. M. Shahverdiev and K. A. Shore. Generalized synchronization in time-delayed systems. Phys. Rev. E, 71:016201, Jan 2005. | ||

In article | CrossRef PubMed | ||