Alphabeta Math
LemmaStatement: 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.

Exact free Gaussian propagation

Statement

Assume the Axiom of Countable Choice and the Axiom of Choice (The Axiom of Countable Choice (ACω), The Axiom of Choice). Fix m>0, ℏ>0, σ>0 and real x0,p0, and let

ψ0(x):=(2πσ2)−1/4exp⁡(−(x−x0)24σ2)exp⁡(ip0(x−x0)ℏ),x∈R.

Let H0 be the free operator of clause (b) of Real multiplication operators and free plus bounded-potential Schrödinger operators with c:=2π2ℏ2/m and b:=0, i.e.

H0:=F2∗ M2π2ℏ2ξ2/m F2,D(H0)={f∈L2(R):ξ2F2f∈L2(R)},

which is self-adjoint with D(H0)=W2,2(R;C) and H0f=−ℏ22mΔf for f∈S(R). Let U be the momentum transform of Momentum Fourier Jacobian: the unitary momentum representation, the surjective isometry L2(R,dx)→L2(R,dp) with U=Dℏ∘F2 on S(R). For t∈R put τ:=ℏt2mσ2 and

ψt(x):=(2πσ2)−1/4(1+iτ)−1/2exp⁡(−(x−x0−p0t/m)24σ2(1+iτ)+ip0(x−x0)ℏ−ip02t2mℏ),

where z−1/2:=exp⁡(−12Log⁡z) is the principal branch of Complex logarithms, the principal logarithm, and principal and multivalued complex powers, legitimate at z=1+iτ because Re⁡(1+iτ)=1>0. Let Qf(x)=xf(x) and Pf(x)=−iℏf′(x) act on S(R). Then:

(a) ψ0∈S(R), ∥ψ0∥2=1, and for every t one has ψt∈S(R)⊆D(H0); the map t↦ψt is continuously differentiable R→L2(R) with ψt∈D(H0) and ddtψt=−iℏH0ψt, and consequently ψt=e−itH0/ℏψ0, where e−itH0/ℏ:=eitT with T:=−H0/ℏ denotes the unitary group generated by the self-adjoint operator T (A self-adjoint operator generates a strongly continuous unitary group).

(b) Position density. ∣ψt(x)∣2=(2πσ2(1+τ2))−1/2exp⁡(−(x−x0−p0t/m)2/(2σ2(1+τ2))); it is a probability density of mean x0+p0t/m and variance σ2(1+τ2).

(c) Momentum density. ∣Uψt(p)∣2=2/π (σ/ℏ)exp⁡(−2σ2(p−p0)2/ℏ2)=∣Uψ0(p)∣2 for almost every p; it is a probability density of mean p0 and variance ℏ2/(4σ2), independent of t.

(d) Standard deviations. With ΔQ:=∥(Q−⟨Qψt,ψt⟩)ψt∥2 and ΔP:=∥(P−⟨Pψt,ψt⟩)ψt∥2 one has ΔQ=σ1+τ2, ΔP=ℏ/(2σ) and

ΔQ ΔP=ℏ21+τ2 ≥ ℏ2,

with equality exactly at t=0; the lower bound ℏ/2 is the one supplied by clause (d) of Weak-form Robertson uncertainty on a common domain with c=ℏ.

(e) Scope. These are exact conditional calculations for the free model defined by H0; no measurement, instrument or laboratory preparation is described, and nothing here asserts that a physical wave packet is Gaussian.

Facts & Assumptions

Given: m,ℏ,σ>0, x0,p0∈R, the vector ψ0 above, the operator H0 and the transform U.

[F1]

S(R) consists of the C∞ functions all of whose polynomial-weighted derivatives are bounded; every Schwartz function lies in L1∩L2 (Schwartz space and its seminorms).

[F2]

U:L2(R,dx)→L2(R,dp) is a surjective complex-linear isometry with U=DℏF2 on S(R), where Dℏg(p)=(2πℏ)−1/2g(p/(2πℏ)); being a surjective isometry, U preserves inner products and U∗=U−1 (Momentum Fourier Jacobian: the unitary momentum representation).

[F3]

If Re⁡z>0 then ∫Re−zx2dx=π z−1/2 for the branch z−1/2=exp⁡(−12Log⁡z), and for G1(z):=∫Re−zx2dx one has G1′(z)=−∫Rx2e−zx2dx (Complex Gaussian continuation(a),(b)).

[F4]

H0 is self-adjoint, D(H0)=W2,2(R;C)={f∈L2:ξ2F2f∈L2} with equivalent norms on this set, and H0f=−ℏ22mΔf for f∈S(R) (Real multiplication operators and free plus bounded-potential Schrödinger operators(b)). Here F2 is the Plancherel transform extending the L1 transform f^(ξ)=∫f(x)e−2πixξdx; it is a surjective isometry preserving inner products (Plancherel theorem, Fourier transform on complex L1 classes).

[F5]

For a self-adjoint operator T the operators V(t):=eitT form a strongly continuous one-parameter unitary group with V(t)D(T)=D(T), TV(t)=V(t)T on D(T), and lim⁡s→01s(V(s)x−x)=iTx exactly for x∈D(T) (A self-adjoint operator generates a strongly continuous unitary group).

[F6]

Dominated convergence: if fn→f pointwise and ∣fn∣≤g for a fixed integrable g, then ∫fn→∫f (Dominated convergence). A continuous real function on an interval has a primitive and ∫abG′=G(b)−G(a) for it (Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G); the integral over R is the limit of the compact-interval integrals when it exists (Improper integrals over unbounded intervals); (fg)′=f′g+fg′ (Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0); for differentiable u,v with integrable derivatives on [a,b], ∫abuv′=u(b)v(b)−u(a)v(a)−∫abu′v (If u,v are differentiable on [a,b] with u′,v′ integrable, then ∫abuv′=u(b)v(b)−u(a)v(a)−∫abu′v). Substitutions: an invertible linear T:R→R scales Lebesgue measure by ∣det⁡T∣ (A linear map T of Rn sends Lebesgue measurable sets to Lebesgue measurable sets, with λn(T[E])=∣det⁡T∣ λn(E) when T is invertible and T[E] Lebesgue null when it is not), and λ1(E+h)=λ1(E) for every Lebesgue measurable E and h∈R (Lebesgue outer measure, Lebesgue measurability and Lebesgue measure are unchanged by translation); hence for every nonnegative measurable H and every real x0, ∫RH(x−x0)dx=∫RH(y)dy, by approximation with simple functions and monotone convergence (Monotone convergence for the integral).

[F7]

For every unit f∈S(R) and every c>0, with Qf(x)=xf(x) and Pcf(x)=−icf′(x), one has ΔQ ΔPc≥c/2, the differences being centered norms (Weak-form Robertson uncertainty on a common domain(d)).

Proof

1.1F1F3F6algebra

[The initial packet is a unit Schwartz vector] The derivative ∂xkψ0 equals Pk(x)ψ0(x) for a polynomial Pk: this follows by induction with the product rule [F6] applied to the three factors (2πσ2)−1/4, e−(x−x0)2/(4σ2) and the bounded factor eip0(x−x0)/ℏ, each differentiation of the last factor multiplying it by the constant ip0/ℏ. Since x↦xjPk(x)e−(x−x0)2/(4σ2) is bounded for all j,k (the exponential decays faster than any polynomial), ψ0∈S(R) [F1]. By [F6] ∫R∣ψ0(x)∣2dx=(2πσ2)−1/2∫Re−y2/(2σ2)dy=(2πσ2)−1/2π (2σ2)1/2=1, using [F3] with z=1/(2σ2) and the substitution x=x0+y.

1.2F3F6algebra

[The complex shift identity] For Re⁡a>0 and b∈R one has ∫Re−ay2eibydy=π/a e−b2/(4a) with π/a:=π a−1/2: setting F(b):=∫Re−ay2+ibydy (absolutely convergent since ∣e−ay2+iby∣=e−Re⁡(a)y2 is integrable) and GR(b):=∫−RRe−ay2+ibydy, the fundamental theorem and [F6] give ddbGR(b)=i∫−RRye−ay2+ibydy, and since ∣y∣e−Re⁡(a)y2 is integrable, GR→F pointwise with the difference quotients converging, so F is differentiable with F′(b)=i∫Rye−ay2+ibydy; integrating ∂y(e−ay2+iby)=(−2ay+ib)e−ay2+iby over [−R,R] gives ∫−RR(−2ay+ib)e−ay2+ibydy=e−aR2+ibR−e−aR2−ibR, of modulus at most 2e−Re⁡(a)R2→0, so [F6] as R→∞ gives −2aFy(b)+ibF(b)=0 with Fy(b):=∫Rye−ay2+ibydy, hence Fy=ib2aF and F′=−b2aF; the product rule [F6] gives ddb[F(b)eb2/(4a)]=0, so the product is constant (fundamental theorem applied to real and imaginary parts [F6]) and its value at b=0 is F(0)=∫Re−ay2dy=π a−1/2 by [F3].

1.3F2F4algebra

[The free operator in momentum coordinates] Let B be the multiplication operator (Bφ)(p)=p22mφ(p) on D(B):={φ∈L2(dp):p2φ∈L2}; then U∗BU=H0 with D(H0)=U∗D(B), and H0U∗=U∗B, because Dℏ∗=Dℏ−1 and (Dℏ−1MgDℏh)(ξ)=g(2πℏξ)h(ξ) for every Borel g and every h in the domain of Dℏ−1MgDℏ, so with g(p)=p2/(2m) one has U∗BU=F2∗Dℏ∗Mp2/(2m)DℏF2=F2∗M2π2ℏ2ξ2/mF2=H0 as operators [F2, F4]; for the domains, f∈D(H0) means ξ2F2f∈L2, which by Uf(p)=(2πℏ)−1/2F2f(p/(2πℏ)) is equivalent to p2Uf∈L2 (the substitution ξ=p/(2πℏ) and the constants (2πℏ)±2 preserve L2-membership), i.e. to Uf∈D(B), so D(H0)=U∗D(B), the identity U∗BU=H0 holds on it, and multiplying on the right by U∗ gives H0U∗=U∗BUU∗=U∗B.

2.1step 1.1step 1.2F2F4F6algebra

[Fourier transform of the initial packet] Since ψ0∈L1∩L2 [step 1.1], F2ψ0 is computed by the L1 integral [F4]; substituting x=x0+y [F6] and using ψ0(x0+y)=(2πσ2)−1/4e−y2/(4σ2)eip0y/ℏ gives F2ψ0(ξ)=(2πσ2)−1/4e−2πix0ξ∫Re−y2/(4σ2)ei(p0/ℏ−2πξ)ydy=(2πσ2)−1/4e−2πix0ξ2σπ e−σ2(p0/ℏ−2πξ)2 by step 1.2 with a=1/(4σ2) and b=p0/ℏ−2πξ; since ψ0∈S(R), Uψ0=DℏF2ψ0 [F2], so φ0:=Uψ0 satisfies φ0(p)=(2πℏ)−1/2F2ψ0(p/(2πℏ))=Ke−ix0p/ℏe−σ2(p−p0)2/ℏ2 with K:=(2πℏ)−1/2(2πσ2)−1/42σπ and K2=2/π σ/ℏ.

2.2step 1.1F1F4algebra

[The packet is Schwartz and its modulus is the spreading Gaussian] Putting zt:=1/(1+iτ) (so Re⁡zt=1/(1+τ2)>0) and ct:=x0+p0t/m, one has ∣ψt(x)∣=(2πσ2)−1/4∣1+iτ∣−1/2e−Re⁡(zt)(x−ct)2/(4σ2)=(2πσ2)−1/4(1+τ2)−1/4e−(x−ct)2/(4σ2(1+τ2)); the remaining x-dependent phase is exp⁡(−iIm⁡(zt)(x−ct)2/(4σ2)+ip0(x−x0)/ℏ), which has modulus one. By the product rule, each x-derivative of ψt is a polynomial in x times ψt, so the displayed decaying modulus makes every polynomial-weighted derivative bounded and ψt∈S(R); in particular ψt∈L1∩L2 and ψt∈D(H0) by [F4].

3.1step 1.2step 2.1step 2.2F2F4F6algebra

[Momentum coordinates of the packet] With At:=(2πσ2)−1/4(1+iτ)−1/2, at:=1/(4σ2(1+iτ)) and bξ:=p0/ℏ−2πξ, substituting x=ct+y [F6] in the L1 integral [F4] of the packet ψt=Ate−ip02t/(2mℏ)e−at(x−ct)2ei(p0/ℏ)(x−x0) [step 2.2] gives F2ψt(ξ)=Atπ/at eip02t/(2mℏ)e−2πiξcte−bξ2/(4at) by step 1.2 with a=at (legal since Re⁡at>0); since ψt∈S(R), Uψt=DℏF2ψt [F2], and with p=2πℏξ, π/at=2σπ(1+iτ)1/2, At(1+iτ)1/2=(2πσ2)−1/4 and 1/(4at)=σ2(1+iτ) one gets Uψt(p)=Keip02t/(2mℏ)e−ipct/ℏe−σ2(p−p0)2/ℏ2e−it(p−p0)2/(2mℏ) with K as in step 2.1.

3.2step 2.2F3F6algebra

[Position density, mean and variance] By step 2.2, ∣ψt(x)∣2=(2πσ2(1+τ2))−1/2e−(x−ct)2/(2dt2) with ct=x0+p0t/m and dt:=σ1+τ2; substituting x=ct+dtu and using ∫Re−u2/2du=2π (that is [F3] with z=1/2) gives ∫R∣ψt∣2=1, while the centering ∫R(x−ct)∣ψt(x)∣2dx=dt∫Rue−u2/22πdu=0 and ∫R(x−ct)2∣ψt(x)∣2dx=dt2∫Ru2e−u2/22πdu=dt2, since ∫Ru2e−u2/2du=−G1′(1/2)=2π by [F3]; hence the density has mean ct=x0+p0t/m and variance σ2(1+τ2), which is clause (b).

4.1step 2.1step 3.1F2algebra

[The packet is the free evolution of the momentum Gaussian] With θt(p):=tp2/(2mℏ) and φt:=e−iθtφ0, steps 2.1 and 3.1 give Uψt(p)=e−itp2/(2mℏ)φ0(p)=φt(p) for every p, because the two exponents differ by i(p02t/(2mℏ)−p(x0+p0t/m)/ℏ−t(p−p0)2/(2mℏ)+px0/ℏ+tp2/(2mℏ))=it(p02−2pp0−(p−p0)2+p2)/(2mℏ)=0; hence Uψt=φt and, U being a surjective isometry [F2], ψt=U∗φt.

4.2step 2.1step 3.1step 3.2F3F6algebra

[Momentum density, mean and variance] By steps 2.1 and 3.1, ∣Uψt(p)∣2=∣φ0(p)∣2=K2e−2σ2(p−p0)2/ℏ2=2/πσℏe−(p−p0)2/(2w2) with w:=ℏ/(2σ) and K2=2/π σ/ℏ; substituting p=p0+wu as in step 3.2 gives total mass 1, mean p0 and variance w2=ℏ2/(4σ2), and the density does not depend on t; this is clause (c).

5.1step 1.3step 2.1step 4.1F2F6algebra

[The momentum path solves the momentum equation] Put q(p):=p2/(2mℏ) and φt(p):=e−itq(p)φ0(p). The Gaussian formula of step 2.1 gives qφ0∈L2(dp) and p2φ0∈L2(dp), so φt∈D(B) for every t. For fixed t and h≠0, φt+h(p)−φt(p)h=φt(p)e−ihq(p)−1h tends pointwise to −iq(p)φt(p)=−(i/ℏ)p22mφt(p). Since ∣e−is−1∣≤∣s∣, the difference between the quotient and its limit has squared modulus at most 4∣q(p)φ0(p)∣2, an integrable function; dominated convergence [F6] therefore gives convergence in L2(dp). Also φt+h→φt and Bφt+h→Bφt in L2(dp) by pointwise convergence and the respective integrable square majorants 4∣φ0∣2 and 4∣Bφ0∣2. Thus t↦φt is continuously differentiable with φt′=−(i/ℏ)Bφt. Consequently ψt=U∗φt [step 4.1] is C1 in L2, lies in D(H0)=U∗D(B), and ddtψt=U∗φt′=−iℏU∗Bφt=−iℏH0U∗φt=−iℏH0ψt, using that U∗ is an isometry [F2] and H0U∗=U∗B from step 1.3; at t=0 the formula of the Statement reduces to ψ0.

6.1step 1.1step 5.1F4F5F6algebra

[Uniqueness of the Cauchy problem] Let V(t):=eitT with T:=−H0/ℏ [F5]; since ψ0∈D(H0)=D(T) [step 1.1, F4], one has V(t)ψ0∈D(T) for all t and t↦V(t)ψ0 is C1 with ddtV(t)ψ0=iTV(t)ψ0=−iℏH0V(t)ψ0, because [F5] gives 1h(V(t+h)−V(t))ψ0=V(t)1h(V(h)ψ0−ψ0)→V(t)(iTψ0) as ψ0∈D(T) and V(t) is bounded, and TV(t)=V(t)T on D(T) by [F5]; by step 5.1, wt:=ψt−V(t)ψ0 is C1 with wt∈D(H0) and wt′=−iℏH0wt, and w0=0 since V(0)=I and the Statement's formula at t=0 is the initial packet; for the real function t↦∥wt∥2 one has ∥wt+h∥2−∥wt∥2h=2Re⁡⟨wt,wt+h−wth⟩+∥wt+h−wt∥2h→2Re⁡⟨wt,wt′⟩, the last term tending to 0; now, since the inner product is conjugate-linear in its second argument, ⟨wt,wt′⟩=⟨wt,−iℏH0wt⟩=iℏ⟨wt,H0wt⟩. The number ⟨wt,H0wt⟩=⟨H0wt,wt⟩ is real because H0 is self-adjoint, hence symmetric, on D(H0) [F4], so Re⁡⟨wt,wt′⟩=0 and ∥wt∥2 has vanishing derivative, hence is constant by the fundamental theorem [F6] and equals ∥w0∥2=0; thus ψt=V(t)ψ0=e−itH0/ℏψ0 for every t, which is clause (a).

7.1step 3.2step 4.2step 6.1F2F6algebra

[Operator standard deviations] For ψ∈S(R) one has U(Pψ)=MpUψ with (Mpg)(p)=pg(p), because F2(Pψ)(ξ)=∫R(−iℏψ′(x))e−2πixξdx=2πℏξ F2ψ(ξ) by integration by parts, the boundary terms [ψ(x)e−2πixξ] vanishing at ±∞ as ψ∈S(R) [F6], and applying Dℏ with p=2πℏξ gives the claim [F2]; the unit vector ψt of step 6.1 is the packet of steps 3.2 and 4.2, and for it this gives bt:=⟨Pψt,ψt⟩=⟨UPψt,Uψt⟩=∫Rp∣Uψt(p)∣2dp=p0, a real number, and ΔP:=∥(P−bt)ψt∥2=∥U(P−bt)ψt∥2=(∫R(p−p0)2∣Uψt(p)∣2dp)1/2=ℏ/(2σ) by step 4.2, while at:=⟨Qψt,ψt⟩=∫Rx∣ψt(x)∣2dx=ct and ΔQ:=∥(Q−at)ψt∥2=(∫R(x−ct)2∣ψt(x)∣2dx)1/2=σ1+τ2 by step 3.2; hence ΔQΔP=ℏ21+τ2.

8.1step 7.1F7algebra∎

[Robertson comparison and scope] By [F7] with c=ℏ, applied to the unit Schwartz vector ψt of step 7.1, one has ΔQΔP≥ℏ/2; the exact value of step 7.1 is ℏ21+τ2, which equals ℏ/2 exactly when τ=0, i.e. exactly at t=0 since τ=ℏt/(2mσ2), and is strictly larger otherwise; thus clause (d) holds, and clause (e) records that all of the above are computations inside the free model defined by H0, with no measurement interpretation.

Depends on

Dependency tree · 2 levels

Sources