Moderate Higher-Order Interactions Enhance Stability While Preserving Basin Structure

Basic Information

• Paper ID: 2510.13321
• Title: Moderate Higher-Order Interactions Enhance Stability While Preserving Basin Structure
• Authors: Zheng Wang, Jinjie Zhu, Xianbin Liu (State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics)
• Classification: nlin.AO (Nonlinear Science - Adaptation and Self-Organizing Systems)
• Publication Date: October 15, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2510.13321

Abstract

Synchronization is a ubiquitous phenomenon in complex systems. The Kuramoto model provides a canonical framework for understanding how coupled oscillators achieve collective rhythms. Traditional approaches focus on pairwise interactions, yet real systems frequently involve higher-order couplings among multiple elements. Previous studies demonstrated that higher-order interactions enrich dynamics but typically shrink the basin of attraction for synchronized states, making synchronization harder to achieve. This paper demonstrates that this picture is incomplete. Through systematic analysis of twisted states on ring networks, we identify a moderate coupling regime where higher-order interactions enhance stability without altering basin structure. The relative distribution among twisted states remains constant, while quasi-potential barriers deepen with increasing coupling strength. By measuring mean first passage times, we demonstrate that pairwise and higher-order couplings synergistically enhance stability, consistent with large deviation theory. These findings provide new insights into the role of higher-order interactions in synchronization.

Research Background and Motivation

Problem to be Addressed

This study aims to resolve a critical theoretical question: Can higher-order interactions enhance synchronization stability without altering the basin of attraction structure? This question challenges existing understanding that higher-order interactions typically shrink the basin of attraction for synchronized states.

Importance of the Problem

1. Theoretical Significance: Deepens understanding of higher-order interaction mechanisms in complex systems, particularly their dual effects on synchronization phenomena
2. Practical Applications: Provides theoretical guidance for designing robust engineering systems (e.g., power grids, neural networks)
3. Natural Phenomena Explanation: Helps explain why higher-order structures are so prevalent in nature

Limitations of Existing Approaches

1. Traditional Research primarily focuses on pairwise interactions, overlooking collective interactions widely present in real systems
2. Existing Higher-Order Studies mainly find that higher-order interactions shrink basins of attraction, lacking deep understanding of stability enhancement mechanisms
3. Lack of Systematic Analysis: Existing research fails to distinguish between basin geometry (basin size) and energy landscape (basin depth)

Research Motivation

The authors pose a fundamental question: Can we separate two complementary aspects of stability—basin structure (determining which initial conditions lead to synchronization) and basin depth (measuring resistance of synchronized states to perturbations)?

Core Contributions

1. Discovery of Moderate Higher-Order Coupling Regime: Identifies a parameter region where higher-order interactions enhance stability while maintaining basin structure invariance
2. Establishment of Quasi-Potential Barrier Framework: First application of Freidlin-Wentzell large deviation theory to higher-order interaction systems, quantifying basin depth through mean first passage time
3. Revelation of Synergistic Enhancement Mechanism: Demonstrates that pairwise and higher-order couplings synergistically enhance stability through energy mechanisms rather than geometric restructuring
4. Systematic Experimental Verification: Validates theoretical predictions through large-scale numerical experiments, including confirmation of exponential scaling relationships

Methodology Details

Task Definition

Study the stability of twisted states in Kuramoto oscillator systems with pairwise and ternary interactions on ring networks. Input: coupling parameters (σ, σ_Δ); Output: basin structure and quasi-potential barrier heights for different twisted states.

Model Architecture

Higher-Order Kuramoto Model

System dynamics described by:

θ̇ᵢ = ω + (σ/2r) Σⱼ sin(θⱼ - θᵢ) + (σ_Δ/[2r(2r-1)]) ΣⱼΣₖ sin(θⱼ + θₖ - 2θᵢ)

Where:

• θᵢ: phase of the i-th oscillator
• σ: pairwise coupling strength
• σ_Δ: ternary coupling strength
• r: coupling radius (r=2)
• n: number of oscillators (n=83)

Twisted State Definition

q-twisted state defined as:

θᵢ⁽ᵠ⁾ = 2πqi/n + Φ

Where q is the winding number and Φ is arbitrary global phase.

Local Coherence Measurement

Zᵢ = (1/(2r+1)) Σⱼ e^(iθⱼ)
R = (1/n) Σᵢ I(|Zᵢ| ≥ ρc)

Where ρc = 0.85 is the empirical threshold.

Technical Innovations

1. Quasi-Potential Barrier Analysis: First application of large deviation theory to high-dimensional non-gradient systems, measuring quasi-potential barriers through noise-induced mean first passage times
2. Dual Stability Framework: Distinguishes geometric stability (basin structure) from energetic stability (basin depth), providing more comprehensive stability analysis
3. Moderate Coupling Regime Identification: Systematically determines parameter regions where higher-order interactions enhance stability without altering basin structure

Experimental Setup

Data Generation

• Random Initial Conditions: 1.6×10⁵ uniformly distributed initial conditions in [-π, π)ⁿ
• Noise Experiments: 1600 sample trajectories with noise intensity D varying to ensure clear exponential scaling
• System Parameters: n=83 oscillators, coupling radius r=2

Evaluation Metrics

1. Basin Proportion: Fraction of initial conditions converging to different twisted states
2. Mean First Passage Time (MFPT): τₑ = Einf{t > 0 : θ(0) ∈ A, θ(t) ∈ B}
3. Quasi-Potential Barrier Height: Extracted through linear fitting of ln(τₑ) ~ Ψ/(2D²)

Comparison Methods

• Pure pairwise coupling systems (σ_Δ = 0)
• Systematic comparison across different coupling strength combinations
• Stability comparison among twisted states with different winding numbers

Implementation Details

• Numerical Integration: Sufficiently long integration time ensures convergence
• Noise Removal Protocol: Periodically remove noise and evolve for 100 time units to determine final state
• Linear Fitting: Use least squares to fit ln(τₑ) vs 1/D² relationship

Experimental Results

Main Results

Basin Structure Analysis

Figure 1 shows that under moderate ternary coupling, twisted states occupy the entire state space. When σ_Δ exceeds a critical threshold, non-twisted states rapidly emerge and eventually dominate. The key finding is that the critical threshold increases with pairwise coupling strength σ.

Stability of Twisted State Distribution

Figure 2 reveals the key finding: within the parameter region where twisted states dominate, the relative distribution of twisted states with different winding numbers remains remarkably stable:

• |q|=1 states consistently occupy the largest proportion (~45%)
• Synchronized states |q|=0 rank second (~25%)
• |q|=2 states account for ~15%
• High winding number states |q|≥4 are negligible

Quasi-Potential Barrier Depth Analysis

Core results from Figures 3-5:

1. Exponential Scaling Verification: ln(τₑ) vs 1/D² exhibits clear linear relationship, validating large deviation theory predictions
2. Synergistic Enhancement Effect: Slope k (∝ quasi-potential barrier height) systematically increases with σ and σ_Δ
3. Universality: All twisted states exhibit the same enhancement pattern

Quantitative Results

• Synchronized state escape: As σ_Δ increases from 0 to 2, slope k increases from 2.197 to 5.236 (Figure 3a)
• Pairwise coupling effect: As σ increases from 1 to 3, slope k increases from 1.708 to 6.636 (Figure 3b)
• Winding number dependence: High winding number states have shallower quasi-potential wells, consistent with physical intuition

Experimental Findings

1. Dual Effect Separation: Successfully separates effects of higher-order interactions on basin structure and basin depth
2. Synergistic Mechanism: Pairwise and higher-order couplings enhance stability through different mechanisms but synergistically
3. Universal Scaling Law: All measurements obey exponential scaling ln(τₑ) ~ 1/D²

Related Work

Main Research Directions

1. Classical Kuramoto Model: Kuramoto (1984), Strogatz (2000) established theoretical foundations for pairwise interactions
2. Network Synchronization: Arenas et al. (2008), Rodrigues et al. (2016) extended models to complex networks
3. Higher-Order Interactions: Skardal & Arenas (2019, 2020), Millán et al. (2020) pioneered higher-order Kuramoto model research

Relationship to Related Work

• Zhang et al. (2024) found higher-order interactions shrink but deepen basins; this paper further separates these two effects
• León et al. (2024) studied anomalous synchronization transitions; this paper focuses on stability enhancement mechanisms
• Muolo et al. (2025) found weak higher-order interactions enhance synchronization; this paper extends to moderate coupling regimes

Advantages of This Work

1. First systematic distinction between geometric and energetic stability
2. Provides rigorous quasi-potential barrier theoretical framework
3. Discovers new mechanisms of higher-order interactions

Conclusions and Discussion

Main Conclusions

1. Moderate Higher-Order Coupling Regime Exists: Within specific parameter ranges, higher-order interactions enhance stability without altering basin structure
2. Synergistic Enhancement Mechanism: Pairwise and higher-order couplings synergistically enhance stability by deepening quasi-potential barriers
3. Theoretical Validation: Large deviation theory successfully describes noise-induced transitions in high-dimensional non-gradient systems

Limitations

1. Parameter Range Restriction: Results hold only in moderate coupling regime; may not apply under strong coupling
2. Network Topology Limitation: Only considers ring networks; other topologies require further investigation
3. System Scale: Choice of n=83 based on computational feasibility; large-scale system behavior needs verification

Future Directions

1. Extend to other network topologies and system scales
2. Study effects of higher-order interactions (quaternary and beyond)
3. Validate theoretical predictions in practical systems (neural networks, power grids)

In-Depth Evaluation

Strengths

1. Theoretical Innovation: First application of quasi-potential barrier theory to higher-order interaction systems, providing new analytical framework
2. Experimental Sufficiency: Large-scale numerical experiments (1.6×10⁵ initial conditions, 1600 noise trajectories) ensure result reliability
3. Result Convincingness: Clear exponential scaling relationships and systematic parameter dependence support theoretical predictions
4. Writing Clarity: Clear logical structure, rigorous mathematical exposition, well-crafted figures

Weaknesses

1. Method Limitations: Quasi-potential barrier analysis requires small noise assumption, potentially limiting practical applications
2. Parameter Selection: Some parameters (e.g., ρc=0.85) based on empirical choice, lacking theoretical justification
3. Computational Complexity: MFPT measurement requires substantial computational resources, limiting system scale

Impact

1. Field Contribution: Provides new perspective for higher-order interaction research, potentially inspiring subsequent studies
2. Practical Value: Provides theoretical guidance for designing robust synchronous systems
3. Reproducibility: Detailed method description and explicit numerical parameters facilitate reproduction

Applicable Scenarios

1. Neural Networks: Understanding synergistic effects of multi-synaptic inputs in the brain
2. Power Systems: Designing more stable synchronous generator networks
3. Biological Systems: Explaining higher-order interaction effects in collective behavior

References

This paper cites 32 important references spanning multiple fields including synchronization theory, network science, and stochastic process theory, providing solid theoretical foundation for the research.

Overall Assessment: This is a high-quality theoretical physics paper that, through rigorous mathematical analysis and large-scale numerical experiments, reveals new mechanisms of higher-order interactions in synchronization. The research methodology is innovative, and results have important theoretical significance and potential practical value, making significant contributions to understanding higher-order interactions in complex systems.