Abstract

We show that dynamical clustering, where a system segregates into distinguishable subsets of synchronized elements, and chimera states, where differentiated subsets of synchronized and desynchronized elements coexist, can emerge in networks of globally coupled robust-chaos oscillators. We describe the collective behavior of a model of globally coupled robust-chaos maps in terms of statistical quantities and characterize clusters, chimera states, synchronization, and incoherence on the space of parameters of the system. We employ the analogy between the local dynamics of a system of globally coupled maps with the response dynamics of a single driven map. We interpret the occurrence of clusters and chimeras in a globally coupled system of robust-chaos maps in terms of windows of periodicity and multistability induced by a drive on the local robust-chaos map. Our results show that robust-chaos dynamics does not limit the formation of cluster and chimera states in networks of coupled systems, as it had been previously conjectured.

1. Introduction

Many smooth nonlinear dynamical systems possess chaotic attractors embedded with a dense set of periodic orbits for any range of parameter values. Therefore, in practical systems operating in chaotic mode, a slight perturbation of a parameter may drive the system out of chaos. Alternatively, there exist dynamical systems that exhibit the property of robust chaos [16]. A chaotic attractor is said to be robust if, for its parameter values, there exists a neighborhood in the parameter space where windows of periodic orbits are absent and the chaotic attractor is unique [1].

Robust chaos constitutes an advantageous feature in applications that require reliable functioning in a chaotic regime, in the sense that the chaotic behavior cannot be removed by arbitrarily small fluctuations of the system parameters. For instance, networks of coupled maps with robust chaos have been efficiently used in communication and encryption algorithms [7] and they have been investigated for information transfer across scales in complex systems [8]. In addition, the existence of robust chaos allows for heterogeneity in the local parameters of a system of coupled oscillators, while guaranteeing the performance of all the oscillators in a chaotic mode.

On the other hand, systems possessing robust chaos may present limitations in the types of collective behaviors that they can achieve, in comparison with systems displaying periodic windows. For example, it has been conjectured that the phenomenon of dynamical clustering in globally coupled networks (where the system segregates into distinguishable subsets of synchronized elements) is only found when stable periodic windows are present in the local elements [911]. Recently, it has also been argued that chimera states (i.e., coexistence of subsets of oscillators with synchronous and asynchronous dynamics) cannot emerge in networks of coupled oscillators having robust chaotic attractors [12, 13].

The phenomenon of dynamical clustering is relevant as it can provide a simple mechanism for the emergence of differentiation, segregation, and ordering of elements in many physical and biological systems [14, 15]. Clustering has been found in systems of globally coupled Rössler oscillators [16], neural networks [17], and biochemical reactions [18] and has been observed experimentally in arrays of globally coupled electrochemical oscillators [19] and globally coupled salt-water oscillators [20]. In addition, the study of chimera states currently attracts much interest (for reviews see [21, 22]). Chimera states have been found in networks of nonlocally coupled phase oscillators [23, 24], in systems with local [2528] and global [2934] interactions, and in networks of time-discrete maps [3538]. These states have been investigated in a diversity of contexts [3947]. Chimera states have been observed in experimental settings, such as populations of chemical oscillators [48], coupled lasers [49], optical light modulators [50], and electronic [51] and mechanical [52, 53] oscillator systems. It has been shown that clustering is closely related to the formation of chimera states in systems of globally coupled periodic oscillators [31].

In this paper, we investigate the occurrence of dynamical clustering and chimera states in systems of coupled robust-chaos oscillators. In Section 2, we describe the characterization of synchronization and cluster and chimera states in globally coupled systems. In Section 3, we consider a network of globally coupled robust-chaos maps and show that cluster and chimera states can actually emerge in this system for several values of parameters. In Section 4, we employ the analogy between the local dynamics of the globally coupled system with the response dynamics of a single driven map. We interpret the occurrence of clusters and chimeras in the globally coupled system in terms of windows of periodicity induced by the drive on the local robust-chaos map. Conclusions are presented in Section 4.

2. Methods

A global interaction in a system can be described as a field or influence acting on all the elements in the system. As a simple model of an autonomous dynamical system subject to a global interaction, we consider a system of maps coupled in the formwhere describes the state variable of the th map in the system at discrete time , the function expresses the local dynamics of the maps, the function represents a global field that depends on the states of the elements in a given subset of the system, at time , and the parameter measures the strength of the coupling of the maps to the field. The form of the coupling in equation (1) is assumed in the commonly used diffusive form. The function may not depend on all the elements, but it must be shared by all the elements of the system to be a global interaction.

A collective state of synchronization or coherence takes place in system equation (1) when for asymptotic times. A desynchronized or incoherent state corresponds to for all times. Dynamical clustering occurs when the system segregates into a number of distinct clusters or subsets of elements such that elements in given subset are synchronized among themselves. In other words, in the th cluster, where denotes the value of in that cluster, with . If is the number of elements belonging to the th cluster, then its relative size is . In general, the number of clusters, their size, and their dynamical evolution (periodic, quasiperiodic, or chaotic) depend on the initial conditions and parameters of the system. A chimera state consists of the coexistence of one or more clusters and a subset of desynchronized elements. If there are clusters, the fraction of elements in the system belonging to clusters is while is the fraction of elements in the desynchronized subset.

In practical applications, we consider that two elements and belong to a cluster at time if the distance between their state variables, defined asis less than a threshold value , i.e., if . The choice of should be appropriate for achieving differentiation between closely evolving clusters. Then, we calculate the fraction of elements that belong to some cluster at time as [16]where for and for . We refer to as the asymptotic time average of . Then, a clustered state in the system can be characterized by the value , while an incoherent state in the system corresponds to . The values characterize a chimera state, where is the minimum cluster size to be taken into consideration.

A synchronization state corresponds to the presence of a single cluster of size and has also the value . To distinguish a synchronization state from a multicluster state, we calculate the asymptotic time average aswhere is the number of discarded transients, is a sufficiently large time, and is the instantaneous standard deviation of the distribution of state variables defined bywhere

Statistically, a synchronization state is characterized by the values and , while a cluster state corresponds to and . Chimera states are characterized by and , and desynchronization is described by , . In this paper, we set and .

Note that, in systems with local or long-range interactions where there is a natural spatial ordering, the synchronized and desynchronized domains for chimera states are localized in space. In contrast, globally coupled systems lack the notion of spatial order. Thus, the chimera and cluster states in our system are characterized in terms of the statistical quantities and , not by the spatial location of the synchronized and the desynchronized domains.

3. Results and Discussion

3.1. Chimeras and Clusters in Globally Coupled Robust-Chaos Maps

Let us consider a network of globally coupled maps described by the equations [14]where the global interaction function is the mean field of the system, and

As local dynamics exhibiting robust chaos, we consider the following smooth, unimodal map defined on the interval [54]:which is chaotic with no periodic windows on the parameter interval . On this interval, the Lyapunov exponent of map equation (9) has the constant value . The bifurcation diagram of map equation (9) in Figure 1 shows the absence of periodicity in the interval .

Figure 2 shows the asymptotic temporal evolution of the states of system equations (7) and (9), for different values of parameters. Since the system is globally coupled, there is no natural spatial ordering. For visualization purposes, the indexes are ordered at time such that if and kept fixed afterwards. The values of the states are represented by distinct color coding; two elements share the same color if . A desynchronized state is displayed in Figure 2(a) and a complete synchronization state occurs in Figure 2(d), while a chimera state and a two-cluster state are visualized in Figures 2(b) and 2(c), respectively.

Figure 3 shows the collective states arising in system equations (7) and (9) on the space of parameters , characterized through the quantities and . Labels indicate the regions where these behaviors occur: CS: complete synchronization; C: cluster states; Q: chimera states, and D: desynchronization.

The linear stability analysis for the complete synchronization state in globally coupled system equation (7) shows that this state is stable if the following condition is satisfied [14]:where is the Lyapunov exponent for the local map . For map equation (9), we obtain that the completely synchronized state is stable for , which agrees with the numerical characterization for this state performed in Figure 3. Figure 3 reveals that both cluster and chimera states can arise in globally coupled map networks for appropriate values of parameters, even when the individual maps lack periodic windows. Clusters and chimera state regions occur adjacent to each other for an intermediate range of values of the coupling parameter on the phase diagram of Figure 3. In fact, chimeras and clusters are closely related collective states in systems subject to global interactions [30]. Chimera states appear to mediate between dynamical clustering and incoherence.

Multicluster chimera states are also possible in systems of globally coupled robust-chaos maps. As an illustration, consider the smooth unimodal map [55]defined on the interval for parameter values . Figure 4 shows the bifurcation diagram of the iterates of map equation (11) as a function of the parameter . The dynamics of the map displays robust chaos with no periodic windows for . The Lyapunov exponent is [55].

Figure 5 shows the temporal evolution of the states of globally coupled system equation (7) with local map equation (11), for different values of parameters. A chimera state with multiple clusters occurs in Figure 5(a), while a two-cluster state is shown in Figure 5(b). Multichimera states or multiheaded chimeras (coexistence of multiple localized domains of incoherence and coherence) have been reported in systems with long-range interactions [56]. However, those states are not equivalent to a chimera state with multiple clusters in a globally coupled system, such as Figure 5(a), where there is no notion of locality.

3.2. Dynamics of Clusters and Chimera States with Global Interactions

Consider a chimera state consisting of clusters and a desynchronized subset in the system of globally coupled map equation (1). The dynamics of this state can be described by the equations

The mean field equation (8) in a chimera state can be expressed as the sum of two contributionswhere

The term is the contribution to the mean field corresponding to elements belonging to clusters, whereas is the average of the states of the elements belonging to the incoherent subset.

Figure 6 shows the temporal behavior of both contributions and in a chimera state for globally coupled autonomous system equations (7) and (9). The time evolution of the cluster contribution is chaotic, similar to that of local map equation (9), but has a smaller amplitude. In general, the form of can be approximated as , where represents a modulation factor reflecting the partition into several clusters. On the other hand, Figure 6 reveals that the time series of fluctuates about a mean value with a small dispersion, corresponding to the superposition of the dynamics of many incoherent chaotic elements. Thus, for the given parameter values, the incoherent contribution to the mean field for a large system size can be expressed approximately as a constant, i.e., .

The dynamics of globally coupled system equation (1), where each map is subject to a feedback field , can be compared to that of a replica system of maps subject to a global external drive in the form

It has been shown that an analogy between autonomous system equation (1) and driven system equation (15) can be established when the time evolution of the field is identical to that of the function [9]. Then, the drive-response dynamics at the local level in both systems is similar, and therefore their corresponding emerging collective states can be equivalent for some appropriate parameter values and initial conditions. In particular, chimera or cluster states in system equation (15) should be induced by an external drive function of the form , with , constants, that imitates the mean field . The realization of these states depends on the parameters and of the drive and on the coupling strength .

Figure 7 shows the temporal evolution of the states of driven system equation (15) with local map equation (9), for some values of parameters. A chimera state with a single cluster takes place in Figure 7(a) for parameter values , where chimera states also occur in the autonomous system equations (7) and (9), as seen in the corresponding phase diagram of Figure 3. Figure 7(b) shows a two-cluster state for values located in the region corresponding to clustered states in Figure 3. The dynamics of the driven system equations (15) shows multistability; depending on initial conditions, chimeras with different partitions may be induced for given parameters values in the region labeled Q in Figure 3. Similarly, different initial conditions may produce cluster states with different partition sizes for fixed parameter values in region C of Figure 3.

System equation (15) can be considered as realizations for different initial conditions of a single driven map

Analogously, each local map in globally coupled system equation (7) can be seen as subject to a field that eventually induces a collective state. Clustering in globally coupled systems of identical elements has been attributed to the existence of periodic windows in the local dynamics [10]. On the other hand, clustering is considered a prerequisite for the occurrence of chimera states in globally coupled systems [31]. Thus, to elucidate the origin of clusters and chimeras in system equation (7) with local robust chaos, one can explore the response dynamics of the single driven map equation (16) with a function of the form and having robust chaos. Then, if periodic windows are induced by the drive on a single map, one may expect that clusters and chimeras should arise in a globally coupled system of those maps.

Even a trivial function can modify the dynamics of a driven robust-chaos map in equation (16) to produce periodic windows. Figure 8(a) shows the bifurcation diagram of in equation (16) versus for the map given by equation (9) with , which is equivalent to a rescaling of . Periodic windows typical of unimodal maps appear in the rescaled map . In general, the driven map equation (16) represents a rescaling of the robust-chaos map that acquires periodic windows. Similarly, the periodic cluster states arising in the globally coupled system equations (7) and (9) are a consequence of the windows of periodicity induced locally by the mean field , in analogy to the periodic windows created by an external drive acting on a single map equation (9). Different initial conditions may lead to different out-of-phase orbits with diverse partitions that appear as clusters in the globally coupled system. A synchronization state in system equations (7) and (9) can be associated to the fixed point interval of the bifurcation diagram of Figure 8(a), while a desynchronization state in the globally coupled system is a manifestation of a chaotic regime as seen in Figure 8(a). Nontrivial forms of the driving function can give rise to multistable behavior besides periodic windows. For example, we have verified that a drive such as in equation (16) induces a region of bistability between chaotic attractors that expresses as chimera states in the associated globally coupled system equations (7) and (9).

These results suggest that the emergence of cluster and chimera states in a globally coupled system of robust-chaos maps can be inferred from the occurrence of periodic windows in the response dynamics of a single map subject to an appropriate drive, as a function of parameters. Figure 8(b) shows the corresponding bifurcation diagram of versus for the map given by equation (11) which also has robust chaos. Again, we see the emergence of periodic windows as the coupling parameter is varied. A globally coupled system of these maps also shows clusters and chimera states, as illustrated in Figure 5. Figure 8(c) presents the bifurcation diagram of versus for the logarithmic map , which possesses robust chaos on the parameter interval and its dynamics is unbounded [2]. In contrast to Figures 8(a)and 8(b), no periodic windows appear on the dynamics of the driven map equation (16) which remains unbounded; only chaotic orbits and a fixed point attractor appear. As a consequence, clusters and chimera states should not be expected in a globally coupled system of logarithmic maps. In fact, only synchronization and nontrivial collective behavior have been observed in such a system [57].

4. Conclusions

Networks of globally coupled identical oscillators are among the simplest symmetric spatiotemporal systems that can display clustering and chimera behavior. Previous works have conjectured that these phenomena cannot occur when the local oscillators possess robust-chaos attractors [913]. We have shown that the presence of global interactions can indeed allow for emergence of both cluster and chimera states in systems of coupled robust-chaos maps. Chimeras appear as partially ordered states between synchronization or clustering and incoherent behavior. We have found that chimera states are associated to the formation of clusters in these systems, a feature that has been observed in other globally coupled systems [31].

The existence of intrinsic periodic windows in the dynamics of local oscillators, such as in logistic maps, is not essential for the emergence of clusters with periodic behavior in a globally coupled system of those oscillators. Windows of periodicity and multistability can be induced in the dynamical response of a robust-chaos map subject to an appropriate external forcing. Because of the analogy between a single driven map and the local dynamics of a globally coupled map system, the global interaction field can also induce periodic windows and multistability on local robust-chaos maps. Those are the essential ingredients for the occurrence of cluster and chimera states in globally coupled systems. Since clustering is a prerequisite for chimeras [31], a single driven robust-chaos map that develops periodic windows on some range of parameters allows us to infer that a globally coupled system of such maps shall also exhibit cluster and chimera states on some range of parameters. Conversely, a robust-chaos map, such as the logarithmic or another singular map, which does not give rise to periodic windows when subject to a drive, implies that a system of globally coupled logarithmic or singular maps does not show clusters nor chimera states.

Further extensions of this work include the investigation of chimera states in networks of globally coupled continuous-time dynamical systems possessing robust chaos or hyperbolic chaotic attractors, the study of interacting populations of robust-chaos elements, and the role of the range of interaction in a network of dynamical robust-chaos units.

Data Availability

The programming codes used in this article are available upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This study was supported by Corporación Ecuatoriana para el Desarrollo de la Investigación y Academia (CEDIA) through project CEPRA-XIII-2019 “Sistemas Complejos.”