Encyclopedia Cosmology Cosmology Fermion Weight Integral Fermi Integral Eq Weight Mul Bose
ARTICLE 3 claims 3 theorems
Cosmology Fermion Weight Integral Fermi Integral Eq Weight Mul Bose
A single number, 7/8, links the energy carried by matter particles to that carried by light in the early universe, and a machine-checked proof now ties that number to a purely mathematical identity.
The 7/8 factor in the cosmic ledger
In the hot, dense early universe, particles of matter and particles of light (radiation) both carry energy. The standard calculation of how much energy each species contributes uses two different integrals. For particles that obey Fermi-Dirac statistics (matter particles like electrons), the energy density integral is over the function t³/(eᵗ+1). For particles that obey Bose-Einstein statistics (light particles like photons), the integral is over t³/(eᵗ−1). The question is how these two integrals compare.
The answer, known in physics since the 1930s, is that the Fermi integral is exactly 7/8 of the Bose integral. This 7/8 factor is a real, measurable consequence of the different quantum statistics: fermions cannot share a quantum state, while bosons can pile into one, so fermions carry slightly less energy at the same temperature. The factor appears in every cosmology textbook as part of the effective number of relativistic degrees of freedom, g*, which sets the expansion rate of the early universe.
The Recognition Science framework's machine-checked library of formal theorems has now closed the mathematical gap between that textbook result and a purely arithmetic identity. The library proves, as a theorem, that the integral of t³/(eᵗ+1) from zero to infinity equals 7/8 times the integral of t³/(eᵗ−1) over the same range. The proof works by expressing each integrand as a geometric series in e⁻ᵗ, then evaluating the resulting sums using the Riemann zeta function ζ(4) and the Dirichlet eta function η(4). The identity η(4) = (7/8)ζ(4) is itself a known result, and the library's contribution is to connect that series identity to the thermodynamic integrals.
The theorem, named fermi_integral_eq_weight_mul_bose, is a statement about pure mathematics: it relates two definite integrals. It does not, by itself, make any physical claim. The physics enters through a separate, explicitly modeled choice: which particles are counted as relativistic in the early universe. That census, the g* bookkeeping, is a modeling decision, not a theorem. The 7/8 factor is the part that is now proved; the list of which particles contribute is not.
What the proof changes is the status of a familiar number. The 7/8 in the entropy-per-photon calculation is no longer an empirical coincidence or a textbook assertion. Within the framework, it is a derived consequence of the definitions of the integrals and the properties of the zeta and eta functions. The remaining freedom in the cosmology is not in the statistics factor, but in the physical inventory of the universe's contents.
THEOREM fermi_div_bose_integral · fermi_integral_eq_weight_mul_bose · IndisputableMonolith/Cosmology/FermionWeightIntegral.lean
/-- **THEOREM (7/8 at the thermodynamic layer).** The Fermi–Dirac energy
integral is exactly 7/8 of the Bose–Einstein one. -/
theorem fermi_div_bose_integral :
(∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t + 1))
/ (∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t - 1)) = 7 / 8 := by
rw [bose_integral_value, fermi_integral_value]
rw [div_eq_iff (by positivity)]
ring
/-- **THEOREM (fermion weight provenance, integral layer).** The `7/8` MODEL
constant of `EntropyPerPhoton.fermionWeight` is the ratio of the actual
thermodynamic integrals: `∫ t³/(eᵗ+1) = fermionWeight · ∫ t³/(eᵗ−1)`.
Together with `FermionWeight.fermionWeight_eq_eta_zeta_ratio` (series layer)
this makes the full mathematical content of the 7/8 factor THEOREM; the
remaining MODEL content of the entropy chain is the relativistic-species
census (`g*`), not the statistics factor. -/
theorem fermi_integral_eq_weight_mul_bose :
(∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t + 1))
= ((EntropyPerPhoton.fermionWeight : ℚ) : ℝ)
* ∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t - 1) := by
rw [bose_integral_value, fermi_integral_value]
unfold EntropyPerPhoton.fermionWeight
push_cast
ring
THEOREM fermi_integral_value · bose_integral_value · IndisputableMonolith/Cosmology/FermionWeightIntegral.lean
/-- **THEOREM (Fermi–Dirac integral).** `∫_{0}^{∞} t³/(eᵗ+1) dt = 7π⁴/120`. -/
theorem fermi_integral_value :
(∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t + 1)) = 7 * π ^ 4 / 120 := by
have h := mellin_fermi_value
rw [mellin_fermi_eq_integral] at h
exact Complex.ofReal_inj.mp h
/-- **THEOREM (Bose–Einstein integral).** `∫_{0}^{∞} t³/(eᵗ−1) dt = π⁴/15`. -/
theorem bose_integral_value :
(∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t - 1)) = π ^ 4 / 15 := by
have h := mellin_bose_value
rw [mellin_bose_eq_integral] at h
exact Complex.ofReal_inj.mp h
THEOREM hasSum_eta_shift · hasSum_zeta_shift · IndisputableMonolith/Cosmology/FermionWeightIntegral.lean
/-- `∑_{n≥0} (−1)ⁿ/(n+1)⁴ = η(4) = (7/8)·(π⁴/90)` (index-shifted
`FermionWeight.hasSum_eta_four`). -/
lemma hasSum_eta_shift :
HasSum (fun n : ℕ => (-1 : ℝ) ^ n / ((n : ℝ) + 1) ^ 4)
(7 / 8 * (π ^ 4 / 90)) := by
have hbase : HasSum (fun m : ℕ => (-1 : ℝ) ^ (m + 1) / (m : ℝ) ^ 4)
(7 / 8 * (π ^ 4 / 90)
+ ∑ i ∈ Finset.range 1, (-1 : ℝ) ^ (i + 1) / (i : ℝ) ^ 4) := by
simpa using FermionWeight.hasSum_eta_four
have h := (hasSum_nat_add_iff
(f := fun m : ℕ => (-1 : ℝ) ^ (m + 1) / (m : ℝ) ^ 4) 1).mpr hbase
exact h.congr_fun fun n => by push_cast [pow_succ]; ring
/-- `∑_{n≥0} 1/(n+1)⁴ = ζ(4) = π⁴/90` (index-shifted `hasSum_zeta_four`). -/
lemma hasSum_zeta_shift :
HasSum (fun n : ℕ => (1 : ℝ) / ((n : ℝ) + 1) ^ 4) (π ^ 4 / 90) := by
have hbase : HasSum (fun m : ℕ => (1 : ℝ) / (m : ℝ) ^ 4)
(π ^ 4 / 90 + ∑ i ∈ Finset.range 1, (1 : ℝ) / (i : ℝ) ^ 4) := by
simpa using hasSum_zeta_four
have h := (hasSum_nat_add_iff
(f := fun m : ℕ => (1 : ℝ) / (m : ℝ) ^ 4) 1).mpr hbase
exact h.congr_fun fun n => by push_cast; ring
What this page does not claim
This theorem does not identify which particles are relativistic; that census is a separate modeling choice. This theorem does not derive the values of ζ(4) or η(4); it uses them as known inputs. This theorem does not make any statement about the physical expansion rate of the universe.
Verify this page
Every tagged claim above names its theorem. To check one yourself rather than trust this page, elaborate the source module with Lean 4 and audit its axiom basis:
$ lake env lean IndisputableMonolith/Cosmology/FermionWeightIntegral.lean
expected axiom basis: [propext, Classical.choice, Quot.sound] (the Lean kernel's standard three; no RS-specific axioms)
A page whose claims cannot be reproduced this way does not ship. In production, every anchor links to the exact declaration in the public source release, and this block carries the build receipt for the page itself.
Derived articles
This page is generated by a question-recursion engine: the questions its answers raise become the next pages. The current agenda, with open targets marked red:
- Which particle species are counted as relativistic in the early universe, and how is that census determined?
- What is the exact value of the effective number of relativistic degrees of freedom, g*, in the standard model?
- How does the 7/8 factor change the predicted primordial helium abundance compared to a universe with only bosons?
- Does the same 7/8 ratio appear for the number density or pressure integrals, not just the energy density?
MACHINE LAYER · GROUNDED CLAIM TABLE · CLICK TO EXPAND
THEOREM fermi_div_bose_integral · fermi_integral_eq_weight_mul_bose · IndisputableMonolith/Cosmology/FermionWeightIntegral.lean
/-- **THEOREM (7/8 at the thermodynamic layer).** The Fermi–Dirac energy integral is exactly 7/8 of the Bose–Einstein one. -/ theorem fermi_div_bose_integral : (∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t + 1)) / (∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t - 1)) = 7 / 8 := by rw [bose_integral_value, fermi_integral_value] rw [div_eq_iff (by positivity)] ring/-- **THEOREM (fermion weight provenance, integral layer).** The `7/8` MODEL constant of `EntropyPerPhoton.fermionWeight` is the ratio of the actual thermodynamic integrals: `∫ t³/(eᵗ+1) = fermionWeight · ∫ t³/(eᵗ−1)`. Together with `FermionWeight.fermionWeight_eq_eta_zeta_ratio` (series layer) this makes the full mathematical content of the 7/8 factor THEOREM; the remaining MODEL content of the entropy chain is the relativistic-species census (`g*`), not the statistics factor. -/ theorem fermi_integral_eq_weight_mul_bose : (∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t + 1)) = ((EntropyPerPhoton.fermionWeight : ℚ) : ℝ) * ∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t - 1) := by rw [bose_integral_value, fermi_integral_value] unfold EntropyPerPhoton.fermionWeight push_cast ringThe integral of t³/(eᵗ+1) from zero to infinity equals 7/8 times the integral of t³/(eᵗ−1) over the same range. fermi_div_bose_integral · fermi_integral_eq_weight_mul_bose · IndisputableMonolith/Cosmology/FermionWeightIntegral.leanTHEOREM fermi_integral_value · bose_integral_value · IndisputableMonolith/Cosmology/FermionWeightIntegral.lean
/-- **THEOREM (Fermi–Dirac integral).** `∫_{0}^{∞} t³/(eᵗ+1) dt = 7π⁴/120`. -/ theorem fermi_integral_value : (∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t + 1)) = 7 * π ^ 4 / 120 := by have h := mellin_fermi_value rw [mellin_fermi_eq_integral] at h exact Complex.ofReal_inj.mp h/-- **THEOREM (Bose–Einstein integral).** `∫_{0}^{∞} t³/(eᵗ−1) dt = π⁴/15`. -/ theorem bose_integral_value : (∫ t in Ioi (0 : ℝ), t ^ 3 / (Real.exp t - 1)) = π ^ 4 / 15 := by have h := mellin_bose_value rw [mellin_bose_eq_integral] at h exact Complex.ofReal_inj.mp hThe Fermi integral evaluates to 7π⁴/120 and the Bose integral evaluates to π⁴/15. fermi_integral_value · bose_integral_value · IndisputableMonolith/Cosmology/FermionWeightIntegral.leanTHEOREM hasSum_eta_shift · hasSum_zeta_shift · IndisputableMonolith/Cosmology/FermionWeightIntegral.lean
/-- `∑_{n≥0} (−1)ⁿ/(n+1)⁴ = η(4) = (7/8)·(π⁴/90)` (index-shifted `FermionWeight.hasSum_eta_four`). -/ lemma hasSum_eta_shift : HasSum (fun n : ℕ => (-1 : ℝ) ^ n / ((n : ℝ) + 1) ^ 4) (7 / 8 * (π ^ 4 / 90)) := by have hbase : HasSum (fun m : ℕ => (-1 : ℝ) ^ (m + 1) / (m : ℝ) ^ 4) (7 / 8 * (π ^ 4 / 90) + ∑ i ∈ Finset.range 1, (-1 : ℝ) ^ (i + 1) / (i : ℝ) ^ 4) := by simpa using FermionWeight.hasSum_eta_four have h := (hasSum_nat_add_iff (f := fun m : ℕ => (-1 : ℝ) ^ (m + 1) / (m : ℝ) ^ 4) 1).mpr hbase exact h.congr_fun fun n => by push_cast [pow_succ]; ring/-- `∑_{n≥0} 1/(n+1)⁴ = ζ(4) = π⁴/90` (index-shifted `hasSum_zeta_four`). -/ lemma hasSum_zeta_shift : HasSum (fun n : ℕ => (1 : ℝ) / ((n : ℝ) + 1) ^ 4) (π ^ 4 / 90) := by have hbase : HasSum (fun m : ℕ => (1 : ℝ) / (m : ℝ) ^ 4) (π ^ 4 / 90 + ∑ i ∈ Finset.range 1, (1 : ℝ) / (i : ℝ) ^ 4) := by simpa using hasSum_zeta_four have h := (hasSum_nat_add_iff (f := fun m : ℕ => (1 : ℝ) / (m : ℝ) ^ 4) 1).mpr hbase exact h.congr_fun fun n => by push_cast; ringThe proof uses the series identity η(4) = (7/8)ζ(4). hasSum_eta_shift · hasSum_zeta_shift · IndisputableMonolith/Cosmology/FermionWeightIntegral.lean