Synchronization of chaos

Last updated June 13, 2024

Synchronization of chaos is a phenomenon that may occur when two or more dissipative chaotic systems are coupled.

Because of the exponential divergence of the nearby trajectories of chaotic systems, having two chaotic systems evolving in synchrony might appear surprising. However, synchronization of coupled or driven chaotic oscillators is a phenomenon well established experimentally and reasonably well-understood theoretically.

The stability of synchronization for coupled systems can be analyzed using master stability. Synchronization of chaos is a rich phenomenon and a multi-disciplinary subject with a broad range of applications.^[1]^[2]^[3]

Synchronization may present a variety of forms depending on the nature of the interacting systems and the type of coupling, and the proximity between the systems.

Identical synchronization

This type of synchronization is also known as complete synchronization. It can be observed for identical chaotic systems. The systems are said to be completely synchronized when there is a set of initial conditions so that the systems eventually evolve identically in time. In the simplest case of two diffusively coupled dynamics is described by

x^{\prime }=F(x)+\alpha (y-x)

y^{\prime }=F(y)+\alpha (x-y)

where $F$ is the vector field modeling the isolated chaotic dynamics and $\alpha$ is the coupling parameter. The regime $x(t)=y(t)$ defines an invariant subspace of the coupled system, if this subspace $x(t)=y(t)$ is locally attractive then the coupled system exhibit identical synchronization.

If the coupling vanishes the oscillators are decoupled, and the chaotic behavior leads to a divergence of nearby trajectories. Complete synchronization occurs due to the interaction, if the coupling parameter is large enough so that the divergence of trajectories of interacting systems due to chaos is suppressed by the diffusive coupling. To find the critical coupling strength we study the behavior of the difference $v=x-y$ . Assuming that $v$ is small we can expand the vector field in series and obtain a linear differential equation - by neglecting the Taylor remainder - governing the behavior of the difference

v^{\prime }=DF(x(t))v-2\alpha v

where $DF(x(t))$ denotes the Jacobian of the vector field along the solution. If $\alpha =0$ then we obtain

u^{\prime }=DF(x(t))u,

and since the dynamics of chaotic we have $\|u(t)\|\leq \|u(0)\|e^{\lambda t}$ , where $\lambda$ denotes the maximum Lyapunov exponent of the isolated system. Now using the ansatz $v=ue^{-2\alpha t}$ we pass from the equation for $v$ to the equation for $u$ . Therefore, we obtain

\|v(t)\|\leq \|u(0)\|e^{(-2\alpha +\lambda )t}

yield a critical coupling strength $\alpha _{c}=\lambda /2$ , for all $\alpha >\alpha _{c}$ the system exhibit complete synchronization. The existence of a critical coupling strength is related to the chaotic nature of the isolated dynamics.

In general, this reasoning leads to the correct critical coupling value for synchronization. However, in some cases one might observe loss of synchronization for coupling strengths larger than the critical value. This occurs because the nonlinear terms neglected in the derivation of the critical coupling value can play an important role and destroy the exponential bound for the behavior of the difference.^[4] It is however, possible to give a rigorous treatment to this problem and obtain a critical value so that the nonlinearities will not affect the stability.^[5]

Generalized synchronization

This type of synchronization occurs mainly when the coupled chaotic oscillators are different, although it has also been reported between identical oscillators. Given the dynamical variables $(x_{1},x_{2},...x_{n})$ and $(y_{1},y_{2},...y_{m})$ that determine the state of the oscillators, generalized synchronization occurs when there is a functional, $\phi$ , such that, after a transitory evolution from appropriate initial conditions, it is $[y_{1}(t),y_{2}(t),...,y_{m}(t)=\phi [x_{1}(t),x_{2}(t),...,x_{n}(t)]$ . This means that the dynamical state of one of the oscillators is completely determined by the state of the other. When the oscillators are mutually coupled this functional has to be invertible, if there is a drive-response configuration the drive determines the evolution of the response, and Φ does not need to be invertible. Identical synchronization is the particular case of generalized synchronization when $\phi$ is the identity.

Phase synchronization

Phase synchronization occurs when the coupled chaotic oscillators keep their phase difference bounded while their amplitudes remain uncorrelated. This phenomenon occurs even if the oscillators are not identical. Observation of phase synchronization requires a previous definition of the phase of a chaotic oscillator. In many practical cases, it is possible to find a plane in phase space in which the projection of the trajectories of the oscillator follows a rotation around a well-defined center. If this is the case, the phase is defined by the angle, φ(t), described by the segment joining the center of rotation and the projection of the trajectory point onto the plane. In other cases it is still possible to define a phase by means of techniques provided by the theory of signal processing, such as the Hilbert transform. In any case, if φ₁(t) and φ₂(t) denote the phases of the two coupled oscillators, synchronization of the phase is given by the relation nφ₁(t)=mφ₂(t) with m and n whole numbers.

Anticipated and lag synchronization

In these cases, the synchronized state is characterized by a time interval τ such that the dynamical variables of the oscillators, $(x_{1},x_{2},...x_{n})$ and $(x'_{1},x'_{2},...x'_{n})$ , are related by $x'_{i}(t)=x_{i}(t+\tau )$ ; this means that the dynamics of one of the oscillators follows, or anticipates, the dynamics of the other. Anticipated synchronization may occur between chaotic oscillators whose dynamics is described by delay differential equations, coupled in a drive-response configuration. In this case, the response anticipates the dynamics of the drive. Lag synchronization may occur when the strength of the coupling between phase-synchronized oscillators is increased.

Amplitude envelope synchronization

This is a mild form of synchronization that may appear between two weakly coupled chaotic oscillators. In this case, there is no correlation between phases nor amplitudes; instead, the oscillations of the two systems develop a periodic envelope that has the same frequency in the two systems.

This has the same order of magnitude than the difference between the average frequencies of oscillation of the two chaotic oscillator. Often, amplitude envelope synchronization precedes phase synchronization in the sense that when the strength of the coupling between two amplitude envelope synchronized oscillators is increased, phase synchronization develops.

All these forms of synchronization share the property of asymptotic stability. This means that once the synchronized state has been reached, the effect of a small perturbation that destroys synchronization is rapidly damped, and synchronization is recovered again. Mathematically, asymptotic stability is characterized by a positive Lyapunov exponent of the system composed of the two oscillators, which becomes negative when chaotic synchronization is achieved.

Some chaotic systems allow even stronger control of chaos, and both synchronization of chaos and control of chaos constitute parts of what's known as "cybernetical physics".

Notes

↑ Arenas, Alex; Díaz-Guilera, Albert; Kurths, Jurgen; Moreno, Yamir; Zhou, Changsong (2008-12-01). "Synchronization in complex networks". Physics Reports. 469 (3): 93–153. arXiv: 0805.2976 . Bibcode:2008PhR...469...93A. doi:10.1016/j.physrep.2008.09.002. S2CID 14355929.
↑ Wu, Chai Wah (2007). Synchronization in Complex Networks of Nonlinear Dynamical Systems. Bibcode:2007scnn.book.....W. doi:10.1142/6570. ISBN 978-981-270-973-8.
↑ Eroglu, Deniz; Lamb, Jeroen S. W.; Pereira, Tiago (2017). "Synchronisation of chaos and its applications". Contemporary Physics. 58 (3): 207–243. Bibcode:2017ConPh..58..207E. doi:10.1080/00107514.2017.1345844. hdl: 10044/1/53479 . ISSN 0010-7514. S2CID 126358436.
↑ Ashwin, Peter (2006-08-09). "Bubbling transition". Scholarpedia. 1 (8): 1725. Bibcode:2006SchpJ...1.1725A. doi: 10.4249/scholarpedia.1725 . ISSN 1941-6016.
↑ Tiago Pereira, Stability of Synchronized Motion in Complex Networks , arXiv:1112.2297v1, 2011.

Related Research Articles

In mathematics, a dynamical system is a system in which a function describes the time dependence of a point in an ambient space, such as in a parametric curve. Examples include the mathematical models that describe the swinging of a clock pendulum, the flow of water in a pipe, the random motion of particles in the air, and the number of fish each springtime in a lake. The most general definition unifies several concepts in mathematics such as ordinary differential equations and ergodic theory by allowing different choices of the space and how time is measured. Time can be measured by integers, by real or complex numbers or can be a more general algebraic object, losing the memory of its physical origin, and the space may be a manifold or simply a set, without the need of a smooth space-time structure defined on it.

<span class="mw-page-title-main">Oscillation</span> Repetitive variation of some measure about a central value

Oscillation is the repetitive or periodic variation, typically in time, of some measure about a central valu or between two or more different states. Familiar examples of oscillation include a swinging pendulum and alternating current. Oscillations can be used in physics to approximate complex interactions, such as those between atoms.

In mathematics and science, a nonlinear system is a system in which the change of the output is not proportional to the change of the input. Nonlinear problems are of interest to engineers, biologists, physicists, mathematicians, and many other scientists since most systems are inherently nonlinear in nature. Nonlinear dynamical systems, describing changes in variables over time, may appear chaotic, unpredictable, or counterintuitive, contrasting with much simpler linear systems.

<span class="mw-page-title-main">Lyapunov exponent</span> The rate of separation of infinitesimally close trajectories

In mathematics, the Lyapunov exponent or Lyapunov characteristic exponent of a dynamical system is a quantity that characterizes the rate of separation of infinitesimally close trajectories. Quantitatively, two trajectories in phase space with initial separation vector $diverge at a rate given by$

Flight dynamics is the science of air vehicle orientation and control in three dimensions. The three critical flight dynamics parameters are the angles of rotation in three dimensions about the vehicle's center of gravity (cg), known as pitch, roll and yaw. These are collectively known as aircraft attitude, often principally relative to the atmospheric frame in normal flight, but also relative to terrain during takeoff or landing, or when operating at low elevation. The concept of attitude is not specific to fixed-wing aircraft, but also extends to rotary aircraft such as helicopters, and dirigibles, where the flight dynamics involved in establishing and controlling attitude are entirely different.

In the study of dynamical systems, the van der Pol oscillator is a non-conservative, oscillating system with non-linear damping. It evolves in time according to the second-order differential equation

The Rabinovich–Fabrikant equations are a set of three coupled ordinary differential equations exhibiting chaotic behaviour for certain values of the parameters. They are named after Mikhail Rabinovich and Anatoly Fabrikant, who described them in 1979.

<span class="mw-page-title-main">Bifurcation theory</span> Study of sudden qualitative behavior changes caused by small parameter changes

Bifurcation theory is the mathematical study of changes in the qualitative or topological structure of a given family of curves, such as the integral curves of a family of vector fields, and the solutions of a family of differential equations. Most commonly applied to the mathematical study of dynamical systems, a bifurcation occurs when a small smooth change made to the parameter values of a system causes a sudden 'qualitative' or topological change in its behavior. Bifurcations occur in both continuous systems and discrete systems.

The Kuramoto model, first proposed by Yoshiki Kuramoto, is a mathematical model used in describing synchronization. More specifically, it is a model for the behavior of a large set of coupled oscillators. Its formulation was motivated by the behavior of systems of chemical and biological oscillators, and it has found widespread applications in areas such as neuroscience and oscillating flame dynamics. Kuramoto was quite surprised when the behavior of some physical systems, namely coupled arrays of Josephson junctions, followed his model.

The Duffing equation, named after Georg Duffing (1861–1944), is a non-linear second-order differential equation used to model certain damped and driven oscillators. The equation is given by

<span class="mw-page-title-main">Hopf bifurcation</span> Critical point where a periodic solution arises

In the mathematical theory of bifurcations, a Hopfbifurcation is a critical point where, as a parameter changes, a system's stability switches and a periodic solution arises. More accurately, it is a local bifurcation in which a fixed point of a dynamical system loses stability, as a pair of complex conjugate eigenvalues—of the linearization around the fixed point—crosses the complex plane imaginary axis as a parameter crosses a threshold value. Under reasonably generic assumptions about the dynamical system, the fixed point becomes a small-amplitude limit cycle as the parameter changes.

In mathematics, delay differential equations (DDEs) are a type of differential equation in which the derivative of the unknown function at a certain time is given in terms of the values of the function at previous times. DDEs are also called time-delay systems, systems with aftereffect or dead-time, hereditary systems, equations with deviating argument, or differential-difference equations. They belong to the class of systems with the functional state, i.e. partial differential equations (PDEs) which are infinite dimensional, as opposed to ordinary differential equations (ODEs) having a finite dimensional state vector. Four points may give a possible explanation of the popularity of DDEs:

Aftereffect is an applied problem: it is well known that, together with the increasing expectations of dynamic performances, engineers need their models to behave more like the real process. Many processes include aftereffect phenomena in their inner dynamics. In addition, actuators, sensors, and communication networks that are now involved in feedback control loops introduce such delays. Finally, besides actual delays, time lags are frequently used to simplify very high order models. Then, the interest for DDEs keeps on growing in all scientific areas and, especially, in control engineering.
Delay systems are still resistant to many classical controllers: one could think that the simplest approach would consist in replacing them by some finite-dimensional approximations. Unfortunately, ignoring effects which are adequately represented by DDEs is not a general alternative: in the best situation, it leads to the same degree of complexity in the control design. In worst cases, it is potentially disastrous in terms of stability and oscillations.
Voluntary introduction of delays can benefit the control system.
In spite of their complexity, DDEs often appear as simple infinite-dimensional models in the very complex area of partial differential equations (PDEs).

In the mathematics of evolving systems, the concept of a center manifold was originally developed to determine stability of degenerate equilibria. Subsequently, the concept of center manifolds was realised to be fundamental to mathematical modelling.

A synchronous frame is a reference frame in which the time coordinate defines proper time for all co-moving observers. It is built by choosing some constant time hypersurface as an origin, such that has in every point a normal along the time line and a light cone with an apex in that point can be constructed; all interval elements on this hypersurface are space-like. A family of geodesics normal to this hypersurface are drawn and defined as the time coordinates with a beginning at the hypersurface. In terms of metric-tensor components $, a synchronous frame is defined such that$

A coupled map lattice (CML) is a dynamical system that models the behavior of nonlinear systems. They are predominantly used to qualitatively study the chaotic dynamics of spatially extended systems. This includes the dynamics of spatiotemporal chaos where the number of effective degrees of freedom diverges as the size of the system increases.

Undulatory locomotion is the type of motion characterized by wave-like movement patterns that act to propel an animal forward. Examples of this type of gait include crawling in snakes, or swimming in the lamprey. Although this is typically the type of gait utilized by limbless animals, some creatures with limbs, such as the salamander, forgo use of their legs in certain environments and exhibit undulatory locomotion. In robotics this movement strategy is studied in order to create novel robotic devices capable of traversing a variety of environments.

In mathematics, the master stability function is a tool used to analyze the stability of the synchronous state in a dynamical system consisting of many identical systems which are coupled together, such as the Kuramoto model.

Supersymmetric theory of stochastic dynamics or stochastics (STS) is an exact theory of stochastic (partial) differential equations (SDEs), the class of mathematical models with the widest applicability covering, in particular, all continuous time dynamical systems, with and without noise. The main utility of the theory from the physical point of view is a rigorous theoretical explanation of the ubiquitous spontaneous long-range dynamical behavior that manifests itself across disciplines via such phenomena as 1/f, flicker, and crackling noises and the power-law statistics, or Zipf's law, of instantonic processes like earthquakes and neuroavalanches. From the mathematical point of view, STS is interesting because it bridges the two major parts of mathematical physics – the dynamical systems theory and topological field theories. Besides these and related disciplines such as algebraic topology and supersymmetric field theories, STS is also connected with the traditional theory of stochastic differential equations and the theory of pseudo-Hermitian operators.

The Dicke model is a fundamental model of quantum optics, which describes the interaction between light and matter. In the Dicke model, the light component is described as a single quantum mode, while the matter is described as a set of two-level systems. When the coupling between the light and matter crosses a critical value, the Dicke model shows a mean-field phase transition to a superradiant phase. This transition belongs to the Ising universality class and was realized in cavity quantum electrodynamics experiments. Although the superradiant transition bears some analogy with the lasing instability, these two transitions belong to different universality classes.

<span class="mw-page-title-main">Heteroclinic channels</span> Robotic control method

Heteroclinic channels are ensembles of trajectories that can connect saddle equilibrium points in phase space. Dynamical systems and their associated phase spaces can be used to describe natural phenomena in mathematical terms; heteroclinic channels, and the cycles that they produce, are features in phase space that can be designed to occupy specific locations in that space. Heteroclinic channels move trajectories from one equilibrium point to another. More formally, a heteroclinic channel is a region in phase space in which nearby trajectories are drawn closer and closer to one unique limiting trajectory, the heteroclinic orbit. Equilibria connected by heteroclinic trajectories form heteroclinic cycles and cycles can be connected to form heteroclinic networks. Heteroclinic cycles and networks naturally appear in a number of applications, such as fluid dynamics, population dynamics, and neural dynamics. In addition, dynamical systems are often used as methods for robotic control. In particular, for robotic control, the equilibrium points can correspond to robotic states, and the heteroclinic channels can provide smooth methods for switching from state to state.

References

Pikovsky, A.; Rosemblum, M.; Kurths, J. (2001). Synchronization: A Universal Concept in Nonlinear Sciences. Cambridge University Press. ISBN 978-0-521-53352-2.
González-Miranda, J. M. (2004). Synchronization and Control of Chaos. An introduction for scientists and engineers. Imperial College Press. ISBN 978-1-86094-488-8.
Fradkov A.L. Cybernetical physics: from control of chaos to quantum control. Springer-Verlag, 2007, (Preliminary Russian version: St.Petersburg, Nauka, 2003).