A coupled map lattice (CML) is a dynamical system that models the behavior of nonlinear systems (especially partial differential equations). 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. [1]
Features of the CML are discrete time dynamics, discrete underlying spaces (lattices or networks), and real (number or vector), local, continuous state variables. [2] Studied systems include populations, chemical reactions, convection, fluid flow and biological networks. More recently, CMLs have been applied to computational networks [3] identifying detrimental attack methods and cascading failures.
CMLs are comparable to cellular automata models in terms of their discrete features. [4] However, the value of each site in a cellular automata network is strictly dependent on its neighbor(s) from the previous time step. Each site of the CML is only dependent upon its neighbors relative to the coupling term in the recurrence equation. However, the similarities can be compounded when considering multi-component dynamical systems.
A CML generally incorporates a system of equations (coupled or uncoupled), a finite number of variables, a global or local coupling scheme and the corresponding coupling terms. The underlying lattice can exist in infinite dimensions. Mappings of interest in CMLs generally demonstrate chaotic behavior. Such maps can be found here: List of chaotic maps.
A logistic mapping demonstrates chaotic behavior, easily identifiable in one dimension for parameter r > 3.57:
In Figure 1, is initialized to random values across a small lattice; the values are decoupled with respect to neighboring sites. The same recurrence relation is applied at each lattice point, although the parameter r is slightly increased with each time step. The result is a raw form of chaotic behavior in a map lattice. However, there are no significant spatial correlations or pertinent fronts to the chaotic behavior. No obvious order is apparent.
For a basic coupling, we consider a 'single neighbor' coupling where the value at any given site is computed from the recursive maps both on itself and on the neighboring site . The coupling parameter is equally weighted. Again, the value of is constant across the lattice, but slightly increased with each time step.
Even though the recursion is chaotic, a more solid form develops in the evolution. Elongated convective spaces persist throughout the lattice (see Figure 2).
Figure 1: An uncoupled logistic map lattice with random seeding over forty iterations. | Figure 2: A CML with a single-neighbor coupling scheme taken over forty iterations. |
CMLs were first introduced in the mid 1980s through a series of closely released publications. [5] [6] [7] [8] Kapral used CMLs for modeling chemical spatial phenomena. Kuznetsov sought to apply CMLs to electrical circuitry by developing a renormalization group approach (similar to Feigenbaum's universality to spatially extended systems). Kaneko's focus was more broad and he is still known as the most active researcher in this area. [9] The most examined CML model was introduced by Kaneko in 1983 where the recurrence equation is as follows:
where and is a real mapping.
The applied CML strategy was as follows:
The CML system evolves through discrete time by a mapping on vector sequences. These mappings are a recursive function of two competing terms: an individual non-linear reaction, and a spatial interaction (coupling) of variable intensity. CMLs can be classified by the strength of this coupling parameter(s).
Much of the current published work in CMLs is based in weak coupled systems [2] where diffeomorphisms of the state space close to identity are studied. Weak coupling with monotonic (bistable) dynamical regimes demonstrate spatial chaos phenomena and are popular in neural models. [10] Weak coupling unimodal maps are characterized by their stable periodic points and are used by gene regulatory network models. Space-time chaotic phenomena can be demonstrated from chaotic mappings subject to weak coupling coefficients and are popular in phase transition phenomena models.
Intermediate and strong coupling interactions are less prolific areas of study. Intermediate interactions are studied with respect to fronts and traveling waves, riddled basins, riddled bifurcations, clusters and non-unique phases. Strong coupling interactions are most well known to model synchronization effects of dynamic spatial systems such as the Kuramoto model.
These classifications do not reflect the local or global (GMLs [11] ) coupling nature of the interaction. Nor do they consider the frequency of the coupling which can exist as a degree of freedom in the system. [12] Finally, they do not distinguish between sizes of the underlying space or boundary conditions.
Surprisingly the dynamics of CMLs have little to do with the local maps that constitute their elementary components. With each model a rigorous mathematical investigation is needed to identify a chaotic state (beyond visual interpretation). Rigorous proofs have been performed to this effect. By example: the existence of space-time chaos in weak space interactions of one-dimensional maps with strong statistical properties was proven by Bunimovich and Sinai in 1988. [13] Similar proofs exist for weakly coupled hyperbolic maps under the same conditions.
For the case where the underlying map is based on a generalised Bernoulli map it can be shown that the complete Lyapunov spectrum for the CML can be evaluated analytically in a range of cases. [14]
CMLs have revealed novel qualitative universality classes in (CML) phenomenology. Such classes include:
The unique qualitative classes listed above can be visualized. By applying the Kaneko 1983 model to the logistic map, several of the CML qualitative classes may be observed. These are demonstrated below, note the unique parameters:
Frozen Chaos | Pattern Selection | Chaotic Brownian Motion of Defect |
Figure 1: Sites are divided into non-uniform clusters, where the divided patterns are regarded as attractors. Sensitivity to initial conditions exist relative to a < 1.5. | Figure 2: Near uniform sized clusters (a = 1.71, ε = 0.4). | Figure 3: Defects exist in the system and fluctuate chaotically akin to Brownian motion (a = 1.85, ε = 0.1). |
Defect Turbulence | Spatiotemporal Intermittency I | Spatiotemporal Intermittency II |
Figure 4: Many defects are generated and turbulently collide (a = 1.895, ε = 0.1). | Figure 5: Each site transits between a coherent state and chaotic state intermittently (a = 1.75, ε = 0.6), Phase I. | Figure 6: The coherent state, Phase II. |
Fully Developed Spatiotemporal Chaos | Traveling Wave | |
Figure 7: Most sites independently oscillate chaotically (a = 2.00, ε = 0.3). | Figure 8: The wave of clusters travels at 'low' speeds (a = 1.47, ε = 0.5). |
Coupled map lattices being a prototype of spatially extended systems easy to simulate have represented a benchmark for the definition and introduction of many indicators of spatio-temporal chaos, the most relevant ones are
Chaos theory is an interdisciplinary area of scientific study and branch of mathematics. It focuses on underlying patterns and deterministic laws of dynamical systems that are highly sensitive to initial conditions. These were once thought to have completely random states of disorder and irregularities. Chaos theory states that within the apparent randomness of chaotic complex systems, there are underlying patterns, interconnection, constant feedback loops, repetition, self-similarity, fractals and self-organization. The butterfly effect, an underlying principle of chaos, describes how a small change in one state of a deterministic nonlinear system can result in large differences in a later state. A metaphor for this behavior is that a butterfly flapping its wings in Brazil can cause a tornado in Texas.
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.
The logistic map is a polynomial mapping of degree 2, often referred to as an archetypal example of how complex, chaotic behaviour can arise from very simple nonlinear dynamical equations. The map, initially utilized by Edward Lorenz in the 1960s to showcase irregular solutions, was popularized in a 1976 paper by the biologist Robert May, in part as a discrete-time demographic model analogous to the logistic equation written down by Pierre François Verhulst. Mathematically, the logistic map is written
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.
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
In descriptive statistics and chaos theory, a recurrence plot (RP) is a plot showing, for each moment in time, the times at which the state of a dynamical system returns to the previous state at , i.e., when the phase space trajectory visits roughly the same area in the phase space as at time . In other words, it is a plot of
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.
In dynamical systems theory, the baker's map is a chaotic map from the unit square into itself. It is named after a kneading operation that bakers apply to dough: the dough is cut in half, and the two halves are stacked on one another, and compressed.
In dynamical systems theory, a period-doubling bifurcation occurs when a slight change in a system's parameters causes a new periodic trajectory to emerge from an existing periodic trajectory—the new one having double the period of the original. With the doubled period, it takes twice as long for the numerical values visited by the system to repeat themselves.
The Lorenz system is a system of ordinary differential equations first studied by mathematician and meteorologist Edward Lorenz. It is notable for having chaotic solutions for certain parameter values and initial conditions. In particular, the Lorenz attractor is a set of chaotic solutions of the Lorenz system. The term "butterfly effect" in popular media may stem from the real-world implications of the Lorenz attractor, namely that tiny changes in initial conditions evolve to completely different trajectories. This underscores that chaotic systems can be completely deterministic and yet still be inherently impractical or even impossible to predict over longer periods of time. For example, even the small flap of a butterfly's wings could set the earth's atmosphere on a vastly different trajectory, in which for example a hurricane occurs where it otherwise would have not. The shape of the Lorenz attractor itself, when plotted in phase space, may also be seen to resemble a butterfly.
Synchronization of chaos is a phenomenon that may occur when two or more dissipative chaotic systems are coupled.
Stellar pulsations are caused by expansions and contractions in the outer layers as a star seeks to maintain equilibrium. These fluctuations in stellar radius cause corresponding changes in the luminosity of the star. Astronomers are able to deduce this mechanism by measuring the spectrum and observing the Doppler effect. Many intrinsic variable stars that pulsate with large amplitudes, such as the classical Cepheids, RR Lyrae stars and large-amplitude Delta Scuti stars show regular light curves.
Valentin Afraimovich was a Soviet, Russian and Mexican mathematician. He made contributions to dynamical systems theory, qualitative theory of ordinary differential equations, bifurcation theory, concept of attractor, strange attractors, space-time chaos, mathematical models of non-equilibrium media and biological systems, travelling waves in lattices, complexity of orbits and dimension-like characteristics in dynamical systems.
The kicked rotator, also spelled as kicked rotor, is a paradigmatic model for both Hamiltonian chaos and quantum chaos. It describes a free rotating stick in an inhomogeneous "gravitation like" field that is periodically switched on in short pulses. The model is described by the Hamiltonian
The dynamical systems approach to neuroscience is a branch of mathematical biology that utilizes nonlinear dynamics to understand and model the nervous system and its functions. In a dynamical system, all possible states are expressed by a phase space. Such systems can experience bifurcation as a function of its bifurcation parameters and often exhibit chaos. Dynamical neuroscience describes the non-linear dynamics at many levels of the brain from single neural cells to cognitive processes, sleep states and the behavior of neurons in large-scale neuronal simulation.
Chaotic scattering is a branch of chaos theory dealing with scattering systems displaying a strong sensitivity to initial conditions. In a classical scattering system there will be one or more impact parameters, b, in which a particle is sent into the scatterer. This gives rise to one or more exit parameters, y, as the particle exits towards infinity. While the particle is traversing the system, there may also be a delay time, T—the time it takes for the particle to exit the system—in addition to the distance travelled, s. In certain systems the two will be equivalent—see below. In a chaotic scattering system, a minute change in the impact parameter, may give rise to a very large change in the exit parameters.
The Rulkov map is a two-dimensional iterated map used to model a biological neuron. It was proposed by Nikolai F. Rulkov in 2001. The use of this map to study neural networks has computational advantages because the map is easier to iterate than a continuous dynamical system. This saves memory and simplifies the computation of large neural networks.
The Frenkel–Kontorova (FK) model is a fundamental model of low-dimensional nonlinear physics.
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 Chialvo map is a two-dimensional map proposed by Dante R. Chialvo in 1995 to describe the generic dynamics of excitable systems. The model is inspired by Kunihiko Kaneko's Coupled map lattice numerical approach which considers time and space as discrete variables but state as a continuous one. Later on Rulkov popularized a similar approach. By using only three parameters the model is able to efficiently mimic generic neuronal dynamics in computational simulations, as single elements or as parts of inter-connected networks.