Critical Behavior at the Onset of Multichimera States in a Coupled-Oscillator Array
Pith reviewed 2026-05-24 20:35 UTC · model grok-4.3
The pith
The transition from global synchronization to multichimera states in a chain of coupled oscillators is continuous and belongs to a new universality class of non-equilibrium critical phenomena.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper establishes that the onset of multichimera states occurs via a continuous transition where the asynchronous fraction ρ_a scales as (α − α_c)^β_a and the spatio-temporal gaps between asynchronous sites exhibit power-law distributions at the critical point. The resulting exponents are shown to lie outside the universality class of (1+1)-dimensional directed percolation and other known absorbing-state transitions, indicating a distinct class of non-equilibrium critical phenomena. Traveling waves rejuvenate asynchronous clusters by mediating non-local interactions between them.
What carries the argument
Traveling waves that rejuvenate asynchronous clusters by mediating non-local interactions between them.
Load-bearing premise
The numerically identified critical point and extracted exponents reflect the infinite-system limit rather than being controlled by finite-size effects or the specific frequency distribution and coupling range chosen for the simulations.
What would settle it
A simulation on arrays orders of magnitude larger that yields exponents converging to those of directed percolation would falsify the claim of a new universality class.
Figures
read the original abstract
We numerically investigate the onset of multi-chimera states in a linear array of coupled oscillators. As the phase delay $\alpha$ is increased, they exhibit a continuous transition from the globally synchronized state to the multichimera state consisting of asynchronous and synchronous domains. Large-scale simulations show that the fraction of asynchronous sites $\rho_a$ obeys the power law $\rho_a \sim (\alpha - \alpha_c)^{\beta_a}$, and that the spatio-temporal gaps between asynchronous sites show power-law distributions at the critical point. The critical exponents are distinct from those of the (1+1)-dimensional directed percolation and other absorbing-state phase transitions, indicating that this transition belongs to a new class of non-equilibrium critical phenomena. Crucial roles are played by traveling waves that rejuvenate asynchronous clusters by mediating non-local interactions between them.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript numerically studies the onset of multichimera states in a one-dimensional array of nonlocally coupled phase oscillators as the phase delay parameter α is increased. It reports a continuous transition from global synchronization to a state with coexisting synchronous and asynchronous domains, with the asynchronous-site density ρ_a obeying a power law ρ_a ∼ (α − α_c)^β_a and the spatio-temporal gaps between asynchronous clusters exhibiting power-law distributions at criticality. The extracted exponents are stated to differ from those of (1+1)-dimensional directed percolation, leading to the claim of a new non-equilibrium universality class, with traveling waves mediating nonlocal interactions between clusters.
Significance. If the numerical identification of α_c and the exponents can be shown to survive finite-size scaling and extrapolation to the thermodynamic limit, the result would establish a novel absorbing-state transition driven by nonlocal coupling and traveling waves. This would be of interest to the study of chimera states and non-equilibrium phase transitions, as it suggests that the combination of oscillator dynamics and nonlocal interactions can produce critical behavior outside known classes such as directed percolation.
major comments (3)
- [§3] §3 (Numerical results): The location of the critical point α_c and the subsequent power-law fits for ρ_a and the gap distributions are performed on finite arrays without any reported finite-size scaling analysis, data-collapse procedure, or extrapolation in system size L. Because the central claim is that the exponents differ from those of (1+1)D directed percolation and therefore define a new class, the absence of these controls leaves open the possibility that the observed power laws are slow crossovers or boundary-induced artifacts rather than true critical behavior.
- [§3] §3, paragraph on gap statistics: The power-law exponents for the spatio-temporal gap distributions are asserted to be distinct from directed percolation, yet no table or figure quantifies the system sizes (L, T), the number of independent realizations, or the fitting range used to extract these exponents. Without convergence checks as L increases, the distinction from DP cannot be established at the level required for a new universality class.
- [§2] §2 (Model definition): The frequency distribution and the range of the nonlocal coupling are fixed parameters whose variation could affect the apparent location of α_c and the measured exponents; no robustness tests against changes in these parameters or against different boundary conditions are presented, even though traveling waves are stated to play a crucial role.
minor comments (2)
- [Abstract] The abstract states that “large-scale simulations” were performed but supplies no numerical values for the largest system sizes; these values should be stated explicitly in the main text or a methods subsection.
- [§2] Notation for the order parameter ρ_a is introduced without an explicit equation; adding a numbered equation in §2 would improve clarity when the power-law scaling is later discussed.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The concerns about finite-size effects and robustness are valid and will be addressed in the revision. Below we respond point by point to the major comments.
read point-by-point responses
-
Referee: [§3] §3 (Numerical results): The location of the critical point α_c and the subsequent power-law fits for ρ_a and the gap distributions are performed on finite arrays without any reported finite-size scaling analysis, data-collapse procedure, or extrapolation in system size L. Because the central claim is that the exponents differ from those of (1+1)D directed percolation and therefore define a new class, the absence of these controls leaves open the possibility that the observed power laws are slow crossovers or boundary-induced artifacts rather than true critical behavior.
Authors: We agree that systematic finite-size scaling (FSS) analysis, data collapse, and extrapolation are necessary to firmly establish the critical point and rule out crossovers. In the revised manuscript we will add simulations for multiple larger system sizes (up to L = 20000), perform data collapse of ρ_a, and extrapolate α_c(L) to the thermodynamic limit. These additions will strengthen the evidence that the observed scaling defines a new universality class. revision: yes
-
Referee: [§3] §3, paragraph on gap statistics: The power-law exponents for the spatio-temporal gap distributions are asserted to be distinct from directed percolation, yet no table or figure quantifies the system sizes (L, T), the number of independent realizations, or the fitting range used to extract these exponents. Without convergence checks as L increases, the distinction from DP cannot be established at the level required for a new universality class.
Authors: We will include a new table (and accompanying text) in the revised §3 that reports the exact values of L and T, the number of independent realizations (typically 100–500), the fitting ranges, and the goodness-of-fit measures for the gap distributions. We will also add a supplementary figure demonstrating convergence of the extracted exponents with increasing L to support the distinction from directed percolation. revision: yes
-
Referee: [§2] §2 (Model definition): The frequency distribution and the range of the nonlocal coupling are fixed parameters whose variation could affect the apparent location of α_c and the measured exponents; no robustness tests against changes in these parameters or against different boundary conditions are presented, even though traveling waves are stated to play a crucial role.
Authors: The frequency distribution (uniform on [−ω,ω]) and coupling range (exponential kernel with fixed decay length) are the standard choices for this nonlocal phase-oscillator model; the transition is driven by the phase delay α. Nevertheless, to address the concern we will add a short robustness subsection (or supplementary material) showing that the critical exponents remain consistent under modest variations of the frequency width and kernel range, as well as for periodic boundary conditions. revision: yes
Circularity Check
Numerical simulation with external benchmark comparison exhibits no circularity
full rationale
The paper reports results from direct large-scale numerical simulations of a coupled-oscillator array, locating the critical value α_c and extracting the power-law exponent β_a for the asynchronous fraction ρ_a along with gap-distribution exponents. These quantities are then compared to the independently known exponents of (1+1)D directed percolation. No equations or procedures reduce a claimed prediction to a fitted input by construction, no self-citations are load-bearing, and no ansatz or uniqueness theorem is smuggled in; the chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Critical Behavior at the Onset of Multichimera States in a Coupled-Oscillator Array
We numer- ically solved (1) with a periodic boundary condition for the system size N up to 2 24 = 16777216. The coupling rangeL = 5 is used unless otherwise stated. By choosing L≪N, we can study the statistical behavior of a large number of asynchronous clusters. A uniformly random and spatially uncorrelated distributions of φx in [0, 2π] is used for the ...
work page internal anchor Pith review Pith/arXiv arXiv 1907
-
[2]
D. M. Abrams and S. H. Strogatz, Phys. Rev. Lett. 93, 174102 (2004)
work page 2004
-
[3]
Y. Kuramoto and D. Battogtokh, Nonlin. Phenom. Com- plex Syst. 5, 380 (2002)
work page 2002
- [4]
- [5]
-
[6]
D. M. Abrams and S. H. Strogatz, Int. Journal of Bifur- cat. Chaos 16, 21 (2006)
work page 2006
-
[7]
D. M. Abrams, R. Mirollo, S. H. Strogatz, and D. A. Wiley, Phys. Rev. Lett. 101, 084103 (2008)
work page 2008
-
[8]
O. E. Omelchenko, Y. L. Maistrenko, and P. A. Tass, Phys. Rev. Lett. 100, 044105 (2008)
work page 2008
-
[9]
G. Bordyugov, A. Pikovsky, and M. Rosenblum, Phys. Rev. E 82, 035205(R) (2010)
work page 2010
-
[10]
O. E. Omel’chenko, Nonlinearity 26, 2469 (2013)
work page 2013
-
[11]
G. C. Sethia and A. Sen, Phys. Rev. Lett. 112, 144101 (2014)
work page 2014
-
[12]
A. Yeldesbay, A. Pikovsky, and M. Rosenblum, Phys. Rev. Lett. 112, 144103 (2014)
work page 2014
-
[13]
M. J. Panaggio and D. M. Abrams, Nonlinearity , R67 (2015)
work page 2015
-
[14]
C. R. Laing, Physica D 238, 1569 (2009)
work page 2009
-
[15]
C. R. Laing, Chaos 19, 013113 (2009)
work page 2009
- [16]
-
[17]
I. Omelchenko, A. Provata, J. Hizanidis, E. Sch¨ oll, and P. H¨ ovel, Phys. Rev. E91, 022917 (2015)
work page 2015
- [18]
-
[19]
J. Hizanidis, N. E. Kouvaris, G. Zamora-L´ opez, A. D´ ıaz- Guilera, and C. G. Antonopoulos, Sci. Rep. 6, 19845 (2016)
work page 2016
- [20]
-
[21]
M. R. Tinsley, S. Nkomo, and K. Showalter, Nat. Phys. 8, 662 (2012)
work page 2012
- [22]
-
[23]
J. F. Totz, J. Rode, M. R. Tinsley, K. Showalter, and H. Engel, Nat. Phys. 14, 282 (2018)
work page 2018
-
[24]
E. A. Martens, S. Thutupalli, A. Fourri` ere, and O. Hal- latschek, Proc. Natl. Acad. Sci. 110, 10563 (2013)
work page 2013
-
[25]
T. Kapitaniak, P. Kuzma, J. Wojewoda, K. Czolczynski, and Y. Maistrenko, Sci. Rep. 4, 6379 (2014)
work page 2014
-
[26]
D. P. Rosin, D. Rontani, N. D. Haynes, E. Sch¨ oll, and D. J. Gauthier, Phys. Rev. E 90, 030902(R) (2014)
work page 2014
-
[27]
L. V. Gambuzza, A. Buscarino, S. Chessari, L. Fortuna, R. Meucci, and M. Frasca, Phys. Rev. E 90, 032905 (2014)
work page 2014
- [28]
- [29]
-
[30]
G. C. Sethia, A. Sen, and F. M. Atay, Phys. Rev. Lett. 100, 144102 (2008)
work page 2008
-
[31]
G. C. Sethia, A. Sen, and G. L. Johnston, Phys. Rev. E 88, 042917 (2013)
work page 2013
-
[32]
I. Omelchenko, O. E. Omel’chenko, P. H¨ ovel, and E. Sch¨ oll, Phys. Rev. Lett.110, 224101 (2013)
work page 2013
- [33]
-
[34]
J. Xie, E. Knobloch, and H.-C. Kao, Phys. Rev. E 90, 022919 (2014)
work page 2014
-
[35]
I. Omelchenko, A. Zakharova, P. H¨ ovel, J. Siebert, and E. Sch¨ oll, Chaos25, 083104 (2015)
work page 2015
-
[36]
J. Xie, E. Knobloch, and H.-C. Kao, Phys. Rev. E 92, 042921 (2015)
work page 2015
- [37]
-
[38]
O. E. Omel’chenko, M. Wolfrum, and C. R. Laing, Chaos 24, 023102 (2014)
work page 2014
-
[39]
M. Wolfrum, S. V. Gurevich, and O. E. Omelchenko, Nonlinearity 29, 257 (2016)
work page 2016
-
[40]
O. E. Omelchenko, Nonlinearity 31, R121 (2018)
work page 2018
- [41]
-
[42]
K. A. Takeuchi, M. Kuroda, H. Chat´ e, and M. Sano, Phys. Rev. Lett. 99, 234503 (2007)
work page 2007
-
[43]
G. Lemoult, L. Shi, K. Avila, S. V. Jalikop, M. Avila, and B. Hof, Nat. Phys. 12, 254 (2016)
work page 2016
- [44]
- [45]
-
[46]
R. M. Ziff, E. Gulari, and Y. Barshad, Phys. Rev. Lett. 56, 2553 (1986)
work page 1986
-
[47]
D. Mollison, J. R. Stat. Soc. Series B Stat. Methodol. 39, 283 (1977). 4 FIG. 1. Spatio-temporal patterns of the phase difference for (a) α = 0.43, (b) α = 0.44, and (c) α = 0.45 in an N = 512 system. They are dominated by synchronized, traveling waves and multichimera states in the late stage, respectively. The spatial profiles of|∆φx| in (d) and (e) are t...
work page 1977
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.