Synchronization of chaos

Last updated October 31, 2022

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. 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 value 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 because 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.

Quantum chaos is a branch of physics which studies how chaotic classical dynamical systems can be described in terms of quantum theory. The primary question that quantum chaos seeks to answer is: "What is the relationship between quantum mechanics and classical chaos?" The correspondence principle states that classical mechanics is the classical limit of quantum mechanics, specifically in the limit as the ratio of Planck's constant to the action of the system tends to zero. If this is true, then there must be quantum mechanisms underlying classical chaos. If quantum mechanics does not demonstrate an exponential sensitivity to initial conditions, how can exponential sensitivity to initial conditions arise in classical chaos, which must be the correspondence principle limit of quantum mechanics?

Various types of stability may be discussed for the solutions of differential equations or difference equations describing dynamical systems. The most important type is that concerning the stability of solutions near to a point of equilibrium. This may be discussed by the theory of Aleksandr Lyapunov. In simple terms, if the solutions that start out near an equilibrium point $stay near forever, then is Lyapunov stable . More strongly, if is Lyapunov stable and all solutions that start out near converge to, then is asymptotically stable. The notion of exponential stability guarantees a minimal rate of decay, i.e., an estimate of how quickly the solutions converge. The idea of Lyapunov stability can be extended to infinite-dimensional manifolds, where it is known as structural stability, which concerns the behavior of different but "nearby" solutions to differential equations. Input-to-state stability (ISS) applies Lyapunov notions to systems with inputs.$

The competitive Lotka–Volterra equations are a simple model of the population dynamics of species competing for some common resource. They can be further generalised to the Generalized Lotka–Volterra equation to include trophic interactions.

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.

The Kuramoto model, first proposed by Yoshiki Kuramoto, is a mathematical model used to 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 gauge covariant derivative is a variation of the covariant derivative used in general relativity, quantum field theory and fluid dynamics. If a theory has gauge transformations, it means that some physical properties of certain equations are preserved under those transformations. Likewise, the gauge covariant derivative is the ordinary derivative modified in such a way as to make it behave like a true vector operator, so that equations written using the covariant derivative preserve their physical properties under gauge transformations.

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

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).

A coupled map lattice (CML) is a dynamical system that models the behavior of non-linear 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.

A geometric stable distribution or geo-stable distribution is a type of leptokurtic probability distribution. Geometric stable distributions were introduced in Klebanov, L. B., Maniya, G. M., and Melamed, I. A. (1985). A problem of Zolotarev and analogs of infinitely divisible and stable distributions in a scheme for summing a random number of random variables. These distributions are analogues for stable distributions for the case when the number of summands is random, independent of the distribution of summand, and having geometric distribution. The geometric stable distribution may be symmetric or asymmetric. A symmetric geometric stable distribution is also referred to as a Linnik distribution. The Laplace distribution and asymmetric Laplace distribution are special cases of the geometric stable distribution. The Mittag-Leffler distribution is also a special case of a geometric stable distribution.

A synchronization network is a network of coupled dynamical systems. It consists of a network connecting oscillators, where oscillators are nodes that emit a signal with somewhat regular frequency, and are also capable of receiving a signal.

In the fields of dynamical systems and control theory, a fractional-order system is a dynamical system that can be modeled by a fractional differential equation containing derivatives of non-integer order. Such systems are said to have fractional dynamics. Derivatives and integrals of fractional orders are used to describe objects that can be characterized by power-law nonlocality, power-law long-range dependence or fractal properties. Fractional-order systems are useful in studying the anomalous behavior of dynamical systems in physics, electrochemistry, biology, viscoelasticity and chaotic systems.

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

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).