Alphabeta Math
Physics theoremStatement: AI-adaptedProof: AI-adaptedprecheck 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.

Finite-time error control for the small-oscillation approximation: the quadratic remainder of the mechanical field and the variation-of-constants bound

Physical scope

Finite-time error control for the small-oscillation model of a mechanical system in a chart at a mechanical equilibrium: for a C^3 kinetic metric field and a C^3 potential with vanishing gradient at q_0, the mechanical field is C^2 near the equilibrium, its Taylor remainder at the linearization is quadratic with constants obtained from compactness, and the finite-time variation-of-constants estimate bounds the difference between the nonlinear motion and the linearized (normal-mode) motion by C delta^2 t exp((||A||+2L)t) on every horizon T with delta exp(LT)<r_0. The statement is a conditional finite-time estimate for the declared chart model under explicit regularity and smallness hypotheses; it asserts no all-time validity of the linear approximation, no stability conclusion, no large-data statement and no empirical claim, and the constants depend on the model data only through its behaviour on the stated ball.

Statement

Fix positive coordinate units si, a time unit t∗ and an energy unit e∗ as in The coupled quadratic oscillator model: mass metric, stiffness, generalized force, modal vocabulary and the small-oscillation model at a mechanical equilibrium. Throughout this theorem qi=Qi/si, t=tphys/t∗ and vi=(t∗/si)Q˙i are dimensionless numerical chart and phase coordinates. M and U are the corresponding normalized metric and potential specified in that definition. In particular (δq,δv) is a vector in one Euclidean numerical space; r0,δ,C,L,A,T below refer to that space and numerical time. All constants are calculated after this normalization. If the numerical error bound is b(t), its physical position and velocity components are bounded respectively by sib(t) and (si/t∗)b(t) at physical time t∗t. Fixed positive linear coordinate/time changes preserve the stated regularity and mechanical equation, so no original chart model is excluded.

Let n≥1, let Ω⊆Rn be open, let M:Ω→Mn(R) be a field of class C3 with M(q) symmetric positive definite for every q∈Ω, and let U:Ω→R be of class C3 with ∇U(q0)=0 at a point q0∈Ω. Let G be the mechanical vector field of The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart on Ω×Rn, let r0>0 satisfy B‾(q0,r0)⊆Ω and define F:{ζ=(δq,δv):q0+δq∈Ω}→R2n by F(δq,δv):=G(q0+δq,δv). Put A:=DF(0)=DG(q0,0). Then:

  1. F is of class C2, F(0)=0, F is locally Lipschitz in the state variable, and DG(q0,0)(δq,δv)=(δv, −M(q0)−1D2U(q0) δq); moreover there are constants C,L≥0 with ∥F(ζ)−Aζ∥2≤C∥ζ∥22,∥F(ζ)∥2≤L∥ζ∥2(∥ζ∥2≤r0); one may take L=sup⁡∥η∥2≤r0∥DF(η)∥ (The operator norm as the least bound and as the unit-sphere or unit-ball supremum) and C=Λ/2 for any Lipschitz constant Λ of DF on the closed ball B‾(0,r0)⊆R2n (Open ball, closed ball and sphere in a metric space).
  2. For every datum z0=(q0+δq0,v0) with δ:=∥(δq0,v0)∥2>0 and every T>0 with δeLT<r0, the maximal solution z=(q,v) of z′=G(z) with z(0)=z0 is defined on [0,T] (Extensions, maximal solutions, maximal intervals, and global solutions) and satisfies, for every t∈[0,T], ∥z(t)−y(t)∥2≤Cδ2 t e(∥A∥+2L)t, where y(t)=(q0+ξ(t),η(t)) and y~=(ξ,η):R→R2n solves y~′=Ay~, y~(0)=(δq0,v0). Thus y is the linear approximation translated back to the original chart.
  3. Writing y~=(ξ,η) for the displacement linear solution of clause 2, the first component ξ is exactly the unique motion of the small-oscillation model at q0 of clause (d) of The coupled quadratic oscillator model: mass metric, stiffness, generalized force, modal vocabulary and the small-oscillation model at a mechanical equilibrium with datum (δq0,v0): ξ is C2 on R with M(q0)ξ′′+D2U(q0)ξ=0, ξ(0)=δq0 and ξ′(0)=v0.

No all-time statement, no statement for data violating the smallness condition δeLT<r0, and no value of any measured quantity is asserted.

Facts & Assumptions

Given: The C3 metric field M, the C3 potential U with ∇U(q0)=0, the mechanical field G, the radius r0, the translated field F, the matrix A=DF(0) and a datum z0=(q0+δq0,v0) with δ>0.

[F1]

Regularity of the mechanical field (The mechanical first-order field, Lagrangian and energy of a configuration-dependent kinetic metric in a chart): G(q,v)=(v,M(q)−1B(q,v)) with B(q,v)=12∇q(vTM(q)v)−(DqM(q)[v])v−∇U(q); the maps M, ∂kM, ∇q(vTM(q)v), (DqM(q)[v])v, ∇U and q↦M(q)−1 are respectively C3, C2, C2, C2, C2 and C3 under the stated hypotheses, because M−1=det⁡(M)−1adj(M) has entries polynomial in the entries of M divided by the nonvanishing continuous function det⁡M (If det⁡(A) is a unit, then A−1=det⁡(A)−1adj⁡(A), If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative, Ck Euclidean maps and diffeomorphisms); hence G and F are C2, so F is locally Lipschitz in the state variable (Local Lipschitz continuity in the state variable, locally uniform in time and parameters).

[F2]

Linearization and equality of derivatives (Linearisation of the mechanical system at an equilibrium: M(q_0) xi'' + D²U(q_0) xi = 0, and the small-oscillation quadratic model clauses 1-2): G is totally differentiable at (q0,0) and DG(q0,0)(δq,δv)=(δv,−M(q0)−1D2U(q0)δq); D2U(q0) is symmetric (The Hessian of a C2 scalar field is symmetric, The Hessian matrix and critical points of a scalar field). Since F is the translation of G to the equilibrium, A=DF(0)=DG(q0,0) and F(0)=G(q0,0)=0 because q0 is an equilibrium.

[F3]

Compactness and mean value (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, The mean value inequality: if f:[a,b]→Rm is continuous and differentiable on (a,b) with ∥f′∥2≤M, then ∥f(b)−f(a)∥2≤M(b−a), Open ball, closed ball and sphere in a metric space): the closed ball B‾(0,r0) is a compact convex subset of the open domain of F; a continuous real function on it attains its maximum, and a continuous map with ∥DF∥≤L on the ball satisfies ∥F(η)−F(η′)∥2≤L∥η−η′∥2 there.

[F4]

Fundamental theorem along rays (If f:[a,b]→Rm is differentiable with integrable f′ then ∫abf′=f(b)−f(a); and a bounded derivative makes f Lipschitz, For a≤b and f:[a,b]→Rm integrable when a<b, ∥∫abf∥2≤∫ab∥f∥2; for a<b, ∥f∥2 is integrable, The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral): for ∥ζ∥2≤r0 the curve t↦F(tζ) on [0,1] is C1 with derivative DF(tζ)ζ, and F(ζ)−F(0)=∫01DF(tζ)ζ dt,∥∫01ϕ(t) dt∥2≤∫01∥ϕ(t)∥2 dt.

[F5]

Finite-time linearization bound (Finite-time linearization error bound: quadratic remainder, exponential growth bound and the variation-of-constants estimate with state dimension 2n): for the open domain of F, the constants r0,C,L and the matrix A of the statement, the maximal solution ζ of ζ′=F(ζ), ζ(0)=ζ0:=(δq0,v0) is defined on [0,T] whenever δeLT<r0, it satisfies ∥ζ(t)∥2≤δeLt there, and ∥ζ(t)−y~(t)∥2≤Cδ2∫0te∥A∥(t−s)e2Ls ds≤Cδ2te(∥A∥+2L)t for t∈[0,T], where y solves y~′=Ay~, y~(0)=ζ0; every initial value problem has a unique maximal solution (Every Picard–Lindelöf initial value problem has one maximal solution on an open interval, Extensions, maximal solutions, maximal intervals, and global solutions).

[F6]

Model dictionary (The coupled quadratic oscillator model: mass metric, stiffness, generalized force, modal vocabulary and the small-oscillation model at a mechanical equilibrium clauses (b), (d)): the small-oscillation model at q0 is the coupled quadratic model with mass metric matrix M(q0) and stiffness D2U(q0), whose phase system is y′=Ay with the displayed block matrix A=(0In−M(q0)−1D2U(q0)0); a C2 curve ξ is a motion of it if and only if M(q0)ξ′′+D2U(q0)ξ=0, with the datum read off (ξ(0),ξ′(0)).

Proof

Proof technique: direct; prove the quadratic-remainder bounds for the translated mechanical field by compactness and the fundamental theorem along rays, then apply the finite-time linearization theorem and identify the linear solution with the small-oscillation motion.

1.1F1F2

Clauses 1 and the derivative identity: by [F1] the translated field F is C2 on its open domain and locally Lipschitz; F(0)=0 and A=DF(0)=DG(q0,0) with the displayed block form are [F2].

1.2F1F3F4algebra

Quantitative remainder and growth bounds: since F is C2, DF is C1, hence continuous on the compact ball B‾(0,r0); thus L:=sup⁡∥η∥2≤r0∥DF(η)∥ is finite, and DF is Lipschitz on the convex ball with some finite constant Λ (each entry of DF has continuous gradient and is Lipschitz on the ball by [F3] with the supremum of that gradient; a finite sum of such bounds bounds ∥DF(η)−DF(η′)∥). For ∥ζ∥2≤r0, [F4] applied to the C1 curve t↦F(tζ)−tAζ, whose derivative is DF(tζ)ζ−Aζ=(DF(tζ)−DF(0))ζ, gives F(ζ)−Aζ=∫01(DF(tζ)−DF(0))ζ dt, hence ∥F(ζ)−Aζ∥2≤∫01Λt∥ζ∥22dt=Λ2∥ζ∥22; and F(ζ)=F(ζ)−F(0) with [F3] gives ∥F(ζ)∥2≤L∥ζ∥2. This proves the bounds of clause 1 with C=Λ/2.

2.1F5step 1.2given

The finite-time bound: apply [F5] to F, the constants r0,C,L, the matrix A and the displaced datum ζ0=(δq0,v0) of norm δ. The maximal solution ζ of ζ′=F(ζ) is defined on [0,T] and obeys the error bound against y~. Translating back gives z(t)=(q0+ζq(t),ζv(t)), the physical phase solution through z0; y(t)=(q0+ξ(t),η(t)) has exactly the same translation, so ∥z(t)−y(t)∥2=∥ζ(t)−y~(t)∥2; this is clause 2.

3.1F2F6step 2.1algebra

Identification of the linear solution: by [F2] the matrix A is the block matrix (δq,δv)↦(δv,−M(q0)−1D2U(q0)δq); the displacement linear solution y~=(ξ,η) of y~′=Ay~, y~(0)=(δq0,v0) therefore satisfies η=ξ′, hence ξ′′=−M(q0)−1D2U(q0)ξ, equivalently M(q0)ξ′′+D2U(q0)ξ=0, with ξ(0)=δq0, ξ′(0)=v0; by [F6] this is exactly the motion of the small-oscillation model with that datum, and it is unique by the model's phase-system uniqueness. This is clause 3.

4.1step 1.1step 1.2step 2.1step 3.1∎

Clauses 1-3 are proved. The estimate is a finite-time statement under the explicit smallness hypothesis δeLT<r0; the constants C,L exist by compactness and are not asserted to be optimal, and no long-time, stability or empirical conclusion is drawn.

Remarks

  • Where the exponential factor comes from. The linearized flow itself may grow exponentially at rate ∥A∥ (an unstable equilibrium included), and the nonlinear remainder is controlled by the Lipschitz growth e2Lt of the solution; the product is the displayed rate ∥A∥+2L. On a fixed horizon this rate is a constant, and the bound is of order δ2 in the initial displacement.

  • Constant versus asymptotic statements. The theorem fixes a horizon T and shrinks the datum; it does not assert that the nonlinear motion stays close to the linear one for all time, nor that a small initial displacement produces a small displacement later. Such all-time statements are separate problems (for the stable case they concern nonlinear stability, for the unstable case nonlinear instability) and are not obtained from this estimate.

Depends on

Dependency tree · 2 levels

Sources