Saddle torus and mutual synchronization of periodic oscillators

1. 1 2 3 4 Saddle torus and mutual synchronization of periodic oscillators N. Janson 1, A. Balanov 2,V. Astakhov 3, P.V.E. McClintock 4

2. Phase and frequency synchronization: definitions 1 2 F1,2 - phases f 1,2 – basic frequencies Frequency synchronization: n, m - integers Phase synchronization:

3. What happens in the phase space? torus limit cycle Poincare section stroboscopic section Poincare section stroboscopic section synchronous nonsynchronous

4. Periodic oscillator f0,fr ff Forced Periodic Oscillators C Uni-directional coupling, or forcing Forcing parameters: C – amplitude ff - forcing frequency Simplest form of forcing!

5. Synchronization Regions of Forced Periodic Oscillators Weak coupling – phase locking Strong coupling – suppression of natural dynamics C Synchronization regions suppression phase locking ff:f0 1:4 1:2 1 2

6. 1st classical mechanism of synchronization frequency locking Van der Pol oscillator under harmonic forcing Evolution of spectra:

7. 1st classical mechanism of synchronization frequency locking Van der Pol oscillator under harmonic forcing Stroboscopic section (in a period of forcing):

8. 2nd classical mechanism of synchronization suppression of natural oscillations Van der Pol oscillator under harmonic forcing Evolution of spectra:

9. 2nd classical mechanism of synchronization suppression of natural oscillations Van der Pol oscillator under harmonic forcing Stroboscopic section (in a period of forcing):

10. Periodic oscillator 1 Periodic oscillator 2 f2 f1 Mutual coupling Mutually Coupled Periodic Oscillators General forms of coupling functions Most popular coupling function: diffusive

11. Knowledge vs Assumption. Question 1 What can be (and is often) assumed: The structure of the phase space in mutually coupled periodic oscillators is similar to the one in forced oscillators. Question 1: Is this true? What was known previously: Multistability inside synchronization region is common in mutually coupled periodic oscillators*. * Aronson 1986; Astakhov et al 1990;J. Rasmussen et al 1996;D. Postnov et al 1999; S. Park et al 1999;Vadivasova et al 2000; E. Mosekilde et al 2003.

12. Linear vs Non-linear Coupling. Question 2. Linear Diffusive Symmetric Coupling in one Equation Non-linearly diffusively coupled Van der Pol oscillators Linearly diffusively coupled Van der Pol oscillators Most earlier studies are devoted to linearly coupled systems. Question 2: Does nonlinear coupling produce substantially different effect to the one by linear coupling?

13. Note When two periodic oscillators are coupled mutually, the smallest dimension of phase space is 4. Visualization of phase space structure is more difficult than with forced systems. One can reasonably well visualize only Poincare map.

14. System under Study Possibly, the simplest case: Two diffusively coupled Van der Pol oscillators Nonlinearities of partial oscillators Close to partial basic frequencies Detuning between oscilators I.e. around 1:1 synchronization region.

15. Established Facts about the System (1) • The dimension of phase space is 4. • 2. There is 1 fixed point at the origin for all parameter values. • 3. There are 6 cycles involved for parameters at least around • synchronization region 1:1. Namely, • 2 stable cycles S1 and S2 • 2 saddle cycles (with one-dimensional unstable manifold) • S1* and S2* • 2 twice saddle cycles (with two-dimensional unstable • manifolds) U1 and U2

16. Bifurcation Diagramfor Linear Coupling (C2=0) Solid: Saddle-Node; Dashed: Torus Birth/Death suppression locking

17. Established Facts about the System (2) 4. In locking region a pair of stable and a pair of saddle cycles lie on the same stable resonant torus. 5. Any saddle cycle can merge with any stable cycle (because they lie on the same torus). 6. Any saddle cycle can merge with any twice saddle cycle. 7. In suppression region there are no saddle or twice saddle cycles.

18. Question 3 In locking region – 1. Two stable cycles coexist. 2. Their basins of attractions are supposedly separated by manifolds of saddle cycles (like in case of a forced oscillator) In suppression region - 1. Two stable cycles coexist. 2. There are no saddle or twice saddle cycles (they have died via SN bifurcation). 3. Therefore, there are no separating manifolds of saddle cycles. Question 3: What separates the basins of attraction of two cycles in suppression region?

19. BDs for Nonlinear Coupling (different C1) Solid: Saddle-Node: Dashed: Torus Birth/Death suppression suppression locking locking suppression suppression locking locking

20. 3-Parameter Bifurcation Diagram C2 C1 p

21. p-C2 Bifurcation Diagram for C1 =0.1 Black solid: Saddle-Node: Dashed: Torus Birth/Death Green solid: line where stable torus vanishes A stable cycle + A stable ergodic torus 2 stable cycles

22. Question 4 ?? • When the stable torus vanishes it - • remains of finite size • is not distorted • its quasiperiod does not tend to infinity • Thus: There are no signs of the torus merging with some saddle • cycle, i.e. no homoclinic-like bifurcation. • 3. For a double check, such a stray cycle was sought for and • not found. Question 4: Why does the torus disappear?

23. Speculation The assumption that the phase space structure is the same as for a forced periodic oscillator does not explain the phenomena observed. We need to reveal the true phase space structure that explains all observed phenomena. We hypothesize the following structure

24. Hypothesis to Fit the Facts Inside locking region where two stable cycles coexist - • Fact: All stable and saddle cycles lie on the same stable torus. • 2. Fact: Any saddle cycle can merge with any twice saddle cycle. Suggestion: Saddle and twice saddle cycles could lie on the same saddle torus. • To imagine a saddle torus, draw analogy with a saddle cycle in 3D. • A saddle cycle is a closed curve in • the original 3D space. • A saddle torus is a closed curve in 3D • Poincare section of a 4D space. • A saddle cycle is an intersection • of two 2D manifolds in 3D space. • A saddle torus is an intersection of two • 2D manifolds in 3D Poincare section.

25. Hypothesis Continued Consequence: A stable and a saddle torus should lie on the same manifold and intersect at two saddle cycles. 3. Fact: trajectories go to the stable cycles along spirals in Poincare section. Cutting the cylinder and squeezing its boundaries to a point stable torus Hypothesized Poincare map of all objects in the phase space saddle torus Stable cycle St Saddle cycle Twice saddle cycle

26. The Hypothesis would Answer Question 3 Question 3: What separates the basins of attraction of two cycles in suppression region? Answer: The stable manifold of the saddle torus (blue one). suppression S1 S2

27. The Hypothesis would Answer Question 4 Question 4: At a non-linear coupling, the just born torus disappears while staying of finite size, smooth, with finite quasiperiod. Why does the torus disappear? S2 S1 Answer: The stable torus merges with saddle torus and vanishes via SN bifurcation.

28. How to Verify the Hypothesis? The hypothesis looks nice, but needs verification. Problem:There seem to be no well established numerical methods for plotting a saddle torus Solution: If the hypothesis is correct, what follows from it is correct, too. We can check if the consequences are true. What follows from the hypothesis (1): • The section of a Saddle torus is intersection of two 2D surfaces. • One of them is the “sphere”. • “Sphere” is an attracting object. • Inside the “sphere” there should be a repelling object – and • there is one! An unstable fixed point at the origin – a repeller!.

29. How to Verify the Hypothesis? • Reverse the time: • repeller becomes stable fixed point • the “sphere” will be the borderline • between the basins of attraction of the fixed point and infinity • Thus: The “sphere” can be calculated via calculation of basins of attraction in reversed time. • What follows from the hypothesis (2): • In locking region two stable cycles coexist • Their basins of attraction are separated by the • “plane” • Thus: The “plane” can be calculated via calculation of basins • of attraction of the two stable cycles in direct time. locking

30. How to Verify the Hypothesis? The saddle torus can be found as intersection of a “sphere” and the “plane”. Note: This method does not require the knowledge of where the saddle or the twice saddle cycles are located. If the saddle torus is found, and if it fits the hypothesis, all 4 cycles should lie on it! This is a test for the saddle torus validity.

31. Stable torus Saddle torus How Saddle Torus was Found Borderline between basins of stable cycles (“plane”) Stable cycle “Sphere” Saddle cycle Twice saddle cycle

32. Phase Space Structure Hypothesized 1 Numerically estimated locking

33. Phase Space Structure S2 Hypothesized 1 T2 T* S1 Numerically estimated 3 2

34. Surface of the “Sphere” (Poincare section)

35. If systems are slightly non-identical, how will their BD change? C1=0.1 C2 p

36. What if we Couple Other Oscillators? Van der Pol mutually coupled oscillators FitzHugh-Nagumo mutually coupled oscillators far from non-local bifurcations

37. Conclusions • The structure of the phase space in two mutually coupled 2D periodic oscillators, each being far from non-local bifurcations, is more complicated than in a forced oscillator. It involves a saddle torus. • Nonlinearity in coupling function influences the structure of synchronization region and bifurcations. E.g. it can induce saddle-node bifurcation of tori. • An approach is proposed to vizualize a saddle torus.

38. p-C1 Bifurcation Diagrams (SN lines) C2 p e1= e2 = 0.5 w01=1

39. p-C2 Bifurcation Diagrams C2 p e1= e2 = 0.5 w01=1

40. Some Basic Publications • V.I. Arnold “Geometrical Methods in the Theory of Ordinary Differential Equations” (New York, Springer, 1993). • P.S. Landa “Self-Oscillations in Systems with Finite Numbers of Degrees of Freedom” (Moscow, Nauka, 1980, in Russian). • A. Pikovsky, M. Rosenblum, J. Kurths “Synchronization: a Universal Concept in Nonlinear Science” (Cambridge University Press, 2001). • E. Mosekilde, Yu. Maistrenko, D. Postnov “Chaotic Synchronization. Application to Living Systems” (World Scientific, Singapore, Series A, v 42, 2002). • S. Boccaletti, J. Kurths, G. Osipov, D.L. Valladares, C.S. Zhou Phys. Rep. 366, 1 (2002). • V.S. Anishchenko, V.V. Astakhov, A.B. Neiman, T.E. Vadivasova, L. Schimansky-Geier “Nonlinear Dynamics of Chaotic and Stochastic Systems” (Springer, Berlin, Heidelberg 2002). • J. Guckenheimer, P. Holmes “Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields” (Springer-Verlag, New York, 1986).