Abstract

The synchronization behaviors of coupled oscillators under time-varying couplings are of both theoretical and practical significance. While recent studies show that synchronization is suppressed by time-varying coupling in general, the underlying mechanism is still not very clear. Here, by the kernel of sinusoidal coupling function, we revisit the effects of periodic coupling on the synchronization of networked phase oscillators. It is found that the suppressed synchronization by periodic coupling is attributed to the formation of synchronization clusters in the transition from desynchronization to global synchronization. The clusters are different in size and frequency but are all locked to the frequency of the periodic coupling. We demonstrate this phenomenon numerically in different network models and conduct a theoretical analysis on the numerical results based on the method of dimension reduction. The findings extend our knowledge on the dynamical responses of a complex network to external drivings, and shed lights on the mechanism of suppressed synchronization in periodically coupled oscillators.

1. Introduction

The study of collective behaviors covers many different fields, which range from physics to chemistry, biology, and social economy [14]. However, realistic systems are unavoidably affected by the surrounding circumstances. Typical examples include the seasonal variations in ecosystems [5], the circadian rhythm of the nervous system [6], the diurnal cycles on gene expressions [7], and the electrical and magnetic stimulations in brain areas [8]. Pharmacological and neurological evidence has revealed that the molecular and cellular pathology of neural disorders can be altered by externally added low-frequency periodic stimuli, yet the underlying mechanism remains unknown [810]. Revealing the mechanism behind these collective behaviors has always been an important research field of nonlinear dynamics and network science which expands the understanding of the macroscopic dynamics in complex systems [3, 4].

Synchronization is a typical collective behavior in network science. The Kuramoto model is a classical model to describe the synchronization phenomenon of coupled phase oscillators [1, 2]. When the coupling strength between oscillators exceeds a certain critical value (i.e., critical coupling strength ), the partial synchronization of coupled phase oscillators can be observed. As the coupling strength further increases, more and more oscillators are added to the synchronized cluster, forming a macroscopic nonzero order parameter. Since it was introduced, the classical Kuramoto model and its generalization have received wide attention aimed at analyzing collective synchronization [3, 4].

The robustness of synchronization to external perturbations is one of the central issues in probing oscillator synchronization. Due to the nonlinear characteristics of system dynamics, many interesting phenomena can be produced by the existence of external perturbations, e.g., the controlling synchronization by periodic signals [1113], the steering synchronization by periodic drivings [1418], the generating synchronization by specifying initial conditions [19], and the enhancing synchronization by random noise [2023]. Similar studies also exist in ecosystems (where the interaction of biological species are changing with seasons) and metabolic systems (where gene expressions are affected by diurnal cycles) [5, 7]. Recently, the field of time-varying coupling has received a great deal of attention [2430]. It should be noted that periodic and even random modulation of the coupling strength does not always enhance synchronization [17, 31]. For example, periodic coupling suppresses synchronization in coupled phase oscillators in general, and the suppressed synchronization by periodic coupling has been attributed to the asymmetric oscillation of the temporal order parameter at the macroscopic scale [31]. However, it is still not very clear how this happened at the microscopic scale. In coupled nonlinear systems, some microscopical synchronous behaviors are exhibited when the inner elements compete or cooperate [29, 32]. In these studies, two questions are naturally raised: How the synchronization behaviors of coupled oscillators are affected by periodic couplings at the microscopic scale? How the critical coupling strength of coupled oscillators are affected by periodic couplings?

Motivated by the questions mentioned previously, we revisit how the synchronization behaviors of networked phase oscillators could be affected by periodic coupling [31]. Using the dimension reduction method, we analyze the impacts of periodic coupling on synchronization and give the formula of synchronization order parameter analytically. In particular, it is proven theoretically that the critical coupling strength remains unchanged by the periodic coupling. Numerical results demonstrate that the adoption of periodic coupling will suppress synchronization in general, showing the decreased synchronization order parameter as compared with the results of constant coupling. The theoretical predictions of the dimension reduction method are in good agreement with the results of numerical simulations. We investigate further the microscopic mechanism underlying the observed phenomena using the bifurcation tree, which shows clearly that the suppressed synchronization by periodic coupling is attributed to the formation of microscopic clustering patterns in the system. The findings shed lights on the interaction between complex nonlinear systems and external stimuli, and it might have implications for the interventions (evaluations) of function (operation) in some realistic systems, for example, the therapeutic interventions of neurological disorders, the control of cellular metabolism by diurnal cycles and circadian rhythm, and the evaluation of the impacts of seasonal variation on ecological processes [68].

The rest of this paper is organized as follows. In the next section, we will present the generalized model with a unimodal distribution of natural frequencies. In Section 3, we will derive a dimension reduction on the effects of periodic coupling on synchronization. Specially, we will provide rigorous proof of the critical strength coupling. In Section 4, we explore the microscopic mechanism of the suppressed synchronization, and discuss the generalizations to different network models, containing Zachary’s karate club network and the neural network of nematode C. elegans. The conclusion will be given in Section 5.

2. The Generalized Model

We consider the dynamics of network-coupled phase oscillators governed by the following equations:withwhere are the oscillator (node) indices, is the phase of the th oscillator, and and are, respectively, natural frequency of the th oscillator and the size of the system. The natural frequency follows a prescribed probability density function . The coupling relationship between oscillators can be described by the adjacency matrix , with if two oscillators are connected by a link, otherwise . represents the number of connections associated with oscillator , that is, its degree. In contrast to the case of the classical Kuramoto model where the coupling strength among oscillators is an invariant constant, the generalized model considered here involves a coupling strength function that depends on the coupling frequency and coupling amplitude of periodic coupling [15]. If , we have and the model is equivalent to the generalized Kuramoto model with constant coupling [4]. Furthermore, when all the oscillators are connected (i.e., for all the nondiagonal elements in matrix ), equation (1) depicts the classical Kuramoto model [4]. The main work of this paper is to investigate how the variations of and will affect the synchronization degree of the oscillators at the microscopic scale.

Following the references [33, 34], we select the natural frequency from the Lorentzian distribution as follows:

Here, and are the scale parameter and the central frequency of the Lorentzian distribution, respectively. To describe the synchronization behaviors of the oscillator, Kuramoto introduces the macroscopic order parameter, which is defined as follows [2]:where denotes the modulus of the complex order parameter , is a collective phase, and is an imaginary unit. describes the degree of synchronization between oscillators, such that denotes the desynchronization state, whereas is the global (perfectly) synchronization state.

3. Theoretical Analysis

To make the model theoretically tractable, we choose and employ the all-to-all coupling structure (all-coupled model). Using the dimension reduction method proposed by Ott and Antonsen [33, 34], we investigate the dynamical mechanism of the network-coupled phase oscillators. Under the thermodynamic limit condition , the state of the system at time can be depicted by a continuous probability function . The evolution of satisfies the continuity equation as follows:with

The complex order parameter is written in the continuous form as follows:

Here, denotes the complex conjugate of . Expanding the Fourier series of with respect to , we can obtain the following:

Here, . stands for the complex conjugate of . Using the Ott-Antonsen ansatz method [33, 34], we getwe can obtain the following equations by substituting equation (8) into equations (5) and (7):

The Lorentzian distribution can be described as . We thus obtain by the residue theorem. Substituting this result into equation (10), we obtain the following nonlinear equation with respect to the order parameter:where is the modulus of .

For the case of coupling frequency , the solution of equation (12) is as follows [33, 34]:where and is the initial value of the order parameter. Equation (13) shows that the value of approaches for and approaches 0 for .

For the case of coupling frequency , the desynchronization state is described by . The linearization of is performed in equation (12). Letting , we obtain

Neglecting the higher order term , we have

Therefore, we get the averaged eigenvalue as follows:

Here, the brackets denote a time average. Hence, will be transversely stable if . It is found that the frequency and amplitude of periodic coupling do not appear on the right-hand side of equation (16). This decoupling process shows that is the eigenvalue connected with the transverse perturbation (In the case where is the incoherent state). If , we can deduce that . Therefore, the critical coupling strength (= ) is independent of periodic coupling, and it depends only on the scale parameter of the Lorentzian distribution.

When , equation (12) cannot be solved analytically. Therefore, we use numerical simulation to solve it. By solving equation (12) numerically, the variation of in the parameter space is plotted in Figure 1. It shows that compared to the constant coupling, the synchronization performance of the system deteriorates with periodic coupling, and the order parameter is gradually decreased with decreasing frequency of periodic coupling. Furthermore, equation (12) indicates that the evolution of is independent of the central frequency, . These studies are in line with the findings in the periodically coupled phase oscillators [31]. Different from the previous study [31], we provide theoretical analysis for the critical strength coupling that is independent of periodic coupling.

4. Numerical Results

In this section, we first consider the all-coupled case and show that suppressed synchronization emerges using periodic coupling. We verify our result using the theoretical analysis. Then, we explore the underlying microscopic mechanism based on the method of bifurcation tree, and finally, we extend the study to complex network cases.

4.1. All-Coupled Case

In numerical simulations, we choose the system size as , and update equations (1) and (2) using the fourth-order Runge–Kutta algorithm with step size and the transient time . In obtaining the order parameter, , we average the instant order parameter, , over a period of . Figure 2 (a1–c1) show the values of versus the coupling amplitude in the all-coupled model. For all the cases, the values of are staying around 0 when . When , the value of for the case of is clearly smaller than that of constant coupling , indicating that synchronization is suppressed. Furthermore, for the case of , we find that the value of is very close to that of constant coupling, signifying that synchronization is hardly affected. One can notice that the synchronization is suppressed for periodic coupling, and the theoretical results in equation (12) are in line with the simulation results in equation (1).

While recent studies show that synchronization is suppressed by time-varying coupling in general, the underlying mechanism is still not very clear [31]. Here, by the method of bifurcation tree, we illustrate the effects of periodic coupling on synchronization. To characterize the underlying mechanism, the microscopic effective frequencies are defined as follows [29, 32, 35]:with being the average time window. If two oscillators and are phase locked, then , otherwise . When all the oscillators have a common value, it implies that a global phase-locking of the system has occurred. However, if behaves like a random form, it indicates that the system is in an incoherent state.

For the effective frequencies , we fix and choose . Figure 2(a2–c2) plot the bifurcation trees of the effective frequencies versus the coupling amplitude for different coupling frequencies . It initially forms a cluster centered at the central frequency for all three cases. For coupling frequency , as coupling amplitude increases from 0, a majority of oscillators move close to , and it finally generates a large cluster (Figure 2(a2)). For coupling frequency , as coupling amplitude increases from 0, a majority of oscillators form several clusters which are centered at (Figure 2(b2)). This decreases the degree of synchronization and indicates that coupling frequency suppresses synchronization in coupled phase oscillators. As coupling amplitude continually increases, more oscillators merge into the synchronous cluster . This phenomenon is different from the case of . For coupling frequency , as coupling amplitude increases from 0, a majority of oscillators form several clusters which are centered at (Figure 2(c2)). This decreases the degree of synchronization and indicates that coupling frequency suppresses synchronization in coupled phase oscillators. However, because are away from the central frequency , the degree of synchronization has hardly affected by the periodic coupling. It indicates that coupling frequency has little effect on the synchronization performance of coupled phase oscillators. As coupling amplitude continually increases, more oscillators merge into the synchronous cluster . In contrast to the case of constant coupling , a majority of oscillators form several clusters which are centered at for . The clusters are different in size and frequency, but are all locked to the frequency of the periodic coupling (with a fairly uniform spacing between adjacent clusters). Furthermore, the closer is the central frequency , the bigger is the cluster size (Figure 2(a2–c2)).

Therefore, it is found that when coupling amplitude is small , although the variation of , the effective frequencies are almost a random form for all the cases. As increases from , the effective frequencies have clearly several common values for the case of periodic coupling, indicating that synchronization is suppressed by periodic coupling. However, for the case of high-frequency coupling , we observe that the effective frequencies are approximate to that of constant coupling for the whole area, signifying that synchronization is less influenced by high-frequency coupling. These studies demonstrate that the adoption of periodic coupling will suppress synchronization in general, and the degree of synchronization gradually deteriorates with decreasing the frequency of periodic coupling. The findings of this study are reported in the periodically coupled phase oscillators [31]. It is important to note that the suppressed synchronization is attributed to the formation of microscopic clustering patterns in our study, which was previously not found in the periodically coupled phase oscillators [31].

When the time scale of the nodal dynamics is almost identical to that of periodic coupling, the synchronization performance could be greatly changed in the periodically coupled chaotic oscillators [15]. In particular, it is shown that at some characteristic frequencies of the periodic coupling, the upper bound of stable synchronization in the parameter space of the coupling strength can be enhanced. For the periodically coupled phase oscillators, the averaged time scale of the nodal dynamics is measured by the central frequency in the Lorentzian distribution. However, it is shown that is irrelevant to the synchronization performance for periodic coupling [31]. This arouses our interest in studying the influence of on synchronization using bifurcation trees.

Changing to and 2, we plot the variations of with respect to for in Figure 3. We see that the clusters are different in size and frequency, but are all locked to . Furthermore, the closer is the central frequency , the bigger is the cluster size (Figures 2 and 3). This signifies that the clusters are dependent on the central frequency and the frequency of the periodic coupling. We also see that the results of and 2 are identical to that of except for the position of clusters. It suggests that the synchronization performance of periodically coupled phase oscillators is irrelevant to the time scale of the nodal dynamics (the central frequency of the Lorentzian distribution). The results of this study are in line with the findings in the periodically coupled phase oscillators [31]. It should be noted that the main phenomenon is independent of the specific form of the frequency distribution. It has been verified from bifurcation trees that Gaussian distribution and the triangular frequency distribution also have similar phenomena.

4.2. Complex Network Cases

We finally test the universality of the observed phenomena in the system with complex network structures. The first system we investigate is the Zachary’s karate club network which consists of 34 nodes and 78 links [36]. It is one of the most commonly used in complex networks, sociological analysis, and other fields. Significantly, it has two clusters in Zachary’s karate club network. Besides the network structure, other settings of Zachary’s karate club network are the same as the globally connected network in Figure 2, such as the coupling forms and the frequency distribution. Figure 4 shows the variations of with respect to for the three cases of coupling frequency , coupling frequency and coupling frequency . We see that before the onset of synchronization, the effective frequencies of for all three cases are very similar; after the onset point, the effective frequencies of for the cases of are different from that of constant coupling . Specifically, we see that several clusters have been formed for . These results are in agreement with the results obtained in the globally connected network in Figure 2. It should be noted that due to the sparse connectivity of Zachary’s karate club network, the synchronization order parameter is not precisely described by equation (12). However, we have revealed that the dynamic mechanism is universal, and is independent of the network structure (Figures 2 and 4). Therefore, we see that compared to the case of constant coupling, synchronization is suppressed by periodic couplings in general, and this is attributed to the formation of clusters. In addition, the spacing between adjacent clusters is identical to the frequency of periodic coupling.

The second system we study is the neural network of the nematode C. elegans, which contains 297 nodes and 2,148 links [37]. The natural frequencies of the oscillators are still selected from the Lorentzian distribution, and the parameters are the same as those of the globally connected network (Figure 2). Figure 5 shows the variations of with respect to for the cases of coupling frequency , coupling frequency , and coupling frequency . The results are similar to the results of Zachary’s karate club network (Figure 4), i.e., synchronization is suppressed by periodic couplings compared to the scenario of constant coupling. Therefore, periodic coupling suppresses synchronization in the general network structure and this is attributed to the formation of clusters.

5. Conclusion

As the dynamical basis of the normal brain function, the synchronous activity of distributed neurons has attracted the interest of neuroscientists [38]. In particular, there have been a variety of stimulation techniques aimed at improving the brain function or restoring the brain function from disorders in the past few decades [8, 9], including deep-brain stimulations and transcranial direct-current stimulations. Although this evidence has revealed that the molecular and cellular pathology of neural disorders can be changed by externally added low-frequency periodic stimuli. However, the underlying mechanism remains unclear.

The collective behavior of a complex system implies the emergence of a microscopic order from the organizations of populations of units with mutual interaction. Usually, the number of degrees of freedom at the macroscopic scale is so large that an exact description of the system at this level is impossible and also unnecessary. Therefore, a microscopic study of the complex system is significant [29, 32]. Here, treating the stimuli as periodic perturbations that modify the strength of couplings, we evaluate how the synchronization behaviors of networked phase oscillators could be affected by periodic coupling at the microscopic scale.

In summary, we have studied the generalized model with periodic coupling in which the coupling strength is changing periodically with time. We theoretically and numerically verify that the critical coupling strength is irrelevant to the periodic coupling. The bifurcation trees have revealed that the system can exhibit rich clustering features in a system of finite size. The synchronization and desynchronization are characterized by the bifurcation tree, and the suppressed synchronization is attributed to the formation of clusters. In addition, these analyses can be generalized to realistic networks. The findings extend our knowledge on the dynamical responses of a complex network to external perturbations, and provide lights on the mechanism of suppressed synchronization in periodically coupled oscillators. Future studies could apply this model to study these periodically coupled biological, physical, social, and technological systems and to better understand their associated behaviors.

Data Availability

The data used to support the findings of this study are available upon request from the corresponding author.

Conflicts of Interest

The author declares that he has no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under the Grant no. 12102004, Natural Science Foundation of Ningxia under the Grant no. 2020AAC03240, Research Foundation of North Minzu University under Grant nos. 2020XYZSX06 and 2020KYQD05, and Construction Project of First-Class Disciplines in Ningxia Higher Education under the Grant no. NXYLXK2017B09.