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.

Gaussian and gamma product tools

1 · Prerequisites

2 · Summary

This page supplies finite Gaussian and gamma product identities, a Gaussian covariance calculation under the stated Countable Choice measure hypothesis, the normalized marginal of a specified spherical-shell preparation under Countable Choice, and fine-entropy conservation under supplied smooth volume-preserving flow data. Physical results then calculate canonical and grand-canonical ideal gases, a finite Ising ring, a one-dimensional hard-rod gas, and a finite spherical-shell gas. Each physical conclusion is conditional on its explicitly adopted preparation, reference measure, units and finite-domain assumptions; the continuous Lebesgue product and Jordan-volume calculations state their Countable Choice conditions. No experiment, equilibration mechanism, or general ensemble-equivalence claim is asserted. Moment differentiation follows the explicitly evaluated moments and Gamma recurrence; finite-ring transfer identities use the published periodic Ising-chain theorem. The continuous models share the explicitly adopted positive action reference, with their original label-counting conventions retained. Draft status and the ordinary independent review obligations remain.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

LemmaStatement: AI-adaptedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Gaussian, gamma and beta product integrals

Statement

For a>0 and an integer r≥0, define Ir(a):=∫Rx2re−ax2 dx as the two-sided improper Riemann integral. Then

Ir(a)=Γ(r+12)ar+1/2,∫Rx2r+1e−ax2 dx=0,

where Γ(s):=∫0∞ts−1e−t dt for s>0. In particular, I0(a)=π/a and I1(a)=π/(2a3/2). For every compact interval K=[a0,a1]⊂(0,∞), Ir is continuously differentiable on K and

Ir′(a)=−Ir+1(a).

For s,t>0, put B(s,t):=∫01us−1(1−u)t−1 du. Then B(s,t) is finite and

Γ(s)Γ(t)=Γ(s+t)B(s,t),Γ(s+1)=sΓ(s),Γ(12)=π.

Consequently, if independent positive random variables have densities xs−1e−x/Γ(s) and xt−1e−x/Γ(t), their sum has density xs+t−1e−x/Γ(s+t).

Facts & Assumptions

Given: a>0, an integer r≥0, and positive real parameters s,t.

[F1]

The improper integral of e−x2 over R is finite and positive: the improper Gaussian integral is finite and positive.

[F2]

The square of the one-dimensional Gaussian integral is the plane Gaussian integral: the square of the one-dimensional Gaussian integral equals the nonnegative plane Gaussian integral.

[F3]

The plane Gaussian integral equals π by polar coordinates: the nonnegative plane Gaussian integral equals π, with the polar substitutions justified on compact half-annuli away from the origin.

[F4]

Change of variables for an injective C1 map on a compact Jordan set: an injective C1 map with invertible derivative on an open neighborhood of a compact Jordan set transforms a bounded integrable function by the absolute Jacobian determinant.

[F5]

Continuity and derivatives of positive-base real powers: for any fixed real exponent γ, the positive-base power aγ is differentiable on (0,∞), with derivative γaγ−1; these positive-base powers are continuous.

[F6]

Measures agreeing on a generating pi-system are equal under an increasing finite-measure exhaustion from that pi-system: two measures agreeing on a generating pi-system are equal if that system contains an increasing common finite-measure exhaustion of the whole space.

Proof

1.1givenalgebra

For s>0, the gamma integrand on (0,1) is bounded by ts−1, whose integral is 1/s. For t≥1, ts−1e−t/2 is bounded, so ts−1e−t≤Cse−t/2 and the tail converges. On a compact positive s-interval, the lower endpoint bounds the near-zero power and the upper endpoint bounds the exponential tail uniformly. Hence Γ(s) is finite and has controlled endpoint tails.

1.2F2F3algebra

Let I=∫Re−x2dx. By [F2], I2 equals the nonnegative plane Gaussian integral, and [F3] evaluates that integral as π. Hence I2=π.

2.1F1step 1.2algebra

By [F1], I is finite and positive. The positive square root in step 1.2 is therefore I=π.

2.2step 1.1algebra

Integration by parts on [ε,R] gives ∫εRtse−t dt=[−tse−t]εR+s∫εRts−1e−t dt. The boundary terms vanish by step 1.1, so Γ(s+1)=sΓ(s); direct integration gives Γ(1)=1.

2.3step 1.1algebra

For s,t>0, the beta integrand is integrable near 0 by comparison with us−1 and near 1 by comparison with (1−u)t−1. On each compact rectangle [1/j,j]2 the continuous product integrand xs−1yt−1e−(x+y) has integral equal to the product of its one-dimensional integrals, as follows by factoring the finite rectangular Riemann sums. These nested rectangles exhaust the positive quadrant and their integrals tend to Γ(s)Γ(t) by step 1.1.

3.1step 1.2step 2.1algebra

Substitution t=x2 on compact half-intervals and passage to the convergent improper limits gives Γ(12)=2∫0∞e−x2 dx=I=π by steps 1.2 and 2.1.

3.2step 1.1step 2.3F4algebra

On Aj=[1/j,j]×[1/j,1−1/j] in coordinates (r,u), the map Φ(r,u)=(ru,r(1−u)) is injective, has determinant −r≠0, and is a C1 diffeomorphism on an open neighborhood. Its image is a compact Jordan set. The images exhaust the positive quadrant cofinally with the rectangles in step 2.3: each compact rectangle stays a positive distance from both axes and has bounded x+y, while each image is compact and stays away from the axes. Applying [F4] on Aj and taking the increasing limits gives Γ(s)Γ(t)=∫0∞rs+t−1e−r dr B(s,t)=Γ(s+t)B(s,t). The limits are finite by step 1.1 and the endpoint bounds in step 2.3.

4.1givenF1F4step 1.1step 1.2step 2.1step 2.2step 3.1algebra

For integer r≥0, x2re−ax2/2 is bounded on R, so x2re−ax2≤Ce−ax2/2 and [F1] gives absolute convergence; odd moments converge absolutely and vanish by reflection. On compact intervals the substitution y=a x is covered by [F4], and its improper limit gives I0(a)=a−1/2I0(1)=π/a by steps 1.2 and 2.1. The substitution t=ax2 on both half-lines gives Ir(a)=a−r−1/2Γ(r+12); finiteness follows from step 1.1, the half-integer starting value from step 3.1, and recurrence from step 2.2.

4.2givenstep 1.1step 3.2F4algebra

Let C=[c,d]⊂(0,∞) and set Qj=[1/j,j]2. The event Aj=Qj∩{(x,y):x+y∈C} is Jordan measurable. On each grid rectangle I×J⊂Qj, independence gives probability P(X∈I)P(Y∈J)=(∫Ifs)(∫Jft); because the marginal densities are continuous on Qj, the finite rectangular Riemann sums factor and converge to the integral of fs(x)ft(y) on that rectangle. The inner and outer grid unions for Aj differ only in cells meeting its finitely many line-segment boundary pieces; their total area tends to zero, and fsft is bounded on Qj. Thus P((X,Y)∈Aj)=∫Ajfsft dx dy. As Qj increases to the positive quadrant, countable additivity gives probabilities increasing to P(X+Y∈C), and the nonnegative improper integrals increase to the integral over x,y>0. The images of C×[1/j,1−1/j] under (r,u)↦(ru,r(1−u)) increase to AC={x,y>0:x+y∈C}. On each such compact set the map is injective with absolute Jacobian r, so [F4] transforms the integral to rs+t−1e−rus−1(1−u)t−1/(Γ(s)Γ(t)). Its rectangular Riemann sums factor; as j→∞ the u-integral tends to B(s,t). By step 3.2, P(X+Y∈C)=∫Crs+t−1e−r/Γ(s+t) dr. A singleton C has probability zero by the same zero-area grid-boundary argument. The density integrates to one by step 1.1.

5.1F5step 2.2step 4.1algebra

Fix the integer r≥0. Step 4.1 evaluated Ir(a)=Γ(r+12)a−r−1/2 for every a>0 using only convergent improper substitutions, and step 2.2 proved the Gamma recurrence by integration by parts; neither argument used parameter differentiation. By [F5], differentiation of this established explicit function gives Ir′(a)=−(r+12)Γ(r+12)a−r−3/2=−Γ(r+32)a−r−3/2=−Ir+1(a). The same positive-base power formula makes this derivative continuous on (0,∞), hence on a neighborhood of every compact K⊂(0,∞). This proves the full parameter clause without interchanging a derivative and an improper integral.

6.1F6step 4.2construct∎

Let μ be the law of X+Y on (0,∞) and let ν be the probability measure with the density in step 4.2. The collection P={∅}∪{[a,b]:0<a≤b<∞} is closed under intersections, allowing a singleton or the empty set. It generates the Borel sigma-algebra of the positive half-line: its members are closed there, and every relatively open interval is a countable union of closed subintervals with positive endpoints; unbounded intervals are obtained by an additional increasing union. Step 4.2 gives equality of μ and ν on every member of P, including the singleton case. The sets Pj=[1/j,j] for integers j≥1 increase to (0,∞), lie in P, and have common measure at most 1. Thus [F6] proves μ=ν on all Borel sets, which establishes the claimed gamma convolution law.

Remarks

The parameter-derivative clause follows from the independently evaluated moments and Gamma recurrence. No research-module token is used as a logical supplier, and these integral identities make no empirical claim. No claim is made for a≤0, s≤0 or t≤0.

LemmaStatement: AI-adaptedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Jacobian identity for a specified smooth diffeomorphic flow

Statement

Let n≥1, let I⊆R be an open interval containing 0, and let U⊆Rn be open. Suppose χ:I×U→Rn is jointly C2, each χt(a):=χ(t,a) is a C2 diffeomorphism from U onto an open set Ut, and χ0(a)=a. Let D:={(t,x):t∈I, x∈Ut} be open and let b:D→Rn be C1 with ∂tχ(t,a)=b(t,χ(t,a)). Define F(t,a):=Daχ(t,a), J(t,a):=det⁡F(t,a), and div⁡b:=∑i=1n∂xibi. Then

∂tF=(Dxb)(t,χ)F,∂tJ=J(div⁡b)(t,χ),J(0,a)=1.

In particular J(t,a)>0 for all (t,a)∈I×U, and if div⁡b=0 on D, then J(t,a)=1. The same identities hold on the flat torus Tn:=Rn/Λ, where Λ=Zℓ1+⋯+Zℓn for a linearly independent basis ℓ1,…,ℓn and the smooth structure and Euclidean metric descend from Rn. In that clause, let χ:I×Tn→Tn be jointly C2, with each χt a diffeomorphism and χ0=id⁡. Let π:TTn→Tn be the tangent-bundle projection, and let b:I×Tn→TTn have C1 chart coefficients and satisfy π(b(t,a))=a and ∂tχ(t,a)=b(t,χ(t,a)). In every translation chart write b=∑ibi∂xi and define div⁡b:=∑i∂xibi. Define F(t,a) to be the matrix of Daχt:TaTn→Tχt(a)Tn in the translation frames and J(t,a):=det⁡F(t,a). Translation-chart changes have derivative I, so these matrices, determinants and divergences agree on overlaps. These are identities for supplied flow data; no existence or global-flow theorem is asserted.

Facts & Assumptions

Given: The interval I, open set U, jointly C2 family χ, diffeomorphism and initial-identity hypotheses, and C1 velocity field b in the Statement.

[F1]

Clairaut--Schwarz theorem for continuous second partial derivatives: the two mixed second partial derivatives of a C2 scalar function on an open Euclidean set agree.

[F2]

The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a): for totally differentiable maps, the derivative of a composite is the composite of the derivatives.

[F3]

The derivative of the determinant along an invertible C¹ matrix curve: (det⁡A)′=det⁡A  tr⁡(A−1A′): if A(t) is a differentiable invertible real square-matrix curve, then (det⁡A)′=det⁡A tr⁡(A−1A′).

[F4]

For A∈Mm×n(F) and B∈Mn×m(F), tr⁡(AB)=tr⁡(BA): for matching finite real square matrices, tr⁡(AB)=tr⁡(BA).

[F5]

The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a): a function continuous on a closed interval and differentiable on its interior has an interior point where its derivative equals the interval's average slope.

Proof

1.1givenF1F2algebra

For each component χi, joint C2 regularity and [F1] give ∂t∂ajχi=∂aj∂tχi. Differentiating the given velocity equation in a and applying [F2] yields ∂tF=LF, where L(t,a):=(Dxb)(t,χ(t,a)).

2.1givenF2F3step 1.1

Each F(t,a) is invertible: differentiating χt−1∘χt=id⁡ by [F2] gives Dχt−1(χt(a))F(t,a)=I. Thus [F3] applies to the C1 matrix curve t↦F(t,a) and gives ∂tJ=Jtr⁡(F−1LF).

3.1F4step 1.1step 2.1algebra

By [F4], tr⁡(F−1LF)=tr⁡(FF−1L)=tr⁡L=∑i∂xibi(t,χ(t,a)). Also F(0,a)=I, so J(0,a)=1.

4.1givenF6step 2.1step 3.1

The determinant J(t,a) is continuous and nonzero on I by step 2.1. If it were negative for some t, [F6] on the closed segment with endpoints 0 and t would give a point where J=0, contradicting invertibility. Hence J>0.

4.2givenF5step 3.1algebra

If div⁡b=0, step 3.1 gives ∂tJ=0. For each fixed a, the case t=0 is already J(0,a)=1. For t≠0, apply [F5] to J(⋅,a) on [min⁡(0,t),max⁡(0,t)], whose endpoints are distinct; since its derivative is zero, its endpoint difference is zero, so J(t,a)=J(0,a)=1.

5.1givenstep 1.1step 2.1step 3.1step 4.1step 4.2∎

For the torus clause, fix a∈Tn and t0∈I. By continuity, choose a time neighborhood of t0 on which a and χt(a) stay in translation charts. Their transition maps are translations, so the coordinate matrices of Daχt, their determinants, and the coordinate divergence of b agree on overlaps. In these charts the Euclidean calculations in steps 1.1–3.1 prove both differential identities, and χ0=id⁡ gives J(0,a)=1. Since t0 was arbitrary, the identities hold throughout I; applying steps 4.1 and 4.2 on the closed segment between 0 and any t∈I gives positivity and the divergence-free conclusion.

Remarks

The hypotheses describe an already supplied family of diffeomorphisms and its velocity field. This lemma does not prove that an ODE generates such a family, that a trajectory is global, or that any physical fluid or Hamiltonian model satisfies the hypotheses.

TheoremStatement: AI-adaptedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Fine entropy under a specified incompressible periodic flow

Statement

Use log⁡ for the natural logarithm of The natural logarithm as the inverse of the exponential function. Let Tn:=Rn/Λ be a flat translation torus, where Λ=Zℓ1+⋯+Zℓn for a linearly independent basis ℓ1,…,ℓn, with the smooth structure and Euclidean metric descended from Rn. Let μ=c da be its fixed translation-invariant volume measure in dimensionless translation coordinates, with c>0. Let I⊆R be an open interval containing 0. Suppose χ:I×Tn→Tn is jointly C2, each χt is a diffeomorphism, and χ0=id⁡. Let π:TTn→Tn be the tangent-bundle projection and let b:I×Tn→TTn have C1 chart coefficients with π(b(t,a))=a and ∂tχ(t,a)=b(t,χ(t,a)). In translation charts write b=∑ibi∂xi and define div⁡b=∑i∂xibi. Assume div⁡b=0 in these charts. Let f0∈C1(Tn) satisfy 0<m≤f0≤M<∞ and ∫f0 dμ=1, and define ft:=f0∘χt−1. Then ft is a probability density for every t∈I, the integral H(ft):=−∫Tnftlog⁡ft dμ is finite, and H(ft)=H(f0). Here H is the dimensionless fine entropy of this density; multiplication by a physical constant such as kB is a separate interpretation.

Facts & Assumptions

Given: The flat torus, fixed measure μ=c da, jointly C2 diffeomorphisms χt, divergence-free velocity, and positive C1 probability density f0 bounded above and below as in the Statement.

[F1]

Jacobian identity for a specified smooth diffeomorphic flow: on the flat translation torus, for a jointly C2 diffeomorphic family χt with χ0=id⁡ and a C1 time-dependent vector field b satisfying ∂tχ(t,a)=b(t,χ(t,a)), the local-chart Jacobian satisfies ∂tJ=Jdiv⁡b, is positive, and equals 1 when the chart divergence ∑i∂xibi vanishes.

[F2]

Change of variables for an injective C1 map on a compact Jordan set: an injective C1 map with invertible derivative on an open neighborhood of a compact Jordan set obeys the Riemann change-of-variables formula for bounded integrable functions.

[F3]

The natural logarithm as the inverse of the exponential function defines log⁡ on (0,∞) as the inverse of the real exponential; Order, continuity, range, and the product, quotient, and reciprocal laws for the natural logarithm makes it continuous there. In particular it is bounded and continuous on every compact interval [m,M]⊂(0,∞).

Proof

1.1givenF1

Apply [F1] to the given family. Its hypotheses are exactly the supplied C2 diffeomorphisms, identity initial map and C1 divergence-free velocity, so every local translation-chart Jacobian is Jt=1.

2.1givenF1F2step 1.1

For fixed t, around each point choose a source chart and a smaller coordinate neighborhood whose closure maps under χt into a target chart. Continuity and compactness give finitely many closed coordinate boxes Qj with interiors covering the torus and with each χt(Qj) inside a target chart. Put E1=Q1 and Ej=Qj∖⋃i<jint⁡(Qi); each Ej is compact, the pieces cover the torus, and any overlap lies on one of the box faces. Each Ej is Jordan because its boundary lies in a finite union of coordinate faces, which are compact C1 hypersurface patches. A C1 graph patch has bounded slope on a compact subpatch, so its mesh-h box cover has total n-volume O(h) and hence content zero. On compact chart neighborhoods χt has bounded derivative, so it maps these zero-content boundaries to zero-content boundaries. Each compact piece and its image therefore have Jordan boundaries and lie in source and target charts; [F2] applies there. Summing over the pieces, overlaps do not affect the integrals, translation chart changes preserve da, and the constant c is the same on both sides. Thus for every continuous φ on the torus, ∫φ∘χt dμ=∫φJt dμ=∫φ dμ.

3.1givenstep 2.1

The function ft=f0∘χt−1 is positive and continuous. The substitution of step 2.1 with φ=ft gives ∫ft dμ=∫f0 dμ=1, so ft is normalized.

4.1givenF3step 2.1algebra∎

The bounds m≤f0≤M also hold for ft; by [F3], log⁡ is continuous on [m,M], so ftlog⁡ft is bounded and continuous on the compact torus and its integral is finite. Applying step 2.1 to φ=ftlog⁡ft gives ∫ftlog⁡ft dμ=∫(ftlog⁡ft)∘χt dμ=∫f0log⁡f0 dμ, since ft∘χt=f0. Negating proves H(ft)=H(f0).

Remarks

This is a statement about supplied smooth mathematical flow data. It does not assert existence of a flow for a vector field, Hamiltonian structure, relaxation, mixing, ergodicity, monotonicity of coarse entropy, or a physical entropy interpretation.

TheoremStatement: AI-adaptedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Gaussian covariance for positive quadratic forms and a boundary counterexample

Statement

Assume the Axiom of Countable Choice ACω as defined in The Axiom of Countable Choice (ACω). Let n≥1, let K∈Rn×n be symmetric positive definite, and let β>0. For x∈Rn put H(x)=12xTKx, and use the Euclidean Lebesgue measure λn to define the density Z−1e−βH(x). Then

Z=∫Rne−βH(x)dx=(2πβ)n/2(det⁡K)−1/2,E[x]=0,E[xxT]=β−1K−1,E[H]=n2β.

Positive definiteness and β>0 are essential for this normalized full-space Gaussian. The unrestricted identity E[x H′(x)]=β−1 is false without a boundary hypothesis: on [−L,L] with L>0 and H≡0, the normalized expectation on the left is 0.

Facts & Assumptions

Given: ACω, n≥1, symmetric positive-definite K, β>0, and H(x)=12xTKx on Rn.

[F1]

Gaussian, gamma and beta product integrals: for a>0, ∫Re−ax2dx=π/a, odd Gaussian moments vanish, and the second moment is π/(2a3/2).

[F2]

Under ACω, a C1 diffeomorphism transforms nonnegative Lebesgue integrals by its absolute Jacobian (A C^1 diffeomorphism satisfies the change-of-variables formula for nonnegative Lebesgue measurable functions).

[F3]

Under ACω, Euclidean Lebesgue measure is the finite product of one-dimensional Lebesgue measures on Borel rectangles, and each box has measure equal to the product of its side lengths (On Borel subsets of R^{m+n}, the product lambda_m times lambda_n agrees with lambda_{m+n}, A box in Rn with parameters ai≤bi is Lebesgue measurable of measure ∏i<n(bi−ai), whichever of its faces are included); each factor is sigma-finite. Tonelli factors nonnegative products, and Fubini factors absolutely integrable signed products (Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Fubini's theorem for L^1 functions on a sigma-finite product).

[F4]

The Gaussian integrals and moments in [F1] are improper Riemann integrals. On bounded intervals they agree with Lebesgue integrals under ACω, and monotone convergence passes their nonnegative truncations to R (Riemann–Lebesgue comparison for distribution test integrands, Monotone convergence for the integral).

Proof

1.1givenconstructalgebra

We construct R with K=RTR by induction on n. For n=1, take R=K11>0. For n>1, write K=(abTbC); positive definiteness gives a>0 and S:=C−bbT/a positive definite, since yTSy=(−bTy/a,y)TK(−bTy/a,y)>0 for y≠0. By induction S=R1TR1; then R=(abT/a0R1) satisfies RTR=K. Its diagonal entries are positive, so det⁡R>0; taking determinants gives det⁡K=(det⁡R)2 and therefore det⁡R=det⁡K.

2.1givenF1F2F3F4step 1.1algebra

Put T(x)=Rx. This is a C1 diffeomorphism of Rn with ∣det⁡DT∣=det⁡R. Since H(x)=∣Rx∣2/2, [F2] gives Z=(det⁡R)−1∫Rne−β∣y∣2/2 dλn(y). By [F3], λn is the finite product of n sigma-finite one-dimensional Lebesgue measures; Tonelli and the factorization e−β∣y∣2/2=∏ie−βyi2/2 give the integral as I0(β/2)n. By [F1] and [F4], I0(β/2)=2π/β. Hence Z=(2π/β)n/2(det⁡K)−1/2, which is finite and positive.

3.1F1F3step 2.1algebra

Let Y have the normalized density proportional to e−β∣y∣2/2 with respect to λn. For each i, [F3] and Tonelli factor the nonnegative second-moment integral into I1(β/2)I0(β/2)n−1, while the normalizer is I0(β/2)n; [F1] gives E[Yi2]=I1(β/2)/I0(β/2)=1/β. Also ∫∣yi∣e−β∣y∣2/2dλn<∞, since ∣yi∣≤1+yi2 and the same Tonelli factorization makes both integrals finite. For i≠j, ∣yiyj∣≤(yi2+yj2)/2, so the cross moments are absolutely integrable. Fubini [F3] factors these signed moments; the one-dimensional first and odd moments vanish by [F1]. Thus E[Y]=0 and E[YYT]=β−1I.

4.1F2step 1.1step 3.1algebra

By [F2], the normalized law of X is the pullback of the law of Y under x=R−1y. The absolute integrability established in step 3.1 permits the same change of variables for signed first and second moments by applying [F2] to their positive and negative parts. Hence E[X]=0 and E[XXT]=R−1(β−1I)R−T=β−1(RTR)−1=β−1K−1. Since H(X)=∣Y∣2/2, summing the diagonal second moments gives E[H]=n/(2β).

5.1givenF3algebra∎

For the stated boundary counterexample, the density on [−L,L] is the constant (2L)−1 because H=0 and the normalizing integral is 2L by [F3]. Since H′(x)=0 everywhere, E[xH′(x)]=0≠β−1 for β>0.

Remarks

If K has a nonzero kernel, the integrand is constant along an unbounded kernel direction and Z=∞; if β≤0, the full-space integral also diverges. The theorem does not assert a named-degree-of-freedom equipartition rule on arbitrary domains or for discrete variables.

TheoremStatement: AI-adaptedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Exact spherical-shell coordinate marginal and fixed-block Gaussian limit

Statement

Assume the Axiom of Countable Choice ACω as defined in The Axiom of Countable Choice (ACω). Let D≥2, 1≤k<D, and R>0. Let Z1,…,ZD be independent real random variables with common Lebesgue probability density ϕ(x):=(2π)−1/2e−x2/2; its integral is one by the Gaussian integral in Gaussian, gamma and beta product integrals, and this is the standard-normal coordinate law. Set Z=(Z1,…,ZD) and define P:=RZ/∣Z∣ when Z≠0, and P:=0 when Z=0. The joint density of Z makes P(Z=0)=0. The distribution of P is rotationally invariant, meaning that for every orthogonal map Q∈O(D) and Borel set A⊆RD, P(QP∈A)=P(P∈A). It is concentrated on the sphere in the precise sense that P(∣P∣=R)=1. Its first k coordinates have Lebesgue density

fD,k(y)=Γ(D/2)Γ((D−k)/2)(πR2)k/2(1−∣y∣2R2)(D−k−2)/2(∣y∣<R),

and zero for ∣y∣≥R. This density is integrable even when it diverges at the boundary. Also E[Pi]=0, E[PiPj]=0 for i≠j, and E[Pi2]=R2/D.

For fixed k and σ>0, set RD=D σ and let D→∞ through integers D>k. If g:Rk→R is bounded and Lipschitz in the Euclidean norm, meaning that finite constants Bg,Lg satisfy ∣g(x)∣≤Bg and ∣g(x)−g(y)∣≤Lg∣x−y∣ for all x,y∈Rk, then

E g((P1,…,Pk))⟶E g(σ(Z1,…,Zk)).

This is a fixed-block bounded-Lipschitz test limit; it makes no claim uniform in k or for unbounded test functions. If k=D, the projection is singular on the sphere and the displayed Lebesgue density is not asserted.

Facts & Assumptions

Given: The Axiom of Countable Choice ACω, integers D,k with D≥2 and 1≤k<D, radius R>0, and the product probability law of D independent coordinates with density ϕ(x)=(2π)−1/2e−x2/2.

[F1]

Gaussian, gamma and beta product integrals: Γ(a) is finite for a>0, Γ(1/2)=π, positive gamma densities with the same scale convolve by adding their shape parameters, and the Gaussian second and fourth moments are 1 and 3 for a standard normal coordinate.

[F2]

A C^1 diffeomorphism satisfies the change-of-variables formula for nonnegative Lebesgue measurable functions: assuming ACω, for a C1 diffeomorphism T:U→V between open subsets of Euclidean space and nonnegative Lebesgue-measurable f on V, the Lebesgue integral transforms by ∣det⁡DT∣.

[F3]

On Borel subsets of R^{m+n}, the product lambda_m times lambda_n agrees with lambda_{m+n}: assuming ACω, product Lebesgue measure and Euclidean Lebesgue measure agree on Borel subsets of Rm+n for positive integers m,n.

[F4]

Tonelli's theorem for nonnegative measurable functions on a sigma-finite product: for a nonnegative product-measurable function on a product of sigma-finite measure spaces, either order of integration agrees with the product integral.

[F5]

The Axiom of Countable Choice (ACω): ACω says every sequence of nonempty sets has a choice function; it is the explicit set-theoretic hypothesis used by [F2].

Proof

1.1givenF2F3F5algebra

Under Countable Choice [F5], the Borel product-measure identity [F3] and independence give the joint Lebesgue density ϕD(z)=(2π)−D/2e−∣z∣2/2. A singleton has Lebesgue measure zero, so Z≠0 almost surely. For any orthogonal Q and Borel A⊆RD, the map Q−1:RD→RD is a C1 diffeomorphism with ∣det⁡Q−1∣=1; applying [F2] to 1AϕD and using ∣Q−1z∣=∣z∣ gives P(QZ∈A)=P(Z∈A). Since Q(RZ/∣Z∣)=RQZ/∣QZ∣ off the null event, the law of P is invariant under every orthogonal map, and ∣P∣=R almost surely.

1.2F1F2algebra

Write Z=(A,B) with A∈Rk and B∈RD−k. The vectors are independent. On each open half-line (0,∞) and (−∞,0), the map z↦z2/2 is a C1 diffeomorphism onto (0,∞); applying [F2] to the Gaussian density times the indicator of any Borel target set on each branch shows that each Zj2/2 has density w−1/2e−w/Γ(1/2). Repeated gamma convolution from [F1] then shows T:=∣B∣2 has density tb−1e−t/2/(2bΓ(b)) for t>0, where b=(D−k)/2, independently of A.

1.3givenconstructalgebra

For fixed t>0, the first k coordinates y=RA/∣A∣2+t range over ∣y∣<R. The inverse is A=t y/R2−∣y∣2. Its derivative has radial eigenvalue t R2/(R2−∣y∣2)3/2 and k−1 tangential eigenvalues t/(R2−∣y∣2)1/2, hence Jacobian tk/2R2/(R2−∣y∣2)k/2+1.

2.1F1F2F3F4step 1.2step 1.3algebra

The map Ψ(y,t):=(A(y,t),t) from U:={(y,t):∣y∣<R, t>0} to V:=Rk×(0,∞) is a C1 diffeomorphism, with inverse (a,t)↦(Ra/∣a∣2+t,t). Apply [F2] to the nonnegative joint density of (A,T) and this diffeomorphism. The product Lebesgue measure on Rk×(0,∞) agrees with Euclidean Lebesgue measure on Borel sets by [F3], and Tonelli's theorem [F4] gives the marginal density by integration in t. Substituting the Gaussian density of A, the gamma density of T, and the Jacobian from step 1.3 yields fD,k(y)=R2(2π)k/22bΓ(b)(R2−∣y∣2)k/2+1∫0∞tD/2−1exp⁡ ⁣(−tR22(R2−∣y∣2))dt. The substitution u=tR2/[2(R2−∣y∣2)] evaluates the integral as Γ(D/2)[2(R2−∣y∣2)/R2]D/2. Since D/2=k/2+b, simplification gives the displayed density. Its boundary singularity, when present, has exponent b−1>−1 and is integrable in the radial variable; the pushforward construction shows total mass one. Because T=∣B∣2>0 almost surely, the marginal has no boundary atom at ∣y∣=R.

2.2step 1.1givenalgebra

Sign changes of one Gaussian coordinate show E[Pi]=0 and E[PiPj]=0 for i≠j. Permutation invariance makes all E[Pi2] equal, while ∑i=1DPi2=R2 almost surely; summing expectations gives E[Pi2]=R2/D.

3.1F1step 1.1algebra∎

For each finite D, take the product Gaussian vector Z=(Z1,…,ZD) in that dimension and set SD:=D−1∑i=1DZi2. By [F1] and independence, ESD=1, E[(SD−1)2]=2/D, and E∣Z1:k∣2=k. Fix 0<δ<1/2 and M>0. On ED={∣Z1:k∣≤M, ∣SD−1∣≤δ}, SD∈[1/2,3/2] and the identity ∣u−1/2−1∣=∣u−1∣/[u(1+u)] gives ∣SD−1/2−1∣≤2δ; hence ∣σZ1:k/SD−σZ1:k∣≤2σMδ. If g is bounded Lipschitz with bound Bg and Lipschitz constant Lg, then M21{∣Z1:k∣>M}≤∣Z1:k∣2 and δ21{∣SD−1∣>δ}≤(SD−1)2 imply P(EDc)≤k/M2+2/(Dδ2). Therefore ∣Eg(P1:k)−Eg(σZ1:k)∣≤Lg2σMδ+2Bg(kM2+2Dδ2). Choose M large, then δ small, then D large. The right side becomes arbitrarily small, proving the stated test convergence without assuming independence between Z1:k and SD.

Remarks

This is a specified rotational Gaussian-direction preparation. It is not a general coarea formula, a proof that every energy-shell measure has this marginal, or a dynamical equilibration result.

PostulateStatement: Literature-sourcedProof: Not applicablejudge pass (gpt-6-luna)Open item page →

Declared example preparations and thermodynamic interpretation

Physical scope

Separate finite classical canonical, grand-canonical, spherical-shell and product preparations in the named model examples; no claim of equilibration, empirical truth, or statistical inference from observations.

Postulate

Here log⁡ denotes the natural logarithm defined in The natural logarithm as the inverse of the exponential function.

For each named finite classical model in this page, adopt the preparation explicitly stated in its theorem. In a canonical branch, take a measurable phase space with a specified reference measure ν, a measurable Hamiltonian H, an inverse energy parameter β>0, and a finite positive partition function Z=∫e−βH dν; the adopted probability measure is dP=Z−1e−βH dν. In a grand-canonical branch, take disjoint particle-number sectors with sector partition functions ZN and a dimensionless activity z>0 for which Ξ=∑N≥0zNZN is finite; sector N has probability weight zNZN/Ξ. The spherical-shell and product preparations are separate explicitly specified probability laws in their consuming theorem statements. These are alternative model assumptions, not consequences of a Hamiltonian flow and not assertions that a laboratory system realizes them.

For a canonical model whose reference measure and Hamiltonian are fixed, adopt the statistical thermodynamic definitions U=EP[H], F=−β−1log⁡Z, and S=kB(log⁡Z+βU) when the displayed integrals are finite. For a volume-dependent model, define the model pressure as −∂VF only along the declared fixed-shape volume family. Here kB>0 has SI units J K−1, temperature T>0 has units K, and β=(kBT)−1 has units J−1. The fixed reference measure is part of the preparation; changing it changes the numerical partition function and entropy convention.

For the continuous particle examples, adopt the operational species and counting convention of Classical particle counting and physical observables. Use one fixed primitive action reference a0>0 in J s, and write h0:=a0 where the consuming model uses the notation h0. In a Cartesian physical phase chart dq dp has action dimension (J s)dN, so the reference dq dp/(∏sNs! a0dN) is dimensionless. Each model retains its declared operational labels: an omitted factorial for resolved harmonic coordinates is not replaced by a particle-exchange quotient. The finite spin examples use their stated discrete counting measure, without a continuous phase-cell factor. Choosing a numerical reference action, or identifying it with Planck's constant, is an additional convention; no such value or identification is imposed here. The Gaussian and gamma identities in Gaussian, gamma and beta product integrals establish the mathematical normalizations, not the physical preparation or empirical equilibration.

Remarks

A finite partition function and normalized density are mathematical conditions of each canonical consumer. No negative-temperature unbounded gas, ergodic hypothesis, observed outcome, or universal ensemble-selection principle is included.

Physics theoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Confined coupled quadratic canonical systems

Physical scope

Finite classical continuous canonical model with positive diagonal mass matrix, strictly positive stiffness matrix, fixed positive action reference and beta>0; the Lebesgue product calculations assume ACω; no unpinned zero mode or empirical equilibration claim.

Statement

Assume ACω as defined in The Axiom of Countable Choice (ACω) for the Euclidean Lebesgue product measure used below. Let n≥1, let M=diag⁡(m1,…,mn) with mi>0 in kg, and let K be a real symmetric positive-definite stiffness matrix with entries in J m−2. For displacement q∈Rn in metres and momentum p∈Rn in kg m s−1, define H(q,p)=12pTM−1p+12qTKq. Adopt the canonical probability density proportional to e−βH relative to dq dp/h0n, where β>0 has units J−1 and the fixed action reference h0>0 has units J s. Under the preparation of Declared example preparations and thermodynamic interpretation, the dimensionless partition function and moments are Z=(2πβh0)ndet⁡Mdet⁡K,E[q]=E[p]=0,E[qqT]=β−1K−1,E[ppT]=β−1M,E[qpT]=0. The model mean energy is U=E[H]=n/β. With kB>0 in J K−1 and T>0 in K from Declared example preparations and thermodynamic interpretation, so that β=(kBT)−1, define the fixed-parameter heat-capacity response by Charm:=(∂U∂T)M,K,h0=nkB. It has units J K−1. This model has no volume variable, so Charm is not a fixed-volume derivative.

Facts & Assumptions

Given: ACω, n≥1, M positive diagonal, K symmetric positive definite, β>0, fixed h0>0, and the displayed continuous-coordinate canonical model.

[F1]

Declared example preparations and thermodynamic interpretation: in each named model with finite positive Z, the canonical probability is Z−1e−βH relative to the specified reference; U=EH and the fixed-reference thermodynamic interpretation is adopted separately from the probability preparation.

[F2]

Gaussian covariance for positive quadratic forms and a boundary counterexample: assuming ACω, for a symmetric positive-definite quadratic form on Rn, the normalized Gaussian integral is (2π/β)n/2(det⁡K)−1/2, with mean zero and covariance β−1K−1.

[F3]

Under ACω, Euclidean phase-space Lebesgue measure is the product of the q and p Lebesgue measures; Tonelli factors nonnegative densities and Fubini factors their integrable moments (On Borel subsets of R^{m+n}, the product lambda_m times lambda_n agrees with lambda_{m+n}, Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Fubini's theorem for L^1 functions on a sigma-finite product).

[A1]

The canonical density is a model preparation, not a claim that an isolated trajectory equilibrates or that a laboratory realizes this law.

Proof

1.1givenF1F2F3algebra

By [F1], the phase-space weight factors as e−βpTM−1p/2e−βqTKq/2. The product-measure identity and Tonelli [F3] justify multiplying the two nonnegative Lebesgue integrals. Applying [F2] to the q factor gives normalization (2π/β)n/2(det⁡K)−1/2 and covariance β−1K−1.

1.2givenF2algebra

The momentum quadratic form has symmetric positive-definite matrix M−1. Applying [F2] with K=M−1 gives its integral (2π/β)n/2det⁡M and covariance of p equal to β−1(M−1)−1=β−1M.

2.1A1F3step 1.1step 1.2algebra

Multiplying the two integrals from steps 1.1–1.2 and dividing by the reference factor h0n gives Z=(2π/(βh0))ndet⁡M/det⁡K. The phase-space Lebesgue measure is the product of its q and p factors by [F3]; the normalized canonical density therefore factors into independent q and p laws. Each factor is even, so its mean is zero, and its covariance is the one computed in steps 1.1–1.2.

3.1A1F1step 1.1step 1.2algebra∎

By [F1], taking the expectation of H and summing the diagonal traces of the two covariance identities gives U=12E[pTM−1p]+12E[qTKq]=n2β+n2β=nβ. The adopted relation β−1=kBT makes U=nkBT. Holding the model parameters M,K,h0 fixed therefore gives Charm=(∂U/∂T)M,K,h0=nkB, with units J K−1.

Remarks

Positive definiteness is a confinement condition. If K has a translation zero mode, the full-space configurational integral diverges and this normalized ensemble does not exist without an additional constraint or confinement.

Physics theoremStatement: Literature-sourcedProof: AI-adaptedprecheck passOpen item page →

Exact canonical and grand ideal-gas preparations

Physical scope

Finite-N canonical sectors of a classical ideal gas of one operationally indistinguishable species on a bounded rectangular box, their grand-canonical sum including N=0, and a separately adopted three-dimensional classical H=c|p| branch; the Lebesgue product and polar calculations assume ACω; excludes quantum statistics, actual equilibration, and identifying the massless branch with m->0.

Statement

Assume ACω as defined in The Axiom of Countable Choice (ACω), for the Euclidean Lebesgue product and polar-coordinate suppliers used below. Let exp⁡(x):=∑n=0∞xn/n! be the real exponential of The real exponential function and the number e by a power series, and let log⁡ be the natural logarithm of The natural logarithm as the inverse of the exponential function, with log⁡′(x)=1/x for x>0 by The natural logarithm has derivative 1/x and equals the integral from 1 to x of 1/t. Let d,N≥1 and let B=∏j=1d[0,Lj]⊂Rd, where each side length Lj>0 is in metres and V:=∏jLj is in md. Let m>0 be the common particle mass in kg, β>0 in J−1, and let h0:=a0>0 be the fixed action reference in J s from Declared example preparations and thermodynamic interpretation. Take one species of N operationally indistinguishable particles as in Classical particle counting and physical observables. On its auxiliary labelled cover, positions qi∈B are in metres and momenta pi∈Rd are in kg m s−1; the Hamiltonian is permutation invariant and the reference measure is dνN=dq dpN!h0dN. The symmetric labelled law represents the quotient preparation; a one-particle marginal below means the law of a uniformly selected particle. For N≥1, define HN(q,p)=∑i=1N∣pi∣22m and use the canonical preparation of Declared example preparations and thermodynamic interpretation. For the vacuum sector N=0, define B0×Rd⋅0 to be a singleton with reference measure ν0=1, set H0=0, and use empty-product conventions V0=0!=h00=1; thus Z0=1. The mass parameter m applies when particles are present; the vacuum sector has no mass coordinate. For all N≥0, set λ:=h0β/(2πm) (metres) and ZN=∫e−βHN dνN=VNN!λdN. For N≥1, the labelled-cover position coordinates are independent uniform points of B and are independent of the momenta; all momentum components are independent centered Gaussians of variance m/β. The joint law is permutation invariant. The energy satisfies U=EHN=dN/(2β) and Var⁡(HN)=dN/(2β2), with variance in J2. With the adopted Boltzmann constant kB>0 in J K−1 and temperature T:=(kBβ)−1>0 in K from Declared example preparations and thermodynamic interpretation, define F=−β−1log⁡ZN and scale all box side lengths uniformly when varying V, holding the shape, N,T,d,m,h0 fixed. The model pressure P:=−(∂F/∂V)N,T,d,m,h0, shape has units J m−d (equivalently N m1−d; Pa when d=3) and obeys PV=N/β=NkBT. The fixed-reference entropy is S=kB(log⁡ZN+βU) and retains the chosen N! and h0 convention.

For d=3, the speed v=∣p∣/m of a uniformly selected particle, in m s−1, has density f(v)=4π(mβ2π)3/2v2e−βmv2/2,v≥0. Its squared-speed expectation, modal speed, and mean speed are E[v2]=3/(mβ) in m2 s−2, 2/(mβ) in m s−1, and E[v]=8/(πmβ) in m s−1, respectively.

For the separate grand-canonical branch of Declared example preparations and thermodynamic interpretation, let μchem be a chemical potential in joules and z=exp⁡(βμchem)>0 the dimensionless activity. Specify its full preparation on Γg:=⨆j≥0({j}×Γj), where Γj=Bj×Rdj for j≥1 and Γ0 is the vacuum singleton. Use the sigma-algebra of sets A whose slices Aj:={x∈Γj:(j,x)∈A} are Borel, and adopt, whenever 0<Ξ:=∑j≥0zjZj<∞, the probability law Pg(A):=Ξ−1∑j≥0zj∫Aje−βHj dνj. The particle count is the map N:Γg→{0,1,…}, N(j,x)=j. This full-law assumption supplements the sector weights; those weights alone do not determine the within-sector law. Then Ξ=∑N=0∞zNZN=exp⁡(zV/λd),P(N=j)=exp⁡(−ν)νjj!,ν:=zV/λd, so EN=Var⁡(N)=ν and the conditional law given N is the canonical law above.

A distinct three-dimensional classical ultra-relativistic model adopts Hur(p)=c∣p∣ for every p∈R3, with c>0 in m s−1 and p in kg m s−1, in the same box volume V in m3. This extends the nonzero massless momentum formula continuously to p=0 only to define a density on all of R3; the point has Lebesgue measure zero and no physical zero-momentum massless state is asserted. Its separate canonical partition function is ZN,ur=VNN!h03N(8π(βc)3)N, with Uur=3N/β J and Var⁡(Hur)=3N/β2 J2. Define its fixed-volume heat capacity by CV,ur:=(∂Uur/∂T)V,N,c,h0=3NkB; it has units J K−1. Its three-dimensional pressure obeys PV=N/β and has units Pa. This branch is not obtained by setting m=0 in the Gaussian model and makes no photon-number or quantum-statistics claim.

Facts & Assumptions

Given: ACω, d,N≥1, the rectangular box B with side lengths in metres and volume V in md, one operationally indistinguishable species, m>0 kg, h0=a0>0 J s, β>0 J−1, and the Hamiltonian and symmetric reference measure in the Statement.

[F1]

Declared example preparations and thermodynamic interpretation: for a specified model with finite positive partition function, the canonical density is Z−1e−βH relative to its fixed reference; grand-sector weights are zNZN/Ξ. This is an adopted preparation, not an equilibration theorem.

[F2]

Gaussian, gamma and beta product integrals: positive Gaussian coefficients have the stated Gaussian moments and finite products; gamma integrals, equal-scale gamma convolution, and moment differentiation follow from its explicit evaluations and Gamma recurrence.

[F3]

Under ACω, product Lebesgue measure agrees with Euclidean Lebesgue measure on Borel rectangles (On Borel subsets of R^{m+n}, the product lambda_m times lambda_n agrees with lambda_{m+n}).

[F5]

The Lebesgue factors are sigma-finite; Tonelli factors nonnegative product integrals and Fubini factors absolutely integrable signed product integrals (Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Fubini's theorem for L^1 functions on a sigma-finite product).

[F6]

Polar coordinates in R3 give ∫g(p) dλ3(p)=∫0∞∫S2g(rω)r2 dσ(ω) dr for nonnegative Borel g, under ACω; σ(S2)=4π by the unit-ball volume formula (Polar coordinates decompose Lebesgue measure into r^{n-1} dr d sigma, The closed form for the volume of the unit n-ball).

[F7]

The Gaussian integrals and moments in [F2] are improper Riemann integrals. Under ACω, their bounded truncations agree with Lebesgue integrals, and monotone convergence passes nonnegative truncations to the full line (Riemann–Lebesgue comparison for distribution test integrands, Monotone convergence for the integral).

[F8]

Massless test-particle action and momentum scale: the adopted massless test sector has nonzero future-null mechanical momentum P and does not identify it with mU at m=0.

[F9]

Mechanical momentum and observer energy: for an observer n with g(n,n)=−1, define spatial momentum qn(P)=P+g(P,n)n and observer energy En(P)=−c g(P,n); in the massless sector P is nonzero future null.

[F10]

Observers, rest spaces and the spatial projection at an event: in a time-oriented Lorentzian background (M,g,N), N is the smooth future timelike field; at an event n=N/−g(N,N) is an observer, whose rest space is n⊥ and whose spatial projection is hn(P)=P+g(n,P)n.

[F11]

The observer rest-space decomposition and the spatial projection: g(hn(P),hn(P))=g(P,P)+g(P,n)2, and g is positive definite on the three-dimensional rest space n⊥.

[F12]

Classical particle counting and physical observables: for one operationally indistinguishable species, the symmetric labelled lift uses the factor 1/N!; label-specific one-particle marginals of a permutation-invariant state are interpreted as those of a uniformly selected particle.

[F13]

The real exponential function and the number e by a power series: exp⁡(x) is the everywhere-convergent series ∑n≥0xn/n!.

[F14]

The exponential is positive and satisfies exp⁡(−x)=1/exp⁡(x): exp⁡(−x)=1/exp⁡(x) and exp⁡(x)>0 for real x.

[F15]

The natural logarithm as the inverse of the exponential function: log⁡ is the real natural logarithm on (0,∞).

[F17]

Order, continuity, range, and the product, quotient, and reciprocal laws for the natural logarithm gives log⁡(xy)=log⁡x+log⁡y and log⁡(x/y)=log⁡x−log⁡y for positive x,y; finite iteration gives log⁡(xk)=klog⁡x for integers k≥0.

[F18]

Under ACω, every point and coordinate hyperplane in R3 is Lebesgue null (A box with a degenerate side is Lebesgue null, and so is every coordinate hyperplane in Rn).

[A1]

The canonical probability laws, grand-sector weights and thermodynamic identifications are adopted as in [F1], with the full grand law specified in [A3]; the calculations establish conditional consequences only. The one-species operational indistinguishability and its permutation-invariant labelled lift are model assumptions as in [F12], not consequences of common mass.

[A2]

The ultra-relativistic branch separately adopts the classical Hamiltonian Hur(p)=c∣p∣ in one observer's orthonormal spatial momentum coordinates. It is not a massless limit of the massive Hamiltonian and is not a quantum radiation model.

[A3]

In the grand branch adopt the full law Pg on the disjoint union and slice sigma-algebra specified in the Statement. In particular, for a Borel E⊆Γj, Pg({j}×E)=Ξ−1zj∫Ee−βHj dνj. This is an additional preparation assumption, not a consequence of sector weights alone.

Proof

Given: ACω, d,N≥1, the rectangular box B with side lengths in metres and volume V in md, one operationally indistinguishable species, m>0 kg, h0=a0>0 J s, β>0 J−1, and the Hamiltonian and symmetric reference measure in the Statement.

1.1givenF1F2F3F4F5F7F12algebra

By [F3], Euclidean Lebesgue measure in the position and momentum coordinates is their product measure; [F5] then factors the nonnegative normalization integral. The box-volume formula [F4] gives λdN(BN)=VN. Each of the dN scalar momentum factors contributes ∫Re−βp2/(2m)dp=2πm/β by [F2] and [F7]. Including the stipulated factor 1/(N!h0dN) from [F12] yields ZN=VN(N!h0dN)−1(2πm/β)dN/2=VN/(N!λdN).

1.2givenA2F8F9F10F11algebra

For the separate massless branch, let Ntime denote the background's smooth future-timelike field N from [F10] and normalize it at the event by n=Ntime/−g(Ntime,Ntime). [F8] makes P future null, so g(P,n)<0. By [F11], ∣qn(P)∣2=g(qn(P),qn(P))=g(P,n)2, hence En(P)=−cg(P,n)=c∣qn(P)∣. Identifying qn with its components in an orthonormal basis of the positive rest space gives the adopted Hur(p)=c∣p∣ with the stated units.

2.1givenF2F3F5step 1.1

After normalization, the canonical density is the product of the uniform position density on BN and the dN scalar Gaussian momentum densities. For measurable coordinate rectangles, [F3] and [F5] show their probabilities factor; hence the positions and all momentum components are independent, each position is uniform on B, and each momentum component is centered with variance m/β by [F2].

2.2A1A3F1F13F14step 1.1algebra

In the grand preparation [A3] with sector weights from [F1], use the vacuum convention Z0=1 and substitute ZN=VN/(N!λdN) for all N≥0 into the sector sum. With ν=zV/λd, [F13] gives Ξ=∑N≥0νN/N!=exp⁡(ν), and [F14] gives the normalized sector weight νN/(N!exp⁡(ν))=exp⁡(−ν)νN/N!. Each summand in [A3] is a nonnegative measure on its sector, so countable additivity and this finite positive total mass make Pg a probability measure.

2.3A1A2F1F2F3F5F6F18step 1.1algebra

By [F6], the momentum integral for one ultra-relativistic particle is ∫R3e−βc∣p∣ dλ3(p)=4π∫0∞r2e−βcrdr=4πΓ(3)(βc)3=8π(βc)3, using [F2]. The point p=0 is Lebesgue null by [F18], so the continuous extension adds no massless zero-momentum state. Product measure and Tonelli [F3, F5] give the N-particle product, and the stated positional/counting reference then gives ZN,ur.

3.1givenF2step 2.1algebra

Each scalar momentum coordinate has second moment m/β and fourth moment 3m2/β2 from [F2]. Thus each term piα2/(2m) has mean 1/(2β) and variance 1/(2β2); independence from step 2.1 makes means and variances add, giving the stated U and Var⁡(HN).

3.2givenF2F6step 2.1algebra

For d=3, write a=β/(2m). The one-particle momentum density is (a/π)3/2e−a∣p∣2. By [F6], for every Borel A⊆[0,∞) the probability that v=∣p∣/m lies in A is 4π(a/π)3/2∫mAr2e−ar2 dr. The substitution r=mv therefore gives the speed density 4π(mβ/(2π))3/2v2e−βmv2/2 on v≥0.

3.3A3F13F14step 1.1step 2.2algebra

Shifting the index in the absolutely convergent series [F13] gives ∑j≥0jνj/j!=νexp⁡(ν) and ∑j≥0j(j−1)νj/j!=ν2exp⁡(ν). After normalization by [F14], EN=ν and E[N(N−1)]=ν2, hence Var⁡(N)=E[N(N−1)]+EN−(EN)2=ν. For each j≥0, step 2.2 gives Pg(N=j)=zjZj/Ξ>0. For any Borel E⊆Γj, [A3] therefore gives Pg({j}×E∣N=j)=(Ξ−1zj∫Ee−βHj dνj)/(Ξ−1zjZj)=Zj−1∫Ee−βHj dνj, exactly the canonical law, including the vacuum singleton.

4.1A1F1F15F16F17step 1.1step 3.1algebra

By [F17] and [A1], F=−β−1[Nlog⁡(V/λd)−log⁡N!] for fixed T,N and fixed box shape. The log derivative [F16] gives −∂VF=N/(βV), with P in J m−d, and therefore PV=N/β. The entropy definition in [F1] gives S=kB(log⁡ZN+βU) with the same fixed h0 and factorial convention.

4.2F2step 3.2algebra

The speed density is proportional to v2e−av2 with a=βm/2. Its derivative has the sign of 2v(1−av2) for v>0, so its unique positive maximum occurs at vmode=a−1/2. The second moment is 3/(2a)=3/(mβ) by the gamma recurrence in [F2]. The mean is Ev=4a3/2π∫0∞v3e−av2dv=4a3/2πΓ(2)2a2=2πa=8πβm. These are model-distribution moments, not values for an individual measured particle.

5.1A1F1F2F5F16step 2.3algebra∎

Under the normalized radial density proportional to r2e−βcr, gamma moments give mean energy 3/β and variance 3/β2 per particle. Independence from the product preparation adds these to Uur=3N/β and Var⁡(Hur)=3N/β2; differentiating at fixed V,N,c,h0 with β−1=kBT gives CV,ur=3NkB J K−1. Since ZN,ur is proportional to VN, [F16] gives −∂VF=N/(βV) and PV=N/β, with P in Pa for this three-dimensional branch.

Remarks

The source M18 Ising-chain supplier is not used here. The massless momentum relation is a consequence of the adopted null-sector convention and the stated observer projection; the classical canonical H=c∣p∣ preparation is a separate model assumption. No empirical data or precision are asserted.

Physics theoremStatement: Literature-sourcedProof: AI-adaptedprecheck passjudge pass (gpt-6-luna)Open item page →

Finite periodic Ising-ring correlations and exact decay

Physical scope

Finite periodic one-dimensional nearest-neighbor Ising ring with N>=2, exactly N indexed bonds, J>=0, zero field, and finite beta>0; no general long-range or higher-dimensional claim.

Statement

Use the real exponential exp⁡(x)=∑n≥0xn/n! of The real exponential function and the number e by a power series and the hyperbolic functions cosh⁡x=(exp⁡x+exp⁡(−x))/2, sinh⁡x=(exp⁡x−exp⁡(−x))/2. Define tanh⁡x:=sinh⁡x/cosh⁡x; the denominator is positive for every real x. Let N≥2, let ΩN={−1,1}N with cyclic indexing sN+1=s1, and count exactly the N indexed periodic bonds (so N=2 has two bonds). For coupling J≥0 in joules, zero field, and β>0 in J−1, adopt the finite canonical law PN(s)=ZN−1exp⁡ ⁣(βJ∑i=1Nsisi+1). Set t=tanh⁡(βJ)∈[0,1). For 1≤r≤N−1, EN[si]=0,EN[sisi+r]=tr+tN−r1+tN. For each fixed positive integer separation r and N→∞, the correlation tends to tr, with absolute error at most tN−r+tN+r. When J>0, define the limiting correlation length ξ>0 in lattice spacings as the unique number satisfying tr=e−r/ξ for every positive integer separation r; equivalently, ξ=−1/log⁡t. The zero-field susceptibility per spin, defined by χN:=N−1∂hEN,h[∑isi]∣h=0 for the field-coupled law with energy − ⁣J∑isisi+1−h∑isi, where the field parameter h∈R is in joules, has limit χ∞=β1+t1−t. Both χN and χ∞ have units J−1 per spin (the spin variables are dimensionless). Define the finite-ring model free energy by FN:=−β−1log⁡ZN and its limiting free energy per spin, when the limit exists, by f∞:=lim⁡N→∞FN/N. Then f∞=−β−1log⁡(2cosh⁡(βJ)), analytic at every finite positive temperature. No statement is made for arbitrary one-dimensional interactions.

Facts & Assumptions

Given: The finite spin space ΩN, the stated bond convention, N≥2, J≥0, and finite β>0.

[F1]

Declared example preparations and thermodynamic interpretation: each finite model uses its explicitly adopted canonical probability on the stated reference/counting space; the preparation is a model assumption.

[F2]

One-dimensional periodic Ising-chain limit: for the periodic nearest-neighbor ring with N≥2, finite real J,h and β>0, the same cosh⁡ and sinh⁡ definitions give Ts,t=exp⁡(βJst+βh(s+t)/2) and ZN=tr⁡(TN). Its explicit eigenvalues satisfy λ+>∣λ−∣, and the logarithmic partition limit and each fixed finite order of parameter derivatives converge uniformly on compact subsets of (0,∞)×R2. This is the current published mathematical supplier for the exact indexed-bond model used here.

[F3]

The exponential is positive and satisfies exp⁡(−x)=1/exp⁡(x): exp⁡(x)>0 for every real x, hence the defined cosh⁡x is positive and tanh⁡x is well-defined.

[F4]

The exponential function is strictly increasing: exp⁡ is continuous and x≤y implies exp⁡(x)≤exp⁡(y).

[F5]

A real-analytic function on an open subset of R is locally represented by a convergent real power series: a real-valued function on an open interval is real analytic when each point has a neighborhood on which the function equals a convergent real power series.

[F6]

The sum of a real power series is real analytic throughout the open interval determined by its radius: the sum of a real power series is real analytic on the interior of its radius of convergence; the exponential series has infinite radius.

[F7]

Real-analytic functions are closed under sums, products and compositions, and under quotients where the denominator is nonzero: sums, products, quotients with nonzero denominator, and compositions of real-analytic functions on their open domains are real analytic.

[F8]

[The power series for log(1+x) on (-1,1], including the Abel endpoint](/item/thm-log-one-plus-x-power-series): for ∣u∣<1, log⁡(1+u)=∑n=1∞(−1)n+1un/n, so u↦log⁡(1+u) is real analytic near u=0.

[F9]

Order, continuity, range, and the product, quotient, and reciprocal laws for the natural logarithm: for positive x,y, log⁡(xy)=log⁡x+log⁡y.

[A1]

For the response clause, the field-coupled finite law is PN,h(s)∝exp⁡(βJ∑isisi+1+βh∑isi) for real field parameter h in joules; this is a finite model family, not an empirical magnetic law.

Proof

1.1givenF2F4algebra

In the zero-field case the transfer matrix in [F2] is Ts,t=eKst with K=βJ, so T=(eKe−Ke−KeK). Direct multiplication gives orthonormal eigenvectors u+=(1,1)/2, u−=(1,−1)/2 with eigenvalues λ+=2cosh⁡K and λ−=2sinh⁡K. Hence ZN=tr⁡(TN)=λ+N+λ−N=λ+N(1+tN), where t=λ−/λ+=tanh⁡K. Since K≥0, [F4] gives eK≥e−K, so λ−=eK−e−K≥0; [F2] gives λ+>∣λ−∣. Thus t∈[0,1).

2.1A1F1F2step 1.1algebra

By the finite preparation [F1], the partition sum is the normalization for the adopted law. Put D0=diag⁡(1,−1); then D0u+=u− and D0u−=u+. Expanding the finite canonical sum over the N cyclic indices with sisi+r inserted gives numerator tr⁡(D0TrD0TN−r). Evaluating this trace in the eigenbasis yields λ−rλ+N−r+λ+rλ−N−r. Dividing by ZN from step 1.1 gives the stated correlation formula.

3.1step 2.1algebra

The map s↦−s preserves every bond product and reverses ∑isi, so at zero field EN[si]=0. For fixed r, subtracting tr from the formula of step 2.1 gives (tN−r−tN+r)/(1+tN); its absolute value is at most tN−r+tN+r because 0≤t<1. Thus the fixed-distance limit is tr. When 0<t<1, writing tr=e−r/ξ gives ξ=−1/log⁡t; for J=0, t=0 and all nonzero-separation correlations vanish.

3.2A1step 2.1algebra

For the field-coupled finite law [A1], differentiating its finite numerator and denominator sums gives ∂hEN,h[M]=β(EN,h[M2]−EN,h[M]2), where M=∑isi. At h=0, spin-flip symmetry gives EM=0 and translation invariance gives 1NE[M2]=1+∑r=1N−1E[s1s1+r]=1+2∑r=1N−1tr1+tN. The finite geometric sum tends to 2t/(1−t) as N→∞, so χN→β(1+t)/(1−t). This also gives χN=β when J=0.

4.1step 1.1step 3.1algebra

From step 1.1, N−1log⁡ZN=log⁡(2cosh⁡(βJ))+N−1log⁡(1+tN). Since 1≤1+tN≤2, the second term tends to zero, so f∞=−β−1log⁡(2cosh⁡(βJ)). The correlation length diverges only in the excluded zero-temperature limit βJ→∞.

5.1F3F4F5F6F7F8F9step 4.1∎

Fix any β0>0 and put A(β):=2cosh⁡(βJ) and u(β):=A(β)/A(β0)−1. By [F5, F6, F7], A and u are real analytic near β0; A(β0)>0 by [F3] and u(β0)=0. Continuity from [F4] allows a smaller neighborhood where ∣u(β)∣<1. By [F8], [F6], and the product law [F9], log⁡A(β)=log⁡A(β0)+log⁡(1+u(β)) is real analytic there, and [F7] makes f∞(β)=−β−1log⁡A(β) real analytic since β≠0. The relation β=(kBT)−1 is real analytic for T>0, so the limit is analytic at every finite positive temperature as well.

Remarks

The transfer-matrix premise uses the published exact periodic-chain theorem. The local spin-insertion and finite-response arguments prove the additional correlation and susceptibility formulas. No Gaussian-tail research token is used as a supplier.

Physics theoremStatement: Literature-sourcedProof: AI-adaptedprecheck passOpen item page →

Exact Tonks hard-rod pressure and convergent virial series

Physical scope

Finite-N one-dimensional classical hard rods of common length a>=0 with L>Na and quadratic momenta; the Euclidean Lebesgue product and Jordan-volume calculations assume ACω; exact thermodynamic and virial statements only for this Tonks model, not an all-fluid cluster expansion.

Statement

Assume ACω as defined in The Axiom of Countable Choice (ACω) for the Euclidean Lebesgue product and Jordan-volume identifications below. Use the natural logarithm log⁡ of The natural logarithm as the inverse of the exponential function. Let N≥1, let a≥0 be a hard-rod length in metres, and let L>Na be the one-dimensional box length in metres. Rod centers are xi∈[a/2,L−a/2] with ∣xi−xj∣≥a for i≠j. Let each rod have mass m>0 kg, momentum pi∈R in kg m s−1, and Hamiltonian H=∑i=1Npi2/(2m) on the allowed configuration set, with infinite energy outside it. Adopt the reference dx dp/(N!h0N), where h0>0 is a fixed action in J s, and canonical inverse temperature β>0 in J−1. Put λ=h0β/(2πm). Then ZN=(L−Na)NN!λN,U=EH=N2β,p1D:=−(∂F∂L)N,a,m,β,h0=Nβ(L−Na),F=−β−1log⁡ZN. Here p1D is the one-dimensional force conjugate to length and has units N; it is not a three-dimensional pressure. With number density ρ=N/L, the exact equation is βp1D=ρ1−aρ=∑j=1∞aj−1ρj,aρ<1. After the first M≥1 terms, the exact nonnegative remainder (strictly positive when a>0) is ρ(aρ)M/(1−aρ), bounded by ρrM/(1−r) whenever aρ≤r<1. At fixed ρ>0 with aρ<1, as N,L→∞ with N/L=ρ, lim⁡βFN=log⁡ ⁣(ρλ1−aρ)−1, with finite-size excess between 0 and (1+log⁡N)/N.

Facts & Assumptions

Given: ACω, N≥1, a≥0, L>Na, m>0, β>0, and the hard-core position set and Hamiltonian in the Statement.

[F1]

Declared example preparations and thermodynamic interpretation: for each specified model with finite positive partition function, the canonical density is Z−1e−βH relative to its fixed reference, and F=−β−1log⁡Z is the adopted free-energy definition.

[F2]

Gaussian, gamma and beta product integrals: for m,β>0, ∫Re−βp2/(2m)dp=2πm/β and the normalized Gaussian has second moment m/β; the positive gamma integral is finite.

[F3]

Change of variables for an injective C1 map on a compact Jordan set: an injective C1 map with nonzero Jacobian on an open neighborhood of a compact Jordan set obeys the compact Jordan change-of-variables formula.

[F4]

The natural logarithm as the inverse of the exponential function defines log⁡ on positive reals; Order, continuity, range, and the product, quotient, and reciprocal laws for the natural logarithm makes it continuous and increasing with the product and quotient laws, and The natural logarithm has derivative 1/x and equals the integral from 1 to x of 1/t gives log⁡′(x)=1/x and log⁡x=∫1xdt/t. Finite iteration of the product law gives the integer-power identity.

[F6]

Tonelli factors nonnegative product integrals and Fubini factors absolutely integrable signed product integrals on sigma-finite measure spaces (Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, Fubini's theorem for L^1 functions on a sigma-finite product).

[F7]

A bounded set is Jordan measurable when its boundary has content zero, and, for m≥1, a continuous graph over a closed nondegenerate rectangle in Rm has content zero in Rm+1 (A bounded set in Rm is Jordan measurable iff its boundary is null, equivalently of content zero, The graph of a continuous function on a closed nondegenerate rectangle in Rm has content zero in Rm+1).

[F8]

The Gaussian moments in [F2] are improper Riemann integrals; on bounded intervals they agree with Lebesgue integrals under ACω, and monotone convergence passes nonnegative truncations to the full line (Riemann–Lebesgue comparison for distribution test integrands, Monotone convergence for the integral).

[A1]

The finite hard-core configuration rule and indistinguishable-classical-particle reference convention are the adopted Tonks model assumptions; no interaction beyond hard exclusion is included.

Proof

1.1givenF2F5F6F8algebra

The single-rod momentum Lebesgue integral is ∫Re−βp2/(2m)dp=2πm/β by [F2] and [F8]. By [F5] and Tonelli [F6], the N momentum factors integrate to the product of these integrals. Dividing by h0N gives λ−N.

1.2givenF2F8algebra

By [F2] and the bounded Riemann-to-Lebesgue agreement in [F8], each normalized scalar momentum Gaussian has second moment m/β. Thus the mean kinetic energy of each rod is 1/(2β) and the total mean is U=N/(2β).

1.3givenA1F3F5F7F9algebra

Put ℓ=L−Na>0. For N=1, the allowed position set is [a/2,L−a/2], whose Lebesgue measure is L−a=ℓ directly by the box formula [F5]; no graph theorem is needed. For N≥2, the allowed labelled position set is a finite union of closed bounded ordering polytopes in [a/2,L−a/2]N, hence compact by [F9]. Its boundary is contained in finitely many box faces and contact hyperplanes xi−xj=±a, each contained in a graph of an affine function over a sufficiently large closed nondegenerate (N−1)-dimensional coordinate rectangle, so [F7] makes it Jordan measurable. Fix one label ordering π and sort its centers as xπ(1)≤⋯≤xπ(N). On that sector set yi=xπ(i)−a/2−(i−1)a ; the hard-core inequalities become the closed bounded, hence compact, ordered simplex Sℓ:={0≤y1≤⋯≤yN≤ℓ} by [F9]. Its boundary is contained in the box faces and finitely many tie hyperplanes, each contained in such an affine graph over a nondegenerate (N−1)-dimensional rectangle; [F7] therefore makes Sℓ Jordan measurable. The affine permutation map has absolute determinant 1, so [F3] preserves Jordan volume. The N! coordinate permutations of Sℓ partition [0,ℓ]N up to the tie hyperplanes, which have content zero by [F7]; their equal volumes sum to ℓN, hence each simplex has volume ℓN/N!. The N! labelled order sectors cover the allowed set and overlap only on the stated content-zero contact hyperplanes. Therefore its Jordan content is ℓN, and [F5] gives the full labelled position Lebesgue measure ℓN.

2.1F1F5F6step 1.1step 1.3algebra

The reference measure divides the position-momentum integral by N!h0N. Product Lebesgue measure and Tonelli [F5, F6] separate the position measure ℓN from step 1.3 and the momentum integral from step 1.1. Thus ZN=ℓN/(N!λN).

3.1F1F4step 2.1algebra

By the product, quotient and integer-power laws in [F4], F=−β−1[Nlog⁡((L−Na)/λ)−log⁡N!], using the logarithm derivative in [F4] at fixed N,a,β gives p1D=N/[β(L−Na)]. Substituting ρ=N/L yields βp1D=ρ/(1−aρ).

4.1step 3.1algebra

At a=0, the virial series has first term ρ and all higher terms zero, so its remainder after any M≥1 terms is zero. For 0≤aρ<1, the finite geometric identity ∑j=1M(aρ)j−1=(1−(aρ)M)/(1−aρ) gives the stated virial partial sum and exact remainder after multiplication by ρ. If aρ≤r<1, then (aρ)M/(1−aρ)≤rM/(1−r), proving the uniform bound.

5.1F4step 2.1algebra∎

At fixed ρ>0, L−Na=N(1/ρ−a)>0, so the laws in [F4] give βFN=−log⁡((L−Na)/λ)+log⁡N!N. For each integer k=1,…,N−1, monotonicity of log⁡ bounds its integral on [k,k+1] between log⁡k and log⁡(k+1). Summing gives ∫1Nlog⁡x dx≤log⁡N!≤∫1Nlog⁡x dx+log⁡N. For N=1 the sums and integral are empty and zero. Since [F4] gives (xlog⁡x−x)′=log⁡x, the integral is Nlog⁡N−N+1. Hence log⁡N!=Nlog⁡N−N+ϵN with 1≤ϵN≤1+log⁡N. Substitution and [F4] cancel log⁡N, yielding βFN=log⁡ ⁣(ρλ1−aρ)−1+ϵNN. The excess is nonnegative and at most (1+log⁡N)/N. Finally, log⁡N=∫1Ndx/x≤∫1Ndx/x=2(N−1), so that upper bound tends to zero and proves the stated thermodynamic limit.

Remarks

The infinite virial series is the geometric expansion of this exact one-dimensional hard-rod equation of state. It is not a formal or convergent cluster expansion for arbitrary fluids, attractions, dimensions, or interactions.

Physics theoremStatement: Literature-sourcedProof: AI-adaptedprecheck passOpen item page →

Finite spherical-shell gas and its exact marginal

Physical scope

Finite-N classical ideal kinetic energy for one operationally indistinguishable species, conditional on Countable Choice for the exact shell marginal, with a specifically adopted rotational Gaussian-direction preparation; no dynamical ergodicity or general coarea claim.

Statement

Assume ACω as defined in The Axiom of Countable Choice (ACω), the mathematical hypothesis required by the spherical-marginal supplier. Let N,d≥1, D=dN≥2, and let B=∏j=1d[0,Lj]⊂Rd have side lengths Lj>0 in metres and volume V=∏jLj in md. Let m>0 kg and E>0 joules, and set R=2mE in kg m s−1. Take one operationally indistinguishable species of N particles as specified by Classical particle counting and physical observables. On the auxiliary labelled cover, take independent uniform positions qi∈B with density V−1 and, independently, the momentum law P=RZ/∣Z∣, where Z∈RD has the independent standard-normal coordinate law defined in Exact spherical-shell coordinate marginal and fixed-block Gaussian limit and P=0 is assigned on the null event Z=0. This product law is invariant under particle permutations, so it defines the stated model preparation on the quotient; a one-particle marginal means the law of a uniformly selected particle. This is a specified spherical-shell preparation: H(P)=∣P∣2/(2m)=E almost surely, with no assumption about a trajectory or equilibration. For any momentum block of k coordinates with 1≤k<D, its density is fD,k from Exact spherical-shell coordinate marginal and fixed-block Gaussian limit. In particular the total-energy variance is zero J2.

Compare the separate canonical ideal-gas preparation of Exact canonical and grand ideal-gas preparations with the same N,d,m,B and one-species convention, at β=D/(2E)>0. Its mean total energy is also E, but its variance is D/(2β2)>0 J2. If N>1, the uniformly selected particle's d-component shell marginal has compact support because d<D, whereas its canonical Gaussian momentum has positive probability outside every finite-radius ball. The two finite preparations therefore differ despite their matched mean energy. For N=1, the full momentum vector is singular on its sphere and the density formula for k=D is not asserted.

Facts & Assumptions

Given: ACω, the rectangular box with side lengths in metres and volume in md, the one-species operational indistinguishability assumption, particle number and dimension, positive mass and energy, shell radius R, and the explicitly adopted independent position/spherical-momentum preparation in the Statement.

[F1]

Declared example preparations and thermodynamic interpretation: each shell or product law is an explicitly adopted probability preparation, separate from a canonical law and not inferred from Hamiltonian dynamics.

[F2]

Exact spherical-shell coordinate marginal and fixed-block Gaussian limit: assuming ACω, for D≥2, 1≤k<D, and R>0, the independent standard-normal coordinate law has an orthogonally invariant direction; RZ/∣Z∣ is concentrated on ∣P∣=R, its k-coordinate marginal has the stated exact density, and its coordinate second moment is R2/D.

[F3]

Exact canonical and grand ideal-gas preparations: for the massive d-dimensional canonical ideal gas with β>0, total energy has mean dN/(2β) and variance dN/(2β2), and each component of the uniformly selected particle's momentum is Gaussian with variance m/β.

[F4]

Classical particle counting and physical observables: for one operationally indistinguishable species, the labelled lift carries the symmetric permutation action and a permutation-invariant one-particle marginal is interpreted as the law of a uniformly selected particle.

[A1]

This permutation-invariant spherical-shell law is a chosen finite classical distribution on one operationally indistinguishable species; it is not a claim that an isolated gas dynamically reaches it or that it is the only microcanonical convention.

Proof

1.1givenF1A1algebra

Under the shell preparation [F1] and its explicit law [A1], ∣P∣=R almost surely by the definition of P; hence H(P)=R2/(2m)=E on every realization outside the null event Z=0. A constant random variable has variance zero.

2.1F2step 1.1

For 1≤k<D, the first k momentum coordinates are the projection in [F2], so their exact Lebesgue density is the displayed fD,k, which vanishes outside the closed radius-R ball. The event that this block has norm R requires the remaining D−k Gaussian coordinates to vanish; their joint density and D−k≥1 imply that this event has probability zero.

2.2givenF3step 1.1algebra

Choose β=D/(2E). By [F3], the canonical mean energy is D/(2β)=E and its variance is D/(2β2)>0. Thus it cannot have the shell law of step 1.1, whose energy variance is zero.

3.1A1F2F3F4step 2.1∎

If N>1, then d<D. Particle permutations act by orthogonal coordinate permutations on RD, so [F2] and [A1] make every particle block marginal identical; by the selected-particle convention [F4], their uniform mixture is the law of a uniformly selected particle. It has a density supported in the closed radius-R ball, with no boundary atom by step 2.1. By [F3], the canonical Gaussian of a uniformly selected particle has strictly positive density at every point of Rd, so its integral over the open annulus R<∣p∣<2R is positive, while the shell probability of that annulus is zero. Therefore the one-particle marginals differ as well. For N=1, d=D and the full shell law is singular, so no k=D density assertion has been used.

Remarks

This result compares two stated finite probability laws. It proves no equivalence or inequivalence theorem for general thermodynamic ensembles and no dynamical equilibration claim.

5 · Examples, experiments and counterexamples

None yet.

Sources