The Universe as a Curve in the Simplex
When someone asks why the universe expands the way it does, there are two possible answers.
The first is the traditional one: because the Friedmann equation says that the square of the Hubble rate is proportional to the energy density, and that density has components that dilute in different ways. Radiation dilutes as a⁻⁴, matter dilutes as a⁻³, and the cosmological constant does not dilute. All of this is true, correct, and operationally indispensable. But taken in isolation, it still feels more like bookkeeping than explanation.
The second answer is: the universe expands the way it does because its background history is a curve in probability space, and that curve obeys a law of selection. Sectors that dilute more slowly gain share; sectors that dilute more rapidly lose share. The speed of this competition is Fisher information.
This article shows how this second formulation reduces, with mathematical rigor, to the first, and why this profoundly changes the way we look at cosmology.
Part I — Discovering the Right Time Variable
The first choice of the program is not a dynamical equation. It is a clock.
When the universe is measured by the scale factor a, expansion is multiplicative: the universe grows by successive factors. This exponential behavior is numerically inconvenient and, more importantly, hides the additive structure inherent in the problem.
The correct choice is:
N := ln a.
Each unit of N is one e-fold. In N, multiplying the size of the universe by a constant means traveling a fixed distance. Cosmological time becomes an additive ruler.
The second step is to write
Z(N) := H²(N).
This looks like a mere renaming, but it is the decisive conceptual point. In statistical mechanics, a partition function is a sum of positive exponential weights. In flat FLRW cosmology with constant-equation-of-state sectors, the square of the Hubble rate has exactly this form:
H²(N) = Z(N) = ∑ᵢ Aᵢ e^(−λᵢN), Aᵢ > 0.
Each term is a cosmic sector. The amplitudes Aᵢ indicate how much of that sector is present; the rates
λᵢ := 3(1 + wᵢ)
indicate how fast it dilutes. For radiation, w = 1/3, hence λ = 4. For matter, w = 0, hence λ = 3. For a cosmological constant, w = −1, hence λ = 0.
The consequence is immediate: if H² is a partition function, then
Φ(N) := ln Z(N) = ln H²(N)
is the logarithmic potential of the system. In statistical language, Φ is the log-partition function. In cosmological language, it is simply the logarithm of H².
From this point onward, much of background cosmology becomes the geometry of this potential.
Part II — The Simplex and the Equation That Governs Everything
Dividing each exponential weight by the total sum, we obtain the density fractions:
Ωᵢ(N) = Aᵢ e^(−λᵢN) / ∑ⱼ Aⱼ e^(−λⱼN).
These fractions satisfy
∑ᵢ Ωᵢ(N) = 1,
and
Ωᵢ(N) > 0.
Therefore, Ω(N) lies in the probability simplex:
Ω(N) ∈ Δⁿ⁻¹.
At each instant, the universe is a point in the simplex. Cosmic history is a curve in that space.
The instantaneous mean dilution rate is
λ̄(N) := ∑ᵢ Ωᵢ(N) λᵢ.
It measures the effective dilution rate of the total cosmic content at that moment.
Differentiating Ωᵢ with respect to N, we obtain:
Ωᵢ′ = −Ωᵢ(λᵢ − λ̄).
This is the replicator equation. It appears in evolutionary dynamics, game theory, and natural selection. Here, it governs the competition among cosmic sectors.
The interpretation is direct: if λᵢ > λ̄, that sector dilutes faster than the average and its fraction decreases. If λᵢ < λ̄, that sector dilutes more slowly than the average and its fraction increases.
Cosmic expansion is natural selection acting on cosmic sectors.
This sentence is not a loose metaphor. It is the direct reading of the equation.
Ωᵢ′ = −Ωᵢ(λᵢ − λ̄).
Radiation diluted faster than matter, so it lost share. Matter dilutes faster than the cosmological constant, so it also loses share. The cosmological constant does not dilute, so in the ΛCDM model it dominates asymptotically.
Part III — Fisher: The Speed of History
The next question is: how fast does this selection occur? The answer is the classical Fisher information of the curve Ω(N):
F_C(N) := ∑ᵢ (Ωᵢ′)² / Ωᵢ.
Substituting the replicator equation,
Ωᵢ′ = −Ωᵢ(λᵢ − λ̄),
we obtain
F_C = ∑ᵢ Ωᵢ(λᵢ − λ̄)².
In other words,
F_C = Var_Ω(λ).
The Fisher information of the universe is the variance of the dilution rates.
When all sectors dilute at the same rate, there is no compositional competition and F_C = 0. When sectors with very different dilution rates coexist, F_C is large.
But the central identity is even stronger. Since
Φ′(N) = Z′/Z = −λ̄,
we have
Φ″ = −λ̄′.
And since
λ̄′ = −F_C,
it follows that
Φ″ = F_C.
Moreover, the deceleration parameter is
q := −1 − H′/H.
Since
H′/H = ½Φ′ = −½λ̄,
we obtain
q = −1 + ½λ̄.
Differentiating,
q′ = ½λ̄′ = −½F_C.
Therefore,
F_C = Var_Ω(λ) = Φ″ = −λ̄′ = −2q′.
This is the central Fisher–FLRW identity.
It states that one single quantity appears simultaneously as:
the variance of the dilution rates;
the curvature of ln H²;
the decline of the mean dilution rate;
the variation of the deceleration parameter.
The last equivalence is physically powerful. The parameter q is cosmographic: it measures whether expansion accelerates or decelerates. The identit
q′ = −½F_C
shows that the decline of q is directly controlled by the variance of the dilution rates.
In the geometric reading of the program, cosmographic measurements are not merely detecting some mysterious dark substance; they are reading the statistical dispersion of the dilution rates of the cosmic sectors.
Since F_C ≥ 0, it follows that
q′ ≤ 0.
The deceleration parameter is monotonically non-increasing in Regime A. Cosmic history is a history of decreasing mean dilution, guided by the selection of the sectors that dilute more slowly.
Part IV — The Trinity of Geometries
There are three natural ways to measure distance in the space of cosmic compositions.
The first is the Fisher–Rao metric:
ds²_FR = ∑ᵢ dΩᵢ² / Ωᵢ.
Along the cosmic trajectory,
(ds_FR/dN)² = F_C.
Therefore,
ds²_FR = F_C dN².
The second comes from the square-root embedding:
|ψ_N⟩ := ∑ᵢ √Ωᵢ(N) |i⟩.
Sinc
⟨ψ_N|ψ_N⟩ = ∑ᵢ Ωᵢ = 1,
this construction places each composition as a real pure state. The quantum Fisher information of this real curve is
F_Q_real = 4(⟨∂_Nψ|∂_Nψ⟩ − |⟨ψ|∂_Nψ⟩|²).
Since
⟨ψ|∂_Nψ⟩ = ½∑ᵢ Ωᵢ′ = 0,
it follows that
F_Q_real = F_C.
The third is the Fubini–Study metric of the projective embedding:
ds²_FS = ¼F_C dN².
Therefore,
F_C dN² = ds²_FR = 4ds²_FS.
Classical Fisher information, Fisher–Rao geometry, and Fubini–Study geometry are the same geometry seen in three languages.
There is also a simple global representation. Define
xᵢ := 2√Ωᵢ.
Then
∑ᵢ xᵢ² = 4.
Thus the Fisher–Rao simplex is isometric to the positive orthant of a sphere of radius 2:
Δⁿ⁻¹_FR ≃ Sⁿ⁻¹_R=2,+.
Cosmic history is not an abstraction: it is a curve on a sphere.
For a complete binary transition, the Fisher–Rao length is
L_FR^binary = ∫_−∞^+∞ √F_C dN = π.
Every complete binary transition travels exactly π in this geometry. Radiation–matter, matter–Λ, or any idealized binary transition between two pure sectors describes half of a great circle in the Fisher–Rao simplex.
Part V — The Cumulants and the Hidden Clock
Once we recognize that
Φ(N) = ln Z(N),
a statistical hierarchy emerges automatically.
Define the cumulant generator of the dilution rates:
K(t; N) := ln ∑ᵢ Ωᵢ(N)e^(tλᵢ).
But
∑ᵢ Ωᵢ(N)e^(tλᵢ) = Z(N − t) / Z(N),
therefore
K(t; N) = Φ(N − t) − Φ(N).
The cumulants are
κᵣ(N) := ∂ʳK/∂tʳ |_{t=0}.
Therefore,
κᵣ(N) = (−1)ʳ Φ⁽ʳ⁾(N).
In particular,
κ₁ = λ̄,
κ₂ = F_C,
κ₃ = −Φ‴,
κ₄ = Φ⁗.
The entire chain obeys
κᵣ′ = −κᵣ₊₁.
Thus,
F_C′ = −κ₃.
The skewness of the dilution distribution decides whether Fisher information increases or decreases. The kurtosis determines the curvature of Fisher information. Higher cosmography is the cumulant hierarchy of cosmic composition.
Measuring q′, q″, jerk, snap, and higher derivatives is, in this language, measuring cumulants of the dilution-rate distribution.
There is no deep separation between cosmography and statistics: cosmography is compositional statistics written in spacetime.
Part VI — Distinguishing Epochs: The Triangle Theorem
How should we measure the difference between two cosmic epochs?
The natural distance between the compositions Ω(N) and Ω(M) is the Kullback–Leibler divergence:
D_KL[Ω(N) ‖ Ω(M)] := ∑ᵢ Ωᵢ(N) ln[Ωᵢ(N)/Ωᵢ(M)].
The ratio between the fractions at two epochs is
ln[Ωᵢ(N)/Ωᵢ(M)] = λᵢ(M − N) + Φ(M) − Φ(N).
Substituting into the Kullback–Leibler divergence,
D_KL[Ω(N) ‖ Ω(M)] = Φ(M) − Φ(N) − Φ′(N)(M − N).
This is exactly the Bregman divergence generated by Φ.
By the integral form of Taylor’s remainder,
Φ(M) = Φ(N) + Φ′(N)(M − N) + ∫_N^M (M − s)Φ″(s) ds.
Since
Φ″ = F_C,
it follows that
D_KL(N ‖ M) = ∫_N^M (M − s)F_C(s) ds.
Every oriented informational distance between two epochs is accumulated Fisher information with a triangular weight.
The reverse divergence is
D_KL(M ‖ N) = ∫_N^M (s − N)F_C(s) ds.
The sum of the two gives the Jeffreys divergence:
J(N, M) = D_KL(N ‖ M) + D_KL(M ‖ N) = (M − N)∫_N^M F_C(s) ds.
The universe becomes informationally distant from itself only when Fisher information is present, that is, when there is diversity in dilution rates. In an epoch dominated by a single sector, expansion may continue, but the composition barely changes.
Transition eras are the epochs that move the universe the most in statistical space.
Part VII — Crossing the Threshold
Cosmic acceleration begins when
q = 0.
But
q = −1 + ½λ̄.
Therefore,
q = 0 ⇔ λ̄ = 2.
Since
w_eff = −1 + ⅓λ̄,
we also have
λ̄ = 2 ⇔ w_eff = −⅓.
Therefore,
q = 0 ⇔ λ̄ = 2 ⇔ w_eff = −⅓.
The universe begins to accelerate exactly when the mean dilution rate crosses the value 2.
In the matter–Λ case, this means that the cosmological constant does not need to be dominant in order to trigger acceleration. Since matter has λ_m = 3 and Λ has λ_Λ = 0, the condition
λ̄ = 3Ω_m + 0·Ω_Λ = 2
implies
Ω_m = ⅔, Ω_Λ = ⅓.
Acceleration begins when Λ reaches one third of the total density, before equality Ω_m = Ω_Λ = ½.
Part VIII — Entering the Branch: Catalan and the Discriminant
We begin with the familiar universe: radiation, matter, and dark energy contributing positive terms to H². In this standard FLRW setting, cosmic composition traces a curve in a probability simplex. After this geometric structure is exposed, the program considers a second step: what happens if horizon information also gravitates as an effective contribution proportional to H⁴? This hypothesis leads to the algebraic extension
H² = H²_bg + αηD H⁴.
Define
y := H²/H²_bg,
and
ξ := αηD H²_bg.
The equation reduces to
ξy² − y + 1 = 0.
The physical branch, continuous with General Relativity when ξ → 0, is
y₋(ξ) = 2/[1 + √(1 − 4ξ)].
Equivalently,
y₋(ξ) = [1 − √(1 − 4ξ)]/(2ξ).
The reality of the branch requires
0 ≤ ξ ≤ ¼.
The informational fraction is
Ω_I = 1 − 1/y.
Since the quadratic equation implies
Ω_I = ξy,
we also obtain
ξ = Ω_I(1 − Ω_I).
The branch is called Catalan because
y₋(ξ) = ∑_{n=0}^∞ C_n ξⁿ, C_n = 1/(n+1) · binom(2n,n).
The radius of convergence of this series is
R_ξ = ¼,
exactly the edge of the discriminant.
Criticality becomes clear through the identity
½ − Ω_I = √(¼ − ξ).
The distance to saturation Ω_I = ½ is the square root of the distance to the discriminant ξ = ¼.
Moreover,
dΩ_I/dξ = 1/√(1 − 4ξ) = 1/[2√(¼ − ξ)].
Thus the susceptibility diverges with critical exponent ½. If Regime B described real physics, small variations in ξ near the boundary would produce large responses in H, Ω_I, and the associated cosmography.
Part IX — The Golden Point
The Catalan branch contains a special interior point.
If
φ := (1 + √5)/2,
then
φ² = φ + 1, φ − 1 = φ⁻¹.
On the branch,
ξ = (y − 1)/y².
Taking
y = φ,
we obtain
ξ = (φ − 1)/φ² = φ⁻¹/φ² = φ⁻³.
Therefore,
y = φ ⇔ ξ = φ⁻³.
At this point,
Ω_I = 1 − 1/y = 1 − φ⁻¹ = φ⁻²,
and
Ω_bg = 1/y = φ⁻¹.
Hence,
Ω_I = φ⁻², Ω_bg = φ⁻¹, Ω_bg/Ω_I = φ.
The golden ratio is not inserted into the model. It appears as a special interior point of the Catalan geometry, where the partition between the background and the informational sector obeys a self-similar relation.
Part X — Massieu Duality: Two Sides of the Same Curvature
Compositional thermodynamics appears when we define the weighted entropy
S_A(Ω) := −∑ᵢ Ωᵢ ln(Ωᵢ/Aᵢ).
Since
Ωᵢ = Aᵢ e^(−λᵢN)/Z,
we have
ln(Ωᵢ/Aᵢ) = −λᵢN − Φ.
Therefore,
S_A = −∑ᵢ Ωᵢ(−λᵢN − Φ) = Nλ̄ + Φ.
Hence,
S_A = Φ + Nλ̄.
This is the Massieu form of compositional cosmology.
Since
dS_A = N dλ̄,
and
dλ̄/dN = −F_C,
it follows that
dN/dλ̄ = −1/F_C.
Therefore,
d²S_A/dλ̄² = −1/F_C.
Fisher information appears as positive curvature in N:
Φ″ = F_C,
but as negative inverse curvature in the dual λ̄-space:
S_A″(λ̄) = −1/F_C.
The same quantity is, at once:
statistical variance;
convex curvature;
geometric speed;
inverse entropic rigidity.
This is the signature of a deep structure: the same mathematical object reappears in independent languages.
Epilogue: What the Program Really Says
The program is not merely a proposal for a cosmological correction. Its strongest core is an exact rewriting of positive FLRW cosmology as statistical geometry.
Cosmic expansion, read in the correct variables, is the selective flow of the least-diluting sectors in the simplex of density fractions. Fisher information measures the speed of that selection:
F_C = Var_Ω(λ) = Φ″ = −λ̄′ = −2q′.
Cosmography — H, q, jerk, snap, and higher-order derivatives — is the cumulant hierarchy of a distribution of dilution rates.
Fisher–Rao geometry turns the simplex into a sphere of radius 2. Aitchison geometry turns the evolution into a straight line of log-ratios. Massieu duality turns Φ and S_A into conjugate potentials. The Catalan branch turns the H⁴ modification into a discriminant algebra with square-root criticality.
The conditional part of the program — associated with Landauer, horizons, operational opacity, and H⁴ — must be judged by microphysics and data. But the algebraic core is independent of that: if H² is a positive sum of exponential modes, then the entire Fisher–FLRW structure follows.
The final statement can be written as:
The background universe is a curve in the simplex.
Expansion selects the sectors that dilute the least. Fisher information measures the speed of that selection.
Cosmography is the geometric shadow of this statistics.
At this level, the universe is not merely a solution of differential equations. It is a moving exponential family.