Cluster singularities

Cluster singularities: the unfolding of clustering behavior in globally coupled oscillatory systems

Felix P. Kemeth Physik-Department, Nonequilibrium Chemical Physics, Technische Universität München, James-Franck-Str. 1, D-85748 Garching, Germany Institute for Advanced Study - Technische Universität München, Lichtenbergstr. 2a, D-85748 Garching, Germany    Sindre W. Haugland Physik-Department, Nonequilibrium Chemical Physics, Technische Universität München, James-Franck-Str. 1, D-85748 Garching, Germany Institute for Advanced Study - Technische Universität München, Lichtenbergstr. 2a, D-85748 Garching, Germany    Katharina Krischer Physik-Department, Nonequilibrium Chemical Physics, Technische Universität München, James-Franck-Str. 1, D-85748 Garching, Germany

The ubiquitous occurrence of cluster patterns in nature still lacks a comprehensive understanding. It is known that the dynamics of many such natural systems is captured by ensembles of Stuart-Landau oscillators. Here, we investigate clustering dynamics in a mean-coupled ensemble of such limit-cycle oscillators. In particular we show how clustering occurs in minimal networks, and elaborate how the observed 2-cluster states crowd when increasing the number of oscillators. Using persistence, we discuss how this crowding leads to a continuous transition from balanced cluster states to synchronized solutions via the intermediate unbalanced 2-cluster states. These cascade-like transitions emerge from what we call cluster singularities. At those points, the bifurcations of all 2-cluster states collapse and the stable balanced cluster state bifurcates into the synchronized solution supercritically. We confirm our results using numerical simulations, and discuss how our conclusions apply to spatially extended systems.

: New J. Phys.

Clustering dynamics in ensembles of coupled oscillators, that is the splitting up of the whole ensemble into two or more groups, the members of which are behaving identically, still raise many open questions, and many aspects of it still remain unsolved. Typically, clustering occurs in system with long-range interactions [1]. Oscillatory systems with global interactions play a crucial role in the understanding of phenomena observed in nature and technology, such as the visual perception in the mammalian brain, the circadian rhythm in the heart or the behavior of coupled Josephson junctions and electrochemical oscillators. Even phenomena such as synchronous chirping of crickets, the flashing of fireflies in unison and the synchronous clapping of an audience can be traced back to the action of a long range coupling between oscillating units (see Refs. [2, 3] and references therein).
If the coupling between such individual units is weak, each oscillator may be represented by its phase value only, and the resulting reduced models can be analyzed using powerful approaches as proposed by Okuda [1], Watanabe and Strogatz [4, 5] or Ott and Antonsen [6, 7]. In many physical systems, however, amplitude effects play a crucial role, and such a reduction is no longer possible. Nevertheless, it has been shown that oscillators close to the onset of oscillations can be reduced to normal forms such as the Stuart-Landau oscillators [8, 9], and thus one is able to draw general conclusion when investigating ensembles of such normal forms [10, 11].
Systems of globally coupled Stuart-Landau oscillators have been investigated intensively in recent years [12], revealing phenomena such as collective chaos [13], aging [14], chimera states [15, 16], oscillation death [17] and clustering [18, 19, 20, 21].
In this paper, we further analyze how clustering arises in systems of mean-coupled Stuart-Landau oscillators. In particular, we first investigate the stability of 2-cluster solutions in minimal networks of just two and then four oscillators. We subsequently elaborate how such cluster states bifurcate from the synchronous solution, and draw connections to previous works on clustering in globally coupled systems. We then show how new cluster states appear in phase space when increasing the number of oscillators in the network. Exploiting properties such as persistence yields new insights into how those states are arranged in phase space. Doing so, we find codimension-2 points which we dub cluster singularities. There, the balanced cluster states with two clusters of equal size supercritically bifurcate off the synchronized solution. We continue our considerations with an investigation of how ensembles of Stuart-Landau oscillators behave close to this point, and what this implies for spatially extended systems. We conclude with some open questions concerning clustering in coupled oscillators, and discuss a few promising directions to address them.

Globally coupled Stuart-Landau oscillators

It has been shown that an oscillatory system close to the onset of oscillations can be reduced to the dynamics of a so-called Stuart-Landau oscillator [8, 9]. Such analysis has also been extended to cases in which there is a field of oscillators which interact with each other. In particular, for such systems which interact linearly through the common mean of some of their variables, one obtains an ensemble of the form


with , the complex amplitude , the shear , and the coupling constant ,  [10, 11]. Hereby, the coupling is diffusive in the sense that it vanishes if  [18].
Since we take all oscillators as identical, Equation (1) is -equivariant with respect to the symmetric group , that is, with respect to permutations of the indices . In addition, the dynamics are unchanged under a rotation in the complex plane, , and under the reflection , with the bar indicating complex conjugation.
Equation (1) can be seen as a normal form for oscillatory systems with a quickly diffusing variable or a coupling through e.g. a gas phase. There is a vast number of different dynamical states that can be observed in such an ensemble, including synchronized motion and splay states [12], cluster states [1, 19, 20], chaotic dynamics [13, 22, 23] and chimera states [15].
For the synchronized solution, in which all oscillators move in phase, and the splay-state, in which all oscillators are frequency-synchronized but with a phase difference of , the respective stability boundaries can be obtained analytically, see for example Ref. [20]. For chaotic or chimera-like dynamics, however, only numerical approaches exist.
In this article, we restrict our analysis to solutions in which the ensemble splits into just two groups with each group being internally synchronized, also called a 2-cluster state [19]. Let denote the fraction of oscillators in cluster and the fraction of oscillators in cluster , then Equation (1) reduces to

with .
Following Ref. [24], the dynamics can be further reduced by exploiting the rotational invariance mentioned above and introducing polar coordinates , and the phase difference . This yields the reduced equations


describing the dynamics in a three-dimensional phase space with . Note that the roots of the reduced model, Equationss. (2) to (4), correspond to sinusoidal oscillations in the full system Eq. 1.

Balanced 2-Cluster States

First, we restrict our analysis to 2-cluster states with an equal number of oscillators in each cluster, , that is for . Then Equationss. (2) to (4) in fact describe the dynamics of just two identical oscillators. If the coupling is weak, that is if , then the amplitudes relax to , on a fast time scale, and the dynamics can be described, to a good approximation, by the phase equation only [8]. Therefore, with the adiabatic approximation the system, Equations. (2) to (4), reduces to


and thus to a sine-coupled system [1]. Notice that it has two fixed points at and , and therefore only the in-phase and anti-phase solutions exist for weak coupling. Important findings for larger ensembles of limit-cycle oscillators with weak coupling, that is for globally coupled phase oscillator systems, are summarized in [25] and [3].
The fixed points of the two-oscillator system, Equations. (2) to (4) with , can be determined as

  • , the synchronized solution in which both oscillators have amplitude equal to one and phase difference zero,

  • , the anti-phase solution corresponding to a splay state for two oscillators,

  • and two symmetry-broken solutions (see Appendix A for their derivation).

Note that the anti-phase solution exists only for parameter values , and that the locations in phase space of both the synchronized and anti-phase solutions are independent of the parameters and . Their stability boundaries can be determined by calculating the eigenvalues of the Jacobian evaluated at these solutions, and are given by the Benjamin-Feir condition [13, 12, 26],


for the synchronized solution and by


for the anti-phase solutions. It is worth mentioning that these stability boundaries are independent of the number of oscillators in the ensemble [13]. The stability diagram of these two solutions for is depicted in Fig. 1(a). There, one can observe that the synchronized solution is stable for positive values, and either loses one stable direction at the solid blue curve and subsequently a second stable direction at the dotted blue curve, or two stable directions through a Hopf bifurcation at the horizontal blue lines. Analogously, the anti-phase solution is unstable for , and gains one stable direction either at the dotted red curve and subsequently becomes stable at the solid red curve, or it becomes directly stable through a Hopf bifurcation at the horizontal red line where it gains two stable eigendirections. Note that there exist two regions in which the two solutions are bistable.
The asymmetric solutions (see Appendix A for their derivation) have the property that the amplitudes of the two oscillators differ and the oscillators have a non-zero phase difference, see also Ref. [21].

Figure 1: (a) Stability of the synchronized solution (blue) and anti-phase solution (red) for . From negative to positive , the synchronized solution , having two unstable directions, gains a stable direction at the dotted blue curve, and eventually becomes stable at the solid blue curve. At the horizontal blue line, the real parts of two complex conjugate eigenvalues of the Jacobian cross the imaginary axis, indicating a Hopf bifurcation. In contrast, the anti-phase solution loses two stable directions at the horizontal solid red line, whereas it looses one unstable direction at the solid red curve and another at the dotted red curve. This fixed point solution eventually disappears at , indicated by the dashed red line, where the amplitudes and vanish. (b) Parameter range for in which the asymmetric solutions exist. Hereby, solution exists for parameter values between the saddle-node bifurcations, indicated by the dotted green lines, and between the dashed magenta curves symbolizing pitchfork bifurcations. The boundaries of the solutions are again the dotted green lines, and the pitchfork indicated by the dashed black lines. The gray lines indicate the one-parameter continuation cuts shown in Fig. 2.

Note that, since our original system is -equivariant, every symmetry-broken state must belong to a subgroup of . In particular, for any 2-cluster states , the isotropy subgroup has the form . Furthermore, let denote the set of operations not included in the isotropy subgroup of , then there exist cluster solutions related to , which form the so-called group orbit of  [27]. For example, considering the 2-cluster solution , then is also a 2-cluster solution and , belong to the same group orbit.

Figure 2: Continuation of the fixed point solutions in the reduced model of four mean-coupled oscillators with AUTO, (a) along the vertical cut shown in Fig. 1(b) with , (b) along the horizontal cut in Fig. 1(b) with . Stable solutions are shown as solid curves, whereas unstable branches are indicated by dashed curves. Note that the branches at symbolize the same solution .

In the case of no shear, , asymmetric solutions bifurcate off the anti-phase solution for large values (see the dotted magenta and black curves for in Fig. 1(b)) through pitchfork bifurcations. For smaller values, they either bifurcate in a saddle-node with each other (see the dotted green line in Fig. 1(b)) or with the synchronized solution through pitchfork bifurcations (see the dotted magenta and black curves for in Fig. 1(b)). That means, (to be more precise, the two elements of the group orbit of ) exists between the dotted magenta curves (pitchforks) and the dotted green lines (saddle-nodes), and the solution exists in the two small regions between the dotted black curves (pitchforks) and the dotted green lines (saddle-nodes).
The continuation of the different solutions along one-parameter cuts (as indicated by the gray lines in Fig. 1(b)) is shown in Fig. 2, where each solution is represented by its variable. This means that the synchronized solution is located at and the anti-phase solution at . In Fig. 2(a), is held fixed and the synchronous solution is continued with increasing . There, one observes a subcritical pitchfork at which the solution branches of are created and gains an additional stable direction. After another pitchfork the synchronized solution becomes stable and the two branches of are born. Those two branches then reach the anti-phase solution where they get destroyed in a pitchfork bifurcation, rendering unstable. Finally, the branches of bifurcate with in another pitchfork, adding the second unstable direction.
Fixing at and starting from negative values, we observe two saddle-node bifurcations, each involving one branch of and , respectively. The branches of are annihilated through a subcritical pitchfork at the synchronized solution. Due to the symmetry , this scenario holds also for positive values. Note that there is no bifurcation at . The crossing of the two branches of is no real crossing but is a consequence of the projection of the solutions onto the variable.
The stability of the asymmetric solutions can further be investigated using the eigenvalues of the Jacobian. Then, we find that is always unstable in the regarded parameter regimes, but is stable for certain ranges of , and . In particular, the number of unstable eigendirections of in the parameter plane is shown in Fig. 3 for different values of .

Figure 3: Stability of the asymmetric solution in an ensemble of two coupled Stuart-Landau oscillators for (a) , (b) , (c) , (d) and (e) . The color encodes the number of unstable eigendirections whereas gray indicates the parameter range where the asymmetric solution does not exist. Since the system is invariant under , plots of only positive values are shown here. A magnification of the orange window is shown in Fig. 4(a).

Note that even for small regions exist in which is stable, but these regions grow for larger absolute values of , see Figs. 3(b)-(e). The stable 2-cluster states lose stability through a Hopf bifurcation, leading to the yellow patches in Fig. 3 with two unstable directions. The Hopf bifurcation is supercritical for small , creating stable periodic orbits [13]. Due to the symmetries in the equations, as discussed above, the bifurcation diagrams are symmetric under a simultaneous exchange , , leading to the left-right symmetry in Fig. 3(a).

Balanced -Cluster States

The reduced model, Equations. (2) to (4), and the eigenvalues of its Jacobian offer only limited information about the stability of 2-cluster states. This becomes obvious when considering the different ways in which 2-cluster solutions can bifurcate: either on or transverse to the cluster manifold. That is, either the two cluster clumps each remain together, but their relative locations change, or one of the two clusters splits up into an arbitrary number of subclusters. The stability properties of the former are fully captured by the reduced equations, but since the latter cannot happen in the reduced model, we can draw no conclusions about stability transverse to the cluster manifold.
In order to rephrase these arguments, we follow Ref. [20] and describe the stability of 2-cluster states by two kinds of Lyapunov exponents, the cluster integrity exponents , and the cluster system orbit stability exponents . The latter describe the stability along the 2-cluster manifold and constitute 3 real numbers. The former, the cluster integrity exponents , describe the internal stability of each cluster . Due to the symmetries of the oscillator ensemble, each is degenerate in the sense that they consist of equal real values, with being the number of oscillators in cluster .
We extend our considerations by regarding four instead of two equally sized clusters. By again introducing polar coordinates, we obtain the dynamics of the four amplitudes , , and three phase differences , . One thus has dynamics in a seven-dimensional phase space , which is equivalent to the dynamics of four coupled Stuart-Landau oscillators, each representing one cluster [16]. Note that the asymmetric solutions now correspond to cluster solutions with two oscillators in each cluster. Due to the new dimensions in phase space, however, the stability of those 2-2 cluster solutions might differ from the stability of the solutions in the two-oscillator ensemble. This can be visualized by evaluating the Jacobian of the four-oscillator system at the 2-2 cluster solutions and investigating the number of positive eigenvalues.

Figure 4: Stability of the balanced cluster solution in an ensemble of (a) two and (b) four mean-coupled Stuart-Landau oscillators for . The parameter range corresponds to the orange window indicated in Fig. 3(e). As in Fig. 3, the color encodes the number of unstable eigendirections, and gray indicates that the balanced cluster solution does not exist for this set of parameters. The bifurcation curves plotted on the right were obtained using AUTO and symbolize a supercritical Hopf bifurcation curve (blue), a pitchfork bifurcation curve (red) and a saddle-node bifurcation curve (green). The position of the cluster singularity (as discussed in the subsequent sections of this article) is indicated by the pink x.

For the parameter window shown in Fig. 4(a) (corresponding to the orange window in Fig. 3(e)), the stability of the cluster solution in the four-oscillator system for is shown in Fig. 4(b). There, one can observe that the parameter range in which the 2-2 cluster solution is stable is smaller than the parameter regime in which the solution is stable in the two-oscillator ensemble. Using the numerical continuation software AUTO, one can further observe that the 2-2 cluster solution may either become unstable through a pitchfork (red curve in Fig. 4(b)), a saddle-node (green curve in Fig. 4(b)), or through a supercritical Hopf bifurcation (blue curve in Fig. 4(b)).
Note that for the 2-2 cluster in the four-oscillator ensemble, we have cluster integrity exponents of one cluster, cluster integrity exponents of the other cluster and 3 cluster system orbit stability exponents, and thus 7 exponents in total. However, the degeneracy of the cluster integrity exponents implies that if a 2-2 cluster state is stable for certain parameters in the four oscillator ensemble, then any cluster state with an even-valued is stable for such parameters. Therefore, the parameter window shown in Fig. 4(b) indicates where balanced 2-cluster states (that is, cluster states with the same number of oscillators in each of the two clusters) for any even number of oscillators (larger than four) are stable.

Figure 5: (a) Continuation of the cluster solutions in ensemble of four coupled Stuart-Landau oscillators for and , as indicated by the dashed orange line in Fig. 4(a) but extended to more negative values, using AUTO. Dotted curves indicate repellors, whereas solid lines represent attracting solutions. (b) Stability of the balanced cluster solution in an ensemble of four mean-coupled Stuart-Landau oscillators for . The parameter range corresponds to the orange window indicated in Fig. 3(e). As in Fig. 3, the color encodes the number unstable eigendirections, and gray indicates that the balanced cluster solution does not exist for this set of parameters. The regions in which the 3-1 cluster states are stable are indicated by the white-shaded patches.

The question remains how this 2-2 cluster state bifurcates from or to the synchronized solution. In order to investigate this, we look at a cut in parameter space (indicated by the dashed orange line in Fig. 4(b)) and continue the 2-2 cluster solution by varying for fixed and . The amplitude of one cluster, , from an exemplary continuation with AUTO is depicted in Fig. 5(a). There, the 2-2 cluster state is indicated by blue curves, which are dotted when this cluster solution is unstable and solid when the 2-2 cluster is stable. Note that there exist two varieties of the 2-2 cluster state, as indicated by the upper and lower blue line in Fig. 5(a). These two solutions differ in their respective value of the amplitude , but are both stable for the same values. Moreover, they belong to the same group orbit (as they can be transformed into one another by interchanging the oscillators). Furthermore one can observe that the 2-2 cluster states bifurcate into 2-1-1 cluster states (green and orange curves in Fig. 5(a)), which are unstable for the parameter values regarded here. These 2-1-1 branches, however, subsequently bifurcate into a 3-1 cluster branch (red curve in Fig. 5(a)), which shows two stable regions. In contrast to the stable 2-2 cluster solutions, the stable cluster solutions do not belong to the same group orbit, but correspond to a state in which the cluster with 3 oscillators has an amplitude (for ) and to a state in which the cluster with three oscillators has an amplitude (for ). Note that at high respective low values the stable 3-1 cluster states become unstable through saddle-node bifurcations. In addition, the cluster states bifurcate with the synchronized solution (black line in Fig. 5(a)) through a transcritical bifurcation, a scenario which has already been described in Refs. [28, 19]. This is in contrast to the ensemble of symmetrically related 2-2 cluster states, which bifurcate off the synchronized solution through equivariant pitchfork bifurcations at the same points. The overlapping stable 2-2 and 3-1 cluster solutions shown in Fig. 5(a) also explain the hysteretic transitions from balanced cluster states to homogeneous oscillations and vice versa. The regime in which the cluster solution is stable is shown as a shaded region in Fig. 5(b).

Figure 6: Cluster solutions in an ensemble of 16 coupled Stuart-Landau oscillators obtained using AUTO. Shown is (a) the amplitude of the larger cluster, , as a function of the bifurcation parameter , as indicated by the dashed orange line in Fig. 5(b), and (b) the amplitude of the smaller cluster, . Other parameter values are and .
Figure 7: The amplitude of one of 16 oscillators as a function of the parameter starting from the stable synchronized solution and increasing (blue curve) and when subsequently reducing (orange curve). Due to the addition of finite noise in the numerical simulations when increasing , not all of the densely located cluster states close to are resolved.

The analysis from above can easily be extended to larger ensembles of oscillators. In particular, we now consider oscillators and investigate the clustering behavior along the parameter cut indicated by the dashed orange line in Fig. 4(b). That is, we fix and and vary . We do this by simulating the full model starting from the balanced cluster solution, and increasing stepwise. Over the course of this increase, the system goes from a stable solution to a stable solution, solution and so forth until it settles on the synchronized solution. In this way we obtain solutions for every cluster distribution, which we then use for continuation with AUTO. The continuation curves we obtain this way are depicted in Fig. 6(a), where again the amplitude of the cluster with the larger number of oscillators, , is shown as a function of the continuation parameter . The amplitudes of the smaller cluster for the same parameter window are shown in Fig. 6(b). Note the correspondence of the cluster solution, shown in blue, with the cluster solution in Fig. 5(a).
There are two values, and , between which the synchronized solution is unstable, see Fig. 6. At the bifurcation points, there are also transcritical bifurcations of all unbalanced cluster states and pitchfork bifurcations of the balanced state. The two solution branches of each transcritical and of the pitchfork bifurcations, respectively, connect the two bifurcation points at and . Each of the unbalanced clusters is born respectively destroyed in a saddle-node bifurcation, the most outer one corresponding to the -1 cluster state being the only one that posses a stable branch. The other cluster states, as well as the balanced one, are stabilized through further (equivariant) pitchfork bifurcations (not shown), compare Fig. 5(a). In this way, two staircases of overlapping stable cluster states are generated, whereby the cluster distribution of the cluster states in subsequent steps differ by just one oscillator. This leads to two cascades of transitions between the two synchronized regions.
If we start from the stable synchronized solution for and slowly increase , the system goes from the synchronized solution to a cluster state, and then two a state and so forth, traversing a cascade up to the balanced 8-8 cluster and back, until it settles again on the synchronized solution, see the blue curve in Fig. 7. Thereby, the originally larger cluster with amplitude (cf. Fig. 6(a)) becomes the smaller one with amplitude (cf. Fig. 6(b)) beyond the balanced cluster state. The second staircase, and with it the hysteretic behavior, can be seen when we subsequently reduce again, cf. the orange curve in Fig. 7.
Furthermore, from Fig. 6(a) it becomes obvious that all unbalanced cluster solutions, that is all cluster solutions except the cluster, bifurcate with the synchronized solution in a transcritical bifurcation. In particular, the unbalanced cluster states exist on both sides of both bifurcation points where the synchronized solution changes stability. However, all cluster solutions lose stability through equivariant pitchfork bifurcations (the resulting branches are not shown in Fig. 6(a)), in which either the cluster with the larger or smaller amplitude breaks up. This is true except for the most unbalanced, the , cluster, for which the smallest cluster cannot break up. There we find that this stable state is destroyed in a saddle-node bifurcation instead. In addition, each kind of unbalanced cluster solution is stable in two different regions in parameter space. Each of these stable regions lies close but slightly shifted to the stable regions of neighboring cluster states.

Persistence and Cluster Singularities

The observed phenomena of slightly shifted bifurcations can be explained with the concept of persistence. Loosely speaking, if two attractors are hyperbolic and close in phase space, then bifurcations of those attractors are also close in parameter space [19]. In addition, one can infer that between the cluster states shown in Fig. 6, which must also exist for larger ensembles, there lie many more cluster states in larger networks, and using persistence, their stability must be similar to that of the ones seen in Fig. 6.
This explains the cascade-like transition from balanced cluster states to the homogeneous solution in large ensembles, where, when changing a parameter, one oscillator after another joins the other cluster until the synchronized solution is reached. We conjecture that for infinitely large ensembles, the cluster attractors are infinitesimal close, and thus this process becomes continuous.
Turning back to Fig. 4(b), we observe that there is a codimension-2 point (pink point in Fig. 4(b)) where the stable 2-2 cluster bifurcates into the synchronized solution in a pitchfork bifurcation. This is in contrast to the phenomena observed in the literature, where the transition to the synchronized solution occurs via the unbalanced cluster solutions [20]. For the Stuart-Landau ensemble, Equation (1), this point can be found analytically as


with the derivation shown in Appendix B.
The characteristics of such a point is that, when starting from a balanced 2-cluster solution and changing the parameters over this bifurcation point, the two clusters approach each other and finally merge and form the synchronized solution. This is what we call a cluster singularity. However, when varying the parameters such that one turns around the codimension-2 point either clock- or anticlockwise (that is, changing and along a path which circumvents the singularity on the left or on the right, cf. pink arrows in Fig. 4(b)), then either of the clusters shrinks and single oscillators join the other cluster until all oscillators finally form the synchronized solution. This scenario can be verified using numerical simulations and is shown in Fig. 8. There, simulations of Stuart-Landau oscillators are shown when avoiding the cluster singularity clockwise (Fig. 8(a)), when directly crossing over the cluster singularity (Fig. 8(b)) and when avoiding the cluster singularity in an anticlockwise manner (Fig. 8(c)).
The cluster singularity serves as an organizing center for nearby unbalanced cluster solutions. Recall that all unbalanced 2-cluster solutions get destroyed in saddle-node bifurcations, cf. Fig. 6. In the cluster singularity, all these saddle-node bifurcations as well as the pitchfork bifurcations that alter the stability properties of the cluster states collapse to a single point in phase space, suggesting the name cluster singularity. Note that when crossing the singularity, the stable balanced cluster solution directly bifurcates into the synchronized solution.

Figure 8: Simulations of the globally coupled Stuart-Landau ensemble close to the cluster singularity for oscillators and , (a) , (b) and (c) . The direction in which is changed is from small to large values.

Clustering in Spatially Extended Systems

Adding a diffusive coupling in addition to the global coupling, one obtains a globally coupled version of the complex Ginzburg-Landau equation [29],


with and . In a sense, such a system can be viewed as an ensemble of infinitely many oscillators, coupled locally and globally. If the local coupling is weak (), then we expect the solutions of the Stuart-Landau ensemble to exist also in the spatially extended system. For infinitely many oscillators, however, the 2-cluster solutions become infinitesimally close in phase space (cf. Fig. 6 for 16 oscillators) and thus infinitesimally small perturbations are sufficient to drive the solution from one cluster state to another. This is also what we observe in numerical simulations: the diffusive coupling leads to the selection of a particular cluster distribution, and the multistability of different cluster solutions, as apparent in Fig. 6 for Stuart-Landau oscillators, seems no longer to exist. What is special about the then globally stable 2-cluster solution, however, still remains unknown.


To summarize our results, we have shown ways how clustering can occur in globally coupled ensembles of Stuart-Landau oscillators. In particular, starting from small ensembles, we described how 2-cluster branches bifurcate, and extended this analysis to larger ensembles of oscillators. Doing so, we found a codimension-2 point which we dubbed a cluster singularity: at this point, the stable balanced cluster solution bifurcates directly into the synchronized solution. In addition, all saddle-node bifurcations generating unbalanced cluster solutions collapse in this point. Using numerical simulations, we showed how ensembles of Stuart-Landau oscillators behave close to this cluster singularity. Since any oscillatory system close to the onset of oscillations can be mapped onto the dynamics of the Stuart-Landau oscillator, we believe that cluster singularities are common in oscillatory systems with global coupling, and that an experimental observation of these should be possible. Concludingly, we discussed how our results extend to spatially extended systems, where the diffusive coupling seems to destroy the multistability.
We believe that our considerations may serve as a further step towards a better understanding of clustering behavior in coupled oscillators. In addition, 2-cluster solutions in the regarded parameter windows may also become unstable through supercritical Hopf bifurcations for smaller values, and even bifurcate into chimera states [13, 16]. How this transition occurs for different cluster distributions is still an open question.

The authors thank Matthias Wolfrum, Oliver Junge, and Munir Salman for fruitful discussions. Financial support from the Deutsche Forschungsgemeinschaft (Grant no. KR 1189/18-1), the Institute of Advanced Study - Technische Universität München, funded by the German Excellence Initiative, and the Studienstiftung des deutschen Volkes is gratefully acknowledged.

Appendix A: Balanced cluster solution

In order to find the solutions of balanced 2-cluster states, it is sufficient to find the fixed points of the reduced two-oscillator system


In addition, these equations can be simplified by introducing the sum and the difference of the squared amplitudes, and , with

This transforms Equations. (11) to (12) into

and using ,

At a fixed point solution, this system of equations must satisfy

Solving the first two equations for and yields


and inserted into the last equation,

Solving for


Using the identity , we can write Equations. (14) and (15), yielding


By inserting Equation (16) into Equation (17), we can solve it for and obtain


Together with Equation (16), this can be used to calculate , and .

Appendix B: Cluster Singularities

The idea is that at the cluster singularity, the saddle-node bifurcations of all unbalanced 2-cluster solutions hit the pitchfork at which the synchronous solution becomes unstable. This must be true for any . Therefore, for simplicity, we take the limit in system Equationss. (2) to (4), yielding

This means and thus leaves

Setting means we are at the point where the cluster solution meets the synchronous solution ( follows from ), and from this the previous expression turns into


which coincides with the curve at which the homogeneous solution becomes unstable. For the saddle-node curve of the cluster solutions, the two solutions of must equal, and thus the discriminant must equal zero

Now use that from Equation (19) above,


Equation (20) gives the saddle-node curve of the cluster with . So for the cluster singularity, this saddle-node bifurcation coincides with the point at which the homogeneous solution becomes unstable, as given by Equation (19), which yields

Subtraction of these two equations yields


and thus

This solution plugged into Equation (21) yields

So, in total, we have at the cluster singularity


This indicates two possible solutions for the cluster singularity. Furthermore, it is worth mentioning that it seems to exist for all values.

Appendix C: Numerical Methods

For the integration of the Stuart-Landau ensemble, an implicit Adams method with a fixed time step of is used. All figures are generated using matplotlib [30].



  • [1] Koji Okuda. Variety and generality of clustering in globally coupled oscillators. Physica D: Nonlinear Phenomena, 63(3-4):424–436, 1993.
  • [2] Steven H. Strogatz. From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators. Physica D: Nonlinear Phenomena, 143(1):1–20, 2000.
  • [3] Arkady Pikovsky and Michael Rosenblum. Dynamics of globally coupled oscillators: Progress and perspectives. Chaos: An Interdisciplinary Journal of Nonlinear Science, 25(9):097616, 2015.
  • [4] Shinya Watanabe and Steven H. Strogatz. Integrability of a globally coupled oscillator array. Phys. Rev. Lett., 70:2391–2394, Apr 1993.
  • [5] Shinya Watanabe and Steven H. Strogatz. Constants of motion for superconducting Josephson arrays. Physica D: Nonlinear Phenomena, 74(3):197–253, 1994.
  • [6] Edward Ott and Thomas M. Antonsen. Low dimensional behavior of large systems of globally coupled oscillators. Chaos: An Interdisciplinary Journal of Nonlinear Science, 18(3):037113, 2008.
  • [7] Edward Ott and Thomas M. Antonsen. Long time evolution of phase oscillator systems. Chaos: An Interdisciplinary Journal of Nonlinear Science, 19(2):023117, 2009.
  • [8] Y. Kuramoto. Chemical Oscillations, Waves and Turbulence, volume 19. Springer-Verlag Berlin Heidelberg, 1984.
  • [9] Vladimir García-Morales and Katharina Krischer. The complex Ginzburg-Landau equation: An introduction. Contemporary Physics, 53(2):79–95, 2012.
  • [10] Vladimir García-Morales and Katharina Krischer. Normal-form approach to spatiotemporal pattern formation in globally coupled electrochemical systems. Phys. Rev. E, 78:057201, Nov 2008.
  • [11] Vladimir García-Morales and Katharina Krischer. Nonlocal complex Ginzburg-Landau equation for electrochemical systems. Phys. Rev. Lett., 100:054101, Feb 2008.
  • [12] Vincent Hakim and Wouter-Jan Rappel. Dynamics of the globally coupled complex Ginzburg-Landau equation. Physical Review A, 46(12):R7347–R7350, 1992.
  • [13] N. Nakagawa and Y. Kuramoto. Collective chaos in a population of globally coupled oscillators. Progress of Theoretical Physics, 89(2):313–323, 1993.
  • [14] Hiroaki Daido and Kenji Nakanishi. Aging and clustering in globally coupled oscillators. Phys. Rev. E, 75:056206, May 2007.
  • [15] Gautam C. Sethia and Abhijit Sen. Chimera states: The existence criteria revisited. Phys. Rev. Lett., 112:144101, Apr 2014.
  • [16] Felix P. Kemeth, Sindre W. Haugland, and Katharina Krischer. Symmetries of chimera states. Phys. Rev. Lett., 120:214101, May 2018.
  • [17] Mauparna Nandan, C. R. Hens, Pinaki Pal, and Syamal K. Dana. Transition from amplitude to oscillation death in a network of oscillators. Chaos: An Interdisciplinary Journal of Nonlinear Science, 24(4):043103, 2014.
  • [18] D.G. Aronson, G.B. Ermentrout, and N. Kopell. Amplitude response of coupled oscillators. Physica D: Nonlinear Phenomena, 41(3):403–449, 1990.
  • [19] Murad Banaji. Clustering in globally coupled oscillators. Dynamical Systems: An International Journal, 17(3):263–285, 2002.
  • [20] Wai Lim Ku, Michelle Girvan, and Edward Ott. Dynamical transitions in large systems of mean field-coupled landau-stuart oscillators: Extensive chaos and cluster states. Chaos: An Interdisciplinary Journal of Nonlinear Science, 25(12):123122, 2015.
  • [21] André Röhm, Kathy Lüdge, and Isabelle Schneider. Bistability in two simple symmetrically coupled oscillators with symmetry-broken amplitude- and phase-locking. Chaos: An Interdisciplinary Journal of Nonlinear Science, 28(6):063114, 2018.
  • [22] Naoko Nakagawa and Yoshiki Kuramoto. From collective oscillations to collective chaos in a globally coupled oscillator system. Physica D: Nonlinear Phenomena, 75(1):74 – 80, 1994.
  • [23] Naoko Nakagawa and Yoshiki Kuramoto. Anomalous Lyapunov spectrum in globally coupled oscillators. Physica D: Nonlinear Phenomena, 80(3):307 – 316, 1995.
  • [24] Arkady Pikovsky, Michael Rosenblum, and Jürgen Kurths. Mutual synchronization of two interacting periodic oscillators, pages 222–235. Synchronization. Cambridge University Press (CUP), 2001.
  • [25] Yoshiki Kuramoto and Ikuko Nishikawa. Statistical macrodynamics of large dynamical systems. case of a phase transition in oscillator communities. Journal of Statistical Physics, 49(3-4):569–605, 1987.
  • [26] T. Brooke Benjamin and J. E. Feir. The disintegration of wave trains on deep water part 1. theory. Journal of Fluid Mechanics, 27(03):417, 1967.
  • [27] Martin Golubitsky and Ian Stewart. The Symmetry Perspective: From Equilibrium to Chaos in Phase Space and Physical Space. Birkäuser Verlag, Basel, Boston, Berlin, 2003.
  • [28] P. Ashwin and J. W. Swift. The dynamics of n weakly coupled identical oscillators. Journal of Nonlinear Science, 2(1):69–108, Mar 1992.
  • [29] D. Battogtokh and A. Mikhailov. Controlling turbulence in the complex Ginzburg-Landau equation. Physica D: Nonlinear Phenomena, 90(1):84–95, 1996.
  • [30] J. D. Hunter. Matplotlib: A 2d graphics environment. Computing In Science & Engineering, 9(3):90–95, 2007.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description