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.

Mixture Dynamics And Shells

1 · Prerequisites

2 · Summary

This page develops conditional mathematical results for mixtures, finite-site shells, continuous ideal gases, regular elastic collision models, a specified Ornstein–Uhlenbeck process, Hamiltonian volume dynamics, discrete quantum dephasing, and compact regular energy shells. Each result states the regularity, integrability, choice, measure, trace-class, or dynamical assumptions it uses.

For mixtures, compactness of the generator space and continuity of the energy give attainment and a finite-support optimizer; noncompact attainment needs additional coercivity or closedness hypotheses. The finite-site microcanonical result uses the exact width condition δN≥q∑j∣ej∣ and δN/N→0. The continuous gas calculation uses the declared classical measure dq dp/(N!hdN) and a normalized-Gaussian preparation on the exact kinetic-energy sphere; it does not establish equivalence for finite-width shells or interacting particles.

The collision entropy inequality is conditional on the symmetric weak collision identity and the stated regular solution class. The identity is a model premise; it is not derived here from reversible many-particle mechanics. Hamiltonian volume preservation gives recurrence only almost everywhere in finite invariant measure. It does not imply ergodicity, and the time-average conclusion requires ergodicity of the particular sampled map and an L1 observable.

The Ornstein–Uhlenbeck result is a mathematical Brownian model under the stated filtration hypothesis. A physical Einstein relation or measured noise spectrum is not asserted. Quantum Cesàro dephasing is restricted to the batch 11 trace-class Gibbs eigenbasis interface and a discrete spectrum with finite multiplicities; its proof remains conditional on that current draft supplier. The shell definition supplies a candidate graph-weighted integral on a regular level whenever its graph partition is supplied. Under AC and the compact regular-band hypotheses, the shell theorem proves chart independence, coarea and positive finite mass on nonempty levels. Invariance additionally requires the stated Hamiltonian diffeomorphism on the full band. Singular levels and normalization on noncompact bands require separate arguments.

The declared earlier page prerequisites include td-math-ensembles-and-limits, with the earlier chain through td-math-nonsmooth-equilibrium to td-math-ordered-composition. The current carriers and owner Step-1 readiness records supply this authored context. Those availability records do not replace independent review of the precise supplier hypotheses and arguments.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

DefinitionDefinition: Literature-sourcedProof: Not applicableverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Compact phase-generator mixture

Definition

Let K be a nonempty compact subset of Rd, d≥1, with its Borel sigma-algebra, and let a:K→R be continuous. The convex hull conv⁡K is the set of finite convex combinations of points of K. A phase generator is a point y∈K with declared energy a(y). A mixture is a Borel probability measure P on K. Its barycenter and mean energy are the coordinatewise Lebesgue integrals yˉ(P):=∫Ky dP(y)∈Rd,aˉ(P):=∫Ka(y) dP(y)∈R. These integrals are finite because each coordinate and a are bounded on compact K. For prescribed mean x∈conv⁡K, the mixture energy is amix(x):=inf⁡{aˉ(P):P is a Borel probability on K, yˉ(P)=x}. The admissible set is nonempty since every finite convex combination of points of K is the barycenter of a finite atomic probability. The definition does not assert that the infimum is attained; the theorem on this page gives a compact continuous sufficient condition.

Scope

The definition concerns a mathematical probability mixture and an assigned scalar energy. It does not assert that physical phases can be prepared, exchanged, or brought to equilibrium.

DefinitionDefinition: Literature-sourcedProof: Not applicableverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Elastic collision model and weak collision identity

Definition

Fix d≥2. The elastic collision set is Qel:={(v,w,v′,w′)∈R4d:v+w=v′+w′, ∣v∣2+∣w∣2=∣v′∣2+∣w′∣2}, with its relative Borel sigma-algebra. A collision measure is a nonnegative countably additive measure ν on this measurable space. It is symmetric if it is invariant under the exchanges (v,w)↔(w,v), (v′,w′)↔(w′,v′), and (v,w)↔(v′,w′). These invariances are hypotheses on ν, not consequences of energy and momentum conservation.

Fix the velocity measure space (Rd,L(Rd),dv), where dv is the Lebesgue set function of Lebesgue measurable sets, the family L(Rn), and the restricted set function λn and its countably additive measure-space structure is assumed as part of the data. This definition makes no existence assertion for that structure or for the collision operator. An input density is a finite-valued positive Lebesgue-measurable function f:Rd→(0,∞) whose four coordinate evaluations f(v),f(w),f(v′),f(w′) on Qel are ν-measurable. Let Dν be a declared set of such inputs. Put Lloc1:={q:Rd→R:q is Lebesgue measurable and ∫K∣q∣ dv<∞ for every compact K⊂Rd}. Its elements here are chosen function representatives, rather than equivalence classes. A collision model includes a specified map Qν:Dν→Lloc1; write Q(f):=Qν(f).

For f∈Dν, put A=f(v)f(w) and B=f(v′)f(w′). For a finite-valued complex Lebesgue-measurable function φ on velocity space, put δφ:=φ(v)+φ(w)−φ(v′)−φ(w′). Define the complete admissible test class Tf to consist of exactly those φ for which δφ is ν-measurable and ∫Rd∣φ(v)Qν(f)(v)∣ dv<∞,∫Qel∣(B−A)δφ∣ dν<∞. Absolute integrability is in the sense of Integrable real and complex functions, and their integrals. The model also declares, for each f∈Dν, a complex linear testing subspace Sf⊆Tf containing at least one nonzero test. Require the following weak collision identity for every f∈Dν and every φ∈Sf. A formulation adopting all admissible tests specifies Sf=Tf; a result using only particular tests names those tests and their linear span explicitly: ∫Rdφ(v)Q(f)(v) dv=14∫Qel(f(v′)f(w′)−f(v)f(w))(φ(v)+φ(w)−φ(v′)−φ(w′)) dν, Both sides are absolutely convergent by the test-class definition. Admissible tests need not be bounded or compactly supported: 1, velocity coordinates, squared speed, and log⁡(f/fref) for a positive constant fref belong to Tf whenever they meet the displayed measurability and integrability conditions. A consumer invoking one of these tests must establish or assume those conditions and require it to belong to the declared testing domain Sf. In particular, the entropy branch uses the span of 1 and log⁡(f/fref), with absolute convergence explicitly assumed. The map Qν, its input domain and its testing domains are adopted model data; the weak identity is their admissibility condition, not an assertion that every positive density and collision measure admit such a map or that these tests alone uniquely determine it. Molecular chaos or a microscopic derivation of the identity is not part of the definition.

Scope

The measure's normalization is chosen so Q(f) has units of f per unit time in an application. No specific cross-section or measured collision rate is asserted.

DefinitionDefinition: Literature-sourcedProof: Not applicableverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Finite-site canonical law and microcanonical shell

Definition

Fix a finite alphabet A={1,…,q}, q≥1, and real single-site energies e1,…,eq in a fixed energy unit. At inverse parameter β∈R, define Z1(β):=∑j=1qe−βej,pj(β):=e−βej/Z1(β),E(β):=∑j=1qpj(β)ej. For N≥1, a configuration i=(i1,…,iN)∈AN has additive energy HN(i):=∑r=1Neir. Given a width δN≥0, the shell is ΓN(β,δN):={i∈AN:∣HN(i)−NE(β)∣≤δN}. When this finite set is nonempty, its microcanonical law means the uniform probability on ΓN, and its dimensionless entropy is log⁡∣ΓN∣. A type is a vector t=(n1/N,…,nq/N), where the nonnegative integer counts sum to N; its multiplicity is N!/∏jnj!.

The convention fixes both the shell center and the law within the shell. It makes no physical claim that an apparatus has resolution δN.

DefinitionDefinition: Literature-sourcedProof: Not applicableverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Ornstein–Uhlenbeck integral equation

Definition

Let (Ω,F,(Ft)t≥0,P) carry a standard Brownian motion B adapted to the filtration, with future increments independent of the present sigma-algebra. Fix constants γ>0 and σ>0, and an F0-measurable real random variable X0 with finite second moment. A continuous adapted solution of the Ornstein–Uhlenbeck equation on [0,T] is a continuous adapted process X satisfying, almost surely for every t∈[0,T], Xt=X0−γ∫0tXs ds+σBt, where the time integral is the ordinary Lebesgue integral along the continuous path. In differential notation this is dXt=−γXt dt+σ dBt, with the stochastic differential interpreted through the displayed integral equation.

The definition asserts neither existence nor uniqueness. The Itô integral in the explicit representation is the square-integrable predictable-process integral of Ito integral for square-integrable predictable processes.

DefinitionDefinition: Literature-sourcedProof: Not applicableverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Regular energy shell and weighted area measure

Definition

Let O⊂Rm, m≥2, be open, let H∈C2(O,R), and let J⊂R be open. An energy E∈J is regular if ∇H(x)≠0 for every x∈ΣE:=H−1({E}). By the implicit-function theorem, a neighborhood of each point of a nonempty regular level is a graph xi=h(z).

In a graph chart, use the local area density 1+∣∇h(z)∣2 dz, where dz is (m−1)-dimensional Lebesgue measure. A graph partition on ΣE is a locally finite graph-chart cover with nonnegative smooth functions χa supported in their respective charts and summing to one on the level. Whenever such a partition is supplied, define the candidate weighted shell integral by μE(g):=∑a∫g(z,ha(z))χa(z,ha(z))∣∂iH(z,ha(z))∣−1 dz for nonnegative Borel g. Its independence from the chosen chart partition, and its equality with ∫ΣEg∣∇H∣−1dAE, are proved in the theorem on this page. If 0<μE(1)<∞, the microcanonical shell probability is μ^E:=μE/μE(1).

This is an explicit chart-based definition, not a formal delta distribution δ(H−E). Critical levels are excluded.

TheoremStatement: Literature-sourcedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Continuous ideal-gas bulk limit and finite-block bridge

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Fix dimension d≥1, particle mass m>0, action reference h>0, inverse energy β>0, and positive cube volumes VN with N/VN→ρ∈(0,∞). For N indistinguishable classical noninteracting particles in a cube of volume VN, use the phase-space reference dq dp/(N!hdN) and Hamiltonian HN(q,p)=∑i=1N∣pi∣2/(2m). Define λβ=hβ/(2πm). Then the dimensionless partition integral and free-energy limits are ZN=VNNN!λβdN,−β−1log⁡ZNN⟶β−1[log⁡(ρλβd)−1],−β−1log⁡ZNVN⟶β−1ρ[log⁡(ρλβd)−1]. There is also the following finite-momentum-block limit. Let M=dN, let G1,…,GM be independent standard normal variables, set AN=M−1∑a=1MGa2, and define canonical momenta Pacan=m/β Ga and shell-prepared momenta Pash=m/β Ga/AN. On the null event AN=0, set the shell-prepared vector to (mM/β,0,…,0). The shell-prepared vector has exactly energy dN/(2β), and for every fixed r and bounded continuous ψ:Rr→R, Eψ(P1sh,…,Prsh)⟶Eψ(P1can,…,Prcan). This is weak convergence for each fixed block in this free-particle model. The shell preparation is the stated normalized-Gaussian construction, and no finite-width shell or interacting-gas equivalence is asserted.

Facts & Assumptions

Given: The free-particle Hamiltonian, phase-space normalization, and normalized-Gaussian shell preparation.

[A1]

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)), as required by the Lebesgue-measure construction and change-of-variables suppliers below. d,N≥1, m,h,β,VN>0, N/VN→ρ∈(0,∞), and the N-particle Hamiltonian is additive kinetic energy in a cubic position box.

[A2]

For the block limit, the Ga are independent with density (2π)−1/2e−x2/2 and r is fixed while N→∞.

[F2]

The product measure of two sigma-finite measure spaces and For sigma-finite factors, the product measure exists, has the rectangle formula, is sigma-finite, and is unique: finite products of the sigma-finite Lebesgue factors carry the rectangular product measure. The Euclidean Lebesgue measure is the completion of the product of the factor Lebesgue measures identifies its completion with Euclidean Lebesgue measure, so the displayed Borel densities have the same integrals in either convention.

[F3]

Tonelli's theorem for nonnegative measurable functions on a sigma-finite product and Fubini's theorem for L^1 functions on a sigma-finite product: nonnegative product integrals factor by Tonelli; integrable products may be evaluated iteratively.

[F4]

The Gaussian integral ∫−∞∞e−x2 dx=π: the improper Riemann Gaussian integral is ∫−∞∞e−x2 dx=π. A nonnegative improper Riemann integral on a half-line agrees with the Lebesgue integral identifies each convergent nonnegative improper integral on [0,∞) with its Lebesgue integral. Apply it to the continuous even Gaussian and its even polynomial moments; reflection gives the negative half-line. Thus the improper integral evaluations below also evaluate their Lebesgue integrals.

[F7]

If u,v are differentiable on [a,b] with u′,v′ integrable, then ∫abuv′=u(b)v(b)−u(a)v(a)−∫abu′v: on each finite interval, differentiable factors with integrable derivatives obey the ordinary integration-by-parts identity.

[F8]

The exponential dominates every fixed nonnegative integer power at +∞: Rke−R2/2→0 for each fixed nonnegative integer k, by setting x=R2/2 and bounding xk/2 by a larger integer power.

Proof

Given: [A1] and, for the second conclusion, [A2].

1.1A1F1F2F3F4F5

The integrand is nonnegative, so [F3] factors the position and momentum integrals: ZN=1N!hdN∫ΛNN∫RdNe−β∑i∣pi∣2/(2m) dp dq=VNNN!hdN(∫Re−βp2/(2m) dp)dN. By [F4]–[F5] the one-coordinate integral is 2πm/β, proving the stated exact formula with λβ=hβ/(2πm). The dimensions are consistent: λβ is a length, so VNN/λβdN and ZN are dimensionless.

1.2A2F1F2F3F4F5F7F8algebra

Let G have the density in [A2]. By [F4]–[F5], it is normalized. The polynomial and Gaussian factors have continuous derivatives on every finite interval, so apply [F7] on [−R,R] first with u(x)=x, v′(x)=xe−x2/2, and then with u(x)=x3, v′(x)=xe−x2/2. In each case v(x)=−e−x2/2. Letting R→∞, the endpoint terms vanish by [F8], giving EG2=1 and EG4=3. Thus E(G2−1)=0 and E(G2−1)2=2. The independence in [A2] and [F3] factor the cross integrals, so expanding the square gives E(AN−1)2=M−2∑a=1M2=2/M. For any η>0, η21{∣AN−1∣>η}≤(AN−1)2, so integration proves P(∣AN−1∣>η)≤2/(Mη2)→0. Hence AN→1 in probability.

2.1A1F6step 1.1algebra

Monotonicity of log⁡ from [F6] gives ∫1Nlog⁡x dx≤log⁡N!≤∫1N+1log⁡x dx, hence log⁡N!=Nlog⁡N−N+O(log⁡(N+1)). Again by [F6], 0≤log⁡x=∫1xdt/t≤2x for x≥1, so log⁡(N+1)/N≤2N+1/N→0. Substituting the exact ZN yields N−1log⁡ZN=log⁡(VN/N)+1−dlog⁡λβ+o(1)→−log⁡ρ+1−dlog⁡λβ. Multiplication by −β−1 proves the per-particle limit; multiplying by N/VN→ρ proves the per-volume limit.

2.2A2step 1.2

Since AN>0 almost surely, the shell-prepared formulas define random variables on a probability-one set. On the null event AN=0, assign the fixed shell vector (mM/β,0,…,0); this measurable extension leaves every expectation unchanged and has the same stated energy. On AN>0, algebra gives ∑a=1dN(Pash)22m=12βAN∑aGa2=M2β=dN2β. This exact energy constraint holds on every outcome where AN>0, an event of probability one since the joint Gaussian density assigns zero mass to the origin.

2.3A2F1F3F9step 1.2

For fixed r, the scalar multiplier AN−1/2 tends to 1 in probability by continuity at 1 and step 1.2. Therefore (P1sh,…,Prsh)−(P1can,…,Prcan) tends to zero in probability: restrict first to the event ∑a=1rGa2≤R2, where its norm is at most m/βR∣AN−1/2−1∣, and then bound the complement by P(∑a=1rGa2>R2)≤r/R2 using the second moment. For bounded continuous ψ, choose R so this tail is small and additionally restrict to ∣AN−1/2−1∣≤1, whose complement has probability tending to zero. On this event both blocks lie in the ball of radius 2m/βR. Use uniform continuity of ψ on the compact ball of radius 2m/βR, and then use convergence in probability of the difference. For any ϵ>0, choose δ>0 by that uniform continuity. The expectation difference is at most ϵ+2∥ψ∥∞[r/R2+P(∣AN−1/2−1∣>1)+P(∥P1:rsh−P1:rcan∥>δ)], with the multiplier assigned arbitrarily on its null exceptional event. Sending first N→∞ and then R→∞ proves the claimed expectation limit.

3.1A1A2F1step 2.3∎

If positions are included, choose them independently and uniformly in the same cube under both preparations. After coordinatewise rescaling by the side length, every fixed position block already has the same uniform law in both models and is independent of the momentum vector. The same bounded-continuous test argument as in 2.3, first on compact sets and then on their small tails, proves the joint fixed-block limit. This does not change the energy computation because the Hamiltonian has no position dependence.

TheoremStatement: Literature-sourcedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Hamiltonian volume, recurrence, and ergodic scope

Statement

Assume the Axiom of Choice. Let (M,ω) be a smooth 2n-dimensional symplectic manifold and H∈C∞(M,R), with Hamiltonian field defined by ιXHω=dH. Wherever its local flow ϕt is defined, it preserves the symplectic form and the volume measure m determined by ωn/n!. No completeness is asserted. Let R be a measurable region on which a nonzero-time map T=ϕτ is defined and is forward invariant, T(R)⊆R. The restriction T:R→R preserves m∣R if T(R)=R, or more generally if m(R∖T(R))=0. This latter condition holds automatically when m(R)<∞. Thus for every finite-measure forward-invariant R and every measurable E⊆R, almost every point of E returns to E infinitely often under forward iteration. Mere forward invariance of an infinite-volume region does not imply relative measure preservation.

If m(R)=1 and T is ergodic, then for every f∈L1(R,m) the averages N−1∑j=0N−1f(Tjx) converge almost everywhere and in L1 to ∫Rf dm. Without ergodicity, the limit is the invariant conditional expectation supplied by Birkhoff's theorem, not necessarily the space mean. Finally, fix any real scalar kB. If each stated time map ϕt is a diffeomorphism of a finite invariant region R, a density p0≥0 with ∫Rp0 dm=1 and ∫R∣p0log⁡p0∣ dm<∞ transports by pt=p0∘ϕ−t with conserved fine-grained entropy −kB∫Rptlog⁡pt dm, using 0log⁡0=0. Volume preservation alone does not imply ergodicity or a time-average/ensemble identity.

Facts & Assumptions

Given: The symplectic/Hamiltonian definitions, the local flow domain, and each separately stated finite invariant probability system.

[A1]

AC is adopted; (M,ω) is smooth symplectic of dimension 2n, H is smooth, and XH follows ιXHω=dH.

[A2]

The local flow is used only where defined. The sampled region is measurable and forward invariant; for an infinite-measure region the restricted-measure assertion additionally requires m(R∖T(R))=0 (in particular, onto invariance suffices). For entropy transport, each stated ϕt is a diffeomorphism of R with inverse ϕ−t. For recurrence and Birkhoff, R is forward invariant, T=ϕτ and every iterate exist on R, and m(R)<∞; for the space-mean conclusion m(R)=1 and T is ergodic in the sense of [F5].

[A3]

The entropy scale kB is a fixed real mathematical parameter, with no physical value or unit adopted here. The transported density is measurable, nonnegative, normalized, and has finite absolute entropy integral, with 0log⁡0=0.

[F1]

The Axiom of Choice, Symplectic form and symplectic manifold, Hamiltonian vector field and Hamiltonian function, and Symplectic manifolds have a canonical orientation and volume form: AC is the inherited measure-theorem convention; ω is closed and nondegenerate, ιXHω=dH, and ωn/n! determines the canonical positive volume measure.

[F2]

Hamiltonian flows preserve the symplectic form: every local Hamiltonian flow preserves ω wherever defined. Liouville volume preservation then gives preservation of the symplectic volume on the same domain.

[F3]

Measure-preserving transformations and systems and Integral invariance under measure-preserving maps: the stated invariant flow map is measure preserving, and integrals of integrable functions are unchanged by its measure-preserving substitution.

[F4]

Poincare recurrence for finite measure-preserving systems: on a finite-measure system, almost every point of each measurable set returns to that set infinitely often.

[F5]

Ergodicity relative to an invariant measure, Birkhoff pointwise ergodic theorem, and Birkhoff ergodic theorem for ergodic finite-measure systems: ergodicity means every invariant event is null or conull; for an integrable observable on a finite measure-preserving probability system, Birkhoff supplies almost-everywhere convergence. Ergodic averages converge in Lp on finite-measure spaces at p=1 supplies L1 convergence on finite measure spaces; Finite-measure identification of the Birkhoff limit under AC identifies the limit as the unique integrable function measurable for I={E:T−1E=E} with ∫Ef∗ dm=∫Ef dm for every E∈I. This is the invariant conditional expectation meant here; the ergodic case has constant limit equal to the space integral.

[F6]

Probability measures and probability spaces and Integrable real and complex functions, and their integrals: normalized densities and the stated absolute-integrability condition define the probability and entropy integrals.

Proof

Given: [A1]–[A3] and the inspected supplier interfaces [F1]–[F6].

1.1A1A2F1F2F3algebra

By [F1], ω is closed and ιXHω=dH. The Hamiltonian-flow theorem [F2], under its local existence domain, gives ϕt∗ω=ω and hence ϕt∗(ωn/n!)=ωn/n!. The volume measure defined by this positive form is therefore preserved by each local flow diffeomorphism on its domain. For a measurable E⊆R, injectivity and volume preservation give m((T∣R)−1E)=m(E∩T(R)). If m(R∖T(R))=0, this equals m(E), so the restriction is the measure-preserving system of [F3]. When m(R)<∞, volume preservation gives m(T(R))=m(R); because T(R)⊆R, finite additivity gives m(R∖T(R))=0. Thus the finite forward-invariant case follows without an onto hypothesis. Also dH(XH)=ω(XH,XH)=0 because a two-form is alternating, so each defined orbit stays on its energy level. This argument grants no flow outside its domain.

2.1A2A3F3F6step 1.1

Apply [F3] to pt=p0∘ϕ−t. It preserves ∫pt dm=∫p0 dm=1. The scalar function u↦ulog⁡u, defined to be 0 at u=0 by [A3], is Borel; its absolute integrability is preserved because the same measure-preserving substitution gives ∫R∣ptlog⁡pt∣ dm=∫R∣p0log⁡p0∣ dm<∞. Applying integral invariance to its positive and negative parts proves ∫Rptlog⁡pt dm=∫Rp0log⁡p0 dm. Multiplication by the fixed real scalar −kB gives the stated scaled entropy conservation.

2.2A2F3F4step 1.1

The map T on R is measure preserving by 1.1 and [A2]. Its total measure is finite, so [F4] gives infinite returns for almost every point in each measurable E⊆R. The conclusion is explicitly almost everywhere and applies to the sampled map at the specified τ.

3.1A2F5step 1.1∎

Under the additional probability and ergodicity assumptions in [A2], the hypotheses of [F5] hold for each f∈L1(R,m). Birkhoff therefore gives almost-everywhere and L1 convergence to the invariant conditional expectation, and its ergodic corollary identifies that limit with ∫Rf dm. If ergodicity is removed, only the invariant conditional expectation conclusion remains. Neither [F2] nor recurrence [F4] supplies ergodicity.

TheoremStatement: Literature-sourcedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Attainment and finite support for compact continuous mixtures

Statement

Let K⊂Rd be nonempty compact, d≥1, and a:K→R continuous. Define amix as in Compact phase-generator mixture. Then for every x∈conv⁡K the infimum defining amix(x) is attained by a Borel probability measure supported on at most d+1 points. The function amix is lower semicontinuous and convex on conv⁡K, and is the greatest convex function on conv⁡K bounded above by a on K. Minimizing over all Borel probabilities on K gives the same attained value as minimizing over finite mixtures. This is a mathematical compactness result; it makes no phase-preparation or equilibrium postulate.

Facts & Assumptions

Given: The compact generator space, continuous energy, and Borel-probability mixture convention.

[A1]

K⊂Rd is nonempty compact, d≥1, and a:K→R is continuous.

[F1]

Compact phase-generator mixture: a mixture is a Borel probability on K; its barycenter and mean energy are coordinatewise Lebesgue integrals; amix is their constrained infimum.

[F2]

Borel probability law on a polish space: P is a countably additive Borel measure with P(K)=1. Integrable real and complex functions, and their integrals defines integrability by ∫∣f∣ dP<∞. For bounded Borel f with ∣f∣≤M, Monotonicity and nonnegative homogeneity of the nonnegative integral and the constant simple integral from The integral of a nonnegative simple function, identified with the nonnegative integral by The nonnegative integral agrees with the simple integral on simple functions, give ∫∣f∣ dP≤∫M dP=MP(K)=M, so f is integrable. The Lebesgue integral is linear on L1(μ) supplies linearity on these integrable functions. In particular, if ∣f∣≤ϵ, integrating the nonnegative functions ϵ+f and ϵ−f and using linearity gives ∣∫f dP∣≤ϵ.

[F3]

Heine-Borel in Rn: with the Euclidean metric a subset of Rn is compact if and only if it is closed and bounded, and the proof by bisection uses no choice principle; the same holds on the real line: in finite-dimensional Euclidean space, a set is compact exactly when it is closed and bounded; in particular, compact sets are closed and bounded.

[F4]

A continuous real-valued function on a nonempty compact metric space is bounded and attains a greatest and a least value: a continuous real function on a nonempty compact metric space attains its minimum and maximum.

Proof

Given: The data in [A1] and the mixture definitions in [F1].

1.1A1F3F4algebra

Define Φ(y)=(y,a(y))∈Rd+1 and G=Φ(K). The map Φ is continuous. Each coordinate of Φ is continuous on compact K, so [F4] bounds it; hence G is bounded. To see G is closed, take (y,r)∉G. If y∉K, closedness of K gives a neighborhood of y disjoint from K. If y∈K, then r≠a(y), and continuity of a gives neighborhoods of y and r whose product misses the graph of a. Thus the complement of G is open. By [F3], G is compact.

2.1A1F3step 1.1algebra

In a finite convex combination, discard zero-weight terms. If more than d+2 positive-weight points gi∈G⊂Rd+1 remain, the vectors (gi,1)∈Rd+2 are linearly dependent; hence there are coefficients ci, not all zero, with ∑ici=0 and ∑icigi=0. Both signs occur. Replacing the positive weights θi by θi−tci, where t=min⁡ci>0θi/ci, preserves their sum and combination, keeps them nonnegative, and makes at least one zero. Repeating gives at most d+2 points. Pad shorter combinations by zero weights and repeat any fixed point of nonempty G. The parameter set Gd+2×Δd+2 is closed and bounded in a finite-dimensional Euclidean space, hence compact by [F3]. The map Θ:Gd+2×Δd+2→Rd+1, Θ((gi),θ)=∑iθigi, is continuous and has image exactly conv⁡G: every open cover of the image pulls back to an open cover of the compact parameter set and therefore has a finite subcover. Thus conv⁡G is compact and closed by [F3].

3.1A1F1F3F4step 2.1

Fix x∈conv⁡K. Its fiber Fx={r:(x,r)∈conv⁡G} is nonempty, since each finite convex combination of generators gives such a point. It is closed and bounded because conv⁡G is closed and bounded by [F3]; hence it is compact by [F3]. Therefore [F4] gives a minimum rx. A convex-combination representation of (x,rx) by points (yi,a(yi)) defines a finite atomic Borel probability Px with barycenter x and energy rx.

4.1A1F1F2F3step 1.1step 2.1step 3.1

Conversely, for any Borel probability P on K, the vector-valued map Φ is bounded and coordinatewise integrable by [F2]. Its integral is in conv⁡G: for each ϵ>0, cover compact G by finitely many open Euclidean balls of radius ϵ/2 and pull the cover back to K by Φ. Disjointifying the finite preimages in order gives a finite Borel partition of K whose pieces have image diameter at most ϵ. Choose a point in each nonempty piece and weight it by that piece's P-mass. The resulting finite convex combination differs coordinatewise from ∫Φ dP by at most ϵ; as conv⁡G is closed by step 2.1, the integral belongs to it. If yˉ(P)=x, then (x,aˉ(P))∈conv⁡G, so aˉ(P)≥rx. Together with step 3.1 this proves that the all-measure infimum equals rx and is attained by Px.

5.1A1F3F4step 2.1step 3.1step 4.1

To prove lower semicontinuity, let xn→x in conv⁡K and put L=lim inf⁡namix(xn). The values lie between min⁡Ka and max⁡Ka, attained by [F4], so L is finite. By the definition of liminf, choose a subsequence nj with amix(xnj)→L. Each pair (xnj,amix(xnj)) lies in the closed set conv⁡G by step 4.1, and the pairs converge to (x,L); hence (x,L)∈conv⁡G. The minimum in the fiber at x satisfies amix(x)≤L=lim inf⁡namix(xn).

5.2F1step 3.1step 4.1

Combine minimizers Px and Pz with weights s and 1−s, 0≤s≤1. The combined finite measure has barycenter sx+(1−s)z and energy samix(x)+(1−s)amix(z). By step 4.1, the minimum at the combined barycenter is no greater than this value. Hence amix is convex.

5.3A1F1step 3.1step 4.1algebra

The support reduction in 2.1 first gives at most d+2 atoms. If an optimal representation has more than d+1 positive weights, the points yi∈Rd have an affine dependence ∑iciyi=0, ∑ici=0. For sufficiently small positive t, both weight lists θi±tci remain nonnegative and preserve the barycenter. If ∑icia(yi)≠0, one of these lists strictly lowers energy, contradicting optimality; hence this sum is zero. Increase t in one sign to the first value at which a positive coefficient becomes zero. The new weights preserve barycenter and energy and still attain the minimum; iteration leaves at most d+1 positive weights.

6.1F1step 5.2step 5.3algebra∎

Finally, let g be any convex function on conv⁡K with g(y)≤a(y) for y∈K. For the finite optimal measure from 4.1, finite Jensen (obtained by iterating the two-point convexity inequality) gives g(x)≤∑iθig(yi)≤∑iθia(yi)=amix(x). Since amix(y)≤a(y) for y∈K by the point mass at y, it is itself such a convex minorant and dominates every other one.

TheoremStatement: Literature-sourcedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Finite-site microcanonical and canonical block equivalence

Statement

Let A={1,…,q}, q≥1, have real energies ej, and fix finite β∈R. Set Z1=∑je−βej, pj=e−βej/Z1>0, and E=∑jpjej. For each N≥1, let δN≥q∑j=1q∣ej∣ and δN/N→0. Define the shell ΓN={i∈AN:∣∑r=1Neir−NE∣≤δN} and the uniform law on that shell. For probability laws μ,ν on the finite set Ak, total variation distance means 12∑x∈Ak∣μ(x)−ν(x)∣. Then ΓN is nonempty for every N; for each fixed k≥1, the first k-site marginal (for N≥k) converges in that distance to p⊗k; and 1Nlog⁡∣ΓN∣⟶H(p):=−∑j=1qpjlog⁡pj. The conclusion is for this finite alphabet, additive energy, exact shell width and uniform shell law. It does not assert finite-size equality or equivalence for interacting or continuous models.

Facts & Assumptions

Given: The finite-site canonical law and uniform shell convention.

[A1]

A is finite, q≥1, ej∈R, β∈R, and δN≥q∑j∣ej∣ with δN/N→0.

[F1]

Finite-site canonical law and microcanonical shell: configurations have additive energies, the canonical probabilities are pj=e−βej/Z1, and the shell law is uniform when the shell is nonempty.

[F2]

Finite sums and finite products, by recursion and Order, continuity, range, and the product, quotient, and reciprocal laws for the natural logarithm: the finite sums and logarithms obey the usual finite distributive and product identities.

[F3]

The real exponential function and the number e by a power series and The natural logarithm as the inverse of the exponential function: exponential weights are positive and logarithm is the inverse of the positive real exponential.

[F4]

The natural logarithm has derivative 1/x and equals the integral from 1 to x of 1/t: log⁡ is continuous on (0,∞), has derivative 1/x, and log⁡x=∫1xdt/t.

[F5]

Heine-Borel in Rn: with the Euclidean metric a subset of Rn is compact if and only if it is closed and bounded, and the proof by bisection uses no choice principle; the same holds on the real line: the finite-dimensional probability simplex is closed and bounded, hence compact; closed bounded subsets of finite-dimensional Euclidean space are compact.

[F6]

A continuous real-valued function on a nonempty compact metric space is bounded and attains a greatest and a least value: a continuous real-valued function on a nonempty compact metric space attains its minimum and maximum.

Proof

Given: [A1] and the conventions in [F1].

1.1A1F2F3F4algebra

For a type tj=nj/N, the number of configurations of that type is M(n)=N!/∏jnj!. The increasing function log⁡x gives ∫1mlog⁡x dx≤log⁡(m!)≤∫1m+1log⁡x dx(m≥1). Evaluating the antiderivative xlog⁡x−x and using log⁡(1+1/m)≤1/m shows log⁡(m!)=mlog⁡m−m+O(log⁡(m+1)+1) with an absolute constant; the case m=0 uses 0!=1 and 0log⁡0=0. Subtracting the q denominator estimates from the numerator estimate gives a constant Cq such that NH(t)−Cqlog⁡(N+1)−Cq≤log⁡M(n)≤NH(t)+Cqlog⁡(N+1)+Cq, where H(t)=−∑jtjlog⁡tj. The sharper upper bound M(n)≤eNH(t) follows by expanding 1=(∑jtj)N and retaining the nonnegative terms with counts n. For later limits, [F4] gives 0≤log⁡x=∫1xdt/t≤2x for x≥1, hence log⁡(N+1)/N≤2N+1/N→0.

1.2A1F1F2F3F4

For any probability vector t, define D(t∥p)=∑j:tj>0tjlog⁡(tj/pj). From [F4], log⁡u≥1−1/u for u>0, with equality only at u=1: apply log⁡z≤z−1 to z=1/u and negate. Therefore D(t∥p)≥∑j:tj>0(tj−pj)=1−∑j:tj>0pj≥0. Since every pj>0, equality forces every index into the support and tj/pj=1 for each j, so t=p. Also log⁡pj=−βej−log⁡Z1, hence D(t∥p)=−H(t)+β∑jtjej+log⁡Z1. In particular, among types of exact mean energy E, entropy is at most H(p), with equality only at p.

2.1A1F1F4step 1.1

Round down Npj for j<q and put the remainder in nq. Then ∣nj/N−pj∣≤1/N for j<q and ∣nq/N−pq∣≤(q−1)/N. Thus ∣∑jnjej−NE∣≤∑j<q∣ej∣+(q−1)∣eq∣≤q∑j∣ej∣≤δN. This type belongs to ΓN, proving nonemptiness. Its frequencies tend to p. For 0<x≤1, [F4] gives 0≤−xlog⁡x=xlog⁡(1/x)≤2x, by integrating 1/t≤1/t on [1,1/x]; hence −xlog⁡x extends continuously by 0 at x=0 and H(t)→H(p). By the lower bound in 1.1, the shell contains at least exp⁡(NH(p)−o(N)) configurations.

3.1A1F1F4F5F6step 1.1step 2.1step 1.2

Fix ϵ>0 and let Δq={t∈[0,1]q:∑jtj=1} and C={t∈Δq:∣t−p∣≥ϵ}. The simplex is closed and bounded, and C is closed in it, so both are compact by [F5]. If C=∅, the claim is vacuous. Otherwise let F={t∈C:∑jtjej=E}, a closed subset of C. If F=∅, the continuous function t↦∣∑jtjej−E∣ has a positive minimum on C by [F6], so choose τ smaller than this minimum and any η>0. If F≠∅, [F6] gives a maximum hF of H on F; step 1.2 gives hF<H(p). Put η=(H(p)−hF)/2 and let U={t∈C:H(t)<H(p)−η}, a relative open neighborhood of F. The compact set C∖U, if nonempty, has a positive minimum m of ∣∑jtjej−E∣ by [F6]; choose 0<τ<m. If C∖U=∅, choose any τ>0. In either case, every t∈C with ∣∑jtjej−E∣≤τ satisfies H(t)≤H(p)−η. For large N, δN/N≤τ; therefore every shell type outside the ϵ-ball has entropy at most H(p)−η. There are at most (N+1)q types, and 1.1 bounds the number of configurations of each by exp⁡[NH(t)]. By the log-growth estimate in 1.1, this type-count factor is subexponential, so division by the shell lower bound in step 2.1 shows that the shell probability outside the ball tends to zero exponentially. Thus the random type converges in probability to p.

4.1A1F1F5step 3.1algebra

Conditional on a type t, a uniform configuration is a uniformly ordered list of its N labeled positions. Couple the first k draws without replacement with k independent uniform position draws: if the independent positions are distinct, keep them; if not, replace them by an independently chosen uniform ordered distinct k-tuple. The resulting tuple is uniform among distinct positions and agrees with the original draws except on the collision event. A union bound over the (k2) pairs gives collision probability at most k(k−1)/(2N). Hence the type-conditional label marginal differs in total variation from t⊗k by at most this quantity. Since t→p in probability, for each fixed tuple (i1,…,ik)∈Ak the random variable XN=∣∏r=1ktir−∏r=1kpir∣ tends to zero in probability and is bounded by 1. For every ϵ>0, EXN≤ϵ+P(XN>ϵ), so EXN→0; summing over the finite set Ak proves the marginal converges to p⊗k in total variation.

5.1A1F1F4F5step 1.1step 1.2step 2.1algebra∎

The rounded shell type and lower bound in step 2.1 give lim inf⁡NN−1log⁡∣ΓN∣≥H(p). For any shell type t, step 1.2 gives the exact identity H(t)=H(p)+β(∑jtjej−E)−D(t∥p)≤H(p)+∣β∣δN/N. There are at most (N+1)q types, and step 1.1 bounds each multiplicity by eNH(t). Therefore ∣ΓN∣≤(N+1)qexp⁡ ⁣(NH(p)+∣β∣δN). By the log-growth estimate in 1.1, qlog⁡(N+1)/N→0. Since also δN/N→0, the limsup is at most H(p). Together with the lower bound this proves the entropy-density limit, including β=0, without invoking a nonexistent compactness step.

TheoremStatement: Literature-sourcedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Ornstein–Uhlenbeck equation, transition law, and fluctuations

Statement

Assume the Axiom of Choice The Axiom of Choice. Assume the filtered probability space, Brownian increment hypothesis, and initial data of Ornstein–Uhlenbeck integral equation, including EX02<∞ and future increments independent of F0. For γ,σ>0, on every finite time interval the process Xt=e−γtX0+σe−γt∫0teγs dBs has a continuous adapted version, solves the stated Ornstein–Uhlenbeck integral equation, and is the unique continuous solution driven by the same (X0,B). Put V=σ2/(2γ). For 0≤s<t, its transition law conditional on Fs is normal with mean e−γ(t−s)Xs and variance V(1−e−2γ(t−s)). Consequently EXt=e−γtEX0,Var⁡(Xt)=e−2γtVar⁡(X0)+V(1−e−2γt). From any such initial law, Xt converges in distribution to N(0,V). If X0∼N(0,V) and is independent of future Brownian increments, the process is stationary Gaussian with covariance C(τ)=Ve−γ∣τ∣, and its transition kernel satisfies detailed balance with the N(0,V) density. Its transition density solves the Fokker–Planck equation for positive time. Under the convention S(ω)=(2π)−1∫RC(τ)eiωτ dτ, the stationary spectrum and time-average bound are S(ω)=Vγπ(γ2+ω2),Var⁡ ⁣(1T∫0TXt dt)≤2VγT. All these statements concern the specified linear stochastic model. An Einstein relation between σ and a physical bath temperature would be an additional physical postulate.

Facts & Assumptions

Given: The Axiom of Choice, the exact Brownian filtration hypothesis, and the OU equation definition.

[A1]

B is adapted and each increment Bt−Bs is independent of Fs with law N(0,t−s); X0 is F0-measurable, has finite second moment, and is independent of future Brownian increments.

[A2]

γ>0, σ>0, and time is restricted to an arbitrary finite interval [0,T].

[F1]

Ornstein–Uhlenbeck integral equation: a solution is a continuous adapted process satisfying Xt=X0−γ∫0tXsds+σBt almost surely for every t∈[0,T] on each finite interval [0,T], with T>0.

[F2]

Brownian motion, Existence of continuous Brownian motion, and Elementary predictable Brownian integrands: under the assumed Axiom of Choice, Brownian paths have a common continuous version and deterministic continuous square-integrable integrands satisfy the standing predictable Brownian hypotheses on finite intervals.

[F3]

Ito integral for square-integrable predictable processes and Ito isometry and linearity in predictable L2: under the assumed Axiom of Choice, the stochastic integral is defined for the deterministic square-integrable integrand and isometry computes its second moments.

[F4]

The Ito integral process has a continuous martingale version and Deterministic Ito integrals are Gaussian: under the assumed Axiom of Choice, the integral process has a continuous version, and deterministic-integrand integrals are centered Gaussian with variance equal to the squared-integrand integral.

[F5]

The exponential function is smooth and (exp⁡)′=exp⁡ and Fubini's theorem for L^1 functions on a sigma-finite product: the ordinary exponential derivative and finite continuous-integral Fubini identity hold.

[F7]

The Gaussian integral ∫−∞∞e−x2 dx=π: the centered Gaussian density with variance v>0 is normalized after the linear change of variable from the unit Gaussian.

Proof

Given: The Axiom of Choice, [A1]–[A2], and [F1]–[F7].

1.1A1A2F2F3F4F5

On [0,T], g(s)=eγs is deterministic, continuous, and square-integrable. Thus Wt=∫0tg(s) dBs is defined by [F3], and [F4] supplies a continuous version. For a partition 0=t0<⋯<tn=t, the deterministic step sums ∑kg(tk)(Btk+1−Btk) converge in L2 to Wt by the Itô isometry and uniform approximation of g. Discrete summation by parts rewrites them as g(t)Bt−g(0)B0−∑k=0n−1Btk+1(g(tk+1)−g(tk)) up to an endpoint term tending to zero. Brownian path continuity makes the last sums converge pathwise to ∫0tg′(s)Bsds. The L2 limit and this almost-sure limit agree in probability, so Wt=eγtBt−γ∫0teγsBs ds almost surely for each fixed t. Both sides have continuous versions; equality on rational t holds on one probability-one event and therefore extends to all t∈[0,T].

1.2A1A2F3F4

Split the stochastic integral at s and use its linearity: Xt=e−γ(t−s)Xs+σ∫ste−γ(t−u) dBu. Step-function approximants to the last deterministic-integrand integral depend only on increments after s, so [A1] makes them independent of Fs. Their L2 limit remains independent of Fs: for bounded Fs-measurable Y and bounded continuous φ, convergence in probability plus boundedness gives convergence of E[Yφ(⋅)], preserving the factorization. By [F4] the limit is centered Gaussian; [F3] computes its variance as σ2∫ste−2γ(t−u)du=V(1−e−2γ(t−s)). This proves the conditional transition law and the Markov property.

2.1A1A2F1F4F5step 1.1

Set Xt=e−γtX0+σe−γtWt. By [F4] it is continuous and adapted. Step 1.1 rewrites it as Xt=e−γtX0+σBt−σγ∫0te−γ(t−s)Bs ds. Integrating this continuous expression in t, ordinary Fubini [F5] gives ∫0tXsds=X0(1−e−γt)/γ+σ∫0tBsds−σ∫0tBu(1−e−γ(t−u)) du. Substitution into X0−γ∫0tXsds+σBt cancels the two ∫Bsds terms and leaves exactly the displayed expression for Xt. Thus [F1] holds.

2.2A1A2F3F4step 1.2

Taking expectations in the explicit solution, the stochastic integral has mean zero by [F4], giving EXt=e−γtEX0. The integral is independent of X0 and centered, so its variance adds to e−2γtVar⁡(X0); [F3] evaluates the added variance as V(1−e−2γt). The independent Gaussian term tends in law to N(0,V) because its variance tends to V and e−γtX0→0 in L2; for any bounded continuous test, restrict to a compact interval, use uniform continuity and the second-moment tail bound, then take the interval radius to infinity.

2.3A1A2F6F7step 1.2

Write r=e−γt and π(x)=(2πV)−1/2e−x2/(2V). The transition density is Kt(x,y)=[2πV(1−r2)]−1/2exp⁡ ⁣(−(y−rx)22V(1−r2)). Multiplying by π(x) gives a symmetric prefactor and exponent proportional to −[x2+y2−2rxy]/[2V(1−r2)], unchanged when x and y are exchanged. Hence π(x)Kt(x,y)=π(y)Kt(y,x), the detailed-balance identity.

2.4A2F5F7step 1.2

For t>0 let μ=e−γtx, v=V(1−e−2γt), and p(t,x,y) be the normal density with mean μ and variance v. Direct differentiation gives ∂tp=p(μ′(y−μ)v+v′2[(y−μ)2v2−1v]),μ′=−γμ,v′=σ2−2γv. Also direct y-differentiation gives γ∂y(yp)+(σ2/2)∂yyp equal to the same expression after writing y=μ+(y−μ). Thus ∂tp=γ∂y(yp)+(σ2/2)∂yyp. The Gaussian formula has total integral one and converges to point mass at x as t↓0; this is the stated positive-time Fokker–Planck solution.

3.1A2F1step 2.1algebra

If X and X~ are two continuous solutions with the same data, their difference z satisfies z(t)=−γ∫0tz(s)ds. Let D(t)=∫0t∣z(s)∣ds. Then D(0)=0, D≥0, and D′(t)=∣z(t)∣≤γD(t) by the equation. The derivative of e−γtD(t) is nonpositive, so this nonnegative function is at most its initial value zero. Hence D=0 and continuity gives z=0 at every time.

3.2A1A2F4step 1.2step 2.2

If X0 is centered normal with variance V, step 1.2 gives centered Gaussian finite-dimensional laws, since each stochastic integral over a disjoint increment interval is an independent Gaussian and X0 is independent of future increments. Its one-time variance is V. For s≤t, the split formula in 2.1 and independence of the later increment yield E[XsXt]=e−γ(t−s)E[Xs2]=Ve−γ(t−s). Thus covariance depends only on ∣t−s∣. A centered Gaussian finite-dimensional law is determined by its covariance matrix; shifting every time by the same amount preserves that matrix and proves stationarity.

4.1A2F6step 3.2

For R>0, the derivative and componentwise FTC in [F6] give ∫0Re−(γ−iω)tdt=(1−e−(γ−iω)R)/(γ−iω). The modulus formula gives ∣e−(γ−iω)R∣=e−γR→0 and ∣e−(γ−iω)t∣=e−γt, whose integral is finite. Thus the stationary covariance is integrable. By evenness, S(ω)=VπRe⁡∫0∞e−(γ−iω)tdt=VπRe⁡1γ−iω=Vγπ(γ2+ω2). The integral converges absolutely because γ>0.

5.1A1A2F5step 3.2∎

For the stationary time average, Fubini applies to the random product because E∣XsXt∣≤(EXs2EXt2)1/2=V, so its absolute integral over the finite square and probability space is at most VT2. Changing variables in the two triangular halves gives Var⁡ ⁣(T−1∫0TXtdt)=2VT2∫0T(T−s)e−γsds≤2VγT. The bound tends to zero, so these averages converge to zero in L2 and in probability. It asserts no almost-sure theorem for arbitrary observables.

TheoremStatement: Literature-sourcedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Conditional Boltzmann entropy inequality and equality form

Statement

Fix periods L1,…,Ld>0, let Ω=∏j=1d(R/LjZ) with quotient Lebesgue measure dr, let d≥2, and let I⊂R be an open time interval. Suppose f:Ω×Rd×I→(0,∞) is C1 in space and time and C2 in velocity. Fix a positive constant reference density fref with the same units as f, and put ℓ(r,v,t)=log⁡(f(r,v,t)/fref). For each t∈I, assume f, f∣v∣2, and f∣ℓ∣ are integrable over Ω×Rd. For every compact interval I0⊂I, assume there is a nonnegative GI0∈L1(Ω×Rd,dr dv) such that ∣(1+ℓ)∂tf∣≤GI0 almost everywhere for every t∈I0. Also assume (1+∣ℓ∣)∣v⋅∇rf∣∈L1(Ω×Rd,dr dv) for each t∈I. Write A=f(r,v,t)f(r,w,t) and B=f(r,v′,t)f(r,w′,t). For each t∈I, assume r↦∫Qel∣(A−B)(ℓ(r,v,t)+ℓ(r,w,t)−ℓ(r,v′,t)−ℓ(r,w′,t))∣ dν is measurable and belongs to L1(Ω,dr). Let Q and the symmetric collision measure ν satisfy the weak identity in Elastic collision model and weak collision identity for tests 1 and ℓ(r,⋅,t) at each (r,t), with absolute convergence, and suppose ∂tf+v⋅∇rf=Q(f). Then for H[f](t)=∫Ω∫Rdfℓ dv dr, ddtH[f](t)=−14∫Ω∫Qel(A−B)log⁡(A/B) dν dr≤0. The integrand is nonnegative because log⁡ is strictly increasing. Equality in the entropy-production integral holds exactly when A=B for ν-almost every collision at almost every r. Under the additional assumptions that ν is locally finite, has full support on Qel, and f is positive and C2 in v, this balance implies for almost every r that log⁡(f(r,v,t)/fref)=α(r,t)+b(r,t)⋅v+λ(r,t)2∣v∣2. If f(r,⋅,t) is integrable over Rd, then λ(r,t)<0, so its velocity profile is Maxwellian. Conversely, every profile of this form has zero collision entropy production. The weak collision identity is an explicit model hypothesis; this is not a derivation from reversible many-particle mechanics or a PDE existence theorem.

Facts & Assumptions

Given: The regular positive solution, the elastic collision identities, and the adopted symmetric weak collision model.

[A1]

d≥2, Ω=∏j(R/LjZ) has its quotient Lebesgue measure, f>0 has the stated Cr,t1 and Cv2 regularity, and fref>0 is a constant reference with matching units.

[A2]

The displayed moments and entropy are integrable. On each compact time interval a common L1(dr dv) majorant bounds ∣(1+ℓ)∂tf∣; for each time, (1+∣ℓ∣)∣v⋅∇rf∣ is integrable. The weak collision identity is absolutely convergent for tests 1 and ℓ, and the nonnegative entropy-production integrand has a measurable spatial integral in L1(dr) at each time.

[A3]

For the equality characterization only, ν is locally finite with full support on the elastic collision set; velocity integrability holds for almost every spatial point.

[F1]

Elastic collision model and weak collision identity: collisions conserve pair momentum and squared speed; the adopted weak identity has factor 1/4 and is symmetric under incoming/outgoing and particle exchange.

[F2]

Lebesgue measurable sets, the family L(Rn), and the restricted set function λn, The product measure of two sigma-finite measure spaces, For sigma-finite factors, the product measure exists, has the rectangle formula, is sigma-finite, and is unique, Tonelli's theorem for nonnegative measurable functions on a sigma-finite product, and Fubini's theorem for L^1 functions on a sigma-finite product: the spatial and velocity Lebesgue factors are sigma-finite, so their product integrals and absolutely convergent rearrangements are legitimate. These product-measure suppliers are not applied to the arbitrary collision measure ν; collision integration is performed first at each fixed spatial point.

[F3]

Dominated convergence: the stated integrable envelope permits differentiation of the entropy integral by dominated convergence.

Proof

Given: [A1]–[A3] and the weak identity and integration interfaces [F1]–[F5].

1.1A1A2F2F3F4algebra

For a fixed time, write ℓ=log⁡(f/fref). By [F3]–[F4], dominated differentiation gives dHdt=∫Ω∫Rd(1+ℓ) ∂tf dv dr. Substitute the transport equation. The spatial contribution is −∫Ω∫Rd(1+ℓ) v⋅∇rf dv dr=−∫Ω∫Rdv⋅∇r(fℓ) dv dr. For each fixed v, the integral of each coordinate derivative over the periodic rectangle is zero by the one-dimensional fundamental theorem of calculus and opposite-face identification; absolute integrability in [A2] permits Fubini [F2]. Thus the transport contribution vanishes.

2.1A1A2F1F2step 1.1

The weak identity with test 1 has zero bracket on every elastic collision, so ∫Q(f) dv=0 at each r,t. Consequently the collision contribution from (1+ℓ)Q(f) equals that from ℓQ(f). Apply the weak identity with φ=ℓ(r,⋅,t) at fixed (r,t). Since fref is constant, the sum ℓ(v)+ℓ(w)−ℓ(v′)−ℓ(w′) is log⁡(A/B). Hence ∫ℓ(v)Q(f)(v) dv=14∫Qel(B−A)log⁡(A/B) dν. The velocity-side function is measurable in (r,v) because Q(f)=∂tf+v⋅∇rf and f>0 have the stated regularity. Integration of its positive and negative parts over the sigma-finite Lebesgue velocity factor gives a measurable function of r; the weak identity makes that function equal to the collision integral at each (r,t). By [A2] this collision integral is in L1(Ω,dr), so integrate the identity in r to obtain the asserted derivative formula. no collision/spatial product measure or interchange of ν with spatial Lebesgue measure is used.

3.1A2F2F5step 2.1algebra

For positive A,B, strict increase of log⁡ [F5] makes (A−B)log⁡(A/B)≥0, with equality exactly when A=B. The derivative in step 2.1 is therefore nonpositive. Let g(r) be the nonnegative inner collision integral; it is measurable and integrable by [A2] and step 2.1. If ∫Ωg(r) dr=0, then g(r)=0 for Lebesgue-almost every r. For each such r, a nonnegative function with zero ν-integral vanishes ν-almost everywhere: its sets where the value is at least 1/n all have ν-measure zero, and their countable union contains its positive set. This argument holds for an arbitrary measure, with no sigma-finiteness hypothesis. Thus A=B for ν-almost every collision at almost every r. The converse follows by substitution in the derivative identity.

4.1A1A3F1step 3.1algebra

Under [A3], fix an r for which collision balance holds ν-almost everywhere. The function (v,w,v′,w′)↦A−B is continuous. If it were nonzero at an allowed collision, it would remain nonzero on a relative open neighborhood; full support of ν gives that neighborhood positive measure, contradicting balance. Hence balance holds for every elastic collision. The set of such r has full Lebesgue measure and is dense in Ω; continuity in r extends the identity to any limit point when needed.

5.1A1A3step 4.1algebra

Fix such an r and put g(v)=log⁡(f(r,v,t)/fref). For any a,u∈Rd and any orthogonal rotation R, the pairs (a+u,a−u) and (a+Ru,a−Ru) have the same total momentum 2a and total squared speed 2∣a∣2+2∣u∣2. Step 4.1 therefore gives g(a+u)+g(a−u)=g(a+Ru)+g(a−Ru). For unit vectors θ,η, put u=sθ and Rθ=η. The second-order Taylor expansion at a gives g(a+sθ)+g(a−sθ)=2g(a)+s2θ⊤D2g(a)θ+o(s2). Equating the expansion for η, dividing by s2 and sending s→0 shows that the quadratic form of the symmetric Hessian is constant on the unit sphere. Testing coordinate vectors makes all diagonal entries equal; testing (ei+ej)/2 makes each off-diagonal entry zero. Thus D2g(a)=λ(a)I.

6.1A1step 5.1algebra

Since ∂ijg=0 for i≠j, the fundamental theorem of calculus along coordinate segments shows ∂ig(v) depends only on vi. Write it as Gi(vi). The diagonal Hessian equalities give Gi′(vi)=Gj′(vj) for every pair of independent real coordinates vi,vj. Because d≥2, varying these coordinates independently forces all these derivatives to be one constant λ. Integrating each first derivative gives g(v)=α+b⋅v+λ∣v∣2/2.

7.1A1A3F1step 6.1algebra∎

If λ>0, eg(v) grows quadratically in every sufficiently large direction and is not integrable. If λ=0, eα+b⋅v is constant when b=0 and grows on a half-line when b≠0, so it is not integrable over Rd. Thus integrability forces λ<0. Completing the square writes f(r,v,t)=C(r,t)exp⁡[−c(r,t)∣v−u(r,t)∣2] for C,c>0, the Maxwellian velocity form. Conversely, when g=α+b⋅v+λ∣v∣2/2, the sum g(v)+g(w) is unchanged by every collision because total momentum and total squared speed are conserved; hence A=B and the production is zero. This equality branch concerns the collision term only; a spatially varying Maxwellian may still evolve under transport.

TheoremStatement: Literature-sourcedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Coarea and invariant measure on compact regular energy shells

Statement

Assume AC. Let n≥1, let O⊂R2n be open with its canonical symplectic form, let H∈C2(O,R), and let J⊂R be open. Suppose ∇H≠0 on H−1(J) and H−1(K) is compact in O for every compact K⊂J. Define μE by Regular energy shell and weighted area measure. Then for every nonnegative Borel f compactly supported in H−1(J), ∫Of(x) dx=∫JμE(f) dE. For every E∈J with nonempty shell ΣE=H−1({E}), 0<μE(1)<∞, and μ^E=μE/μE(1) is a probability measure. Let ϕt be a Hamiltonian flow map defined as a diffeomorphism on a full compact energy band containing ΣE and preserving that band. Then μ^E is invariant under ϕt for every such regular E. Critical levels and noncompact bands are outside this statement.

Facts & Assumptions

Given: The regular-band data, the weighted-area definition, and the complete Hamiltonian flow on the particular band.

[A1]

AC is adopted; O⊂R2n is open, H∈C2, ∇H≠0 on H−1(J), and compact energy intervals have compact inverse images in O.

[A2]

The nonnegative Borel function is compactly supported in H−1(J). For invariance, ϕt is a diffeomorphism on a compact band, maps it onto itself, and is the Hamiltonian flow of H.

[F1]

The Axiom of Choice and Regular energy shell and weighted area measure: AC is the standing hypothesis of the change-of-variables supplier; the candidate shell integral is the sum over a supplied graph partition, and graph area in a chart has density 1+∣∇h∣2 dz.

[F2]

The Euclidean implicit function theorem with derivative formula: if one partial derivative of H is nonzero, a neighborhood of that level is a C1 graph xi=h(z,E) with ∂Eh=1/∂iH.

[F3]

A C^1 diffeomorphism satisfies the change-of-variables formula for nonnegative Lebesgue measurable functions: a C1 injective change of variables with nonzero Jacobian transforms a nonnegative Lebesgue integral by the absolute Jacobian.

[F5]

Write x=(q,p)∈Rn×Rn. For the canonical form ω=∑jdqj∧dpj, the C1 field associated to H∈C2 is XH=(∇pH,−∇qH), so ιXHω=dH. Here its flow means the C1 local family solving ∂sϕs(x)=XH(ϕs(x)), ϕ0(x)=x, with inverse ϕ−s on its flow domains. The volume and energy invariance needed here are proved locally in step 1.2 at this C2 regularity, rather than taken from a smooth-only supplier.

[F6]

Dominated convergence: a common integrable chart bound justifies continuity in E of compactly supported shell integrals.

[F8]

Finite measures agreeing on a generating pi-system and on the whole space are equal: two finite measures that agree on a generating pi-system and have equal total mass agree on the generated sigma-algebra.

Proof

Given: [A1]–[A2] and the geometric and measure interfaces [F1]–[F8].

1.1A1A2F2F4construct

Let C be a compact energy band containing the support of f. At each x∈H−1(C), some coordinate derivative ∂iH(x) is nonzero. By [F2] a neighborhood is parameterized by the remaining coordinates z and energy E, with xi=h(z,E). Compactness gives a finite cover by such charts. Refine the cover by finitely many balls whose closures lie in chart domains and whose smaller concentric balls cover H−1(C). Choose standard smooth bump functions ψa supported in the larger balls and positive on the smaller balls, and let Ψ=∑aψa. Its minimum on H−1(C) is positive. Choose a smooth scalar cutoff η that is zero below one quarter of this minimum and one above one half of it. Then χa=η(Ψ)ψa/Ψ, extended by zero where Ψ=0, is smooth, supported in its chart, and the sum is one near H−1(C).

1.2A1A2F3F5algebra

At the stated C2 regularity, XH in [F5] has continuous first derivatives. Along every flow curve, ddsH(ϕs(x))=∑j(HqjHpj−HpjHqj)(ϕs(x))=0. Differentiating its integral equation ϕs(x)=x+∫0sXH(ϕu(x)) du in x gives Js=I+∫0sDXH(ϕu(x))Ju du, where Js=Dxϕs(x); on each compact flow segment the C1 derivatives are continuous, so difference quotients pass through this finite integral. Thus ∂sJs=DXH(ϕs(x))Js. Differentiating the determinant by multilinearity in its columns gives ∂sdet⁡Js=(tr⁡DXH)(ϕs(x))det⁡Js. Indeed Js is invertible by the inverse flow, and the sum of determinants with one column replaced by DXH times that column is det⁡Js tr⁡(Js−1DXHJs)=det⁡Js tr⁡DXH. Equality of the continuous mixed second derivatives of H gives tr⁡DXH=∑j(∂qj∂pjH−∂pj∂qjH)=0. Since det⁡J0=1, det⁡Js=1 throughout the flow domain; [F3] then gives Lebesgue-volume preservation. Both invariances therefore hold for the entire C2 scope of this item.

2.1A1A2F1F2F3F4step 1.1algebra

In each chart, the map (z,E)↦(z,ha(z,E)) is injective after shrinking the chart, and its absolute Jacobian is ∣∂Eha∣=∣∂iH∣−1. By [F3] and Tonelli [F4], ∫Of(x)χa(x) dx=∫J∫f(z,ha(z,E))χa(z,ha(z,E))∣∂iH∣−1 dz dE. Implicit differentiation gives ∂zjha=−∂jH/∂iH, so the graph area density is 1+∣∇zha∣2 dz=∣∇H∣∣∂iH∣ dz. The inner integral is therefore the corresponding chart contribution to the candidate μE(f) in [F1]. Sum the finitely many charts and use ∑aχa=1 on the support to obtain the asserted coarea identity.

2.2A1F1F2step 1.1algebra

On overlaps, if one graph parameterization is obtained from another by a C1 coordinate change ψ, its Gram matrix is transformed by G↦(Dψ)⊤G(Dψ). Taking determinants shows the area density changes by ∣det⁡Dψ∣, exactly the ordinary chart Jacobian. Implicit differentiation gives 1+∣∇h∣2/∣∇H∣=∣∂iH∣−1, so the chart contributions define the same measure independent of the partition and equal ∫ΣEg∣∇H∣−1dAE. Each shell is a closed subset of compact H−1({E}), hence compact. A finite collection of graph charts with ∇H bounded away from zero on that compact level gives finite area and hence μE(1)<∞. If the shell is nonempty, one graph chart contains an open patch; its coordinate domain has positive (2n−1)-dimensional measure and the continuous density ∣∇H∣−1dAE/dz is positive there, so μE(1)>0.

3.1A1A2F1F5F6step 1.2step 2.1

Fix a compact band and a compactly supported continuous test g on it. For a continuous compactly supported scalar function q on J, the coarea identity from step 2.1 and the Hamiltonian volume/energy invariance from step 1.2 give ∫Jq(E)μE(g∘ϕt) dE=∫Oq(H(x))g(ϕtx) dx=∫Oq(H(x))g(x) dx=∫Jq(E)μE(g) dE. The second equality is the change of variables by the volume-preserving diffeomorphism ϕt; H∘ϕt=H keeps the energy weight fixed. For each test g, both functions of E in the outer integrals are continuous locally: use the finite chart representation in step 2.1, compact chart support, the positive lower bound for ∣∇H∣, and dominated convergence [F6]. Their difference integrates to zero against every compactly supported continuous q. If the continuous difference were nonzero at one energy, it would keep one strict sign on a neighborhood; a nonnegative continuous bump q supported there would make the integral nonzero. Thus μE(g∘ϕt)=μE(g) for every E in the band.

4.1F5F6F7step 2.2step 3.1

Fix a closed subset F of a nonempty shell. If F=∅, both measures assign it zero. Otherwise let χ be a continuous compactly supported cutoff equal to one on the compact shell; the same finite-ball bump construction used in 1.1 supplies it. In the ambient band set gn(x)=χ(x)max⁡{0,1−n dist⁡(x,F)}. By [F7], gn is continuous, 0≤gn≤1, and its restriction to ΣE converges pointwise to 1F. Because ϕt maps the shell onto itself, gn∘ϕt converges there to 1ϕt−1F. Apply step 3.1 to each gn and use dominated convergence [F6] for the finite shell measure from step 2.2; this gives μE(ϕt−1F)=μE(F) for every closed F⊆ΣE. Closed sets form a pi-system generating the shell Borel sigma-algebra, and both μE and its pushforward have the same finite total mass because ϕt maps the shell onto itself.

5.1F8step 2.2step 4.1∎

By [F8], agreement on closed sets and equality of total masses imply that μE equals its pushforward on every Borel subset of ΣE. Divide by 0<μE(1)<∞ from step 2.2. Hence the normalized regular-shell probability is invariant for each such E, not merely almost every energy.

TheoremStatement: Literature-sourcedProof: AI-adaptedverified 2026-10-15 (gpt-6-luna, gpt-6.1-sol, and delegated Codex review/repair agents; see completed run evidence)Open item page →

Trace-norm Cesàro dephasing for a discrete Gibbs spectrum

Statement

Assume the Axiom of Choice and the exact first-clause interface of Trace-class Gibbs normalization for a semibounded self-adjoint operator: a nonzero complex Hilbert space has a finite or countable complete orthonormal energy eigenbasis (uj) for a self-adjoint semibounded A, with each energy eigenspace finite-dimensional. The supplier is a current draft whose statement and complete argument have been inspected; this theorem is conditional on that interface and does not claim its independent acceptance. Let ℏ>0, define Utuj=e−iEjt/ℏuj, and extend linearly and continuously. For any positive trace-class density operator ρ with Tr⁡ρ=1, set ρt=UtρUt∗. If PE denotes the projection onto the full eigenspace of energy E, then the trace-norm Riemann average ρˉT:=1T∫0Tρt dt converges as T→∞ to ρdiag:=∑EPEρPE in trace norm. The full degenerate energy blocks are retained. Thus expectations of every bounded observable converge to their diagonal-ensemble values. This asserts Cesàro convergence, not convergence of ρt at individual late times, and is restricted to the stated discrete trace-class spectrum.

Facts & Assumptions

Given: The inspected discrete eigenbasis interface, a positive trace-one trace-class density, and the unitary phase convention.

[A1]

AC is adopted; the actual draft supplier statement Trace-class Gibbs normalization for a semibounded self-adjoint operator assumes a nonzero complex Hilbert space, self-adjoint A with spectral support in [a,∞), β>0, and positive trace-class Gibbs operator Kβ. Its first clause supplies a finite or countable complete orthonormal basis of eigenvectors Auj=Ejuj, finite multiplicity for each energy, and a positive finite partition trace. Its complete proof is present as draft; the dependency is used conditionally.

[A2]

ρ is positive trace class with trace one; ℏ>0 is fixed and the phase is e−iEjt/ℏ.

[F1]

Trace class operator: trace-class operators have finite trace norm and form the space on which the trace and trace norm are defined.

[F2]

Trace class is a two sided Banach operator ideal: trace class is a Banach space, and ∥ARB∥1≤∥A∥∥R∥1∥B∥ for bounded A,B.

[F3]

Finite-rank orthogonal compressions converge in trace norm: if finite-rank orthogonal projections converge strongly to the identity on a separable Hilbert space, their two-sided compressions of a trace-class operator converge in trace norm.

[F4]

Trace of a positive operator is the sum of its eigenvalues and Trace is absolutely convergent and basis independent: positive trace is the sum of diagonal eigenbasis entries, and the trace is basis independent.

[F5]

Cyclicity of the trace: for trace-class S and bounded B, Tr⁡(SB)=Tr⁡(BS).

[F7]

Real and complex inner-product spaces and their induced length: the Hilbert inner product is linear in the first argument.

Proof

Given: [A1]–[A2] and [F1]–[F7].

1.1A1A2F6F7

Each multiplier e−iEjt/ℏ has modulus one by [F6]. Thus Ut preserves the norm of finite linear combinations of the complete orthonormal basis; it extends by density to an isometry on the Hilbert space. The similarly defined U−t is its inverse, so Ut is unitary and depends continuously on t on each finite-dimensional energy span. The arguments Ejt/ℏ are dimensionless.

2.1A2F2F5F7step 1.1

By [F2], ρt=UtρUt∗ is trace class. It is positive because ⟨ρtx,x⟩=⟨ρUt∗x,Ut∗x⟩≥0, and cyclicity [F5] gives Tr⁡ρt=Tr⁡ρ=1.

2.2A1A2F2F3step 1.1

Enumerate the distinct energy eigenspaces and let PN project onto the first N of them, with PN=I after the list is exhausted if it is finite. Their union spans the eigenbasis, so PN→I strongly. The Hilbert space is separable by the countable complete basis in [A1]. Therefore [F3] gives ∥PNρPN−ρ∥1→0. Unitary invariance of trace norm follows from [F2] in both directions: ∥UtRUt∗∥1≤∥R∥1 and R=Ut∗(UtRUt∗)Ut gives the reverse inequality. Consequently ∥UtPNρPNUt∗−ρt∥1=∥PNρPN−ρ∥1 uniformly in t.

3.1F2step 2.2

For fixed N, the matrix of UtPNρPNUt∗ on its finite-dimensional range has continuous entries, hence is continuous in trace norm. Step 2.2 makes ρt a uniform trace-norm limit on compact time intervals of these continuous maps; therefore t↦ρt is trace-norm continuous. On [0,T] it is uniformly continuous. Riemann sums for a uniformly continuous Banach-valued function are Cauchy: two sums are compared on a common refinement and their difference is bounded by T times the modulus of continuity at the larger mesh. Completeness in [F2] supplies the trace-class integral and its usual norm bound.

4.1A1A2F6F7step 3.1

In the eigenbasis, the (i,j) matrix entry of UtPNρPNUt∗ is ⟨ρuj,ui⟩ei(Ej−Ei)t/ℏ. Its time-average factor is 1 if Ei=Ej; if Ei≠Ej, direct integration using [F6] gives absolute value at most 2ℏ/(T∣Ej−Ei∣), which tends to zero. Since the compression is finite dimensional, entrywise convergence implies trace-norm convergence. Hence the finite compressed average tends to ∑E among the first NPEρPE.

5.1A1A2F1F4step 4.1

Each block PEρPE is positive trace class. The sum of their traces is ∑E∑uj∈ran⁡PE⟨ρuj,uj⟩=Tr⁡ρ=1 by [F4] and the complete basis in [A1]. Therefore their finite partial sums are Cauchy in trace norm: for a positive finite block tail its trace norm equals its trace, and the scalar trace tail tends to zero. Their limit ρdiag exists in trace norm, is positive, and has trace one.

6.1F2step 2.2step 3.1step 4.1step 5.1

Given ϵ>0, choose N so the compression error in step 2.2 and the diagonal tail in step 5.1 are each below ϵ/3. The norm bound for the Riemann integral in step 3.1 makes the average of the compression error at most ϵ/3, uniformly in T. For this fixed N, step 4.1 makes the finite compressed average differ from its finite diagonal limit by less than ϵ/3 for all sufficiently large T. The triangle inequality then gives ∥ρˉT−ρdiag∥1<ϵ.

7.1A2F2F4step 6.1∎

For bounded observable B, the trace-class ideal and trace estimates in [F2] and [F4] give ∣Tr⁡((ρˉT−ρdiag)B)∣≤∥ρˉT−ρdiag∥1∥B∥. Step 6.1 sends the right-hand side to zero. The conclusion is only convergence of time-averaged states in the stated trace norm.

5 · Examples, experiments and counterexamples

None yet.

Sources