Mesoscopic model reduction for the collective dynamics of sparse coupled oscillator networks

TL;DR

This paper develops two mesoscopic model reductions based on collective coordinates to accurately and efficiently describe the bifurcation from global to partial synchronization in finite sparse oscillator networks, overcoming limitations of standard methods.

Contribution

The paper introduces two novel mesoscopic reductions that improve modeling of bifurcation dynamics in finite sparse networks, with calibration and reduced computational complexity.

Findings

Accurate modeling of bifurcation dynamics in sparse networks.

Reduction of computational complexity from O(N^2) to O(1).

Excellent agreement with numerical simulations near bifurcation.

Abstract

The behavior at bifurcation from global synchronization to partial synchronization in finite networks of coupled oscillators is a complex phenomenon, involving the intricate dynamics of one or more oscillators with the remaining synchronized oscillators. This is not captured well by standard macroscopic model reduction techniques which capture only the collective behavior of synchronized oscillators in the thermodynamic limit. We introduce two mesoscopic model reductions for finite sparse networks of coupled oscillators to quantitatively capture the dynamics close to bifurcation from global to partial synchronization. Our model reduction builds upon the method of collective coordinates. We first show that standard collective coordinate reduction has difficulties capturing this bifurcation. We identify a particular topological structure at bifurcation consisting of a main synchronized cluster, the oscillator that desynchronizes at bifurcation, and an intermediary node connecting them. Utilizing this structure and ensemble averages we derive an analytic expression for the mismatch between the true bifurcation from global to partial synchronization and its estimate calculated via the collective coordinate approach. This allows to calibrate the standard collective coordinate approach without prior knowledge of which node will desynchronize. We introduce a second mesoscopic reduction, utilizing the same particular topological structure, which allows for a quantitative dynamical description of the phases near bifurcation. The mesoscopic reductions significantly reduce the computational complexity of the collective coordinate approach, reducing from $\mathcal{O}(N^2)$ to $\mathcal{O}(1)$. We perform numerical simulations for Erd\H{o}s-R\'enyi networks and for modified Barab\'asi-Albert networks demonstrating excellent quantitative agreement at and close to bifurcation.