Dynamics of matrix coupled Kuramoto oscillators on modular networks: excitable behavior and global decoherence

Guilherme S. Costa,Marcus A. M. de Aguiar

2024-06-13

Abstract:Synchronization is observed in many natural systems, with examples ranging from neuronal activation to walking pedestrians. The models proposed by Winfree and Kuramoto stand as the classic frameworks for investigating these phenomena. The Kuramoto model, in particular, has been extended in different ways since its original formulation to account for more general scenarios. One such extension replaces the coupling parameter with a coupling matrix, describing a form of generalized frustration with broken rotational symmetry. A key feature of this model is the existence of {\it phase tuned states}, characterized by having the phase of the order parameter pointing in the direction of the dominant eigenvector of the coupling matrix. Here we investigate the matrix coupled Kuramoto model on networks with two modules, such that one module is in the phase tuned state and the other in a state where the order parameter rotates. We identified different regimes in which one or the other module dominates the dynamics. We found, in particular, that the phase tuned module can create a bottleneck for the oscillation of the rotating module, leading to a behavior similar to the charge and fire regimes of excitable systems. We also found an extended region in the parameter space where motion is globally disordered, even though one of the modules presented high levels of synchronization when uncoupled.

Physics and Society,Pattern Formation and Solitons

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the dynamic behavior of matrix - coupled Kuramoto oscillators in modular networks. In particular, when the network contains two modules, one of which is in a phase - tuning state and the ordered parameter of the other module rotates, the dynamic behavior of the system. Specifically, the researchers focus on the following aspects:

1. **Identification of different dynamic regions**: By changing the inter - module coupling strength \(K_{12}\) and the phase - lag parameter \(\alpha_1\), the researchers identify the dynamic behaviors of the system under different parameters. These behaviors include:
- **Independent module region**: The two modules operate almost independently.
- **G1 - dominated region**: The G1 module dominates, and the two modules rotate at similar frequencies.
- **G2 - dominated region**: The G2 module dominates, and the phases of the two modules are tuned to a specific direction.
- **Asynchronous region**: The interaction between the modules leads to global disorder.
- **Exciting non - linear oscillation region**: The modules exhibit charge - discharge behaviors similar to excitatory neurons.

2. **Boundary between synchronization and asynchronization**: By analyzing the time - averaged detuning \(\langle \dot{\xi}\rangle\) and the time - averaged order parameter \(\langle p\rangle\), the researchers determine the boundary between synchronization and asynchronization. In particular, in the case of \(J_2 = 0\), the boundary curve of the asynchronous region is obtained analytically.

3. **Effect of inter - module coupling**: The researchers explore the effects of the inter - module coupling strength \(K_{12}\) and the intra - module coupling strength \(K_2\) on the system dynamics. The results show that as \(K_2\) increases, the asynchronous region gradually shrinks and eventually disappears.

4. **Simulation of actual systems**: To verify the results of the theoretical analysis, the researchers conduct direct simulations and use natural frequencies with a Gaussian distribution for numerical integration. The simulation results show that the behavior of the system is consistent with the theoretical analysis, including the existence of exciting non - linear oscillations and asynchronous islands.

In summary, this paper aims to gain an in - depth understanding of the dynamic behavior of matrix - coupled Kuramoto oscillators in modular networks, especially the synchronization and asynchronization characteristics of the system under different parameters, and the influence of inter - module coupling on these characteristics.