Experimental observation of chimera and cluster states in a minimal globally coupled network

I. INTRODUCTION

• Chimeras were initially thought to require large, nonlocally coupled networks and special initial conditions, but later studies found them in small and globally coupled networks and from random starts.They have also been observed experimentally in several physical and electronic systems.Source: §I · p. 1
• Chimeras often appear in multistable regions alongside other synchronous patterns, and recent work proposed a four-laser network linking chimera emergence to multistability.Cluster synchrony and chimera research had largely developed independently, despite recent links between them in globally coupled networks.Source: §I · p. 1
• The paper experimentally observes chimera and other partially synchronous states in four globally coupled opto-electronic oscillators and applies symmetry-based methods to calculate their linear stability.It also discusses the importance of multistability among partially synchronous states.Source: §I · p. 2

A. Experimental apparatus

• Each node uses a laser diode, Mach–Zehnder modulator, photoreceiver, and DSP board to filter, delay, amplify, and feed back its electrical signal.Coupling signals are optically combined and processed as a second electrical signal, with coupling delay τc generally distinct from feedback delay τf.Source: A. Experimental apparatus · p. 2
• The experiment fixes round-trip gain at β = 3.8 and feedback delay at τf = 1.4 ms, choosing phase bias so an uncoupled node behaves chaotically.The filter uses a two-pole Butterworth bandpass with cutoff frequencies ωH/2π = 100 Hz and ωL/2π = 2.5 kHz.Source: A. Experimental apparatus · p. 2
• For identical global coupling, each of four nodes receives input from the other three; the study varies coupling strength ε and coupling delay τc.Coupling is not generally Laplacian and becomes Laplacian only when τc = τf; trials initialize from noise, allow uncoupled transients to decay, then classify the final 400 ms.Source: A. Experimental apparatus · p. 2

B. Observation of partial synchronization

• The five four-node outcomes are global synchrony, doublet-doublet, triplet-singlet, doublet-singlet-singlet (DSS), and desynchrony.The authors call doublet-doublet and triplet-singlet cluster states, and DSS a chimera state.Source: B. Observation of partial… · p. 2
• The experiment observes all possible synchronous states, including a chimera that persists for many delay times and appears stable; the four-node network is the minimal globally coupled network that can support one.Different initial conditions place nodes in different clusters, indicating the patterns are not due to parameter mismatch.Source: B. Observation of partial… · pp. 2–3
• The network sketch distinguishes each node’s self-feedback delay τf from its coupling delay τc to every other node.The experimental schematic distinguishes optical connections from electronic ones.Source: Figure 1 · p. 3

III. SYNCHRONIZATION PATTERNS AND NETWORK SYMMETRIES

• Figure 2 presents the four-node synchronous patterns, with nodes of the same color grouped into the same cluster.The illustrated patterns are global, doublet-doublet, triplet-singlet, and doublet-singlet-singlet synchrony.Source: Figure 2 · p. 3, §III · p. 3
• The reported examples use ε = 0.40 and τc = 1.8 ms for global synchrony, ε = 0.45 and τc = 1.8 ms for both cluster states, and ε = 0.40 and τc = 2.3 ms for the chimera.Numerical simulations of the governing equations produce similar time series.Source: §III · p. 3
• Subgroups of the global permutation symmetry allow any partition of N oscillators, including one synchronized cluster and N − Ns desynchronized singlets; stability analysis determines whether such a chimera can be observed.Subgroup orbits provide the partial-symmetry clusters that can emerge through symmetry breaking.Source: §III · p. 4

• A symmetry-preserving coordinate transformation block-diagonalizes the adjacency matrix and decouples manifold motion from transverse directions, reducing the stability calculation.Source: A. Stability analysis · p. 4
• For the four-node DSS chimera, only the transverse IRR basis vector ∆v⊥ = ∆v4 is needed after transformation.Source: A. Stability analysis · p. 4
• At regularly spaced ε–τc points, at least 20 random-start experimental trials were run; shaded regions mark negative-LLE stability, while dots and crosses mark observed and unobserved states.Source: A. Stability analysis · p. 5, Figure 3 · p. 5
• Experiments and calculated stability agree quite well for all four synchronous states; slight mismatches are attributed to finite random-initial-condition sampling.Source: A. Stability analysis · p. 5
• For a 10-node global network, stability analysis of a chimera with one cluster of five and five singlets agrees with direct simulations.Source: A. Stability analysis · p. 5, Figure 4 · p. 6

IV. CHIMERA STATES AND MULTISTABILITY

• The authors consider a globally synchronized pattern stable if any of its dynamical realizations is stable, without distinguishing chaotic from quasiperiodic behavior.The globally synchronized state can appear chaotic or quasiperiodic.Source: §IV · p. 6
• Stable calculated DSS chimeras coincide well with multistable regions, and experiments also find chimeras at parameter values exhibiting multistability.Multistability means at least two of global synchrony, doublet-doublet, and triplet-singlet states are stable.Source: Figure 5 · p. 6, §IV · p. 6
• In this system, mismatch between coupling and feedback delays breaks maximal symmetry; when the delays match, only global synchrony is observed.The matching-delay result is also supported by prior work on these opto-electronic oscillator networks [27].Source: §IV · p. 6
• Different initial conditions and basins of attraction yield fragmented stability regions, including global synchrony and complete desynchrony.The parameter space can be multistable both within one synchronous pattern and between different patterns.Source: §IV · p. 7
• Multistability is presented as a requirement for chimeras, though its physical mechanism can differ across systems.Time-delay coupling induces multistability in this system, whereas amplitude-phase coupling does so in the cited laser system [21].Source: §IV · p. 7

V. DISCUSSION

• Four globally coupled oscillators support experimentally observed chimera states from random initial conditions, and stability calculations show they are stable physical states, not long transients.The study also reports all possible partially synchronous states and agreement between calculated linear stability and experiments.Source: §V · p. 7
• Isolated desynchronization lets some clusters separate from synchrony while other clusters remain synchronized through the network’s subgroup symmetries.In the DSS chimera, the two doublet oscillators receive the same total signal from the desynchronized singlets.Source: §V · p. 7
• The group-theoretical analysis and isolated-desynchronization mechanism extend to networks with identically synchronized coherent populations, including nonlocal and star networks.The analysis is not restricted to globally coupled networks.Source: §V · p. 7
• The amount of experimental heterogeneity compatible with the group-theoretical stability analysis remains an open question.Persistent chimera and cluster states were observed despite small experimental heterogeneities, although simulations and stability analysis assume identical oscillators and coupling.Source: §V · p. 7
• The authors identify two different time delays as the mechanism breaking the symmetry present in Laplacian coupling and generating multistability in their system.They note that the mechanism generating multistability can differ between systems.Source: §V · p. 7

• Computing T generally requires finding IRRs for the synchronous state's symmetry subgroup, but for the DSS chimera it can be determined by inspection.Source: APPENDIX: CALCULATION OF THE IRR… · pp. 7–8
• The DSS chimera has a three-dimensional synchronization manifold and, with four nodes, only one transverse direction.The manifold comprises one dimension for the synchronized doublet and one for each singlet cluster.Source: APPENDIX: CALCULATION OF THE IRR… · p. 8
• Stacking normalized row basis vectors with manifold directions first makes T block-diagonalize the adjacency matrix into manifold and transverse perturbations.The upper block corresponds to synchronization-manifold perturbations and the lower block to transverse perturbations.Source: APPENDIX: CALCULATION OF THE IRR… · p. 8