Alphabeta Math
Pipeline-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.

Cm Math Geometry — Examples

1 · Prerequisites

2 · Summary

Finite verifications of the definitions and lemmas of cm-math-geometry. Worked change of adapted chart with a rotation, a boost and a spatial shift works a change of adapted chart with a rotation, a boost and a spatial shift through explicit coordinates, checks R∈SO(3) and the transformation law x′=Rx+vt+a, and verifies the velocity and acceleration transformation on a stated curve. The rotating-frame acceleration identity for a uniform rotation about an axis computes Ω=RTR′ for a uniform rotation and checks the acceleration identity of Angular velocity of a moving frame and the rotating-frame acceleration identity on two explicit curves, isolating the Coriolis-type and centrifugal-type terms. The mass metric of a particle constrained to the unit sphere in spherical coordinates establishes the smooth two-dimensional sphere and the spherical chart away from the poles and longitude seam, then evaluates the mass metric of the inclusion S2↪R3 in spherical coordinates, obtaining m(θ˙2+sin⁡2θ φ˙2) on the chart domain and noting the chart degeneration at the poles. A collinear labeling admits a nontrivial rotation refutes the claim that every nonempty labeling has trivial SO(3)-stabilizer by exhibiting the half-turn about the axis of {0,e3}, complementing A labeling containing two independent vectors has trivial SO(3)-stabilizer. A closed field on a star-shaped domain with its radial primitive computed differentiates a closed field on the plane, exhibits its potential and evaluates the radial primitive of the published Poincare lemma, with the punctured-plane obstruction recorded as a separate published counterexample and not used as a premise. The stroboscopic map of the plane oscillator computes the time-P map of the undamped plane oscillator, shows DP=R2 for every P, and identifies the identity case P=2π/ω and the half-turn case P=π/ω. All items are mathematics and nothing outside this pair depends on them.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

None yet.

5 · Examples, experiments and counterexamples

ExampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

Worked change of adapted chart with a rotation, a boost and a spatial shift

Example

Take V=R4 with coordinates (t,x1,x2,x3), τ(t,x1,x2,x3)=t, so that K={0}×R3 with the standard inner product g, and take A=V as an affine space over itself. Let the first adapted chart be the standard one, o=0, u=e0, (e1,e2,e3) the standard basis of K. Let the second adapted chart be

o′=e2,u′=e0+e1,e1′=e2,e2′=−e1,e3′=e3,

which is adapted: τ(u′)=1, and (e1′,e2′,e3′) is a positively oriented g-orthonormal basis of K.

For a point with first coordinates (t,x1,x2,x3), equating the two expansions gives t′=t, x2=x1′+1 and x1=t′−x2′, that is

x1′=x2−1,x2′=t−x1,x3′=x3,

which is the transformation x′=Rx+vt+a of Transformation between adapted charts and its kinematic consequences with b=0,

R=(010−100001),v=(010),a=(−100).

The matrix R satisfies RTR=I3 and det⁡R=1, so R∈SO(3).

Facts & Assumptions

Given: The two adapted charts displayed above, a point p with first coordinates (t,x1,x2,x3), and the map γ:R→A, γ(s)=(s,s,2s,0) in the first coordinates.

[L1]

Adapted coordinates are unique: the coefficients in p=o+tu+∑ixiei and in the primed expansion are uniquely determined (Existence and uniqueness of adapted coordinates, and the affine structure of simultaneous slices).

[L2]

Two adapted charts are related by t′=t+b and x′=Rx+vt+a with R∈O(3), and R∈SO(3) when both bases are positively oriented (Transformation between adapted charts and its kinematic consequences, The orthogonal group O(3) and the rotation group SO(3)).

Verification

technique · direct; substitute the chart data and check the displays
1.1givenalgebra

The matrix R has rows (0,1,0), (−1,0,0), (0,0,1), so RTR has entries: first diagonal 02+(−1)2+02=1, second diagonal 12+02+02=1, third diagonal 1, and off-diagonal entries 0⋅1+(−1)⋅0+0⋅0=0, 0⋅0+(−1)⋅0+0⋅1=0, 1⋅0+0⋅0+0⋅1=0; hence RTR=I3, and det⁡R=0⋅(0⋅1−0⋅0)−1⋅((−1)⋅1−0⋅0)+0=1, so R∈SO(3).

2.1step 1.1L1L2algebra

For the point p with first coordinates (t,x1,x2,x3) the direct substitution p=o′+t′u′+∑ixi′ei′ reads (t,x1,x2,x3)=(t′, t′−x2′, 1+x1′, x3′), whose unique solution is t′=t, x1′=x2−1, x2′=t−x1, x3′=x3; comparison with x′=Rx+vt+a gives Rx+vt+a=(x2,−x1,x3)+(0,t,0)+(−1,0,0)=(x2−1, t−x1, x3), so the transformation law of the lemma holds in this chart pair with the displayed b,v,a,R.

3.1step 2.1L2algebra∎

For the curve γ(s)=(s,s,2s,0), with first coordinate representative x(s)=(s,2s,0), the second representative is y(s)=(2s−1, s−s, 0)=(2s−1,0,0) by step 2.1, so y′(s)=(2,0,0) and Rx′(s)+v=R(1,2,0)+(0,1,0)=(2,−1,0)+(0,1,0)=(2,0,0); also x′′=0 and y′′=0, and R⋅0=0. This verifies y′=Rx′+v and y′′=Rx′′ in the worked case.

Remarks

  • The shift a=(−1,0,0) is the coordinate vector of o−o′=−e2 in the primed spatial basis. The change from u to u′=u+e1 supplies the boost v=(0,1,0); since b=0, it adds no constant shift.

  • The matrix R is the change-of-basis matrix from (e1,e2,e3) to (e1′,e2′,e3′), and its determinant is 1 because both bases are positively oriented. Had one of the two bases been negatively oriented, the same computation would give det⁡R=−1 and R∈O(3)∖SO(3).

ExampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

The rotating-frame acceleration identity for a uniform rotation about an axis

Example

Let ν>0 and define R:R→SO(3) by the rotation about the third standard axis through angle νt:

R(t)=(cos⁡νt−sin⁡νt0sin⁡νtcos⁡νt0001)(t∈R).

Then R(t)∈SO(3) for every t, and Ω=RTR′ is

Ω(t)=ν(0−10100000).

The angular-velocity vector of Angular velocity of a moving frame and the rotating-frame acceleration identity is ω(t)=(0,0,ν), independent of t. For any twice continuously differentiable y:R→R3, put x(t):=R(t)y(t) and let c:R→R3 be identically zero. That lemma gives

R(t)Tx′′(t)=y′′(t)+2 ω×y′(t)+ω×(ω×y(t))(x=Ry).

Facts & Assumptions

Given: A number ν>0 and the rotation matrix R(t) displayed above, defined for t∈R; write ct:=cos⁡νt and st:=sin⁡νt.

[L1]

The identity ct2+st2=1 holds for every real t; sine and cosine have derivatives (sin⁡)′=cos⁡ and (cos⁡)′=−sin⁡, so repeated differentiation gives continuous derivatives of all orders. Vector-valued formulas are differentiated componentwise (Sine and cosine defined by their real power series, Parity and the Pythagorean identity for sine and cosine, The derivatives of sine and cosine are cosine and minus sine, The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral).

[L2]

For twice continuously differentiable R:I→SO(3) and y, with Ω=RTR′ and ω(t) the unique vector with Ωz=ω×z for all z, Angular velocity of a moving frame and the rotating-frame acceleration identity gives RTx′′=y′′+2 ω×y′+ω′×y+ω×(ω×y)+RTc′′ when x=Ry+c.

[L3]

The cross product in R3 of the standard basis vectors is e1×e2=e3, e2×e3=e1, e3×e1=e2, and it is bilinear and alternating (The cross product in R3).

Verification

technique · direct; compute $\Omega$ and check the identity on two explicit curves defined on $\mathbb R$
1.1L1givenalgebra

By [L1], R′=ν(−st−ct0ct−st0000), and the product RTR′ equals ν(ctst0−stct0001)(−st−ct0ct−st0000)=ν(−ctst+stct−ct2−st20st2+ct2stct−ctst0000)=ν(0−10100000), using ct2+st2=1; hence reading off the entries of Ω as in clause 2 of [L2] gives ω(t)=(0,0,ν) and Ωz=ω×z for every z. Also RTR=I3 and det⁡R=ct2+st2=1, so R(t)∈SO(3); [L1] gives its required C2 regularity.

2.1step 1.1L1L2L3algebra

Take y(t)=e1 for t∈R, so y′=y′′=0, x(t)=R(t)e1=(ct,st,0) and x′′=−ν2(ct,st,0) by [L1]. Then RTx′′=−ν2RT(ct,st,0)=−ν2e1; the right-hand side is ω×(ω×e1)=(0,0,ν)×(0,ν,0)=(−ν2,0,0) by [L3] and step 1.1. The two sides agree.

3.1step 1.1step 2.1L1L2L3algebra∎

Take y(t)=te1 for t∈R, so y′=e1, y′′=0, x(t)=t(ct,st,0) and x′′=(−2νst−tν2ct, 2νct−tν2st, 0) by [L1]. Multiplying by RT gives RTx′′=(−tν2, 2ν, 0); the right-hand side is 2(0,ν,0)+t(−ν2,0,0)=(−tν2,2ν,0) by [L3] and step 1.1, using 2 ω×e1=2(0,ν,0) and ω×(ω×te1)=t(−ν2,0,0). The two sides agree, exhibiting the two rotational terms separately.

Remarks

  • Both checks are finite computations with the convention Ω=RTR′ and angular velocity ω=(0,0,ν). The other convention R′RT gives the same vector in this fixed-axis example, so the example does not test the distinction between conventions.

  • The identity is kinematic. The terms are named "Coriolis" and "centrifugal" only in a mechanical interpretation, which requires the force and mass hypotheses of the consuming physical items.

ExampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

The mass metric of a particle constrained to the unit sphere in spherical coordinates

Example

Let Q:=S2={x∈R3:∥x∥2=1} be the unit sphere (Euclidean spheres and closed balls as subspaces of Rn), with the smooth two-dimensional embedded structure supplied by A regular level set is an embedded submanifold applied to F(x)=⟨x,x⟩2 at the regular value 1 (verified below), let X:Q↪R3 be the inclusion, and let m>0 be the mass of one particle, so that the mass matrix of Constrained configurations, the mass matrix and the mass metric is M=mI3.

The mass metric of the inclusion is Gq(u,v)=m⟨u,v⟩2 for u,v∈TqS2, hence positive definite at every point by The mass metric is symmetric, positive semidefinite, and positive definite exactly when the configuration differential is injective. On D:=(0,π)×(0,2π) consider the parametrization

Ψ(θ,φ):=(sin⁡θcos⁡φ, sin⁡θsin⁡φ, cos⁡θ).

Its image is the open subset U:=S2∖{x∈S2:x2=0, x1≥0}, which excludes the poles and the longitude seam. The map Ψ:D→U is a diffeomorphism; in the manifold-to-coordinate convention the chart is (U,Ψ−1). Its differential is injective on D, and in the coordinate basis (∂θ,∂φ) the mass metric has the classical spherical form

G(θ˙∂θ+φ˙∂φ, θ˙∂θ+φ˙∂φ)=m(θ˙2+sin⁡2θ φ˙2).

Equivalently, the kinetic energy of a tangent vector with coordinate velocity (θ˙,φ˙) is T=m2(θ˙2+sin⁡2θ φ˙2).

Facts & Assumptions

Given: The unit sphere Q=S2, the inclusion X:Q↪R3, a mass m>0, and the parametrization Ψ(θ,φ)=(sin⁡θcos⁡φ,sin⁡θsin⁡φ,cos⁡θ) on the domain 0<θ<π, 0<φ<2π.

[L1]

The mass metric of (Q,X,m) is Gq(u,v)=m⟨u,v⟩2 for the inclusion and tangent vectors u,v, and it is positive definite (Constrained configurations, the mass matrix and the mass metric, The mass metric is symmetric, positive semidefinite, and positive definite exactly when the configuration differential is injective).

[L2]

Finite products and sums of smooth Euclidean maps are smooth, and continuous partial derivatives give the total derivative with its Jacobian matrix (Ck Euclidean maps and diffeomorphisms, Ck Euclidean maps are closed under componentwise algebra and composition, If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative). For a smooth map Ψ:D→S2, the chain rule gives d(i∘Ψ)=di∘dΨ, where the inclusion i:S2↪R3 has injective differential because it is a smooth embedding; hence independence of the ambient derivative columns implies that dΨ is injective (The differential of a smooth map, The chain rule for differentials of smooth maps, The inclusion of an embedded submanifold is a smooth embedding, Immersions, submersions, and constant-rank maps).

[L3]

The inner product is bilinear, cos⁡2+sin⁡2=1, and sin⁡′=cos⁡, cos⁡′=−sin⁡; repeated differentiation makes sine and cosine smooth (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, Parity and the Pythagorean identity for sine and cosine, The derivatives of sine and cosine are cosine and minus sine).

[L4]

A nonempty regular level set of a smooth map Mm→Nn is an embedded submanifold of dimension m−n (A regular level set is an embedded submanifold). An open subset of a smooth manifold has its restricted smooth structure (An open subset of a smooth manifold has a canonical restricted smooth structure); a map into an embedded submanifold is smooth if its ambient composite is smooth (Smoothness into an embedded submanifold is an initial property); and a smooth map with an isomorphic differential has a smooth local inverse (The smooth inverse function theorem on manifolds).

[L5]

The quarter-turn values are sin⁡(π/2)=1, cos⁡0=1, cos⁡π=−1 and sin⁡0=sin⁡π=0 (Quarter-turn values and shifts by pi/2 and pi, The derivatives of sine and cosine are cosine and minus sine). The sine monotonicity below therefore gives sin⁡θ>0 on (0,π). Cosine is strictly decreasing from 1 to −1 on [0,π], and φ↦(cos⁡φ,sin⁡φ) is a bijection from [0,2π) onto the unit circle (Signs, monotonicity intervals, and ranges of sine and cosine, t↦(cos⁡t,sin⁡t) is a bijection from [0,2π) onto the real unit circle).

Verification

technique · direct; differentiate the parametrization and evaluate the inner products of the coordinate vector fields
1.1givenL1L2L3L4algebra

The map F:R3→R, F(x)=⟨x,x⟩2, is smooth because it is a finite sum of coordinate squares and [L2] applies. For h∈R3, F(x+h)−F(x)−2⟨x,h⟩2=⟨h,h⟩2=∥h∥22, whose ratio to ∥h∥2 tends to 0; hence DFx(h)=2⟨x,h⟩2. At every x∈F−1(1), the vector h=x/2 has DFx(h)=1, so this differential is surjective. The fibre is nonempty since e3∈F−1(1); [L4] therefore makes S2=F−1(1) a smooth embedded submanifold of dimension 3−1=2. Its inclusion is a smooth embedding by [L2], so [L1] applies and gives Gq(u,v)=m⟨u,v⟩2 and positive definiteness at every point of the sphere.

2.1givenL2L3L5step 1.1algebra

The partial derivatives of Ψ are ∂θΨ=(cos⁡θcos⁡φ,cos⁡θsin⁡φ,−sin⁡θ) and ∂φΨ=(−sin⁡θsin⁡φ,sin⁡θcos⁡φ,0). For 0<θ<π these are linearly independent: a linear combination α∂θΨ+β∂φΨ=0 has third component −αsin⁡θ and first two components combining as cos⁡θ(αcos⁡φ,αsin⁡φ)+β(−sin⁡θsin⁡φ,sin⁡θcos⁡φ)=0; the first relation gives α=0 because sin⁡θ≠0, and then β=0 because sin⁡θ≠0. By [L2] this proves that dΨ is injective.

3.1step 1.1step 2.1L3L4L5algebra

The image of Ψ lies in S2 by [L3]. For x∈U, the exclusion of the poles gives −1<x3<1, so [L5] gives a unique θ∈(0,π) with cos⁡θ=x3 and sin⁡θ=1−x32>0. The pair (x1,x2)/sin⁡θ lies on the unit circle and is not (1,0); hence [L5] gives a unique φ∈(0,2π) producing x. Conversely an angle in that interval cannot give (1,0) by the same circle bijection. Thus Ψ:D→U is bijective. The removed meridian is relatively closed, so U is open in S2. Smoothness follows from [L3] and the embedded-target criterion [L4]. Its differential is an isomorphism between two-dimensional tangent spaces by steps 1.1 and 2.1, so [L4] supplies smooth local inverses everywhere; bijectivity makes these restrictions of its single global inverse. Thus Ψ is a diffeomorphism onto U, and (U,Ψ−1) is a smooth chart.

3.2step 2.1L1L3algebra

The three inner products of the coordinate vector fields are ⟨∂θΨ,∂θΨ⟩2=cos⁡2θcos⁡2φ+cos⁡2θsin⁡2φ+sin⁡2θ=1, ⟨∂θΨ,∂φΨ⟩2=−cos⁡θsin⁡θcos⁡φsin⁡φ+cos⁡θsin⁡θsin⁡φcos⁡φ=0 and ⟨∂φΨ,∂φΨ⟩2=sin⁡2θsin⁡2φ+sin⁡2θcos⁡2φ=sin⁡2θ, using cos⁡2+sin⁡2=1 for each angle.

4.1step 2.1step 3.1step 3.2L1algebra∎

Since Ψ is the parametrization of the inclusion chart, the mass metric in these coordinates is Gij=m⟨∂iΨ,∂jΨ⟩2; by step 3.2 its matrix is m(100sin⁡2θ), so for a coordinate velocity (θ˙,φ˙) the kinetic energy is T=12m(θ˙2+sin⁡2θ φ˙2), and the metric is positive definite on the chart domain because sin⁡θ≠0.

Remarks

  • The coordinate expression degenerates at θ=0 and θ=π, but this is a property of the chart, not of the mass metric: the chart domain also excludes the longitude seam and excludes the poles, where the parametrization is not an immersion. The metric of the inclusion is positive definite on the whole sphere, as recorded above, and any consumer working at a pole must use a different chart.

  • The sphere here is a mathematical constraint set. No claim is made that any physical particle is constrained to a sphere, or that the displayed kinetic energy is measured; a physical consumer states the constraint model and its scope in its own items.

CounterexampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

A collinear labeling admits a nontrivial rotation

Statement refuted

Every nonempty labeling S⊆R3 determines the rotation: if R∈SO(3) satisfies Rs=s for every s∈S, then R=I3.

Facts & Assumptions

Given: Write e1=(1,0,0) and e3=(0,0,1), where the tuple lists the values at the imported indices 0,1,2. Matrix-vector action is (Av)(k)=∑j<3A(k,j)v(j), the sole column of the rectangular product with Cv(j,0)=v(j). The set S={0,e3}⊆R3 and the matrix R=diag⁡(−1,−1,1).

[L2]

Transposition is involutive and reverses matrix products (Transpose is linear and involutive, and (AB)T=BTAT). The determinant is the finite Leibniz sum det⁡A=∑π∈S3sgn⁡(π)∏i<3A(π(i),i) (For n≥1, the determinant over a commutative ring by the Leibniz formula, and ∣det⁡A∣ for a real matrix).

Counterexample

technique · direct; exhibit a nontrivial rotation fixing a nonempty labeling
1.1givenalgebra

Take S={0,e3}⊆R3 and R=(−1000−10001). Then Re3=e3 and R0=0, while R≠I3 because Re1=−e1.

2.1step 1.1L1L2algebra∎

The matrix R lies in SO(3): its transpose equals itself, so RTR=R2=diag⁡(1,1,1)=I3, and in the Leibniz sum of [L2] every nonidentity permutation has an index i with π(i)≠i, giving a zero off-diagonal factor. Only the identity permutation contributes, with sign 1, so det⁡R=(−1)⋅(−1)⋅1=1.

Remarks

  • The set S={0,e3} is collinear: all its points lie on the line Re3. Its pointwise stabilizer is defined in A labeling containing two independent vectors has trivial SO(3)-stabilizer. The witness includes 0 because every linear map fixes 0; that lemma requires only two independent vectors in S and makes no separate origin assumption.

  • The displayed matrix fixes the axis Re3 and sends e1 to −e1; this verifies its nonidentity and its action on S directly. No physical-body claim is made.

ExampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

A closed field on a star-shaped domain with its radial primitive computed

Example

On U=R2, which is star-shaped about the origin (Star-shaped open subsets of Euclidean space), consider the C1 field F(x,y):=(2xy, x2) (Exact and closed C1 vector fields). Its potential

ϕ(x,y):=x2y

has gradient ∇ϕ=(2xy,x2)=F, so F is exact with potential ϕ. The closed field condition ∂1F0=∂0F1 holds as well, and the radial formula of the published Poincare lemma returns exactly ϕ:

∫01⟨F(t(x,y)),(x,y)⟩2 dt=∫013x2y t2 dt=x2y.

Facts & Assumptions

Given: The open set U=R2 with centre 0, the field F(x,y)=(2xy,x2) and the function ϕ(x,y)=x2y.

[L1]

R2 is open and star-shaped with respect to 0: for every x∈R2 and 0≤t≤1 the point tx lies in R2 (Star-shaped open subsets of Euclidean space).

[L2]

A C1 field F=(F0,F1) on an open subset of R2 is closed when ∂1F0=∂0F1, and exact when there is a C2 function ϕ with ∇ϕ=F (Exact and closed C1 vector fields).

[L3]

Poincare's lemma: on a star-shaped open domain U with centre a, every closed C1 field is exact, and a C2 potential is ϕ(x)=∫01⟨F(a+t(x−a)),x−a⟩2 dt (Poincare's lemma on a star-shaped domain: every closed C1 field is exact).

Verification

technique · direct; compute the partial derivatives and evaluate the radial integral
1.1givenL2algebra

The component functions F0(x,y)=2xy and F1(x,y)=x2 are C1, with ∂1F0=2x=∂0F1; hence F is closed by [L2].

2.1step 1.1L2algebra

The function ϕ(x,y)=x2y is C2 with ∂0ϕ=2xy=F0 and ∂1ϕ=x2=F1, so ∇ϕ=F and F is exact with potential ϕ by [L2].

3.1step 2.1L1L3L4algebra∎

For a=0 the radial formula of [L3] reads ϕrad(x,y)=∫01⟨F(tx,ty),(x,y)⟩2 dt=∫01(2t2x2y+t2x2y)dt=∫013x2yt2 dt. The integrand is a continuous polynomial on [0,1], hence Riemann integrable by [L4]. By [L4] applied to G(t):=x2yt3, whose derivative is 3x2yt2, the integral equals G(1)−G(0)=x2y, so ϕrad=ϕ.

Remarks

  • The example exercises the exactness side of Poincare's lemma on the whole plane. Exactness is not a consequence of closedness on domains with holes: the standard angular field on R2∖{0} is closed but has no global primitive, a result established separately in the published counterexample cex-the-closed-angular-form-on-the-punctured-plane-is-not-exact, which is not used as a premise here.

  • Both displays agree on the potential because the radial formula is an explicit construction, not an independent existence argument; the differentiation check of step 2.1 is what identifies ϕ as a potential.

ExampleConstruction: AI-adaptedVerification: AI-adaptedprecheck passOpen item page →

The stroboscopic map of the plane oscillator

Example

Let ω>0 and let F:R×R2→R2, F(t,z1,z2):=(z2,−ω2z1), the vector field of the undamped plane oscillator. Since F does not depend on t, it is P-periodic for every P>0; it is continuous and locally Lipschitz in the state variable, so The stroboscopic (time-P) map of a P-periodic first-order system applies for every P>0.

For z=(z1,z2)∈R2 the curve

xz(t)=(z1cos⁡ωt+z2ωsin⁡ωt, −ωz1sin⁡ωt+z2cos⁡ωt)

is defined on all of R and solves x′=F(t,x) with x(0)=z. Hence DP=R2 for every P>0, and the stroboscopic map is

θP(z)=(z1cos⁡ωP+z2ωsin⁡ωP, −ωz1sin⁡ωP+z2cos⁡ωP).

In particular, for P=2π/ω the map θP is the identity of R2, so every initial state is a fixed point; for P=π/ω it is θP(z)=−z, whose only fixed point is 0.

Facts & Assumptions

Given: A number ω>0 and the field F(t,z1,z2)=(z2,−ω2z1).

[L1]

The stroboscopic map of a continuous P-periodic field, locally Lipschitz in the state variable, is defined on DP={z:[0,P]⊆Jz} by θP(z)=xz(P), where xz is the unique maximal solution through (0,z); every solution through (0,z) is a restriction of xz (The stroboscopic (time-P) map of a P-periodic first-order system, Every Picard–Lindelöf initial value problem has one maximal solution on an open interval).

[L2]

The derivatives of sine and cosine are (sin⁡)′=cos⁡, (cos⁡)′=−sin⁡, with sin⁡0=0, cos⁡0=1. The shift identities sin⁡(x+π)=−sin⁡x, cos⁡(x+π)=−cos⁡x give sin⁡π=sin⁡2π=0, cos⁡π=−1, cos⁡2π=1 (Quarter-turn values and shifts by pi/2 and pi). Vector-valued curves are differentiated componentwise (Sine and cosine defined by their real power series, The derivatives of sine and cosine are cosine and minus sine, The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral).

[L3]

The system is autonomous, so F(t,z)=F(0,z) for all t. For z,w∈R2,

∥F(t,z)−F(t,w)∥22=(z2−w2)2+ω4(z1−w1)2≤max⁡{1,ω4}∥z−w∥22,

so F is globally Lipschitz in the state variable with constant max⁡{1,ω2} and is continuous; it satisfies the Picard-Lindelof hypotheses (Every Picard–Lindelöf initial value problem has one maximal solution on an open interval, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

Verification

technique · direct; verify the explicit solution and evaluate the time-$P$ map
1.1givenL2algebra

Differentiating the two components of the displayed xz with respect to t using [L2] gives xz′(t)=(−ωz1sin⁡ωt+z2cos⁡ωt, −ω2z1cos⁡ωt−ωz2sin⁡ωt), whose first component is the second component of xz and whose second component is −ω2 times the first component of xz; hence xz′=F(t,xz) on all of R, and xz(0)=(z1,z2)=z because sin⁡0=0 and cos⁡0=1.

2.1step 1.1L1L3

The curve of step 1.1 is a solution on all of R; by [L1] and [L3], every solution of the initial value problem with data (0,z) is a restriction of the unique maximal solution xz, whose domain is an open interval. Since a solution with domain R exists, that maximal domain is all of R; hence [0,P]⊆Jz for every P and DP=R2. Substituting t=P in step 1.1 gives the displayed formula for θP(z) for every P>0.

3.1step 2.1L2algebra∎

If P=2π/ω then cos⁡ωP=1 and sin⁡ωP=0, so θP(z)=(z1,z2)=z for every z; if P=π/ω then cos⁡ωP=−1 and sin⁡ωP=0, so θP(z)=(−z1,−z2)=−z, which equals z only for z=0. This proves the two special cases.

Remarks

  • The example verifies the iterate identity of Iterates of the stroboscopic map, and equilibria as fixed points implicitly: since DP=R2 here, every iterate is defined and θPk(z)=xz(kP); for P=2π/ω this makes every iterate the identity.

  • The map θP is a linear map of R2 and preserves the area element; no symplectic or stability conclusion is drawn here, since those notions belong to the symplectic items and to the stability definitions of this collection.

Sources