Alphabeta Math
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Stability, perturbations and chaos — Examples

1 · Prerequisites

2 · Summary

These examples exercise the stability, linearisation, chaos and recurrence items of cm16-stability-perturbations-and-chaos on declared model data and probe the boundaries of the page's statements.

The ideal pendulum on declared data: strict minimum and Lyapunov stability at the downward equilibrium, negative Hessian and instability at the upward equilibrium, the exact linearised eigenvalues and oscillation frequency, and failure of attraction works the ideal pendulum in the angle chart with m=1 kg, ℓ=1 m and g=9.81 m s−2: the downward equilibrium is a strict minimum, hence Lyapunov stable and not attracting, with the purely imaginary linearisation frequencies ±i 3109/10 s−1; the upward equilibrium has the negative Hessian −mgℓ and is unstable, with the real rates ±3109/10 s−1. Two declared linear models worked through: the harmonic flow has zero finite-time growth in every direction, while the repulsive linear model has exact rates ±1 and the positive Lyapunov exponent 1 at its equilibrium computes the flows and growth rates of two dimensionless one-label Cartesian models in R3: the harmonic force F(q)=−q has an orthogonal rotation flow, so every finite-time growth rate vanishes and every existing Lyapunov exponent is 0; the repulsive force F(q)=q has exact rates ±1 along its eigendirections and positive Lyapunov exponent 1 at equilibrium. The three-dimensional isotropic oscillator: the stroboscopic map at half period is minus the identity, every state is periodic, the periodic points are dense, and the map is neither mixing nor sensitive shows the three-dimensional isotropic force model's stroboscopic map at half period to be minus the identity on its six-dimensional phase space: every state is periodic, the periodic points are dense, but the map is neither topologically mixing nor sensitive, so continuous maps with dense periodic points need not be chaotic in the page's stated sense. The pulsed horseshoe on declared data: a=100, the radicand box [74,126], the contraction factor 1/74, the limiting quadratic fixed points 1±101, and the physical bounds of the coded states at ℓ0=1 m, m=1 kg, ω0=2π s−1 evaluates the constants of the pulsed horseshoe on its invariant axis at a=100: the coding box R=13, the radicands in [74,126], the contraction factor 1/74, the limiting quadratic fixed points 1±101, the limiting kick composition, and the axial physical bounds 13 m and 26π kg m s−1 at ℓ0=1 m, m=1 kg, ω0=2π s−1.

The counterexamples delimit the statements. A positive Lyapunov exponent at a hyperbolic equilibrium does not produce a chaotic invariant set of the time-one map: the three-dimensional repulsive model has rate 1 and no compact invariant set with more than one point exhibits the repulsive linear model, whose equilibrium has a positive Lyapunov exponent and whose time-one map has no nonempty compact invariant set beyond {0}; that singleton fails sensitivity, so a positive exponent at an equilibrium does not produce a chaotic invariant set of that map. Lyapunov stability does not imply attraction: the harmonic model is stable at its equilibrium while every nonzero datum keeps its distance and never converges exhibits the harmonic model, Lyapunov stable at its equilibrium while the datum (δ,0) keeps its distance: stability is not attraction. A trapped one-label model: the compact energy sublevel of the radial quartic potential is invariant with finite phase volume, so the recurrence theorem applies, while the free model's energy sublevels are unbounded and the hypothesis fails verifies the recurrence theorem on the trapped one-label, three-dimensional radial quartic model of the flat Cartesian presentation, whose compact invariant energy sublevels have finite six-volume, and contrasts it with the free model, whose unbounded energy sublevels violate the finite-volume hypothesis. Unit restoration and scale checks for the two-rotor KAM statement: the scaling identities at I=E0=1, the physical angular frequencies, the exclusion of rational frequency labels, and the equal-rescaling invariance restores units for the two-rotor KAM statement at I=1 kg m2, E0=1 J, checks the exact scaling identities, the physical angular frequencies ω0νi, the exclusion of rational frequency labels and the invariance under equal rescaling of inertia and energy; it does not evaluate the KAM constant c. All data on this page are declared model data; no measurement or empirical comparison is asserted anywhere.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

None yet.

5 · Examples, experiments and counterexamples

ExampleConstruction: AI-adaptedVerification: AI-adaptedOpen item page →

The ideal pendulum on declared data: strict minimum and Lyapunov stability at the downward equilibrium, negative Hessian and instability at the upward equilibrium, the exact linearised eigenvalues and oscillation frequency, and failure of attraction

Physical scope

A declared-data check of the page's energy-stability theorem on the ideal planar pendulum: with mass 1 kg, length 1 m and g=9.81 m s^-2 in the angle chart, the downward equilibrium is a strict minimum, hence Lyapunov stable but not attracting, with purely imaginary linearised eigenvalues and a positive angular oscillation frequency; the upward equilibrium has a strictly negative Hessian, hence is unstable, with real linearised eigenvalues. All data are declared model data; no measurement and no empirical comparison is made.

Example

Declared model data. Take one label on the real line as the angle chart of the ideal planar pendulum: the chart domain is Ω:=R, the configuration coordinate is the angle q (radians, identified modulo 2π at the end), the kinetic metric is the constant matrix M:=mℓ2 with m:=1 kg, ℓ:=1 m, and the potential is U(q):=mgℓ (1−cos⁡q),g:=9.81 m s−2,mgℓ=9.81 J. This is the mechanical chart datum of The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart with Ω=R, M constant and U of class C∞. Its equations are mℓ2q¨=−mgℓsin⁡q, equivalently q¨=−(g/ℓ)sin⁡q.

Assertions to check. (i) q=0 and q=π are mechanical equilibria; (ii) at q=0 the potential is a strict local minimum and Energy stability and instability of the mechanical model at an equilibrium: the strict minimum is Lyapunov stable, a negative Hessian eigenvalue gives linear instability, a negative definite potential Hessian gives nonlinear instability, and exact energy conservation rules out attraction clause 1 makes (0,0) Lyapunov stable, while the linearised eigenvalues are ±ig/ℓ=±i 3109/10 s−1 (so the angular oscillation frequency is 3109/10 s−1) and the equilibrium is not asymptotically stable; (iii) at q=π the Hessian is −mgℓ<0 and clause 3 makes (π,0) unstable, with linearised eigenvalues ±3109/10 s−1.

Facts & Assumptions

Given: the declared chart data M=1 kg m2, U(q)=9.81(1−cos⁡q) J on Ω=R and the two equilibria q=0, q=π.

[L1]

The energy-stability theorem (Energy stability and instability of the mechanical model at an equilibrium: the strict minimum is Lyapunov stable, a negative Hessian eigenvalue gives linear instability, a negative definite potential Hessian gives nonlinear instability, and exact energy conservation rules out attraction): under the chart hypotheses, a strict local minimum of U makes (q0,0) Lyapunov stable; a negative Hessian eigenvalue gives an unbounded solution eκt(w,κw) of the associated vector equation Y′=DG(q0,0)Y, the action on an initial deviation of the normalized matrix equation (linear instability); a negative definite Hessian D2U(q0) makes the equilibrium unstable; and exact energy conservation prevents asymptotic stability whenever states of different energy accumulate at the equilibrium.

[L2]

Trigonometric facts: (sin⁡)′=cos⁡, (cos⁡)′=−sin⁡, sin⁡0=0 and cos⁡0=1 (The derivatives of sine and cosine are cosine and minus sine); sin⁡π=0, cos⁡π=−1 and the half-turn shifts are those of Quarter-turn values and shifts by pi/2 and pi. Both functions are 2π-periodic, and sin⁡q=0 exactly when q∈πZ (The zero sets of sine and cosine and the least positive common period 2 pi). Cosine is strictly decreasing on [2jπ,(2j+1)π] and strictly increasing on [(2j+1)π,(2j+2)π] for every integer j (Signs, monotonicity intervals, and ranges of sine and cosine); with its value 1 at −2π,0,2π, these inequalities give cos⁡q<1 whenever 0<∣q∣<2π.

Verification

technique · direct
1.1givenL2algebra

Equilibria and second derivative: by [L2], U′(q)=mgℓsin⁡q and U′′(q)=mgℓcos⁡q; hence U′(0)=U′(π)=0, so q=0 and q=π are mechanical equilibria (The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart), with U′′(0)=mgℓ=981/100 J>0 and U′′(π)=−mgℓ<0.

2.1step 1.1givenL1L2algebra

Downward equilibrium: U(q)−U(0)=mgℓ(1−cos⁡q)>0 for 0<∣q∣<2π, because cos⁡q<1 there; thus q=0 is a strict local minimum (take r0:=π), and [L1] first clause gives Lyapunov stability of (0,0) in the chart. Since g/ℓ=981/100 s−2, the linearised first-order matrix is A0=(01−g/ℓ0), whose characteristic polynomial is λ2+g/ℓ and whose eigenvalues are ±ig/ℓ=±i 3109/10 s−1. The coefficient ω:=g/ℓ in η¨+ω2η=0 is the angular oscillation frequency of this linearised model, so the equilibrium is not hyperbolic and no exponential conclusion follows.

2.2step 1.1givenL1algebra

Upward equilibrium: U′′(π)=−981/100 J<0, and since n=1 the Hessian is negative definite at q=π; [L1] third clause therefore makes (π,0) unstable. The linearised first-order matrix is Aπ=(01g/ℓ0), with characteristic polynomial λ2−g/ℓ and real eigenvalues ±g/ℓ=±3109/10 s−1, one positive, matching the linear-instability statement of [L1] second clause.

3.1step 2.1givenL1algebra

Failure of attraction: for every 0<ε<π the datum z0:=(ε,0) has energy E(z0)=mgℓ(1−cos⁡ε)>0=E(0,0), so letting ε↓0 supplies states of different energy accumulating at the equilibrium; by [L1] fourth clause (and the strict-minimum clause of step 2.1) (0,0) is Lyapunov stable but not asymptotically stable. In particular no datum of positive energy can converge to (0,0), because the energy is conserved.

4.1step 2.1step 2.2givenalgebra∎

All quantities above are declared model data substituted into the chart formulas; M has dimension kg m2, the potential and the energies share the dimension J=kg m2 s−2, and the frequencies have dimension s−1 (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework). No pendulum, material or measurement is involved.

Remarks

ExampleConstruction: AI-adaptedVerification: AI-adaptedOpen item page →

Two declared linear models worked through: the harmonic flow has zero finite-time growth in every direction, while the repulsive linear model has exact rates ±1 and the positive Lyapunov exponent 1 at its equilibrium

Physical scope

A declared-data check of two physical one-label Cartesian force models with m0=1 kg, ell0=1 m and t0=1 s: F_1,phys=-(m0/t0^2)X and F_2,phys=(m0/t0^2)X, in newtons. The formulas use normalized q=X/ell0, v=t0 V/ell0 and t=T/t0, mass 1 and potentials U_1(q)=||q||^2/2 and U_2(q)=-||q||^2/2 with energy reference E0=m0 ell0^2/t0^2. In this numerical state metric the harmonic derivative is orthogonal, so all finite-time rates and Lyapunov exponents are zero; the repulsive derivative has exact rates plus and minus 1 along its eigendirections, with positive exponent 1 at equilibrium. The corresponding physical rates divide by t0. All data are declared, with no measurement and no claim of chaos.

Example

Declared model data and normalization. Take one label with Q:=R3 and declare physical reference data m0:=1 kg, ℓ0:=1 m, t0:=1 s. For σ1=−1, σ2=1 define physical force maps Fi,phys:R×Q×R3→R3 by Fi,phys(T,X,V)=σi(m0/t02)X, with T in seconds, X in metres, V in metres per second, and outputs in newtons. These are smooth globally state-Lipschitz force maps of Adopted force model: force maps, the equation of motion, additive composition of force contributions and frame covariance with physical mass m0. Put t=T/t0, q=X/ℓ0, v=t0V/ℓ0 and E0=m0ℓ02/t02. Then m0X¨=Fi,phys becomes q′′=σiq, with q′=v, since X¨=(ℓ0/t02)q′′. The normalized force maps are Fi(t,q,v)=Fi,phys(t0t,ℓ0q,(ℓ0/t0)v)/(m0ℓ0/t02)=σiq, with normalized mass m=1. The physical potentials are Ui,phys(X)=−σim0∥X∥22/(2t02); normalizing their energies by E0 gives, in the constant-metric chart Ω:=R3 with M:=I3, they have the potentials U1(q):=12∥q∥22 and U2(q):=−12∥q∥22 of The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart. The associated fields and complete motion maps are listed below; e1,e2,e3 denote the coordinate basis, explicitly ej(k)=1 for k=j−1 and 0 otherwise, for k∈{0,1,2} and j∈{1,2,3}. All state and tangent norms below apply to the dimensionless numerical coordinates (q,v).

  1. The harmonic model has field G1(q,v)=(v,−q) and motion map Φt1(q0,v0)=(q0cos⁡t+v0sin⁡t, −q0sin⁡t+v0cos⁡t).
  2. The repulsive linear model has field G2(q,v)=(v,q) and motion map Φt2(q0,v0)=(12(q0+v0)et+12(q0−v0)e−t, 12(q0+v0)et−12(q0−v0)e−t).

For a state a and a nonzero tangent vector u with DaΦt(a)u≠0, the finite-time tangent growth rate is t−1log⁡(∥DaΦt(a)u∥2/∥u∥2); when the forward lifetime is unbounded and its finite limit as t→∞ exists, that limit is the Lyapunov exponent (Finite-time tangent growth rates and Lyapunov exponents of the model's flow: the rate per unit time of a transported deviation, with its existence convention and its limits clauses (i)-(ii)). Both displayed flows are complete, so this lifetime condition holds for each state.

Assertions to check. (i) Both displayed maps are the local flows of their fields; (ii) for the harmonic model DaΦt1(a) is the orthogonal block rotation (cos⁡t I3sin⁡t I3−sin⁡t I3cos⁡t I3), so every finite-time tangent growth rate vanishes and every existing Lyapunov exponent is 0; (iii) for the repulsive model the rates along the directions (e1,−e1) and (e1,e1) are −1 and +1, so the Lyapunov exponent of the equilibrium along the unstable eigenvector is 1.

Facts & Assumptions

Given: the two force maps and their equivalent chart data with M=I3, the displayed candidate flows and the finite-time growth-rate definition of Finite-time tangent growth rates and Lyapunov exponents of a C¹ flow as recorded in Finite-time tangent growth rates and Lyapunov exponents of the model's flow: the rate per unit time of a transported deviation, with its existence convention and its limits.

[L1]

Uniqueness of force-model motions (Local existence and uniqueness for the Newtonian initial-value problem, with the quantitative Picard-Lindelöf cylinder): for the adopted physical force maps, every initial datum has a unique maximal model motion; consequently a globally defined curve that solves the equation with the given initial data is the model's maximal motion. The invertible state and time normalization transfers this uniqueness and maximality to the numerical systems.

[L3]

Pythagorean identity (Parity and the Pythagorean identity for sine and cosine): for every real t, sin⁡2t+cos⁡2t=1.

Verification

technique · direct
1.1givenalgebra

The two fields: for the constant metric M=I3 the mechanical field of The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart is G(q,v)=(v,−∇U(q)), giving G1(q,v)=(v,−q) and G2(q,v)=(v,q). These are also the normalized force-model fields of the declared physical F1,phys,F2,phys.

2.1givenstep 1.1L1L2algebra

Harmonic flow: differentiating the displayed Φt1 with respect to t by [L2] gives (q˙,v˙)=(v,−q); at t=0 the map is the identity. By [L1] the displayed map is the force model's maximal flow.

2.2givenstep 1.1L1L2algebra

Repulsive flow: differentiating the displayed Φt2 gives q˙=v(t) and v˙=q(t); at t=0 the map is the identity, so by [L1] it is the force model's maximal flow of G2.

3.1step 2.1givenL2L3algebra

Harmonic growth rates: Φt1 is linear in a=(q0,v0), with derivative DaΦt1(a)=(cos⁡t I3sin⁡t I3−sin⁡t I3cos⁡t I3); writing c=cos⁡t, s=sin⁡t, its transpose times itself has diagonal blocks (c2+s2)I3=I3 by [L3] and off-diagonal blocks csI3−scI3=0, so it is orthogonal, whence ∥DaΦt1(a)u∥2=∥u∥2 for every u≠0. The explicit flow exists for every real time, so the forward lifetime is unbounded, as required by the exponent definition. Therefore r(t,a,u)=t−1log⁡1=0 for every t>0, and the limit defining the Lyapunov exponent exists and equals 0; the variational equation of Linearisation along a motion of the mechanical model: the flow derivative solves the variational equation, with Lipschitz dependence, the constant-metric second-order form, and the conditional finite-time error bound of the linear approximation is satisfied by this derivative matrix, whose columns solve Y˙=DG1Y by step 2.1.

3.2step 2.2givenL2algebra

Repulsive growth rates: for an initial tangent vector (δq,δv) put u:=δq+δv and w:=δq−δv. The derivative flow satisfies δq(t)+δv(t)=etu and δq(t)−δv(t)=e−tw. For each j∈{1,2,3}, input (ej,ej) gives output et(ej,ej) and rate +1, while input (ej,−ej) gives output e−t(ej,−ej) and rate −1. The explicit flow exists for every real time, so the lifetime hypothesis holds and these Lyapunov exponents exist. At the equilibrium a=(0,0) this gives a positive exponent along each unstable coordinate eigenvector.

4.1step 3.1step 3.2givenL1L2algebra∎

Scope check: the chain rule transfers each numerical motion to X(T)=ℓ0q(T/t0) satisfying its physical Newton equation. A rate r per numerical time t becomes r/t0 per physical time T; thus the physical rates are 0 and ±1/t0, and the positive physical exponent is 1/t0. All steps substitute declared normalized data into the chart formulas and the finite-time-rate definition; no step produces a chaotic invariant set, a mixing property or a measurement, and the positive exponent of step 3.2 concerns the linear model at its equilibrium, as the page's counterexample discusses.

Remarks

  • What is genuinely checked. The example verifies that the variational equation holds with the displayed derivative matrices and that the growth rate is invariant under multiplication of the deviation by nonzero scalars (the logarithm of a positive factor cancels against the norm ratio), so the rate is a function of the direction, not of the length of u. Both flows are complete on the declared chart; no domain-escape phenomenon occurs.

  • The unstable direction. For the repulsive model the directions (ej,ej) and (ej,−ej) are exactly the eigendirections of A=DG2(0,0)=(0I3I30), with eigenvalues +1 and −1, respectively, so the rates ±1 are the eigenvalues of the linearisation, in agreement with Hyperbolic linearisation of an equilibrium of the model: exponential dichotomy of the linearised flow, local exponential asymptotic stability in the stable case, and instability in the saddle case clause 1 with c=1.

ExampleConstruction: AI-adaptedVerification: AI-adaptedOpen item page →

The three-dimensional isotropic oscillator: the stroboscopic map at half period is minus the identity, every state is periodic, the periodic points are dense, and the map is neither mixing nor sensitive

Physical scope

A declared-data check on the physical one-label force model with mass m0=1 kg, force F_phys(T,X,V)=-(m0/t0^2)X and potential U_phys(X)=m0||X||^2/(2t0^2), using ell0=1 m and t0=1 s. The formulas use numerical coordinates q=X/ell0, v=t0 V/ell0 and t=T/t0, normalized mass 1 and potential ||q||^2/2. The stroboscopic map at numerical half-period pi is minus the identity and at 2pi is the identity; the first-return map to q_y=0, v_y>0 is the identity. All states return, including the constant zero state. Periodic points of the map are dense, but the map is neither mixing nor sensitive and is not chaotic in the stated sense. Physical periods multiply by t0. All parameters are declared model data; no measurement is reported.

Example

Declared model data and normalization. Fix one label, Q:=R3, m0:=1 kg, ℓ0:=1 m and t0:=1 s. The physical force map is Fphys:R×Q×R3→R3, Fphys(T,X,V)=−(m0/t02)X, in newtons, with T in seconds, X in metres and V in metres per second. This smooth globally state-Lipschitz map and mass m0 specify a force model of Adopted force model: force maps, the equation of motion, additive composition of force contributions and frame covariance. Put t=T/t0, q=X/ℓ0 and v=t0V/ℓ0. Newton's equation becomes q′′=−q since X¨=(ℓ0/t02)q′′. The normalized force is F(t,q,v):=Fphys(t0t,ℓ0q,(ℓ0/t0)v)/(m0ℓ0/t02)=−q; m=1 is its normalized mass. With Uphys(X)=m0∥X∥22/(2t02) and E0=m0ℓ02/t02, U=Uphys/E0 gives, in the chart Ω:=R3 the constant kinetic metric is M:=I3 and the potential is U(q):=12∥q∥22 of The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart. The field and motion through (q0,v0) are G(q,v)=(v,−q),q(t)=q0cos⁡t+v0sin⁡t,v(t)=−q0sin⁡t+v0cos⁡t. Regard the autonomous force model as P-periodic for any P>0, and let ΨP be its stroboscopic map of The stroboscopic map of a periodic model and the Poincare return map of a transversal section: domains, iterates, fixed points and the local C1 diffeomorphism property clause (ii). The complete flow makes each ΨP defined on all of R6. Numerical period P means physical period t0P; all state norms below use the dimensionless coordinates (q,v). Use the section Σ+:={(q,v):qy=0, vy>0}; its first-return map sends each section state to the state at its least positive return time, with the local C1 construction supplied by The stroboscopic map of a periodic model and the Poincare return map of a transversal section: domains, iterates, fixed points and the local C1 diffeomorphism property clause (iii).

Assertions to check. (i) The displayed curves are the model's motions; (ii) Ψπ=−id and Ψ2π=id, so every state of the model is periodic under the flow with period 2π and every point is a periodic point of Ψπ with period dividing 2 (the zero state has a constant motion; only nonzero states give nonconstant periodic motions); (iii) the first-return map of Σ+ is the identity; (iv) with the Euclidean metric of the phase chart, Ψπ on X:=R6 is continuous and has dense periodic points, but is not topologically mixing and has no sensitivity constant, hence is not chaotic in the stated sense of Chaos in the stated sense for maps of the model: topological mixing, dense periodic points and sensitive dependence, as one precise list among several.

Facts & Assumptions

Given: the declared one-label force model Q=R3, m=1, F(T,q,v)=−q, its equivalent chart data M=I3, U(q)=∥q∥22/2, the field G and the displayed motion formulas.

[L1]

Unique force-model motions (Local existence and uniqueness for the Newtonian initial-value problem, with the quantitative Picard-Lindelöf cylinder): for the adopted physical smooth force map, every initial datum has a unique maximal motion, and every other solution through that datum is its restriction. The invertible state and time normalization above transfers existence, uniqueness and maximality to the displayed numerical system.

[L3]

The chaos list (Topological mixing, dense periodic points, sensitive dependence, and the stated chaos list): a map is chaotic in the stated sense exactly when it is continuous, topologically mixing, has dense periodic points and has sensitive dependence with a positive constant; the negation of topological mixing or of sensitivity therefore suffices to fail the list.

[L4]

The local first-return proposition (A transversal return of a C¹ flow yields a local C¹ Poincaré map, with inverse given by the reversed crossing): for a C1 autonomous field and a C1 section, two transversal crossings define a local return map, and exclusion of intervening hits from the section closure makes it a first-return map.

[L5]

Trigonometric values and zeros (Quarter-turn values and shifts by pi/2 and pi, The zero sets of sine and cosine and the least positive common period 2 pi): sin⁡π=0, cos⁡π=−1, and sine and cosine have period 2π, so sin⁡(2π)=sin⁡0=0, cos⁡(2π)=cos⁡0=1 using [L2]. For every real t, sin⁡t=0 exactly when t=kπ for an integer k. In particular, the only zero in 0<t<2π is t=π.

Verification

technique · direct
1.1givenL1L2algebra

Field and motion: with M=I3 and U(q)=∥q∥22/2 the mechanical field is G(q,v)=(v,−q) (The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart), which is also the normalized field of the declared physical force map. Differentiating the displayed vector formulas gives q˙=v and v˙=−q, and at t=0 they reproduce the initial datum; by [L1] these are the unique maximal model motions.

2.1step 1.1givenL2L5algebra

Stroboscopic maps: evaluating step 1.1 at t=π uses cos⁡π=−1, sin⁡π=0 and gives Ψπ(q,v)=(−q,−v); at t=2π it gives the identity. Thus every state returns after 2π, and Ψπ is an involution, so every point of X=R6 is a periodic point of Ψπ with period 1 or 2; in particular the periodic points are all of X and are dense.

2.2step 1.1givenL2L4L5algebra

First return to Σ+: take the open state neighbourhood W:={(q,v):vy>0} and the section function h:W→R, h(q,v):=qy, so h−1(0)=Σ+ within W. At every point of Σ+, Dh G=vy>0. A state with qy0=0 and vy0>0 has qy(t)=vy0sin⁡t and vy(t)=vy0cos⁡t, so the next hit of Σ+ is at t=2π; at that time the full six-dimensional state returns to its initial value by the formulas of step 1.1. The open segment 0<t<2π avoids Σ+‾={qy=0,vy≥0}: by [L5] its only intermediate zero is at t=π, where vy=−vy0<0. Applying [L4] at each reference point gives the local first-return map, and the explicit calculation identifies it with the identity on all of Σ+.

3.1step 2.1L3algebra

Failure of mixing: let Z∈X=R6 be a state with ∥Z∥2>1 and let U:=B(Z,1) be the open Euclidean ball (Open ball, closed ball and sphere in a metric space); then U∩(−U)=∅, and step 2.1 gives Ψπn(U)=U for even n and Ψπn(U)=−U for odd n. Taking V:=U in the mixing condition of [L3], the intersection Ψπn(U)∩V is empty for every odd n; consequently no integer N has Ψπn(U)∩V≠∅ for all n≥N, and Ψπ is not topologically mixing on X.

4.1step 2.1step 3.1L3algebra

Failure of sensitivity: for all states Z,Z′∈R6 and all n≥1, step 2.1 gives ∥Ψπn(Z)−Ψπn(Z′)∥2=∥Z−Z′∥2 because each iterate is ±id, an isometry; hence for no δ>0 can one find, for every Z and every ball around Z, a nearby state whose iterate is δ-separated. The sensitivity condition of [L3] fails, and with it the chaos list; the periodic-point and continuity clauses alone do not make Ψπ chaotic in the stated sense.

5.1step 2.2step 4.1givenL1L2algebra∎

Reading: the maps use the dimensionless numerical state and time variables defined above. By the chain rule, X(T)=ℓ0q(T/t0) obeys m0X¨=−(m0/t02)X; uniqueness transfers through this invertible normalization. The physical first-return time is 2πt0 and the physical half-period is πt0. No measurement, sampling or empirical comparison is performed.

Remarks

ExampleConstruction: AI-adaptedVerification: AI-adaptedOpen item page →

The pulsed horseshoe on declared data: a=100, the radicand box [74,126], the contraction factor 1/74, the limiting quadratic fixed points 1±101, and the physical bounds of the coded states at ℓ0=1 m, m=1 kg, ω0=2π s−1

Physical scope

A declared-data check of the constants and bounds in the page's pulsed-horseshoe theorem: at amplitude a=100 the coding box is R=13, the radicands lie in [74,126], the contraction factor is 1/sqrt(74), the limiting quadratic map has fixed points 1 plus or minus sqrt(101), the limiting kick satisfies K composed with the rotation = F, and with length scale 1 m, mass 1 kg and frequency scale 2 pi per second the coded states on the invariant axis are bounded by 13 m in axial position and 26 pi kg m per second in axial momentum, with period (pi/2+delta)/(2 pi) seconds. All quantities are declared model data; the pulse-width threshold is not numerically evaluated and no measurement is involved.

Example

Declared model data. Take the full physical pulsed oscillator model of A smooth periodically forced one-label Cartesian model with a transverse homoclinic horseshoe: the physical stroboscopic map has a compact invariant set conjugate to the shift, chaotic in the stated sense clause (i), and restrict its coded states to the invariant axis, where transverse coordinates and momenta vanish. Use the declared data a:=100,ℓ0:=1 m,m:=1 kg,ω0:=2π s−1, and a pulse width 0<δ<δ0 with δ0 the supplier threshold (not evaluated numerically here). Let R:=a+3, L:=a−2R and let r± be the fixed points of the limiting quadratic map.

Assertions to check. (i) R=13 and the coding radicands a+xn−1+xn+1 lie in [74,126]=[a−2R,a+2R], with 126<13=R; (ii) the contraction factor of the coding map is 1/L=1/74 and r±=1±101; (iii) for the limiting pulse, K∘R0=F with R0(q,p)=(p,−q), K(q,p)=(q,p+q2−a) and F(x,y)=(y,y2−a−x); (iv) every coded state satisfies ∣Q∣≤13 m and ∣Π∣≤26π kg m s−1, and the forcing period is Tp=(π/2+δ)/(2π) s.

Facts & Assumptions

Given: the constants a=100, R=a+3, L=a−2R and the scales ℓ0,m,ω0 above; the dimensionless model, its rotation R0, pulse endpoint maps Kδ and stroboscopic map Gδ=Kδ∘R0 of The smooth P-periodic pulsed oscillator: the normalized pulse, its Hamiltonian and its stroboscopic (period) map; the limiting kick K.

[L1]

The horseshoe theorem and its bounds (A smooth periodically forced one-label Cartesian model with a transverse homoclinic horseshoe: the physical stroboscopic map has a compact invariant set conjugate to the shift, chaotic in the stated sense clauses (ii)-(iv)): for 0<δ<δ0 the full physical stroboscopic map Ψδ has a compact invariant coded set on the invariant axis, with transverse coordinates and momenta zero; each state satisfies ∣Q∣≤Rℓ0 and ∣Π∣≤mℓ0ω0R, and its coded trajectories exist for all times.

[L2]

The coding contraction (The square-root recurrence Ts(x)n=sna+xn−1+xn+1 is a contraction of BR with a unique coded solution) and the limiting pulse (The rescaled pulse endpoint map converges to the kick K(q,p)=(q,p+q2−a) in C1, and the period map is smooth, symplectic and C1-close to F): the coding map Tσ on the ball of bounded sequences with ∣x∣∞≤R is a contraction with factor L−1 when the radicands stay positive and the derivative bound holds, and Kδ→K in C1 as δ→0 with K(q,p)=(q,p+q2−a).

[L3]

Elementary arithmetic: R2−a−2R=4a+3>0 for a=100, so R=a+3 satisfies a+2R<R; the fixed points of x↦ the constant-code recurrence x2−a=2x solve x2−2x−a=0.

Verification

technique · direct
1.1givenL2L3algebra

Arithmetic of the coding box: for a=100, R=100+3=13; hence a−2R=100−26=74 and a+2R=100+26=126, so the radicands of [L2] lie in [74,126] whenever the coded sequence satisfies ∣xn∣≤R. Moreover 126<13 because 126<169, and 74>0, so all square roots in the coding recurrence are defined with the derivative bound required by [L2]. Also R2−a−2R=169−100−26=43=4100+3>0, which is the supplier's inequality ensuring Tσ maps the ball into itself.

1.2givenL2algebra

Limiting pulse composition: with R0(q,p)=(p,−q) and K(q,p)=(q,p+q2−a) one computes K(R0(q,p))=K(p,−q)=(p,−q+p2−a) and F(q,p)=(p,p2−a−q), which agree; this is the algebraic core of the C1-closeness statement Gδ→F of [L2].

2.1step 1.1givenL2L3algebra

Contraction and fixed points: the radicand bound of step 1.1 gives the derivative bound 1274 for the square root, hence the contraction factor (1274)⋅2=174=L−1 for the coding map, as in [L2] with L=74. Constant codes satisfy xn+1+xn−1=xn2−a with xn≡r, that is 2r=r2−a or r2−2r−a=0, whose solutions are r=1±1+a=1±101. The branch conditions hold with exact inequalities: 101>8 because 101>64, so ∣r−∣=101−1>7; 101>1+74 because 101>75+274 reduces to 13>74, that is 169>74; so ∣r±∣>L; and r+≤13, ∣r−∣≤13 because 101≤12 and 101−1≤13, both from 101≤144.

2.2step 1.1givenL1algebra

Physical bounds: by [L1] clause (iv) every coded physical state satisfies ∣Q∣≤Rℓ0=13 m and ∣Π∣≤mℓ0ω0R=13⋅2π kg m s−1=26π kg m s−1; the forcing period is Tp=P/ω0=(π/2+δ)/(2π) s and the pulse duration is δ/(2π) s.

3.1step 2.1step 1.2step 2.2givenalgebra∎

Reading: every number above is declared model data or an exact surd; the threshold δ0 is not numerically evaluated, no decimal approximation is substituted into a proof step, and no measurement or apparatus is involved (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework).

Remarks

  • Status of the numerical values. The values 13, 74, 126, 1/74, 1±101 and 26π are exact consequences of the declared data; they are quoted as a check of the supplier's constants, not as an empirical estimate of any oscillation. In particular the threshold δ below which the horseshoe exists is left symbolic because the supplier does not compute it.

  • Unit consistency. [ℓ0]=ℓ, [Π]=mℓ0ω0=μℓτ−1, [Tp]=τ, and the coded box is bounded in both position and canonical momentum (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework).

CounterexampleConstruction: AI-adaptedVerification: AI-adaptedOpen item page →

A positive Lyapunov exponent at a hyperbolic equilibrium does not produce a chaotic invariant set of the time-one map: the three-dimensional repulsive model has rate 1 and no compact invariant set with more than one point

Physical scope

A declared-model counterexample in the one-label three-dimensional force branch, refuting the implication from a positive Lyapunov exponent at an equilibrium to a chaotic compact invariant set of the time-one map: the dimensionless repulsive model with potential U(q)=-||q||^2/2 has a complete flow, exponent 1 at its equilibrium along each unstable direction, and a nonlinearly unstable equilibrium, yet every nonempty compact invariant set of its time-one map is {0}; this singleton fails the sensitivity clause of the page's chaos list. The counterexample is a model computation; it reports no measurement.

Statement refuted

False statement. Let G:S0→R6N be the C1 associated field of a complete autonomous stability model in the sense of The model vocabulary of stability and chaos: the associated flow of a declared force model, equilibria, periodic and homoclinic orbits, and invariant sets, with flow Φ:R×S0→S0. If an equilibrium has a positive real eigenvalue λ of DG and Lyapunov exponent λ along a corresponding eigenvector, then the time-one map f:=Φ1:S0→S0 possesses a compact invariant set with more than one point on which its restriction is chaotic in the stated sense of Chaos in the stated sense for maps of the model: topological mixing, dense periodic points and sensitive dependence, as one precise list among several (continuous, topologically mixing, dense periodic points, sensitive dependence).

Facts & Assumptions

Given: the dimensionless one-label Cartesian force model with Q=R3, mass m=1, kinetic matrix M=I3, potential U(q)=−12∥q∥22, force F(q,v)=q, field G(q,v)=(v,q), equilibrium (0,0), and complete flow Φt(q,v)=(12(q+v)et+12(q−v)e−t, 12(q+v)et−12(q−v)e−t).

[L1]

The growth computation (Two declared linear models worked through: the harmonic flow has zero finite-time growth in every direction, while the repulsive linear model has exact rates ±1 and the positive Lyapunov exponent 1 at its equilibrium): for u:=q+v and w:=q−v, the derivative flow acts as (u,w)↦(etu,e−tw); along (e1,e1) the finite-time rate at the equilibrium is 1 for every t>0, so the Lyapunov exponent is 1, and the linearisation has positive eigenvalue 1.

[L2]

The chaos list (Topological mixing, dense periodic points, sensitive dependence, and the stated chaos list): sensitivity requires δ>0 such that for every x and every ε>0 there are y and n≥1 with d(x,y)<ε and d(fn(x),fn(y))≥δ; on a one-point space the only admissible y is x, so no positive δ satisfies the condition.

[L3]

Time-one invariance and instability (The model vocabulary of stability and chaos: the associated flow of a declared force model, equilibria, periodic and homoclinic orbits, and invariant sets clause (iv), Hyperbolic linearisation of an equilibrium of the model: exponential dichotomy of the linearised flow, local exponential asymptotic stability in the stable case, and instability in the saddle case clause 3): for the complete flow, f=Φ1 is a bijection with inverse Φ−1; a compact invariant set of f satisfies f(K)=K, and a hyperbolic equilibrium with a positive-real-part eigenvalue is unstable.

[L4]

Norms (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, Open ball, closed ball and sphere in a metric space): ∥(q,v)∥22=12(∥u∥22+∥w∥22), so boundedness in the state chart implies boundedness of both transformed vectors.

Counterexample

technique · direct
1.1givenL1

The antecedent holds: by [L1] the model has equilibrium (0,0), the linearisation has eigenvector (e1,e1) with eigenvalue 1>0, and the Lyapunov exponent along that direction exists and equals 1. The explicit flow is defined for every real time. This satisfies the antecedent of the refuted statement.

1.2givenL1L3L4algebra

Classification of compact invariant sets of f: let K be compact with f(K)=K, and let a state of K have transformed coordinate u≠0. By [L1], ∥fn(q,v)∥2≥12∥u∥2en for n≥0, which is unbounded; this contradicts boundedness of K. Hence u=0 on K. Since f is a bijection and f(K)=K, applying the same argument to f−n(K)=K shows w=0 on K. Therefore every nonempty compact invariant set of f is the singleton {(0,0)}.

2.1step 1.2L2algebra

The one-point invariant set fails the chaos list: for X={(0,0)} every admissible y in the sensitivity condition of [L2] coincides with the base point, so d(fn(x),fn(y))=0<δ for every n and every δ>0; the sensitivity clause fails, and with it the four-clause list of Chaos in the stated sense for maps of the model: topological mixing, dense periodic points and sensitive dependence, as one precise list among several.

3.1step 1.1step 1.2step 2.1L3∎

Conclusion: the antecedent of the refuted statement holds for the declared model (step 1.1), while its time-one map has no compact invariant set with more than one point (step 1.2) and its only compact invariant set is not chaotic in the stated sense (step 2.1); the implication is therefore false. Nonlinear instability of the equilibrium from [L3] does not repair it: instability of a single state and chaos of a compact invariant set are different properties.

Remarks

CounterexampleConstruction: AI-adaptedVerification: AI-adaptedOpen item page →

Lyapunov stability does not imply attraction: the harmonic model is stable at its equilibrium while every nonzero datum keeps its distance and never converges

Physical scope

A declared-model counterexample refuting the implication from Lyapunov stability to asymptotic stability for the page's mechanical models: the one-dimensional harmonic model with unit mass and unit stiffness has the equilibrium 0, which is Lyapunov stable by the strict-minimum theorem, while the datum (delta,0) has the periodic motion delta cos(t) whose distance from the equilibrium returns to delta at every period, so it never converges. The counterexample is a model computation; it reports no measurement.

Statement refuted

False statement. Every Lyapunov stable equilibrium of a one-dimensional mechanical chart model attracts every sufficiently close initial state; equivalently, Lyapunov stability implies asymptotic stability.

Facts & Assumptions

Given: the dimensionless harmonic chart data: the one-dimensional model of The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart with constant kinetic metric M=1, potential U(q)=12q2, field G(q,v)=(v,−q) and equilibrium (0,0); the comparison data (δ,0) for real δ>0.

[L2]

The harmonic motions (The motions of the harmonic model are the sinusoids Acos⁡(ω(t−t0))+Bsin⁡(ω(t−t0)): amplitude and phase, the least period, the conserved energy kR2/2 and the elliptical phase curves clauses 1-3): for the harmonic model with m=k=1 and ω=k/m=1, every motion is x(t)=Acos⁡(t−t0)+Bsin⁡(t−t0) with A=x(t0), B=v(t0); if R=A2+B2>0 then x(t)=Rcos⁡((t−t0)−φ) and the least positive period is 2π; the datum (x(t0),v(t0))=(δ,0) has A=δ, B=0.

Counterexample

technique · direct
1.1givenL1

The equilibrium is Lyapunov stable: U(q)=12q2 has the strict local minimum q=0 (U(q)>U(0) for q≠0), so [L1] makes (0,0) Lyapunov stable in the chart model.

1.2givenL2algebra

The nonzero datum does not converge: by [L2] the datum (δ,0) at time t0=0 has the maximal motion x(t)=δcos⁡t, with A=δ, B=0, R=δ>0; the least positive period is 2π and x(2πn)=δcos⁡(2πn)=x(0)=δ for every integer n. Hence the sequence tn=2πn tends to infinity while ∥(x(2πn),x′(2πn))−(0,0)∥2≥δ; the state does not converge to the equilibrium.

2.1step 1.1step 1.2given∎

Conclusion: the equilibrium is Lyapunov stable (step 1.1) but for every δ>0 the datum (δ,0), arbitrarily close to the equilibrium, fails to converge (step 1.2); the equilibrium is therefore not asymptotically stable, and the implication of the refuted statement is false for the declared model. The stronger statement of Energy stability and instability of the mechanical model at an equilibrium: the strict minimum is Lyapunov stable, a negative Hessian eigenvalue gives linear instability, a negative definite potential Hessian gives nonlinear instability, and exact energy conservation rules out attraction clause 4 explains the failure: the energy 12(v2+q2) is conserved, equals δ2/2 on this datum, and no state of different energy can converge to the equilibrium.

Remarks

ExampleConstruction: AI-adaptedVerification: AI-adaptedOpen item page →

A trapped one-label model: the compact energy sublevel of the radial quartic potential is invariant with finite phase volume, so the recurrence theorem applies, while the free model's energy sublevels are unbounded and the hypothesis fails

Physical scope

A declared-data check in dimensionless normalized coordinates of the page's recurrence theorem on a trapped one-label, three-dimensional Cartesian model: for the unit-mass radial quartic potential U(x)=||x||^4/4 the six-dimensional phase-space energy sublevel {E<=E1}, E1>0, is compact with finite dimensionless coordinate volume bounded by 512 E1^(9/4), invariant under the flow by energy conservation and compact continuation, so the recurrence theorem gives almost-everywhere infinite returns; in contrast the dimensionless free model U=0 has unbounded energy sublevels and each nonzero-velocity trajectory eventually exits every bounded set, showing the finite-volume hypothesis is not automatic. No SI quantity, measurement or empirical comparison is asserted.

Example

Declared model data. Use dimensionless normalized position, velocity, mass, time and energy coordinates; this is a dimensionless model computation and no SI units are assigned. Take one label in the flat Cartesian presentation of Poincare recurrence for the model's Hamiltonian flow on a bounded invariant set of finite phase volume: almost every state of any Borel subset returns to it infinitely often with N=1: the configuration domain is Q:=R3, the mass is m:=1 so that M=I3 and W=Q×R3 with coordinates (x,v), and the radial quartic potential is U(x):=14∥x∥24, a C∞ datum of The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart with field G(x,v)=(v,−∥x∥22x) and energy E(x,v)=12∥v∥22+14∥x∥24. For E1>0 put AE1:={(x,v)∈W:E(x,v)≤E1}.

Assertions to check. (i) AE1 is a nonempty compact Borel subset of W with λ6(AE1)≤(2(4E1)1/4)3(22E1)3=512 E19/4<∞; (ii) the flow of the model maps AE1 onto itself for every real time, so the hypotheses of Poincare recurrence for the model's Hamiltonian flow on a bounded invariant set of finite phase volume: almost every state of any Borel subset returns to it infinitely often are met and for every Borel E⊆AE1 almost every state of E returns to E infinitely often; (iii) for the flat model U≡0 the energy sublevels are unbounded, the phase measure of W is infinite, and each nonzero-velocity free motion eventually exits every bounded set (hence visits a fixed bounded set at only finitely many integer times), so the finite-volume hypothesis is not automatic.

Facts & Assumptions

Given: the data N=1, M=I3, U(x)=∥x∥24/4, its field and energy; the sublevel AE1; the flat comparison data U≡0 with field G0(x,v)=(v,0).

[L1]

Recurrence (Poincare recurrence for the model's Hamiltonian flow on a bounded invariant set of finite phase volume: almost every state of any Borel subset returns to it infinitely often): for an autonomous potential presentation with a bounded invariant Borel set A whose compact closure lies in the phase domain and whose flow is defined for all real times on A, phase volume is preserved on its Borel subsets by finite local covers and time subdivision; almost every state of every Borel E⊆A returns to E infinitely often at integer times. The theorem assumes ACω, inherited from the phase-volume result.

[L2]

Energy conservation and compact continuation (Energy is conserved along solutions of the mechanical system with a configuration-dependent kinetic metric, Compact continuation for locally Lipschitz ODEs: a solution whose graph stays in a compact subset of the domain near a finite endpoint extends past it): the energy of the mechanical system is constant along every maximal solution, and a solution whose graph stays in a compact subset of the state domain cannot reach a finite forward or backward endpoint.

[L3]

Borel and Lebesgue vocabulary (The Borel sigma-algebra of a topological space, Lebesgue measurable sets, the family L(Rn), and the restricted set function λn): closed sets are Borel, and the Lebesgue measure of a box in R6 is the product of its side lengths.

Verification

technique · direct
1.1givenL3algebra

Compactness and phase volume: E is continuous and AE1=E−1((−∞,E1]) is closed; if (x,v)∈AE1 then 14∥x∥24≤E1 and 12∥v∥22≤E1, so ∥x∥2≤(4E1)1/4 and ∥v∥2≤2E1; hence AE1 is contained in the box of side lengths 2(4E1)1/4 (three position coordinates) and 22E1 (three velocity coordinates), is bounded and closed, therefore compact, and by [L3] has λ6(AE1)≤(2(4E1)1/4)3(22E1)3=512 E19/4<∞; it is nonempty because (0,0)∈AE1.

2.1step 1.1givenL2algebra

Invariance: along a maximal solution of the quartic model the energy is constant by [L2], so the solution stays in AE1; because AE1 is compact and contained in the state domain W, [L2] also shows the solution is defined for all times, and applying the same argument to the reversed motion gives Φt(AE1)⊇AE1 as well; hence Φt(AE1)=AE1 for every t∈R.

3.1step 1.1step 2.1L1

Recurrence conclusion: AE1 is a bounded Borel set of finite measure with closure in the flow domain and is invariant by step 2.1; [L1] therefore gives, for every Borel E⊆AE1, that almost every state of E returns to E at infinitely many positive integer times.

4.1step 1.1givenL1algebra∎

Failure of the hypothesis for the flat model: for U≡0, E0(x,v)=12∥v∥22 and the sublevel {E0≤E1}=R3×{v:∥v∥2≤2E1} is unbounded, indeed it meets every translate of the velocity ball in the position directions; the phase domain W=R6 has infinite Lebesgue measure, so [L1]'s boundedness and finite-measure hypotheses fail. For the motion through (0,e1), x(t)=te1: if a fixed bounded phase set has ∥x∥2≤R on it, this trajectory can be in the set only for 0≤t≤R. It eventually exits and makes only finitely many integer-time visits; this does not claim that no finite return occurs. The finite-volume hypothesis is a genuine restriction.

Remarks

ExampleConstruction: AI-adaptedVerification: AI-adaptedOpen item page →

Unit restoration and scale checks for the two-rotor KAM statement: the scaling identities at I=E0=1, the physical angular frequencies, the exclusion of rational frequency labels, and the equal-rescaling invariance

Physical scope

A declared-data check of the units and labels in the page's two-rotor KAM theorem: with moment of inertia 1 kg m^2 and energy scale 1 J the scaling constants are J0=1 kg m^2 s^-1, t0=1 s and omega0=1 s^-1; the physical and dimensionless Hamilton equations correspond exactly; a surviving torus with dimensionless label nu carries physical angular frequencies omega0 nu_i; every rational label of the frequency cube is excluded from the Diophantine set; and equal rescaling of inertia and energy leaves the dimensionless model and the threshold unchanged. The constant c in the smallness condition is not numerically evaluated and no measured value is asserted.

Example

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)), as in the cited two-rotor KAM theorem.

Declared model data. Take the two-rotor model of The two-rotor model: a positive-measure family of invariant Lagrangian tori with Diophantine frequencies for sufficiently small coupling, with the explicit units and the stated scope with I:=1 kg m2,E0:=1 J,ε dimensionless, so that J0=IE0=1 kg m2 s−1, t0=J0/E0=1 s and ω0=E0/I=1 s−1.

Assertions to check. (i) The dimensions [Hphys]=J, [Πi]=kg m2 s−1, [I]=kg m2, [ε]=1 and [ω0]=s−1 (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework). (ii) With Ii=Πi/J0, s=t/t0 the physical Hamilton equations θ˙i=Πi/I, Π˙i=E0εsin⁡θicos⁡θ3−i are equivalent to θi′=Ii, Ii′=εsin⁡ϑicos⁡ϑ3−i. (iii) A surviving torus with dimensionless label ν carries the physical angular frequencies ω0νi, here νi s−1. (iv) Every rational label in the cube C∗ is excluded from Ωα,τ for every α>0 and τ>1. (v) Replacing (I,E0) by (λI,λE0) with λ>0 leaves ω0, the dimensionless model and the smallness condition unchanged, while rescaling J0 by λ.

Facts & Assumptions

Given: the declared data I=1 kg m2, E0=1 J and the scaling constants J0,t0,ω0; the frequency cube C∗ and the Diophantine set Ωα,τ of The two-rotor model: a positive-measure family of invariant Lagrangian tori with Diophantine frequencies for sufficiently small coupling, with the explicit units and the stated scope.

[L1]

The two-rotor KAM statement and its scaling (The two-rotor model: a positive-measure family of invariant Lagrangian tori with Diophantine frequencies for sufficiently small coupling, with the explicit units and the stated scope): the physical Hamiltonian is Π12/(2I)+Π22/(2I)+E0εcos⁡θ1cos⁡θ2 with J02=IE0 and t0=J0/E0; under Ii=Πi/J0, s=t/t0 it becomes 12(I12+I22)+εcos⁡ϑ1cos⁡ϑ2, and each surviving torus carries φ(t)=φ0+ω0νt in its torus parameter, with the physical angle given by the embedding.

[L2]

Diophantine labels (Analytic KAM data: angular tori, action-angle neighbourhoods, Fourier cutoffs and Diophantine frequency labels): a frequency ν is (α,τ)-Diophantine exactly when ∣k⋅ν∣≥α∣k∣1−τ for every nonzero integer vector k.

Verification

technique · direct
1.1givenL1algebra

Dimensions: the physical Hamiltonian has [Πi2/I]=kg2m4s−2/(kg m2)=kg m2 s−2=J and [E0ε]=1⋅J with ε dimensionless; ω0=E0/I has [ω0]=kg m2 s−2/(kg m2)=s−1 (Dimensions and SI unit conventions of the non-relativistic classical mechanics framework). With the declared values, J0=1 kg m2 s−1, t0=1 s, ω0=1 s−1.

1.2givenL1algebra

Scaling of the equations, first equation: with Πi=J0Ii and t=t0s the chain rule gives θ˙i=t0−1θi′ and Π˙i=t0−1J0Ii′; substituting into θ˙i=Πi/I and using J0t0=I yields t0−1θi′=J0Ii/I=Ii/t0, that is the dimensionless equation θi′=Ii.

1.3givenL2algebra

Rational labels are non-Diophantine: let ν=(p1/q,p2/q)∈C∗ with integers p1,p2 and q≥1. If p1=p2=0, choose k=(1,0); otherwise set k:=(p2,−p1)≠0; in either case k⋅ν=0<α∣k∣1−τ, so ν∉Ωα,τ by [L2]. Hence the Diophantine set of the theorem contains no rational point of the cube; the positive measure it retains is carried by irrational labels.

2.1step 1.2givenL1algebra

Scaling of the equations, second equation: substituting into Π˙i=E0εsin⁡θicos⁡θ3−i gives t0−1J0Ii′=E0εsin⁡θicos⁡θ3−i; since J0/t0=J0E0/J0=E0, this is equivalent to Ii′=εsin⁡θicos⁡θ3−i, the dimensionless Hamilton equation of the potential εcos⁡ϑ1cos⁡ϑ2. With the declared values J0/t0=1 J and E0=1 J the substitution is an identity.

2.2step 1.1givenL1algebra

Physical frequencies: by [L1] a surviving torus with label ν carries φ(t)=φ0+ω0νt in its torus parameter; the physical angle is its embedded angular component; with ω0=1 s−1 the components ω0νi=νi s−1 are the physical angular frequencies, of dimension s−1 as in step 1.1.

3.1step 1.1step 1.2givenL1algebra∎

Equal rescaling: replace (I,E0) by (λI,λE0), λ>0; then ω0=E0/I is unchanged, J0=IE0 is multiplied by λ and t0=J0/E0 is unchanged, so the dimensionless Hamiltonian, the frequency cube, the Diophantine set and the smallness condition are all unchanged, while the physical momenta Πi=J0Ii are rescaled by λ. The constant c of the theorem is not numerically evaluated in this check; only the exact scaling behaviour is asserted.

Remarks

Sources