Alphabeta Math
PropositionStatement: AI-adaptedProof: AI-adaptedPipeline-generatedprecheck pass
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.

The derivative of the Hamiltonian flow satisfies YTJY=J and det⁡Y=1

Statement

Let n≥1, let U be a phase domain with standard two-form ωU and matrix J, let H∈C2(U), and let Φ be the local flow of the Hamiltonian field XH=J∇H (The Hamiltonian field of a C2 function on a phase domain is J∇H, and the Poisson bracket is ∑i(∂qif ∂pig−∂pif ∂qig), The local flow of an autonomous first-order system). Let z0∈U, let I⊆R be the open interval of times t with (t,z0) in the flow domain, and for t∈I put

Y(t):=DzΦt(z0)∈M2n(R).

Then t↦Y(t) is of class C1 on I and, for every t∈I:

  1. (variational equation) Y′(t)=J Hess H(Φt(z0)) Y(t) and Y(0)=I2n;
  2. (canonical matrix identity) Y(t)TJ Y(t)=J;
  3. (unit determinant) det⁡Y(t)=1.

The same conclusions hold for a time-dependent Hamiltonian H:T×U→R of class C2, where T is an open real interval. For t0∈T and z0∈U, write E(t,t0,z) for the maximal solution of z′=J∇zH(t,z) with z(t0)=z; its domain consists of the triples where this solution exists. On the maximal open time interval It0,z0 set Y(t)=DzE(t,t0,z0). Then Y′=J Hess⁡zH(t,E(t,t0,z0))Y,Y(t0)=I2n,YTJY=J,det⁡Y=1. The evolution E is C1 on its open domain, and each fixed-endpoint map has inverse E(t0,t,⋅) on its image. Autonomous H is recovered by H(t,z)=H(z) and E(t,t0,z)=Φ(t−t0,z).

Facts & Assumptions

Given: The phase domain U, the C2 Hamiltonian H, the local flow Φ, the point z0∈U and the derivative Y(t)=DzΦt(z0); for the additional case the open interval T, jointly C2 Hamiltonian H, and initial time t0∈T are given.

[F2]

The Hessian Hess H(z) is the matrix (∂i∂jH(z)), and for H∈C2 it is symmetric: Hess H(z)T=Hess H(z) (The Hessian matrix and critical points of a scalar field, The Hessian of a C2 scalar field is symmetric).

[F3]

Block identities JT=−J, J2=−I2n, J−1=JT and ωz(v,w)=vTJw (The standard two-form of a phase domain is closed and nondegenerate, with contraction ιvω=∑iai dpi−∑ibi dqi).

[F4]

Determinant calculus: the determinant is a polynomial in the matrix entries (For every fixed finite size at least one, the determinant of a real square matrix is a polynomial in its matrix entries), and a finite sum of finite products of continuous real functions is continuous (Sums, products, absolute values, finite maxima and minima, and quotients of continuous real-valued maps on a topological space are continuous where defined, with the metric and topological notions of continuity identified by Metric continuity characterisations, with countable choice for the sequential converse); the entry functions t↦Yij(t) are continuous as components of the continuous map Y (Continuity of a map between metric spaces, at a point and globally, in the ε-δ form), so t↦det⁡Y(t) is continuous on I. The Leibniz formula gives det⁡I2n=1 (For n≥1, the determinant over a commutative ring by the Leibniz formula, and ∣det⁡A∣ for a real matrix); moreover det⁡(AB)=det⁡(A)det⁡(B) and det⁡(AT)=det⁡(A) for square matrices (For same-sized finite square matrices over a commutative ring, det⁡(AB)=det⁡(A)det⁡(B), For every square matrix over a commutative ring, det⁡(AT)=det⁡(A)), and det⁡J≠0 because J is invertible with J−1=JT=−J and over the field R a square matrix is invertible exactly when its determinant is nonzero (A positive-sized square matrix over a commutative ring is invertible if and only if its determinant is a unit, Invertible Euclidean linear maps).

[F5]

Sign rigidity on an interval (Intermediate value theorem, by bisection with a canonical left-half rule: a continuous function on [a,b] takes every value between f(a) and f(b), Intervals of R: the nine order-convex forms, nondegeneracy, and length): a continuous g:I→R on an interval with g(t)2=1 for all t∈I and g(0)=1 equals 1; indeed, if g(t)=−1 for some t>0, the intermediate value theorem on [0,t] gives c with g(c)=0, contradicting g(c)2=1, and the case t<0 is identical on [t,0].

Proof

technique · direct; the variational equation gives the derivative, the antisymmetry of $J$ cancels the quadratic form, and continuity upgrades the squared determinant
1.1F1algebra

Clause 1: by [F1] the initial-state derivative Y exists on I. On a small cylinder the local variational theorem makes it C1 with the stated initial condition. For any compact reference orbit segment from 0 to an admissible t, finitely many such local cylinders cover the segment. Flow composition expresses the initial-state derivative as a product of the local derivatives; on each successive patch the chain rule differentiates this product in the last time variable and gives the same variational equation with the preceding product as its initial value. This extends the C1 derivative and its equation to the entire segment, for negative times using the reversed field. Since t is arbitrary, Y is C1 on I and solves Y′=DF(Φt(z0))Y with Y(0)=I2n; since XH=J∇H and J is constant, the chain rule gives DF(z)=J Hess H(z), so Y′(t)=J Hess H(Φt(z0))Y(t).

1.2F1givenconstruct

For the time-dependent case use the autonomous C1 field F~(s,z)=(1,J∇zH(s,z)) on the open set T×U. Its maximal flow Φ~ has open domain and is C1 by The flow domain of a C¹ autonomous field is open and the flow is C¹ on it. Its first coordinate is s+r at elapsed time r, so define E(t,t0,z) as the second coordinate of Φ~(t−t0,(t0,z)). The preimage of the open flow domain under this change of variables is open, and E is C1 by composition. The composition and inverse laws for the extended autonomous flow give E(t0,t,E(t,t0,z))=z and the reverse identity, with both domains included; thus fixed-endpoint maps are C1 diffeomorphisms onto open images.

2.1F1F2F3step 1.1algebra

Clause 2: put A(t):=J Hess H(Φt(z0)). By [F2] and [F3], AT=(Hess H)TJT=(Hess H)(−J)=−Hess H J and JA=J2Hess H=−Hess H, so ATJ+JA=−Hess H JJ−Hess H=(Hess H)(−J2)−Hess H=0; hence ddt(YTJY)=Y′TJY+YTJY′=YT(ATJ+JA)Y=0, and Y(t)TJY(t)=Y(0)TJY(0)=J for all t∈I.

3.1step 2.1F4F5algebra

Clause 3: taking determinants in clause 2 and using [F4] gives det⁡(YTJY)=det⁡(YT)det⁡(J)det⁡(Y)=(det⁡Y)2det⁡J, so (det⁡Y(t))2det⁡J=det⁡J; since det⁡J≠0, (det⁡Y(t))2=1 for all t∈I; by [F4] the function t↦det⁡Y(t) is continuous with det⁡Y(0)=det⁡I2n=1, so the sign rigidity [F5] gives det⁡Y≡1.

4.1F1F2F3F4F5step 2.1step 3.1step 1.2∎

Apply the autonomous smooth-dependence/variational theorem to F~ on small balls along an extended solution, and use the flow composition law to concatenate the initial-state derivatives over finitely many such local patches on each compact subinterval of its lifetime. This gives the variational equation for the derivative with respect to the original z-initial coordinates on the full maximal interval. The first coordinate t does not depend on those coordinates, so its derivative block is zero; the remaining block is Y, and its equation is Y′=JHess⁡zH(t,E(t,t0,z0))Y with Y(t0)=I. The z-Hessian is symmetric at each time. The cancellation in step 2.1 therefore applies with this time-dependent symmetric Hessian, giving ddt(YTJY)=0 and YTJY=J. The determinant argument of step 3.1 applies on the interval containing t0 with the initial sign det⁡Y(t0)=1, yielding det⁡Y=1. All conclusions, including the nonautonomous case, follow.

Remarks

Depends on

Dependency tree · 2 levels

Sources