---
title: Whiskered Tori in Hamiltonian Dynamics
url: https://www.emergentmind.com/topics/whiskered-tori
type: topic
---

# Whiskered Tori in Hamiltonian Dynamics

Searching arXiv for relevant papers on whiskered tori and related methodologies.
I’m looking up recent and foundational arXiv papers on whiskered tori to ground the article.
Whiskered tori are invariant quasi-periodic tori endowed with hyperbolic directions transverse to the torus. In the standard Hamiltonian setting, an $\ell$-dimensional invariant torus carries $\ell$ neutral directions tangent to the torus together with stable and unstable directions in the normal dynamics; the associated stable and unstable manifolds are the “whiskers.” In integrable models these whiskers may coincide along homoclinic separatrices, whereas under perturbation they typically split, creating transverse homoclinic or heteroclinic structures that organize transport and instability. The modern literature treats whiskered tori as a unifying object across exact symplectic maps and flows, nearly-integrable Hamiltonians, conformally symplectic systems, coupled map lattices, Hamiltonian PDEs, and celestial-mechanical models [1306.0728] [1004.5231].

## 1. Geometric notion and invariant splitting

In coupled map lattices, a whiskered torus is described by an analytic embedding
\[
K:\mathbb T^\ell \to \mathcal M
\]
together with an invariant splitting
\[
T_{K(\theta)}\mathcal M
=
E^s_{K(\theta)}\oplus E^c_{K(\theta)}\oplus E^u_{K(\theta)},
\]
where $E^s$, $E^c$, and $E^u$ are respectively exponentially contracting, neutral, and expanding under the linearized cocycle [1307.2374]. In exact symplectic finite-dimensional dynamics, the same structure appears as
\[
T_{K(\theta)}M=E^s_{K(\theta)}\oplus E^c_{K(\theta)}\oplus E^u_{K(\theta)},
\]
with $\dim E^c=2\ell$, the center bundle being generated by $DK(\theta)$ and its symplectic conjugate $J^{-1}DK(\theta)$ [1004.5231].

A particularly concrete realization arises in the periodically perturbed planar circular restricted three-body problem. There, in the 4D stroboscopic map, a one-dimensional invariant torus $K:\mathbb T\to\mathbb R^4$ is whiskered when
\[
F_\varepsilon(K(\theta))=K(\theta+\omega)
\]
and the derivative admits one tangent direction, one symplectic-conjugate center direction, and hyperbolic stable and unstable directions with multipliers $\lambda_s(\theta)$ and $\lambda_u(\theta)=1/\lambda_s(\theta)$ [2105.11100]. In extended phase space for periodically perturbed PCRTBP, autonomous periodic orbits of the unperturbed problem typically persist as two-dimensional invariant tori whose linearized stroboscopic map has two neutral directions and one contracting plus one expanding multiplier [2109.14814].

The designation “whiskered” therefore refers not to the torus itself but to the existence of invariant stable and unstable manifolds attached to it. A common simplification is to identify whiskered tori only with finite-dimensional, exact-symplectic Hamiltonian systems. The literature instead includes localized tori in infinite lattices [1307.2374], one-dimensional invariant circles in 4D symplectic maps [2105.11100], and normally parabolic tori at infinity in the planar $(n+1)$-body problem [1812.01286].

## 2. Invariance equations, persistence, and nondegeneracy

The basic functional equation for a whiskered torus is the conjugacy relation
\[
F(K(\theta))=K(\theta+\omega)
\]
for maps, or
\[
\partial_\omega K(\theta)=X(K(\theta))
\]
for flows [1004.5231]. In coupled map lattices this is written as $F\circ K=K\circ T_\omega$, with $T_\omega(\theta)=\theta+\omega$ [1307.2374]. In Hamiltonian PDEs the same structure appears as
\[
\partial_\omega K(\theta)=\mathcal L K(\theta)+\mathcal N(K(\theta)),
\]
where the linearization along the torus admits a stable-center-unstable splitting with infinite-dimensional stable and unstable bundles and a finite-dimensional center bundle of dimension $2\ell$ [1602.03775].

Persistence results are typically formulated in an a-posteriori KAM form. One starts from an approximate embedding $K_0$ with small defect, an approximately invariant splitting, and explicit nondegeneracy data, and concludes the existence of a true nearby whiskered torus. In exact symplectic dynamics this requires Diophantine frequency, twist, and hyperbolicity bounds computable from the approximate solution [1004.5231]. In conformally symplectic systems one seeks $(K,\mu)$ satisfying
\[
f_\mu\circ K(\theta)=K(\theta+\omega),
\]
with $f_\mu^\ast\Omega=\lambda \Omega$, and the Newton scheme simultaneously corrects the torus embedding, the drift parameter, and the invariant splitting under an exponential trichotomy [1901.07483]. In Hamiltonian PDEs, the corresponding theorem combines Diophantine estimates on $\omega$, smoothing semigroups in the hyperbolic directions, and a finite-dimensional twist condition on the center reduction [1602.03775].

These formulations make precise that whiskered tori are not merely normally hyperbolic invariant manifolds in an abstract sense. Their persistence theory is tied to cohomological equations over quasi-periodic rotations, small-divisor estimates, and explicit rate conditions on the linearized cocycle. In exact and conformally symplectic contexts, the center dynamics is further constrained by the geometry of the symplectic or conformally symplectic form [1004.5231] [1901.07483].

## 3. Whiskers, parameterization methods, and fast computation

Once the torus is known, one parameterizes its stable and unstable manifolds by an invariance equation of the form
\[
F(W(\theta,s))=W(\theta+\omega,\lambda s),
\]
where $\lambda$ is the stable or unstable multiplier [2105.11100]. In the PCRTBP literature this is expanded as a Fourier-Taylor series
\[
W(\theta,s)=K(\theta)+\sum_{k=1}^d W_k(\theta)s^k,
\]
with $W_1(\theta)$ equal to the stable or unstable direction and higher coefficients obtained order by order [2105.11100]. In the coupled-map-lattice setting the local stable manifold is parameterized by
\[
W:\mathbb T^\ell\times B^s_\rho\to\mathcal M,
\qquad
F(W(\theta,s))=W(\theta+\omega,P(\theta,s)),
\]
where $P(\theta,s)=A^s(\theta)s+\text{(higher-order terms)}$ encodes the internal stable dynamics [1307.2374].

The dominant computational methodology is the parameterization method combined with Newton or quasi-Newton iteration. For exact symplectic maps, one reduces the Newton equation by splitting the correction into stable, center, and unstable components and using automatic reducibility on the center bundle. If the torus is discretized by $N$ elements, one Newton step requires only $O(N)$ storage and $O(N\ln N)$ operations, and only functions of $\ell$ variables need to be computed, independently of the phase-space dimension [1004.5231]. In the periodically perturbed RTBP, a quasi-Newton method computes simultaneously the torus, the bundle frame $P(\theta)$, and the block-triangular cocycle $\Lambda(\theta)$, with each step reducing both the torus defect and the bundle defect quadratically [2105.11100].

A more recent flow-map formulation computes the torus and its one-dimensional whiskers simultaneously in autonomous Hamiltonian systems. There the time-$T$ map $\phi_T$ is used to solve
\[
\phi_T\circ K(\theta)=K(\theta+\omega),
\qquad
\phi_T(W(\theta,s))=W(\theta+\omega,\lambda s),
\]
with Fourier-Taylor or grid-Taylor discretization and a Newton-like step based on an adapted symplectic frame and block-upper-triangular reduction of the linearized map [2507.06123]. This suggests a methodological trend: rather than computing the torus first and the whiskers afterward, one exploits their coupled invariance equations to improve efficiency and numerical stability.

For global connections between whiskered tori, local parameterizations are combined with globalization and search procedures. In periodically perturbed PCRTBP, 2D manifolds are meshed on $(\theta,s)$ grids, globalized by iterating the map, and searched for heteroclinic intersections by broad-phase 4D grid partitioning, axis-aligned bounding-box rejection, and Möller’s quick triangle-triangle test projected to $(x,y,p_x)$ [2109.14814]. The resulting approximate intersections are then refined by damped Newton iteration applied to
\[
f(\theta_u,s_u,\theta_s,s_s)=W_1^u(\theta_u,s_u)-W_2^s(\theta_s,s_s)=0,
\]
reducing residuals from $10^{-2}$ to $10^{-7}$ or better [2109.14814].

## 4. Splitting of separatrices and arithmetic dependence

In nearly-integrable Hamiltonian systems with a pendulum-like hyperbolic part, whiskered tori are central to the theory of exponentially small splitting of separatrices. For
\[
H_0(x,y,\varphi,I)=P(x,y)+\langle \omega, I\rangle+\tfrac12\langle \Lambda I,I\rangle,
\]
with a hyperbolic pendulum degree of freedom and fast rotors, the torus
\[
T_0=\{x=0,y=0,I=0\}
\]
is whiskered, and in the integrable case its stable and unstable manifolds coincide along homoclinic orbits [1306.0728]. After perturbation $H=H_0+\mu H_1$, the splitting is measured by a periodic splitting function $\mathcal M(\theta)$ or, equivalently, by a splitting potential whose gradient is $\mathcal M$ [1306.0728].

The first-order approximation is provided by the Poincaré-Melnikov method. In the quadratic and cubic frequency settings considered by Delshams, Gonchenko, and Gutiérrez, the Melnikov potential is
\[
L(\theta)= -\int_{-\infty}^{\infty}[h(x_0(t))-h(0)]\,f(\theta+\omega_\varepsilon t)\,dt,
\qquad
\mathcal M(\theta)=M(\theta)+O(\mu^2),
\]
with $M(\theta)=\nabla L(\theta)$ [1306.0728]. Because the homoclinic excursion is exponentially long in fast time and because the Fourier coefficients involve small denominators $\langle k,\omega\rangle$, the maximal splitting distance is exponentially small.

For $\ell=2$ and $\ell=3$, if $\mu=\varepsilon^p$ with $p>3$, the maximal splitting distance satisfies
\[
\Delta(\varepsilon)\asymp
\mu\,\varepsilon^{-1/\ell}
\exp\!\left(
-\frac{C_0\,h_1(\varepsilon)}{\varepsilon^{1/(2\ell)}}
\right),
\]
where $C_0$ depends on $\rho$ and on a positive limit numerator $\gamma^\ast$ determined by the Diophantine properties of the frequency vector [1306.0728]. In the quadratic case $\beta=1/4$ and
\[
\Delta(\varepsilon)\asymp
(\mu/\varepsilon^{1/2})
\exp\!\left(-C_0 h_1(\varepsilon)/\varepsilon^{1/4}\right),
\]
with $h_1(\varepsilon)$ periodic in $\ln\varepsilon$ and bounded by $1\le h_1(\varepsilon)\le A_1$ [1306.0728]. In the cubic golden case $\beta=1/6$ and
\[
\Delta(\varepsilon)\asymp
(\mu/\varepsilon^{1/3})
\exp\!\left(-C_0 h_1(\varepsilon)/\varepsilon^{1/6}\right),
\]
but $h_1(\varepsilon)$ is bounded rather than strictly periodic [1306.0728].

This arithmetic dependence is developed in several directions. For quadratic irrational ratios, the dominant harmonics are organized by resonant sequences derived from the continued fraction of $\Omega$, and the functions governing the exponent are periodic in $\ln\varepsilon$ with period $4\ln\lambda$ [1407.6524] [1507.07397]. For the silver ratio $\Omega=\sqrt2-1$, the splitting function has exactly four simple zeros for all sufficiently small $\varepsilon$, yielding four transverse homoclinic orbits, and the maximal splitting and transversality inherit the same periodic modulation in $\ln\varepsilon$ [1409.4944]. For irrational numbers of constant type, continued-fraction boundedness yields exponentially small lower bounds of order
\[
\max_{\theta\in\mathbb T^2}\|\Delta(\theta)\|
\succeq
\frac{\mu}{\varepsilon^{1/2}}
\exp\!\left(-\frac{C(\Omega,\rho)}{\varepsilon^{1/4}}\right),
\]
with $C(\Omega,\rho)$ explicitly controlled by Diophantine constants and continued-fraction data [1402.1654].

A common simplification is to treat the logarithmic modulation in the exponent as universally periodic. The cubic-frequency theory shows that this is not the case. For 3D whiskered tori with cubic frequencies, the lack of a one-dimensional continued-fraction algorithm is replaced by a unimodular Koch matrix, resonant sequences in $\mathbb Z^3$, and a function $h_1(\varepsilon)$ that is quasiperiodic, not periodic, with respect to $\ln\varepsilon$ [1906.01439].

## 5. Extensions beyond exact symplectic finite-dimensional tori

Whiskered tori persist in dissipative geometries. For conformally symplectic maps satisfying
\[
f_\mu^\ast\Omega=\lambda\Omega,
\]
one fixes a Diophantine frequency $\omega$, an approximate embedding $K_0$, an approximate parameter $\mu_0$, and an approximately invariant trichotomy
\[
T_{K(\theta)}M=E^s(\theta)\oplus E^c(\theta)\oplus E^u(\theta),
\]
with $\dim E^c=2d$, and proves the existence of a nearby exact whiskered torus and drift parameter under explicit twist and rate conditions [1901.07483]. The center correction is solved through automatic reducibility and a $2d\times 2d$ averaged nondegeneracy matrix [1901.07483]. This suggests that the essential KAM mechanism survives even when the symplectic form is only preserved up to a scalar factor.

Infinite-dimensional and spatially extended settings require different analytic frameworks. In Hamiltonian PDEs, whiskered tori are embedded quasi-periodic solutions
\[
K:\mathbb T^\ell\supset D_\rho\to X
\]
for which the linearized equation admits analytic stable, center, and unstable bundles, with $\dim X^c_\theta=2\ell$ and $\dim X^s_\theta=\dim X^u_\theta=\infty$ [1602.03775]. The stable and unstable dynamics are controlled by smoothing semigroups, while the center equation is again reduced to cohomological equations via a moving frame [1602.03775]. In coupled Hamiltonian lattices and coupled map lattices, the decisive device is the use of decay functions and weighted Banach spaces, which retain locality and allow hyperbolic splittings and invariant manifolds to remain localized near excited sites [1307.2374] [1406.1091].

The term also extends to normally parabolic regimes. In the planar $(n+1)$-body problem, one encounters parabolic tori at infinity whose normal linearization is trivial and whose first nonzero normal term is of higher order. Such a torus is called whiskered when it admits one-dimensional stable and unstable manifolds of dimension $1+d$ tangent to the parabolic directions [1812.01286]. The corresponding semiconjugacy
\[
F\circ K=K\circ R
\]
is solved by Fourier-Taylor recursion and an a-posteriori theorem, and the resulting whiskered parabolic tori yield solutions tending to parabolic escape while the remaining bodies execute bounded motion [1812.01286].

## 6. Celestial mechanics, transport, and large-scale dynamics

In celestial mechanics, whiskered tori function as organizing centers for global transport. In periodically perturbed PCRTBP and related RTBP models, unstable periodic orbits of the autonomous problem persist as whiskered tori, and intersections between their stable and unstable manifolds create heteroclinic pathways that enable spacecraft to greatly modify their orbits without using propellant [2109.14814]. High-order Fourier-Taylor parameterizations substantially enlarge the usable local manifold domain: in the Jupiter-Europa planar elliptic RTBP, degree-$5$ expansions produced fundamental domains $50$–$200\times$ larger than linear domains [2105.11100].

The computational literature has emphasized both efficiency and scalability. For torus and manifold computation in periodically perturbed RTBP, each quasi-Newton iteration has cost $O(N\log N)$ for $N$ Fourier modes, whereas prior collocation-based methods cost $O(N^3)$ per step [2105.11100]. For heteroclinic search in the Jupiter-Europa planar elliptic RTBP, GPU implementation in Julia and OpenCL reduced the mean kernel time to $\approx 0.03\,\mathrm{s}$ and total time to $\approx 10\,\mathrm{s}$ for checking $112$ half-layer pairs, compared with mean kernel time $\approx 0.49\,\mathrm{s}$ and total time $\approx 61\,\mathrm{s}$ in a CPU-only version; the reported speedup was $\approx 16\times$ in kernel and $\approx 5.8\times$ overall on a modern laptop, with older hardware showing $5$–$7\times$ program speedup [2109.14814].

In the planetary three-body problem, whiskered tori coexist with maximal elliptic KAM tori in the same phase-space region. One rigorous result establishes a Cantor-like family of real-analytic, $3$-dimensional whiskered KAM tori, each with $4$-dimensional stable and unstable manifolds, coexisting with real-analytic, $4$-dimensional maximal KAM tori near the outer-retrograde-coplanar equilibrium [1611.00508]. A later quantitative KAM approach proves continuation of maximal and whiskered quasi-periodic motions in properly degenerate systems and applies it to the planar three-body problem, yielding co-existing families of $5$-dimensional maximal tori and $3$-dimensional whiskered tori with $2$-dimensional stable and unstable manifolds [2302.00279].

Recent work places whiskered tori within a broader theory of robust homoclinic complexity. For any $C^s$ symplectic diffeomorphism having a one-dimensional whiskered torus with a homoclinic orbit, an arbitrarily $C^s$-small perturbation can create a symplectic blender; this, in turn, implies persistence phenomena for saddle-center homoclinics and extends to the corresponding continuous-time settings [2603.20830]. A plausible implication is that whiskered tori should be viewed not only as KAM remnants or transport channels, but also as seeds for higher-dimensional hyperbolic mechanisms that survive perturbation in a robust manner.

Source: https://www.emergentmind.com/topics/whiskered-tori