Experimental study of the irrational phase synchronization of coupled nonidentical mechanical metronomes.

It has recently been observed in numerical simulations that the phases of two coupled nonlinear oscillators can become locked into an irrational ratio, exhibiting the phenomenon of irrational phase synchronization (IPS) [Phys. Rev. E 69, 056228 (2004)]. Here, using two coupled nonidentical periodic mechanical metronomes, we revisit this interesting phenomenon through experimental studies. It is demonstrated that under suitable couplings, the phases of the metronomes indeed can become locked into irrational ratios. Numerical simulations confirm the experimental observations and also reveal that in the IPS state, the system dynamics are chaotic. Our studies provide a solid step toward further studies of IPS.

Introduction

As a universal concept in nonlinear science, synchronization has been widely observed in both natural and manmade systems [1]. One well-known example is the synchronization of hand clapping after a show performance, in which audiences are found to begin clapping their hands with the same rhythm after a short transient period [2]. Additional examples of synchronization are observed in biological systems [3,4], psychiatric diseases [5,6], and many other contexts [7,8,9]. In general, synchronization refers to the coherent motion of coupled oscillators, which typically occurs when the coupling strength exceeds some critical value. Depending on the specific form of the coherent motion, synchronization can be classified into various types, including complete synchronization, phase synchronization [10,11,12], generalized synchronization, and lag synchronization [1].

Phase synchronization is defined as the locking of the phases of coupled oscillators [13,14]. Specifically, let ϕ 1,2(t) be the instantaneous phases of the coupled oscillators; phase synchronization is said to be achieved if |ϕ 1 (t) – (m / n) ϕ 2(t)| < c, where c ∊ [0,π] is a constant and m and n are two integers. Because any 2π increment in ϕ 1 is accompanied by an (m / n)2π increment in ϕ 2, this type of synchronization is also called m:n phase synchronization. Phase synchronization is trivial when the oscillators are periodic because | ϕ 1(t) – (m / n)| ϕ 2(t) < c even in the absence of coupling, given that their natural frequencies (the frequencies of the isolated oscillators) f 10,20 satisfy f 10 / f 20 = m / n. As such, phase synchronization is primarily of interest for chaotic systems [15,16]. Because only the phases of the oscillators are correlated, phase synchronization represents a weaker type of coherent motion than complete synchronization and is typically observed at a lower coupling strength than the latter. However, unlike complete synchronization, which requires the coupled oscillators to exhibit identical dynamics, phase synchronization is observable for nonidentical oscillators, making it more relevant to practical situations.

Previous studies of phase synchronization have been primarily focused on rational phase locking, i.e., cases in which the ratio r = m:n is a rational number. Recently, it has been shown numerically that under certain circumstances, the phases of two coupled chaotic oscillators can also become locked into an irrational ratio, i.e., |Δ ϕ − c| < π, where Δ ϕ (t) = | ϕ 1 − r ϕ 2| is the phase difference and r is an irrational ratio. In Ref. [16], this phenomenon is called irrational phase synchronization (IPS). The numerical discovery of IPS represents a significant extension to the current knowledge of phase synchronization, as it is common sense that for coupled periodic oscillators, the resonance tongues that correspond to the irrational frequency ratios are of zero measure. The question that naturally arises, therefore, is the following: can IPS be observed in a realistic system?

Motivated by the numerical discovery of IPS in Ref. [16], in the present work, we study IPS through experiments. Specifically, by coupling two nonidentical periodic metronomes, we demonstrate that IPS is, indeed, observable in realistic systems. Moreover, through numerical simulations of the corresponding model, it is found that in the IPS state, the system dynamics are weakly chaotic, i.e., the largest Lyapunov exponent is positive.

Materials and Experimental Results

The experiment setup employed in our studies is depicted in Fig. 1. Two metronomes are placed on a platform supported by two aluminum pipes. The specific metronome used is the Taktell Piccolino (Series 890) manufactured by Wittner GmbH & Co. KG, Germany, which has a mass of 94 g and whose energy is supplied by a manual winding spring. By adjusting the position of the bob on the pendulum, the frequency of the metronome can be changed continuously from 40 (largo) to 208 (prestissimo) beats per minute (BPM). The platform underneath the metronomes is formed of two sheets of A4 paper folded in a zigzag pattern. The mass of the platform is approximately 4.366 g. Two identical aluminum pipes with inner and external diameters of 39 mm and 41 mm, respectively, are placed in parallel below the platform. Because the motion of the pendulum bob is perpendicular to the pipe axes, a bidirectional coupling between the metronomes is established through the platform. Because the platform is light and the diameter of the pipes is large, the coupling between the metronomes is relatively strong. Underneath the pipes is placed a horizontal support that serves as the base of the experiment.

Fig 1

The experimental setup.

(a) The experimental equipment and the CCD acquisition system. (b) The details of the experimental equipment.

The motions of the metronomes are monitored by a charge-coupled device (CCD) acquisition system, and the recorded data are analyzed using computer software [17,18]. To facilitate the recording process, we affix a red wafer to the top of the pendulum on each metronome. This slightly increases the mass of the metronome (by approximately 0.01 g) and, consequently, slightly decreases the frequency of the metronome. The metronome device, the acquisition system, and the data analysis software have been developed and improved through a series of studies and have proven to be stable and reliable in various experiments [17,18,19]. The metronome device was initially designed by Zou et al. in Ref. [19] and was later improved on by several other researchers [20,21]. In Ref. [20], it is shown that by varying the coupling strength and the frequency mismatch, the metronomes can achieve phase or envelope synchronization. In Ref. [21], it is further demonstrated that by varying the system parameters, the basin of the phase synchronization can be significantly modified.

In our experiments, we fix the frequency of the first metronome to 160 BPM while varying the frequency of the second metronome gradually from 152 to 168 BPM. Note that the attachment of the wafer shifts the natural frequency of the first metronome to f 10 = 158.47 BPM. With the natural frequency of the second metronome chosen to be f 20 = 156.52 BPM, we plot in Fig. 2A the time evolution of the two phases, ϕ 1(t) and ϕ 2(t). The envelopes of ϕ 1 and ϕ 2 exhibit the phenomenon of anti-phase synchronization [20]. To provide further details on the dynamics of the phases, in Fig. 2B, we present the system trajectory in the spaces. It is observed that the trajectory is irregular in both spaces. Using the Poincaré surface-of-section method (i.e., by recording ϕ 2 at the moments when and ), we replot in Fig. 2C the discrete evolution of the system trajectory in the space, which reveals that the system dynamics are chaotic. A careful evaluation of the time evolution of ϕ 1,2 also indicates that |Δ ϕ − c| < π, with r = π / 3.1 (the details of how to determine this irrational ratio are presented in the next section). That is, the two metronomes achieve the IPS state. To validate the experimental findings, we gradually vary the natural frequency of the second metronome from 152 to 168 BPM and investigate the synchronization behaviors of the two phases. The experimental results are presented in Fig. 2D. It is observed that as f 20 increases, the phase-locking ratio, r, varies continuously. In particular, near the testing frequency f 20 = 156.52, r continuously decreases as f 20 increases. According to Ref. [16], given that r is continuously varying over some interval of the parameter space, IPS can always be achieved through careful tuning of the parameter within this interval.

Fig 2

The experimental results.

For f 10 = 158.47 BPM and f 20 = 156.52 BPM, we present (a) the time evolutions of ϕ 1,2, where envelopes of ϕ 1,2 exhibit anti-phase synchronization; (b) the phase trajectory in the space; and (c) the Poincaré maps. (d) The variation of the phase-locking ratio, r, as a function of the natural frequency of the second metronome, f 20. The arrow denotes the point at which r = π / 3.1.

Methods and Numerical Results

To confirm the experimental findings, we next study IPS using numerical simulations based on a simplified model. By treating the metronomes as self-sustained periodic oscillators and regarding the platform as being driven by a liner spring, the experimental setup depicted in Fig. 1 can be simplified to the model schematically illustrated in Fig. 3 [22,23]. Based on this simplified model, the phase dynamics of the metronomes, ϕ 1,2(t), and the displacement of the platform, x(t), can be calculated using the following Euler-Lagrange equations (please refer to Ref. [18] for further details regarding how to derive the equations from the simplified model): Where i = 1,2 is the metronome index. In Eq. (1) (the equation for the phase dynamics), m is the mass of the pendulum bob, l is the pendulum length, c is the damping coefficient, and g ≈ 9.8 is the gravitational acceleration. M is the external energy supplied by the metronome escapement, which is used to compensate for the energy dissipated by friction. The external energy has the form where θ is a small angle (π / 36 rad) below which the escapement is activated. In Eq. (2) (the equation for the platform displacement), M is the total mass of the two metronomes and the paper platform, c is the friction of the pipes, and k is the linear damping of the beam.

Fig 3

A schematic illustration of the simplified model.

For consistency with the experimental setup, we set the (dimensionless) parameters in the modeling equations as follows: m 1 = m 2 = m = 1, c = 0.01, M = 30, θ = π / 36, c = 1.5, and k = 1. By setting l 1 = 0.014248, we specify f 10 = 160 BPM for the first metronome. The natural frequency of the second metronome can be modified continuously by adjusting l 2, following the relation l 2 ≈ g(60/πf 20)2. For the simulations, we set the initial conditions to and and choose ϕ 1,2(t = 0) randomly within the range (−π / 2,π / 2). Eqs. (1) and (2) are integrated using the 4th-order Runge-Kutta method with a time step of Δt = 0.0001. In our simulations, the working frequencies, f 1,2, of the metronomes are obtained by averaging over a period of t = 1 × 104. (Please note the difference between the working frequencies, f 1,2, and the natural frequencies, f 10,20, of the metronomes. The former characterizes the coupled metronomes, whereas the latter describes the isolated metronomes.) The first question we are interested in is the following: can the working frequencies of the metronomes become locked into an irrational ratio, i.e., r′ = f 1 / f 2, where r′ is an irrational number? (The locking of the working frequencies into an irrational ratio is a necessary condition for the observation of IPS.)

To construct a global picture of the locking of the working frequencies, we first consider the variation of r′ as a function of f 20. The results are presented in Fig. 4A. It is observed that with the exception of the interval [f 10 − 2.4,f 10 + 2.4], where r′ = f 1 / f 2 = 1, the value of r′ is continuously decreasing as the increase of f 20. This behavior of r′ is consistent with the experimental results presented in Fig. 2D. To acquire further details on the locking of the two working frequencies, we then investigate the variations of f 1,2 as functions of f 20. The results are plotted in Fig. 4B. It is observed that as f 20 increases, both f 1 and f 2 vary continuously. This finding confirms the continuous variation in r′ observed in Fig. 3A. Interestingly, in Fig. 4B, it is also observed that within the interval [f 10 − 2.4,f 10 + 2.4], the locked frequencies (with the locking ratio r = 1) are largely divergent from the natural frequency of the first metronome, f 10. Specifically, f 1,2 > f 10,20 holds in the interval [f 10 − 0.4,f 10 + 0.4], and f 1,2 < f 10,20 holds in the intervals [f 10 − 2.4,f 10 − 0.4] and [f 10 + 0.4,f 10 + 2.4]. This frequency fluctuation behavior may be induced by the resonance of the metronome oscillation with the platform motion, and this possibility will be investigated in our future studies. To exclude the possibility of quasiperiodic motion, we also calculate the variations of the Lyapunov exponents, λ 1,2,3, as functions of f 20, based on the algorithm introduced in Refs. [24,25]. The results are plotted in Fig. 4C. It is observed that the largest Lyapunov exponent, λ 1, is equal to 0 within the region [f 10 − 2.4,f 10 + 2.4] and is positive in the other regions. Fig. 4 thus implies that in the regions f 20 < f 10 − 2.4 and f 20 > f 10 + 2.4, the working frequencies of the metronomes can become locked into irrational ratios and the system dynamics are chaotic.

Fig 4

The numerical results.

(a) The variation in the ratio of the working frequencies, r′ = f 1 / f 2, as a function of f 20 for f 10 = 160 BPM. As a reference, the ratio of the natural frequencies, r 0 = f 10 / f 20, is also plotted. (b) The variations in the working frequencies, f 1,2, as functions of f 20. Resonance-like behavior is observed in the region f 20 ∊ [160 − 2.4,160 + 2.4]. (c) The variations in the Lyapunov exponents, λ 1,2,3, as functions of f 20. At f 20 = 157.08915, the working frequencies are locked into the ratio r′ = π / 3.1, and the largest Lyapunov exponent is λ 1 ≈ 0.001.

Setting f 20 = 157.08915, we now investigate the evolutional dynamics of the phases in detail. From Fig. 4A and B, we know that in this case, the phases will be locked into the ratio r ≈ π / 3.1, and from Fig. 4, we also know that the system dynamics are chaotic, with λ ≈ 0.001. In Fig. 5A, we plot the time evolution of the two phases ϕ 1,2. Consistent with the experimental results [Fig. 2A], the envelopes of ϕ 1 and ϕ 2 exhibit the phenomenon of anti-phase synchronization. Applying the same method used in the experimental analysis, in Fig. 5B, 5C and 5D, we present the phase trajectory in the space, the Poincaré maps and a magnified view of the latter. Consistent with the observation of the positive Lyapunov exponent, the phases exhibit chaotic motion [26].

Fig 5

Examples of IPS observed in numerical simulations.

For f 10 = 160 BPM and f 20 = 157.08915, we present (a) the time evolution of the phases, ϕ 1,2, where anti-phase synchronization is observed in the phase envelopes; (b) the phase trajectory in the space; (c) the Poincaré maps; and (d) a magnified view of the Poincaré maps. (e) The variation of the phase difference Δ ϕ at t = 1 × 104 as a function of the phase ratio r. Δ ϕ reaches a minimum as r approaches r 1 = π / 3.1. (f) The time evolution of Δ ϕ for r = r 1 (red) and r = r 1 − 4.17 × 10−4 (blue). (g) The variation of the phase difference Δ ϕ at t = 1 × 104 as a function of r near . (f) The time evolution of Δ ϕ for r = r 2 (red) and r = r 1 − 7.93 × 10−4 (blue).

Can the phases also become locked into an irrational ratio? Because of the chaotic nature of the system dynamics, the locking of the working frequencies does not guarantee the same locking of the phases. Furthermore, because of the limited computational precision, it is difficult to judge whether the phases could become locked exactly into a truly irrational ratio or merely into a decimal ratio with a finite number of digits. To address this question, an effective approach is to investigate the variation of the phase difference Δ ϕ (t) = |ϕ − rϕ | as a function of r in the neighboring region of the irrational ratio of interest. If the phase difference Δ ϕ gradually approaches 0 as r approaches the irrational ratio r 1, then one can conclude that at r 1, the two metronomes achieve the state of IPS [16]. Using this method, we plot in Fig. 5E the variation in Δ ϕ (obtained at the instant t = 1×104 of the system evolution) as a function of r in the range r ∊ [1.1032,1.0140]. It is clearly observed that Δ ϕ approaches 0 as r approaches r = π / 3.1 ≈ 1.013417. As depicted in Fig. 5F, the value of Δ ϕ remains equal to 0 for r = r 1 but is monotonically increasing when r = r 1 − 4.17×10−4.

Fig. 5G and H illustrate another example of IPS. In Fig. 5G, we change the natural frequency of the second metronome to f 20 = 157.06244 and plot the variation of the phase difference Δ ϕ as a function of r within the same range used in Fig. 5E. As r approaches , Δ ϕ monotonically decreases to 0. Again, in Fig. 5H, it is observed that Δ ϕ remains equal to 0 when r = r 2 and monotonically increases for r = r 2 − 7.93×10−4. By calculating the largest Lyapunov exponent, it is shown that the system is chaotic for f 20 = 157.06244, as λ 1 > 0. We thus conclude from Fig. 5E-5H that IPS can indeed be achieved between the coupled metronomes, given that the natural frequencies of the metronomes are properly chosen, e.g., (f 10,f 20) = (160,157.08915) (r = π / 3.1) or (f 10,f 20) = (160,157.06244) .

Discussion

A few remarks regarding the present work are in order. First, although it is impossible for the phases of two coupled oscillators to be locked exactly into an irrational ratio in practice, our studies show that the impossibility of this phenomenon is attributable only to limited experimental and numerical precision. In other words, through improvements to the experimental and numerical precision, IPS can be gradually approached in realistic systems. This fact can be partially understood from the results presented in Fig. 5E and 5F, where the phase difference Δ ϕ gradually decreases to 0 as r approaches r . Second, we would like to note that in the present work, although the IPS state is chaotic, the motion of the isolated metronomes is, in fact, periodic. The underlying mechanism of these chaotic dynamics is still unknown to us. Therefore, it will be important to confirm whether IPS can be experimentally realized in coupled chaotic oscillators. Third, in numerical studies, we have observed resonance-like behavior in the locked working frequencies, i.e., r′ = 1 near f 20. The dynamical mechanism is also unknown to us. Finally, considering the fact that many realistic systems consist of many coupled oscillators, it will be intriguing to investigate whether IPS can be observed in spatially extended systems, e.g., three oscillators coupled in a ring structure. We hope to address these questions in our future studies.

To summarize, we have experimentally studied IPS in coupled nonidentical metronomes. We find that when the natural frequencies of the metronomes are properly chosen, the phases of the two metronomes can become locked into irrational ratios. The experimental findings are verified by numerical simulations based on a simplified model. Our studies suggest that IPS is not merely of theoretical interest but also of practical concern.