Abstract:We present an adaptive coupling strategy to induce hysteresis/explosive synchronization (ES) in complex networks of phase oscillators Sakaguchi-Kuramoto model). The coupling strategy ensures explosive synchronization with significant explosive width enhancement. Results show the robustness of the strategy and the strategy can diminish (by inducing enhanced hysteresis loop) the contrarian impact of phase frustration in the network, irrespective of network structure or frequency distributions. Additionally, we design a set of frequency for the oscillators which eventually ensure complete in-phase synchronization behavior among these oscillators (with enhanced explosive width) in the case of adaptive-coupling scheme. Based on a mean-field analysis, we develop a semi-analytical formalism, which can accurately predict the backward transition of synchronization order parameter.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to achieve explosive synchronization (ES) in the phase oscillator model (Sakaguchi - Kuramoto model) in complex networks and significantly enhance its explosive width. Specifically, the researchers proposed an adaptive coupling strategy. By adjusting the parameter \(z\) in the coupling function, the explosive synchronization phenomenon is induced and enhanced. The paper explores the effectiveness of this strategy under different frequency distributions and network structures, and analyzes its suppression effect on the influence of phase frustration.
### Main problems
1. **Induction and enhancement of explosive synchronization**: How to induce explosive synchronization in complex networks through an adaptive coupling strategy and significantly enhance its explosive width?
2. **Influence of phase frustration**: How to suppress the negative influence of phase frustration on network synchronization behavior through the proposed adaptive coupling strategy?
3. **Theoretical and numerical verification**: How to verify the effectiveness of the adaptive coupling strategy through semi - analytical methods and numerical simulations?
### Solutions
1. **Adaptive coupling strategy**: The researchers proposed an adaptive coupling function \(l_i = R^{z - 1}\), where \(R\) is the global synchronization order parameter and \(z\) is a positive real number. By adjusting the value of \(z\), the width of explosive synchronization can be significantly enhanced.
2. **Numerical simulation**: Through numerical simulations on different types of complex networks (such as scale - free networks, small - world networks), the effectiveness of the adaptive coupling strategy is verified.
3. **Semi - analytical method**: Based on the Ott - Antonsen ansatz, semi - analytical equations describing the system behavior are derived, and the backward transition points of the synchronization order parameter are predicted through these equations and compared with the numerical simulation results.
### Key findings
- **Enhancement of explosive width**: By increasing the value of \(z\), the width of explosive synchronization can be significantly enhanced, regardless of the network structure or frequency distribution.
- **Suppression of phase frustration**: The adaptive coupling strategy can effectively suppress the negative influence of phase frustration on network synchronization behavior.
- **Consistency between theory and experiment**: The results obtained by the semi - analytical method are highly consistent with the numerical simulation results, verifying the accuracy of the theoretical model.
### Formula summary
- **Dynamics equation of phase oscillator**:
\[
\frac{d\theta_i}{dt}=\omega_i+\lambda l_i\sum_{j = 1}^N A_{ij}\sin(\theta_j-\theta_i-\alpha)
\]
where \(\omega_i\) is the natural frequency of the \(i\)-th oscillator, \(A_{ij}\) is an element of the adjacency matrix, \(\alpha\) is the phase frustration parameter, \(\lambda\) is the coupling strength, and \(l_i = R^{z - 1}\) is the adaptive coupling term.
- **Global synchronization order parameter**:
\[
Re^{i(\psi-\alpha)}=\frac{1}{N}\sum_{j = 1}^N e^{i(\theta_j-\alpha)}
\]
where \(R\) represents the synchronization degree and \(\psi\) is the average phase.
- **Semi - analytical equation**:
\[
\frac{\partial\beta}{\partial t}+i\omega\beta+\frac{\lambda k}{2}R^{z - 1}(C\beta^2 - C^*) = 0
\]
where \(\beta\) is the Fourier coefficient of the oscillator density function and \(C\) is the global synchronization order parameter in complex form.
Through these methods and findings, the paper successfully solves the problem of inducing and enhancing explosive synchronization in complex networks and provides theoretical and experimental support.