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 ; 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 — and the map knows everything qualitative about the flow while living in one dimension fewer. Periodic orbits become fixed points; stability becomes eigenvalues of (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 , 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.
- A Poincaré section is a codimension-one surface transverse to the flow (); successive same-direction crossings define the first-return map , an -dimensional discrete system.
- Fixed points of are periodic orbits of the flow; the eigenvalues of are the Floquet multipliers, with the trivial along-flow multiplier automatically removed.
- Stability: all ; the three unit-circle exits are the bifurcation dictionary — saddle-node, period-doubling, complex pair Neimark–Sacker (torus birth).
- Area bookkeeping transfers exactly: — 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 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.
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-, for period-, 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.
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 — 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 are exactly the deviations that matter, and the eigenvalues of 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 collides orbits (saddle-node), through 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 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.
The instrument, formalized
The setup: an autonomous flow and a well-placed pane of glass,
where: — any smooth scalar whose zero set is the section surface; the transversality condition says the flow crosses 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 , and with it the Poincaré map
where: — the flow map (solution operator); maps the -dimensional surface to itself. Note that the return time varies from point to point — the stroboscopic map of driven systems, with fixed , is the special case where the drive phase supplies a section with constant return time. The dictionary begins with periodic orbits:
where: a -cycle of likewise corresponds to a closed orbit crossing the pane times. Stability of the orbit is stability of the fixed point: linearize the map,
where: — the characteristic (Floquet) multipliers of the periodic orbit. Their relation to the flow's own linearization is the monodromy matrix over one period,
where: — 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 always has the eigenvalue along the flow — useless for stability, and removed by the section: the spectrum of is (Derivation A). Volume bookkeeping crosses over intact. Liouville's formula for the variational determinant gives, per return,
where: the integral runs along the orbit. Hamiltonian flow: , so — 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: — dots are drawn inward, and fixed points, invariant circles, or fractal attractor sections await. The driven pendulum's was equation (6) with constant divergence. Finally, the bifurcation dictionary — the classification this chapter has used piecemeal, now in one line:
where: the exits through 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 want | use | notes |
|---|---|---|
| build the map | transversal + return (2) | driven system: strobe at the drive period |
| find/certify a periodic orbit | fixed point of , eq. (3) | Newton on beats orbit-hunting |
| orbit stability | multipliers of , eq. (4) | trivial multiplier already removed |
| which bifurcation just happened | exit dictionary (7) | , , complex pair |
| Hamiltonian vs dissipative section | , eq. (6) | area-preserving; attractors |
| read a computed section | dots / curve / dust chart | Figure 2; islands + sea = Hamiltonian mix |
We are deriving the instrument's licence: that transversality makes 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 whose orbit returns to at time . Consider ; we need to solve for near . The time-derivative is
by transversality — so the implicit function theorem yields a unique smooth near , and is a composition of smooth maps. Why transversality is non-negotiable: where the flow grazes the section (), crossings merge or vanish, 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 with respect to : the vector satisfies — it solves the variational equation. Over one period the orbit closes, so
— the velocity is an eigenvector of the monodromy with eigenvalue exactly . 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 into . Deviations in evolve, per return, by the composition of with the projection back to along the flow (the return-time adjustment): this composition is precisely . Its spectrum is the monodromy spectrum with the trivial removed:
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 . 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 is smooth in states and parameters, Newton's method on 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 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
The goal is equation (6), , 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: , hence over one period
Step 2. Factor the spectrum: . 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 gives — the number used twice in that article, now derived in general.
Step 3. Specialize to the plane. A 2D flow has a 1D section: is a scalar map along a small transversal through the cycle , and stability is one number,
with the cycle stable iff , i.e. iff the divergence integrated around the loop is negative. Why no separate trivial multiplier appears: in the plane the monodromy is ; 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 : , limit cycle , divergence . Around the loop: , so
— stable, with perturbations shrinking by per lap (about at ), 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, , areas preserved — multipliers in reciprocal pairs (: 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 recomputed; .
Go deeper · Derivation C · 2 minThe 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: , , both angles mod . (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 , transverse whenever . Return time: needs to wrap once. During that lap, advances by :
where 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 rational, every point of the circle is periodic with period : the section shows exactly dots, and the flow orbit is a closed curve winding times one way and times the other. If is irrational, no point ever returns exactly; the iterates 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 . Irrational- curves with good (Diophantine) irrationality survive as slightly deformed invariant curves — these are the KAM tori in section; rational- 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); 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
| Case | Section behaviour | Meaning |
|---|---|---|
| periodic orbit, period | dots forever | fixed point / -cycle of |
| quasiperiodic 2-torus, irrational | dots fill a smooth closed curve | rigid rotation; dense by equidistribution |
| rational | exactly dots | mode-locked; Poincaré–Birkhoff chains under perturbation |
| dissipative chaos | fractal layered dust, | strange-attractor cross-section |
| Hamiltonian chaos | fat dust + islands + curves, | area-preserving mix; nothing attracts |
| flow tangent to somewhere | map tears; crossings merge | transversality violated — move the pane |
| driven system | stroboscopic map, constant | the special case with a built-in clock |
Worked examples
1 · Reading a stroboscopic ledger (confidence builder)
A mechanical oscillator driven at Hz is sampled once per drive cycle. The section shows 6 dots, and the fitted return-map Jacobian at the associated cycle has eigenvalues and . Classify the motion, its stability, and predict what happens if the second multiplier drifts to .
Step 1. Six dots = a 6-cycle of the stroboscopic map: the response repeats every drive periods — a period-6 subharmonic at Hz, period s.
Step 2. Stability: both multipliers inside the unit circle (, ): the subharmonic is stable. The negative sign of means transverse deviations alternate sides at each return — harmless, but a warning of which exit is nearest.
Step 3. If the multiplier reaches : the 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 , 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 and at , 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 in the small- regime (using ).
Step 2. At : — perturbations halve each lap. At : — near-deadbeat return to the cycle in one lap (though at the approximation is already fraying; treat as an estimate).
Step 3. Decay laps at : need , i.e. — 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 and (rad/s). (i) Describe its section at . (ii) Repeat for . (iii) How would you distinguish the two cases from 200 measured returns?
Step 1. : rational. The return map is rotation by ; every orbit closes after 8 returns. Section: exactly dots, visited in the fixed order (multiples of 5 mod 8) — a mode-locked resonance.
Step 2. : the golden mean, maximally irrational. Dots never repeat; after returns they fill the invariant circle to spacing 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 . Quantitatively, plot the number of distinct dot positions (at measurement resolution) versus : saturation at 8 versus linear growth is unambiguous by . 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
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.
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 via transversality and the implicit function theorem follows Guckenheimer–Holmes §1.5 and Strogatz §8.7. Monodromy structure, trivial multiplier , and 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 matches Strogatz §7.1; the van der Pol evaluation is recomputed directly. The bifurcation dictionary at 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: , where 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 , return-time function ) reconstructs the flow completely, up to smooth equivalence, via the suspension construction: flow along a "roof" of height over the section and glue by . Quantities that need the clock use explicitly: the period of an orbit through a -cycle is the sum of the return times; flow Lyapunov exponents are map exponents divided by the average return time (); 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 with '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 lie on both the stable and unstable manifolds of the saddle (a homoclinic point), with the crossing transverse. Both manifolds are invariant: maps each to itself. Hence , , … are ALL homoclinic points too — marching toward 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 , 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 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 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 flips, then locate the crossing properly: either interpolate the last step to (a root-solve on the dense output), or — the elegant classical trick (Hénon's) — switch the independent variable from to for one final partial step, so the integrator lands exactly on the section. Sloppy crossing location smears dots by , which is enough to blur island chains and fake fractality. (2) Respect orientation and transversality: record only crossings (or only ), and check that no recorded crossing has small — 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: , where is the period and 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 define map Lyapunov exponents per return, and the flow exponents are recovered by dividing by the mean return time: — 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 matched its exponent sum times , and why chaos certified on a section (positive map exponent) certifies chaos in the flow with no further work.