Second Law of Thermodynamics (Wikipedia) I: Carnot's Theorem and the Thermodynamic Temperature ScaleTextbook
Motivation
The second law of thermodynamics is a physical law based on universal empirical observation concerning heat and energy interconversions. Two of its classical statements are impossibility statements about machines. Rudolf Clausius (1854): "Heat can never pass from a colder to a warmer body without some other change, connected therewith, occurring at the same time." Lord Kelvin: "It is impossible for any system to operate in a thermodynamic cycle and deliver a net amount of energy by work to its surroundings while receiving energy by heat transfer from a single thermal reservoir." The Wikipedia article Second law of thermodynamics records that these two have been shown to be equivalent.
The law is not only prohibitive. Sadi Carnot's analysis of heat engines (1824) yields an upper bound on the efficiency of any engine running between two heat reservoirs, attained by reversible engines and strictly missed by irreversible ones, and — crucially — the same for all reversible engines whatever the working substance. That substance-independence is what makes a thermodynamic temperature possible: as the article puts it, "The second law allows a distinguished temperature scale, which defines an absolute, thermodynamic temperature, independent of the properties of any particular reference thermometric body."
This mission formalizes that chain: from the two impossibility statements, through Carnot's theorem, to the existence and uniqueness of the absolute temperature scale normalized by the Kelvin convention at the triple point of water.
Setting
A cyclically operating device is recorded by the heat it exchanges with thermal reservoirs over one complete cycle. Reservoirs are labelled by empirical temperature, an element of a linearly ordered type Θ (s<t: the reservoir at s is colder). A cycle is a finitely supported function c:Θ→R, where c(t) is the net heat absorbed by the device from the reservoir at t (negative if rejected). Because the device returns to its initial state, the first law makes the net work delivered equal to the total heat absorbed:
W(c)=t∑c(t).
The cycles that can actually be realized form a set C, assumed to be a convex cone: the null cycle is realizable; two realizable devices can be run side by side, their heat exchanges adding; and a realizable device can be scaled by any factor r≥0. A cycle is reversible when c∈C and −c∈C, and operates betweenh and k when it exchanges heat with no other reservoir. A Carnot engine between a hot h and a cold k is a reversible cycle operating between them with c(h)>0, c(k)<0 and W(c)>0. The efficiency relative to the hot reservoir is η=W(c)/c(h).
In this language:
Clausius statement. For k<h: every realizable cycle operating between h and k with W(c)=0 satisfies c(k)≤0.
Kelvin statement. Every realizable cycle exchanging heat with a single reservoir satisfies W(c)≤0.
Cadmits Carnot engines when a Carnot engine exists between every pair of reservoirs at different temperatures.
Formalization targets
Goal: the thermodynamic temperature scale
Assume C is a convex cone satisfying the Clausius statement and admitting Carnot engines, and fix a reference reservoir T0. Then there is exactly one T:Θ→R with
T>0,T(T0)=273.16,∣e(h)∣∣e(k)∣=T(h)T(k)for every Carnot engine e between k<h.
The content is that the left-hand ratio, a property of an engine, depends on the engine only through its two reservoirs, and does so as a ratio of numbers attached to the reservoirs. The normalization 273.16 at T0 is the Kelvin convention; as the article notes, "any reference temperature and any positive numerical value could be used".
The two impossibility statements
Clausius(C)⟺Kelvin(C)
for a convex cone admitting Carnot engines over at least two reservoir temperatures.
Carnot's theorem
For k<h and a Carnot engine e between them, every realizable engine c between the same two reservoirs satisfies η(c)≤η(e), with strict inequality when c is irreversible (−c∈/C); and any two Carnot engines between h and k have equal efficiency.
Significance
The goal makes explicit which features of the realizable processes are used when a textbook defines absolute temperature. The three hypotheses — convex cone, Clausius statement, existence of Carnot engines — are exactly what the derivation needs; the working substance, the equation of state, and the mechanism of a device never appear. Uniqueness is the part that is usually left implicit: the Kelvin normalization at one reservoir pins the scale completely.
The milestones follow the article's own order: equation (1) for the efficiency of an arbitrary engine; the equivalence of the Clausius and Kelvin statements; the three clauses of Carnot's theorem; the composition law f(T1,T3)=f(T1,T2)f(T2,T3) for the reversible heat ratio across an intermediate reservoir; the conversion of that law into a scale; and the resulting Carnot efficiency η=1−TC/TH.
Difficulty
The mathematics is elementary — convex cones in a space of finitely supported real functions, and one functional equation — but the bookkeeping is where faithfulness is won or lost. Carnot's theorem is proved by scaling a candidate engine against a reversed Carnot engine so that the composite does no net work, which makes the Clausius statement applicable; the equality case has to identify the candidate as a positive multiple of the Carnot engine before concluding that it is reversible. The composition law requires the intermediate reservoir to be exactly balanced, so that the composite is again a two-reservoir engine — the article's remark that "the heat q2 flown from the engine 2 to the intermediate reservoir must be equal to the heat q2∗ flown out from the reservoir to the engine 3". The functional equation must be solved without assuming the reference reservoir is the hottest or the coldest, so the scale is defined by f(T0,⋅) below T0 and by 1/f(⋅,T0) above it, and the composition law is applied in four different orderings.
Formalization scope
Reservoir temperatures form an arbitrary linearly ordered type; the order is the empirical hotter/colder relation and nothing else is assumed of it.
Heat is positive into the device. Work is the sum of the heats, which is the first law for a cycle; the first law is therefore built into the model rather than stated.
A device is recorded only by its heat exchanges. Working substance, the states traversed, quasi-staticity and the internal construction do not appear; "reversible" means exactly that the reversed heat exchange pattern is also realizable.
The Clausius statement is formalized for a cycle with zero net work exchanging heat with exactly two reservoirs: that is the article's "without some other change, connected therewith, occurring at the same time".
Out of scope for this mission: the Clausius inequality ∮δQ/Tsurr≤0, the definition of entropy as a state function, Carathéodory's and Planck's formulations, the statistical-mechanical derivations, and the Clausius–Duhem inequality for continua.
Einstein Field Equations I: Local Energy–Momentum Conservation in D DimensionsTextbook
Motivation
The Einstein field equations (EFE) are the field equations of general relativity. They relate the curvature of spacetime, encoded in the metric gμν, to its matter–energy content, encoded in the stress–energy tensor Tμν:
Gμν+Λgμν=κTμν,κ=c48πG.
They were published by Einstein in 1915. A structural requirement on any such equation is that it be compatible with the local conservation of energy and momentum, ∇βTαβ=0. For the EFE this compatibility is automatic: the left-hand side is identically divergence-free by the contracted Bianchi identity, so the field equations force conservation of the right-hand side. This mission formalizes that implication together with the standard algebraic reformulations of the equations, following the Wikipedia article Einstein field equations (revision 1374672993).
Setting
Everything is written in one coordinate chart. A point is x∈RD, where D is the spacetime dimension (not fixed to 4). A 2-tensor field is a matrix-valued function x↦Tab(x). A metric on an open set U is a 2-tensor field g whose components are C∞ on U and which is symmetric and nondegenerate at each point of U; gab denotes the inverse matrix. Partial derivatives ∂a are coordinate (Fréchet) derivatives. Using the Misner–Thorne–Wheeler (+++) convention adopted by the source:
Ricci tensorRμν=Rαμαν, scalar curvatureR=gμνRμν, Einstein tensorGμν=Rμν−21Rgμν;
covariant derivative∇λTμν=∂λTμν−ΓλμσTσν−ΓλνσTμσ and divergence(divT)μ=gνλ∇λTμν, the index-lowered form of ∇νTμν;
the field equations on U: Gμν(x)+Λgμν(x)=κTμν(x) for all x∈U.
Formalization targets
Goal: local energy–momentum conservation
For a smooth nondegenerate symmetric metric g on an open U⊆RD, constants G,c>0, κ=8πG/c4, any Λ, and any T such that the field equations hold on U:
gβλ∇λTαβ=0on U.
Milestones
Symmetry of the Einstein tensor, Gμν=Gνμ.
Metric compatibility, ∇λgμν=0.
The differential Bianchi identity Rμα[βγ;ε]=0.
The contracted Bianchi identity, Gαβ;β=0.
The trace of the EFE, R−2DR+DΛ=κT.
Equivalence with the trace-reversed form Rμν−D−22Λgμν=κ(Tμν−D−21Tgμν) for D=2.
The vacuum equations Rμν=Λgμν/(2D−1) when T=0, D=2.
The cosmological term as vacuum stress–energy Tμν(vac)=−κΛgμν, with ρvac=−pvac=Λ/κ.
Minkowski space diag(−1,1,1,1) is a vacuum solution.
Significance
Conservation of Tμν is a consistency condition that any relativistic field equation for gravity has to satisfy; it is also the starting point for deriving the motion of matter (for dust, the geodesic equation) from the field equations alone. The trace-reversed and vacuum forms are the ones used in practice to study vacuum solutions such as Schwarzschild and Kerr, and Einstein manifolds in geometry.
All of these results are classical and fully proved in the literature. What is not yet available is a machine-checked development of the Levi-Civita curvature calculus in coordinates suitable for these statements: Mathlib currently has no Riemann or Ricci tensor of a pseudo-Riemannian metric. The mission produces that layer and the proofs on top of it.
Difficulty
The algebraic milestones (trace, trace reversal, vacuum form, vacuum energy) only need matrix algebra and the symmetry of g; without symmetry the identity gμνgμν=D fails, so symmetry is essential. The analytic core is the contracted Bianchi identity: it needs the symmetry of mixed second and third partial derivatives (Schwarz's theorem on an open set), derivatives of the matrix inverse g−1, and lengthy index bookkeeping with sums over FinD. Since all objects are total functions, every derivative step must be justified by differentiability on the open set U, not taken for granted.
Formalization scope
Coordinates: RD is represented as Fin D → ℝ; tensors are Matrix (Fin D) (Fin D) ℝ-valued functions; no manifold structure. D is arbitrary, and no Lorentzian signature is imposed (the results hold for any nondegenerate symmetric metric, which includes the Lorentzian case).
Regularity: the metric components are C∞ on the open set U; T carries no assumptions, as the field equations determine it on U.
Total-function conventions: the inverse of a singular matrix is 0 and the derivative of a non-differentiable function is 0; nondegeneracy on U and openness of U make these junk values irrelevant at points of U.
Divergence is taken on the second index; for symmetric tensors such as G this is the same as the first.
Not formalized: the Newtonian limit (it rests on approximations), the Einstein–Maxwell system, and the polynomial form.
The electroweak interaction is the unified description of electromagnetism and the weak nuclear force in the Standard Model of particle physics. Glashow, Salam and Weinberg shared the 1979 Nobel Prize in Physics for this unification, and its central predictions — neutral weak currents (Gargamelle, 1973) and the massive W± and Z0 bosons (UA1/UA2, 1983) — have been confirmed experimentally. The algebraic core of the theory is short: an SU(2)×U(1) gauge theory whose symmetry is broken by the vacuum value of a Higgs doublet, after which the neutral gauge fields are rotated by the weak mixing angleθW into the photon and the Z0, and the masses satisfy mZ=mW/cosθW.
This mission formalizes those tree-level algebraic facts as they are stated in the Wikipedia article Electroweak interaction (revision 1360331872), sections Formulation and Lagrangian.
Timeline (as recounted in the article). After the Wu experiment (1956) exposed parity violation, Glashow combined a chiral and an achiral symmetry and predicted a Z boson; Salam and Ward (1964) proposed a massless photon and three massive gauge bosons with manually broken symmetry; Weinberg (1967) used spontaneous symmetry breaking (Englert–Brout, Higgs, Guralnik–Hagen–Kibble, 1964) and predicted rough W and Z masses; 't Hooft (1971) proved that spontaneously broken gauge theories are renormalizable. Neutral currents were observed in 1973 and the W, Z in 1983.
Setting
Fix the weak isospin couplingg>0 and the weak hypercharge couplingg′>0. The gauge fields at a spacetime point, for one Lorentz component, are real numbers W1,W2,W3 (weak isospin) and B (weak hypercharge).
Weak mixing angle.θW=arctan(g′/g), the angle of the right triangle with legs g, g′ and hypotenuse g2+g′2 (the figure Weinberg's weak mixing angle in the article). The electromagnetic coupling is the altitude of that triangle, e=gg′/g2+g′2.
Physical neutral fields.γ=cosθB+sinθW3 and Z0=−sinθB+cosθW3.
Charged fields.W±=(W1∓iW2)/2.
Higgs sector. The Higgs field is a complex doublet h∈C2 with ∣h∣2=∣h1∣2+∣h2∣2 and potential V(h)=λ(∣h∣2−v2/2)2. Its vacuum is h0=(0,v/2).
Covariant derivative.Dμ=∂μ−i2g′YBμ−i2gσjWμj with σj the Pauli matrices. On the constant vacuum the derivative term vanishes, so ∣Dμh0∣2=∣Mh0∣2 with M=2g′YB1+2g(W1σ1+W2σ2+W3σ3) and Y=1 for the Higgs.
Charge.Q=T3+21Y.
Masses.mW=gv/2 and mZ=g2+g′2v/2.
Formalization targets
Goal: the gauge boson mass terms
For all real v,W1,W2,W3,B,
∣Mh0∣2=mW2W+W−+21mZ2(Z0)2,
which is the mass part mW2Wμ+W−μ+21mZ2ZμZμ of LK in the article. The photon field γ does not appear: it stays massless.
Milestones
The Weinberg triangle: cosθW=g/g2+g′2, sinθW=g′/g2+g′2.
e=gsinθW=g′cosθW.
(γ,Z0) is a rotation of (B,W3): it can be inverted and γ2+(Z0)2=B2+W32.
h0 minimizes the Higgs potential, and the zero set of V is the sphere ∣h∣2=v2/2.
Q=T3+21Y annihilates h0, and (for v=0) it is the only combination aT3+bY/2, up to scale, that does.
The neutral-current coupling: gT3W3+g′2YB=eQγ+cosθWg(T3−sin2θWQ)Z0.
Significance
These identities are the content behind the article's statements that "the axes ... have been rotated by the angle θW", that mZ=mW/cosθW, that e=gsinθW=g′cosθW, and that U(1)em is unbroken because it does not couple to the Higgs. The formal library produced here (Pauli matrices, the Higgs doublet, the mixing rotation, the mass formulas) can be reused by later missions on the full electroweak Lagrangian, on the Yukawa sector, or on the custodial relation ρ=mW2/(mZ2cos2θW)=1.
All statements are classical tree-level results. None is open; the work is to formalize them.
Difficulty
The mathematics is finite-dimensional linear algebra and trigonometry. The main pitfalls are conventions: the factor 21 in front of the Pauli matrices, the hypercharge normalization Q=T3+Y/2 (so the Higgs has Y=1), and the sign conventions in the rotation matrix. The goal also requires handling complex doublets, 2 factors and the inverse trigonometric definition of θW at the same time.
Formalization scope
All fields are real numbers at one spacetime point and for one Lorentz component; the Lorentz contraction and the kinetic (derivative) terms are not modeled. The Higgs doublet is a function Fin 2 → ℂ, and its squared norm is written out as ∣h1∣2+∣h2∣2 (not the sup norm of the product type). The Pauli matrices are defined explicitly. The angle is θW=arctan(g′/g); with g>0 this agrees with the triangle. The standing convention g,g′>0 is a hypothesis wherever θW appears. The goal is stated as an equation in C so that the product W+W− appears exactly as in LK.
All declarations live in the namespace ElectroweakWiki.
Wheeler–DeWitt Equation: Superspace, Semiclassical Limit and Lattice Constraint AlgebraResearch Paper
Motivation
The Wheeler–DeWitt equation is the quantum version of the Hamiltonian constraint of general relativity. In the canonical formulation, a spacetime is described by the evolution of the metric hab of a spatial slice, and Einstein's equations include a constraint: the Hamiltonian density vanishes at every point of space. Quantizing this constraint, as DeWitt did (Phys. Rev. 160, 1967), gives a second-order functional differential equation for a wavefunctional Ψ[h] with no external time parameter.
Three structural facts make the equation distinctive:
its kinetic term is built from an indefinite metric on the space of metrics, the DeWitt supermetric, so the equation is hyperbolic;
the timelike direction of that supermetric is the local volume, which serves as an intrinsic time;
in the limit ℏ→0, a WKB ansatz turns the equation into the Einstein–Hamilton–Jacobi equation, recovering classical general relativity.
A fourth issue is not settled: consistency. In the continuum, the Poisson bracket of two Hamiltonian constraints is a combination of the momentum constraints, so imposing all constraints is consistent (Dirac's first-class condition). Any regularization of the quantum equation, a lattice for instance, must reproduce this. Generic lattice discretizations do not: the constraints become second class (Bander 1987; Dittrich 2008; Bahr–Dittrich 2009). Whether, and with which anomaly, the Hamiltonian-constraint commutators close as the lattice is refined was singled out as an open question by Loll (Living Rev. Rel. 1998, §2.12).
The first three facts are the standard entry points of quantum geometrodynamics, as reviewed by Kiefer (arXiv 0812.0295): Section 2 for the constraints and their quantization, eqs. (3)–(8), and Section 3 for the Hamilton–Jacobi equation, eq. (10), and the semiclassical approximation, eq. (11).
Setting
A spatial metric at a point is a symmetric positive-definite 3×3 matrix h=(hab). Write h=deth. The DeWitt supermetric is
Gabcd(h)=2hhachbd+hadhbc−habhcd.
It is a quadratic form G(p,p)=Gabcdpabpcd on symmetric momenta. Its inverse Gabcd=21h(hachbd+hadhbc−2habhcd) acts on velocities.
Space is represented by a finite set X of sites; this is a lattice regularization. A configuration assigns to each site x a matrix h(x). All nine entries are coordinates, so configuration space is R9∣X∣, and physical configurations are those that are symmetric and positive definite at every site. A wavefunction is a function Ψ from configurations to C. The functional derivative δ/δhab(x) is the partial derivative ∂ab(x) with respect to the coordinate hab(x), with the cell volume normalized to one.
The Wheeler–DeWitt operator at site x, with gravitational coupling κ, cosmological constant Λ, and a supplied potential R(h,x) playing the role of the spatial scalar curvature, is
Write Vz(h)=−2κh(z)(R(h,z)−2Λ) for the potential at site z. On the classical phase space of pairs (h,p), with pab(z) conjugate to hab(z), the lattice constraints are
Hz(h,p)=2κGabcd(h(z))pab(z)pcd(z)+Vz(h),
with the canonical Poisson bracket. They are first class if each bracket {Hx,Hy} vanishes wherever all constraints vanish.
The supermetric stands to the left of both derivatives (Kiefer, eq. (7), vacuum case; the coupling is κ=8πG, so 2κ=16πG as in Kiefer's eq. (10)). The kinetic term involves only derivatives at the same site, as in the continuum; it is ultralocal. The potential may couple neighbouring sites.
Formalization targets
Milestones: geometry of the supermetric
Signature.G(p,p)>0 for every nonzero symmetric p with tr(hp)=0, and G(h−1,h−1)<0.
Inverse.GabcdGcdefpef=pab for symmetric p.
Intrinsic time. With h=sh^, deth^=1 and τ=32/3s3/4, the velocity form of the supermetric is −dτ2+323τ2tr((h^−1dh^)2).
The ℏ-independent part is the Einstein–Hamilton–Jacobi functional of Kiefer's eq. (10), Hamiltonian part, with S the eikonal of the WKB form (11) (amplitude C=1). This is the weakest statement exhibiting the classical limit, and it holds for every potential.
If every Vz depends only on h(z), then {Hx,Hy}=0 identically.
Goal: lattice constraints are second class
If the symmetric part of ∂Vy/∂hab(x) is nonzero at a physical h, for some x=y, then there is a point of the constraint surface over h at which {Hx,Hy}=0.
Every finite-difference curvature couples neighbouring sites, so the goal applies to every such discretization. It is the weakest statement that the discretized constraint algebra is inconsistent, and it holds uniformly in the discretization.
Stronger targets
The limit form ℏ→0, as a Tendsto statement.
The same results with the Laplace–Beltrami ordering of the kinetic term on the six-dimensional fibre of symmetric matrices, and an explicit formula for the difference between the two orderings.
A potential realized as the scalar curvature of a smooth metric field on a coordinate patch, computed from Christoffel symbols, with its scaling law R[sh]=R[h]/s for constant s>0.
The momentum constraints (Kiefer, eq. (8)).
Open as far as the drafter could determine. Fix a cubic lattice of spacing ε and a standard finite-difference scalar curvature. Determine the limit ε→0 of the suitably rescaled commutator [Hx,Hy], and decide whether it reproduces the continuum Dirac algebra {H(x),H(y)}=(habHb(x)+habHb(y))∂aδ(x,y) or leaves an anomaly. Give the anomaly explicitly if there is one. This is the metric-variable analogue of the question Loll raised.
Significance
The signature theorem is the mathematical content behind two physical facts: the hyperbolic character of the Wheeler–DeWitt equation, and the instability of the conformal mode (the "conformal factor problem" of Euclidean quantum gravity). DeWitt's cone form identifies the local volume as an intrinsic clock. On this rests the reading of the equation as a wave equation in superspace, with τ as time. The conformal-sector reduction shows explicitly how the timelike direction appears in the operator, with the sign opposite to the traceless modes. The WKB identity is the precise sense in which the quantum constraint contains classical general relativity: its leading order is the Hamilton–Jacobi form of the Einstein constraint (Peres 1962; Kiefer, eq. (10)), and its next order is the one used in derivations of the Schrödinger equation from quantum gravity.
All of these are classical computations. The drafter knows of no machine-checked version. A formalization fixes conventions that the literature handles loosely:
the normalization of Gabcd;
the treatment of symmetric versus independent metric components;
the ordering of operators.
It does so with genuine derivatives in place of formal symbols.
The constraint-algebra part turns a phenomenon described qualitatively in the literature into a theorem with exact hypotheses. Lattice constraints are second class as soon as curvature couples sites, and the quantum commutator is given in closed form. It also isolates precisely what a consistent discretization would have to cancel. The continuum-limit target addresses an open question in canonical quantum gravity.
Difficulty
The supermetric statements are finite-dimensional, but none is a single computation. The signature requires a positivity argument for tr((h1/2ph1/2)2) on an h-traceless subspace. The cone form requires separating the trace part of a velocity in a non-orthonormal frame.
The analytic targets require second derivatives of compositions through the determinant and the matrix inverse. The derivative of det at h in the direction of a single entry is deth(h−1)ba, and the derivative of h−1 produces terms like h−1Eabh−1; a sign or index slip changes the constant 7. Treating the nine entries as independent coordinates, rather than the six symmetric ones, changes the numerical coefficients; the conventions below fix this.
For the constraint algebra, the quantum commutator involves up to fourth derivatives of Ψ, which cancel only by the symmetry of mixed partials. The classical goal needs a point on the constraint surface where the bracket is nonzero. Each constraint fixes the value of an indefinite quadratic form in the momenta, and the bracket is linear in them, so the proof must use the Lorentzian signature to show that the bracket is not constant on that quadric.
A naive formalization that replaces ∂ab by an abstract linear map reduces every target to ring arithmetic, and is ruled out by the definitions.
Formalization scope
Definitions (module Def_WheelerDeWittSuperspace):
metric components as Fin 3 → Fin 3 → ℝ, so that configuration space X → Fin 3 → Fin 3 → ℝ with X a Fintype is a normed space;
matrices via Matrix.of;
∂ab(x) as fderiv ℝ Ψ h applied to the coordinate direction of hab(x);
h as Real.sqrt of the determinant;
the operator with the supermetric on the left.
Conventions:
all nine metric entries are independent coordinates, and statements that need symmetry assume it at the evaluation point (Matrix.PosDef includes symmetry);
lattice cell volume equal to 1;
vacuum (no matter term).
Not formalized; the potential is an interface: the potential R is supplied data on configurations and sites, and every target holds for all such R. The continuum limit of the lattice, a Hilbert-space domain for H, regularization of coincident derivatives in the continuum, and existence of solutions are outside the mission and must not be inferred from it.
Trivial cases to rule out:
the WKB milestone must hold for every C2 phase S, not only for Ψ=0 or constant S;
the goal must exhibit a point on the full constraint surface, with every Hz vanishing, not merely a point where the bracket is nonzero;
the dependence hypothesis is on the symmetric part, because antisymmetric metric directions are invisible to G.
Phase-space convention: the bracket is the canonical one on all nine components per site. The symmetric reduction changes off-diagonal normalizations but not the first-class property.
Welcome contributions: the stronger targets, a matter field with its energy density, and minisuperspace reductions (e.g. homogeneous isotropic metrics h=a2δ) derived from the operator here rather than postulated.
R. Loll, Discrete approaches to quantum gravity in four dimensions, Living Rev. Relativ. 1 (1998), 13, Section 2.12. https://arxiv.org/abs/gr-qc/9805049
R. Loll, On the diffeomorphism-commutators of lattice quantum gravity, Class. Quantum Grav. 15 (1998), 799–809. https://arxiv.org/abs/gr-qc/9708025
P. Renteln, Some results of SU(2) spinorial lattice gravity, Class. Quantum Grav. 7 (1990), 493–502.
Naculich 2020: Regge colour basis of the four-gluon amplitude and the all-loop coefficient B_ℓℓResearch Paper
Motivation
At high energies (s≫−t, the Regge limit) scattering amplitudes of gauge theories simplify and organise into towers of logarithms log∣s/t∣. For the four-gluon amplitude of N=4 super-Yang–Mills theory with gauge group SU(N), the colour structure is what makes the nonplanar (subleading in 1/N) amplitude hard to control. Naculich (2020) builds a basis of colour factors Cik adapted to the Regge limit at every loop order, writes the ℓ-loop amplitude in that basis, and uses one of its coefficients, Bℓℓ(ℓ), to motivate a conjectured all-loop relation with the Regge limit of the N=8 supergravity four-graviton amplitude. The physics conjecture itself is not a mathematical statement; what is exact and checkable is the colour algebra underlying it — finite-dimensional linear algebra over polynomials in N.
Setting
Trace basis. Four-gluon colour factors are expanded in six single and double traces c[1],…,c[6] of fundamental generators (eq. (2.1)). A colour vector is a column v=(v1,…,v6) of polynomials in N, with complex coefficients.
Colour operators. The t-channel Casimir Tt2 and the crossing-odd combination Ts−u2=21(Ts2−Tu2) act on colour vectors as the explicit 6×6 matrices of eq. (4.15) (entries in Z[21][N]). The tree-level t-channel colour factor is C00=c[1]−c[3]. Crossing (exchange of legs 2 and 3) is the permutation P: c[1]↔c[3], c[4]↔c[6].
each with exactly i operators. The admissible pairs are (0,0), (1,1), (2,1),(2,2) and (i,1),(i,i−1),(i,i) for i≥3.
Extended trace basis. At ℓ loops, tλ+6k(ℓ)=Nℓ−2kc[λ] for λ=1,2,3 and Nℓ−2k−1c[λ] for λ=4,5,6, giving 3ℓ+3 elements (eq. (2.3)). The colour-ordered amplitudesAλ(ℓ) of an amplitude are its coordinates in this basis (eq. (2.4)).
Regge-basis amplitude. Eq. (4.24):
A(ℓ)=A1(0)a~ℓi=0∑ℓk∑Bik(ℓ)Nℓ−iCik.
Null vectors. Four vectors per loop order ℓ≥2 (eqs. (2.10), (2.12), extended by prepending zeros) encode the SU(N) group-theory relations among the Aλ(ℓ).
Formalization targets
Goal: the all-loop coefficient Bℓℓ(ℓ) (eqs. (4.26)–(4.27))
For every even ℓ,
A1(0)a~ℓBℓℓ(ℓ)=2⋅3ℓ/21(A3ℓ+1(ℓ)−A3ℓ+3(ℓ)),
and for ℓ=1, A1(0)a~B11(1)=−(A1(1)+A3(1)). The paper obtains the even-ℓ formula by guessing from explicit results through ℓ=9; a proof for all ℓ is not given there.
The Regge colour factors satisfy the four null-vector constraints and have no t2(ℓ) component, for all ℓ≥2.
The Regge colour factors are linearly independent, and for ℓ≥2 span all constrained Regge-limit amplitudes (the "basis" claim of Sec. 4).
Significance
The result. Eq. (4.26) is the input to the paper's analysis of the infrared behaviour of Bℓℓ(ℓ) (at most 1/ϵℓ, manifest for even ℓ because only most-subleading-colour amplitudes appear), and hence to the conjectured all-orders relation between Bℓℓ(ℓ) and the Regge limit of N=8 supergravity (Sec. 6). An all-ℓ proof replaces an extrapolation from nine loop orders by a theorem.
Formalizing it. The low-loop tables (4.16), (4.20), (4.18), (4.22) are finite computations; the goal and milestones 3, 4, 7, 8 are statements for every loop order and need genuine arguments (induction on ℓ, structure of the operators on the extended basis).
Difficulty
The coefficient Bℓℓ(ℓ) is defined through the whole (3ℓ−2)-element basis: a formula for it must vanish on every other Nℓ−iCik, including the nested commutators Ci1 and Ci,i−1 whose components grow in complexity with i. Checking the formula loop order by loop order does not give a proof; one needs to control how Tt2 and Ts−u2 act on the two most-subleading slots of the extended basis uniformly in ℓ.
Formalization scope
Colour vectors are Fin 6 → ℂ[X], with N the polynomial variable X; operators are Matrix (Fin 6) (Fin 6) ℂ[X] acting by Matrix.mulVec. Coefficients Bik and the prefactor κ=A1(0)a~ℓ are arbitrary complex numbers; B is a function ℕ × ℕ → ℂ of which only the admissible pairs matter.
Extended-basis indices λ are 1-based natural numbers, 1≤λ≤3ℓ+3; the coordinate function returns 0 outside this range. The colour-ordered amplitudes are defined as the coordinates of κ∑BikNℓ−iCik, so the goal says: whatever the coefficients, the stated combination of coordinates returns κBℓℓ. This is equivalent to the paper's reading for amplitudes lying in the span of the Regge basis (milestone 8).
The odd-ℓ formula (4.28) for ℓ>1 is not included. It agrees with the model at ℓ=3 (it is the B33 line of eq. (4.22)), but an exploratory exact-rational computation (not machine-checked) found that at ℓ=5 no linear combination of the six amplitudes A3ℓ−8,A3ℓ−6,A3ℓ−5,A3ℓ−3,A3ℓ−2,A3ℓ reproduces B55(5) in this model. Including it would risk a false milestone.
Local checks of the statements used Lean v4.28.0 with Mathlib; only standard Mathlib API (polynomials, matrices, LinearIndependent, finite sums) is used.
Welcome contributions: lemmas describing Tt2 and Ts−u2 as operators on extended-basis coordinates (shift of N-degree), and general tools for iterated commutators.
Selected references
S. G. Naculich, All-loop-orders relation between Regge limits of N=4 SYM and N=8 supergravity four-point amplitudes, arXiv:2012.00030 (2020), https://arxiv.org/abs/2012.00030.
Witten 1998: Anti de Sitter Space and Holography — the massive scalar boundary-value problemResearch Paper
Motivation
In February 1998 Witten's paper Anti de Sitter Space and Holography made precise Maldacena's proposal that large-N conformal field theories in d dimensions are described by gravity on (d+1)-dimensional anti de Sitter spaceAdSd+1. The dictionary it proposes identifies boundary values of bulk fields with sources of boundary operators, and reads correlation functions off the dependence of the bulk action on those boundary values. Its most quoted quantitative consequence is the mass–dimension relation: a scalar of mass m on AdSd+1 corresponds to a boundary operator of conformal dimension
Δ=21(d+d2+4m2),Δ(Δ−d)=m2.
Behind this physics statement sits a piece of classical analysis on hyperbolic space: the Dirichlet-type problem for the operator Δg−m2 with prescribed boundary asymptotics, solved by an explicit Poisson-type integral (Section 2 of the paper, eqs. (2.15)–(2.44)). This mission formalizes that analytic core.
Timeline. Maldacena 1997 (hep-th/9711200) conjectured the correspondence; Gubser–Klebanov–Polyakov (hep-th/9802109) and Witten (hep-th/9802150, both February 1998) proposed the precise dictionary; the Euclidean black-hole thermodynamics used in Section 3 goes back to Hawking–Page 1983.
Setting
Euclidean AdSd+1 is modelled as the upper half space{(x0,x):x0>0,x∈Rd} with the hyperbolic metric
A massive scalar of squared mass m2∈R solves Δgu=m2u (eqs. (2.34), (2.37)). The conformal dimension is Δ=21(d+d2+4m2) (eq. (2.44)), the bulk-to-boundary kernel is
and the Poisson-type integral of a boundary datum φ0 is ∫RdKΔ(x0,x−x′)φ0(x′)dx′ (eqs. (2.20), (2.41)).
Formalization targets
Goal (eq. (2.41) with (2.43)–(2.44))
For m2>−d2/4 there is a constant c>0 such that for every bounded continuous φ0 on Rd,
φ(x0,x)=c∫Rd(x02+∣x−x′∣2)Δx0Δφ0(x′)dx′
is a C2 solution of Δgφ=m2φ on the upper half space and x0Δ−dφ(x0,x)→φ0(x) as x0→0+, i.e. φ∼x0d−Δφ0 at the boundary (eq. (2.36) with f=x0).
Milestones
x0Δ solves the massive equation when Δ(Δ−d)=m2 (eqs. (2.16)–(2.17), (2.37)).
The inverted kernel KΔ solves it too (eqs. (2.18)–(2.19), (2.38)).
Scaling: ∫x02Δ−d(x02+∣x∣2)−Δdx is independent of x0 (eq. (2.40)).
Approximate identity: these kernels tend to a multiple of δ(x) (eq. (2.39)).
The massless case Δ=d with φ→φ0 (eq. (2.20)).
Δ=21(d+d2+4m2) is the larger root of Δ(Δ−d)=m2 (eqs. (2.43)–(2.44)).
The gauge-field profile x0d−2 solves the reduced Maxwell equation (eqs. (2.24)–(2.25)).
Hawking–Page, d=3: β(r+)=12πr+/(1+3r+2) has a maximum and I(r+)=πr+2(1−r+2)/(1+3r+2) is negative for large r+ (eqs. (3.15)–(3.16)).
Significance
The goal is the analytic statement underlying the mass–dimension relation Δ(Δ−d)=m2 of the AdS/CFT dictionary: it shows that the boundary datum φ0 of a massive field is naturally a density of weight d−Δ, which is how the paper identifies the dimension of the dual operator. The same kernel is the starting point for the two-point function (2.42). The result is classical mathematics (scattering theory on hyperbolic space), proved informally in the paper; this mission produces a machine-checked version together with reusable infrastructure for the hyperbolic Laplacian in half-space coordinates.
Difficulty
The pointwise identities (milestones 1, 2, 6, 7, 8) are calculus. The substance is analytic: differentiating the Poisson-type integral twice under the integral sign with only boundedness of φ0, and proving the approximate-identity limit with the right normalisation; the scaling integral requires Δ>d/2, and the endpoint m2=−d2/4 (where Δ=d/2 and the kernel is not integrable) must be excluded.
Formalization scope
Points are pairs (x0,x)∈R×EuclideanSpaceR(Find); partial derivatives are one-variable derivs, so the solution predicate also demands C2 regularity on {x0>0}, which rules out junk values of deriv. Real powers are Real.rpow; integrals are Bochner integrals for Lebesgue measure. Boundary data are bounded continuous functions. All declarations live in the namespace WittenAdSHolography, with the shared definitions in one definition file.
J. Maldacena, The Large N Limit of Superconformal Field Theories and Supergravity, Adv. Theor. Math. Phys. 2 (1998) 231–252. https://arxiv.org/abs/hep-th/9711200
S. S. Gubser, I. R. Klebanov, A. M. Polyakov, Gauge Theory Correlators from Non-Critical String Theory, Phys. Lett. B 428 (1998) 105–114. https://arxiv.org/abs/hep-th/9802109
S. W. Hawking, D. N. Page, Thermodynamics of Black Holes in Anti-de Sitter Space, Commun. Math. Phys. 87 (1983) 577–588. https://doi.org/10.1007/BF01208266
Coupling Constant II: The Beta Function in the Energy Scale and the Running DichotomyTextbook
Motivation
A coupling constant is the number that fixes the strength of an interaction: in the
modern formulation it is the parameter weighting the interaction part of a Lagrangian
against its kinetic part. In a quantum field theory such a parameter is not constant.
Probing a process at a higher energy resolves more of the virtual activity that dresses
the interaction, and the value one measures depends on the energy scale μ at which
the measurement is made. This scale dependence is called the running of the coupling,
and it is the organising phenomenon behind two of the best-known qualitative facts about
the Standard Model: the electromagnetic coupling grows with energy (α≈1/137
at low energy, α≈1/127 at the Z mass, about 90 GeV) and apparently
diverges at a finite energy — the Landau pole — while the strong coupling decreases
logarithmically at high energy, the asymptotic freedom found by Gross, Politzer and
Wilczek and recognised by the 2004 Nobel Prize.
The source treats the running through the beta function, defined there by
β(g)=μ∂μ∂g=∂lnμ∂g,
and through the one-loop strong coupling
αs(k2)≈β0ln(k2/Λ2)1.
Everything in this mission is the mathematics of those two formulas: an autonomous scalar
ordinary differential equation in the logarithm of the energy, and the behaviour of its
solutions according to the sign of the one-loop coefficient.
Setting
Write μ>0 for the energy scale of a process and g, or α, for a real-valued
coupling as a function of that scale. The beta function in the energy scale is the
logarithmic derivative used by the source,
β(g)(μ)=μ⋅dμdg(μ),
a function of μ built from g; a theory whose beta function vanishes identically is
scale-invariant, since its coupling then takes the same value at every scale.
A coupling is said to be one-loop running with coefficient b on a set S of scales
when it satisfies, at every μ∈S,
μdμdα(μ)=bα(μ)2,
equivalently dα/dμ=bα2/μ. Here b<0 is the asymptotically free sign
(QCD) and b>0 the sign the source attributes to QED. The equation is the one-loop
truncation written in μ rather than in t=lnμ: both forms appear in the source,
and the first milestone is precisely their agreement.
The one-loop strong coupling with one-loop coefficient β0>0 and QCD scaleΛ>0 is the explicit function
α(Q)=β0ln(Q2/Λ2)1,Q>Λ,
the formula the source quotes for αs(k2). It is defined by an existing platform
definition, which this mission reuses rather than restates.
Formalization targets
Goal — the running dichotomy in the energy scale
For every reference scale μ0>0 and every positive reference value α0 of the
coupling at that scale:
b<0:α(μ)=1−bα0ln(μ/μ0)α0(μ≥μ0),α(μ)μ→∞0;b>0:no coupling obeys μdμdα=bα2 on all of [μ0,μ0e1/(bα0)] with α(μ0)=α0.
The negative-coefficient half is asymptotic freedom together with the uniqueness of the
solution; the positive-coefficient half is the Landau pole, stated as the non-existence of
any solution up to the pole scale rather than as a divergence of one chosen formula. The
pole scale μ0e1/(bα0) is the one the one-loop equation itself predicts.
Supporting levels
The milestones below range from the definitional identity
μ∂μg=∂lnμg, through the three qualitative readings of the
sign of β that the source states in words (vanishing beta function means scale
invariance; positive beta function means the coupling grows with energy; negative means it
falls), to the identification of the quoted QCD formula as an exact solution of the
one-loop equation and the determination of Λ from a single measured value.
Significance
The dichotomy is the mathematical content of two physical statements that the source
presents side by side. Asymptotic freedom is what makes perturbative QCD predictive at high
energy; the Landau pole is the reason the source warns that the perturbative beta function
cannot be trusted at strong coupling, since the one-loop equation has no solution reaching
the pole scale at all. Formalizing the uniqueness half also pins down what "the" running
coupling means: the explicit formula is not merely one solution of the renormalization-group
equation but the only one through a given data point.
An earlier public mission on the same source
("Coupling Constant: One-Loop Running, Asymptotic Freedom and the Landau Pole") formalized
the one-loop flow in the logarithmic variable t=lnμ with a cubic beta function
β(g)=bg3, and its statements are proved. This mission is stated in the other
variable the source uses — the energy scale μ itself, with the quadratic one-loop
equation for α — and adds the parts of the source that the earlier mission does not
cover: the defining identity between the two forms of the beta function, the sign-to-monotonicity
statements for an arbitrary (not necessarily one-loop) coupling, the explicit pole scale, and
the determination of the QCD scale Λ from a measurement. Where the earlier mission
already settles a statement in the quoted formula's variables, this proposal references the
published theorem instead of restating it.
Difficulty
The equation is elementary, and the difficulty is entirely in what the statements quantify
over. The negative-coefficient half asserts that every function satisfying the equation on
[μ0,∞) is the explicit one, so it needs a uniqueness argument for a nonlinear
scalar ordinary differential equation rather than a verification that a formula works; the
substitution u=1/α that makes the equation linear is legitimate only while α
is known not to vanish, which is part of what has to be proved. The positive-coefficient half
is a non-existence statement on a closed interval reaching the pole, so no local solution can
be extended to it; the naive attempt — differentiate the explicit formula and observe that it
blows up — proves nothing about the other candidate solutions.
Formalization scope
Couplings are total functions R→R, and every statement restricts to
the scales it is about: the positive half-line for the beta-function statements, (Λ,∞)
for the explicit QCD formula, [μ0,∞) or a closed interval for the one-loop flow.
Outside those ranges nothing is claimed — in particular at Q≤Λ, where the logarithm
in the quoted formula vanishes or is negative and the reciprocal returns a junk value.
Differentiability is a hypothesis where it is needed, written as a two-sided derivative at each
scale of the set, so that a statement cannot be satisfied by a one-sided artefact at an endpoint.
The one-loop property is phrased as an exact differential equation, never as an asymptotic
approximation, and the Landau-pole milestone is a non-existence claim rather than a divergence
claim, so it cannot be discharged by exhibiting a formula with a singularity.
A complete development needs the mean value theorem, Gronwall- or uniqueness-type arguments for
scalar ordinary differential equations, and the standard logarithm and exponential API; the
uniqueness lemma for u˙=bu2 and the reciprocal substitution are reusable beyond this
mission. Contributions of a general scalar-ODE uniqueness lemma, usable by both this mission
and the earlier one, are especially welcome.
A. Deur, S. J. Brodsky, G. F. de Téramond, The QCD running coupling, Progress in Particle and Nuclear Physics 90 (2016) 1-74. https://arxiv.org/abs/1604.08082
ATLAS Collaboration, A precise determination of the strong-coupling constant from the recoil of Z bosons with the ATLAS experiment at s=8 TeV, 2023. https://arxiv.org/abs/2309.12986
P. A. Zyla et al. (Particle Data Group), Review of Particle Physics - Chapter 9, Quantum Chromodynamics, PTEP 2020 (8). https://doi.org/10.1093/ptep/ptaa104
12 thms3 active usersReviewed
🏆Completed
Captain: Lucas
Polyakov Action: Symmetries and Classical Equivalence with Nambu–GotoTextbook
Motivation
The Polyakov action is the standard action for the worldsheet of a bosonic string. It was introduced by Deser and Zumino and independently by Brink, Di Vecchia and Howe in 1976, and is associated with Polyakov after he used it to quantize the string in 1981. Its value over the Nambu–Goto (area) action is that it is quadratic in the embedding fields, which makes quantization tractable, while being classically equivalent to the area action. This mission formalizes the classical facts collected in the Wikipedia article Polyakov action.
Setting
A string sweeps out a two-dimensional worldsheet with coordinates σ=(τ,σ). Its embedding in a D-dimensional target is Xμ(σ); the target carries a metric gμν(X), and the worldsheet carries an independent (auxiliary) metric hab(σ) with inverse hab and determinant h. With the induced metric Gab=gμν(X)∂aXμ∂bXν and string tension T, the Polyakov and Nambu–Goto actions are
SP=2T∫d2σ−hhabGab,SNG=T∫d2σ−detGab.
The stress–energy tensor is Tab=T(Gab−21habhcdGcd). Signature convention: timelike +, spacelike −.
Formalization targets
Goal: on-shell equivalence with Nambu–Goto
If the auxiliary metric satisfies its equation of motion Tab=0 (with hcdGcd≥0), then
SP=SNG.
Milestones
Poincaré invariance of the target: translations X→X+b (exact) and infinitesimal Lorentz transformations X→X+εωX (to O(ε2)), for a constant target metric.
Worldsheet diffeomorphism invariance and Weyl invariance h→Λh.
Tracelessness of the stress–energy tensor, Taa=0.
The variations δ−h=−21−hhabδhab and δS/δhab=2T−h(Gab−21habhcdGcd).
On shell, detG=41h(hcdGcd)2.
Conformal gauge: S=2T∫(X˙2−X′2), the constraints X˙⋅X′=0, X˙2+X′2=0, the wave equation □Xμ=0, and its light-cone solution X=X+(τ+σ)+X−(τ−σ).
Significance
These are the classical foundations of bosonic string theory: the symmetries (Poincaré, diffeomorphism, Weyl) are what make the conformal gauge reachable, and the equivalence with Nambu–Goto justifies using the quadratic action. All statements are textbook results; the work here is to formalize them with explicit hypotheses (differentiability, integrability, positivity of hcdGcd) that the physics presentation leaves implicit.
Difficulty
The algebraic statements are 2×2 matrix identities. The diffeomorphism invariance needs the change-of-variables formula; the infinitesimal Lorentz statement needs care with integrability; the wave-equation milestone needs the fundamental lemma of the calculus of variations; and the light-cone solution is d'Alembert's theorem.
Formalization scope
Worldsheet =R2 (coordinates Fin 2 → ℝ, index 0=τ), target =RD. Actions are Bochner integrals over a region U⊆R2. The inverse metric is Mathlib's matrix inverse. Quantization, Virasoro operators, the Liouville action and string boundary conditions are out of scope.
Selected references
Wikipedia, Polyakov action, revision 1292212311.
S. Deser, B. Zumino, A complete action for the spinning string, Phys. Lett. B 65 (1976) 369–373.
L. Brink, P. Di Vecchia, P. Howe, A locally supersymmetric and reparametrization invariant action for the spinning string, Phys. Lett. B 65 (1976) 471–474.
A. M. Polyakov, Quantum geometry of bosonic strings, Phys. Lett. B 103 (1981) 207–210.
Drezet 2023: GHZ Bell-Nonlocality in the Everett MultiverseResearch Paper
Motivation
Whether the many-worlds (Everett) interpretation of quantum mechanics is local in the sense of Bell is a recurring question in the foundations of physics. A number of authors have argued that Everett's theory escapes Bell-type nonlocality theorems because it drops the assumption that each measurement has a single outcome. A. Drezet's paper arXiv:2306.07794 (2023) argues against this claim by revisiting the Greenberger–Horne–Zeilinger (GHZ) nonlocality theorem, which, unlike Bell's inequality, uses only probability-one (perfect) correlations. Its central mathematical step, stated in Section II, is that no locally causal beable model can reproduce the quantum predictions for the GHZ state. The paper's argument is that this step makes no assumption about completeness or hidden variables, and so applies to Everett's theory as well.
Timeline:
1964–1976: J. S. Bell introduces Bell inequalities and, in "The theory of local beables" (1976), the notion of local causality phrased with beables λ (collected in Speakable and Unspeakable in Quantum Mechanics, 2004).
1989: Greenberger, Horne and Zeilinger give a nonlocality proof for three spin-½ particles that uses no inequality ("Going beyond Bell's theorem").
2015: L. Vaidman discusses the GHZ argument in the many-worlds setting ("Bell inequality and Many-Worlds interpretation").
2023: Drezet's paper, the source of this mission.
Setting
Three spin-½ particles are sent to three distant observers, Alice, Bob and Charlie. Each one measures the spin of their particle along a direction n^j∈{x^,y^} and records an outcome α,β,γ∈{+1,−1}. The particles are prepared in the GHZ state
∣ψ⟩=21(∣+z⟩1∣+z⟩2∣+z⟩3−∣−z⟩1∣−z⟩2∣−z⟩3),
and quantum mechanics predicts the joint outcome probabilities through the Born rule, P(α,β,γ∣n^1,n^2,n^3,ψ)=∣⟨αn^1,βn^2,γn^3∣ψ⟩∣2. Here ∣an^⟩ denotes the eigenvector of σn^ with eigenvalue a. In Lean this is bornProb n₁ n₂ n₃ α β γ.
A locally causal (Bell-local) model consists of
a space Λ of beables with a probability distribution ρ(λ) that does not depend on the settings (Eq. (18), statistical independence), and
local response probabilities Pj(α∣λ,n^j), where party j's outcome depends only on λ and its own setting (Eq. (17)).
In Lean the model is a LocalCausalModel Λ and its prediction is M.predict. The beables are arbitrary: λ may be a hidden variable, or the quantum state itself as in the Beltrametti–Bugajski representation, Eq. (15).
Formalization targets
Goal (Section II, Eqs. (10)–(27))
No locally causal model reproduces the quantum predictions on the four GHZ settings (x^,x^,x^), (x^,y^,y^), (y^,x^,y^), (y^,y^,x^):
Eqs. (2)–(5).σx(1)σx(2)σx(3)∣ψ⟩=−∣ψ⟩, and the three σxσyσy-type operators fix ∣ψ⟩.
Born probabilities (Section III, Eqs. (30)–(32)). For (x^,x^,x^) each outcome triple has probability 1/4 if αβγ=−1 and 0 otherwise. For the other three GHZ settings the same holds with αβγ=+1.
Eqs. (10)–(13). No assignment Aj,Bj∈{±1} satisfies all four GHZ constraints.
Eqs. (19)–(27). If a product of three probability distributions on {±1} vanishes off the set {αβγ=s}, then each factor is a point mass, and the three atoms multiply to s.
Section II, summary (ii). A locally causal model that reproduces the GHZ predictions is, for ρ-almost every λ, deterministic, with values satisfying (10)–(13).
Section III (B). Local causality implies non-signalling.
Eq. (29). Quantum non-signalling for the GHZ state, with Alice's marginal P1(±1∣x^,ψ)=21.
Appendix 3, Eqs. (39)–(45). Keeping Eq. (17) but letting the beable density depend on the settings (a "superdeterministic" model) does reproduce the GHZ predictions.
Significance
The goal theorem is the mathematical core of the paper's conclusion that any beable account of the GHZ correlations, including one where λ is just the quantum state, violates Bell's local causality. Milestones 6–8 formalize the paper's separation of notions of locality: non-signalling holds both in quantum mechanics and in every locally causal model, so it cannot tell them apart; and dropping statistical independence (Eq. (18)) is enough to reproduce the correlations. The GHZ theorem itself has been known since 1989. As far as the drafter is aware it has no machine-checked formalization on this platform, and this mission provides one, stated for general measurable beable spaces.
Difficulty
The combinatorial core (milestone 3) is a parity argument. The substance lies in the steps that connect it to probability: showing that perfect correlations force determinism almost surely when ρ is an arbitrary probability measure (null sets, measurability, integrals of nonnegative functions), and computing the Born probabilities exactly in Lean from explicit matrices and eigenvectors.
Formalization scope
Outcomes are elements of ℤˣ (exactly {+1,−1}). Settings form a two-element type Setting with values x and y. Qubit index 0 is ∣+z⟩ and index 1 is ∣−z⟩. Three-qubit operators are Kronecker products.
A LocalCausalModel Λ has a probability measure ρ on a measurable space Λ and measurable response functions P j λ n α with values in [0,1] that sum to one over α. The fixed state ψ is left implicit. Predictions are Bochner integrals, and all integrands are bounded and measurable.
The goal quantifies over all such models on all measurable spaces. Its hypothesis is the equality of all 4×8 predicted and Born probabilities. Because it uses the actual Born probabilities rather than an arbitrary list of constraints, it cannot hold vacuously.
Needed infrastructure: finite-dimensional matrix computations over C with 2, and basic measure theory (a.e. vanishing of nonnegative integrands).
Selected references
A. Drezet, An Elementary Proof That Everett's Quantum Multiverse Is Nonlocal: Bell-Locality and Branch-Symmetry in the Many-Worlds Interpretation, arXiv:2306.07794, 2023. https://arxiv.org/abs/2306.07794
D. M. Greenberger, M. A. Horne, A. Zeilinger, Going beyond Bell's theorem, in Bell's Theorem, Quantum Theory and Conceptions of the Universe, Kluwer, 1989, pp. 69–72. https://arxiv.org/abs/0712.0921
Verlinde 2016: Apparent Dark Matter and the Baryonic Tully-Fisher RelationResearch Paper
Motivation
In emergent gravity proposals, spacetime and gravity are macroscopic descriptions of an underlying quantum system. In Verlinde (2016) it is argued that in de Sitter space (a universe with positive dark energy) the entanglement entropy has, besides the usual area-law part, a thermal volume-law part, and that matter displaces this volume-law entropy. The elastic response of the medium to this displacement produces an extra gravitational pull, which the paper proposes as the origin of the effects usually ascribed to dark matter in galaxies and clusters.
The physical derivation is heuristic, but it ends in explicit mathematical relations between the baryonic (visible) mass profile MB(r) and the apparent dark matter mass profile MD(r) of a spherically symmetric, isolated system. The central one is the integral relation (7.40),
∫0rr′2GMD(r′)2dr′=6MB(r)a0r,
where G is Newton's constant and a0=cH0 is the Hubble acceleration scale. For a point mass it reproduces Milgrom's empirical acceleration relation and the baryonic Tully-Fisher relationvf4=aMGMB with aM=a0/6.
Setting
All quantities are real numbers or real functions of the radius r>0. G>0 and a0>0 are constants. MB(r) and MD(r) are the baryonic and apparent dark masses inside radius r. The Newtonian accelerations are gB(r)=GMB(r)/r2 and gD(r)=GMD(r)/r2 (7.42). The averaged density inside radius r is ρˉ(r)=3M(r)/(4πr3) (7.45), and the slope parameter is β(r)=−dlogρˉ/dlogr=−rρˉ′(r)/ρˉ(r) (7.46). A local density ρ gives M(r)=∫0rρ(r′)4πr′2dr′ (7.39).
The paper also uses two purely mathematical facts about elastic strain in d−1 spatial dimensions: a bound (7.30) on the largest principal strain of a traceless symmetric strain tensor, and the integration-by-parts identity (7.31) between the squared Hessian and the squared Laplacian.
Formalization targets
Goal: Milgrom relation for a central point mass (7.43)
If MB>0 is constant and (7.40) holds for every r>0, then for every r>0
gD(r)=aMgB(r),aM=6a0.
Milestones
(4.44) the strain-squared integral outside a spherical inclusion: ∫c∞(c/V)2dV=c.
(7.30) ε2≤d−1d−2∑ijεij′2 for a traceless symmetric ε′ with eigenvalue ε.
(7.31) ∫(∇i∇jχ)2dV=∫(∇2χ)2dV.
(7.47) ρˉD(r)2=(4−βB(r))8πGra0ρˉB(r) for a general baryonic profile.
(7.48) ρD(r)=(1−31βD(r))ρˉD(r).
(7.51) ΩD2=34ΩB when (7.47) is applied at the Hubble radius with βB=0.
Significance
The point-mass relation (7.43) is what connects the paper to the observed flat rotation curves and the baryonic Tully-Fisher law; (7.47) is its extension to extended mass profiles, which the paper uses to discuss galaxy clusters. The formal content checked here is the mathematics that turns the stated relation (7.40) into these observable consequences, plus the two elasticity lemmas used on the way to (7.40). The physical steps that lead to (7.40) (the entropy-displacement picture, the elastic interpretation, the identification (7.35)) are not mathematical statements and are outside the scope of this mission. None of these relations is known to have been machine-checked before.
Difficulty
The algebra is short once one may differentiate; the work is in the analysis. (7.40) is an integral identity with a singular weight 1/r′2 at the origin, so recovering pointwise information requires the fundamental theorem of calculus on (0,∞) with only continuity there, and care with integrability at 0. (7.31) is a double integration by parts in Rn that must be justified for the chosen decay class. (7.30) needs the spectral theorem for symmetric matrices, or an equivalent argument.
Formalization scope
All statements live in the namespace Verlinde2016. Radial profiles are functions R→R, with hypotheses only on (0,∞). The integral in (7.40) is Lean's interval integral over [0,r]; since it would be 0 for a non-integrable integrand, the positivity of the right-hand side forces integrability. Continuity of MD (and of densities) on (0,∞) is taken as a standing physical convention. In (7.31) the paper's condition that "χ falls off rapidly enough" is formalized as χ smooth with compact support. Spatial dimension is d−1 in (7.30); (7.40)–(7.51) use d=4 as in Section 7.2. Definitions: avgDensity, slopeParam, partialDeriv.
M. Milgrom, A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis, Astrophys. J. 270 (1983) 365. https://doi.org/10.1086/161130
8 thms3 active usersReviewed
🏆Completed
Captain: Lucas
Liouville Field Theory I: Spectrum and Kinematic SymmetriesTextbook
Motivation
Liouville field theory is a two-dimensional conformal field theory whose classical equation of motion generalizes Liouville's equation. It is one of the few interacting conformal field theories with a continuous spectrum that has been solved: its three-point function on the sphere is given by the DOZZ formula (Dorn–Otto 1994; Zamolodchikov–Zamolodchikov 1996), proved rigorously from a probabilistic construction of the path integral by Kupiainen–Rhodes–Vargas (2017), with the conformal bootstrap established by Guillarmou–Kupiainen–Rhodes–Vargas (2020). It appears in two-dimensional quantum gravity, non-critical string theory, and log-correlated random energy models.
This mission is the first in a series following the Wikipedia article on Liouville field theory. It formalizes the kinematic layer of the theory: the algebraic relations between the coupling constant, the background charge, the central charge, the conformal dimensions and the reflection coefficient, on which every later statement (the spectrum, the DOZZ formula, the bootstrap) is phrased.
Setting
Fix a complex number b, the coupling constant. Define
the background chargeQ=b+b1;
the central chargec=1+6Q2;
the conformal dimension of the primary field with momentumα∈C, Δ(α)=α(Q−α);
the reflection coefficient with normalization parameter λ,
RP=±λ−2iPΓ(−2ibP)Γ(−2ib−1P)Γ(2ibP)Γ(2ib−1P),
with sign +1 if c∈(−∞,1) and −1 otherwise.
Two symmetries act on these data: the reflectionα→Q−α and the dualityb→1/b. The spectrum of the theory is the set of dimensions Δ∈24c−1+R+, parametrized by momenta α=2Q+iP, P∈R+.
Formalization targets
Goal: the spectrum parametrization
For every b∈C, the map P↦Δ(2Q+iP) is a bijection
[0,∞)∼24c−1+[0,∞)⊂C.
Milestones
Reflection invariance: Δ(Q−α)=Δ(α).
Duality: c and Δ are invariant under b→1/b.
Marginality: for b=0, b(Q−b)=1⟺Q=b+1/b.
Dimension in terms of momentum: Δ(2Q+iP)=24c−1+P2.
Imaginary coupling: real nonzero β=ib produce exactly the real central charges c≤1.
Reflection relation consistency: RPR−P=1 for λ=0 and 2ibP,2ib−1P∈/Z.
Significance
These identities are the dictionary between the parameters used in the physics literature (b, Q, c, α, P, Δ) and fix the conventions for later missions in the series: the DOZZ formula, the pole structure of correlation functions, and the reflection and duality covariance of structure constants are all stated in these variables. None of these statements is open; the value of the mission is a vetted, reusable definition layer with machine-checked consistency of the conventions.
Difficulty
The individual statements are elementary algebra over C, together with the non-vanishing of the Gamma function away from its poles. The care needed is in conventions: total functions in Lean (0−1=0, principal-branch complex powers, Gamma at its poles) must be handled by explicit hypotheses exactly where the source implicitly excludes degenerate values.
Formalization scope
All quantities are complex-valued functions of b∈C, defined in the namespace LiouvilleFieldTheory. 1/b is Lean's inverse (so b=0 gives junk values Q=0, c=1; statements that need b=0 assume it). R+ is taken to be [0,∞). The condition c∈(−∞,1) is encoded as Imc=0∧Rec<1. λ−2iP is the principal-branch complex power. The analytic content of the article (Virasoro representations, conformal blocks, the special function Υb, the path integral, the DOZZ formula) is out of scope for this mission.
The entropy photon-number inequalityResearch Paper
Motivation
Shannon's entropy power inequality bounds the entropy of a sum of independent random vectors below by the entropies of the summands; it is a basic tool for classical Gaussian channel capacities. In optical communication, addition is replaced by passive mixing of bosonic modes at a beam splitter. Guha, Erkmen and Shapiro conjectured the sharp bosonic analogue, the entropy photon-number inequality (EPnI), stated in terms of the photon number of the thermal state with the same entropy (ITA 2008). If true, it settles the minimum output entropy of thermal attenuators at fixed input entropy and the capacity regions of bosonic broadcast and wiretap channels, which had been reduced to this inequality (PRA 2007).
This mission asks for a formal proof of the EPnI for finite-energy inputs in any finite number of modes, as stated in an OpenAI preprint dated September 24, 2026 (source). The preprint has not been peer reviewed and its claims have not been independently verified; on the platform the Lean goal is open.
Background
1948–1965 — Shannon introduces entropy power; Stam (1959) and Blachman (IEEE IT 1965) prove the entropy power inequality.
2004 — Capacity of the pure-loss bosonic channel (PRL 2004) and the minimum-output-entropy conjecture (PRA 2004).
2007–2008 — Guha, Shapiro and Erkmen relate broadcast capacity to a minimum-output-entropy conjecture; Guha, Erkmen and Shapiro state the EPnI.
2011 — Special cases (one-mode number-diagonal inputs) by Das, Sharma and Muthukrishnan (arXiv:1107.2365).
2014 — König and Smith prove a quantum entropy power inequality, including S(ρC)≥ηS(ρA)+(1−η)S(ρB) (IEEE IT 2014); De Palma, Mari and Giovannetti prove the exponential form eS(ρC)/n≥ηeS(ρA)/n+(1−η)eS(ρB)/n (Nat. Photonics 2014).
2015 — Giovannetti, Holevo and García-Patrón solve the Gaussian optimizer conjecture for unconstrained minimum output entropy (CMP 2015).
2016–2018 — De Palma, Trevisan and Giovannetti settle one-mode constrained problems (IEEE IT 2017; PRL 2017); De Palma and Trevisan give a rigorous finite-energy conditional EPI (CMP 2018).
September 2026 — The OpenAI preprint claims the two-input EPnI for all finite-energy inputs (Theorem 1.1, p. 2).
Setting
An n-mode bosonic system has Fock space ℓ2(Nn) with number basis ∣k⟩, k∈Nn. A state is a positive trace-one operator ρ; it has finite energy if Trρ∑jaj†aj=∑k(k1+⋯+kn)⟨k∣ρ∣k⟩<∞. The von Neumann entropy is S(ρ)=−Trρlogρ (natural log). The thermal entropy per mode is
g(t)=(t+1)log(t+1)−tlogt,t≥0,
increasing from 0 to ∞, and the entropy photon number of ρ is N(ρ)=g−1(S(ρ)/n), the mean photon number per mode of the product thermal state with the same entropy. For independent inputs ρA,ρB and transmissivity 0≤η≤1, the beam-splitter outputρC is the reduced state of the modes cj=ηaj+1−ηbj after passive mixing of ρA⊗ρB. No independence among the modes within either input is assumed.
Formalization targets
Goal: Theorem 1.1 (p. 2)
For every finite n≥1, every pair of finite-energy n-mode states ρA,ρB and every 0≤η≤1,
g−1(nS(ρC))≥ηg−1(nS(ρA))+(1−η)g−1(nS(ρB)).
Thermal inputs with a common mean photon number within each port give equality.
Significance
The result itself. The EPnI is sharper than the exponential quantum entropy power inequality when the two input entropies differ, which is exactly the regime needed for constrained minimum output entropy. The preprint derives the exact minimum output entropy at fixed input entropy for every tensor power of an identical thermal attenuator, with inputs entangled across modes (Corollary 1.2, p. 2), and the exact capacity region of the degraded two-receiver pure-loss broadcast channel (Corollary 7.1, p. 26), completing the converse of Guha, Shapiro and Erkmen.
Formalizing it. The statement involves infinite-dimensional Fock space, trace-class positive operators, continuous functional calculus for the entropy, and the partial trace after a passive unitary. The Lean development builds these concretely in the number basis. A formal proof would give machine-checked infrastructure for quantum entropy on infinite-dimensional spaces, which Mathlib currently lacks. No machine-checked version of any quantum entropy power inequality is known.
Difficulty
Heat-flow (quantum Fisher information) proofs give the linear and exponential bounds, but these are attained only for equal entropies and do not recover the thermal comparison at unequal entropy levels. The preprint instead minimizes a regularized functional over pairs of inputs and must show that minimizers are thermal. This requires an interpolation theorem for diagonal form comparisons (Theorem 3.1, p. 8), attainment and Gibbs regularity of minimizers in infinite dimensions (Lemma 4.1, Proposition 4.2, pp. 15–16), a component entropy-Hessian estimate, and a stationarity identity whose cancellation forces convergence to thermal inputs (Section 6). Compactness relies on the finite-energy hypothesis; without it entropies can be infinite.
Formalization scope
Fock n is lp (fun _ : Fin n → ℕ => ℂ) 2. A State is a positive bounded operator whose diagonal in the number basis sums (HasSum) to 1.
FiniteEnergy ρ is summability of ∑k∣k∣⟨k∣ρ∣k⟩.
entropy ρ is the sum of diagonal entries of cfc(−tlogt)(ρ); for finite-energy states it is finite, so the tsum is the true entropy.
gInv s = sInf {t ≥ 0 | s ≤ g t}, the inverse of g on [0,∞).
IsBeamSplitterOutput η ρA ρB ρC requires every number-basis entry of ρC to be the (convergent) partial-trace sum of the explicitly given beam-splitter coefficients applied to ρA⊗ρB; the inputs are arbitrary density matrices, so internal entanglement is allowed.
The hypotheses are satisfiable (e.g. vacuum or thermal inputs) and the conclusion is the inequality as stated.
S. Guha, B. I. Erkmen, J. H. Shapiro, The entropy photon-number inequality and its consequences, ITA Workshop 2008, 128–130. https://doi.org/10.1109/ITA.2008.4601037
S. Guha, J. H. Shapiro, B. I. Erkmen, Classical capacity of bosonic broadcast communication and a minimum output entropy conjecture, Phys. Rev. A 76 (2007), 032303. https://doi.org/10.1103/PhysRevA.76.032303
G. De Palma, A. Mari, V. Giovannetti, A generalization of the entropy power inequality to bosonic quantum systems, Nat. Photonics 8 (2014), 958–964. https://doi.org/10.1038/nphoton.2014.252
G. De Palma, D. Trevisan, The conditional entropy power inequality for bosonic quantum systems, Comm. Math. Phys. 360 (2018), 639–662. https://doi.org/10.1007/s00220-017-3082-8
V. Giovannetti, A. S. Holevo, R. García-Patrón, A solution of Gaussian optimizer conjecture for quantum channels, Comm. Math. Phys. 334 (2015), 1553–1571. https://doi.org/10.1007/s00220-014-2150-6
V. Giovannetti, S. Guha, S. Lloyd, L. Maccone, J. H. Shapiro, H. P. Yuen, Classical capacity of the lossy bosonic channel: the exact solution, Phys. Rev. Lett. 92 (2004), 027902. https://doi.org/10.1103/PhysRevLett.92.027902
Varying Constants I: Dimensional Analysis, Natural Units and the Buckingham π TheoremResearch Paper
Motivation
Uzan's review Varying Constants, Gravitation and Cosmology (Living Rev. Relativity 14 (2011) 2) is about testing whether the fundamental constants of physics change in space or time. Before any experiment can be read, §2.1 of the review settles which constants it even makes sense to ask about. Two claims from that section carry the whole programme:
p. 14 (§2.1.2): "only the variation of dimensionless constants can be measured and in case such a variation is detected, it is impossible to determine, which dimensional constant is varying";
pp. 15–17 (§2.1.1): three independent constants (e.g. Stoney's G,e,c or Planck's c,G,ℏ) can be used to define the three mechanical units, after which "all other constants are dimensionless quantities" whose values do not depend on the units, and "any variation of constants that will leave these numbers unaffected is actually just a redefinition of units".
Both claims are informal statements of a classical result of dimensional analysis, the Buckingham π theorem (E. Buckingham, Phys. Rev. 4 (1914) 345). This mission states them precisely and asks for formal proofs.
Setting
Fix d base units U1,…,Ud (for mechanics d=3: L,M,T) and n dimensional constants. The dimension matrixD∈Rn×d says that constant i has dimension ∏jUjDij; real exponents are allowed, because Gaussian units give the charge the dimension M1/2L3/2T−1 (Uzan, p. 15). A configurationx∈R>0n lists the numerical values of the constants in some system of units.
A change of units is a vector s∈R>0d. It acts by (s⋅x)i=xi∏jsjDij.
An observablef:Rn→R is unit invariant if f(s⋅x)=f(x) for all positive x,s.
The dimensionless exponents are ZD={a∈Rn:∑iaiDij=0∀j}, and πa(x)=∏ixiai is the corresponding power product.
Given d chosen constants e1,…,ed, a system of natural units for x is a positive s with (s⋅x)ej=1 for all j.
The Planck dimension matrixDP has rows [c]=LT−1, [G]=L3M−1T−2 and [ℏ]=L2MT−1.
Target
Goal (Buckingham π). Put k=n−rankD. There are linearly independent a(1),…,a(k)∈ZD such that every unit-invariant f factors as
f(x)=F(πa(1)(x),…,πa(k)(x))(x∈R>0n).
Milestones, in order:
Dimensionless power products are unit invariant: a∈ZD⇒πa(s⋅x)=πa(x).
dimZD=n−rankD.
If det(Dej,k)=0, every positive x has exactly one system of natural units.
Values in natural units do not depend on the system of units you started from: s′⋅(r⋅x)=s⋅x.
Planck units: c=G=ℏ=1 holds exactly when s=(1/ℓP,1/mP,1/tP), where ℓP=Gℏ/c3, mP=ℏc/G and tP=Gℏ/c5 (Uzan, p. 16).
Two positive configurations agree on every πa with a∈ZD⟺ they differ by a change of units (Uzan, p. 17).
Only dimensionless variations are measurable: a unit-invariant f takes equal values on configurations that agree on all dimensionless πa (Uzan, p. 14).
Significance
The result. Milestones 6–7 justify the review's working convention of quoting all bounds as bounds on dimensionless numbers (αEM, μ=mp/me, αG=Gmp2/ℏc, …). Milestone 2 counts how many independent numbers there are. Milestones 3–5 say that natural units are well defined.
Formalizing it. The mathematics is classical and the proofs are known. What remains is to formalize it: a reusable, general dimensional-analysis layer (unit actions, dimensionless monomials, natural units) that later missions on Uzan's review can build on.
Difficulty
Taking logarithms turns the multiplicative problem into linear algebra: changes of units act on logx by translations through the column space of D, and ZD is the annihilator of that space. Two steps carry the content. The first is the duality step (annihilator of the annihilator) behind milestone 6 and the goal. The second is the bookkeeping with real powers (Real.rpow), which only behaves well on positive bases. The goal also needs a choice of F with no regularity, defined through the quotient Rn/col(D).
Formalization scope
Lean namespace VaryingConstants; all definitions are in one definition file. Indices are Fin n and Fin d, and exponents are real.
Positivity is an explicit hypothesis (IsPositive), not a subtype. Unit invariance constrains f only on positive configurations.
ZD is the kernel of a↦aTD, and the rank is Matrix.rank over R.
Natural units are indexed by a map e : Fin d → Fin n. The independence hypothesis is detDe=0.
There is no trivializing reading: the goal fixes the number of π-groups to n−rankD and requires them to be linearly independent and dimensionless, so the identity map does not qualify.
Contributions such as general lemmas on Real.rpow products and on annihilators of matrix column spaces are reusable beyond this mission.
E. Buckingham, On physically similar systems; illustrations of the use of dimensional equations, Phys. Rev. 4 (1914) 345. https://doi.org/10.1103/PhysRev.4.345
Mathematical Methods in Quantum Mechanics V: Stone's Theorem and the RAGE TheoremTextbook
Motivation
In quantum mechanics the state of a system at time t is ψ(t)=e−itAψ, where the HamiltonianA is a self-adjoint operator. The spectral theorem splits the Hilbert space into subspaces according to the type of the spectral measure of a vector (absolutely continuous, singularly continuous, pure point). That splitting is defined through measure theory, and it is not evident what it has to do with the motion of the system. The RAGE theorem, named after Ruelle (1969), Amrein and Georgescu (1973) and Enß (1978), answers this. Vectors in the continuous subspace are exactly the states that, on time average, leave every bounded region: the scattering states. Vectors in the pure point subspace are exactly the states that stay in a bounded region uniformly in time: the bound states. This dynamical characterization is the starting point of the geometric approach to scattering theory and of the modern theory of quantum dynamics and transport.
This mission formalizes Chapter 5 of Teschl's Mathematical Methods in Quantum Mechanics (AMS GSM 99, 2009). The chapter covers the time evolution e−itA as a unitary group, Wiener's theorem on Fourier transforms of measures, and the abstract RAGE theorem with its Heisenberg-picture corollaries.
Setting
Let H be a complex Hilbert space with inner product ⟨⋅,⋅⟩, conjugate linear in the first argument. An operatorA is a linear map defined on a subspace D(A), its domain. A is self-adjoint if D(A) is dense and A=A∗. L(H) denotes the bounded, everywhere defined operators.
A projection-valued measureP assigns to each Borel set Ω⊆R an orthogonal projection P(Ω), with P(R)=I and P(⋃nΩn)ψ=∑nP(Ωn)ψ for disjoint Ωn. The spectral measure of ψ is the finite Borel measure μψ(Ω)=⟨ψ,P(Ω)ψ⟩. The spectral theorem says that every self-adjoint A is A=∫λdP(λ) for a unique P=PA. Then D(A)={ψ∣∫λ2dμψ<∞} and ⟨ψ,Aψ⟩=∫λdμψ(λ). For a bounded Borel function f, f(A) is the bounded operator with ⟨ψ,f(A)ψ⟩=∫fdμψ. The time evolution is U(t)=e−itA, the case f(λ)=e−itλ.
The spectral subspaces are
Hac={ψ∣μψ≪Lebesgue},Hsc={ψ∣μψ singular, without atoms},Hpp={ψ∣μψ supported on a countable set},
and the continuous subspace is Hc=Hac⊕Hsc. The projectors Pac,Pc,Ppp are the orthogonal projections onto these subspaces. σp(A) is the set of eigenvalues of A.
A finite rank operator is a K∈L(H) with finite dimensional range. The compact operatorsC(H) are the norm closure of the finite rank operators. For z in the resolvent set ρ(A), RA(z)=(A−z)−1∈L(H). An operator K is relatively compact with respect to A if KRA(z)∈C(H) for one z∈ρ(A); in particular D(A)⊆D(K). The Fourier transform of a finite complex Borel measure μ on R is μ^(t)=∫Re−itλdμ(λ).
Formalization targets
Goal: the RAGE theorem (Theorem 5.7)
Let A be self-adjoint and let Kn∈L(H) be relatively compact operators converging strongly to I. Then
Theorem 5.1.U(t)=e−itA is a strongly continuous one-parameter unitary group. limt→0t1(U(t)ψ−ψ) exists iff ψ∈D(A), and then it equals −iAψ. U(t)D(A)=D(A) and AU(t)=U(t)A.
Theorem 5.4 (Wiener).T→∞limT1∫0T∣μ^(t)∣2dt=λ∈R∑∣μ({λ})∣2, and the sum is finite.
Lemma 5.5.C(H) is a closed ∗-ideal in L(H).
Theorem 5.6. For relatively compact K and ψ∈D(A): T1∫0T∥Ke−itAPcψ∥2dt→0 and ∥Ke−itAPacψ∥→0. For bounded K this holds for all ψ.
Theorem 5.8.T→∞limT1∫0TeitAKe−itAψdt=λ∈σp(A)∑PA({λ})KPA({λ})ψ for ψ∈D(A), and for all ψ if K is bounded.
Corollary 5.9. Under the RAGE assumptions, the iterated time averages of eitAKne−itAψ and of eitA(I−Kn)e−itAψ converge to Pppψ and to Pcψ.
Significance
The RAGE theorem connects the spectral types of the Hamiltonian with observable long-time behaviour. For Schrödinger operators −Δ+V on L2(Rn), taking Kn to be multiplication by the indicator of the ball of radius n gives a sequence of relatively compact operators. The theorem then says that states in Hc spend, on time average, a vanishing fraction of time in every ball, while states in Hpp stay uniformly localized. Enß's proof of asymptotic completeness for short-range potentials rests on this characterization. The same argument, through Wiener's theorem, is used to relate the fractal dimensions of spectral measures to rates of quantum transport. Theorem 5.8 and Corollary 5.9 describe the asymptotic observables of the Heisenberg picture.
These results are classical and proved in the literature. They are not formalized. Mathlib has neither a spectral theorem nor a functional calculus for unbounded self-adjoint operators, and no decomposition of a Hilbert space by spectral type. What the platform has are Stone's theorem in the separate BookProof development (for weakly measurable groups, with its own structure of unbounded self-adjoint operators) and the discrete-time von Neumann mean ergodic theorem for contractions. The mean ergodic theorem is related to Theorem 5.8 with K=I but is a different statement. This mission formalizes the continuous-time, spectral-type-resolved statements.
Difficulty
The two sides of the RAGE characterization are defined in unrelated terms. Hc and Hpp are defined by properties of the scalar measures μψ, while the conditions in (5.14) concern norms of vectors moving under a unitary group and seen through operators that are only relatively compact. Wiener's theorem controls scalar quantities ⟨φ,e−itAψ⟩. Passing to norms ∥Ke−itAψ∥ needs compactness, and relative compactness only gives compactness after composing with a resolvent, which is where the restriction to D(A) in Theorems 5.6 and 5.8 comes from. The pure point characterization needs control that is uniform in t≥0, so a limit over n has to be exchanged with a supremum over all times. The limits in (5.14) and (5.19)–(5.20) are iterated, not joint, and a naive exchange of the two limits is not justified. Summing over σp(A) in (5.18) needs an unconditionally convergent sum over a possibly uncountable index set. Wiener's theorem itself needs a Fubini argument for a complex measure against its conjugate.
Formalization scope
The formalization uses Lean 4 with Mathlib. Operators are partially defined linear maps (LinearPMap, H →ₗ.[ℂ] H) on a complex Hilbert space. Self-adjointness is Mathlib's IsSelfAdjoint, which includes density of the domain. The Hilbert space is not assumed separable.
The spectral theorem is not part of this mission. Every statement about A takes a projection-valued measureP on the Borel sets of R as data, together with the hypothesis that A=∫λdP(λ): D(A)={ψ∣∫λ2dμψ<∞} and ⟨ψ,Aψ⟩=∫λdμψ there. This is the form the spectral theorem guarantees, and it makes PA=P. The time evolution e−itA is the bounded operator P(λ↦e−itλ), determined by its quadratic form, and eitA is its value at −t. Nothing else the book derives is assumed.
The spectral subspaces are sets of vectors defined only from the spectral measures, as above, and Hc is the set of sums of a vector of Hac and a vector of Hsc. A formalization that defined Hc or Hpp by the dynamical conditions of (5.14) would make the goal a tautology; that is ruled out here.Pc, Pac, Ppp are the orthogonal projections onto the closed linear spans of these sets. Compact operators are the norm closure of the finite rank operators (the book's definition), not Mathlib's IsCompactOperator. Relative compactness asks for one z∈ρ(A), its resolvent R, and a compact C with RanR⊆D(K) and KR=C. "K bounded" means an everywhere defined K∈L(H).
Iterated limits limnlimT are stated as: for every n the inner limit exists, and the sequence of inner limits converges. The supremum over t≥0 is taken in [0,∞]. Time averages use interval integrals, with Bochner integrals for vector-valued integrands; all integrands are continuous in t. Wiener's theorem uses Mathlib's ComplexMeasure ℝ and its vector-measure integral. The sum over λ∈R is an unconditional sum whose summability is part of the conclusion, and the sum over σp(A) is an unconditional sum in H. The Fourier transform has the book's normalization.
A complete development needs a functional calculus for projection-valued measures (multiplicativity, strong continuity in the function argument), the reduction of A by the spectral subspaces (Lemma 3.19), the Riemann–Lebesgue lemma for absolutely continuous measures, and the approximation of compact operators by finite rank ones. These pieces are reusable across the rest of the book, notably for scattering theory. Proofs of any milestone, and of such supporting lemmas, are welcome.
Selected references
G. Teschl, Mathematical Methods in Quantum Mechanics, With Applications to Schrödinger Operators, Graduate Studies in Mathematics 99, American Mathematical Society, 2009. doi:10.1090/gsm/099. Chapter 5, pp. 123–132.
D. Ruelle, A remark on bound states in potential-scattering theory, Il Nuovo Cimento A 61 (1969), 655–662. doi:10.1007/BF02819607.
W. O. Amrein and V. Georgescu, On the characterization of bound states and scattering states in quantum mechanics, Helvetica Physica Acta 46 (1973), 635–658 (archive: e-periodica.ch).
V. Enss, Asymptotic completeness for quantum mechanical potential scattering. I. Short range potentials, Communications in Mathematical Physics 61 (1978), 285–291. doi:10.1007/BF01940771.
M. H. Stone, On one-parameter unitary groups in Hilbert space, Annals of Mathematics 33 (1932), 643–648. doi:10.2307/1968538.
Mathematical Methods in Quantum Mechanics X: The Uncertainty Principle and the Harmonic OscillatorTextbook
Motivation
The quantum harmonic oscillator is the Schrödinger operator of a particle in a quadratic potential. It describes small vibrations around an equilibrium, the modes of the electromagnetic field, and phonons in a crystal. It is also one of the very few Schrödinger operators whose spectrum and eigenfunctions are known in closed form. The eigenvalues are equally spaced, and the eigenfunctions are Hermite functions. Its algebraic treatment by creation and annihilation operators goes back to Dirac (The Principles of Quantum Mechanics, 1930). That treatment became the template for second quantization and Fock space.
The same chapter of G. Teschl's Mathematical Methods in Quantum Mechanics (AMS Graduate Studies in Mathematics 99, 2009, Chapter 8, AMS page) contains two further results that use commutation relations. The first is the Heisenberg uncertainty principle (Heisenberg 1927; the general form for two observables is due to Robertson, Phys. Rev. 34 (1929)). The second is an abstract statement that A∗A and AA∗ have the same spectrum away from 0. Physicists call this supersymmetric quantum mechanics, and Teschl uses it later for the hydrogen atom. This mission follows Teschl's Sections 8.1–8.4.
Setting
Let H be a complex Hilbert space. A linear operatorA is a linear map defined on a subspace D(A)⊆H, its domain. A is symmetric if D(A) is dense and ⟨φ,Aψ⟩=⟨Aφ,ψ⟩ for all φ,ψ∈D(A). For a closed operator A and z∈C, the resolventRA(z)=(A−z)−1 exists when A−z maps D(A) bijectively onto H with a bounded inverse. The resolvent setρ(A) is the set of such z, and the spectrum is σ(A)=C∖ρ(A). An operator is essentially self-adjoint on its domain if its closure A is self-adjoint.
For ψ∈D(A) with ∥ψ∥=1 (a state), the expectation is Eψ(A)=⟨ψ,Aψ⟩ and the mean-square deviation is Δψ(A)=∥(A−Eψ(A))ψ∥. The commutator is [A,B]ψ=ABψ−BAψ for ψ∈D(AB)∩D(BA), where D(AB)={ψ∈D(B)∣Bψ∈D(A)}.
Fix a frequency ω>0. On R3 with Lebesgue measure, let
Dω=span{xαe−ω∣x∣2/2∣α∈N03}⊆L2(R3),
where xα=x1α1x2α2x3α3 and ∣x∣ is the Euclidean norm. With Δ=∑j∂2/∂xj2 and Teschl's units (ℏ=1, mass 1/2, so the free Hamiltonian is H0=−Δ), the harmonic oscillator is the operator
H=H0+ω2x2,(Hf)(x)=−Δf(x)+ω2∣x∣2f(x),D(H)=Dω.
In one dimension the ladder operators on span{xke−ωx2/2} are
A±=21(ωx∓ω1dxd).
The Hermite polynomials are Hn(x)=(−1)nex2dxndne−x2. The Hermite functions are
ψn(x)=2nn!1(πω)1/4Hn(ωx)e−ωx2/2,
and on R3, ψn1,n2,n3(x)=ψn1(x1)ψn2(x2)ψn3(x3).
Formalization targets
Goal: the harmonic oscillator (Theorem 8.5)
For every ω>0 the operator H above exists. It is essentially self-adjoint on Dω, the functions ψn1,n2,n3, (n1,n2,n3)∈N03, form an orthonormal basis of L2(R3) consisting of eigenvectors of H, and
σ(H)={(2n+3)ω∣n∈N0}.
Milestones
Density (Lemma 8.3).span{xαe−∣x∣2/2∣α∈N0n} is dense in L2(Rn) for every n.
Canonical commutation relation (Eq. (8.36)).[A−,A+]=1 on span{xke−ωx2/2}.
Factorization (Eq. (8.37)). On the same space, −dx2d2+ω2x2=ω(2N+1) with N=A+A−, and the space is invariant under A±.
Hermite functions (Eqs. (8.40)–(8.42)).ψn=A+nψ0/n! with ψ0(x)=(ω/π)1/4e−ωx2/2. Each ψn is a normalized eigenfunction of N for the eigenvalue n, and span{ψn}=span{xke−ωx2/2}.
Heisenberg uncertainty principle (Theorem 8.2). For symmetric A, B and a state ψ∈D(AB)∩D(BA),
Δψ(A)Δψ(B)≥21∣Eψ([A,B])∣,
with equality if (B−Eψ(B))ψ=iλ(A−Eψ(A))ψ for a real λ=0, or if ψ is an eigenvector of A or B.
Abstract commutation (Theorem 8.6). For A closed and densely defined, A∗A∣Ker(A)⊥ and AA∗∣Ker(A∗)⊥ are unitarily equivalent. A maps eigenvectors of A∗A (eigenvalue E) to eigenvectors of AA∗ with ∥Aψ0∥=E∥ψ0∥. For z=0,
The last two milestones do not enter the goal. They are the chapter's other two theorems and use the same operator-theoretic layer.
Significance
Theorem 8.5 is the rigorous form of the most used exactly solvable model of quantum mechanics. Once it is available, the Hermite functions become an orthonormal basis of L2(R3) that diagonalizes a concrete differential operator. That basis underlies Mehler's formula and the Fourier transform's action Fψn=(−i)nψn (Eq. (8.46)). Perturbation arguments around the oscillator also start from it, and so does the construction of Fock space. The uncertainty principle is the basic quantitative constraint on simultaneous measurement. Theorem 8.6 gives an eigenvalue transfer mechanism that Teschl reuses to compute the hydrogen spectrum (Section 10.4).
Formally, the mission asks for a machine-checked proof about an unbounded differential operator on L2(R3), stated on a concrete core. Mathlib has unbounded operators (LinearPMap, closure, adjoint), Lp spaces, Gaussian integrals and the probabilists' Hermite polynomials. It has no statement identifying the spectrum of a concrete Schrödinger operator. These results are proved in every textbook, but none of them has been formalized in this form.
Difficulty
The algebra of Eqs. (8.36)–(8.41) is a finite computation with polynomials times a Gaussian. The work in the goal is analytic. Essential self-adjointness is a statement about the closure of an operator on a non-closed domain. Symmetry of H together with a family of eigenvectors in Dω does not give it by definition, and the relevant criterion for unbounded operators is not in Mathlib. Completeness of the eigenfunctions in L2(R3) is a density statement about a concrete function space (Lemma 8.3), and it has to be transferred from one dimension to three. The spectrum is asked for the closure H, not for H on Dω: the eigenvalues (2n+3)ω are easy to exhibit, but showing that no other point of C lies in σ(H) requires control of (H−z)−1 on all of L2. Finally, the Rodrigues formula (8.42) and the ladder construction (8.40) define the same functions only after an identity for all derivatives of e−x2.
Formalization scope
H is a complex inner product space with CompleteSpace; separability is not assumed, and no statement needs it. Operators are Mathlib LinearPMaps H →ₗ.[ℂ] H. A∗ is LinearPMap.adjoint, and every statement using it assumes a dense domain. Self-adjointness is Mathlib's IsSelfAdjoint, and essential self-adjointness is IsSelfAdjoint H.closure. ρ and σ are defined in the mission for LinearPMaps (bounded two-sided inverse of A−z); Mathlib's Banach-algebra spectrum is not used. L2(Rn) is Lp ℂ 2 volume on EuclideanSpace ℝ (Fin n). The Laplacian is the sum of second Fréchet derivatives along the unit vectors. Hermite polynomials are defined by the Rodrigues formula, not as Mathlib's probabilists' Polynomial.hermite. No projection-valued measure, functional calculus or Fourier transform is taken as data.
The oscillator is not defined as a construction. IsHarmonicOscillator ω H states that D(H)=Dω and that H acts on a representative f as −Δf+ω2∣x∣2f. This determines H uniquely, and the goal also asserts that such an H exists. Trivializing formalization ruled out: defining H as the diagonal operator ∑n(2∣n∣+3)ω⟨ψn,⋅⟩ψn would make the spectrum immediate. Here H is the differential operator on Dω, and its eigenvalues must be derived.
Two choices differ from the printed text. First, Teschl's (8.34) writes the domain with e−x2/2, but his (8.39)–(8.41) and the conclusion span{ψn}=D use e−ωx2/2. The mission uses Dω, which equals (8.34) for ω=1. Lemma 8.3 is stated for (8.20)'s D itself. Second, Theorem 8.6 prints "H1ψ1=ψ1", and the mission states the intended H1ψ1=Eψ1. The uncertainty principle adds ∥ψ∥=1 explicitly; the book assumes it through its definition of a state.
Reusable pieces are the product of LinearPMaps, spectrum and resolvent for unbounded operators, the Gauss–polynomial core and its density, and the Hermite-function API. Contributions of general lemmas are welcome, for example "an operator with an orthonormal eigenbasis in its domain is essentially self-adjoint" and the spectrum of such an operator.
Selected references
G. Teschl, Mathematical Methods in Quantum Mechanics, With Applications to Schrödinger Operators, Graduate Studies in Mathematics 99, AMS, 2009. bookstore.ams.org/gsm-99
H. P. Robertson, "The Uncertainty Principle", Physical Review 34 (1929), 163–164. doi:10.1103/PhysRev.34.163
W. Heisenberg, "Über den anschaulichen Inhalt der quantentheoretischen Kinematik und Mechanik", Zeitschrift für Physik 43 (1927), 172–198. doi:10.1007/BF01397280
In quantum field theory a Ward–Takahashi identity is an identity between correlation functions that follows from a global or gauge symmetry of the theory. In quantum electrodynamics it was used by Ward (1950) and Takahashi (1957) to relate the wave-function renormalization of the electron to its vertex renormalization, and it is described in the source article as "the quantum version of classical current conservation associated to a continuous symmetry by Noether's theorem". In the path-integral formulation the identities are a reflection of the invariance of the functional measure under the symmetry; when the measure is not invariant one obtains the anomalous Ward–Takahashi identity (e.g. the chiral anomaly).
The continuum path integral of interacting QED is not a mathematically defined object, so the identities cannot be formalized literally. This mission formalizes the path-integral derivation in the setting where every step is rigorous: a finite lattice, where the functional measure is Lebesgue measure on a finite-dimensional space.
Setting
A field configuration of a complex scalar field on N sites is φ∈CN. The global U(1) symmetry acts by (Rθφ)y=eiθφy. The local generator at site x is (gxφ)y=δxyiφx, and the local variation of a real-differentiable functional F is δxF(φ)=DF(φ)[gxφ]. An action S:CN→R is U(1)-invariant if S∘Rθ=S for all θ. The (unnormalised, Euclidean) path integral of an observable F is
ZS[F]=∫CNF(φ)e−S(φ)dφ.
For invariant S, jx=δxS is the lattice divergence of the Noether current.
Formalization targets
Goal: two-point Ward–Takahashi identity
For a C1, U(1)-invariant action with the moment conditions (1+∥φ∥3)e−S∈L1 and ∥φ∥3∥DS∥e−S∈L1: current conservation ∑xjx=0, and for all sites x,y,z
ZS[jxφyφz]=i(δxy−δxz)ZS[φyφz].
Milestones
Invariance of Lebesgue measure under Rθ.
Lattice Noether theorem: ∑xδxS=0 for invariant S.
Infinitesimal measure invariance: ∫δxGdφ=0.
Local Ward–Takahashi identity: ZS[δxF]=ZS[FδxS].
Global Ward identity: ∑xZS[δxF]=0 for invariant S.
Anomalous identity for a linear generator A: ∫(DF[Aφ]−FDS[Aφ])e−S=−tr(A)ZS[F].
Significance
The goal is the lattice version of the statement that the divergence of the current inserted into a charged correlation function produces only contact terms — the structure of the QED identity in the source, whose right-hand side consists of amplitudes with shifted external momenta. The milestones isolate the three ingredients of the textbook derivation: invariance of the measure, Noether current conservation, and integration by parts; the anomalous version identifies the anomaly with the Jacobian (trace) of the transformation. The results are classical; to the drafter's knowledge they are not available in Mathlib in this form.
Difficulty
The identities are elementary on paper; the work is analytic bookkeeping. Integration by parts must be justified on the non-compact space CN with non-constant vector fields φ↦gxφ or Aφ, so the "surface terms vanish" step of the physics derivation has to be replaced by explicit integrability hypotheses, and products with complex conjugates require real (not complex) differentiability throughout.
Formalization scope
Fields are Fin N → ℂ with the sup norm and the product Lebesgue measure; derivatives are real Fréchet derivatives (fderiv ℝ); integrals are Bochner integrals, which take the value 0 on non-integrable integrands. The weight is Euclidean, e−S, and all identities are stated for unnormalised integrals (normalised expectations are ratios by ZS[1]). Hypotheses are satisfiable: e.g. any Gaussian action S(φ)=∑y∣φy∣2 is U(1)-invariant and satisfies all moment conditions, so the goal is not vacuous. Continuum QED, renormalization and the S-matrix Ward identity are out of scope.
Susskind 2022: Separation of Scales in Double-Scaled SYKResearch Paper
Motivation
A recent conjecture proposes a holographic duality between a static patch of Jackiw–Teitelboim (JT) gravity with positive cosmological constant, a two-dimensional reduction of three-dimensional de Sitter gravity, and the double-scaled SYK model at infinite temperature (DSSYK∞). The conjecture comes from Susskind (arXiv:2109.14104, arXiv:2205.00315), Lin–Susskind (arXiv:2206.01083) and Rahman. In the paper arXiv:2209.09999, Susskind looks at the semiclassical limit of de Sitter space, where the horizon entropy diverges. There physics splits into a cosmic sector, sensitive to the de Sitter radius Lc, and a microscopic (string) sector, which does not see it. The paper argues that the N→∞ limit of DSSYK∞ shows the same split, with the two sectors separated by a factor of q. This parallel is offered as evidence for the duality.
Most of the paper's arguments are heuristic scaling estimates. Several of them rest on explicit closed-form expressions: the epidemic scrambling function, the large-q correlators, and the second moments of the Hamiltonians. These can be stated and proved as precise limit theorems. This mission collects them.
Setting
The SYK model has N Majorana fermions χi with χi2=1, coupled by q-body interactions with Gaussian random couplings. The cosmic Hamiltonian is Hc=∑∣I∣=qjIχI, where the (qN) couplings jI have variance ⟨jj⟩=q!J2/Nq−1 (eq. (3.4)), and J>0 is a fixed energy. The string Hamiltonian is Hs=Hc/q (eq. (3.7)). Cosmic timetc and string timets are related by ts=qtc.
The double-scaled limit is a sequence of models (Nn,qn) with
Nn→∞,Nnqn2→λ>0(eq. (3.3)).
In particular qn→∞. Infinite temperature means expectations are taken in the maximally mixed state.
In DSSYK∞ scrambling is described by the epidemic scrambling probability (eq. (8.1))
PN,q(tc)=1−(1+Nqe(q−1)Jtc)−q−11.
The infinite-temperature second moments of the two Hamiltonians are
Early scrambling is ordinary fast scrambling in string units (8.3). Late scrambling is a quasinormal-mode decay in cosmic units (8.6).
Milestones
Eq. (8.3): NP(ts/q)→eJts (fast scrambling in string units).
Eq. (8.6): 1−P(tc)→e−Jtc for tc>0 (hyperfast scrambling in cosmic units).
Eqs. (8.4)–(8.5): at Δtc=logq/(qJ), qNP→λlog(1+λ), so fast scrambling ends at P∼q/N.
Eqs. (4.1)–(4.3): ∫Rcosh−2(Jqt)dt=2/(Jq), the emergent string scale Ls=1/(Jq) as IR regulator.
Eqs. (9.3)–(9.4): (cosh−2(Jqtc))1/q→e−2∣J∣∣tc∣, the fermion correlator decaying on the cosmic scale.
Eq. (11.1): ⟨Hc2⟩/N→J2e−λ/2 and ⟨Hc2⟩→∞.
Eq. (11.2): ⟨Hs2⟩→J2e−λ/2/λ (finite).
Significance
The goal theorem is the paper's clearest quantitative instance of the separation of scales. The same function (8.1) shows a finite-rate exponential onset only in string time, and a finite-rate exponential saturation only in cosmic time. Neither regime survives in the other's units. Milestones 6–7 are the quantitative input to the paper's claim, attributed to H. Lin, that the cosmic Hamiltonian is not in the cosmic algebra while the string Hamiltonian passes the corresponding fluctuation test. Milestones 4–5 make precise the emergence of the string scale and the confinement of fermions near the horizon.
On formalization: the paper derives these limits by informal asymptotics, and some are stated only up to "∼". The formal statements fix the exact constants. Two of them differ from the paper's displayed formulas: the factor e−λ/2 in (11.1)–(11.2), which comes from (qN)q!/Nq→e−λ/2 when q2/N→λ, and the factor 2 in the exponent of (9.4). The paper's qualitative conclusions are unchanged. At the time of drafting, a search of the platform found no existing formal statements of these results.
Difficulty
Every statement is a limit in which two large parameters, N and q∼λN, compete inside an exponent. The naive move of expanding (1+x)−1/(q−1) to first order is valid only when x→0. In the goal's late-time regime x=Nqe(q−1)Jtc diverges, and at the end of the fast window x→λ stays finite. Each regime needs its own uniform control of log(1+x)/(q−1). The Hamiltonian milestones need the product asymptotics ∏k<q(1−k/N)→e−λ/2 under q2/N→λ, which is a second-order effect that first-order Stirling-type estimates miss.
Formalization scope
Nn,qn are natural-number sequences. The double-scaled limit is the predicate λ>0∧Nn→∞∧qn2/Nn→λ (real division). There is no parity or q≤N assumption.
PN,q uses a real power with base ≥1. The degenerate values at q∈{0,1} or N=0 affect only finitely many terms.
⟨Hc2⟩ and ⟨Hs2⟩ are defined by their closed forms. The operator-level identity (Wick contraction in the maximally mixed state) is not part of this mission. The integral in milestone 4 is the Lebesgue integral over R, and its nonzero value rules out the "non-integrable, hence 0" reading.
All statements are ordinary limits of real sequences. None is vacuous: the double-scaled predicate is satisfiable, e.g. by Nn=n, qn=⌊λn⌋.
Needed infrastructure: elementary real asymptotics (Mathlib's Filter.Tendsto, Real.rpow, Real.log, Real.cosh) and an improper integral of sech2. Reusable lemmas about (qN)q!/Nq in the regime q2/N→λ are welcome.
Selected references
L. Susskind, De Sitter Space, Double-Scaled SYK, and the Separation of Scales in the Semiclassical Limit, 2022. arXiv:2209.09999
L. Susskind, Scrambling in Double-Scaled SYK and De Sitter Space, 2022. arXiv:2205.00315
H. Lin, L. Susskind, Infinite Temperature's Not So Hot, 2022. arXiv:2206.01083
J. Maldacena, D. Stanford, Remarks on the Sachdev-Ye-Kitaev model, Phys. Rev. D 94 (2016) 106002. arXiv:1604.07818
D. A. Roberts, D. Stanford, A. Streicher, Operator growth in the SYK model, JHEP 06 (2018) 122. arXiv:1802.02633
V. Chandrasekaran, R. Longo, G. Penington, E. Witten, An Algebra of Observables for de Sitter Space, 2022. arXiv:2206.10780
String Theory (Tong) III: Modular Invariance of the One-Loop Partition FunctionTextbook
Motivation
At one loop, the worldsheet of a closed string is a torus, and the one-loop vacuum amplitude of the bosonic string is an integral over all conformally inequivalent tori. Two tori with modular parameters τ and τ′ are equivalent exactly when τ′=(aτ+b)/(cτ+d) for an integer matrix of determinant one, so a consistent string amplitude must be built from SL(2,Z)-invariant data. In David Tong's Cambridge Part III lecture notes (Section 6.4), the partition function of the bosonic string is computed in lightcone gauge and written in terms of the Dedekind eta function; the notes then observe that the measure and the integrand are individually modular invariant. This invariance is the mechanism by which string theory avoids the ultraviolet divergences of point-particle loops: the integration region is a fundamental domain that never reaches Imτ→0 (Polchinski 1986).
Setting
Write H={τ∈C:Imτ>0}. A modular transformation is the map
τ⟼γτ=cτ+daτ+b,a,b,c,d∈Z,ad−bc=1,(6.19)
generated by T:τ↦τ+1 and S:τ↦−1/τ. With q=e2πiτ the Dedekind eta function is
η(τ)=q1/24n=1∏∞(1−qn)=e2πiτ/24n=1∏∞(1−e2πinτ),
convergent and non-vanishing on H. The partition function of a single free scalar field on the torus is, up to constant factors, (6.20)
where F={∣τ∣≥1,∣Reτ∣≤21} is the standard fundamental domain.
Formalization targets
Goal: the string integrand is modular invariant
For every τ∈H and all a,b,c,d∈Z with ad−bc=1,
F(cτ+daτ+b)=F(τ).
Milestones (Section 6.4, in the order of the notes)
S and T generate SL(2,Z).
Every τ∈H is mapped into F by some modular transformation.
Im(γτ)=Imτ/∣cτ+d∣2 and d(γτ)/dτ=1/(cτ+d)2 (so d2τ↦d2τ/∣cτ+d∣4).
The measure d2τ/(Imτ)2 on H is SL(2,Z)-invariant.
η(τ+1)=e2πi/24η(τ).
η(−1/τ)=−iτη(τ).
Zscalar is modular invariant.
Significance
The result. Modular invariance of the integrand is what makes the one-loop string amplitude well defined as an integral over the moduli space of tori, rather than over the whole upper half-plane: restricting to F is only consistent because nothing depends on the representative. It is the first instance of a consistency condition that, in the superstring, constrains the allowed gauge groups and spectra.
Formalizing it. The eta transformation laws are classical (Apostol 1990, Ch. 3); Mathlib contains the definition of the eta function and results on SL(2,Z) and its fundamental domain, but, to the mission author's knowledge, not the S-transformation law of η. The mission produces a self-contained formal statement of the physics claim "both the measure and the integrand are individually modular invariant" and asks for its proof; milestones 1–2 are expected to follow from existing Mathlib results, milestone 6 is the analytic core.
Difficulty
The T-law is immediate from the product. The S-law η(−1/τ)=−iτη(τ) is not visible from the product formula at all: it requires either Poisson summation / the Jacobi triple product for theta functions, or a Mellin-transform / Eisenstein-series argument (logη is an Eichler integral of E2), together with care about the branch of the square root. Passing from generators to all of SL(2,Z) requires an induction on words in S and T in which the multiplier ∣cτ+d∣ must be tracked.
Formalization scope
Points of H are represented as complex numbers τ with the hypothesis 0<Imτ; a modular transformation is given by four integers with ad−bc=1 and acts by TongString.modularAction a b c d τ = (aτ+b)/(cτ+d). The eta function is defined directly by Tong's product formula as TongString.dedekindEta, using Mathlib's unconditional infinite product tprod. The square root in the S-law is the principal branch, (−iτ)1/2, which is the one intended since Re(−iτ)>0. The invariant measure is formalized as the measure on C with density (Imτ)−2 with respect to Lebesgue measure restricted to H, and invariance as equality of push-forward measures. Constant prefactors that Tong explicitly drops (powers of α′, 4π2) are omitted. Because every statement is restricted to Imτ>0, where η has no zeros, the reciprocal 1/∣η∣ is never a division by zero, so no statement is vacuous.
The Mathematics of Water I: Poiseuille's LawTextbook
Motivation
How much fluid does a pipe carry for a given pressure drop? For slow, viscous flow the answer is Poiseuille's law, Q=πΔpd4/(128ηL). It is the standard first quantitative model of blood flow in vessels and appears in Lecture 15 ("The Mathematics of Water") of the Drexel biophysics course PHYS 461/561 (Fall 2011), where it is derived from the Navier–Stokes equation and then used to estimate flow speeds in capillaries. The fourth-power dependence on the diameter is what makes small changes in vessel radius matter so much physiologically.
Setting
A Newtonian fluid of viscosity η>0 flows steadily through a straight cylindrical pipe of diameter d>0 and length L>0 under a pressure drop Δp=p(0)−p(L). By cylindrical symmetry the velocity is axial, v=v(r)ez, where r is the distance to the axis. The force balance on a thin cylindrical shell (slide 14), integrated along the pipe (slide 15), gives the radial equation
r1drd(rdrdv)=−ηLΔp,0<r<2d,
with boundary conditions: v finite at the axis, and no slip, v(d/2)=0, at the wall. The flow rate is Q=∫0d/2v(r)2πrdr and the average velocity is ⟨v⟩=Q/(πd2/4).
Target
Milestones, in order:
(slide 15) Every solution of the radial equation on (0,R) has the form v(r)=−ηLΔp4r2+C1lnr+C2.
(slide 16) Boundedness at the axis and no slip force v(r)=4ηLΔp(4d2−r2) on (0,d/2].
(slide 16) For this profile, Q=128ηLπΔpd4.
(slides 16–17) For this profile, ⟨v⟩=32ηLΔpd2.
Goal (Poiseuille's law): for every pipe flow v as above,
Q=∫0d/2v(r)2πrdr=128ηLπΔpd4.
Significance
The result is classical and completely known; the value of the mission is a clean, machine-checked version of the full chain from the radial ODE and its boundary conditions to the flow-rate formula, rather than just the final integral. The ODE-uniqueness step (milestones 1–2), which uses boundedness at the axis to exclude the logarithmic solution, is the part lecture notes usually wave through.
Difficulty
Milestones 3 and 4 are calculus exercises. The real work is milestones 1–2: turning "integrate twice" into a rigorous statement on an open interval (a function with zero derivative on an interval is constant), and then showing that boundedness near r=0 forces C1=0 (because lnr→−∞) and that continuity plus no-slip at the wall fixes C2.
Formalization scope
All quantities are real. The ODE is imposed pointwise on the open interval (0,d/2) with explicit differentiability; "v(r=0)<∞" is encoded as boundedness of v on (0,d/2), and the wall condition as v(d/2)=0 plus continuity from inside. η,L,d are assumed positive in every theorem; Δp is an arbitrary real. The pressure field itself is not modelled: the statements start from the integrated radial equation of slide 15. Two typos on the slides are corrected: slide 15 writes the homogeneous term as C1/r2 (correct: C1lnr), and slide 16 prints ⟨v⟩=Δpd4/(128ηL) (correct, as used on slide 17: Δpd2/(32ηL)). All definitions are in the single definition file MathematicsOfWater_PipeFlow; all declarations use the namespace MathematicsOfWater.
L. D. Landau and E. M. Lifshitz, Fluid Mechanics, 2nd ed., Pergamon, 1987, §17 (Poiseuille flow).
6 thms3 active usersReviewed
🏆Completed
Captain: Lucas
Introduction to Neutrino Physics I: Three-Flavour Vacuum OscillationsTextbook
Motivation
Neutrinos are produced and detected in charged-current weak interactions as flavour eigenstatesνe,νμ,ντ, but they propagate as superpositions of mass eigenstatesν1,ν2,ν3. When the masses differ, the mass components acquire different quantum-mechanical phases and the flavour content of the beam changes periodically along its path: neutrino flavour oscillations, first proposed by Pontecorvo. Oscillations of atmospheric and solar neutrinos are the experimental evidence that neutrinos are massive and that leptons mix, and the vacuum oscillation formula is what experiments fit to extract mass-squared differences, mixing angles and the CP-violating phase.
The formulas are elementary quantum mechanics, but the standard three-flavour form requires careful bookkeeping of phases, of unitarity relations, and of the sign conventions for the CP-odd term. The source for this mission, Section 4 of P. Lipari's lecture notes Introduction to Neutrino Physics, derives the whole chain in a few pages (Eqs. (44)–(73)). The mission turns that chain into machine-checkable statements.
Setting
Fix n flavours and n mass eigenstates. The mixing matrixU is a complex n×n matrix, rows indexed by flavour α and columns by mass eigenstate j, with ∣να⟩=∑jUαj∗∣νj⟩; physically U is unitary. Let mj2 be the squared masses, L the distance travelled and E the neutrino energy. In the relativistic approximation (L≃t, Ej≃E+mj2/2E) the transition amplitude and transition probability are
The antineutrino probabilityP(νˉα→νˉβ) is the same expression with U replaced by its entrywise complex conjugate U∗. Write Δmjk2=mk2−mj2 and define the CP-even amplitudes and the CP/T-odd coefficients
For two flavours the mixing matrix is the rotation (cosθ−sinθsinθcosθ); for three flavours the standard parametrization is U=R23U13(δ)R12 with mixing angles θ12,θ13,θ23 and CP phase δ (Eq. (57)).
Formalization targets
Goal: three-flavour transition probability (Eq. (71))
p. 135: no CP violation in survival probabilities, P(νˉα→νˉα)=P(να→να).
Significance
The goal is the formula used to interpret long-baseline accelerator and reactor oscillation data: the three amplitudes Aαβjk govern the CP-conserving oscillation pattern, and the single product term isolates all CP- and T-violating effects in vacuum. Together with milestones 3–6 it shows that one real number, the Jarlskog parameter, controls CP violation, and that such violation requires a nontrivial phase, three nonzero mixing angles and three active oscillation frequencies.
These results are textbook material and their proofs are known; what this mission adds is a machine-checked derivation with every convention explicit. Two points in the printed source are resolved by the formal statements: the constant term of Eq. (61) is printed with exponent 4 where squaring the amplitude gives exponent 2, and with the source's own Eqs. (57) and (69), the value of J12eμ has the opposite overall sign from the one printed in Eq. (70). Milestone 6 is stated with the sign that follows from (57) and (69).
Difficulty
No step is conceptually deep. The work is in the algebra: expanding a squared modulus of a sum of three complex exponentials, and then using unitarity to remove the constant term and to identify the sine coefficients with one another. The unitarity relations among the products UαjUβj∗ (footnote 8 of the source) have to be combined in the right order, and it is easy to lose a sign or a conjugate. Eq. (70) needs an explicit expansion of a product of three 3×3 complex matrices.
Formalization scope
The mixing matrix is Matrix (Fin n) (Fin n) ℂ; flavours e,μ,τ and mass labels 1,2,3 are the indices 0, 1, 2. Unitarity is membership in Matrix.unitaryGroup (Fin 3) ℂ, assumed only where the source uses it (Eqs. (67), (68), (71)).
Squared masses are an arbitrary function m2 : Fin n → ℝ; L and E are arbitrary reals. No positivity of E is imposed: all statements remain true at E=0, where Lean's convention x/0=0 makes every phase vanish.
Probabilities are Complex.normSq of the amplitude. Antineutrino probabilities use the entrywise conjugate U.map star, not the conjugate transpose.
The approximate formulas of §4.4 (one-mass-scale dominance, averaged oscillations) and the unit conversions of Eqs. (53)–(54) are out of scope.
Definitions live in two files: the oscillation objects (amplitude, probability, antineutrino probability, Aαβjk, Jjkαβ) and the mixing matrices (two-flavour rotation, standard PMNS parametrization). They are reusable for later missions on matter effects (Section 6 of the source).
Selected references
P. Lipari, Introduction to Neutrino Physics, lecture notes; Section 4 "Neutrino oscillations in vacuum", pp. 130–137 of the proceedings pagination (source of every statement in this mission).
B. Pontecorvo, original proposal of neutrino oscillations (reference [1] of the source).
Weinberg 1965: Infrared Photons and Gravitons — the soft exponents A and BResearch Paper
Motivation
When a charged particle scatters, it can radiate photons of arbitrarily low energy. In perturbative quantum electrodynamics this shows up as infrared divergences: loop integrals over very soft virtual photons diverge logarithmically, and so do the rates for emitting soft real photons. The standard resolution, due to Bloch and Nordsieck and systematized by Yennie, Frautschi and Suura, is that the divergences cancel in any rate that sums over undetectable soft emissions below an energy resolution E.
S. Weinberg's paper Infrared Photons and Gravitons (Phys. Rev. 140, B516 (1965)) shows that the same mechanism works for gravitons, derives explicit formulas for the exponents that control the soft factors, and proves that the extra "collinear" divergences produced by massless hard particles cancel for gravitons but not for photons. The soft-graviton factor (2.5)–(2.7) used there is the factor that appears in Weinberg's soft graviton theorem, and the spectrum formula is used in the paper to estimate the gravitational radiation from thermal collisions in the sun.
Timeline.
1937 — Bloch and Nordsieck: soft-photon divergences cancel between virtual and real emission.
1961 — Yennie, Frautschi and Suura: all-orders exponentiation of soft-photon factors in QED.
1965 — Weinberg: the analogous treatment for gravitons, closed forms (2.16) and (2.26) for the exponents A and B, cancellation of collinear divergences for gravitons (Sec. III), and the Coulomb-phase result conjectured by Dalitz (Sec. V).
Setting
Work in units ℏ=c=1 with the paper's metric: for four-vectors p=(p,p0) and q=(q,q0), p⋅q=p⋅q−p0q0. A process α→β has finitely many external linesn; line n has mass mn, three-momentum pn∈R3, on-shell energy En=∣pn∣2+mn2, charge en, and sign ηn=+1 (outgoing) or −1 (incoming). The relative velocity of lines n,m is
βnm=[1−(pn⋅pm)2mn2mm2]1/2.
For a unit vector q^ (the direction of a soft quantum) define the photon and graviton angular functions
and the infrared exponentsA=∫d2ΩA(q^) and B=∫d2ΩB(q^), integrals over the unit sphere against surface measure. With infrared cutoff λ and dividing point Λ, virtual soft quanta multiply a rate by (λ/Λ)A or (λ/Λ)B, and summing over real soft emission of total energy at most E gives the cutoff-free factors (E/Λ)Ab(A) and (E/Λ)Bb(B) of Eqs. (2.51)–(2.52).
Eq. (2.24) — B≥0 under energy–momentum conservation.
Eqs. (3.4)–(3.5) — if one line becomes massless, B stays finite (the lnm1 terms cancel by energy–momentum conservation).
Eq. (3.3) — for photons the corresponding limit is A→+∞ when the massless line is charged.
Eq. (4.9) — f(β)=1+611β2+4063β4+O(β6).
Significance
The exponents A and B are the only process-dependent quantities in the soft factors: once they are known in closed form, the infrared behaviour of any rate, the shape EA or EB of the soft spectrum, and the soft power spectrum EdΓ=BΓ0dE (2.53) follow. The Sec. III cancellation is the statement that gravitation, unlike massless electrodynamics or Yang–Mills theory, is free of unremovable collinear divergences at this level; the paper uses the contrast to argue that charged massless particles are problematic. The expansion (4.9) feeds the nonrelativistic quadrupole formula (4.13) behind the solar estimate.
The paper's arguments are informal physics derivations. This mission makes the mathematical content of those derivations precise and machine-checkable: the angular integrals, the sign statements, the limiting behaviour as a mass goes to zero, and the combinatorial factorization.
Difficulty
The closed forms need the solid-angle integral of a product of two denominators [En−pn⋅q^]−1[Em−pm⋅q^]−1 for two arbitrary, non-collinear momenta; the one-denominator case is a textbook integral in polar coordinates, but the two-denominator case does not reduce to it by an obvious choice of axis. The positivity statements are false without the conservation laws, so any proof must use them in an essential way. The massless-limit statements concern the interplay of a logarithmically divergent term with a coefficient that vanishes only at the limiting configuration, so a termwise limit fails.
Formalization scope
Three-momenta live in EuclideanSpace ℝ (Fin 3); energies are the on-shell values ∣p∣2+m2 and are not free variables. Solid-angle integrals are Bochner integrals on the unit sphere against Mathlib's volume.toSphere (total mass 4π). External lines are indexed by an arbitrary Fintype; the signs ηn are real numbers constrained to ±1 by hypothesis. Sums over n,m run over all ordered pairs, including n=m, where βnn=0; the kernels β−1ln1−β1+β and f are extended by their limits 2 and 1 at β=0, so the formulas are not trivialized by Lean's convention x/0=0. The printed Eq. (4.6) has an exponent 1/2 on the logarithm's argument that contradicts (2.26), (4.5) and (4.9); the definition of f drops it. The −iηϵ prescriptions are suppressed. The δ(q2) integration in (2.13) is carried out, leaving an integral over the shell λ≤∣q∣≤Λ. Positivity is stated as ≥0, because A=B=0 in degenerate configurations. Sec. IV beyond (4.9) (the velocity expansion (4.8), the quadrupole formula (4.13), the solar estimate) and the Coulomb phases of Sec. V are out of scope for this mission.
A complete development needs solid-angle integration on S2 (a spherical-coordinates change of variables for volume.toSphere), integration in polar coordinates on R3, and elementary inequalities for Minkowski products of timelike and null vectors. These are reusable beyond this mission, and contributions of them as separate lemmas are welcome.
D. R. Yennie, S. C. Frautschi and H. Suura, The infrared divergence phenomena and high-energy processes, Ann. Phys. 13, 379 (1961). https://doi.org/10.1016/0003-4916(61)90151-8
10 thms3 active usersReviewed
🏆Completed
Captain: Lucas
Monopoles, Instantons and Confinement I: Kink Solitons and the Bogomol'nyi BoundTextbook
Motivation
Solitons (localized, finite-energy, non-dissipating solutions of classical field equations) are the simplest examples of the topological objects that dominate the non-perturbative physics of gauge theories: vortices, magnetic monopoles and instantons. G. 't Hooft's lecture notes Monopoles, Instantons and Confinement (arXiv:hep-th/0010225) open with the two textbook examples in one space and one time dimension, the φ⁴ kink and the sine-Gordon soliton, because every later construction in the notes (the Abrikosov–Nielsen–Olesen vortex, the 't Hooft–Polyakov monopole, the BPST instanton) repeats the same pattern: a degenerate vacuum, a static solution interpolating between two vacua, and an energy (or action) bounded below by a topological quantity — the Bogomol'nyi bound — which the soliton saturates.
This mission formalizes that pattern in its cleanest setting, Chapter 1 (Sections 1.1–1.2) together with Exercise (i) of the notes.
Setting
A static real scalar field is a function φ:R→R. Its energy (eq. (1.5)) for a potential V is
EV[φ]=∫−∞∞(21(∂xφ)2+V(φ))dx∈[0,∞].
Two potentials are considered, with parameters λ,A,F>0:
Case (a), the "Mexican hat" (eq. (1.2)): Va(φ)=4!λ(φ2−F2)2, with the two vacua φ=±F and particle mass m given by m2=λF2/3.
Case (b), sine-Gordon (eq. (1.3)): Vb(φ)=A(1−cos(2πφ/F)), with vacua φ=nF, n∈Z, particle mass m=2πA/F and quartic coupling λ=16π4A/F4.
The static field equation is ∂x2φ=V′(φ). The notes give the explicit solutions
which interpolate between −F and F (case (a)) and between 0 and F (case (b)).
Formalization targets
Goal
For every x0: in case (a) the kink has energy 2m3/λ and no differentiable configuration with φ(−∞)=−F, φ(+∞)=F has smaller energy; in case (b) the soliton has energy 8m3/λ and no differentiable configuration with φ(−∞)=0, φ(+∞)=F has smaller energy:
The φ⁴ kink solves the static equation with the stated boundary values.
The sine-Gordon soliton solves the static equation with the stated boundary values.
The φ⁴ kink mass 2m3/λ (p. 5).
The sine-Gordon soliton mass 8m3/λ (p. 5).
Bogomol'nyi bound, case (a) (eq. (1.7), Exercise (i)).
Bogomol'nyi bound, case (b) (eq. (1.7), Exercise (i)).
Significance
The result identifies the soliton mass with a boundary term, E=W(φ)φ−∞φ+∞ with W′=2V, so the mass depends only on the pair of vacua joined, and shows the kink is a global energy minimizer in its topological sector — the classical stability statement behind treating the kink as a particle of mass ∝m3/λ, heavy in the weak-coupling regime.
The physics is classical and fully known. What this mission adds is a machine-checked account: explicit verification of the closed-form solutions, exact evaluation of improper energy integrals, and a rigorous Bogomol'nyi argument for configurations that are merely differentiable, with limits at ±∞ rather than compact support. No existing platform theorem matching these statements was found by a library search at drafting time.
Difficulty
The pointwise inequality 21φ′2+V(φ)≥φ′2V(φ) is elementary; the difficulty lies in turning it into a statement about improper integrals over R for arbitrary differentiable fields. One must handle configurations of infinite energy, pass from integrals over [a,b] to R using only the boundary limits, and deal with 2V being non-smooth at the vacua of Va and at every vacuum of Vb (a configuration may overshoot a vacuum). The explicit energy computations require exact evaluation of ∫sech4 and ∫sech2 type integrals as Lebesgue integrals.
Formalization scope
Fields are functions R→R; derivatives are Mathlib's deriv. Every competitor in the bounds is assumed differentiable everywhere, and boundary values are expressed as limits at −∞ and +∞.
The energy is a lower Lebesgue integral valued in [0,∞] of the (nonnegative) energy density, so configurations with infinite energy are allowed and satisfy the bounds trivially; no integrability hypothesis is imposed.
The parameters are strictly positive: λ>0, A>0, F>0. The masses m and the sine-Gordon coupling are defined from λ,A,F by the relations of Section 1.1, not left as free parameters.
The statements do not assume the competitor solves the field equation, and do not assume finite energy — assuming either would weaken the bound.
All declarations live in the namespace MonopolesInstantonsConfinement, intended to be shared by later missions drawn from the same notes.
Selected references
G. 't Hooft, Monopoles, Instantons and Confinement, lecture notes (Saalburg 1999) written by F. Bruckmann, 2000. arXiv:hep-th/0010225
E. B. Bogomol'nyi, The stability of classical solutions, Sov. J. Nucl. Phys. 24 (1976) 449.
R. Rajaraman, Solitons and Instantons, North-Holland, 1982.
Missing baryons from stacked SZ filaments: the model and statistics of de Graaff et al. (2019)Research Paper
Motivation
Big-bang nucleosynthesis and the acoustic peaks of the cosmic microwave background fix the baryon density of the universe to about 5% of its total energy density, but in the low-redshift universe only some 10% of those baryons are seen in galaxies, with another 10% in the circumgalactic and intracluster medium. Simulations place the remainder — the missing baryons — in a diffuse warm-hot intergalactic medium (WHIM) at 105–107 K, spread along the filaments of the cosmic web. Absorption-line studies and X-ray emission from individual filaments probe the cool and the hot ends of that range respectively, leaving the bulk poorly constrained.
de Graaff, Cai, Heymans & Peacock (2019) attack the problem statistically, by stacking the Planck Compton-y map over 1020334 pairs of CMASS galaxies and measuring the residual signal in the bridge between each pair. Interpreting that residual requires an analytic model of a filament, and that model — a cylinder with a Gaussian cross-section, convolved with the instrument beam — is what this mission formalizes. The same model, applied to the CMB lensing convergence map, converts the measured signals into a gas density and temperature, and hence into a baryon fraction.
Setting
A filament is modelled as an infinite cylinder seen side-on. In the plane perpendicular to the filament axis, write ℓ for the line-of-sight coordinate and r⊥ for the transverse coordinate on the sky. The electron density is taken to be a two-dimensional Gaussian
ne(ℓ,r⊥)=n0exp(−2σ2ℓ2)exp(−2(σ2+σB2)r⊥2),
with central densityn0, intrinsic widthσ and beam widthσB; the transverse direction, unlike the line of sight, is smoothed by the instrument beam, which is why σB appears only there.
The Compton y parameter of a line of sight is the optical-depth-weighted temperature of the scattering electrons,
y=mec2kBTeσT∫nedℓ,
with kB the Boltzmann constant, σT the Thomson cross-section, mec2 the electron rest energy and Te the electron temperature, assumed constant across the filament. The convergenceκ of CMB lensing is the analogous projection of the total matter density contrast δ=ρ/ρˉ−1,
κ=2c23H02Ωm∫0DSDSDL(DS−DL)aδdDL,
with H0 the Hubble constant, Ωm the matter density parameter, a the scale factor and DL,DS the comoving distances of lens and source. For a structure thin compared with the lensing kernel, the geometric factor is constant across it and κ reduces to a line-of-sight integral of δ against a fixed prefactor.
The significance of the measurement is assessed with a second, statistical layer. A stack of N binned profiles yik (k the map index, i the bin index) has mean profileyˉi=N1∑kyik and covariance estimatorCi,j=N1∑k(yik−yˉi)(yjk−yˉj); the deviation of the mean profile from zero is measured by χ2=∑i,jyˉi(C−1)i,jyˉj. Because the stacked maps overlap on the sky, the paper replaces C by a jackknife covariance built from Nsub sky sub-samples, Ci,jJK=NsubNsub−1∑k(yik−yˉi)(yjk−yˉj), and multiplies the inverse by the Hartlap factor(Nsub−n−2)/(Nsub−1), with n the number of bins.
Measurements enter through two numbers: the mean Compton parameter yˉ and the mean convergence κˉ over the boxed filament region, together with the empirical observation that the peak of each profile is close to 1/0.9 times its mean. Every statement in this mission imposes that calibration as a hypothesis on the model.
This is the quantity the paper's baryon budget is computed from: given the measured yˉ, an assumed temperature Te and the beam width, it fixes the number of electrons in a filament of length L.
The four displayed identities are the entire inferential chain from two stacked maps to a baryon fraction. Eq. (A.2) says what the model predicts for the observable; Eq. (A.3) inverts it at the filament axis; Eq. (A.4) turns the inversion into a total electron count, from which the paper obtains gas at (5.5±2.9)ρˉb and T=(2.7±1.7)×106 K, accounting for 11±7% of the cosmic baryon budget; Eq. (A.5) supplies the independent lensing constraint that breaks the density–temperature degeneracy of the SZ measurement alone. The statistical milestones cover the other half of the analysis: the covariance estimators of Eqs. (3) and (5)-(6) that turn a stack into an error bar, the nonnegativity of the χ2 of Eq. (4) that is converted into the quoted 2.9σ, and the Hartlap rescaling of Eq. (7). The beam-dominance milestone quantifies the paper's claim that the result barely depends on the one free shape parameter, the intrinsic width σ, which is taken from simulations rather than measured.
Formalizing them contributes a machine-checked derivation layer for a widely used observational technique: the projection of a Gaussian cylinder onto a Compton-y or convergence map, and the inversion of that projection, recur throughout stacked-SZ and stacked-lensing analyses. What this mission adds over the paper is a statement of each identity with all its hypotheses exposed — which positivity conditions are needed, which parameters genuinely drop out — rather than a new physical result. The physics is not re-derived and the measurement is not re-analysed; the empirical inputs yˉ, κˉ and the peak-to-mean ratio 1/0.9 enter as hypotheses, not as claims.
Difficulty
The identities are Gaussian integrals, so the mathematical depth is modest; the difficulty is bookkeeping and faithfulness. Three specific traps. First, the line-of-sight and transverse directions have different widths — σ against σ2+σB2 — and the asymmetry is what makes the final answer depend on the beam; a formalization that symmetrizes them is a different theorem. Second, the paper's n0 is not free: it is pinned by the calibration ypeak=yˉ/0.9, so the goal must carry that equation as a hypothesis rather than substituting a closed form for n0 silently. Third, the total electron count is an iterated integral over the plane, and the inner and outer variables must be kept in the paper's order.
Formalization scope
Everything is stated over R with no units, so every physical constant is an explicit real variable; integrals are Bochner integrals over the whole real line with Lebesgue measure, which return 0 on a non-integrable integrand. The definition bundle fixes the model once — the SZ prefactor, the Gaussian cylinder profile, the Compton parameter, the thin-lens convergence prefactor, the convergence and the total electron content — and every statement is phrased against it. The convergence definition is the thin-lens specialization of Eq. (2): the geometric kernel is a constant prefactor, not an integral over lens distance, and the mission does not claim the reduction from the full Eq. (2) to that form.
Degenerate readings are ruled out as follows. The calibration hypotheses are satisfiable for every admissible choice of constants, so no statement is vacuous; conversely they genuinely constrain the central amplitude, so no statement holds for a trivial reason. The positivity hypotheses are the minimal ones that make the denominators nonzero, and the widths enter as σ>0 with σB unrestricted, so the beam-free case σB=0 is included rather than excluded. The statistical layer is finite-dimensional: profiles are real arrays indexed by a finite bin type, covariances are real matrices, and matrix inversion is the nonsingular inverse, which returns the zero matrix on a singular argument — hence the positive-definiteness hypothesis in the χ2 milestone. Counts are natural numbers but all arithmetic in the prefactors is performed after casting to the reals, so Nsub−n−2 is a genuine real subtraction and may be negative. The statistical statements assert the structural properties of the estimators as written; they do not assert unbiasedness, nor the reduction C→C/N for the covariance of a mean, nor any property of the resampling scheme. Only Mathlib's Gaussian-integral, measure-theory and positive-semidefinite-matrix API is needed; no new infrastructure is required, and the resulting projection lemmas are reusable for any other stacked-profile analysis.
Selected references
A. de Graaff, Y.-C. Cai, C. Heymans, J. A. Peacock, Probing the missing baryons with the Sunyaev-Zel'dovich effect from filaments, Astronomy & Astrophysics 624, A48 (2019). doi:10.1051/0004-6361/201935159
R. A. Sunyaev, Y. B. Zeldovich, The Observations of Relic Radiation as a Test of the Nature of X-Ray Radiation from the Clusters of Galaxies, Comments on Astrophysics and Space Physics 4, 173 (1972).
Planck Collaboration XXII, Planck 2015 results. XXII. A map of the thermal Sunyaev-Zeldovich effect, A&A 594, A22 (2016). doi:10.1051/0004-6361/201525826
J. Hartlap, P. Simon, P. Schneider, Why your model parameter confidences might be too optimistic: unbiased estimation of the inverse covariance matrix, A&A 464, 399 (2007). doi:10.1051/0004-6361:20066170
H. Tanimura et al., A search for warm/hot gas filaments between pairs of SDSS Luminous Red Galaxies, MNRAS 483, 223 (2019). doi:10.1093/mnras/sty3118
11 thms3 active usersReviewed
🏆Completed
Captain: Lucas
Rydberg constant: Bohr-model derivation of the Rydberg formulaTextbook
Motivation
The Rydberg constantR∞ is the constant of spectroscopy that sets the scale of the hydrogen spectrum. It first appeared as an empirical fitting parameter in J. Rydberg's formula for the hydrogen spectral series; N. Bohr later showed that its value can be computed from more fundamental constants within his model of the atom. Before the 2019 revision of the SI it was, together with the electron spin g-factor, among the most accurately measured physical constants (relative standard uncertainty about 1.1×10−12, CODATA 2022). This mission formalizes the algebraic content of the standard account of R∞ as presented in the Wikipedia article Rydberg constant: its closed form, its reduced-mass correction, its alternative expressions through the fine-structure constant, and the Bohr-model derivation of the Rydberg formula.
Setting
Fix five strictly positive real numbers: the electron rest mass me, the elementary charge e, the vacuum permittivity ε0, the Planck constant h and the speed of light c. From them define
R∞=8ε02h3cmee4,ℏ=2πh,α=4πε01ℏce2,
the Compton wavelengthλe=h/(mec), Compton frequencyfC=mec2/h, Compton angular frequencyωC=2πfC, Bohr radiusa0=4πε0ℏ2/(e2me), classical electron radiusre=4πε01mec2e2, and the Rydberg unit of energyRy=hcR∞.
For a nucleus of mass M>0 the reduced mass is μ=1/(1/me+1/M) and the corrected Rydberg constant is RM=(μ/me)R∞; for hydrogen M=mp and RM=RH.
A Bohr orbit with principal quantum number n of a particle of mass m and charge −e about a fixed charge +e is a circular orbit of radius r>0 and speed v>0 such that the Coulomb force supplies the centripetal force and the angular momentum is quantized:
rmv2=4πε0r2e2,mvr=nℏ.
Its energy is E=21mv2−4πε0re2. The infinite-nuclear-mass model uses m=me; the reduced-mass model uses m=μ.
Formalization targets
Goal: the Rydberg formula with reduced mass
For distinct positive integers n1,n2 and Bohr orbits of mass μ with quantum numbers n1,n2, the wavenumber 1/λ=(En2−En1)/(hc) of the photon emitted in the transition satisfies
λ1=RM(n121−n221).
Milestones
Equivalent forms of the reduced-mass correction: RM=meμR∞=1+me/MR∞=me+MMR∞.
Isotopic shift: RM is strictly increasing in M and stays below R∞.
Rydberg unit of energy: hcR∞=α2mec2/2.
Alternative expressions: R∞=2hα2mec=2λeα2=4πa0α, and 1/R∞=(4π/α)a0.
Bohr energy levels (infinite nuclear mass): En=−hcR∞/n2.
Rydberg formula (infinite nuclear mass): λ1=Ry⋅hc1(n121−n221)=8ε02h3cmee4(n121−n221).
Bohr energy levels with reduced mass: En=−hcRM/n2.
Significance
These identities are the standard bridge between the empirical Rydberg formula and the constants me,e,ε0,h,c: they explain why a single constant governs all hydrogen series, why isotopes such as deuterium show shifted lines (the shift led to the discovery of deuterium), and how R∞ relates to α, a0 and the Compton scales. The results are classical and not open; the mission's contribution is a machine-checked, reusable layer of definitions (Bohr orbits, reduced mass, the named atomic length and energy scales) on which later atomic-physics formalizations can build.
Difficulty
The content is elementary algebra with positivity side conditions. The main point to handle carefully is that the Bohr-orbit conditions determine v and r only implicitly; the energy must be computed from the two defining equations rather than from explicit formulas, and every division must be justified by positivity of the constants.
Formalization scope
All constants are real numbers bundled in a structure with strict positivity fields; no SI numerical values (CODATA figures, 1.09678×107m−1, etc.) are formalized, and the precision-measurement and QED discussion of the source is out of scope. The Bohr model is formalized for a single particle of mass m orbiting a fixed charge +e (hydrogen-like, Z=1); the reduced-mass correction is modelled, as in the source, by substituting μ for me. The photon wavenumber is taken to be the orbital energy difference divided by hc. The Bohr-orbit hypotheses are satisfiable for every n≥1 and every m>0, so the Bohr-model statements are not vacuous. Contributions of alternative proofs and of generalizations (nuclear charge Z, explicit orbit radii and speeds) are welcome.
B. H. Bransden, C. J. Joachain, Quantum Mechanics, 2nd ed., Prentice Hall, 2000.
10 thms3 active usersReviewed
🏆Completed
Captain: Lucas
Elementary Charge: the SI defining relations and charge quantizationTextbook
Motivation
Since the 2019 revision of the International System of Units (SI), the elementary chargee is no longer a measured quantity: it is one of the seven defining constants of the SI and is fixed by definition at
e=1.602176634×10−19C.
Before that revision, e had to be extracted from experiment, and the metrology literature accumulated several independent routes to it: Faraday's laws of electrolysis combined with the Avogadro constant, Millikan and Fletcher's oil-drop experiment (1909), shot-noise analysis, and — since the 1980s — the combination of the Josephson effect and the quantum Hall effect. Each route is a short algebraic recipe that turns other measured constants into e. The 2019 redefinition did not delete those recipes; it inverted their role. With e, the Planck constant h, the Avogadro constant NA and the speed of light c all fixed exactly, the recipes become exact arithmetic identities between defining constants, and the residual measurement uncertainty migrates to the constants that are no longer fixed (the vacuum magnetic permeability μ0, equivalently the fine-structure constant α).
This mission formalizes exactly those identities, as they are stated in the source article, together with the two quantization statements the same article records: Dirac's 1931 argument that a magnetic monopole forces charge quantization, and the closure property behind the observation that although quarks carry charges in multiples of e/3, isolatable particles carry integer multiples of e.
Setting
All quantities are real numbers carrying SI units implicitly; the Lean development works in ℝ and fixes the numerical values as exact rationals, never as floating-point approximations.
The defining constants used here are e=1.602176634×10−19 (coulomb), NA=6.02214076×1023 (per mole), h=6.62607015×10−34 (joule second) and c=299792458 (metre per second).
From these the article builds four derived constants, each of which is a definition in this mission rather than an axiom:
the Faraday constantF=NAe, the charge of one mole of electrons;
the Josephson constantKJ=2e/h, measurable through the Josephson effect;
the von Klitzing constantRK=h/e2, measurable through the quantum Hall effect;
the fine-structure constant in the form used by CODATA, α=μ0ce2/(2h), with μ0 the vacuum magnetic permeability.
A fifth definition captures the natural unit of chargeq0=4πε0ℏc of those natural unit systems in which e=q0α, with ε0 the electric constant and ℏ the reduced Planck constant.
Target
The goal theorem asserts that the three determination routes recorded in the source return one and the same number, the fixed SI value of e, and that the electrolysis route's intermediate constant has its exact SI value:
The milestones are the individual identities, each stated for arbitrary positive values of the constants rather than only at the SI values, plus the two quantization statements:
F/NA=e for any NA=0;
the exact value of F at the SI values of NA and e;
2/(KJRK)=e for e=0, h=0;
the CODATA relation e=2hα/(μ0c);
e=q0α for q0=4πε0ℏc and α=e2/(4πε0ℏc);
Dirac quantization: if a monopole charge g=0 satisfies qg∈2ℏZ for every charge q in a collection, then every such q is an integer multiple of the fixed quantum ℏ/(2g);
if every charge in a collection S is an integer multiple of a quantum q0, then so is every charge in the additive subgroup of R generated by S.
Significance
The result itself. The content is the coherence of the SI's electrical sector: the same e comes out of an electrochemical measurement chain, a quantum-electrical one and the CODATA relation, and the arithmetic that links them is exact. Milestones 6 and 7 record the two quantization statements in the source: the Dirac argument is a genuine implication (a monopole forces the charge spectrum into a discrete subgroup), while milestone 7 isolates the algebraic half of "stable groupings of quarks carry integer charge" — closure of a quantized spectrum under composition.
Formalizing it. What this mission produces is a small, faithful, reusable encoding of the 2019 SI defining constants and of the standard derived constants, with the exact rational values rather than floating-point ones, and with the degenerate cases made explicit. Nothing here is an open problem, and the mission should not be read as one: the value is a checked reference layer for physics-flavoured formalization, not a research advance.
Difficulty
Low, and stated as such deliberately. Milestones 1–3 are field arithmetic; 4 and 5 are field arithmetic plus Real.sqrt of a square, which needs the positivity hypotheses that are carried explicitly; 6 and 7 are elementary manipulations of ℤ-multiples and of AddSubgroup.closure. The only places a solver can go wrong are the degenerate ones: division by a constant that has not been assumed nonzero, and Real.sqrt of a quantity not known to be nonnegative — which is why every statement carries positivity or nonvanishing hypotheses rather than leaving Lean's junk values to decide the outcome.
Formalization scope
Conventions this mission commits to:
Units are implicit; every constant is a bare real number in the SI unit named in its docstring.
The fixed values are exact rationals (for example eSI = 1602176634 / 10 ^ 28), so the numeric milestone is an exact identity, provable by norm_num, not a floating-point comparison.
The derived constants are functions of their arguments, so that each identity can be stated for general positive values and then instantiated at the SI values. This rules out a trivializing reading in which a "relation" holds only because both sides are the same closed numeral.
Every division carries a nonvanishing hypothesis and every Real.sqrt a nonnegativity or positivity hypothesis, so no statement is discharged by a Lean junk value.
Dirac's condition is stated as membership of qg in 2ℏZ for each charge in a given set, and its conclusion as membership of q in 2gℏZ; the physical derivation of the condition is not part of the mission.
No physics is assumed as an axiom: everything is a definition plus arithmetic over ℝ, using only Mathlib.
Contributions welcome beyond the current list: the shot-noise and oil-drop routes, which the source describes qualitatively rather than by a formula, would need a probabilistic or mechanical model before they can be stated, and that model would be a worthwhile addition.
R. A. Millikan, "The isolation of an ion, a precision measurement of its charge, and the correction of Stokes's law", Science 32 (822), 436–448 (1910). DOI: 10.1126/science.32.822.436
J. Preskill, "Magnetic Monopoles", Annual Review of Nuclear and Particle Science 34, 461–530 (1984). DOI: 10.1146/annurev.ns.34.120184.002333