No Arabic abstract
Globally coupled ensembles of phase oscillators serve as useful tools for modeling synchronization and collective behavior in a variety of applications. As interest in the effects of simplicial interactions (i.e., non-additive, higher-order interactions between three or more units) continues to grow we study an extension of the Kuramoto model where oscillators are coupled via three-way interactions that exhibits novel dynamical properties including clustering, multistability, and abrupt desynchronization transitions. Here we provide a rigorous description of the stability of various multicluster states by studying their spectral properties in the thermodynamic limit. Not unlike the classical Kuramoto model, a natural frequency distribution with infinite support yields a population of drifting oscillators, which in turn guarantees that a portion of the spectrum is located on the imaginary axes, resulting in neutrally stable or unstable solutions. On the other hand, a natural frequency distribution with finite support allows for a fully phase-locked state, whose spectrum is real and may be linearly stable or unstable.
In the context of the celebrated Kuramoto model of globally-coupled phase oscillators of distributed natural frequencies, which serves as a paradigm to investigate spontaneous collective synchronization in many-body interacting systems, we report on a very rich phase diagram in presence of thermal noise and an additional non-local interaction on a one-dimensional periodic lattice. Remarkably, the phase diagram involves both equilibrium and non-equilibrium phase transitions. In two contrasting limits of the dynamics, we obtain exact analytical results for the phase transitions. These two limits correspond to (i) the absence of thermal noise, when the dynamics reduces to that of a non-linear dynamical system, and (ii) the oscillators having the same natural frequency, when the dynamics becomes that of a statistical system in contact with a heat bath and relaxing to a statistical equilibrium state. In the former case, our exact analysis is based on the use of the so-called Ott-Antonsen ansatz to derive a reduced set of nonlinear partial differential equations for the macroscopic evolution of the system. Our results for the case of statistical equilibrium are on the other hand obtained by extending the well-known transfer matrix approach for nearest-neighbor Ising model to consider non-local interactions. The work offers a case study of exact analysis in many-body interacting systems. The results obtained underline the crucial role of additional non-local interactions in either destroying or enhancing the possibility of observing synchrony in mean-field systems exhibiting spontaneous synchronization.
We study the dynamics of coupled oscillator networks with higher-order interactions and their ability to store information. In particular, the fixed points of these oscillator systems consist of two clusters of oscillators that become entrained at opposite phases, mapping easily to information more commonly represented by sequences of 0s and 1s. While $2^N$ such fixed point states exist in a system of $N$ oscillators, we find that a relatively small fraction of these are stable, as chosen by the network topology. To understand the memory selection of such oscillator networks, we derive a stability criterion to identify precisely which states are stable, i.e., which pieces of information are supported by the network. We also investigate the process by which the system can switch between different stable states when a random perturbation is applied that may force the system into the basin of attraction of another stable state.
We study a variant of Kuramoto-Sakaguchi model in which oscillators are divided into two groups, each characterized by its coupling constant and phase lag. Specifically, we consider the case that one coupling constant is positive and the other negative, and calculate numerically the traveling speed of two clusters emerging in the system and average separation between them as well as the order parameters for positive and negative oscillators, as the two coupling constants, phase lags, and the fraction of positive oscillators are varied. An expression explaining the dependence of the traveling speed on these parameters is obtained and observed to fit well the numerical data. With the help of this, we describe the conditions for the traveling state to appear in the system.
We study the effects of Janus oscillators in a system of phase oscillators in which the coupling constants take both positive and negative values. Janus oscillators may also form a cluster when the other ones are ordered and we calculate numerically the traveling speed of three clusters emerging in the system and average separations between them as well as the order parameters for three groups of oscillators, as the coupling constants and the fractions of positive and Janus oscillators are varied. An expression explaining the dependence of the traveling speed on these parameters is obtained and observed to fit well the numerical data. With the help of this, we describe how Janus oscillators affect the traveling of the clusters in the system.
Collective behavior plays a key role in the function of a wide range of physical, biological, and neurological systems where empirical evidence has recently uncovered the prevalence of higher-order interactions, i.e., structures that represent interactions between more than just two individual units, in complex network structures. Here, we study the optimization of collective behavior in networks with higher-order interactions encoded in clique complexes. Our approach involves adapting the Synchrony Alignment Function framework to a new composite Laplacian matrix that encodes multi-order interactions including, e.g., both dyadic and triadic couplings. We show that as higher-order coupling interactions are equitably strengthened, so that overall coupling is conserved, the optimal collective behavior improves. We find that this phenomenon stems from the broadening of a composite Laplacians eigenvalue spectrum, which improves the optimal collective behavior and widens the range of possible behaviors. Moreover, we find in constrained optimization scenarios that a nontrivial, ideal balance between the relative strengths of pair-wise and higher-order interactions leads to the strongest collective behavior supported by a network. This work provides insight into how systems balance interactions of different types to optimize or broaden their dynamical range of behavior, especially for self-regulating systems like the brain.