Dynamics of phase oscillators with generalized frequency-weighted coupling

Higher-order interactions induce stepwise explosive phase transitions

In recent studies, it has been established that higher-order interactions in coupled oscillators can induce a process from continuous to explosive phase transition. In this study, we identify a phase transition, termed the stepwise explosive phase transition, characterized by the emergence of multiple critical phase plateaus in a globally frequency-weighted coupled pendulum model. This transition bridges the continuous and explosive phase transitions, arising from a delicate balance between attractive higher-order interactions and repulsive pairwise interactions. Specifically, the stepwise explosive phase transition occurs when the higher-order coupling is moderate, neither large nor small, while the pairwise interactions remain repulsive. Our analysis shows that stronger attractive higher-order interactions necessitate weaker repulsive pairwise interactions, leading to partial frequency locking among oscillators and triggering the stepwise transition. We construct an analytical framework using self-consistent equations to provide an approximation of the steady-state behavior. This study uncovers an alternative pathway to the desynchronization, and it provides additional insights into the phase transitions in coupled dynamical networks.

Solvable dynamics of the three-dimensional Kuramoto model with frequency-weighted coupling

The presence of coupling heterogeneity is deemed to be a natural attribute in realistic systems comprised of many interacting agents. In this work, we study dynamics of the 3D Kuramoto model with heterogeneous couplings, where the strength of coupling for each agent is weighted by its intrinsic rotation frequency. The critical coupling strength for the instability of incoherence is rigorously derived to be zero by carrying out a linear stability analysis of an incoherent state. For positive values of the coupling strength, at which the incoherence turns out to be unstable, a self-consistency approach is developed to theoretically predict the degree of global coherence of the model. Our theoretical analyses match well with numerical simulations, which helps us to deepen the understanding of collective behaviors spontaneously emerged in heterogeneously coupled high-dimensional dynamical networks.

Relaxation dynamics of phase oscillators with generic heterogeneous coupling

The coupled phase oscillator model serves as a paradigm that has been successfully used to shed light on the collective dynamics occurring in large ensembles of interacting units. It was widely known that the system experiences a continuous (second-order) phase transition to synchronization by gradually increasing the homogeneous coupling among the oscillators. As the interest in exploring synchronized dynamics continues to grow, the heterogeneous patterns between phase oscillators have received ample attention during the past years. Here, we consider a variant of the Kuramoto model with quenched disorder in their natural frequencies and coupling. Correlating these two types of heterogeneity via a generic weighted function, we systematically investigate the impacts of the heterogeneous strategies, the correlation function, and the natural frequency distribution on the emergent dynamics. Importantly, we develop an analytical treatment for capturing the essential dynamical properties of the equilibrium states. In particular, we uncover that the critical threshold corresponding to the onset of synchronization is unaffected by the location of the inhomogeneity, which, however, does depend crucially on the value of the correlation function at its center. Furthermore, we reveal that the relaxation dynamics of the incoherent state featuring the responses to external perturbations is significantly shaped by all the considered effects, thereby leading to various decaying mechanisms of the order parameters in the subcritical region. Moreover, we untangle that synchronization is facilitated by the out-coupling strategy in the supercritical region. Our study is a step forward in highlighting the potential importance of the inhomogeneous patterns involved in the complex systems, and could thus provide theoretical insights for profoundly understanding the generic statistical mechanical properties of the steady states toward synchronization.

Route to synchronization in coupled phase oscillators with frequency-dependent coupling: Explosive or continuous?

Interconnected dynamical systems often transition between states of incoherence and synchronization due to changes in system parameters. These transitions could be continuous (gradual) or explosive (sudden) and may result in failures, which makes determining their nature important. In this study, we abstract dynamical networks as an ensemble of globally coupled Kuramoto-like phase oscillators with frequency-dependent coupling and investigate the mechanisms for transition between incoherent and synchronized dynamics. The characteristics that dictate a continuous or explosive route to synchronization are the distribution of the natural frequencies of the oscillators, quantified by a probability density function g(ω), and the relation between the coupling strength and natural frequency of an oscillator, defined by a frequency-coupling strength correlation function f(ω). Our main results are conditions on f(ω) and g(ω) that result in continuous or explosive routes to synchronization and explain the underlying physics. The analytical developments are validated through numerical examples.

Partial locking in phase-oscillator populations with heterogenous coupling

We consider a variant of the mean-field model of coupled phase oscillators with uniform distribution of natural frequencies. By establishing correlations between the quenched disorder of intrinsic frequencies and coupling strength with both in- and out-coupling heterogeneities, we reveal a generic criterion for the onset of partial locking that takes place in a domain with the coexistence of phase-locked oscillators and drifters. The critical points manifesting the instability of the stationary states are obtained analytically. In particular, the bifurcation mechanism of the equilibrium states is uncovered by the use of frequency-dependent version of the Ott-Antonsen reduction consistently with the analysis based on the self-consistent approach. We demonstrate that both the manner of coupling heterogeneity and correlation exponent have influence on the emergent patterns of partial locking. Our research could find applicability in better understanding the phase transitions and related collective phenomena involving synchronization control in networked systems.

Universal scaling and phase transitions of coupled phase oscillator populations

The Kuramoto model, which serves as a paradigm for investigating synchronization phenomena of oscillatory systems, is known to exhibit second-order, i.e., continuous, phase transitions in the macroscopic order parameter. Here we generalize a number of classical results by presenting a general framework for analytically capturing the critical scaling properties of the order parameter at the onset of synchronization. Using a self-consistent approach and constructing a characteristic function, we identify various phase transitions toward synchrony and establish scaling relations describing the asymptotic dependence of the order parameter on the coupling strength near the critical point. We find that the geometric properties of the characteristic function, which depends on the natural frequency distribution, determine the scaling properties of order parameter above criticality.

Dynamics of the generalized Kuramoto model with nonlinear coupling: Bifurcation and stability

We systematically study dynamics of a generalized Kuramoto model of globally coupled phase oscillators. The coupling of modified model depends on the fraction of phase-locked oscillators via a power-law function of the Kuramoto order parameter r through an exponent α, such that α=1 corresponds to the standard Kuramoto model, α<1 strengthens the global coupling, and the global coupling is weakened if α>1. With a self-consistency approach, we demonstrate that bifurcation diagrams of synchronization for different values of α are thoroughly constructed from two parametric equations. In contrast to the case of α=1 with a typical second-order phase transition to synchronization, no phase transition to synchronization is predicted for α<1, as the onset of partial locking takes place once the coupling strength K>0. For α>1, we establish an abrupt desynchronization transition from the partially (fully) locked state to the incoherent state, whereas there is no counterpart of abrupt synchronization transition from incoherence to coherence due to that the incoherent state remains linearly neutrally stable for all K>0. For each case of α, by performing a standard linear stability analysis for the reduced system with Ott-Antonsen ansatz, we analytically derive the continuous and discrete spectra of both the incoherent state and the partially (fully) locked states. All our theoretical results are obtained in the thermodynamic limit, which have been well validated by extensive numerical simulations of the phase-model with a sufficiently large number of oscillators.