Coexistence of synchrony and incoherence in oscillatory media under nonlinear global coupling

Collective dynamics of coupled Lorenz oscillators near the Hopf boundary: Intermittency and chimera states

We study collective dynamics of networks of mutually coupled identical Lorenz oscillators near a subcritical Hopf bifurcation. Such systems exhibit induced multistable behavior with interesting spatiotemporal dynamics including synchronization, desynchronization, and chimera states. For analysis, we first consider a ring topology with nearest-neighbor coupling and find that the system may exhibit intermittent behavior due to the complex basin structures and dynamical frustration, where temporal dynamics of the oscillators in the ensemble switches between different attractors. Consequently, different oscillators may show a dynamics that is intermittently synchronized (or desynchronized), giving rise to intermittent chimera states. The behavior of the intermittent laminar phases is characterized by the characteristic time spent in the synchronization manifold, which decays as a power law. Such intermittent dynamics is quite general and is also observed in an ensemble of a large number of oscillators arranged in variety of network topologies including nonlocal, scale-free, random, and small-world networks.

Nonequilibrium thermodynamic signatures of collective dynamical states around chimera in a chemical reaction network

Different dynamical states ranging from coherent, incoherent to chimera, multichimera, and related transitions are addressed in a globally coupled nonlinear continuum chemical oscillator system by implementing a modified complex Ginzburg-Landau equation. Besides dynamical identifications of observed states using standard qualitative metrics, we systematically acquire nonequilibrium thermodynamic characterizations of these states obtained via coupling parameters. The nonconservative work profiles in collective dynamics qualitatively reflect the time-integrated concentration of the activator, and the majority of the nonconservative work contributes to the entropy production over the spatial dimension. It is illustrated that the evolution of spatial entropy production and semigrand Gibbs free-energy profiles associated with each state are connected yet completely out of phase, and these thermodynamic signatures are extensively elaborated to shed light on the exclusiveness and similarities of these states. Moreover, a relationship between the proper nonequilibrium thermodynamic potential and the variance of activator concentration is established by exhibiting both quantitative and qualitative similarities between a Fano factor like entity, derived from the activator concentration, and the Kullback-Leibler divergence associated with the transition from a nonequilibrium homogeneous state to an inhomogeneous state. Quantifying the thermodynamic costs for collective dynamical states would aid in efficiently controlling, manipulating, and sustaining such states to explore the real-world relevance and applications of these states.

Chimeras in globally coupled oscillators: A review

The surprising phenomenon of chimera in an ensemble of identical oscillators is no more strange behavior of network dynamics and reality. By this time, this symmetry breaking self-organized collective dynamics has been established in many networks, a ring of non-locally coupled oscillators, globally coupled networks, a three-dimensional network, and multi-layer networks. A variety of coupling and dynamical models in addition to the phase oscillators has been used for a successful observation of chimera patterns. Experimental verification has also been done using metronomes, pendula, chemical, and opto-electronic systems. The phenomenon has also been shown to appear in small networks, and hence, it is not size-dependent. We present here a brief review of the origin of chimera patterns restricting our discussions to networks of globally coupled identical oscillators only. The history of chimeras in globally coupled oscillators is older than what has been reported in nonlocally coupled phase oscillators much later. We elaborate the story of the origin of chimeras in globally coupled oscillators in a chronological order, within our limitations, and with brief descriptions of the significant contributions, including our personal experiences. We first introduce chimeras in non-locally coupled and other network configurations, in general, and then discuss about globally coupled networks in more detail.

Chimera resonance in networks of chaotic maps

We explore numerically the impact of additive Gaussian noise on the spatiotemporal dynamics of ring networks of nonlocally coupled chaotic maps. The local dynamics of network nodes is described by the logistic map, the Ricker map, and the Henon map. 2D distributions of the probability of observing chimera states are constructed in terms of the coupling strength and the noise intensity and for several choices of the local dynamics parameters. It is shown that the coupling strength range can be the widest at a certain optimum noise level at which chimera states are observed with a high probability for a large number of different realizations of randomly distributed initial conditions and noise sources. This phenomenon demonstrates a constructive role of noise in analogy with the effects of stochastic and coherence resonance and may be referred to as chimera resonance.

Transient chimera states emerging from dynamical trapping in chaotic saddles

Nonlinear systems possessing nonattracting chaotic sets, such as chaotic saddles, embedded in their state space may oscillate chaotically for a transient time before eventually transitioning into some stable attractor. We show that these systems, when networked with nonlocal coupling in a ring, are capable of forming chimera states, in which one subset of the units oscillates periodically in a synchronized state forming the coherent domain, while the complementary subset oscillates chaotically in the neighborhood of the chaotic saddle constituting the incoherent domain. We find two distinct transient chimera states distinguished by their abrupt or gradual termination. We analyze the lifetime of both chimera states, unraveling their dependence on coupling range and size. We find an optimal value for the coupling range yielding the longest lifetime for the chimera states. Moreover, we implement transversal stability analysis to demonstrate that the synchronized state is asymptotically stable for network configurations studied here.

Detecting partial synchrony in a complex oscillatory network using pseudovortices

Partial synchronization is an important dynamical process of coupled oscillators on various natural and artificial networks, which can remain undetected due to the system complexity. With an analogy between pairwise asynchrony of oscillators and topological defects, i.e., vortices, in the two-dimensional XY model, we propose a robust and data-driven method to identify the partial synchronization on complex networks. The proposed method is based on an integer matrix whose element is pseudovorticity that discretely quantifies asynchronous phase dynamics in every two oscillators, which results in graphical and entropic representations of partial synchrony. As a first trial, we apply our method to 200 FitzHugh-Nagumo neurons on a complex small-world network. Partially synchronized chimera states are revealed by discriminating synchronized states even with phase lags. Such phase lags also appear in partial synchronization in chimera states. Our topological, graphical, and entropic method is implemented solely with measurable phase dynamics data, which will lead to a straightforward application to general oscillatory networks including neural networks in the brain.

Short-lived chimera states

In the classic Kuramoto system of coupled two-dimensional rotators, chimera states characterized by the coexistence of synchronous and asynchronous groups of oscillators are long-lived because the average lifetime of these states increases exponentially with the system size. Recently, it was discovered that, when the rotators in the Kuramoto model are three-dimensional, the chimera states become short-lived in the sense that their lifetime scales with only the logarithm of the dimension-augmenting perturbation. We introduce transverse-stability analysis to understand the short-lived chimera states. In particular, on the unit sphere representing three-dimensional (3D) rotations, the long-lived chimera states in the classic Kuramoto system occur on the equator, to which latitudinal perturbations that make the rotations 3D are transverse. We demonstrate that the largest transverse Lyapunov exponent calculated with respect to these long-lived chimera states is typically positive, making them short-lived. The transverse-stability analysis turns the previous numerical scaling law of the transient lifetime into an exact formula: the "free" proportional constant in the original scaling law can now be precisely determined in terms of the largest transverse Lyapunov exponent. Our analysis reinforces the speculation that in physical systems, chimera states can be short-lived as they are vulnerable to any perturbations that have a component transverse to the invariant subspace in which they live.

Chimera patterns in conservative Hamiltonian systems and Bose-Einstein condensates of ultracold atoms

Experimental realizations of chimera patterns, characterized by coexisting regions of phase coherence and incoherence, have so far been achieved for non-conservative systems with dissipation and exclusively in classical settings. The possibility of observing chimera patterns in quantum systems has rarely been studied and it remains an open question if chimera patterns can exist in closed, or conservative quantum systems. Here, we tackle these challenges by first proposing a conservative Hamiltonian system with nonlocal hopping, where the energy is well-defined and conserved. We show explicitly that such a system can exhibit chimera patterns. Then we propose a physical mechanism for the nonlocal hopping by using an additional mediating channel. This leads us to propose a possible experimentally realizable quantum system based on a two-component Bose-Einstein condensate (BEC) with a spin-dependent optical lattice, where an untrapped component serves as the matter-wave mediating field. In this BEC system, nonlocal spatial hopping over tens of lattice sites can be achieved and simulations suggest that chimera patterns should be observable in certain parameter regimes.

Transition from chimera/solitary states to traveling waves

We study numerically the spatiotemporal dynamics in a ring network of nonlocally coupled nonlinear oscillators, each represented by a two-dimensional discrete-time model of the classical van der Pol oscillator. It is shown that the discretized oscillator exhibits richer behavior, combining the peculiarities of both the original system and its own dynamics. Moreover, a large variety of spatiotemporal structures is observed in the network of discrete van der Pol oscillators when the discretization parameter and the coupling strength are varied. Regimes, such as the coexistence of a multichimera state/a traveling wave and a solitary state are revealed for the first time and are studied in detail. It is established that the majority of the observed chimera/solitary states, including the newly found ones, are transient toward a purely traveling wave mode. The peculiarities of the transition process and the lifetime (transient duration) of the chimera structures and the solitary state are analyzed depending on the system parameters, the observation time, initial conditions, and the influence of external noise.

From Turing patterns to chimera states in the 2D Brusselator model

The Brusselator has been used as a prototype model for autocatalytic reactions and, in particular, for the Belousov-Zhabotinsky reaction. When coupled at the diffusive limit, the Brusselator undergoes a Turing bifurcation resulting in the formation of classical Turing patterns, such as spots, stripes, and spirals in two spatial dimensions. In the present study, we use generic nonlocally coupled Brusselators and show that in the limit of the coupling range R→1 (diffusive limit), the classical Turing patterns are recovered, while for intermediate coupling ranges and appropriate parameter values, chimera states are produced. This study demonstrates how the parameters of a typical nonlinear oscillator can be tuned so that the coupled system passes from spatially stable Turing structures to dynamical spatiotemporal chimera states.

Synchronization of Belousov-Zhabotinsky oscillators with electrochemical coupling in a spontaneous process

A passive electrochemical coupling approach is proposed to induce spontaneous synchronization between chemical oscillators. The coupling exploits the potential difference between a catalyst redox couple in the Belousov-Zhabotinsky (BZ) reaction, without external feedback, to induce surface reactions that impact the kinetics of the bulk system. The effect of coupling in BZ oscillators under batch condition is characterized using phase synchronization measures. Although the frequency of the oscillators decreases nonlinearly over time, by a factor of 2 or more within 100 cycles, the coupling is strong enough to maintain synchronization. In such a highly drifting system, the Gibbs-Shannon entropy of the cyclic phase difference distribution can be used to quantify the coupling effect. We extend the Oregonator BZ model to account for the drifting natural frequencies in batch condition and for electrochemical coupling, and numerical simulations of the effect of acid concentration on synchronization patterns are in agreement with the experiments. Because of the passive nature of coupling, the proposed coupling scheme can open avenues for designing pattern recognition and neuromorphic computation systems using chemical reactions in a spontaneous process.

A global bifurcation organizing rhythmic activity in a coupled network

We study a system of coupled phase oscillators near a saddle-node on invariant circle bifurcation and driven by random intrinsic frequencies. Under the variation of control parameters, the system undergoes a phase transition changing the qualitative properties of collective dynamics. Using Ott-Antonsen reduction and geometric techniques for ordinary differential equations, we identify heteroclinic bifurcation in a family of vector fields on a cylinder, which explains the change in collective dynamics. Specifically, we show that heteroclinic bifurcation separates two topologically distinct families of limit cycles: contractible limit cycles before bifurcation from noncontractibile ones after bifurcation. Both families are stable for the model at hand.

Disordered quenching in arrays of coupled Bautin oscillators

In this work, we study the phenomenon of disordered quenching in arrays of coupled Bautin oscillators, which are the normal form for bifurcation in the vicinity of the equilibrium point when the first Lyapunov coefficient vanishes and the second one is nonzero. For particular parameter values, the Bautin oscillator is in a bistable regime with two attractors-the equilibrium and the limit cycle-whose basins are separated by the unstable limit cycle. We consider arrays of coupled Bautin oscillators and study how they become quenched with increasing coupling strength. We analytically show the existence and stability of the dynamical regimes with amplitude disorder in a ring of coupled Bautin oscillators with identical natural frequencies. Next, we numerically provide evidence that disordered oscillation quenching holds for rings as well as chains with nonidentical natural frequencies and study the characteristics of this effect.

Higher-order interactions promote chimera states

Since the discovery of chimera states, the presence of a nonzero phase lag parameter turns out to be an essential attribute for the emergence of chimeras in a nonlocally coupled identical Kuramoto phase oscillators' network with pairwise interactions. In this Letter, we report the emergence of chimeras without phase lag in a nonlocally coupled identical Kuramoto network owing to the introduction of nonpairwise interactions. The influence of added nonlinearity in the coupled system dynamics in the form of simplicial complexes mitigates the requisite of a nonzero phase lag parameter for the emergence of chimera states. Chimera states stimulated by the reciprocity of the pairwise and nonpairwise interaction strengths and their multistable nature are characterized with appropriate measures and are demonstrated in the parameter spaces.

Nonequilibrium thermodynamic characterization of chimeras in a continuum chemical oscillator system

The emergence of the chimera state as the counterintuitive spatial coexistence of synchronous and asynchronous regimes is addressed here in a continuum chemical oscillator system by implementing a relevant complex Ginzburg-Landau equation with global coupling. This study systematically acquires and characterizes the evolution of nonequilibrium thermodynamic entities corresponding to the chimera state. The temporal evolution of the entropy production rate exhibits a beat pattern with a series of equidistant spectral lines in the frequency domain. Symmetric profiles associated with the incoherent regime appear in descriptions of the dynamics and thermodynamics of the chimera. It is shown that identifying the semigrand Gibbs free energy of the state as the Gabor elementary function can reveal the guiding role of the information uncertainty principle in shaping the chimera energetics.

A new class of chimeras in locally coupled oscillators with small-amplitude, high-frequency asynchrony and large-amplitude, low-frequency synchrony

Chimeras are surprising yet important states in which domains of decoherent (asynchronous) and coherent (synchronous) oscillations co-exist. In this article, we report on the discovery of a new class of chimeras, called mixed-amplitude chimera states, in which the structures, amplitudes, and frequencies of the oscillations differ substantially in the decoherent and coherent regions. These mixed-amplitude chimeras exhibit domains of decoherent small-amplitude oscillations (phase waves) coexisting with domains of stable and coherent large-amplitude or mixed-mode oscillations (MMOs). They are observed in a prototypical bistable partial differential equation with oscillatory dynamics, spatially homogeneous kinetics, and purely local, isotropic diffusion. They are observed in parameter regimes immediately adjacent to regimes in which common large-amplitude solutions exist, such as trigger waves, spatially homogeneous MMOs, and sharp-interface solutions. Also, key singularities, folded nodes, and folded saddles arising commonly in multi-scale, bistable systems play important roles, and these have not previously been studied in systems with chimeras. The discovery of these mixed-amplitude chimeras is an important advance for understanding some processes in neuroscience, pattern formation, and physics, which involve both small-amplitude and large-amplitude oscillations. It may also be of use for understanding some aspects of electroencephalogram recordings from animals that exhibit unihemispheric slow-wave sleep.

Attracting Poisson chimeras in two-population networks

Chimera states, i.e., dynamical states composed of coexisting synchronous and asynchronous oscillations, have been reported to exist in diverse topologies of oscillators in simulations and experiments. Two-population networks with distinct intra- and inter-population coupling have served as simple model systems for chimera states since they possess an invariant synchronized manifold in contrast to networks on a spatial structure. Here, we study dynamical and spectral properties of finite-sized chimeras on two-population networks. First, we elucidate how the Kuramoto order parameter of the finite-sized globally coupled two-population network of phase oscillators is connected to that of the continuum limit. These findings suggest that it is suitable to classify the chimera states according to their order parameter dynamics, and therefore, we define Poisson and non-Poisson chimera states. We then perform a Lyapunov analysis of these two types of chimera states, which yields insight into the full stability properties of the chimera trajectories as well as of collective modes. In particular, our analysis also confirms that Poisson chimeras are neutrally stable. We then introduce two types of "perturbation" that act as small heterogeneities and render Poisson chimeras attracting: A topological variation via the simplest nonlocal intra-population coupling that keeps the network symmetries and the allowance of amplitude variations in the globally coupled two-population network; i.e., we replace the phase oscillators by Stuart-Landau oscillators. The Lyapunov spectral properties of chimera states in the two modified networks are investigated, exploiting an approach based on network symmetry-induced cluster pattern dynamics of the finite-size network.

Between synchrony and turbulence: intricate hierarchies of coexistence patterns

Coupled oscillators, even identical ones, display a wide range of behaviours, among them synchrony and incoherence. The 2002 discovery of so-called chimera states, states of coexisting synchronized and unsynchronized oscillators, provided a possible link between the two and definitely showed that different parts of the same ensemble can sustain qualitatively different forms of motion. Here, we demonstrate that globally coupled identical oscillators can express a range of coexistence patterns more comprehensive than chimeras. A hierarchy of such states evolves from the fully synchronized solution in a series of cluster-splittings. At the far end of this hierarchy, the states further collide with their own mirror-images in phase space - rendering the motion chaotic, destroying some of the clusters and thereby producing even more intricate coexistence patterns. A sequence of such attractor collisions can ultimately lead to full incoherence of only single asynchronous oscillators. Chimera states, with one large synchronized cluster and else only single oscillators, are found to be just one step in this transition from low- to high-dimensional dynamics.

Birhythmicity, intrinsic entrainment, and minimal chimeras in an electrochemical experiment

The coexistence of limit cycles in a phase space, so called birhythmicity, is a phenomenon known to exist in many systems in various disciplines. Yet, detailed experimental investigations are rare, as are studies on the interaction between birhythmic components. In this article, we present experimental evidence for the existence of birhythmicity during the anodic electrodissolution of Si in a fluoride-containing electrolyte using weakly illuminated n-type Si electrodes. Moreover, we demonstrate several types of interaction between the coexisting limit cycles, in part resulting in peculiar dynamics. The two limit cycles exhibit vastly different sensitivities with respect to a small perturbation of the electrode potential, rendering the coupling essentially unidirectional. A manifestation of this is an asymmetric 1:2 intrinsic entrainment of the coexisting limit cycles on an individual uniformly oscillating electrode. In this state, the phase-space structure mediates the locking of one of the oscillators to the other one across the separatrix. Furthermore, the transition scenarios from one limit cycle to the other one at the borders of the birhythmicity go along with different types of spatial symmetry breaking. Finally, the master-slave type coupling promotes two (within the experimental limits) identical electrodes initialized on the different limit cycles to adopt states of different complexity: one of the electrodes exhibits irregular, most likely chaotic, motion, while the other one exhibits period-1 oscillations. The coexistence of coherence and incoherence is the characteristic property of a chimera state, the two coupled electrodes constituting an experimental example of a smallest chimera state in a minimal network configuration.

Connecting minimal chimeras and fully asymmetric chaotic attractors through equivariant pitchfork bifurcations

Highly symmetric networks can exhibit partly symmetry-broken states, including clusters and chimera states, i.e., states of coexisting synchronized and unsynchronized elements. We address the S_{4} permutation symmetry of four globally coupled Stuart-Landau oscillators and uncover an interconnected web of solutions with different symmetries. Among these are chaotic 2-1-1 minimal chimeras that arise from 2-1-1 periodic solutions in a period-doubling cascade, as well as fully asymmetric chaotic states arising similarly from periodic 1-1-1-1 solutions. A backbone of equivariant pitchfork bifurcations mediate between the two cascades, culminating in equivariant pitchforks of chaotic attractors.

Amplitude-mediated spiral chimera pattern in a nonlinear reaction-diffusion system

Formation of diverse patterns in spatially extended reaction-diffusion systems is an important aspect of study that is pertinent to many chemical and biological processes. Of special interest is the peculiar phenomenon of chimera state having spatial coexistence of coherent and incoherent dynamics in a system of identically interacting individuals. In the present article, we report the emergence of various collective dynamical patterns while considering a system of prey-predator dynamics in the presence of a two-dimensional diffusive environment. Particularly, we explore the observance of four distinct categories of spatial arrangements among the species, namely, spiral wave, spiral chimera, completely synchronized oscillations, and oscillation death states in a broad region of the diffusion-driven parameter space. Emergence of amplitude-mediated spiral chimera states displaying drifted amplitudes and phases in the incoherent subpopulation is detected for parameter values beyond both Turing and Hopf bifurcations. Transition scenarios among all these distinguishable patterns are numerically demonstrated for a wide range of the diffusion coefficients which reveal that the chimera states arise during the transition from oscillatory to steady-state dynamics. Furthermore, we characterize the occurrence of each of the recognizable patterns by estimating the strength of incoherent subpopulations in the two-dimensional space.

Dynamics in two interacting subpopulations of nonidentical phase oscillators

Chimera states refer to the dynamical states in which the inherent symmetry of the system is broken. The system composed of two interacting identical subpopulations of phase oscillators provides a platform to study chimera states. In this system, different types of chimera states have been identified and the transitions between them have been investigated. However, the parameter space is not fully explored in this system. In this work, we study a system comprised of two interacting subpopulations of nonidentical phase oscillators. Through numerical simulations and theoretical analyses, we find three symmetry-reserving states, including incoherent state, in-phase synchronous state, and antiphase synchronous state, and three types of symmetry-breaking states, including in-phase chimera states, antiphase chimera states, and weak chimera states. The stability diagrams of these dynamical states are explored on different parameter planes and transition scenarios amongst these states are investigated. We find that the weak chimera states act as the bridge between in-phase and antiphase chimera states. We also observe the existence of a period-two chimera state, chaotic chimera state, and drifting chimera states.

Self-Organized Multifrequency Clusters in an Oscillating Electrochemical System with Strong Nonlinear Coupling

We study the spatiotemporal dynamics of the oscillatory photoelectrodissolution of n-type Si in a fluoride-containing electrolyte under anodic potentials using in situ ellipsometric imaging. When lowering the illumination intensity stepwise, we successively observe uniform oscillations, modulated amplitude clusters, and the coexistence of multifrequency clusters, i.e., regions with different frequencies, with a stationary domain. We argue that the multifrequency clusters emerge due to an adaptive, nonlinear, and nonlocal coupling, similar to those found in the context of neural dynamics.

Mechanism for Strong Chimeras

Chimera states have attracted significant attention as symmetry-broken states exhibiting the unexpected coexistence of coherence and incoherence. Despite the valuable insights gained from analyzing specific systems, an understanding of the general physical mechanism underlying the emergence of chimeras is still lacking. Here, we show that many stable chimeras arise because coherence in part of the system is sustained by incoherence in the rest of the system. This mechanism may be regarded as a deterministic analog of noise-induced synchronization and is shown to underlie the emergence of strong chimeras. These are chimera states whose coherent domain is formed by identically synchronized oscillators. Recognizing this mechanism offers a new meaning to the interpretation that chimeras are a natural link between coherence and incoherence.

Dynamics of multilayer networks with amplification

We study the dynamics of a multilayer network of chaotic oscillators subject to amplification. Previous studies have proven that multilayer networks present phenomena such as synchronization, cluster, and chimera states. Here, we consider a network with two layers of Rössler chaotic oscillators as well as applications to multilayer networks of the chaotic jerk and Liénard oscillators. Intra-layer coupling is considered to be all to all in the case of Rössler oscillators, a ring for jerk oscillators and global mean field coupling in the case of Liénard, inter-layer coupling is unidirectional in all these three cases. The second layer has an amplification coefficient. An in-depth study on the case of a network of Rössler oscillators using a master stability function and order parameter leads to several phenomena such as complete synchronization, generalized, cluster, and phase synchronization with amplification. For the case of Rössler oscillators, we note that there are also certain values of coupling parameters and amplification where the synchronization does not exist or the synchronization can exist but without amplification. Using other systems with different topologies, we obtain some interesting results such as chimera state with amplification, cluster state with amplification, and complete synchronization with amplification.

Quasi phase reduction of all-to-all strongly coupled λ-ω oscillators near incoherent states

The dynamics of an ensemble of N weakly coupled limit-cycle oscillators can be captured by their N phases using standard phase reduction techniques. However, it is a phenomenological fact that all-to-all strongly coupled limit-cycle oscillators may behave as "quasiphase oscillators," evidencing the need of novel reduction strategies. We introduce, here, quasi phase reduction (QPR), a scheme suited for identical oscillators with polar symmetry (λ-ω systems). By applying QPR, we achieve a reduction to N+2 degrees of freedom: N phase oscillators interacting through one independent complex variable. This "quasi phase model" is asymptotically valid in the neighborhood of incoherent states, irrespective of the coupling strength. The effectiveness of QPR is illustrated in a particular case, an ensemble of Stuart-Landau oscillators, obtaining exact stability boundaries of uniform and nonuniform incoherent states for a variety of couplings. An extension of QPR beyond the neighborhood of incoherence is also explored. Finally, a general QPR model with N+2M degrees of freedom is obtained for coupling through the first M harmonics.

Remote pacemaker control of chimera states in multilayer networks of neurons

Networks of coupled nonlinear oscillators allow for the formation of nontrivial partially synchronized spatiotemporal patterns, such as chimera states, in which there are coexisting coherent (synchronized) and incoherent (desynchronized) domains. These complementary domains form spontaneously, and it is impossible to predict where the synchronized group will be positioned within the network. Therefore, possible ways to control the spatial position of the coherent and incoherent groups forming the chimera states are of high current interest. In this work we investigate how to control chimera patterns in multiplex networks of FitzHugh-Nagumo neurons, and in particular we want to prove that it is possible to remotely control chimera states exploiting the multiplex structure. We introduce a pacemaker oscillator within the network: this is an oscillator that does not receive input from the rest of the network but is sending out information to its neighbors. The pacemakers can be positioned in one or both layers. Their presence breaks the spatial symmetry of the layer in which they are introduced and allows us to control the position of the incoherent domain. We demonstrate how the remote control is possible for both uni- and bidirectional coupling between the layers. Furthermore we show which are the limitations of our control mechanisms when it is generalized from single-layer to multilayer networks.

Structural anomalies in brain networks induce dynamical pacemaker effects

Dynamical effects on healthy brains and brains affected by tumor are investigated via numerical simulations. The brains are modeled as multilayer networks consisting of neuronal oscillators whose connectivities are extracted from Magnetic Resonance Imaging (MRI) data. The numerical results demonstrate that the healthy brain presents chimera-like states where regions with high white matter concentrations in the direction connecting the two hemispheres act as the coherent domain, while the rest of the brain presents incoherent oscillations. To the contrary, in brains with destructed structures, traveling waves are produced initiated at the region where the tumor is located. These areas act as the pacemaker of the waves sweeping across the brain. The numerical simulations are performed using two neuronal models: (a) the FitzHugh-Nagumo model and (b) the leaky integrate-and-fire model. Both models give consistent results regarding the chimera-like oscillations in healthy brains and the pacemaker effect in the tumorous brains. These results are considered a starting point for further investigation in the detection of tumors with small sizes before becoming discernible on MRI recordings as well as in tumor development and evolution.

Chimera states in coupled map lattices: Spatiotemporally intermittent behavior and an equivalent cellular automaton

We construct an equivalent cellular automaton (CA) for a system of globally coupled sine circle maps with two populations and distinct values for intergroup and intragroup coupling. The phase diagram of the system shows that the coupled map lattice can exhibit chimera states with synchronized and spatiotemporally intermittent subgroups after evolution from random initial conditions in some parameter regimes, as well as to other kinds of solutions in other parameter regimes. The CA constructed by us reflects the global nature and the two population structure of the coupled map lattice and is able to reproduce the phase diagram accurately. The CA depends only on the total number of laminar and burst sites and shows a transition from co-existing deterministic and probabilistic behavior in the chimera region to fully probabilistic behavior at the phase boundaries. This identifies the characteristic signature of the transition of a cellular automaton to a chimera state. We also construct an evolution equation for the average number of laminar/burst sites from the CA, analyze its behavior and solutions, and correlate these with the behavior seen for the coupled map lattice. Our CA and methods of analysis can have relevance in wider contexts.

Control of inter-layer synchronization by multiplexing noise

We study the synchronization of spatio-temporal patterns in a two-layer network of coupled chaotic maps, where each layer is represented by a nonlocally coupled ring. In particular, we focus on noisy inter-layer communication that we call multiplexing noise. We show that noisy modulation of inter-layer coupling strength has a significant impact on the dynamics of the network and specifically on the degree of synchronization of spatio-temporal patterns of interacting layers initially (in the absence of interaction) exhibiting chimera states. Our goal is to develop control strategies based on multiplexing noise for both identical and non-identical layers. We find that for the appropriate choice of intensity and frequency characteristics of parametric noise, complete or partial synchronization of the layers can be observed. Interestingly, for achieving inter-layer synchronization through multiplexing noise, it is crucial to have colored noise with intermediate spectral width. In the limit of white noise, the synchronization is destroyed. These results are the first step toward understanding the role of noisy inter-layer communication for the dynamics of multilayer networks.

A Brief Review of Chimera State in Empirical Brain Networks

Understanding the human brain and its functions has always been an interesting and challenging problem. Recently, a significant progress on this problem has been achieved on the aspect of chimera state where a coexistence of synchronized and unsynchronized states can be sustained in identical oscillators. This counterintuitive phenomenon is closely related to the unihemispheric sleep in some marine mammals and birds and has recently gotten a hot attention in neural systems, except the previous studies in non-neural systems such as phase oscillators. This review will briefly summarize the main results of chimera state in neuronal systems and pay special attention to the network of cerebral cortex, aiming to accelerate the study of chimera state in brain networks. Some outlooks are also discussed.

Chimera-like states induced by additional dynamic nonlocal wirings

We investigate the existence of chimera-like states in a small-world network of chaotically oscillating identical Rössler systems with an addition of randomly switching nonlocal links. By varying the small-world coupling strength, we observe no chimera-like state either in the absence of nonlocal wirings or with static nonlocal wirings. When we give an additional nonlocal wiring to randomly selected nodes and if we allow the random selection of nodes to change with time, we observe the onset of chimera-like states. Upon increasing the number of randomly selected nodes gradually, we find that the incoherent window keeps on shrinking, whereas the chimera-like window widens up. Moreover, the system attains a completely synchronized state comparatively sooner for a lower coupling strength. Also, we show that one can induce chimera-like states by a suitable choice of switching times, coupling strengths, and a number of nonlocal links. We extend the above-mentioned randomized injection of nonlocal wirings for the cases of globally coupled Rössler oscillators and a small-world network of coupled FitzHugh-Nagumo oscillators and obtain similar results.

Emergence of second coherent regions for breathing chimera states

Chimera states in one-dimensional nonlocally coupled phase oscillators are mostly assumed to be stationary, but breathing chimeras can occasionally appear, branching from the stationary chimeras via Hopf bifurcation. In this paper, we demonstrate two types of breathing chimeras: The type I breathing chimera looks the same as the stationary chimera at a glance, while the type II consists of multiple coherent regions with different average frequencies. Moreover, it is shown that the type I changes to the type II by increasing the breathing amplitude. Furthermore, we develop a self-consistent analysis of the local order parameter, which can be applied to breathing chimeras, and numerically demonstrate this analysis in the present system.

Symmetry breakings in two populations of oscillators coupled via diffusive environments: Chimera and heterosynchrony

We consider two diffusively coupled populations of identical oscillators, where the oscillators in each population are coupled with a common dynamic environment. Existence and stability of a variety of stationary states are analyzed on the basis of the Ott-Antonsen reduction method, which reveals that the chimera state occurs under the diffusive coupling scheme. Furthermore, we find an exotic symmetry-breaking behavior, the so-called the heterosynchronous state, in which each population exhibits in-phase coherence, while the order parameters of two populations rotate at different phase velocities. The chimera and heterosynchronous states emerge from bistabilities of distinct states for decoupled population and occur as a unique continuation for weak diffusive couplings. The heterosynchronous state is caused by an indirect coupling scheme via dynamic environments and could occur for a finite-size system as well, even for the system that consists of one oscillator per population.

Post-canard symmetry breaking and other exotic dynamic behaviors in identical coupled chemical oscillators

We analyze a model of two identical chemical oscillators coupled through diffusion of the slow variable. As a parameter is varied, a single oscillator undergoes a canard explosion-a transition from small amplitude, nearly harmonic oscillations to large-amplitude, relaxation oscillations over a very small parameter interval. In the coupled system, if the two oscillators have the same initial conditions, then the oscillators remain synchronized and exhibit the same canard behavior observed for the single oscillator. If the oscillators are separated initially, then in the region of the canard they display a variety of complex behaviors, including intermittent spiking, mixed-mode oscillation, and quasiperiodicity. Further variation of the parameter leads to a return to synchronized large-amplitude oscillation followed by a post-canard symmetry-breaking, in which one oscillator shows small-amplitude, complex behavior (mixed-mode oscillation, quasiperiodicity, chaos,...) while the other undergoes essentially periodic large amplitude behavior, resembling a master-slave scenario. We analyze the origins of this behavior by looking at several modified coupling schemes.

Robust partial synchronization of delay-coupled networks

Networks of coupled systems may exhibit a form of incomplete synchronization called partial synchronization or cluster synchronization, which refers to the situation where only some, but not all, systems exhibit synchronous behavior. Moreover, due to perturbations or uncertainties in the network, exact partial synchronization in the sense that the states of the systems within each cluster become identical, cannot be achieved. Instead, an approximate synchronization may be observed, where the states of the systems within each cluster converge up to some bound, and this bound tends to zero if (the size of) the perturbations tends to zero. In order to derive sufficient conditions for this robustified notion of synchronization, which we refer to as practical partial synchronization, first, we separate the synchronization error dynamics from the network dynamics and interpret them in terms of a nonautonomous system of delay differential equations with a bounded additive perturbation. Second, by assessing the practical stability of this error system, conditions for practical partial synchronization are derived and formulated in terms of linear matrix inequalities. In addition, an explicit relation between the size of perturbation and the bound of the synchronization error is provided.

Lyapunov spectra and collective modes of chimera states in globally coupled Stuart-Landau oscillators

Oscillatory systems with long-range or global coupling offer promising insight into the interplay between high-dimensional (or microscopic) chaotic motion and collective interaction patterns. Within this paper, we use Lyapunov analysis to investigate whether chimera states in globally coupled Stuart-Landau (SL) oscillators exhibit collective degrees of freedom. We compare two types of chimera states, which emerge in SL ensembles with linear and nonlinear global coupling, respectively, the latter introducing a constraint that conserves the oscillation of the mean. Lyapunov spectra reveal that for both chimera states the Lyapunov exponents split into several groups with different convergence properties in the limit of large system size. Furthermore, in both cases the Lyapunov dimension is found to scale extensively and the localization properties of covariant Lypunov vectors manifest the presence of collective Lyapunov modes. Here, however, we find qualitative differences between the two types of chimera states: Whereas the ones in the system under nonlinear global coupling exhibit only slow collective modes corresponding to Lyapunov exponents equal or close to zero, those which experience the linear mean-field coupling exhibit also faster collective modes associated with Lyapunov exponents with large positive or negative values. Furthermore, for the fastest collective mode we showed that it spreads across both synchonous and incoherent oscillators.

Blinking chimeras in globally coupled rotators

In globally coupled ensembles of identical oscillators so-called chimera states can be observed. The chimera state is a symmetry-broken regime, where a subset of oscillators forms a cluster, a synchronized population, while the rest of the system remains a collection of nonsynchronized, scattered units. We describe here a blinking chimera regime in an ensemble of seven globally coupled rotators (Kuramoto oscillators with inertia). It is characterized by a death-birth process, where a long-term stable cluster of four oscillators suddenly dissolves and is very quickly reborn with a new reshuffled configuration. We identify three different kinds of rare blinking events and give a quantitative characterization by applying stability analysis to the long-lived chaotic state and to the short-lived regular regimes that arise when the cluster dissolves.

Controlling chimera states via minimal coupling modification

We propose a method to control chimera states in a ring-shaped network of nonlocally coupled phase oscillators. This method acts exclusively on the network's connectivity. Using the idea of a pacemaker oscillator, we investigate which is the minimal action needed to control chimeras. We implement the pacemaker choosing one oscillator and making its links unidirectional. Our results show that a pacemaker induces chimeras for parameters and initial conditions for which they do not form spontaneously. Furthermore, the pacemaker attracts the incoherent part of the chimera state, thus controlling its position. Beyond that, we find that these control effects can be achieved with modifications of the network's connectivity that are less invasive than a pacemaker, namely, the minimal action of just modifying the strength of one connection allows one to control chimeras.

Chimera states in coupled logistic maps with additional weak nonlocal topology

We demonstrate the occurrence of coexisting domains of partially coherent and incoherent patterns or simply known as chimera states in a network of globally coupled logistic maps upon addition of weak nonlocal topology. We find that the chimera states survive even after we disconnect nonlocal connections of some of the nodes in the network. Also, we show that the chimera states exist when we introduce symmetric gaps in the nonlocal coupling between predetermined nodes. We ascertain our results, for the existence of chimera states, by carrying out the recurrence quantification analysis and by computing the strength of incoherence. We extend our analysis for the case of small-world networks of coupled logistic maps and found the emergence of chimeralike states under the influence of weak nonlocal topology.

Weak multiplexing in neural networks: Switching between chimera and solitary states

We investigate spatio-temporal patterns occurring in a two-layer multiplex network of oscillatory FitzHugh-Nagumo neurons, where each layer is represented by a nonlocally coupled ring. We show that weak multiplexing, i.e., when the coupling between the layers is smaller than that within the layers, can have a significant impact on the dynamics of the neural network. We develop control strategies based on weak multiplexing and demonstrate how the desired state in one layer can be achieved without manipulating its parameters, but only by adjusting the other layer. We find that for coupling range mismatch, weak multiplexing leads to the appearance of chimera states with different shapes of the mean velocity profile for parameter ranges where they do not exist in isolation. Moreover, we show that introducing a coupling strength mismatch between the layers can suppress chimera states with one incoherent domain (one-headed chimeras) and induce various other regimes such as in-phase synchronization or two-headed chimeras. Interestingly, small intra-layer coupling strength mismatch allows to achieve solitary states throughout the whole network.

Self-adaptation of chimera states

Chimera states in spatiotemporal dynamical systems have been investigated in physical, chemical, and biological systems, and have been shown to be robust against random perturbations. How do chimera states achieve their robustness? We uncover a self-adaptation behavior by which, upon a spatially localized perturbation, the coherent component of the chimera state spontaneously drifts to an optimal location as far away from the perturbation as possible, exposing only its incoherent component to the perturbation to minimize the disturbance. A systematic numerical analysis of the evolution of the spatiotemporal pattern of the chimera state towards the optimal stable state reveals an exponential relaxation process independent of the spatial location of the perturbation, implying that its effects can be modeled as restoring and damping forces in a mechanical system and enabling the articulation of a phenomenological model. Not only is the model able to reproduce the numerical results, it can also predict the trajectory of drifting. Our finding is striking as it reveals that, inherently, chimera states possess a kind of "intelligence" in achieving robustness through self-adaptation. The behavior can be exploited for the controlled generation of chimera states with their coherent component placed in any desired spatial region of the system.

Hierarchical transitions in multiplex adaptive networks of oscillatory units

In this work, we consider two-layer multiplex networks of coupled Stuart-Landau oscillators. The first layer contains oscillators with amplitude heterogeneity and all-to-all adaptive links, while the second layer contains identical oscillators all-to-all coupled by links with constant weights. The links between different layers are adaptive and organized in a one-to-one manner. We study the evolution of one-layer and two-layer networks depending on intra- and interlayer coupling strengths and show hierarchical transitions between oscillatory and quenched regimes.

Networks of coupled oscillators: From phase to amplitude chimeras

We show that amplitude-mediated phase chimeras and amplitude chimeras can occur in the same network of nonlocally coupled identical oscillators. These are two different partial synchronization patterns, where spatially coherent domains coexist with incoherent domains and coherence/incoherence referring to both amplitude and phase or only the amplitude of the oscillators, respectively. By changing the coupling strength, the two types of chimera patterns can be induced. We find numerically that the amplitude chimeras are not short-living transients but can have a long lifetime. Also, we observe variants of the amplitude chimeras with quasiperiodic temporal oscillations. We provide a qualitative explanation of the observed phenomena in the light of symmetry breaking bifurcation scenarios. We believe that this study will shed light on the connection between two disparate chimera states having different symmetry-breaking properties.

Existence and stability of chimera states in a minimal system of phase oscillators

We examine partial frequency locked weak chimera states in a network of six identical and indistinguishable phase oscillators with neighbour and next-neighbour coupling and two harmonic coupling of the form g ( ϕ ) = - sin ⁡ ( ϕ - α ) + r sin ⁡ 2 ϕ . We limit to a specific partial cluster subspace, reduce to a two dimensional system in terms of phase differences, and show that this has an integral of motion for α = π / 2 and r = 0 . By careful analysis of the phase space, we show that there is a continuum of neutrally stable weak chimera states in this case. We approximate the Poincaré return map for these weak chimera solutions and demonstrate several results about the stability and bifurcation of weak chimeras for small β = π / 2 - α and r that agree with numerical path-following of the solutions.

Robustness of chimera states in nonlocally coupled networks of nonidentical logistic maps

We investigate the dynamics of nonlocally coupled time-discrete maps with emphasis on the occurrence and robustness of chimera states. These peculiar, hybrid states are characterized by a coexistence of coherent and incoherent regions. We consider logistic maps coupled on a one-dimensional ring with finite coupling radius. Domains of chimera existence form different tongues in the parameter space of coupling range and coupling strength. For a sufficiently large coupling strength, each tongue refers to a wave number describing the structure of the spatial profile. We also analyze the period-adding scheme within these tongues and multiplicity of period solutions. Furthermore, we study the robustness of chimeras with respect to parameter inhomogeneities and find that these states persist for different widths of the parameter distribution. Finally, we explore the spatial structure of the chimera using a spatial correlation function.

Symmetries of Chimera States

Symmetry broken states arise naturally in oscillatory networks. In this Letter, we investigate chaotic attractors in an ensemble of four mean-coupled Stuart-Landau oscillators with two oscillators being synchronized. We report that these states with partially broken symmetry, so-called chimera states, have different setwise symmetries in the incoherent oscillators, and in particular, some are and some are not invariant under a permutation symmetry on average. This allows for a classification of different chimera states in small networks. We conclude our report with a discussion of related states in spatially extended systems, which seem to inherit the symmetry properties of their counterparts in small networks.

Three-dimensional chimera patterns in networks of spiking neuron oscillators

We study the stable spatiotemporal patterns that arise in a three-dimensional (3D) network of neuron oscillators, whose dynamics is described by the leaky integrate-and-fire (LIF) model. More specifically, we investigate the form of the chimera states induced by a 3D coupling matrix with nonlocal topology. The observed patterns are in many cases direct generalizations of the corresponding two-dimensional (2D) patterns, e.g., spheres, layers, and cylinder grids. We also find cylindrical and "cross-layered" chimeras that do not have an equivalent in 2D systems. Quantitative measures are calculated, such as the ratio of synchronized and unsynchronized neurons as a function of the coupling range, the mean phase velocities, and the distribution of neurons in mean phase velocities. Based on these measures, the chimeras are categorized in two families. The first family of patterns is observed for weaker coupling and exhibits higher mean phase velocities for the unsynchronized areas of the network. The opposite holds for the second family, where the unsynchronized areas have lower mean phase velocities. The various measures demonstrate discontinuities, indicating criticality as the parameters cross from the first family of patterns to the second.