Phase reduction beyond the first order: The case of the mean-field complex Ginzburg-Landau equation

Rotating clusters in phase-lagged Kuramoto oscillators with higher-order interactions

The effect of phase-lag in pairwise interactions has been a topic of great interest for a while. However, real-world systems often have interactions that are beyond pairwise and can be modeled using simplicial complexes. We show that the inclusion of higher-order interactions in phase-lagged coupled Kuramoto oscillators shifts the critical point at which first-order transition from a cluster synchronized state to an incoherent state takes place. Considering the polar coordinates, we obtain the rotation frequency of the clusters, which turns out to be a function of the phase-lag parameter. In turn, the phase- lag can be used as a control parameter to achieve a desired cluster frequency. Moreover, in the thermodynamic limit, by employing the Ott-Antonsen approach we derive a reduced equation for the order parameter measuring cluster synchronization and progress further through the self-consistency method to obtain a closed form of the order parameter measuring global synchronization which was lacking in the Ott-Antonsen approach.

Prolonged hysteresis in the Kuramoto model with inertia and higher-order interactions

The inclusion of inertia in the Kuramoto model has long been reported to change the nature of a phase transition, providing a fertile ground to model the dynamical behaviors of interacting units. More recently, higher-order interactions have been realized as essential for the functioning of real-world complex systems ranging from the brain to disease spreading. Yet analytical insights to decipher the role of inertia with higher-order interactions remain challenging. Here, we study the Kuramoto model with inertia on simplicial complexes, merging two research domains. We develop an analytical framework in a mean-field setting using self-consistent equations to describe the steady-state behavior, which reveals a prolonged hysteresis in the synchronization profile. Inertia and triadic interaction strength exhibit isolated influence on system dynamics by predominantly governing, respectively, the forward and backward transition points. This paper sets a paradigm to deepen our understanding of real-world complex systems such as power grids modeled as the Kuramoto model with inertia.

Determining bifurcations to explosive synchronization for networks of coupled oscillators with higher-order interactions

We determine bifurcations from gradual to explosive synchronization in coupled oscillator networks with higher-order coupling using self-consistency analysis. We obtain analytic bifurcation values for generic symmetric natural frequency distributions. We show that nonsynchronized, drifting, oscillators are non-negligible and play a crucial role in bifurcation. As such, the entire natural frequency distribution must be accounted for, rather than just the shape at the center. We verify our results for Lorentzian- and Gaussian-distributed natural frequencies.

Insights into oscillator network dynamics using a phase-isostable framework

Networks of coupled nonlinear oscillators can display a wide range of emergent behaviors under the variation of the strength of the coupling. Network equations for pairs of coupled oscillators where the dynamics of each node is described by the evolution of its phase and slowest decaying isostable coordinate have previously been shown to capture bifurcations and dynamics of the network, which cannot be explained through standard phase reduction. An alternative framework using isostable coordinates to obtain higher-order phase reductions has also demonstrated a similar descriptive ability for two oscillators. In this work, we consider the phase-isostable network equations for an arbitrary but finite number of identical coupled oscillators, obtaining conditions required for the stability of phase-locked states including synchrony. For the mean-field complex Ginzburg-Landau equation where the solutions of the full system are known, we compare the accuracy of the phase-isostable network equations and higher-order phase reductions in capturing bifurcations of phase-locked states. We find the former to be the more accurate and, therefore, employ this to investigate the dynamics of globally linearly coupled networks of Morris-Lecar neuron models (both two and many nodes). We observe qualitative correspondence between results from numerical simulations of the full system and the phase-isostable description demonstrating that in both small and large networks, the phase-isostable framework is able to capture dynamics that the first-order phase description cannot.

Higher-order interactions induce anomalous transitions to synchrony

We analyze the simplest model of identical coupled phase oscillators subject to two-body and three-body interactions with permutation symmetry and phase lags. This model is derived from an ensemble of weakly coupled nonlinear oscillators by phase reduction, where the first and second harmonic interactions with phase lags naturally appear. Our study indicates that the higher-order interactions induce anomalous transitions to synchrony. Unlike the conventional Kuramoto model, higher-order interactions lead to anomalous phenomena such as multistability of full synchronization, incoherent, and two-cluster states, and transitions to synchrony through slow switching and clustering. Phase diagrams of the dynamical regimes are constructed theoretically and verified by direct numerical simulations. We also show that similar transition scenarios are observed even if a small heterogeneity in the oscillators' frequency is included.

High-order phase reduction for coupled 2D oscillators

Phase reduction is a general approach to describe coupled oscillatory units in terms of their phases, assuming that the amplitudes are enslaved. The coupling should be small for such reduction, but one also expects the reduction to be valid for finite coupling. This paper presents a general framework, allowing us to obtain coupling terms in higher orders of the coupling parameter for generic two-dimensional oscillators and arbitrary coupling terms. The theory is illustrated with an accurate prediction of Arnold's tongue for the van der Pol oscillator exploiting higher-order phase reduction.

Synchronization transitions in Kuramoto networks with higher-mode interaction

Synchronization is an omnipresent collective phenomenon in nature and technology, whose understanding is still elusive for real-world systems in particular. We study the synchronization transition in a phase oscillator system with two nonvanishing Fourier-modes in the interaction function, hence going beyond the Kuramoto paradigm. We show that the transition scenarios crucially depend on the interplay of the two coupling modes. We describe the multistability induced by the presence of a second coupling mode. By extending the collective coordinate approach, we describe the emergence of various states observed in the transition from incoherence to coherence. Remarkably, our analysis suggests that, in essence, the two-mode coupling gives rise to states characterized by two independent but interacting groups of oscillators. We believe that these findings will stimulate future research on dynamical systems, including complex interaction functions beyond the Kuramoto-type.

Tiered synchronization in Kuramoto oscillators with adaptive higher-order interactions

Phase transitions widely occur in natural systems. Incorporation of higher-order interactions in coupled dynamics is known to cause first-order phase transition to synchronization in an otherwise smooth second-order in the presence of only pairwise interactions. Here, we discover that adaptation in higher-order interactions restores the second-order phase transition in the former setup and notably produces additional bifurcation referred as tiered synchronization as a consequence of combination of super-critical pitchfork and two saddle node bifurcations. The Ott-Antonsen manifold underlines the interplay of higher-order interactions and adaptation in instigating tiered synchronization, as well as provides complete description of all (un)stable states. These results would be important in comprehending dynamics of real-world systems with inherent higher-order interactions and adaptation through feedback coupling.

Higher-order interactions shape collective dynamics differently in hypergraphs and simplicial complexes

Higher-order networks have emerged as a powerful framework to model complex systems and their collective behavior. Going beyond pairwise interactions, they encode structured relations among arbitrary numbers of units through representations such as simplicial complexes and hypergraphs. So far, the choice between simplicial complexes and hypergraphs has often been motivated by technical convenience. Here, using synchronization as an example, we demonstrate that the effects of higher-order interactions are highly representation-dependent. In particular, higher-order interactions typically enhance synchronization in hypergraphs but have the opposite effect in simplicial complexes. We provide theoretical insight by linking the synchronizability of different hypergraph structures to (generalized) degree heterogeneity and cross-order degree correlation, which in turn influence a wide range of dynamical processes from contagion to diffusion. Our findings reveal the hidden impact of higher-order representations on collective dynamics, highlighting the importance of choosing appropriate representations when studying systems with nonpairwise interactions.

Solvable Dynamics of Coupled High-Dimensional Generalized Limit-Cycle Oscillators

We introduce a new model consisting of globally coupled high-dimensional generalized limit-cycle oscillators, which explicitly incorporates the role of amplitude dynamics of individual units in the collective dynamics. In the limit of weak coupling, our model reduces to the D-dimensional Kuramoto phase model, akin to a similar classic construction of the well-known Kuramoto phase model from weakly coupled two-dimensional limit-cycle oscillators. For the practically important case of D=3, the incoherence of the model is rigorously proved to be stable for negative coupling (K<0) but unstable for positive coupling (K>0); the locked states are shown to exist if K>0; in particular, the onset of amplitude death is theoretically predicted. For D≥2, the discrete and continuous spectra for both locked states and amplitude death are governed by two general formulas. Our proposed D-dimensional model is physically more reasonable, because it is no longer constrained by fixed amplitude dynamics, which puts the recent studies of the D-dimensional Kuramoto phase model on a stronger footing by providing a more general framework for D-dimensional limit-cycle oscillators.

Synchronization of phase oscillators on complex hypergraphs

We study the effect of structured higher-order interactions on the collective behavior of coupled phase oscillators. By combining a hypergraph generative model with dimensionality reduction techniques, we obtain a reduced system of differential equations for the system's order parameters. We illustrate our framework with the example of a hypergraph with hyperedges of sizes 2 (links) and 3 (triangles). For this case, we obtain a set of two coupled nonlinear algebraic equations for the order parameters. For strong values of coupling via triangles, the system exhibits bistability and explosive synchronization transitions. We find conditions that lead to bistability in terms of hypergraph properties and validate our predictions with numerical simulations. Our results provide a general framework to study the synchronization of phase oscillators in hypergraphs, and they can be extended to hypergraphs with hyperedges of arbitrary sizes, dynamic-structural correlations, and other features.

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.

Efficient moment-based approach to the simulation of infinitely many heterogeneous phase oscillators

The dynamics of ensembles of phase oscillators are usually described considering their infinite-size limit. In practice, however, this limit is fully accessible only if the Ott-Antonsen theory can be applied, and the heterogeneity is distributed following a rational function. In this work, we demonstrate the usefulness of a moment-based scheme to reproduce the dynamics of infinitely many oscillators. Our analysis is particularized for Gaussian heterogeneities, leading to a Fourier-Hermite decomposition of the oscillator density. The Fourier-Hermite moments obey a set of hierarchical ordinary differential equations. As a preliminary experiment, the effects of truncating the moment system and implementing different closures are tested in the analytically solvable Kuramoto model. The moment-based approach proves to be much more efficient than the direct simulation of a large oscillator ensemble. The convenience of the moment-based approach is exploited in two illustrative examples: (i) the Kuramoto model with bimodal frequency distribution, and (ii) the "enlarged Kuramoto model" (endowed with nonpairwise interactions). In both systems, we obtain new results inaccessible through direct numerical integration of populations.

Tiered synchronization in coupled oscillator populations with interaction delays and higher-order interactions

We study synchronization in large populations of coupled phase oscillators with time delays and higher-order interactions. With each of these effects individually giving rise to bistability between incoherence and synchronization via subcriticality at the onset of synchronization and the development of a saddle node, we find that their combination yields another mechanism behind bistability, where supercriticality at onset may be maintained; instead, the formation of two saddle nodes creates tiered synchronization, i.e., bistability between a weakly synchronized state and a strongly synchronized state. We demonstrate these findings by first deriving the low dimensional dynamics of the system and examining the system bifurcations using a stability and steady-state analysis.

Enlarged Kuramoto model: Secondary instability and transition to collective chaos

The emergence of collective synchrony from an incoherent state is a phenomenon essentially described by the Kuramoto model. This canonical model was derived perturbatively, by applying phase reduction to an ensemble of heterogeneous, globally coupled Stuart-Landau oscillators. This derivation neglects nonlinearities in the coupling constant. We show here that a comprehensive analysis requires extending the Kuramoto model up to quadratic order. This "enlarged Kuramoto model" comprises three-body (nonpairwise) interactions, which induce strikingly complex phenomenology at certain parameter values. As the coupling is increased, a secondary instability renders the synchronized state unstable, and subsequent bifurcations lead to collective chaos. An efficient numerical study of the thermodynamic limit, valid for Gaussian heterogeneity, is carried out by means of a Fourier-Hermite decomposition of the oscillator density.

Collective dynamics of phase oscillator populations with three-body interactions

Many-body interactions between dynamical agents have caught particular attention in recent works that found wide applications in physics, neuroscience, and sociology. In this paper we investigate such higher order (nonadditive) interactions on collective dynamics in a system of globally coupled heterogeneous phase oscillators. We show that the three-body interactions encoded microscopically in nonlinear couplings give rise to added dynamic phenomena occurring beyond the pairwise interactions. The system in general displays an abrupt desynchronization transition characterized by irreversible explosive synchronization via an infinite hysteresis loop. More importantly, we give a mathematical argument that such an abrupt dynamic pattern is a universally expected effect. Furthermore, the origin of this abrupt transition is uncovered by performing a rigorous stability analysis of the equilibrium states, as well as by providing a detailed description of the spectrum structure of linearization around the steady states. Our work reveals a self-organized phenomenon that is responsible for the rapid switching to synchronization in diverse complex systems exhibiting critical transitions with nonpairwise interactions.

Dead zones and phase reduction of coupled oscillators

A dead zone in the interaction between two dynamical systems is a region of their joint phase space where one system is insensitive to the changes in the other. These can arise in a number of contexts, and their presence in phase interaction functions has interesting dynamical consequences for the emergent dynamics. In this paper, we consider dead zones in the interaction of general coupled dynamical systems. For weakly coupled limit cycle oscillators, we investigate criteria that give rise to dead zones in the phase interaction functions. We give applications to coupled multiscale oscillators where coupling on only one branch of a relaxation oscillation can lead to the appearance of dead zones in a phase description of their interaction.

Collective behaviors of mean-field coupled Stuart-Landau limit-cycle oscillators under additional repulsive links

We study the collective behaviors of a large population of Stuart-Landau limit-cycle oscillators that coupled diffusively and equally with all of the others via the conjugate of the mean field, where the underlying interaction is shown to break the rotational symmetry of the coupled system. In the model, an ensemble of Stuart-Landau oscillators are in fact diffusively coupled via the mean field in the real parts, whereas additional repulsive links are present in the imaginary parts. All the oscillators are linked via the similar state variables, which distinctly differs from the conjugate coupling through dissimilar variables in the previous studies. We show that depending on the strength of coupling and the distribution of natural frequencies, the coupled system exhibits three qualitatively different types of collective stationary behaviors: amplitude death (AD), oscillation death (OD), and incoherent state. Our goal is to analytically characterize the onset of the above three typical macrostates by performing the rigorous linear stability analyses of the corresponding mean-field coupled system. We prove that AD is able to occur in the coupled system with identical frequencies, where the stable coupling interval of AD is found to be independent on the system's size N. When the natural frequencies are distributed according to a general density function, we obtain the analytic equations that govern the exact stability boundaries of AD, OD, and the incoherence for a coupled system in the thermodynamic limit N→∞. All the theoretical predictions are well confirmed via numerical simulations of the coupled system with a specific Lorentzian frequency distribution.

A universal route to explosive phenomena

Critical transitions are observed in many complex systems. This includes the onset of synchronization in a network of coupled oscillators or the emergence of an epidemic state within a population. "Explosive" first-order transitions have caught particular attention in a variety of systems when classical models are generalized by incorporating additional effects. Here, we give a mathematical argument that the emergence of these first-order transitions is not surprising but rather a universally expected effect: Varying a classical model along a generic two-parameter family must lead to a change of the criticality. To illustrate our framework, we give three explicit examples of the effect in distinct physical systems: a model of adaptive epidemic dynamics, for a generalization of the Kuramoto model, and for a percolation transition.

Traveling waves in non-local pulse-coupled networks

Traveling phase waves are commonly observed in recordings of the cerebral cortex and are believed to organize behavior across different areas of the brain. We use this as motivation to analyze a one-dimensional network of phase oscillators that are nonlocally coupled via the phase response curve (PRC) and the Dirac delta function. Existence of waves is proven and the dispersion relation is computed. Using the theory of distributions enables us to write and solve an associated stability problem. First and second order perturbation theory is applied to get analytic insight and we show that long waves are stable while short waves are unstable. We apply the results to PRCs that come from mitral neurons. We extend the results to smooth pulse-like coupling by reducing the nonlocal equation to a local one and solving the associated boundary value problem.

Discrete Polynomial Optimization with Coherent Networks of Condensates and Complex Coupling Switching

Gain-dissipative platforms consisting of lasers, optical parametric oscillators and nonequilibrium condensates operating at the condensation or coherence threshold have been recently proposed as efficient analog simulators of the two-local spin Hamiltonians with continuous or discrete degrees of freedom. We show that nonequilibrium condensates above the threshold arranged in an interacting network may realize k-local Hamiltonians with k>2 and lead to nontrivial phase configurations. Similarly, many gain-dissipative systems that can be manipulated by optical means can bring about the ground state of the k-local Hamiltonians and solve higher-order binary optimization problems. We show how to facilitate the search for the global solution by invoking complex couplings in the system and demonstrate the efficiency of the method on the sets of complex problems. This approach offers a highly flexible new kind of computation based on gain-dissipative simulators with complex coupling switching.

Comment on "The Winfree model with non-infinitesimal phase-response curve: Ott-Antonsen theory" [Chaos 30, 073139 (2020)]

In a recent paper [Chaos 30, 073139 (2020)], we analyzed an extension of the Winfree model with nonlinear interactions. The nonlinear coupling function Q was mistakenly identified with the non-infinitesimal phase-response curve (PRC). Here, we assess to what extent Q and the actual PRC differ in practice. By means of numerical simulations, we compute the PRCs corresponding to the Q functions previously considered. The results confirm a qualitative similarity between the PRC and the coupling function Q in all cases.

Quasi phase reduction of all-to-all strongly coupled λ-ω oscillators near incoherent states

The dynamics of an ensemble of N weakly coupled limit-cycle oscillators can be captured by their N phases using standard phase reduction techniques. However, it is a phenomenological fact that all-to-all strongly coupled limit-cycle oscillators may behave as "quasiphase oscillators," evidencing the need of novel reduction strategies. We introduce, here, quasi phase reduction (QPR), a scheme suited for identical oscillators with polar symmetry (λ-ω systems). By applying QPR, we achieve a reduction to N+2 degrees of freedom: N phase oscillators interacting through one independent complex variable. This "quasi phase model" is asymptotically valid in the neighborhood of incoherent states, irrespective of the coupling strength. The effectiveness of QPR is illustrated in a particular case, an ensemble of Stuart-Landau oscillators, obtaining exact stability boundaries of uniform and nonuniform incoherent states for a variety of couplings. An extension of QPR beyond the neighborhood of incoherence is also explored. Finally, a general QPR model with N+2M degrees of freedom is obtained for coupling through the first M harmonics.

Understanding the dynamics of biological and neural oscillator networks through exact mean-field reductions: a review

Many biological and neural systems can be seen as networks of interacting periodic processes. Importantly, their functionality, i.e., whether these networks can perform their function or not, depends on the emerging collective dynamics of the network. Synchrony of oscillations is one of the most prominent examples of such collective behavior and has been associated both with function and dysfunction. Understanding how network structure and interactions, as well as the microscopic properties of individual units, shape the emerging collective dynamics is critical to find factors that lead to malfunction. However, many biological systems such as the brain consist of a large number of dynamical units. Hence, their analysis has either relied on simplified heuristic models on a coarse scale, or the analysis comes at a huge computational cost. Here we review recently introduced approaches, known as the Ott-Antonsen and Watanabe-Strogatz reductions, allowing one to simplify the analysis by bridging small and large scales. Thus, reduced model equations are obtained that exactly describe the collective dynamics for each subpopulation in the oscillator network via few collective variables only. The resulting equations are next-generation models: Rather than being heuristic, they exactly link microscopic and macroscopic descriptions and therefore accurately capture microscopic properties of the underlying system. At the same time, they are sufficiently simple to analyze without great computational effort. In the last decade, these reduction methods have become instrumental in understanding how network structure and interactions shape the collective dynamics and the emergence of synchrony. We review this progress based on concrete examples and outline possible limitations. Finally, we discuss how linking the reduced models with experimental data can guide the way towards the development of new treatment approaches, for example, for neurological disease.

Phase-amplitude reduction far beyond the weakly perturbed paradigm

While phase reduction is a well-established technique for the analysis of perturbed limit cycle oscillators, practical application requires perturbations to be sufficiently weak thereby limiting its utility in many situations. Here, a general strategy is developed for constructing a set of phase-amplitude reduced equations that is valid to arbitrary orders of accuracy in the amplitude coordinates. This reduction framework can be used to investigate the behavior of oscillatory dynamical systems far beyond the weakly perturbed paradigm. Additionally, a patchwork phase-amplitude reduction method is suggested that is useful when exceedingly large magnitude perturbations are considered. This patchwork method incorporates the high-accuracy phase-amplitude reductions of multiple nearby periodic orbits that result from modifications to nominal parameters. The proposed method of high-accuracy phase-amplitude reduction can be readily implemented numerically and examples are provided where reductions are computed up to fourteenth order accuracy.