Higher-order interactions promote chimera states

Bifurcations and collective states of Kuramoto oscillators with higher-order interactions and rotational symmetry breaking

We study a network of identical Kuramoto oscillators with higher-order interactions that also break the rotational symmetry of the system. To gain analytical insights into this model, we use the Watanabe-Strogatz Ansatz, which allows us to reduce the dimensionality of the original system of equations. The study of stability and bifurcations of the reduced system reveals a codimension two Bogdanov-Takens bifurcation and several other associated bifurcations. Such analysis is corroborated by numerical simulations of the associated Kuramoto system, which, in turn, unveils a variety of collective behaviors such as synchronized motion, oscillation death, chimeras, incoherent states, and traveling waves. Importantly, this system displays a case where alternating chimeras emerge in an indistinguishable single population of oscillators, which may offer insights into the unihemispheric slow-wave sleep phenomenon observed in mammals and birds.

Explosive transition in adaptive Stuart-Landau oscillators with higher-order interactions

We investigate the nature of the transition to synchronization and death states in nonidentical Stuart-Landau oscillators adaptively coupled with pairwise and higher-order interactions. Our findings reveal that nonidentical Stuart-Landau oscillators globally coupled with higher-order interactions exhibit asynchronous oscillations. In contrast, systems coupled with pairwise interactions exhibit a continuous transition to synchronization. Introducing adaptation in pairwise interactions results in an explosive transition to synchronization. Furthermore, when the system is coupled with both pairwise and higher-order interactions, adaptation leads to an explosive transition to synchronization. For higher values of coupling strength, the coupled system also exhibits an explosive transition to the death state. Extending our analysis to random networks, we observe similar results. These findings align with the behavior observed in coupled Kuramoto oscillators, offering new insights into synchronization dynamics in complex systems with adaptive coupling and higher-order interactions.

Deeper but smaller: Higher-order interactions increase linear stability but shrink basins

A key challenge of nonlinear dynamics and network science is to understand how higher-order interactions influence collective dynamics. Although many studies have approached this question through linear stability analysis, less is known about how higher-order interactions shape the global organization of different states. Here, we shed light on this issue by analyzing the rich patterns supported by identical Kuramoto oscillators on hypergraphs. We show that higher-order interactions can have opposite effects on linear stability and basin stability: They stabilize twisted states (including full synchrony) by improving their linear stability, but also make them hard to find by markedly reducing their basin size. Our results highlight the importance of understanding higher-order interactions from both local and global perspectives.

How combined pairwise and higher-order interactions shape transient dynamics

Understanding how species interactions shape biodiversity is a core challenge in ecology. While much focus has been on long-term stability, there is rising interest in transient dynamics-the short-lived periods when ecosystems respond to disturbances and adjust toward stability. These transitions are crucial for predicting ecosystem reactions and guiding effective conservation. Our study introduces a model that uses convex combinations to blend pairwise and higher-order interactions (HOIs), offering a more realistic view of natural ecosystems. We find that pairwise interactions slow the journey to stability, while HOIs speed it up. Employing global stability analysis and numerical simulations, we establish that as the proportion of HOIs increases, mean transient times exhibit a significant reduction, thereby underscoring the essential role of HOIs in enhancing biodiversity stabilization. Our results reveal a robust correlation between the most negative real part of the eigenvalues of the Jacobian matrix associated with the linearized system at the coexistence equilibrium and the mean transient times. This indicates that a more negative leading eigenvalue correlates with accelerated convergence to stable coexistence abundances. This insight is vital for comprehending ecosystem resilience and recovery, emphasizing the key role of HOIs in promoting stabilization. Amid growing interest in transient dynamics and its implications for biodiversity and ecological stability, our study enhances the understanding of how species interactions affect both transient and long-term ecosystem behavior. By addressing a critical gap in ecological theory and offering a practical framework for ecosystem management, our work advances knowledge of transient dynamics, ultimately informing effective conservation strategies.

Chimera-inspired dynamics: When higher-order interactions are expressed differently

The exploration of chimera-inspired dynamics in nonlocally coupled networks of Kuramoto oscillators with higher-order interactions is still in its nascent stages. Concurrently, the investigation of collective phenomena in higher-order interaction networks is gaining attraction. Here, we observe that hypergraph networks tend to synchronize through lower-order interactions, whereas simplicial complex networks exhibit a preference for higher-order interactions. This observation suggests that higher-order representations manifest substantial differences in chimera-inspired synchronization regions. Moreover, we introduce an explicit expression for identifying the chimera state. With a comprehensive basin stability analysis and the interplay of pairwise and higher-order interaction strengths, the emergence of the chimera state is inherent in high-order interaction networks. Our findings contribute to the understanding of chimera-inspired dynamics in higher-order interaction networks.

Tipping in Stuart-Landau oscillators induced by higher-order repulsive interactions

Tipping phenomena in complex systems represent abrupt transitions in the system behavior due to incremental changes in parameters. Here, we report the emergence of an abrupt transition from an oscillatory to a death state in coupled limit cycle oscillators with higher-order repulsive interactions. This transition contrasts with the typical continuous transitions observed with pairwise repulsive links. It is notable that higher-order interactions introduce more complexity to the system dynamics, thereby leading to the transition to the death state at lower coupling strengths compared to pairwise coupled systems. Importantly, there exists a sudden revival of the oscillation from the death state with increasing the pairwise coupling strength despite the detrimental effect of pairwise couplings. The results are also robust for nonidentical systems. Furthermore, the analytical determination of the critical coupling strength for the tipping point aligns closely with the results obtained from numerical simulations.

Higher-order interaction induced chimeralike state in a bipartite network

We report higher-order coupling induced stable chimeralike state in a bipartite network of coupled phase oscillators without any time-delay in the coupling. We show that the higher-order interaction breaks the symmetry of the homogeneous synchronized state to facilitate the manifestation of symmetry breaking chimeralike state. In particular, such symmetry breaking manifests only when the pairwise interaction is attractive and higher-order interaction is repulsive, and vice versa. Further, we also demonstrate the increased degree of heterogeneity promotes homogeneous symmetric states in the phase diagram by suppressing the asymmetric chimeralike state. We deduce the low-dimensional evolution equations for the macroscopic order parameters using Ott-Antonsen ansatz and obtain the bifurcation curves from them using the software xppaut, which agrees very well with the simulation results. We also deduce the analytical stability conditions for the incoherent state, in-phase and out-of-phase synchronized states, which match with the bifurcation curves.

Chimeric states induced by higher-order interactions in coupled prey-predator systems

Higher-order interactions have been instrumental in characterizing the intricate complex dynamics in a diverse range of large-scale complex systems. Our study investigates the effect of attractive and repulsive higher-order interactions in globally and non-locally coupled prey-predator Rosenzweig-MacArthur systems. Such interactions lead to the emergence of complex spatiotemporal chimeric states, which are otherwise unobserved in the model system with only pairwise interactions. Our model system exhibits a second-order transition from a chimera-like state (mixture of oscillating and steady state nodes) to a chimera-death state through a supercritical Hopf bifurcation. The origin of these states is discussed in detail along with the effect of the higher-order non-local topology which leads to the rise of a distinct and dynamical state termed as "amplitude-mediated chimera-like states." Our study observes that the introduction of higher-order attractive and repulsive interactions exhibit incoherence and promote persistence in consumer-resource population dynamics as opposed to susceptibility shown by synchronized dynamics with only pairwise interactions, and these results may be of interest to conservationists and theoretical ecologists studying the effect of competing interactions in ecological networks.

Complete synchronization of three-layer Rulkov neuron network coupled by electrical and chemical synapses

This paper analyzes the complete synchronization of a three-layer Rulkov neuron network model connected by electrical synapses in the same layers and chemical synapses between adjacent layers. The outer coupling matrix of the network is not Laplacian as in linear coupling networks. We develop the master stability function method, in which the invariant manifold of the master stability equations (MSEs) does not correspond to the zero eigenvalues of the connection matrix. After giving the existence conditions of the synchronization manifold about the nonlinear chemical coupling, we investigate the dynamics of the synchronization manifold, which will be identical to that of a synchronous network by fixing the same parameters and initial values. The waveforms show that the transient chaotic windows and the transient approximate periodic windows with increased or decreased periods occur alternatively before asymptotic behaviors. Furthermore, the Lyapunov exponents of the MSEs indicate that the network with a periodic synchronization manifold can achieve complete synchronization, while the network with a chaotic synchronization manifold can not. Finally, we simulate the effects of small perturbations on the asymptotic regimes and the evolution routes for the synchronous periodic and the non-synchronous chaotic network.

Phase chimera states on nonlocal hyperrings

Chimera states are dynamical states where regions of synchronous trajectories coexist with incoherent ones. A significant amount of research has been devoted to studying chimera states in systems of identical oscillators, nonlocally coupled through pairwise interactions. Nevertheless, there is increasing evidence, also supported by available data, that complex systems are composed of multiple units experiencing many-body interactions that can be modeled by using higher-order structures beyond the paradigm of classic pairwise networks. In this work we investigate whether phase chimera states appear in this framework, by focusing on a topology solely involving many-body, nonlocal, and nonregular interactions, hereby named nonlocal d-hyperring, (d+1) being the order of the interactions. We present the theory by using the paradigmatic Stuart-Landau oscillators as node dynamics, and we show that phase chimera states emerge in a variety of structures and with different coupling functions. For comparison, we show that, when higher-order interactions are "flattened" to pairwise ones, the chimera behavior is weaker and more elusive.

Effect of higher-order interactions on chimera states in two populations of Kuramoto oscillators

We investigate the effect of the fraction of pairwise and higher-order interactions on the emergent dynamics of the two populations of globally coupled Kuramoto oscillators with phase-lag parameters. We find that the stable chimera exists between saddle-node and Hopf bifurcations, while the breathing chimera lives between Hopf and homoclinic bifurcations in the two-parameter phase diagrams. The higher-order interaction facilitates the onset of the bifurcation transitions at a much lower disparity between the inter- and intra-population coupling strengths. Furthermore, the higher-order interaction facilitates the spread of breathing chimera in a large region of the parameter space while suppressing the spread of the stable chimera. A low degree of heterogeneity among the phase-lag parameters promotes the spread of both stable chimera and breathing chimera to a large region of the parameter space for a large fraction of the higher-order coupling. In contrast, a large degree of heterogeneity is found to decrease the spread of both chimera states for a large fraction of the higher-order coupling. A global synchronized state is observed above a critical value of heterogeneity among the phase-lag parameters. We have deduced the low-dimensional evolution equations for the macroscopic order parameters using the Ott-Antonsen Ansatz. We have also deduced the analytical saddle-node and Hopf bifurcation curves from the evolution equations for the macroscopic order parameters and found them to match with the bifurcation curves obtained using the software XPPAUT and with the simulation results.

The complementary contribution of each order topology into the synchronization of multi-order networks

Higher-order interactions improve our capability to model real-world complex systems ranging from physics and neuroscience to economics and social sciences. There is great interest nowadays in understanding the contribution of higher-order terms to the collective behavior of the network. In this work, we investigate the stability of complete synchronization of complex networks with higher-order structures. We demonstrate that the synchronization level of a network composed of nodes interacting simultaneously via multiple orders is maintained regardless of the intensity of coupling strength across different orders. We articulate that lower-order and higher-order topologies work together complementarily to provide the optimal stable configuration, challenging previous conclusions that higher-order interactions promote the stability of synchronization. Furthermore, we find that simply adding higher-order interactions based on existing connections, as in simple complexes, does not have a significant impact on synchronization. The universal applicability of our work lies in the comprehensive analysis of different network topologies, including hypergraphs and simplicial complexes, and the utilization of appropriate rescaling to assess the impact of higher-order interactions on synchronization stability.

First-order transition to oscillation death in coupled oscillators with higher-order interactions

We investigate the dynamical evolution of Stuart-Landau oscillators globally coupled through conjugate or dissimilar variables on simplicial complexes. We report a first-order explosive phase transition from an oscillatory state to oscillation death, with higher-order (2-simplex triadic) interactions, as opposed to the second-order transition with only pairwise (1-simplex) interactions. Moreover, the system displays four distinct homogeneous steady states in the presence of triadic interactions, in contrast to the two homogeneous steady states observed with dyadic interactions. We calculate the backward transition point analytically, confirming the numerical results and providing the origin of the dynamical states in the transition region. The results are robust against the application of noise. The study will be useful in understanding complex systems, such as ecological and epidemiological, having higher-order interactions and coupling through conjugate variables.

Collective dynamics of nonlocally coupled Hindmarsh-Rose neurons modified by magnetic flux

We investigate the dynamics of nonlocally coupled Hindmarsh-Rose neurons, modified by coupling the induced magnetic flux to the membrane potential with a quadratic memristor of strength k. The nonlocal coupling consists of the interaction of each neuron with its neighbors within a fixed radius, which influence the membrane potential of the neuron with coupling intensity σ. For such local dynamics and network of interactions, we investigate how variations of k and σ affect the collective dynamics. We find that when increasing k as well as when increasing σ, coherence typically increases, except for small ranges of these parameters where the opposite behavior can occur. Besides affecting coherence, varying k also affects the pattern of bursts and spikes, namely, for large enough k, burst frequency is augmented, the number and amplitude of the spikes are reduced, and quiescent periods become longer. Results are displayed for an intermediate range of interactions with radius 1/4 of the network size, but we also varied the range of interactions, ranging from first-neighbor to all-to-all couplings, observing in all cases a qualitatively similar impact of induction.

Higher-order interactions in Kuramoto oscillators with inertia

How do higher-order interactions influence the dynamical landscape of a network of the second-order phase oscillators? We address this question using three coupled Kuramoto phase oscillators with inertia under pairwise and higher-order interactions, in search of various collective states, and new states, if any, that show marginal presence or absence under pairwise interactions. We explore this small network for varying phase lag in the coupling and over a range of negative to positive coupling strength of pairwise as well as higher-order or group interactions. In the extended coupling parameter plane of the network we record several well-known states such as synchronization, frequency chimera states, and rotating waves that appear with distinct boundaries. In the parameter space, we also find states generated by the influence of higher-order interactions: The 2+1 antipodal point and the 2+1 phase-locked states. Our results demonstrate the importantance of the choices of the phase lag and the sign of the higher-order coupling strength for the emergent dynamics of the network. We provide analytical support to our numerical results.

Multistability in coupled oscillator systems with higher-order interactions and community structure

We study synchronization dynamics in populations of coupled phase oscillators with higher-order interactions and community structure. We find that the combination of these two properties gives rise to a number of states unsupported by either higher-order interactions or community structure alone, including synchronized states with communities organized into clusters in-phase, anti-phase, and a novel skew-phase, as well as an incoherent-synchronized state. Moreover, the system displays strong multistability with many of these states stable at the same time. We demonstrate our findings by deriving the low dimensional dynamics of the system and examining the system's bifurcations using stability analysis and perturbation theory.

Multiple first-order transitions in simplicial complexes on multilayer systems

The presence of higher-order interactions (simplicial complexes) on globally coupled systems yield abrupt first-order transitions to synchronization. We discover that simplicial complexes on multilayer systems can yield multiple basins of attraction, leading to multiple abrupt first-order transitions to (de)synchronization for associated coupled dynamics. Using the Ott-Antonsen approach, we develop an analytical framework for simplicial complexes on multilayer systems, reducing the high-dimensional evolution equation to a low-dimensional manifold, which thoroughly explains the origin and stability of all possible dynamical states, including multiple synchronization transitions. The study illustrating rich dynamical behaviors could be pivotal in comprehending the impacts of higher-order interactions on dynamics of complex real-world networks, such as brain, social, and technological, which have inherent multilayer architectures.

Stability of synchronization in simplicial complexes with multiple interaction layers

Understanding how the interplay between higher-order and multilayer structures of interconnections influences the synchronization behaviors of dynamical systems is a feasible problem of interest, with possible application in essential topics such as neuronal dynamics. Here, we provide a comprehensive approach for analyzing the stability of the complete synchronization state in simplicial complexes with numerous interaction layers. We show that the synchronization state exists as an invariant solution and derive the necessary condition for a stable synchronization state in the presence of general coupling functions. It generalizes the well-known master stability function scheme to the higher-order structures with multiple interaction layers. We verify our theoretical results by employing them on networks of paradigmatic Rössler oscillators and Sherman neuronal models, and we demonstrate that the presence of group interactions considerably improves the synchronization phenomenon in the multilayer framework.