Alphabeta Math
TheoremStatement: AI-adaptedProof: AI-adaptedPipeline-generated
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.

Instability of spin about the intermediate principal axis of the torque-free top

Statement

Let I1<I2<I3 be positive reals, let I:R3→R3 be the linear map with Iej=Ijej on the standard orthonormal basis, and consider the torque-free Euler system l′=F(l), F(l)=l×(I−1l), on U=R3 (Euler's equations Iω′+ω×Iω=RTτ under an adopted torque balance, and the torque-free invariants and principal-axis equations, Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system). Let L≠0. Then l0:=Le2 is an equilibrium of F and is not Lyapunov stable: with ε0:=∣L∣2,a:=1/I1−1/I21/I2−1/I3>0,c0:=a(1I2−1I3)32∣L∣>0, for every δ>0 there is a solution l with ∥l(0)−l0∥2<δ and some t≥0 with ∥l(t)−l0∥2≥ε0; Moreover, with δ0:=∣L∣/(21+a2), the escape can be arranged at a time t≤c0−1log⁡(δ0/η) where η=l1(0)>0 may be chosen as small as desired.

Facts & Assumptions

Given: Positive reals I1<I2<I3, the diagonal map I, the field F(l)=l×(I−1l) on U=R3, and L≠0.

[L1]

The system l′=F(l) is read as First-order systems, initial value problems, and solutions on intervals. For ∥x∥2,∥y∥2≤r, bilinearity gives F(x)−F(y)=(x−y)×(I−1x)+y×(I−1(x−y)); since I−1 is diagonal with entries at most 1/I1 in absolute value and ∥u×v∥2≤∥u∥2∥v∥2, this yields ∥F(x)−F(y)∥2≤(2r/I1)∥x−y∥2. This bound makes F locally Lipschitz in the state variable and also continuous at each state point by the ε-δ definition (Local Lipschitz continuity in the state variable, locally uniform in time and parameters, The cross product in R3, The squared cross-product norm is the Gram determinant of two vectors). Picard-Lindelof gives local existence and uniqueness, hence a unique maximal solution through every datum (Picard-Lindelöf local existence and uniqueness for first-order systems, Extensions, maximal solutions, maximal intervals, and global solutions, Every Picard–Lindelöf initial value problem has one maximal solution on an open interval); a solution trapped in a compact subset of U=R3 near a finite endpoint extends past it (First-order systems, initial value problems, and solutions on intervals, A non-extendable solution leaves every compact subset of the domain near its finite endpoint; an autonomous solution trapped in a compact subset of its state domain cannot break down there).

[L2]

Along every solution the invariants C(l)=∥l∥22 and E(l)=12⟨l,I−1l⟩ are constant; in the eigenbasis Iej=Ijej one has C=∑jlj2, E=∑jlj2/(2Ij), and the component equations l1′=(1/I3−1/I2)l2l3, l2′=(1/I1−1/I3)l1l3, l3′=(1/I2−1/I1)l1l2 The same algebra as in Euler's equations Iω′+ω×Iω=RTτ under an adopted torque balance, and the torque-free invariants and principal-axis equations applies directly to the present diagonal map: I−1 is symmetric, so C′=2⟨l×(I−1l),l⟩=0 and E′=⟨l×(I−1l),I−1l⟩=0 by cross-product orthogonality. Both scalar functions are continuous and differentiable on the solution interval; A function continuous on an interval I whose derivative vanishes at every interior point of I is constant on I; consequently two such functions with the same derivative differ by a constant therefore makes them constant there. These computations require only the given positive diagonal entries, and no realization as a mass-distribution inertia matrix.

[L3]

The Euler field is even: F(−l)=(−l)×(I−1(−l))=F(l), by bilinearity and linearity (The cross product in R3). Put σ:=L/∣L∣∈{−1,1}, so σ2=1 and L=σ∣L∣.

[L8]

For r>0, log⁡r is the inverse value of the strictly increasing exponential: exp⁡(log⁡r)=r (The natural logarithm as the inverse of the exponential function, The exponential is a continuous bijection from R onto (0,∞), The exponential function is strictly increasing).

[L7]

The Euclidean norm satisfies the triangle inequality, so ∣∥x∥2−∥y∥2∣≤∥x−y∥2; hence each sphere {z:∥z∥2=r} is closed (Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation).

Proof

technique · direct
1.1givenL1L3

F(Le2)=Le2×(L/I2)e2=0, so l0 is an equilibrium. Put σ=L/∣L∣ as in [L3]. Choose the second and third coordinates of the initial datum with signs σ and −σ, respectively.

2.1step 1.1givenL2L3L4algebra

Let 0<η0<δ0:=∣L∣/(21+a2) and let 0<η<η0. Put s(η):=∣L∣2−(1+a2)η2>0 and lη:=(η, σs(η), −σaη), so that C(lη)=L2 and E(lη)=L2/(2I2): indeed C(lη)=η2+∣L∣2−(1+a2)η2+a2η2=L2 and E(lη)=12(η2/I1+(∣L∣2−(1+a2)η2)/I2+a2η2/I3)=12L2/I2+12η2(1/I1−1/I2−a2(1/I2−1/I3))=L2/(2I2) by the definition of a. Also ∥lη−l0∥2→0 as η↓0, since s(η)→∣L∣.

3.1step 2.1givenL1L2L4algebraL6L7

Let l be the unique maximal solution with l(0)=lη; it is defined for all t≥0, because its invariant sphere C=L2 is closed by [L7] and bounded, hence compact by Heine-Borel in Rn: with the Euclidean metric a subset of Rn is compact if and only if it is closed and bounded, and the proof by bisection uses no choice principle; the same holds on the real line, and [L1] rules out breakdown inside a compact set. Let τ:=inf⁡{t>0:l1(t)≥δ0 or l1(t)≤0}, so 0<τ≤+∞, and on (0,τ) one has 0<l1<δ0. On that interval l1(t)2+l3(t)2≤l12+a2l12<δ02(1+a2)=L2/4 by [L2] (which gives l3=±al1 on the invariant level set E=L2/(2I2), since E−C/(2I2)=0 and (1/I2−1/I3)>0), so l2(t)2=L2−l12−l32>3L2/4 and l2 has the constant sign σ on [0,τ), its initial sign. Likewise l3=±al1 cannot change sign on [0,τ) without vanishing, which would force l1=l3=0 contrary to l1>0; since l3(0)=−σaη, one has l3(t)=−σal1(t) there.

4.1step 3.1L2algebra

On (0,τ), the first component equation of [L2] and step 3.1 give l1′=(1/I3−1/I2)l2l3=a(1/I2−1/I3)(σl2)l1 ≥c0l1 with c0=a(1/I2−1/I3)32∣L∣>0, because σl2≥32∣L∣.

5.1step 3.14.1givenL5L8algebra

By step 4.1 and the product and chain rules [L5], (e−c0tl1(t))′=e−c0t(l1′(t)−c0l1(t))≥0 on (0,τ); since e−c0tl1(t) is continuous on [0,τ), it is nondecreasing there by [L5], so l1(t)≥ηec0t for t∈[0,τ). If τ=+∞, this would give l1(t)≥ηec0t>δ0 for t>c0−1log⁡(δ0/η), contradicting l1<δ0 on (0,τ); hence τ<+∞.

6.1step 3.15.1givenL4algebra

At the time τ one has l1(τ)=δ0: indeed l1(τ)≠0 by step 5.1 and continuity, so by the definition of τ and continuity it is the alternative l1(τ)≥δ0 that occurs, and on [0,τ] the bound l1<δ0 holds on (0,τ) with l1 continuous, so l1(τ)=δ0; moreover τ≤c0−1log⁡(δ0/η). Therefore ∥l(τ)−l0∥22≥l1(τ)2+l3(τ)2=(1+a2)δ02=L2/4, so ∥l(τ)−l0∥2≥∣L∣/2=ε0.

7.1step 2.16.1givenalgebra∎

It remains to make the initial point as close to l0 as prescribed: for any δ>0 choose η>0 with η<min⁡{δ0,δ} and ∥lη−l0∥2<δ, possible because ∥lη−l0∥2→0 by step 2.1. The corresponding solution then has ∥l(0)−l0∥2<δ and, by step 6.1, escapes the ball of radius ε0 at time τ≤c0−1log⁡(δ0/η). This is instability of l0 in the sense of Equilibria, Lyapunov stability and asymptotic stability of an autonomous first-order system, with the explicit time bound of the statement.

Depends on

Dependency tree · 2 levels

Sources