Henon Map: Initial Conditions & Attractor Shape

The Hénon map $(x, y) \mapsto (1 - a x^{2} + y, b x)$ is a two-dimensional invertible quadratic diffeomorphism whose classical parameters $a = 1.4, b = 0.3$ produce a strange attractor of fractal dimension $\approx 1.26$ . Key result: for almost every initial condition inside the basin of attraction, orbits converge to a single, geometrically fixed attractor whose shape is independent of the starting point; only points in the escaping set diverge to infinity. Applies to dissipative chaos, dimension theory, and dynamical systems pedagogy. Introduced by Michel Hénon in 1976.

The Hénon map is the canonical example of a low-dimensional dissipative dynamical system exhibiting a strange attractor. Defined on $R^{2}$ by a quadratic nonlinearity and a constant Jacobian determinant $- b$ , it compresses areas by a factor $∣ b ∣$ each iteration while stretching and folding the plane, the mechanism that generates chaos. Its importance lies in being the simplest invertible map with a verified chaotic attractor — a discrete-time analog of the Lorenz flow — making it a testbed for fractal dimension, Lyapunov exponents, and rigorous proofs of chaos. The surprising fact: although different initial conditions follow wildly divergent individual trajectories (sensitive dependence), the set they fill in the limit is the same attractor with the same shape, so the attractor's geometry is a global invariant insensitive to initial conditions, whereas the orbit on it is exquisitely sensitive.

The Core Insight

There is a sharp dichotomy between the orbit and the attractor. The choice of initial condition $(x_{0}, y_{0})$ determines which infinite sequence of points you visit and how fast you approach the limit set — but provided $(x_{0}, y_{0})$ lies in the basin of attraction, it does not determine the shape, location, or fractal dimension of the limiting attractor. Initial conditions partition the plane into the basin (orbits captured by the attractor) and the escaping set (orbits diverging to $\infty$ ); the attractor shape is fixed by the parameters $(a, b)$ alone.

Formal Statement

The Hénon Map

The Hénon map is the diffeomorphism $H_{a, b} : R^{2} \to R^{2}$ given by

H_{a, b} (x, y) = (1 - a x^{2} + y, b x)

⟳This is the standard Hénon (1976) form; the alternative

x_{n + 1} = y_{n} + 1 - a x_{n}^{2}, y_{n + 1} = b x_{n}

is identical componentwise.

Here $a > 0$ controls the strength of the nonlinear (folding) term and $b$ with $0 < ∣ b ∣ < 1$ controls dissipation. The orbit of an initial condition $(x_{0}, y_{0})$ is the sequence $(x_{n}, y_{n}) = H_{a, b}^{n} (x_{0}, y_{0})$ , i.e.

x_{n + 1} = 1 - a x_{n}^{2} + y_{n}, y_{n + 1} = b x_{n} .

⟳Apply [1] componentwise;

y_{n + 1} = b x_{n}

is the linear contraction,

x_{n + 1}

carries the quadratic nonlinearity.

The Jacobian matrix is

D H = (- 2 a x b 10), det D H = - b .

⟳Partial-differentiate the components of [1]; the off-diagonal

b

and

1

give

det = (- 2 a x) (0) - (1) (b) = - b

The constant determinant $- b$ means every iteration multiplies oriented area by $- b$ : contraction (since $∣ b ∣ < 1$ ) plus an orientation reversal.

Inverse and Invertibility

Because $det D H = - b \neq = 0$ for $b \neq = 0$ , $H_{a, b}$ is a diffeomorphism with explicit inverse

H_{a, b}^{- 1} (x, y) = (\frac{y}{b}, x - 1 + a \frac{y ^{2}}{b ^{2}}) .

⟳From

y = b x

solve

x = y / b

; substitute into

x_{old new} = 1 - a x^{2} + y_{old}

and solve for the preimage's

y

-coordinate.

Invertibility distinguishes Hénon from the (non-invertible) logistic map and makes it a genuine 2D analog of a Poincaré return map.

Existence of the Attractor — Benedicks–Carleson

For a positive-measure set of parameters near $b = 0$ and $a$ close to $2$ , the Hénon map possesses a topologically transitive attractor $Λ$ supporting an SRB (Sinai–Ruelle–Bowen) measure, with a positive Lyapunov exponent along the unstable direction. Consequently orbits exhibit sensitive dependence on initial conditions while statistically converging to a unique invariant measure on $Λ$ .

Basin and Escape Dichotomy

For the classical parameters $a = 1.4, b = 0.3$ there exists a bounded trapping region $R$ (a quadrilateral) with $H (R) \subset int R$ . Every initial condition in the basin $B = ⋃_{n \geq 0} H^{- n} (R)$ has $ω$ -limit set contained in the attractor $Λ = ⋂_{n \geq 0} H^{n} (R)$ ; every initial condition outside $\overline{B}$ escapes to infinity. The attractor shape $Λ$ is independent of the chosen initial condition within $B$ .

Why It Holds / Derivation

Why the Attractor Shape Is Initial-Condition Independent

Step 1 — Trapping. Construct a quadrilateral

R

whose image

H (R)

is a horseshoe-like bent strip lying strictly inside

R

. Existence is verified by checking the four corner images numerically; this gives

H (R) \subset int R

. Step 2 — Nested intersection. Define

R_{n} = H^{n} (R)

. Trapping implies

R_{n + 1} \subset R_{n}

. The attractor is the decreasing intersection

Λ = n \geq 0 ⋂ H^{n} (R) .

⟳A nested family of nonempty compact sets has nonempty compact intersection;

H

-invariance follows since

H (Λ) = Λ

Step 3 — Area collapse. By [3],

area (R_{n}) = ∣ b ∣^{n} area (R) \to 0

. Thus

Λ

has zero Lebesgue measure: it is a thin (Cantor-like) set transverse to the unstable foliation. Step 4 — Attraction. Any

(x_{0}, y_{0}) \in B

eventually enters

R

and thereafter stays; its iterates lie in

R_{n}

for ever larger

n

, so

dist (H^{n} (x_{0}, y_{0}), Λ) \to 0

. The limit set is contained in

Λ

for every such start. Step 5 — Independence.

Λ

in [5] is defined purely through

H_{a, b}

and

R

; no initial condition appears. Different starts merely yield different orbits dense in (essentially all of, for transitive

Λ

) the same fixed set.

■

Fixed Points and Their Stability

Fixed points solve $x = 1 - a x^{2} + b x, y = b x$ , giving

x_{\pm} = \frac{( b - 1 ) \pm ( 1 - b ) ^{2} + 4 a}{2 a}, y_{\pm} = b x_{\pm} .

⟳Set

x_{n + 1} = x_{n} = x

y_{n + 1} = y_{n} = y = b x

; substitute to get

a x^{2} + (1 - b) x - 1 = 0

and apply the quadratic formula.

Stability follows from the eigenvalues $λ$ of $D H$ at $(x_{\pm}, y_{\pm})$ , roots of

λ^{2} + 2 a x_{\pm} λ - b = 0.

⟳Characteristic polynomial

det (D H - λ I) = λ^{2} - (tr) λ + det = λ^{2} + 2 a x λ - b

using

tr = - 2 a x, det = - b

For $a = 1.4, b = 0.3$ , $(x_{+}, y_{+}) \approx (0.6314, 0.1894)$ is a saddle whose unstable manifold's closure is the visible attractor; $(x_{-}, y_{-})$ is also a saddle lying off the attractor. The attractor is thus the closure of the unstable manifold of $(x_{+}, y_{+})$ — again a construction with no reference to initial conditions.

Michel Hénon | 1931–2013

A French astronomer and mathematician at the Nice Observatory, Hénon worked on stellar dynamics and the restricted three-body problem. While studying the Lorenz system he sought the simplest discrete map that would reproduce its stretching-and-folding chaos. By collapsing the Lorenz flow's Poincaré section into a two-line recurrence, he produced his 1976 map — deliberately minimal, analytically tractable, and explicitly invertible. His earlier work on the "Hénon–Heiles" Hamiltonian (1964) already revealed his taste for simple models exposing deep chaotic structure.

Building Intuition

Picture the trapping quadrilateral as a sheet of dough. Each iteration does three things encoded in [1]: the quadratic term $- a x^{2}$ folds the sheet over (a parabola-shaped bend), the $+ y$ term shears, and the $y^{'} = b x$ rule squashes it to $0.3$ of its width. Fold–shear–squash, repeated forever, lays down infinitely many parallel sheets nested at every scale — the Cantor cross-section that gives fractal dimension $\approx 1.26$ . Now drop two nearby drops of ink (two nearby initial conditions). They both get pulled onto the dough's accumulating layers (same shape!) but their individual paths diverge exponentially because the folding repeatedly separates neighbors. That is the whole story: same canvas, unpredictable brushstrokes.

Hénon attractor: orbit converges to a shape fixed by (a,b), not by the start

Worked Examples

Locating the Fixed Point and Verifying the Saddle

Given:

a = 1.4, b = 0.3

. Find the "physical" fixed point

(x_{+}, y_{+})

and confirm it is a saddle. Step 1 — Apply [6]. Discriminant

= (1 - b)^{2} + 4 a = (0.7)^{2} + 5.6 = 0.49 + 5.6 = 6.09

6.09 = 2.4678

x_{+} = \frac{( 0.3 - 1 ) + 2.4678}{2 ( 1.4 )} = \frac{- 0.7 + 2.4678}{2.8} = \frac{1.7678}{2.8} = 0.6313.

y_{+} = b x_{+} = 0.3 (0.6313) = 0.1894.

Step 2 — Eigenvalues via [7].

2 a x_{+} = 2 (1.4) (0.6313) = 1.7676

. Solve

λ^{2} + 1.7676 λ - 0.3 = 0

λ = \frac{- 1.7676 \pm 1.767 6 ^{2} + 1.2}{2} = \frac{- 1.7676 \pm 4.3244}{2} = \frac{- 1.7676 \pm 2.0795}{2} .

λ_{u} = 0.1560, λ_{s} = - 1.9236.

Step 3 — Classify. Wait: check magnitudes —

∣ λ_{s} ∣ = 1.9236 > 1

and

∣ λ_{u} ∣ = 0.156 < 1

, so one expanding, one contracting direction.

(x_{+}, y_{+}) = (0.6313, 0.1894), λ \approx {- 1.924, 0.156} \Rightarrow saddle

Two Initial Conditions, Same Limit Set

Given: Iterate [2] from

(x_{0}, y_{0}) = (0, 0)

and from

(0.1, 0.1)

a = 1.4, b = 0.3

. Step 1 — Start A $(0, 0)$ :

(x_{1}, y_{1}) = (1 - 0 + 0, 0) = (1, 0)

;

(x_{2}, y_{2}) = (1 - 1.4 + 0, 0.3) = (- 0.4, 0.3)

;

(x_{3}, y_{3}) = (1 - 1.4 (0.16) + 0.3, - 0.12) = (1.076, - 0.12)

. Step 2 — Start B $(0.1, 0.1)$ :

(x_{1}, y_{1}) = (1 - 1.4 (0.01) + 0.1, 0.03) = (1.086, 0.03)

;

(x_{2}, y_{2}) = (1 - 1.4 (1.179) + 0.03, 0.3258) = (- 0.621, 0.326)

. Step 3 — Observe divergence then capture. After

\sim

5 steps the two orbits are macroscopically different (sensitive dependence), yet after discarding

\sim

1000 transient iterations, the point clouds coincide — both fill the same crescent attractor bounded roughly by

x \in [- 1.28, 1.27], y \in [- 0.38, 0.38]

Different transients \Rightarrow identical limiting attractor shape; D_{0} \approx 1.26

Common Misconceptions

"Changing the initial condition changes the attractor." False inside the basin. The attractor's geometry, box-counting dimension (

\approx 1.26

), and SRB measure are functions of

(a, b)

only. What the initial condition changes is the transient and the specific point sequence. The one case where the start matters qualitatively is whether it lies inside the basin at all — start outside and the orbit diverges to

\infty

(no attractor reached). "Sensitive dependence means the long-term statistics are unpredictable." The pointwise future is unpredictable, but the statistical future is not: histograms of long orbits converge to the SRB measure regardless of initial condition (for a.e. start). Predictability of averages coexists with unpredictability of trajectories.

An Initial Condition That Does NOT Reach the Attractor

Take $(x_{0}, y_{0}) = (10, 10)$ with $a = 1.4, b = 0.3$ . Then $x_{1} = 1 - 1.4 (100) + 10 = - 129$ , $x_{2} = 1 - 1.4 (12 9^{2}) + 3 = - 23297$ , and $∣ x_{n} ∣ \to \infty$ super-exponentially. This orbit never approaches $Λ$ — it lies in the escaping set. Hence the claim "every initial condition converges to the Hénon attractor" is false; the correct statement restricts to the basin $B$ .

Periodic Windows Hide the Strange Attractor

At $a = 1.4, b = 0.3$ the attractor is strange, but at, e.g., $a = 1.0, b = 0.3$ the attractor is a stable period-7 cycle — a finite point set, dimension $0$ . Picking initial conditions there shows orbits converging to seven points, not a fractal. Thus "Hénon $\Rightarrow$ strange attractor" depends entirely on parameters, not on initial conditions.

Connections

The Hénon map is the discrete sibling of the Lorenz System: Hénon distilled the Lorenz Poincaré-section folding into two lines. Its fold-and-stretch mechanism is a smooth realization of the Smale Horseshoe, which provides the symbolic dynamics underlying chaos. The cross-sectional Cantor structure ties it to Fractal Dimension and the Box-Counting Dimension; its exponential separation of orbits is quantified by Lyapunov Exponents. The statistical convergence to a fixed measure is governed by SRB Measures, and the basin/escape split exemplifies Basins of Attraction. The 1D shadow of Hénon as $b \to 0$ is the Logistic Map.

Real-World Applications

Plasma & particle accelerators — Hénon-type maps model the transverse phase-space dynamics of charged-particle beams; the trapping region predicts the dynamic aperture (region of stable, non-escaping initial conditions).
Secure communications / chaos encryption — the attractor's sensitive dependence and reproducible SRB statistics are exploited in pseudo-random keystream generators based on Hénon iteration.
Stellar dynamics — Hénon's original motivation; the map approximates Poincaré sections of orbits in galactic potentials, distinguishing regular from chaotic stellar trajectories.
Time-series benchmarking — Hénon data is the standard test set for delay-coordinate embedding, attractor reconstruction (Takens), and nonlinear prediction algorithms.
Climate/toy attractors — used pedagogically to validate SRB-measure estimators and ergodic-average convergence before deployment on high-dimensional geophysical models.
Neural and circuit oscillators — discrete maps of Hénon form arise as Poincaré returns of forced nonlinear oscillators (Chua circuits, FitzHugh–Nagumo discretizations).
Image/data compression research — fractal-dimension measurements on Hénon attractors calibrate multifractal-spectrum software.

Historical Development

1963 — Edward Lorenz publishes "Deterministic Nonperiodic Flow," exhibiting the first strange attractor in a flow and motivating discrete analogs.
1964 — Hénon & Heiles study a 2-DOF Hamiltonian, revealing chaotic Poincaré sections in celestial mechanics.
1967 — Stephen Smale formalizes the horseshoe, the topological skeleton of folding chaos.
1976 — Michel Hénon publishes "A two-dimensional mapping with a strange attractor" (Commun. Math. Phys.), introducing [1] with $a = 1.4, b = 0.3$ .
1980 — Numerical estimates pin the attractor's fractal dimension near $1.26$ and Lyapunov exponent $λ_{1} \approx 0.42$ .
1985 — Benedicks–Carleson begin the program proving abundance of chaos in 1D unimodal maps, the technical precursor to the 2D case.
1991 — Benedicks & Carleson prove, for a positive-measure parameter set near $b = 0$ , that the Hénon map has a genuine strange attractor with positive Lyapunov exponent.
1993 — Benedicks & Young construct the SRB measure for these Hénon attractors, establishing statistical regularity.
2000s — Rigorous computer-assisted dimension and topological-entropy bounds (interval arithmetic) refine the classical-parameter picture.

Subtleties & Limits

At the classical parameters $a = 1.4, b = 0.3$ it is not rigorously proven that the observed object is a strange attractor in the Benedicks–Carleson sense — those theorems hold for small $b$ near a positive-measure parameter set, which may not include $(1.4, 0.3)$ . The set could in principle contain a long-period stable sink with a microscopically thin basin. Open questions: (i) Determine the exact Hausdorff (vs. box-counting) dimension of the classical attractor. (ii) Characterize the full measure-theoretic structure of the basin boundary, which appears fractal. (iii) Map the precise parameter regions where the attractor is strange versus periodic — the "Hénon parameter zoo" is only partially understood, with intricate interleaving of chaotic and periodic windows.