Three-Body Interactions Can Split Identical Oscillators Into Chimeras

When Three Bodies Are All You Need

Imagine a crowd where everyone speaks the same language, moves at the same natural pace, and follows the same simple rules—yet somehow, half the crowd falls into lockstep synchronization while the other half churns in chaotic disorder. No leader. No difference in the individuals. No pairwise influence whatsoever. Just three-body interactions, rippling through the system like whispers in a room where nobody is directly talking to anyone else.

This is the chimera state: a striking phenomenon in which symmetry breaks spontaneously in a perfectly homogeneous system. First described theoretically in 2002 by Yoshiki Kuramoto and later popularized by Daniel Abrams and Steven Strogatz, chimera states have been observed in lasers, chemical reactions, power grids, and—just recently—in the synchronized flashing of firefly swarms. They seem to violate our intuition about how identical elements should behave. If everyone is the same, shouldn't they all do the same thing?

For two decades, the answer researchers gave was: not necessarily—but you still need pairwise interactions. Some oscillators must directly influence others through simple one-to-one connections. The coupling might be nonlocal, spanning across the system, but it was always fundamentally pairwise. Two oscillators, talking to each other.

A new paper published on arXiv in July 2026 flips this assumption on its head. Sudo Yi, Gugyoung Kim, Mi Jin Lee, S.-W. Son, and B. Kahng, working across three Korean institutions, have demonstrated that chimera states can arise from purely triadic interactions—where three oscillators influence each other simultaneously, and zero pairwise coupling exists. The work establishes triadic coupling as what the authors call a "minimal mechanism" for chimera formation, opening a new chapter in our understanding of how complex collective behavior emerges from simple many-body rules.

The Science: Building a World Without Pairwise Coupling

The researchers set out to answer a deceptively simple question: Can higher-order interactions alone, without any pairwise coupling, generate chimera states? By "higher-order," they mean interactions that involve three or more elements at once—triadic interactions, in this case. In many physical systems, such interactions are present alongside pairwise ones, making it impossible to disentangle their intrinsic effects. The researchers wanted to isolate that effect completely.

Their model consists of 2N oscillators divided into two distinct groups. Each oscillator's phase evolves according to an equation that captures triadic coupling and nothing else. The key term in the dynamics involves a sine function with three phases: two oscillators from potentially different groups influence a third, with a phase lag parameter α that modulates how the synchronization unfolds. The coupling strength depends on whether the oscillators belong to the same group or different groups, quantified by a parameter μ that measures inter-group interaction strength.

This setup explicitly excludes pairwise coupling. No term in the equation describes direct two-oscillator interaction. The only influence any oscillator feels comes through triadic interactions—three-body encounters that reshape its phase.

The coupling tensor K determines how strongly triadic interactions proceed. When all three oscillators belong to the same group, the interaction strength is set to unity. When one oscillator comes from the opposite group, the strength drops to μ. When two oscillators come from the opposite group, it drops further to μ². This creates a hierarchy: intra-group triadic interactions dominate, while inter-group interactions are suppressed by the parameter μ ∈ [0, 1).

The researchers then applied a modified version of the Ott–Antonsen (OA) ansatz—a powerful analytical technique that reduces infinite-dimensional dynamics to a handful of macroscopic variables. The standard OA ansatz works well when oscillator phases cluster into a single unimodal distribution. But triadic interactions, the researchers discovered, naturally generate bimodal distributions: two clusters of synchronized oscillators separated by a phase difference of π (180 degrees). This π-symmetry, inherent to the triadic coupling structure, cannot be captured by the standard approach.

To solve this, the researchers decomposed each group's phase density into two anti-phase subpopulations, introducing an asymmetry parameter η that captures the relative weight of each subpopulation. Applying the OA ansatz separately to each subpopulation, they derived a reduced description parameterized by the complex amplitude a_σ = r_σ e^{-iψ_σ}, where r_σ measures the synchronization within each peak of the bimodal distribution and ψ_σ captures the collective phase. This "bimodal Ott–Antonsen reduction" yields an exact low-dimensional description of the full dynamics.

Numerical simulations were conducted with N₁ = N₂ = 10,000 oscillators in each group, placing the system firmly in the thermodynamic limit where finite-size effects become negligible. The initial conditions played a crucial role: Group 1 was prepared in a symmetry-broken synchronized state (η₁ ≠ 0.5, Q₁ = 1), while Group 2 was initialized in a symmetric bimodal state (η₂ = 0.5, 0 < Q₂ < 1). This asymmetry between the groups acts as a dynamical trigger, driving the emergence of chimera behavior.

What They Found: The Anatomy of a Triadic Chimera

The reduced dynamics reduces to two governing equations for the amplitude r (essentially the synchronization level in Group 2) and the phase difference φ between the two groups:

$$\dot{r}=\frac{1}{2r}(1-r^4){\mu }^2(2\eta -1{)}^2\cos (2\varphi +\alpha )$$$$\dot{\varphi }=-\frac{1}{2r^2}(1+r^4){\mu }^2(2\eta -1{)}^2\sin (2\varphi +\alpha )+(2\eta -1{)}^2\sin \alpha$$

These equations admit steady-state solutions that reveal the critical threshold separating synchronized behavior from chimera states. The boundary occurs at:

$$\mu =\sqrt{\sin \alpha }$$

Above this line (μ ≥ √sin α), the system converges to full synchronization: both groups achieve Q₁ = Q₂ = 1, meaning all oscillators lock into coherent bimodal patterns. Below this line (μ < √sin α), chimera states emerge—one group remains fully synchronized (Q₁ = 1) while the other exhibits partial synchronization (Q₂ < 1, equivalent to r < 1).

Chimera boundary in (α, μ) space at μ = 0.5

Label	Value
α = 0.2	0.95 Q₂ value
α = 0.3	0.6 Q₂ value
α = 0.5	0.45 Q₂ value
α = 0.7	0.35 Q₂ value
α = 0.9	0.3 Q₂ value

μ = 0.5 fixed. Below the theoretical boundary μ = √sin α, the system exhibits breathing chimera behavior with persistent oscillations in Q₂.

The chimera states that emerge are not static. They are breathing chimeras, characterized by persistent periodic oscillations in the partially synchronized group. The oscillation amplitude and period depend on both the system parameters and the initial conditions.

For a representative case with α = π/4 and μ = 0.5, the researchers found:

• Oscillation amplitude: ΔQ₂ ≈ 0.266 (theory) vs. 0.266 (simulation)
• Oscillation period: T ≈ 475.0 (theory) vs. T ≈ 474.9 (simulation)

The agreement between theory and simulation is striking—differences on the order of 10⁻⁴ for amplitude and 10⁻² for period. This precision validates the bimodal OA reduction as an exact description of the macroscopic dynamics.

Breathing dynamics of Q₂ over time

Label	Value
t = 0	0.6 Q₂
t = 100	0.75 Q₂
t = 200	0.85 Q₂
t = 300	0.8 Q₂
t = 400	0.6 Q₂
t = 500	0.75 Q₂
t = 600	0.85 Q₂
t = 700	0.8 Q₂

Breathing chimera with persistent oscillations: amplitude ΔQ₂ ≈ 0.266, period T ≈ 475 (theory). Simulation matches to within 10⁻⁴ for amplitude.

The reduced dynamics can be rewritten in an even more compact form by combining the amplitude and phase into a single complex order parameter Z = r² e^{iΦ}. The resulting equation is a Riccati-type equation:

$$\dot{Z}=A(1-Z^2)+2iBZ$$

where A = μ²(1 - 2η)² and B = sin α (1 - 2η)². This equation admits a conserved invariant that confines the dynamics to a closed orbit in the complex plane. The breathing chimera, in geometric terms, is a circular trajectory in this abstract space—a periodic orbit that neither grows nor decays, but cycles endlessly.

Phase portraits: fixed point vs breathing chimera orbits

Label	Value
α = 0.2	1 Orbit type
α = 0.3	0 Orbit type
α = 0.5	0 Orbit type
α = 0.7	0 Orbit type
α = 0.9	0 Orbit type

μ = 0.5. The closed orbits in the complex plane represent the geometric interpretation of breathing chimera as persistent periodic motion.

The period of this orbit is given by:

$$T/\pi =(1-2\eta {)}^{-2}({\sin }^2\alpha -{\mu }^4{)}^{-1/2}$$

This expression reveals that the period diverges as μ → √sin α, indicating critical slowing down near the bifurcation boundary. The oscillation amplitude, by contrast, depends on initial conditions through the conserved invariant:

$$\Delta Q\sim \frac{{\mu }^2}{\sin \alpha }$$

As the inter-group coupling μ increases, the amplitude grows—more cross-group triadic interactions disrupt the synchronization in Group 2, pushing it further from lockstep. As the phase lag α increases (within the chimera region), the amplitude shrinks, reflecting the stabilizing role of larger phase delays.

Why This Changes Things: The Minimal Mechanism

The significance of this work extends far beyond the elegant mathematics of Riccati equations and complex-plane orbits. At its core, it answers a fundamental question about the minimal requirements for symmetry breaking in collective systems.

For twenty years, researchers assumed that chimera states required pairwise interactions as a necessary ingredient. The logic was intuitive: for some oscillators to synchronize while others don't, there must be some mechanism allowing differential influence—some way for local interactions to create patches of coherence within a sea of disorder. Pairwise coupling, even nonlocal pairwise coupling, seemed essential.

Yi and colleagues have dismantled this assumption. Their model contains zero pairwise coupling. The only interactions are triadic—three oscillators influencing each other simultaneously. Yet chimera states emerge robustly, with one group falling into full synchronization while the other oscillates perpetually between partial synchronization and partial disorder.

This matters for several reasons.

First, it establishes triadic coupling as a sufficient mechanism for chimera formation. Sufficient, not merely helpful. The researchers have proven that you don't need pairwise interactions to get this form of spontaneous symmetry breaking—you only need three-body encounters.

Second, it identifies a new pathway to chimera states that bypasses mechanisms previously considered essential. Classical chimera theory emphasizes the role of phase lag and nonlocal coupling topology. The triadic mechanism operates differently: the π-symmetry inherent to three-body interactions naturally generates bimodal phase distributions, and an asymmetry in the initial conditions (one group more balanced than the other) triggers the coexistence of synchronized and unsynchronized populations.

Third, the work provides an analytically tractable framework for studying symmetry-broken collective dynamics in systems dominated by many-body interactions. The bimodal OA reduction yields exact low-dimensional equations—a rare achievement in nonlinear dynamics. Researchers can now explore breathing chimera behavior analytically, deriving conditions for emergence, stability, and bifurcation with mathematical rigor rather than relying solely on simulation.

The geometric interpretation deserves particular attention. By recasting the reduced dynamics as a Riccati equation, the researchers have given breathing chimeras a vivid spatial representation: closed periodic orbits in the complex plane. This isn't merely aesthetic—it offers intuition. The chimera state is not a fixed point but a trajectory, an orbit that the system traverses endlessly. The invariant of motion that confines the trajectory to a closed curve ensures the breathing continues indefinitely, neither dying out nor exploding into full synchronization.

For the broader study of collective phenomena, this work suggests that higher-order interactions may support symmetry-broken states that pairwise coupling cannot easily produce. In neural systems, where collective activity often depends on interactions among many neurons rather than pairwise synapses, such mechanisms could underlie interesting macroscopic patterns. In ecological systems governed by group-level dynamics—predator-prey interactions involving multiple species, for instance—triadic interactions might generate collective states with no analog in pairwise models.

What's Next: Extensions, Caveats, and Open Questions

The paper is careful to acknowledge its limitations. The analysis focuses on identical oscillators with all-to-all purely triadic interactions and prescribed asymmetric initial conditions. These assumptions make the model analytically tractable but leave open questions about robustness.

Frequency heterogeneity, for instance, remains unexplored. Real-world oscillators rarely share exactly the same natural frequency. Introducing disorder in the ω parameter could destabilize the chimera states or modify the breathing dynamics. The researchers note that addressing this extension would clarify the relevance of triadic chimera mechanisms to more realistic systems.

Noise is another frontier. Biological and neural systems are inherently noisy; stochastic fluctuations might disrupt the delicate periodic orbits that characterize breathing chimeras. Whether robust breathing behavior survives the presence of noise remains an open question.

The topology of interactions also warrants investigation. The current model assumes all-to-all coupling—all oscillators participate in all triadic interactions. Real systems are often sparse, with triadic interactions occurring only among connected nodes in a hypergraph. How sparsity affects the emergence and stability of triadic chimeras is not yet known.

Perhaps most interesting is the question of mixed coupling. Real systems typically exhibit both pairwise and higher-order interactions simultaneously. The researchers' finding that triadic coupling alone suffices for chimera formation raises the question of what happens when both mechanisms are present. Do they compete? Reinforce? Generate novel hybrid states?

The breathing chimera's sensitivity to initial conditions, while analytically elegant, also hints at practical challenges. The asymmetric initialization—one group in a symmetry-broken state, the other in a symmetric bimodal state—may be difficult to achieve in experimental realizations. Whether nature routinely provides such asymmetry, or whether additional mechanisms are needed to initialize triadic chimeras, remains to be seen.

Despite these open questions, the work represents a significant conceptual advance. It establishes a new minimal mechanism for chimera formation, provides a tractable mathematical framework for analyzing it, and opens a pathway for understanding symmetry-broken collective dynamics in systems where many-body interactions are intrinsic rather than incidental.

The fireflies that synchronize their flashing, the neurons that coordinate their firing, the neurons in a sleeping bird that allow one hemisphere to rest while the other stays alert—these phenomena may share a common thread that runs through triadic interactions. Three-body encounters, rippling through a system, may be enough to split a homogeneous crowd into synchronized and unsynchronized factions. The math, at least, says so.

As the researchers acknowledge, the road from theoretical model to biological relevance is long. But the destination has become clearer: somewhere in the space of many-body interactions, symmetry breaks in ways that pairwise thinking alone cannot explain.