Attractor

For other uses, see Attractor (disambiguation).

"Strange attractor" redirects here. For other uses, see Strange Attractor (disambiguation).

Visual representation of a strange attractor. Another visualization of the same 3D attractor is this video. Code capable of rendering this is available.

In the mathematical field of dynamical systems, an attractor is a set of states toward which a system tends to evolve, ^[2] for a wide variety of starting conditions of the system. System values that get close enough to the attractor values remain close even if slightly disturbed.

In finite-dimensional systems, the evolving variable may be represented algebraically as an n-dimensional vector. The attractor is a region in n-dimensional space. In physical systems, the n dimensions may be, for example, two or three positional coordinates for each of one or more physical entities; in economic systems, they may be separate variables such as the inflation rate and the unemployment rate.^{[ not verified in body ]}

If the evolving variable is two- or three-dimensional, the attractor of the dynamic process can be represented geometrically in two or three dimensions, (as for example in the three-dimensional case depicted to the right). An attractor can be a point, a finite set of points, a curve, a manifold, or even a complicated set with a fractal structure known as a strange attractor (see strange attractor below). If the variable is a scalar, the attractor is a subset of the real number line. Describing the attractors of chaotic dynamical systems has been one of the achievements of chaos theory.

A trajectory of the dynamical system in the attractor does not have to satisfy any special constraints except for remaining on the attractor, forward in time. The trajectory may be periodic or chaotic. If a set of points is periodic or chaotic, but the flow in the neighborhood is away from the set, the set is not an attractor, but instead is called a repeller (or repellor).

Motivation of attractors

A dynamical system is generally described by one or more differential or difference equations. The equations of a given dynamical system specify its behavior over any given short period of time. To determine the system's behavior for a longer period, it is often necessary to integrate the equations, either through analytical means or through iteration, often with the aid of computers.

Dynamical systems in the physical world tend to arise from dissipative systems: if it were not for some driving force, the motion would cease. (Dissipation may come from internal friction, thermodynamic losses, or loss of material, among many causes.) The dissipation and the driving force tend to balance, killing off initial transients and settle the system into its typical behavior. The subset of the phase space of the dynamical system corresponding to the typical behavior is the attractor, also known as the attracting section or attractee.

Invariant sets and limit sets are similar to the attractor concept. An invariant set is a set that evolves to itself under the dynamics. ^[3] Attractors may contain invariant sets. A limit set is a set of points such that there exists some initial state that ends up arbitrarily close to the limit set (i.e. to each point of the set) as time goes to infinity. Attractors are limit sets, but not all limit sets are attractors: It is possible to have some points of a system converge to a limit set, but different points when perturbed slightly off the limit set may get knocked off and never return to the vicinity of the limit set.

For example, the damped pendulum has two invariant points: the point x₀ of minimum height and the point x₁ of maximum height. The point x₀ is also a limit set, as trajectories converge to it; the point x₁ is not a limit set. Because of the dissipation due to air resistance, the point x₀ is also an attractor. If there was no dissipation, x₀ would not be an attractor. Aristotle believed that objects moved only as long as they were pushed, which is an early formulation of a dissipative attractor.

Some attractors are known to be chaotic (see strange attractor), in which case the evolution of any two distinct points of the attractor result in exponentially diverging trajectories, which complicates prediction when even the smallest noise is present in the system. ^[4]

Mathematical definition

Let ${\displaystyle t}$ represent time and let ${\displaystyle f(t,\cdot )}$ be a function which specifies the dynamics of the system. That is, if ${\displaystyle a}$ is a point in an ${\displaystyle n}$-dimensional phase space, representing the initial state of the system, then ${\displaystyle f(0,a)=a}$ and, for a positive value of ${\displaystyle t}$, ${\displaystyle f(t,a)}$ is the result of the evolution of this state after ${\displaystyle t}$ units of time. For example, if the system describes the evolution of a free particle in one dimension then the phase space is the plane ${\displaystyle \mathbb {R} ^{2}}$ with coordinates ${\displaystyle (x,v)}$, where ${\displaystyle x}$ is the position of the particle, ${\displaystyle v}$ is its velocity, ${\displaystyle a=(x,v)}$, and the evolution is given by

Attracting period-3 cycle and its immediate basin of attraction for a certain parametrization of the Julia set, which iterates the function f(z) = z + c. The three darkest points are the points of the 3-cycle, which lead to each other in sequence, and iteration from any point in the basin of attraction leads to (usually asymptotic) convergence to this sequence of three points.

${\displaystyle f(t,(x,v))=(x+tv,v).\ }$

An attractor is a subset ${\displaystyle A}$ of the phase space characterized by the following three conditions:

• ${\displaystyle A}$ is forward invariant under ${\displaystyle f}$: if ${\displaystyle a}$ is an element of ${\displaystyle A}$ then so is ${\displaystyle f(t,a)}$, for all ${\displaystyle t>0}$.
• There exists a neighborhood of ${\displaystyle A}$, called the basin of attraction for ${\displaystyle A}$ and denoted ${\displaystyle B(A)}$, which consists of all points ${\displaystyle b}$ that "enter" ${\displaystyle A}$ in the limit ${\displaystyle t\to \infty }$. More formally, ${\displaystyle B(A)}$ is the set of all points ${\displaystyle b}$ in the phase space with the following property:

For any open neighborhood ${\displaystyle N}$ of ${\displaystyle A}$, there is a positive constant ${\displaystyle T}$ such that ${\displaystyle f(t,b)\in N}$ for all real ${\displaystyle t>T}$.

• There is no proper (non-empty) subset of ${\displaystyle A}$ having the first two properties.

Since the basin of attraction contains an open set containing ${\displaystyle A}$, every point that is sufficiently close to ${\displaystyle A}$ is attracted to ${\displaystyle A}$. The definition of an attractor uses a metric on the phase space, but the resulting notion usually depends only on the topology of the phase space. In the case of ${\displaystyle \mathbb {R} ^{n}}$, the Euclidean norm is typically used.

Many other definitions of attractor occur in the literature. For example, some authors require that an attractor have positive measure (preventing a point from being an attractor), others relax the requirement that ${\displaystyle B(A)}$ be a neighborhood. ^[5]

Types of attractors

Attractors are portions or subsets of the phase space of a dynamical system. Until the 1960s, attractors were thought of as being simple geometric subsets of the phase space, like points, lines, surfaces, and simple regions of three-dimensional space. More complex attractors that cannot be categorized as simple geometric subsets, such as topologically wild sets, were known of at the time but were thought to be fragile anomalies. Stephen Smale was able to show that his horseshoe map was robust and that its attractor had the structure of a Cantor set.

Two simple attractors are a fixed point and the limit cycle. Attractors can take on many other geometric shapes (phase space subsets). But when these sets (or the motions within them) cannot be easily described as simple combinations (e.g. intersection and union) of fundamental geometric objects (e.g. lines, surfaces, spheres, toroids, manifolds), then the attractor is called a strange attractor .

Fixed point

Weakly attracting fixed point for a complex number evolving according to a complex quadratic polynomial. The phase space is the horizontal complex plane; the vertical axis measures the frequency with which points in the complex plane are visited. The point in the complex plane directly below the peak frequency is the fixed point attractor.

A fixed point of a function or transformation is a point that is mapped to itself by the function or transformation. If we regard the evolution of a dynamical system as a series of transformations, then there may or may not be a point which remains fixed under each transformation. The final state that a dynamical system evolves towards corresponds to an attracting fixed point of the evolution function for that system, such as the center bottom position of a damped pendulum, the level and flat water line of sloshing water in a glass, or the bottom center of a bowl containing a rolling marble. But the fixed point(s) of a dynamic system is not necessarily an attractor of the system. For example, if the bowl containing a rolling marble was inverted and the marble was balanced on top of the bowl, the center bottom (now top) of the bowl is a fixed state, but not an attractor. This is equivalent to the difference between stable and unstable equilibria. In the case of a marble on top of an inverted bowl (a hill), that point at the top of the bowl (hill) is a fixed point (equilibrium), but not an attractor (unstable equilibrium).

In addition, physical dynamic systems with at least one fixed point invariably have multiple fixed points and attractors due to the reality of dynamics in the physical world, including the nonlinear dynamics of stiction, friction, surface roughness, deformation (both elastic and plasticity), and even quantum mechanics. ^[6] In the case of a marble on top of an inverted bowl, even if the bowl seems perfectly hemispherical, and the marble's spherical shape, are both much more complex surfaces when examined under a microscope, and their shapes change or deform during contact. Any physical surface can be seen to have a rough terrain of multiple peaks, valleys, saddle points, ridges, ravines, and plains. ^[7] There are many points in this surface terrain (and the dynamic system of a similarly rough marble rolling around on this microscopic terrain) that are considered stationary or fixed points, some of which are categorized as attractors.

Finite number of points

In a discrete-time system, an attractor can take the form of a finite number of points that are visited in sequence. Each of these points is called a periodic point. This is illustrated by the logistic map, which depending on its specific parameter value can have an attractor consisting of 1 point, 2 points, 2ⁿ points, 3 points, 3×2ⁿ points, 4 points, 5 points, or any given positive integer number of points.

Limit cycle

Main article: Limit cycle

A limit cycle is a periodic orbit of a continuous dynamical system that is isolated. It concerns a cyclic attractor. Examples include the swings of a pendulum clock, and the heartbeat while resting. The limit cycle of an ideal pendulum is not an example of a limit cycle attractor because its orbits are not isolated: in the phase space of the ideal pendulum, near any point of a periodic orbit there is another point that belongs to a different periodic orbit, so the former orbit is not attracting. For a physical pendulum under friction, the resting state will be a fixed-point attractor. The difference with the clock pendulum is that there, energy is injected by the escapement mechanism to maintain the cycle.

Van der Pol phase portrait: an attracting limit cycle

Limit torus

There may be more than one frequency in the periodic trajectory of the system through the state of a limit cycle. For example, in physics, one frequency may dictate the rate at which a planet orbits a star while a second frequency describes the oscillations in the distance between the two bodies. If two of these frequencies form an irrational fraction (i.e. they are incommensurate), the trajectory is no longer closed, and the limit cycle becomes a limit torus. This kind of attractor is called an N_t -torus if there are N_t incommensurate frequencies. For example, here is a 2-torus:

A time series corresponding to this attractor is a quasiperiodic series: A discretely sampled sum of N_t periodic functions (not necessarily sine waves) with incommensurate frequencies. Such a time series does not have a strict periodicity, but its power spectrum still consists only of sharp lines.^{[ citation needed ]}

Strange attractor

A plot of Lorenz's strange attractor for values ρ = 28, σ = 10, β = 8/3

An attractor is called strange if it has a fractal structure, that is if it has non-integer Hausdorff dimension. This is often the case when the dynamics on it are chaotic, but strange nonchaotic attractors also exist. If a strange attractor is chaotic, exhibiting sensitive dependence on initial conditions, then any two arbitrarily close alternative initial points on the attractor, after any of various numbers of iterations, will lead to points that are arbitrarily far apart (subject to the confines of the attractor), and after any of various other numbers of iterations will lead to points that are arbitrarily close together. Thus a dynamic system with a chaotic attractor is locally unstable yet globally stable: once some sequences have entered the attractor, nearby points diverge from one another but never depart from the attractor. ^[8]

The term strange attractor was coined by David Ruelle and Floris Takens to describe the attractor resulting from a series of bifurcations of a system describing fluid flow. ^[9] Strange attractors are often differentiable in a few directions, but some are like a Cantor dust, and therefore not differentiable. Strange attractors may also be found in the presence of noise, where they may be shown to support invariant random probability measures of Sinai–Ruelle–Bowen type. ^[10]

Examples of strange attractors include the double-scroll attractor, Hénon attractor, Rössler attractor, and Lorenz attractor.

Attractors characterize the evolution of a system

Bifurcation diagram of the logistic map. The attractor(s) for any value of the parameter ${\displaystyle r}$ are shown on the ordinate in the domain ${\displaystyle 0<x<1}$. The colour of a point indicates how often the point ${\displaystyle (r,x)}$ is visited over the course of 10 iterations: frequently encountered values are coloured in blue, less frequently encountered values are yellow. A bifurcation appears around ${\displaystyle r\approx 3.0}$, a second bifurcation (leading to four attractor values) around ${\displaystyle r\approx 3.5}$. The behaviour is increasingly complicated for ${\displaystyle r>3.6}$, interspersed with regions of simpler behaviour (white stripes).

The parameters of a dynamic equation evolve as the equation is iterated, and the specific values may depend on the starting parameters. An example is the well-studied logistic map, ${\displaystyle x_{n+1}=rx_{n}(1-x_{n})}$, whose basins of attraction for various values of the parameter ${\displaystyle r}$ are shown in the figure. If ${\displaystyle r=2.6}$, all starting ${\displaystyle x}$ values of ${\displaystyle x<0}$ will rapidly lead to function values that go to negative infinity; starting ${\displaystyle x}$ values of ${\displaystyle x>1}$ will also go to negative infinity. But for ${\displaystyle 0<x<1}$ the ${\displaystyle x}$ values rapidly converge to ${\displaystyle x\approx 0.615}$, i.e. at this value of ${\displaystyle r}$, a single value of ${\displaystyle x}$ is an attractor for the function's behaviour. For other values of ${\displaystyle r}$, more than one value of ${\displaystyle x}$ may be visited: if ${\displaystyle r}$ is 3.2, starting values of ${\displaystyle 0<x<1}$ will lead to function values that alternate between ${\displaystyle x\approx 0.513}$ and ${\displaystyle x\approx 0.799}$. At some values of ${\displaystyle r}$, the attractor is a single point (a "fixed point"), at other values of ${\displaystyle r}$ two values of ${\displaystyle x}$ are visited in turn (a period-doubling bifurcation), or, as a result of further doubling, any number ${\displaystyle k\times 2^{n}}$ values of ${\displaystyle x}$; at yet other values of ${\displaystyle r}$, any given number of values of ${\displaystyle x}$ are visited in turn; finally, for some values of ${\displaystyle r}$, an infinitude of points are visited. Thus one and the same dynamic equation can have various types of attractors, depending on its parameters.

Basins of attraction

An attractor's basin of attraction is the region of the phase space, over which iterations are defined, such that any point (any initial condition) in that region will asymptotically be iterated into the attractor. For a stable linear system, every point in the phase space is in the basin of attraction. However, in nonlinear systems, some points may map directly or asymptotically to infinity, while other points may lie in a different basin of attraction and map asymptotically into a different attractor; other initial conditions may be in or map directly into a non-attracting point or cycle. ^[11]

Linear equation or system

An univariate linear homogeneous difference equation ${\displaystyle x_{t}=ax_{t-1}}$ diverges to infinity if ${\displaystyle |a|>1}$ from all initial points except 0; there is no attractor and therefore no basin of attraction. But if ${\displaystyle |a|<1}$ all points on the number line map asymptotically (or directly in the case of 0) to 0; 0 is the attractor, and the entire number line is the basin of attraction.

Likewise, a linear matrix difference equation in a dynamic vector ${\displaystyle X}$, of the homogeneous form ${\displaystyle X_{t}=AX_{t-1}}$ in terms of square matrix ${\displaystyle A}$ will have all elements of the dynamic vector diverge to infinity if the largest eigenvalues of ${\displaystyle A}$ is greater than 1 in absolute value; there is no attractor and no basin of attraction. But if the largest eigenvalue is less than 1 in magnitude, all initial vectors will asymptotically converge to the zero vector, which is the attractor; the entire ${\displaystyle n}$-dimensional space of potential initial vectors is the basin of attraction.

Similar features apply to linear differential equations. The scalar equation ${\displaystyle dx/dt=ax}$ causes all initial values of ${\displaystyle x}$ except zero to diverge to infinity if ${\displaystyle a>0}$ but to converge to an attractor at the value 0 if ${\displaystyle a<0}$, making the entire number line the basin of attraction for 0. And the matrix system ${\displaystyle dX/dt=AX}$ gives divergence from all initial points except the vector of zeroes if any eigenvalue of the matrix ${\displaystyle A}$ is positive; but if all the eigenvalues are negative the vector of zeroes is an attractor whose basin of attraction is the entire phase space.

Nonlinear equation or system

Equations or systems that are nonlinear can give rise to a richer variety of behavior than can linear systems. One example is Newton's method of iterating to a root of a nonlinear expression. If the expression has more than one real root, some starting points for the iterative algorithm will lead to one of the roots asymptotically, and other starting points will lead to another. The basins of attraction for the expression's roots are generally not simple—it is not simply that the points nearest one root all map there, giving a basin of attraction consisting of nearby points. The basins of attraction can be infinite in number and arbitrarily small. For example, ^[12] for the function ${\displaystyle f(x)=x^{3}-2x^{2}-11x+12}$, the following initial conditions are in successive basins of attraction:

A Newton fractal showing basins of attraction in the complex plane for using Newton's method to solve x − 1 = 0. Points in like-colored regions map to the same root; darker means more iterations are needed to converge.

2.35287527 converges to 4;

2.35284172 converges to −3;

2.35283735 converges to 4;

2.352836327 converges to −3;

2.352836323 converges to 1.

Newton's method can also be applied to complex functions to find their roots. Each root has a basin of attraction in the complex plane; these basins can be mapped as in the image shown. As can be seen, the combined basin of attraction for a particular root can have many disconnected regions. For many complex functions, the boundaries of the basins of attraction are fractals.

Partial differential equations

Parabolic partial differential equations may have finite-dimensional attractors. The diffusive part of the equation damps higher frequencies and in some cases leads to a global attractor. The Ginzburg–Landau, the Kuramoto–Sivashinsky, and the two-dimensional, forced Navier–Stokes equations are all known to have global attractors of finite dimension.

For the three-dimensional, incompressible Navier–Stokes equation with periodic boundary conditions, if it has a global attractor, then this attractor will be of finite dimensions. ^[13]

See also

• Cycle detection
• Hyperbolic set
• Stable manifold
• Steady state
• Wada basin
• Hidden oscillation
• Rössler attractor
• Stable distribution
• Convergent evolution

References

1. The image and video show the attractor of a second order 3-D Sprott-type polynomial, originally computed by Nicholas Desprez using the Chaoscope freeware (cf. http://www.chaoscope.org/gallery.htm and the linked project files for parameters).
2. Weisstein, Eric W. "Attractor". MathWorld . Retrieved 30 May 2021.
3. Carvalho, A.; Langa, J.A.; Robinson, J. (2012). Attractors for infinite-dimensional non-autonomous dynamical systems. Vol. 182. Springer. p. 109.
4. Kantz, H.; Schreiber, T. (2004). Nonlinear time series analysis. Cambridge university press.
5. John Milnor (1985). "On the concept of attractor". Communications in Mathematical Physics. 99 (2): 177–195. Bibcode:1985CMaPh..99..177M. doi:10.1007/BF01212280. S2CID   120688149.
6. Greenwood, J. A.; J. B. P. Williamson (6 December 1966). "Contact of Nominally Flat Surfaces". Proceedings of the Royal Society. 295 (1442): 300–319. Bibcode:1966RSPSA.295..300G. doi:10.1098/rspa.1966.0242. S2CID   137430238.
7. Vorberger, T. V. (1990). Surface Finish Metrology Tutorial (PDF). U.S. Department of Commerce, National Institute of Standards (NIST). p. 5.
8. Grebogi Celso, Ott Edward, Yorke James A (1987). "Chaos, Strange Attractors, and Fractal Basin Boundaries in Nonlinear Dynamics". Science. 238 (4827): 632–638. Bibcode:1987Sci...238..632G. doi:10.1126/science.238.4827.632. PMID   17816542. S2CID   1586349.
9. Ruelle, David; Takens, Floris (1971). "On the nature of turbulence". Communications in Mathematical Physics. 20 (3): 167–192. Bibcode:1971CMaPh..20..167R. doi:10.1007/bf01646553. S2CID   17074317.
10. Chekroun M. D.; Simonnet E. & Ghil M. (2011). "Stochastic climate dynamics: Random attractors and time-dependent invariant measures". Physica D. 240 (21): 1685–1700. Bibcode:2011PhyD..240.1685C. CiteSeerX   10.1.1.156.5891 . doi:10.1016/j.physd.2011.06.005.
11. Strelioff, C.; Hübler, A. (2006). "Medium-Term Prediction of Chaos". Phys. Rev. Lett. 96 (4): 044101. Bibcode:2006PhRvL..96d4101S. doi:10.1103/PhysRevLett.96.044101. PMID   16486826.
12. Dence, Thomas, "Cubics, chaos and Newton's method", Mathematical Gazette 81, November 1997, 403–408.
13. Geneviève Raugel, Global Attractors in Partial Differential Equations, Handbook of Dynamical Systems, Elsevier, 2002, pp. 885–982.

Further reading

• John Milnor (ed.). "Attractor". Scholarpedia .
• David Ruelle; Floris Takens (1971). "On the nature of turbulence". Communications in Mathematical Physics. 20 (3): 167–192. Bibcode:1971CMaPh..20..167R. doi:10.1007/BF01646553. S2CID   17074317.
• D. Ruelle (1981). "Small random perturbations of dynamical systems and the definition of attractors" (PDF). Communications in Mathematical Physics. 82 (1): 137–151. Bibcode:1981CMaPh..82..137R. doi:10.1007/BF01206949. S2CID   55827557.
• David Ruelle (1989). Elements of Differentiable Dynamics and Bifurcation Theory. Academic Press. ISBN   978-0-12-601710-6.
• Ruelle, David (August 2006). "What is...a Strange Attractor?" (PDF). Notices of the American Mathematical Society . 53 (7): 764–765. Retrieved 16 January 2008.
• Celso Grebogi; Edward Ott; Pelikan; Yorke (1984). "Strange attractors that are not chaotic". Physica D. 13 (1–2): 261–268. Bibcode:1984PhyD...13..261G. doi:10.1016/0167-2789(84)90282-3.
• Chekroun, M. D.; E. Simonnet; M. Ghil (2011). "Stochastic climate dynamics: Random attractors and time-dependent invariant measures". Physica D. 240 (21): 1685–1700. Bibcode:2011PhyD..240.1685C. CiteSeerX   10.1.1.156.5891 . doi:10.1016/j.physd.2011.06.005.
• Edward N. Lorenz (1996) The Essence of Chaos ISBN   0-295-97514-8
• James Gleick (1988) Chaos: Making a New Science ISBN   0-14-009250-1

External links

• Basin of attraction on Scholarpedia
• A gallery of trigonometric strange attractors
• Double scroll attractor Chua's circuit simulation
• A gallery of polynomial strange attractors
• Chaoscope, a 3D Strange Attractor rendering freeware
• Research abstract and software laboratory Archived 28 June 2022 at the Wayback Machine
• Online strange attractors generator
• Interactive trigonometric attractors generator
• Economic attractor