chaos · dynamical systems · universality

Bifurcation Diagram & Feigenbaum Constants

The logistic map xn+1 = r·xn(1−xn) transitions from order to chaos through an infinite cascade of period-doublings, governed by the universal constant δ ≈ 4.6692.

Preset
Colour

The Logistic Map

The logistic map is the simplest one-dimensional dynamical system that exhibits chaos. Originally introduced by ecologist Robert May (1976) to model population dynamics, it became the paradigm for the period-doubling route to chaos:

x_{n+1} = r · x_n · (1 − x_n) x_n ∈ [0,1], r ∈ [0,4]

Despite its seeming simplicity — a single quadratic equation — the long-run behaviour depends sensitively on r and fills a bifurcation diagram of remarkable complexity.

Period-Doubling Cascade

As r increases from 0 to 4, the logistic map undergoes successive bifurcations where a stable periodic orbit loses stability and splits into an orbit of twice the period:

PeriodBifurcation at r ≈Width ΔrRatio Δrₙ/Δrₙ₊₁
1 → 23.000 0000.449 490
2 → 43.449 4900.096 2574.669
4 → 83.544 0900.020 6164.669
8 → 163.564 4070.004 4134.669
16 → 323.568 7590.000 9454.669
∞ (chaos onset)3.569 945 7…0δ ≈ 4.6692…

The ratio of successive interval widths converges to Feigenbaum's constant δ = 4.6692040890097185…

Feigenbaum's Constants

δ = lim_{n→∞} (r_n − r_{n−1}) / (r_{n+1} − r_n) ≈ 4.6692040890097185... α = lim_{n→∞} d_n / d_{n+1} ≈ 2.502907875095892... (vertical scaling)

In 1978 Mitchell Feigenbaum discovered that both constants are universal: they appear in any smooth one-dimensional map with a quadratic maximum — sin(πx), x(1−x²), and countless others all exhibit the same ratio. This universality is explained by a fixed-point equation in the space of maps, analogous to the renormalisation group in statistical physics.

🔮

Universality

δ and α appear in every smooth unimodal map. Period-doubling cascades have been observed in fluid turbulence, electronic circuits, and optical lasers.

🌀

Self-Similarity

Zoom into any chaotic band and you see miniature copies of the whole diagram, scaled by δ horizontally and α vertically — a fractal structure.

🪟

Periodic Windows

Within the chaotic regime, narrow windows of stable periodicity appear (period 3, 5, 7, …) before bifurcating into their own period-doubling cascades.

📐

Sharkovskii Ordering

Period 3 implies all periods: 3 ▷ 5 ▷ 7 ▷ … ▷ 6 ▷ 10 ▷ … ▷ 2 ▷ 4 ▷ 8 ▷ … The period-3 window around r≈3.831 is the most prominent.

The Logistic Map in Detail

Fixed Points

x* = r·x*(1−x*) gives fixed points x* = 0 and x* = 1 − 1/r. The non-trivial fixed point is stable (|df/dx|<1) for r ∈ (1, 3). At r=3 the stability eigenvalue equals −1, triggering the first bifurcation.

Lyapunov Exponent

λ = lim_{N→∞} (1/N) Σ ln|r(1 − 2x_n)|

λ < 0 indicates a stable periodic orbit; λ > 0 indicates chaos. At the period-doubling accumulation point λ = 0. In the chaotic regime λ fluctuates around +0.69.

Chaos Onset r∞ ≈ 3.5699456718…

The accumulation point of period doublings is r∞ ≈ 3.5699456718789691. Beyond this value, for most r, the orbit densely fills one or more intervals (chaotic bands). The number of bands halves at each reverse-bifurcation as r increases further.

Zooming the Diagram

The simulation supports three zoom methods:

Moving your mouse over the diagram displays the live orbit (time series x_n vs n) for the r value under the cursor in the strip below the diagram.

Summary Table

QuantityValue / FormulaSignificance
Feigenbaum δ4.6692040890097185…Horizontal period-doubling ratio
Feigenbaum α2.5029078750958928…Vertical attractor scaling
First bifurcationr₁ ≈ 3.000Period 1→2
Chaos onset r∞3.5699456718…Accumulation of period doublings
Period-3 windowr ≈ 3.8284 – 3.8571Sharkovskii: period 3 ⟹ all periods
Lyapunov λ (chaotic)≈ +0.69Exponential sensitivity to initial conditions
Hausdorff dimension≈ 0.538Fractal dimension of attractor at r∞
Universality classAll unimodal mapsδ, α identical for all quadratic maxima

Frequently Asked Questions

What is a bifurcation diagram?
A bifurcation diagram plots the long-run attractors of a dynamical system versus a control parameter. For the logistic map x_{n+1}=r·x_n(1−x_n), each vertical slice at a fixed r shows the set of x-values the orbit eventually settles on (after discarding a transient). A single point indicates a fixed point; two points indicate a period-2 cycle; a filled region indicates chaos.
What is Feigenbaum's constant δ ≈ 4.6692?
δ = lim (r_n − r_{n−1})/(r_{n+1} − r_n) ≈ 4.6692040890…, where r_n is the r-value of the nth period-doubling bifurcation. It measures how rapidly the bifurcations accelerate as the chaos onset is approached. Remarkably, this constant is universal — it is the same for any smooth one-dimensional map with a single quadratic maximum. Feigenbaum derived it from a functional fixed-point equation in 1978.
What is the period-3 window?
Around r ≈ 3.8284 a stable period-3 cycle appears within the chaotic regime. By Sharkovskii's theorem (1964), a dynamical system with a period-3 orbit must also have orbits of every other period. Li and Yorke's 1975 paper "Period Three Implies Chaos" coined the word "chaos" in mathematics and showed that a period-3 orbit implies the existence of an uncountable set of initial conditions with chaotic (non-periodic) behaviour.
Why does this simulation matter beyond the logistic map?
The period-doubling route to chaos has been observed experimentally in fluid dynamics (Ruelle-Takens scenario), Josephson junctions, heart arrhythmias, and electronic oscillators. The Feigenbaum constants appear in each case. This universality — predicted by renormalisation group theory — was one of the first hints that chaos is a fundamental and broadly applicable phenomenon in nature.

Explore More Chaos & Dynamical Systems

Dive deeper into chaotic attractors and emergent complexity.