Amplitude death in globally coupled oscillators with time-scale diversity

Wei Zou, Meng Zhan, and Jürgen Kurths

Phys. Rev. E 98, 062209 – Published 12 December 2018

Abstract

We analyze amplitude death (AD) in large systems of globally coupled oscillators with randomly distributed time scales. We show that the distribution of characteristic eigenvalues of a large but finite system is well approximated by the continuous and discrete spectra of the system in the thermodynamic limit. The stability analysis from the continuous and discrete spectra of the infinite system provides a fairly accurate prediction for the onset of AD in the large finite system. We prove the argument by examining the stability of AD in a paradigmatic system of coupled Stuart-Landau limit cycles with mismatched time scales. The proposed technique is extended to systems of globally coupled nonlinear oscillators of a general form with time-scale diversity, which is verified in coupled chaotic Rössler oscillators. Our study provides analytical insight into the understanding of the emergence of AD in populations of globally coupled dynamical systems.

Received 10 September 2018
Revised 9 November 2018

DOI:https://doi.org/10.1103/PhysRevE.98.062209

Physics Subject Headings (PhySH)

Coupled oscillators

Nonlinear Dynamics

Authors & Affiliations

Wei Zou^1,*, Meng Zhan², and Jürgen Kurths^3,4,5

¹School of Mathematical Sciences, South China Normal University, Guangzhou 510631, People's Republic of China
²State Key Laboratory of Advanced Electromagnetic Engineering and Technology, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
³Potsdam Institute for Climate Impact Research, Telegraphenberg, Potsdam D-14415, Germany
⁴Institute of Physics, Humboldt University Berlin, Berlin D-12489, Germany
⁵Saratov State University, 4410012 Saratov, Russia

^*weizou83@gmail.com

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand

Issue

Vol. 98, Iss. 6 — December 2018

Reuse & Permissions

Access Options

Author publication services for translation and copyediting assistance advertisement

Images

Figure 1
Distribution of eigenvalues of the coupled system (2) with $N = 30$ , $K = 5$ , $w_{j} = 10$ , and $σ = 0.4$ . The time scales $τ_{j}$ 's are uniformly distributed on $[1 - σ, 1 + σ]$ . The black squares denote the eigenvalues numerically obtained from the characteristic equation (4), which can be nicely predicted by the discrete and continuous spectra in Eqs. (8) and (9) represented by the red open circle and the red line, respectively. All eigenvalues have a negative real part implying that AD is stable.
Reuse & Permissions
Figure 2
Stability regions for AD of the coupled system (1) with $N = 30$ in the parameter space of $(σ, K)$ for (a) $w = 10$ and (b) $w = 5$ . The black squares are numerical results for the AD boundary obtained from directly integrating the coupled system (1), where the lower and upper black lines are the theoretical predictions by Eqs. (10) and (11), respectively. The time scales $τ_{j}$ 's are sampled from a uniform distribution over $[1 - σ, 1 + σ]$ .
Reuse & Permissions
Figure 3
(a) Stability region for AD of the coupled system (1) with $N = 30$ in the parameter space of $(K, w)$ for a uniform distribution of $g (τ)$ with $σ = 0.4$ . The coupling interval for AD decreases monotonically as $w$ decreases, and completely disappears if $w < w_{min} = 4.95$ . The black squares and the black solid line, respectively, denote the numerical simulations and the theoretical predictions as in Fig. 2. (b) The dependence of $w_{min}$ on $σ$ . The black square points of $w_{min}$ are obtained as in (a). The red solid line is a power-law fit $w_{min} \propto σ^{γ}$ with $γ = - 0.76$ , which is clearly demonstrated in the upper-right inset with the log-log plot.
Reuse & Permissions
Figure 4
Effect of the distribution of $w_{j}$ on AD induced by time-scale diversity in the coupled system (1). The stability regions for AD are depicted in the parameter space of $(σ, K)$ with $w_{j}$ and $τ_{j}$ independently and uniformly distributed over $[1 - σ, 1 + σ]$ and $[w - α, w + α]$ , respectively. $N = 100$ and $w = 5$ are fixed. For $α > 0$ , a lot of different combinations for $τ_{j}$ and $w_{j}$ are possible, where the two extreme cases are (i) $τ_{j} = 1 + σ (2 j - N - 1) / (N - 1)$ and $w_{j} = w + α (2 j - N - 1) / (N - 1)$ , and (ii) $τ_{j} = 1 + σ (2 j - N - 1) / (N - 1)$ and $w_{(N - j + 1)} = w + α (2 j - N - 1) / (N - 1)$ . The AD region expands for $α = 0.5$ with the case (i) (enclosed by the black lines) and shrinks for $α = 0.5$ with the case (ii) (enclosed by the blue lines) compared with that for $α = 0$ (enclosed by the red lines), as the level of diversity for the effective time scales is enhanced in the case (i) and is weakened in the case (ii) by the distribution of $w_{j}$ .
Reuse & Permissions
Figure 5
(a) Bifurcation diagram by plotting the local maxima of the mean of the $x_{j} (j = 1, \dots, N)$ for the globally coupled chaotic Rössler oscillators with $N = 100$ and a uniform distribution of $g (τ)$ with $σ = 0.4$ . With increasing the value of $K$ , AD is stabilized within a certain coupling interval of $0.1 < K < 0.72$ . (b) Stability region for AD in the parameter space of $(σ, K)$ . The black squares are numerical results for the AD boundary from direct integrations of the $N = 100$ coupled chaotic Rössler oscillators as in (a), where the lower and upper solid lines are the theoretical predictions from Eqs. (15) and (16) with the chaotic Rössler oscillator model.
Reuse & Permissions
Figure 6
Comparison of the bifurcation boundaries of the AD region by the method of order-parameter expansion in [40] with our results for (a) coupled Stuart-Landau oscillators and (b) coupled chaotic Rössler oscillators. The time scales $τ_{j}$ are uniformly distributed over $[1 - σ, 1 + σ]$ . The prediction of AD from the method of order-parameter expansion in [40] is indicated by the red solid lines. Our theoretical prediction is marked by the black solid lines, and the black squares denote the numerical observation of AD, which are the same as in Figs. 2 and 5, respectively. Our theory gives rise to a better prediction of the boundaries of the AD region compared with the method of order-parameter reduction.
Reuse & Permissions

Physical Review E

covering statistical, nonlinear, biological, and soft matter physics