Poincaré Sections and Return Maps

Every article in this chapter has quietly leaned on one instrument. The driven pendulum's chaos became visible as a stroboscopic dot-cloud; the double pendulum's KAM transition was read off dot patterns at θ1=0\theta_1 = 0; the period-doubling cascade was defined through multipliers of a "return map" whose construction was deferred. This article pays the debt: it builds the Poincaré section properly — the single most important diagnostic instrument in dynamics — and proves the facts the chapter has been borrowing. The idea is Poincaré's characteristic economy. A trajectory of a flow is a tangled curve that is hard to interrogate; so stop watching the curve. Stretch a surface across the flow, transverse to it, and record only the sequence of points where the trajectory pokes through. Those piercings define a map — the first-return map PP — and the map knows everything qualitative about the flow while living in one dimension fewer. Periodic orbits become fixed points; stability becomes eigenvalues of DPDP (the Floquet multipliers, with the useless along-flow direction already quotiented away); quasiperiodic tori become smooth invariant curves; chaos becomes structured dust. The entire fixed-point technology of the early M10 articles — Jacobians, classification, bifurcations — re-enters, now applied to PP, and the three ways a multiplier can leave the unit circle reproduce exactly the bifurcation dictionary of the cascade articles. The construction also carries the chapter's deepest historical cargo. It was on a section of the three-body problem, in 1889, that Poincaré followed the stable and unstable manifolds of a saddle point, realized they must cross — and cross again, infinitely often, in a figure he declared he would not even attempt to draw. That homoclinic tangle, discovered while correcting a prize-winning error, is the first sighting of chaos in mathematics, and the section is the instrument that made it visible.

Key takeaways
  • A Poincaré section is a codimension-one surface Σ\Sigma transverse to the flow (∇g⋅f≠0\nabla g \cdot \mathbf{f} \neq 0); successive same-direction crossings define the first-return map P:Σ→ΣP: \Sigma \to \Sigma, an (n−1)(n-1)-dimensional discrete system.
  • Fixed points of PP are periodic orbits of the flow; the eigenvalues of DPDP are the Floquet multipliers, with the trivial along-flow multiplier 11 automatically removed.
  • Stability: all ∣μi∣<1|\mu_i| < 1; the three unit-circle exits are the bifurcation dictionary — +1+1 saddle-node, −1-1 period-doubling, complex pair Neimark–Sacker (torus birth).
  • Area bookkeeping transfers exactly: det⁡DP=exp⁡∮∇⋅f dt\det DP = \exp\oint \nabla\cdot\mathbf{f}\,\mathrm{d}t — equal to 1 for Hamiltonian flows (area-preserving sections, KAM structure) and less than 1 in modulus for dissipative ones (attractors possible).
  • Reading chart: finite dots = periodic; dots filling a smooth closed curve = quasiperiodic torus; structured fractal dust = chaos.
  • For driven systems the natural section is the stroboscopic map — sample once per drive period — with constant return time.
  • Key idea: sample a flow only where it crosses a transverse surface: the resulting first-return map PP has the same periodic orbits, stability, and bifurcations as the flow, but lives in one fewer dimension, so periodicity, quasiperiodicity, and chaos become visually and quantitatively distinct.

Where this lives today

The Poincaré section is arguably the most heavily industrialized object in this chapter. In accelerator physics it is not an analysis tool but the native description: a stored particle's turn-by-turn coordinates at a fixed azimuth of the ring are literally iterates of a return map, and storage-ring design is the engineering of that map — tune (rotation number) placement away from low-order rationals, dynamic-aperture surveys that are section portraits, island chains and chaotic layers appearing exactly per Derivation C. In magnetic-confinement fusion, field-line Poincaré plots are the daily diagnostic of stellarator and tokamak design: nested invariant curves are good magnetic surfaces, island chains are resonant field errors, and stochastic (dust-filled) regions leak heat — a KAM textbook rendered in plasma. Celestial mechanics maps stability of exoplanetary systems and designs low-energy spacecraft trajectories through the tangle-generated transport channels of the three-body problem — the "interplanetary superhighway" is a homoclinic web used as infrastructure. Within this course, the section is the connective instrument of the whole chaos block: it defined the maps whose multipliers ran the cascade, displayed the driven pendulum's strange attractor and the double pendulum's KAM transition, and carries Hamiltonian chaos (via area-preserving maps like the standard map) to the KAM and Chirikov material that closes the track. One article remains: the chapter's simplest object — a one-dimensional map studied on its own terms — where everything the section technology reveals can be computed to arbitrary depth.

The physical picture

A long-exposure photograph of an orbit is a blur; the section is the fix. Stretch a pane of glass across phase space where the flow punches through it cleanly, and mark a dot at every piercing (counting only one crossing direction, so each lap contributes once). The blur becomes a ledger of returns: the section discards the one boring fact — that the trajectory keeps moving forward — and keeps everything about whether, where, and how it comes back. In three dimensions, the pane is two-dimensional and the flow has become a planar map: the entire arsenal built for planar systems earlier in this chapter applies again, one level up.

Figure
The Poincaré section — stop watching the curve, watch only where it comes back through the glass
section Σ (transverse) flow trajectory record only the piercings: x(k+1) = P(x(k)) an n-dimensional flow becomes an (n−1)-dimensional map — the along-flow direction is quotiented out, nothing else is lost

The dots organize into a reading chart that the whole chapter has been using informally. A periodic orbit pierces the pane at the same spot every lap: finitely many dots, forever — one for period-TT, kk for period-kTkT, and the doubling cascade is dots splitting in two. A quasiperiodic orbit winds on a torus; the torus meets the pane in a circle, and the dots hop around that circle by a fixed angle each return, never repeating (irrational rotation) but filling the curve densely — the nested rings of the double pendulum's low-energy sections. Chaos scatters dots into structured dust: fractal layering if the system is dissipative (the driven pendulum's attractor cross-section), a fat chaotic sea with island chains if it is Hamiltonian. One glance at a section and the dynamical taxonomy is read off.

Figure
Reading a section at a glance — dots for periodicity, curves for tori, structured dust for chaos
periodic orbit — finite dots period-3: 3 dots quasiperiodic — invariant curve torus section: dots fill a curve Hamiltonian mix — islands + sea three orbits of one return map

Why is the map sharper than the flow it samples? Because of what it silently removes. Perturbations along the orbit direction always carry the marginal multiplier 11 — they merely shift the clock — and in the flow's monodromy matrix this trivial eigenvalue sits in the way of every stability statement. The section quotients it out by construction: deviations within Σ\Sigma are exactly the deviations that matter, and the eigenvalues of DPDP are exactly the nontrivial Floquet multipliers. Better still, the return map converts the continuous bifurcation zoo into the discrete dictionary already mastered in the cascade articles: a multiplier through +1+1 collides orbits (saddle-node), through −1-1 doubles the period (flip), and a complex pair exiting the unit circle births an invariant curve — a torus of the flow — in the Neimark–Sacker bifurcation. And the section is where chaos was discovered. On it, a saddle fixed point carries one-dimensional stable and unstable manifolds — curves, drawable. Poincaré, mapping the restricted three-body problem, asked how the unstable manifold returns and found it must cross the stable one transversally; but the image of a crossing point under PP is another crossing point, so one intersection begets infinitely many, and the unstable manifold must fold back and forth across the stable one in ever-wilder oscillations without ever self-intersecting. The resulting homoclinic tangle forces stretch-and-fold dynamics — Smale later showed a horseshoe map, with all its chaos, lives inside every transverse tangle. The figure below is the picture Poincaré described and declined to draw.

Figure
The homoclinic tangle on the section — stable and unstable manifolds of one saddle crossing infinitely often
saddle point of P stable manifold W s unstable manifold W u — folded across W s infinitely often each amber crossing is a homoclinic point; one crossing forces infinitely many — the tangle Poincaré refused to draw

The instrument, formalized

The setup: an autonomous flow and a well-placed pane of glass,

x˙=f(x),Σ={x:g(x)=0, ∇g⋅f≠0}(1)\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x}), \qquad \Sigma = \{\mathbf{x} : g(\mathbf{x}) = 0,\ \nabla g\cdot\mathbf{f} \neq 0\} \qquad (1)

where: gg — any smooth scalar whose zero set is the section surface; the transversality condition ∇g⋅f≠0\nabla g \cdot \mathbf{f} \neq 0 says the flow crosses Σ\Sigma at a nonzero angle, never grazing it. Transversality is the entire licence for what follows: it guarantees (Derivation A) a smooth, locally unique first-return time τ(x)\tau(\mathbf{x}), and with it the Poincaré map

xk+1=P(xk)=ϕτ(xk)(xk)(2)\mathbf{x}_{k+1} = P(\mathbf{x}_k) = \boldsymbol{\phi}_{\tau(\mathbf{x}_k)}(\mathbf{x}_k) \qquad (2)

where: ϕt\boldsymbol{\phi}_t — the flow map (solution operator); PP maps the (n−1)(n-1)-dimensional surface to itself. Note that the return time varies from point to point — the stroboscopic map of driven systems, P=ϕTP = \boldsymbol{\phi}_T with fixed T=2π/ωdriveT = 2\pi/\omega_{\text{drive}}, is the special case where the drive phase supplies a section with constant return time. The dictionary begins with periodic orbits:

P(x∗)=x∗  ⟺  the flow has a periodic orbit through x∗ of period T=τ(x∗)(3)P(\mathbf{x}^*) = \mathbf{x}^* \iff \text{the flow has a periodic orbit through } \mathbf{x}^* \text{ of period } T = \tau(\mathbf{x}^*) \qquad (3)

where: a kk-cycle of PP likewise corresponds to a closed orbit crossing the pane kk times. Stability of the orbit is stability of the fixed point: linearize the map,

ξk+1=DP(x∗) ξk,stable  ⟺  ∣μi∣<1 for all eigenvalues μi of DP(4)\boldsymbol{\xi}_{k+1} = DP(\mathbf{x}^*)\,\boldsymbol{\xi}_k, \qquad \text{stable} \iff |\mu_i| < 1 \ \text{for all eigenvalues } \mu_i \text{ of } DP \qquad (4)

where: μ1,…,μn−1\mu_1, \ldots, \mu_{n-1} — the characteristic (Floquet) multipliers of the periodic orbit. Their relation to the flow's own linearization is the monodromy matrix over one period,

M(T)=Texp⁡(∫0TDf(ϕt(x∗)) dt),M(T) f(x∗)=f(x∗)(5)M(T) = \mathcal{T}\exp\Big(\int_0^T D\mathbf{f}\big(\boldsymbol{\phi}_t(\mathbf{x}^*)\big)\,\mathrm{d}t\Big), \qquad M(T)\,\mathbf{f}(\mathbf{x}^*) = \mathbf{f}(\mathbf{x}^*) \qquad (5)

where: Texp⁡\mathcal{T}\exp — the time-ordered exponential (the propagator of the variational equation from the Lyapunov article, evaluated over one period). The displayed identity is the key structural fact: the velocity vector itself solves the variational equation, so MM always has the eigenvalue 11 along the flow — useless for stability, and removed by the section: the spectrum of MM is {1}∪{μi}\{1\} \cup \{\mu_i\} (Derivation A). Volume bookkeeping crosses over intact. Liouville's formula for the variational determinant gives, per return,

det⁡DP(x∗)=exp⁡(∫0T∇⋅f dt)(6)\det DP(\mathbf{x}^*) = \exp\Big(\int_0^{T} \nabla\cdot\mathbf{f}\,\mathrm{d}t\Big) \qquad (6)

where: the integral runs along the orbit. Hamiltonian flow: ∇⋅f=0\nabla\cdot\mathbf{f} = 0, so det⁡DP=1\det DP = 1 — the return map is area-preserving (indeed symplectic), multipliers come in reciprocal pairs, and nothing can attract: the KAM world of curves, islands, and seas. Dissipative flow: ∣det⁡DP∣<1|\det DP| < 1 — dots are drawn inward, and fixed points, invariant circles, or fractal attractor sections await. The driven pendulum's det⁡DP=e−2βT≈0.009\det \mathrm{D}P = e^{-2\beta T} \approx 0.009 was equation (6) with constant divergence. Finally, the bifurcation dictionary — the classification this chapter has used piecemeal, now in one line:

μ=+1:saddle-node;μ=−1:period-doubling;μ=e±iω:Neimark–Sacker(7)\mu = +1: \text{saddle-node}; \qquad \mu = -1: \text{period-doubling}; \qquad \mu = e^{\pm i\omega}: \text{Neimark–Sacker} \qquad (7)

where: the exits through ±1\pm 1 were the working machinery of the cascade articles; the complex exit creates an attracting invariant curve around the fixed point — the section of a newborn torus, the discrete sibling of the Hopf bifurcation that opened this chapter.

Which statement answers which question

you wantusenotes
build the maptransversal Σ\Sigma + return (2)driven system: strobe at the drive period
find/certify a periodic orbitfixed point of PP, eq. (3)Newton on P(x)−xP(\mathbf{x}) - \mathbf{x} beats orbit-hunting
orbit stabilitymultipliers of DPDP, eq. (4)trivial multiplier already removed
which bifurcation just happenedexit dictionary (7)+1+1, −1-1, complex pair
Hamiltonian vs dissipative sectiondet⁡DP\det DP, eq. (6)=1= 1 area-preserving; <1< 1 attractors
read a computed sectiondots / curve / dust chartFigure 2; islands + sea = Hamiltonian mix
Derivation AWhy the return map exists, is smooth, and carries exactly the right multipliers

We are deriving the instrument's licence: that transversality makes PP well-defined and smooth, and that the section automatically discards the flow's trivial multiplier while keeping all the others. The strategy is the implicit function theorem for the first claim and a direct inspection of the monodromy spectrum for the second. The step to watch is the identification of the trivial eigenvector — one line of calculus that explains why map eigenvalues, not flow eigenvalues, are the correct stability currency for periodic orbits.

Step 1. Existence and smoothness of the return time. Fix x0∈Σ\mathbf{x}_0 \in \Sigma whose orbit returns to Σ\Sigma at time τ0\tau_0. Consider F(t,x)=g(ϕt(x))F(t, \mathbf{x}) = g(\boldsymbol{\phi}_t(\mathbf{x})); we need to solve F=0F = 0 for tt near (τ0,x0)(\tau_0, \mathbf{x}_0). The time-derivative is

∂F∂t=∇g⋅f(ϕt(x))≠0\frac{\partial F}{\partial t} = \nabla g\cdot\mathbf{f}\big(\boldsymbol{\phi}_t(\mathbf{x})\big) \neq 0

by transversality — so the implicit function theorem yields a unique smooth τ(x)\tau(\mathbf{x}) near x0\mathbf{x}_0, and P(x)=ϕτ(x)(x)P(\mathbf{x}) = \boldsymbol{\phi}_{\tau(\mathbf{x})}(\mathbf{x}) is a composition of smooth maps. Why transversality is non-negotiable: where the flow grazes the section (∇g⋅f=0\nabla g\cdot\mathbf{f} = 0), crossings merge or vanish, τ\tau jumps, and the "map" tears. Practical section-placing is the art of keeping the dynamics of interest away from tangencies.

Step 2. The trivial multiplier of the flow. Along the periodic orbit, differentiate ϕ˙t=f(ϕt)\dot{\boldsymbol\phi}_t = \mathbf{f}(\boldsymbol\phi_t) with respect to tt: the vector v(t)=f(ϕt(x∗))\mathbf{v}(t) = \mathbf{f}(\boldsymbol\phi_t(\mathbf{x}^*)) satisfies v˙=Df v\dot{\mathbf{v}} = D\mathbf{f}\,\mathbf{v} — it solves the variational equation. Over one period the orbit closes, so

M(T) f(x∗)=f(x∗)M(T)\,\mathbf{f}(\mathbf{x}^*) = \mathbf{f}(\mathbf{x}^*)

— the velocity is an eigenvector of the monodromy with eigenvalue exactly 11. Why this eigenvalue is dynamically empty: a perturbation along the orbit is the same orbit with a shifted clock; it neither converges nor diverges. Any stability verdict that keeps it (e.g. naively asking for ALL monodromy eigenvalues inside the unit circle) would declare every periodic orbit unstable — the classic beginner's paradox.

Step 3. The section as quotient. Decompose tangent space at x∗\mathbf{x}^* into TΣ⊕span(f)T\Sigma \oplus \mathrm{span}(\mathbf{f}). Deviations in Σ\Sigma evolve, per return, by the composition of M(T)M(T) with the projection back to Σ\Sigma along the flow (the return-time adjustment): this composition is precisely DPDP. Its spectrum is the monodromy spectrum with the trivial 11 removed:

spec M(T)={1}∪spec DP\mathrm{spec}\,M(T) = \{1\} \cup \mathrm{spec}\,DP

Step 4. Conclude the dictionary entries (3)–(4): the orbit is stable iff all map multipliers lie inside the unit circle, and every bifurcation of the orbit is a bifurcation of the fixed point of PP. The flow problem and the map problem are not analogous — they are the same problem, with the map holding the cleaner ledger.

Step 5. Note the smooth dependence bonus: because PP is smooth in states and parameters, Newton's method on P(x)−x=0P(\mathbf{x}) - \mathbf{x} = 0 locates periodic orbits (stable or unstable) and continuation tracks them across parameter space — the computational backbone of the bifurcation studies cited throughout this chapter. Verified: IFT construction per Strogatz §8.7 / Guckenheimer–Holmes §1.5; the Mf=fM\mathbf{f} = \mathbf{f} identity by direct differentiation; spectrum splitting per Floquet theory (Ott ch. 3).

Go deeper · Derivation B · 2 minArea bookkeeping on the section — and the one-number stability of a planar limit cycle
Derivation BArea bookkeeping on the section — and the one-number stability of a planar limit cycle

The goal is equation (6), det⁡DP=exp⁡∮∇⋅f dt\det DP = \exp\oint\nabla\cdot\mathbf{f}\,\mathrm{d}t, and its sharpest payoff: for a planar flow the section is one-dimensional and the single multiplier of a limit cycle is computable as a divergence integral around the loop. The strategy is Liouville's formula for the variational determinant, then the factor-out of the trivial multiplier. Watch how the abstract formula collapses, in the plane, to a checkable number for the van der Pol cycle.

Step 1. Liouville's formula along the orbit: ddtdet⁡M(t)=(∇⋅f) det⁡M(t)\frac{\mathrm{d}}{\mathrm{d}t}\det M(t) = \big(\nabla\cdot\mathbf{f}\big)\,\det M(t), hence over one period

det⁡M(T)=exp⁡(∫0T∇⋅f dt)\det M(T) = \exp\Big(\int_0^T \nabla\cdot\mathbf{f}\,\mathrm{d}t\Big)

Step 2. Factor the spectrum: det⁡M(T)=1⋅∏iμi=det⁡DP\det M(T) = 1 \cdot \prod_i \mu_i = \det DP. The trivial multiplier contributes the factor 1, so the flow's volume budget passes to the section undiluted — equation (6). Why this is the master consistency check: it is the periodic-orbit version of the Lyapunov sum rule, and any computed multiplier set violating it is wrong. For the driven pendulum, constant divergence −2β-2\beta gives det⁡DP=e−2βT\det DP = e^{-2\beta T} — the number used twice in that article, now derived in general.

Step 3. Specialize to the plane. A 2D flow has a 1D section: PP is a scalar map sk+1=P(sk)s_{k+1} = P(s_k) along a small transversal through the cycle γ\gamma, and stability is one number,

μ=P′(0)=exp⁡(∮γ∇⋅f dt)\mu = P'(0) = \exp\Big(\oint_\gamma \nabla\cdot\mathbf{f}\,\mathrm{d}t\Big)

with the cycle stable iff μ<1\mu < 1, i.e. iff the divergence integrated around the loop is negative. Why no separate trivial multiplier appears: in the plane the monodromy is 2×22\times 2; one eigenvalue is the trivial 1, so the determinant IS the other multiplier — the formula falls out of Step 2 with nothing left over.

Step 4. Check it on the chapter's favourite cycle. Van der Pol at small μ\mu: x¨−μ(1−x2)x˙+x=0\ddot x - \mu(1 - x^2)\dot x + x = 0, limit cycle x≈2cos⁡tx \approx 2\cos t, divergence ∇⋅f=μ(1−x2)\nabla\cdot\mathbf{f} = \mu(1 - x^2). Around the loop: ∮μ(1−4cos⁡2t) dt=μ(2π−4π)=−2πμ\oint \mu(1 - 4\cos^2 t)\,\mathrm{d}t = \mu(2\pi - 4\pi) = -2\pi\mu, so

μcycle=e−2πμ<1\mu_{\text{cycle}} = e^{-2\pi\mu} < 1

— stable, with perturbations shrinking by e−2πμe^{-2\pi\mu} per lap (about 0.530.53 at μ=0.1\mu = 0.1), consistent with the amplitude relaxation derived by averaging in the van der Pol article.

Step 5. Read the two structural regimes off (6). Hamiltonian: exponent zero, det⁡DP=1\det DP = 1, areas preserved — multipliers in reciprocal pairs (μ,1/μ\mu, 1/\mu: saddles or rotation pairs, never attraction), the section geometry of the double pendulum. Dissipative: negative average exponent — contraction, attractors, and (via strong contraction) the near-1D one-hump maps that made the cascade universal. Verified: Liouville/Jacobi formula per Ott §4.2; planar multiplier formula matches Strogatz §7.1's stability-of-cycles criterion; vdP integral ∮(1−4cos⁡2t)dt=−2π\oint(1-4\cos^2 t)\mathrm{d}t = -2\pi recomputed; e−0.628=0.534e^{-0.628} = 0.534.

Go deeper · Derivation C · 2 minThe torus on the section — circle maps, rotation numbers, and why irrationality means dense curves
Derivation CThe torus on the section — circle maps, rotation numbers, and why irrationality means dense curves

We derive the middle panel of the reading chart: why quasiperiodic motion on a 2-torus appears on the section as dots hopping around a closed invariant curve, filling it densely when the frequencies are incommensurate. The strategy is to take the torus flow in angle coordinates — the natural end-state of the action-angle machinery from the W4 articles — and compute its return map exactly: a rigid rotation of a circle. The step to watch is the rational/irrational dichotomy, which is also the door through which KAM theory and the small divisors enter.

Step 1. Motion on an invariant 2-torus in action-angle form: θ˙1=ω1\dot\theta_1 = \omega_1, θ˙2=ω2\dot\theta_2 = \omega_2, both angles mod 2π2\pi. (For the double pendulum's regular regions, this is the Liouville–Arnold normal form of the W4-4 article.)

Step 2. Choose the natural section Σ={θ1=0}\Sigma = \{\theta_1 = 0\}, transverse whenever ω1≠0\omega_1 \neq 0. Return time: θ1\theta_1 needs τ=2π/ω1\tau = 2\pi/\omega_1 to wrap once. During that lap, θ2\theta_2 advances by ω2τ\omega_2\tau:

P(θ2)=θ2+2πω2ω1(mod2π)≡θ2+2πρP(\theta_2) = \theta_2 + 2\pi\frac{\omega_2}{\omega_1} \pmod{2\pi} \equiv \theta_2 + 2\pi\rho

where ρ=ω2/ω1\rho = \omega_2/\omega_1 is the rotation number — the entire dynamics of the return map in one number. Why the map is a RIGID rotation: on the unperturbed torus the angles advance uniformly; all the potential complexity of the flow has been compressed into the single ratio of frequencies. The section has performed its dimensional reduction perfectly: 4D phase space → 3D energy surface → torus family → one circle map per torus.

Step 3. The dichotomy. If ρ=p/q\rho = p/q rational, every point of the circle is periodic with period qq: the section shows exactly qq dots, and the flow orbit is a closed curve winding pp times one way and qq times the other. If ρ\rho is irrational, no point ever returns exactly; the iterates θ2+2πnρ\theta_2 + 2\pi n\rho are all distinct and — by Weyl's equidistribution theorem — fill the circle densely and uniformly. The section shows dots accumulating into a smooth closed curve: the quasiperiodic signature of Figure 2.

Step 4. Perturb, and meet KAM on the section. Adding a small coupling makes the return map a circle map θ↦θ+2πρ+εh(θ)\theta \mapsto \theta + 2\pi\rho + \varepsilon h(\theta). Irrational-ρ\rho curves with good (Diophantine) irrationality survive as slightly deformed invariant curves — these are the KAM tori in section; rational-ρ\rho circles shatter into alternating elliptic/hyperbolic periodic points (Poincaré–Birkhoff), the island chains of the double pendulum's middle panel, with tangles and chaotic layers seeded at the hyperbolic points. Why the section is the natural home of this story: the entire KAM/Birkhoff phenomenology — which tori die, which survive, what replaces the dead ones — is a statement about circle maps and their rotation numbers, objects that exist only after the section has reduced the flow.

Step 5. Extract the measurable: from a computed or experimental section, the rotation number is estimated as (average angular advance per return)/2π/2\pi; its continued-fraction quality predicts robustness of the torus (the golden mean being the famous last survivor in the standard map). Rational lockings appear as finite dot counts — the mode-locking steps of driven oscillators. Verified: torus normal form per the Liouville–Arnold article; equidistribution per Weyl's theorem; Poincaré–Birkhoff island structure as in the double-pendulum sections; circle-map/KAM phenomenology per Ott ch. 6 and Strogatz §8.6.

Limiting cases

CaseSection behaviourMeaning
periodic orbit, period kTreturnkT_{\text{return}}kk dots foreverfixed point / kk-cycle of PP
quasiperiodic 2-torus, irrational ρ\rhodots fill a smooth closed curverigid rotation; dense by equidistribution
rational ρ=p/q\rho = p/qexactly qq dotsmode-locked; Poincaré–Birkhoff chains under perturbation
dissipative chaosfractal layered dust, ∣det⁡DP∣<1\lvert\det DP\rvert < 1strange-attractor cross-section
Hamiltonian chaosfat dust + islands + curves, det⁡DP=1\det DP = 1area-preserving mix; nothing attracts
flow tangent to Σ\Sigma somewheremap tears; crossings mergetransversality violated — move the pane
driven systemstroboscopic map, constant τ=T\tau = Tthe special case with a built-in clock

Worked examples

1 · Reading a stroboscopic ledger (confidence builder)

A mechanical oscillator driven at fd=50f_d = 50 Hz is sampled once per drive cycle. The section shows 6 dots, and the fitted return-map Jacobian at the associated cycle has eigenvalues 0.30.3 and −0.8-0.8. Classify the motion, its stability, and predict what happens if the second multiplier drifts to −1.02-1.02.

Step 1. Six dots = a 6-cycle of the stroboscopic map: the response repeats every 66 drive periods — a period-6 subharmonic at fd/6=8.33f_d/6 = 8.33 Hz, period 0.120.12 s.

Step 2. Stability: both multipliers inside the unit circle (∣0.3∣<1|0.3| < 1, ∣−0.8∣<1|-0.8| < 1): the subharmonic is stable. The negative sign of −0.8-0.8 means transverse deviations alternate sides at each return — harmless, but a warning of which exit is nearest.

Step 3. If the multiplier reaches −1.02-1.02: the −1-1 exit is the flip — the 6-cycle destabilizes and a stable 12-cycle is born; the section's 6 dots each split in two. By the cascade articles, expect further doublings at parameter intervals shrinking by δ=4.669\delta = 4.669, with chaos an accumulation point away.

2 · One number decides a limit cycle (exam level)

Using the planar multiplier formula, evaluate the stability of the van der Pol cycle at μ=0.1\mu = 0.1 and at μ=0.5\mu = 0.5, and give the number of laps needed for a 10% amplitude perturbation to decay below 1%.

Step 1. From Derivation B, the cycle's single multiplier is μcyc=exp⁡(∮∇⋅f dt)=e−2πμ\mu_{\text{cyc}} = \exp\big(\oint\nabla\cdot\mathbf{f}\,\mathrm{d}t\big) = e^{-2\pi\mu} in the small-μ\mu regime (using x≈2cos⁡tx \approx 2\cos t).

Step 2. At μ=0.1\mu = 0.1: μcyc=e−0.628=0.53\mu_{\text{cyc}} = e^{-0.628} = \boxed{0.53} — perturbations halve each lap. At μ=0.5\mu = 0.5: e−3.14=0.043e^{-3.14} = \boxed{0.043} — near-deadbeat return to the cycle in one lap (though at μ=0.5\mu = 0.5 the x=2cos⁡tx = 2\cos t approximation is already fraying; treat as an estimate).

Step 3. Decay laps at μ=0.1\mu = 0.1: need 0.53n<0.10.53^n < 0.1, i.e. n>ln⁡10/0.628=3.7n > \ln 10/0.628 = 3.7 — 4\boxed{4} laps. The section converts "how stable is this oscillation" from a vague adjective into a per-lap contraction factor.

3 · Rotation numbers and dot counts (applied)

A two-frequency system has ω1=1.000\omega_1 = 1.000 and ω2=0.625\omega_2 = 0.625 (rad/s). (i) Describe its section at θ1=0\theta_1 = 0. (ii) Repeat for ω2=0.618034\omega_2 = 0.618034. (iii) How would you distinguish the two cases from 200 measured returns?

Step 1. ρ=0.625=5/8\rho = 0.625 = 5/8: rational. The return map is rotation by 2π⋅5/82\pi\cdot 5/8; every orbit closes after 8 returns. Section: exactly 8\boxed{8} dots, visited in the fixed order 0,5,2,7,4,1,6,30, 5, 2, 7, 4, 1, 6, 3 (multiples of 5 mod 8) — a mode-locked 8:58{:}5 resonance.

Step 2. ρ=0.618034≈(5−1)/2\rho = 0.618034 \approx (\sqrt5 - 1)/2: the golden mean, maximally irrational. Dots never repeat; after nn returns they fill the invariant circle to spacing ∼2π/n\sim 2\pi/n with famous uniformity. Section: a smooth closed curve materializing dot by dot.

Step 3. With 200 returns: case (i) shows 8 point-clusters whose widths are set only by noise — additional returns revisit clusters, never fill gaps; case (ii) shows 200 distinct points with maximum gap shrinking like 1/n1/n. Quantitatively, plot the number of distinct dot positions (at measurement resolution) versus nn: saturation at 8 versus linear growth is unambiguous by n≈30n \approx 30. This dot-count census is exactly how mode-locking tongues are mapped in driven-oscillator experiments.

Origins · 1 minHenri Poincaré and how this was discovered
HP
Henri Poincaré
1854–1912 · Paris

Poincaré appears twice in this chapter's gallery — the Poincaré–Bendixson article told the planar story — and the section is his second, larger bequest. It was forged in the fire of the most famous error in mathematical physics. His 1888 entry for King Oscar II's prize proved the restricted three-body problem stable — using precisely the section-and-return-map method he had invented for the purpose. After the prize was awarded, the young editor Lars Phragmén pressed him on an obscure step; Poincaré discovered the flaw ran deep: the stable and unstable manifolds of a periodic orbit's saddle point on the section, which he had implicitly assumed to rejoin smoothly, must instead cross transversally — and then cross infinitely often, forced by the map itself. The already-printed memoir was recalled and destroyed (Poincaré paid the printing costs, exceeding the prize), and the corrected 1890 version contains the first description of a homoclinic tangle: "one is struck by the complexity of this figure, which I shall not even attempt to draw." In the Méthodes nouvelles he went further — near a tangle lie infinitely many periodic orbits interlaced with motions of arbitrary complexity: deterministic chaos, seen and named as a phenomenon of the return map, seventy years before computers drew the pictures. The instrument and the discovery are inseparable: chaos was found on a section.

History
1881
— Poincaré begins the memoir series on curves defined by differential equations, inventing the qualitative theory and the transversal-section idea.
1890
— The corrected prize memoir on the three-body problem describes homoclinic intersections — the tangle, and with it the first mathematical chaos.
1899
— Les Méthodes nouvelles de la mécanique céleste, volume III, systematizes sections, return maps, and the infinite web of periodic orbits near a tangle.
1963
— Hénon and Heiles compute sections of a two-degree-of-freedom system, turning the instrument into the standard numerical diagnostic.
1967
— Smale's horseshoe paper proves every transverse homoclinic tangle contains chaotic symbolic dynamics — Poincaré's undrawable figure, domesticated.

Going further

References. Strogatz, Nonlinear Dynamics and Chaos, 2nd ed., §8.7 (Poincaré maps) · Guckenheimer & Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, §1.5 · Ott, Chaos in Dynamical Systems, 2nd ed., ch. 3–4 · Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., ch. 11 · Lichtenberg & Lieberman, Regular and Chaotic Dynamics, 2nd ed. · Hénon & Heiles (1964), Astron. J. 69, 73

Provenance. Construction of PP via transversality and the implicit function theorem follows Guckenheimer–Holmes §1.5 and Strogatz §8.7. Monodromy structure, trivial multiplier Mf=fM\mathbf{f}=\mathbf{f}, and det⁡DP=exp⁡∫∇⋅f dt\det DP=\exp\int\nabla\cdot\mathbf{f}\,\mathrm{d}t are standard Floquet/Liouville results (Ott ch. 3–4); a verifier should re-derive the trivial multiplier by differentiating the flow equation and the determinant formula from Jacobi's identity. The planar cycle multiplier exp⁡∮∇⋅f dt\exp\oint\nabla\cdot\mathbf{f}\,\mathrm{d}t matches Strogatz §7.1; the van der Pol evaluation ∮μ(1−4cos⁡2t) dt=−2πμ\oint\mu(1-4\cos^2t)\,\mathrm{d}t=-2\pi\mu is recomputed directly. The bifurcation dictionary at μ=+1,−1,e±iω\mu=+1,-1,e^{\pm i\omega} is textbook (Guckenheimer–Holmes ch. 3). Torus section = rigid rotation and the rational/irrational dichotomy follow the action-angle normal form with Weyl equidistribution.

Builds on: Phase Portraits and the Phase Plane · Fixed Points and Equilibria · Linear Stability Analysis and the Jacobian · Classification of Fixed Points: Nodes, Saddles, Spirals, Centers · Conservative Systems and Hamiltonian Phase Flow Leads to: Poincaré–Bendixson Theorem · Limit Cycles and Self-Sustained Oscillations · Hopf Bifurcation and Birth of Limit Cycles · Period-Doubling Cascade and the Route to Chaos · Bifurcations and the Loss of Stability

Does the choice of section surface change the results?

Not the results that matter — and this is a theorem, not a convention. Return maps built on two different transversals to the same orbit region are smoothly conjugate: P2=h∘P1∘h−1P_2 = h \circ P_1 \circ h^{-1}, where hh maps one surface to the other along the flow. Conjugacy preserves everything qualitative and most things quantitative: the count and periods of periodic orbits, their multipliers (eigenvalues are conjugation-invariant), the bifurcation structure, Lyapunov exponents, rotation numbers, and the topology of invariant curves and dust. What DOES change is the picture's parametrization — dots land at different coordinates, curves are deformed — and, near a bad choice, the practical quality: a pane close to tangency has wildly varying return times and distorted, hard-to-read portraits. So the honest statement is: section choice is a gauge choice; pick the gauge in which your structure is legible (for driven systems the drive phase is nature's suggestion), and never attribute physical meaning to features that a re-slicing would move.

The return time varies from point to point — is that information just thrown away?

It is set aside, not destroyed, and it can be reattached whenever needed. The section deliberately quotients out when things happen to expose what happens — that is its power. But the pair (return map PP, return-time function τ\tau) reconstructs the flow completely, up to smooth equivalence, via the suspension construction: flow along a "roof" of height τ(x)\tau(\mathbf{x}) over the section and glue by PP. Quantities that need the clock use τ\tau explicitly: the period of an orbit through a kk-cycle is the sum of the kk return times; flow Lyapunov exponents are map exponents divided by the average return time (λflow=λmap/⟨τ⟩\lambda_{\text{flow}} = \lambda_{\text{map}}/\langle\tau\rangle); power spectra need the actual crossing times. One caution runs the other way: for driven systems the stroboscopic clock is rigid, so subharmonics are exact integer multiples of the drive period — while for autonomous systems the "period" read from dot counts is k⟨τ⟩k\langle\tau\rangle with τ\tau's variation contributing phase jitter that a naive spectral analysis of the dots misses.

How do sections work above three dimensions — does the picture advantage die?

The construction is dimension-blind — one flow dimension is always removed — but the display advantage needs supplementing. A 4D flow gives a 3D map: already hard to draw. The standard tricks: (1) For Hamiltonian systems, fix the energy first (3D surface), THEN section: a 2D area-preserving map — this is exactly the double-pendulum pipeline, and it works for any two-degree-of-freedom system. (2) For higher-dimensional dissipative systems, strong contraction usually collapses the section dynamics onto a low-dimensional attractor within it, so 2D projections of the section remain informative (the driven pendulum in higher-dimensional guise). (3) When neither applies, one keeps the map algebraically — Newton for periodic orbits, multipliers for stability, all of Derivations A–B — and forgoes pictures; the section's analytic content never degrades. (4) Double sections (sectioning the section) can peel further dimensions for special structures. The instrument's soul — trading a flow for a map with the trivial direction removed — is dimension-independent; only the poster-friendliness is a low-dimensional luxury.

What is so terrible about one homoclinic crossing — why does it force infinitely many and then chaos?

Because crossings breed under iteration. Let q≠x∗\mathbf{q} \neq \mathbf{x}^* lie on both the stable and unstable manifolds of the saddle x∗\mathbf{x}^* (a homoclinic point), with the crossing transverse. Both manifolds are invariant: PP maps each to itself. Hence P(q)P(\mathbf{q}), P2(q)P^2(\mathbf{q}), … are ALL homoclinic points too — marching toward x∗\mathbf{x}^* along the stable manifold — as are the preimages, marching in along the unstable one. Infinitely many crossings, from one. Now watch the geometry strain: between consecutive crossings near x∗\mathbf{x}^*, the unstable manifold must connect points that are exponentially close along the stable direction, while manifolds may never self-intersect; the only resolution is that WuW^u folds into wilder and longer oscillations — exponentially stretched lobes crammed into shrinking windows. That lobe structure IS stretch-and-fold: Smale's theorem makes it exact by locating a horseshoe map in some iterate of PP near the tangle, whence symbolic dynamics on two symbols — uncountably many aperiodic orbits, a dense infinity of periodic ones, positive entropy, sensitive dependence. This is also the practical moral: prove one transverse homoclinic intersection (Melnikov's method does it perturbatively) and chaos is certified, no simulation required.

What is the right way to compute a section numerically?

Three rules cover most of it. (1) Never take the nearest stored step. Integrate until the sign of g(x)g(\mathbf{x}) flips, then locate the crossing properly: either interpolate the last step to g=0g = 0 (a root-solve on the dense output), or — the elegant classical trick (Hénon's) — switch the independent variable from tt to gg for one final partial step, so the integrator lands exactly on the section. Sloppy crossing location smears dots by O(Δt)O(\Delta t), which is enough to blur island chains and fake fractality. (2) Respect orientation and transversality: record only ∇g⋅f>0\nabla g\cdot\mathbf{f} > 0 crossings (or only <0< 0), and check that no recorded crossing has small ∣∇g⋅f∣|\nabla g\cdot\mathbf{f}| — near-tangencies deserve either a moved pane or special care. (3) Match the integrator to the structure: for Hamiltonian sections use a symplectic integrator, or the numerical dissipation will slowly spiral your KAM curves (the double-pendulum thread on this is instructive); discard transients for dissipative systems; and for return-map Jacobians, integrate the variational equations between crossings and apply the section projection of Derivation A rather than finite-differencing dots.

Are the multipliers of the map related to the Lyapunov exponents of the flow?

Directly — they are the same information at two levels of averaging. For a periodic orbit, the relation is exact and per-orbit: μi=eλi(orbit)T\mu_i = e^{\lambda_i^{(orbit)} T}, where TT is the period and λi(orbit)\lambda_i^{(orbit)} are the orbit's Floquet exponents; stability "inside the unit circle" and "negative real part" are the same statement through the exponential. For chaotic trajectories, the section's iterated Jacobians DP(xk)⋯DP(x1)DP(\mathbf{x}_k)\cdots DP(\mathbf{x}_1) define map Lyapunov exponents per return, and the flow exponents are recovered by dividing by the mean return time: λflow=λmap/⟨τ⟩\lambda^{\text{flow}} = \lambda^{\text{map}}/\langle\tau\rangle — with the flow's trivial zero exponent being exactly the direction the section removed (a map has no obligatory zero). This bookkeeping is why the driven pendulum's per-flash contraction det⁡DP=e−2βT\det DP = e^{-2\beta T} matched its exponent sum −2β-2\beta times TT, and why chaos certified on a section (positive map exponent) certifies chaos in the flow with no further work.

Before this

Where this leads