Feigenbaum Constant

DOWNLOAD Mathematica Notebook Contribute to this entry FeigenbaumConstantBifurcation

The Feigenbaum constant delta is a universal constant for functions approaching chaos via period doubling. It was discovered by Feigenbaum in 1975 (Feigenbaum 1979) while studying the fixed points of the iterated function

 f(x)=1-mu|x|^r,
(1)

and characterizes the geometric approach of the bifurcation parameter to its limiting value as the parameter mu is increased for fixed x. The plot above is made by iterating equation (1) with r=2 several hundred times for a series of discrete but closely spaced values of mu, discarding the first hundred or so points before the iteration has settled down to its fixed points, and then plotting the points remaining.

FeigenbaumConstantIteration

A similar plot that more directly shows the cycle may be constructed by plotting f^n(x)-x as a function of mu. The plot above (Trott, pers. comm.) shows the resulting curves for n=1, 2, and 4.

Let mu_n be the point at which a period 2^n-cycle appears, and denote the converged value by mu_infty. Assuming geometric convergence, the difference between this value and mu_n is denoted

 lim_(n->infty)mu_infty-mu_n=Gamma/(delta^n),
(2)

where Gamma is a constant and delta>1 is a constant now known as the Feigenbaum constant. Solving for delta gives

 delta=lim_(n->infty)(mu_(n+1)-mu_n)/(mu_(n+2)-mu_(n+1))
(3)

(Rasband 1990, p. 23; Briggs 1991). An additional constant alpha, defined as the separation of adjacent elements of period doubled attractors from one double to the next, has a value

 alpha=lim_(n->infty)(d_n)/(d_(n+1)),
(4)

where d_n is the value of the nearest cycle element to 0 in the 2^n cycle (Rasband 1990, p. 37; Briggs 1991).

For equation (1) with r=2, the onsets of bifurcations occur at mu=0.75, 1.25, 1.368099, 1.39405, 1.399631, ..., giving convergents to delta for n=1, 2, 3, ... of 4.23374, 4.5515, 4.64617, ....

For the logistic map,

delta=4.669201609102990...
(5)
Gamma=2.637...
(6)
mu_infty=3.569945672...
(7)
alpha=-2.502907875....
(8)

(OEIS A006890, A098587, and A006891; Broadhurst 1999; Wolfram 2002, p. 920), where delta is known as the Feigenbaum constant and alpha is the associated "reduction parameter."

Briggs (1991) calculated delta to 84 digits, Briggs (1997) to 576 decimal places (of which 344 were correct), and Broadhurst (1999) to 1018 decimal places. It is not known if the Feigenbaum constant delta is algebraic, or if it can be expressed in terms of other mathematical constants (Borwein and Bailey 2003, p. 53).

Briggs (1991) calculated alpha to 107 digits, Briggs (1997) to 576 decimal places (of which 346 were correct), and Broadhurst (1999) to 1018 decimal places.

Amazingly, the Feigenbaum constant delta and associated reduction parameter alpha are "universal" for all one-dimensional maps f(x) if f(x) has a single locally quadratic maximum. This was conjecture by Feigenbaum, and demonstrated rigorously by Lanford (1982) for the case r=2, and by Epstein (1985) for all r<14.

More specifically, the Feigenbaum constant is universal for one-dimensional maps if the Schwarzian derivative

 D_(Schwarzian)=(f^(''')(x))/(f^'(x))-3/2[(f^('')(x))/(f^'(x))]^2
(9)

is negative in the bounded interval (Tabor 1989, p. 220). Examples of maps which are universal include the Hénon map, logistic map, Lorenz attractor, Navier-Stokes truncations, and sine map x_(n+1)=asin(pix_n). The value of the Feigenbaum constant can be computed explicitly using functional group renormalization theory. The universal constant also occurs in phase transitions in physics.

The value of alpha for a universal map may be approximated from functional group renormalization theory to the zeroth order by solving

 1-alpha^(-1)=(1-alpha^(-2))/([1-alpha^(-2)(1-alpha^(-1))]^2),
(10)

which can be rewritten as the quintic equation

 alpha^5+2alpha^4-2alpha^3-alpha^2+2alpha-1=0.
(11)

Solving numerically for the smallest real root gives alpha=-2.48634..., only 0.7% off from the actual value (Feigenbaum 1988).

For an area-preserving two-dimensional map with

x_(n+1)=f(x_n,y_n)
(12)
y_(n+1)=g(x_n,y_n),
(13)

the Feigenbaum constant is delta=8.7210978... (Tabor 1989, p. 225).

For a function of the form (1), the Feigenbaum constant for various r is given in the following table (Briggs 1991, Briggs et al. 1991, Finch 2003), which updates the values in Tabor (1989, p. 225).

rdeltaalpha
35.9679687038...1.9276909638...
47.2846862171...1.6903029714...
58.3494991320...1.5557712501...
69.2962468327...1.4677424503...

Broadhurst (1999) considered additional Feigenbaum constants. Let g(x) and f(x) be even functions with g(0)=f(0)=1 and

(g(alphax))/alpha=g(g(x))
(14)
(deltaf(alphax))/alpha=g^'(g(x))f(x)+f(g(x))
(15)

and delta as large as possible. Let (b,c,d) be positive numbers with

 g(b)=0=1/(g(c+id))
(16)

and (b,c^2+d^2) as small as possible. Also let kappa be the order of the nearest singularity, with

 1/(g(c+id+z))=O(z^kappa)
(17)

as z tends to zero. The values of these constants are summarized in the following table.

constantSloanevalue
bA1192770.83236723690531642484...
cA1192781.8312589849371314853...
dA1192792.6831509004740718014...
kappaA1192801.3554618047064087438...

Wolfram Web Resources

Mathematica »

The #1 tool for creating Demonstrations and anything technical.

Wolfram|Alpha »

Explore anything with the first computational knowledge engine.

Wolfram Demonstrations Project »

Explore thousands of free applications across science, mathematics, engineering, technology, business, art, finance, social sciences, and more.

Computerbasedmath.org »

Join the initiative for modernizing math education.

Online Integral Calculator »

Solve integrals with Wolfram|Alpha.

Step-by-step Solutions »

Walk through homework problems step-by-step from beginning to end. Hints help you try the next step on your own.

Wolfram Problem Generator »

Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.

Wolfram Education Portal »

Collection of teaching and learning tools built by Wolfram education experts: dynamic textbook, lesson plans, widgets, interactive Demonstrations, and more.

Wolfram Language »

Knowledge-based programming for everyone.