Abstract
The spatially flat Friedmann equation describes the expansion of the Universe as a sum of energetic contributions. In its simplest binary form, it can be written as
H² = H_B² + H_X²,
where B represents a background sector and X a complementary sector. This article shows that every positive decomposition of this type admits an exact reparametrization by the generating function of the Catalan numbers. Defining
y = H²/H_B²,
Ω_X = H_X²/H²,
ξ_X = Ω_X(1 − Ω_X),
one obtains identically
y = 1 + ξ_X y².
This is the same functional equation satisfied by the generating function of the Catalan numbers,
C(ξ) = 1 + ξ C(ξ)².
On the regular branch, corresponding to the case in which X is subdominant, namely
0 ≤ Ω_X ≤ 1/2,
it follows that
H²/H_B² = C[Ω_X(1 − Ω_X)].
The reparametrization provides a normal form for cosmological dominance transitions. The variable ξ_X measures the competition between sectors, its maximum value ξ_X = 1/4 marks the equality H_X² = H_B², and the series expansion is rigidly fixed by the Catalan numbers. Cosmological equal dominance therefore coincides with the branch point ξ = 1/4 of the Catalan generating function. The construction does not modify the dynamics of the ΛCDM model; it reorganizes the Friedmann equation in natural transition coordinates.
Keywords: FLRW cosmology; Friedmann equation; Catalan numbers; dominance transitions; ΛCDM; dark energy.

1. Introduction
One of the central goals of cosmology is to describe how the Universe expands. This expansion is measured by the Hubble parameter,
H(t) = ȧ(t)/a(t),
where a(t) is the scale factor.
In the standard cosmological model, called ΛCDM, the content of the Universe is described, to a good approximation, by three main sectors:
radiation, matter, and the cosmological constant Λ.
In a spatially flat FLRW geometry, the Friedmann equation is
H²(a) = H₀² [Ω_r0 a⁻⁴ + Ω_m0 a⁻³ + Ω_Λ0].
Each term dominates in a different era:
radiation dominates in the primordial Universe;
matter dominates at intermediate times;
Λ dominates at late times.
Therefore, cosmic history may be viewed as a sequence of dominance transitions:
radiation → matter → Λ.
Traditionally, these transitions are located by equalities such as
ρ_r = ρ_m,
ρ_m = ρ_Λ.
These equalities are correct, but they do not fully display the algebraic structure of the transition. The purpose of this article is to show that there is a more organized way of writing these passages.
Central idea: every positive decomposition of the Friedmann equation into two sectors admits a Catalan normal form.
In other words, the Friedmann equation does not merely add densities. It also allows a reparametrization that reveals the algebraic grammar of its transitions.

2. The Friedmann Equation in Two Sectors
Let us begin with a simple decomposition:
H² = H_B² + H_X².
Here,
B = background sector,
X = complementary sector.
For example, during the matter era, when one wishes to study the approach to Λ domination, one may take
B = m,
X = Λ.
Then
H_B² = H_m²,
H_X² = H_Λ².
In another regime, if radiation is still relevant, one may choose
B = m + r,
X = Λ.
The decomposition is only a choice of language. It does not change the Friedmann equation.
We now define the ratio between the total expansion and the background expansion:
y = H²/H_B².
Since
H² = H_B² + H_X²,
we have
y = (H_B² + H_X²)/H_B²,
that is,
y = 1 + H_X²/H_B².
Thus, y measures how much the total expansion differs from the expansion produced by the background sector alone.
If X is negligible, then H_X² ≪ H_B² and y ≈ 1.
If X becomes comparable to B, then y departs from 1.
If X strongly dominates, while B is kept fixed, then y can become very large.

3. The Fraction of the Complementary Sector
We define the fraction of sector X relative to the total:
Ω_X = H_X²/H².
This quantity is analogous to the usual density fractions in cosmology.
If X is small, then Ω_X ≈ 0.
If X dominates, then Ω_X ≈ 1.
From the definition,
H_X² = Ω_X H².
Since
H² = H_B² + H_X²,
we have
H_B² = H² − H_X².
Substituting H_X² = Ω_X H², we obtain
H_B² = H²(1 − Ω_X).
Therefore,
H²/H_B² = 1/(1 − Ω_X).
That is,
y = 1/(1 − Ω_X).
This relation will be the key to the reparametrization.

4. The Transition Variable
The fraction Ω_X measures how much of the total lies in sector X. But to study a dominance transition, we want to measure something else: we want to know when two sectors are competing.
For this purpose, we introduce
ξ_X = Ω_X(1 − Ω_X).
This variable has three simple properties.
First, if Ω_X ≈ 0, then
ξ_X ≈ 0.
In this case, X is negligible.
Second, if Ω_X ≈ 1, then again
ξ_X ≈ 0.
In this case, X clearly dominates.
Third, when
Ω_X = 1/2,
we have
ξ_X = (1/2)(1 − 1/2) = 1/4.
Hence,
0 ≤ ξ_X ≤ 1/4.
The variable ξ_X is maximal exactly when the two sectors have equal weight.
Ω_X measures abundance. ξ_X measures dominance competition.
This distinction is the physical motivation for the reparametrization.

5. Derivation of the Catalan Form
We have already seen that
y = 1/(1 − Ω_X).
Now compute ξ_X y².
Since
ξ_X = Ω_X(1 − Ω_X),
we have
ξ_X y² = Ω_X(1 − Ω_X) · 1/(1 − Ω_X)².
Canceling one factor of 1 − Ω_X, we obtain
ξ_X y² = Ω_X/(1 − Ω_X).
Then
1 + ξ_X y² = 1 + Ω_X/(1 − Ω_X).
Putting the terms over a common denominator,
1 + ξ_X y² = (1 − Ω_X)/(1 − Ω_X) + Ω_X/(1 − Ω_X),
therefore
1 + ξ_X y² = 1/(1 − Ω_X).
But
1/(1 − Ω_X) = y.
Hence,
y = 1 + ξ_X y².
This is the central equation of the article.
It is not an approximation. It is an algebraic identity derived directly from the positive decomposition
H² = H_B² + H_X².

6. Appearance of the Catalan Numbers
The Catalan numbers appear in many areas of mathematics. They count, for example, certain ways of parenthesizing expressions, binary trees, and lattice paths that do not cross a boundary.
The generating function of the Catalan numbers is defined by
C(ξ) = 1 + ξ C(ξ)².
Comparing this with the cosmological equation
y = 1 + ξ_X y²,
we see that the structure is the same.
Therefore, on the regular branch,
y = C(ξ_X).
Since
y = H²/H_B²,
we obtain
H²/H_B² = C(ξ_X) = C[Ω_X(1 − Ω_X)].
The explicit solution of the Catalan function is
C(ξ) = [1 − √(1 − 4ξ)]/(2ξ),
with continuous extension
C(0) = 1.
More precisely,
C(ξ) = [1 − √(1 − 4ξ)]/(2ξ), for 0 < ξ ≤ 1/4,
and
C(0) = 1.
The continuous extension C(0) = 1 must be understood on the regular branch, that is, in the limit Ω_X → 0. The same value ξ_X → 0 also occurs when Ω_X → 1, but in that case we are on the other branch of the quadratic equation, corresponding to the dominance of X over B.
The series expansion is
C(ξ) = 1 + ξ + 2ξ² + 5ξ³ + 14ξ⁴ + 42ξ⁵ + ⋯.
Thus,
H² = H_B² [1 + ξ_X + 2ξ_X² + 5ξ_X³ + 14ξ_X⁴ + 42ξ_X⁵ + ⋯].
The coefficients
1, 1, 2, 5, 14, 42, …
are the Catalan numbers.

7. Why Is This Useful?
A generic expansion could have the form
y = 1 + b₁ξ + b₂ξ² + b₃ξ³ + ⋯.
In that case, the coefficients b₁, b₂, b₃, … would be arbitrary.
But the Friedmann equation, when written in the variable ξ_X, fixes these coefficients:
b₁ = 1,
b₂ = 2,
b₃ = 5,
b₄ = 14,
and so on.
That is,
bₙ = Cₙ,
where Cₙ are the Catalan numbers.
The reparametrization turns the cosmological transition into a universal series, without arbitrary phenomenological coefficients.
This is the formal advantage of the construction.
The important point is that the universality does not come from a new physical hypothesis. It comes from the algebra of the positive decomposition
H² = H_B² + H_X²
together with the choice of the competition variable
ξ_X = Ω_X(1 − Ω_X).

8. The Radius of Convergence and Cosmological Equality
The Catalan generating function has a special point at
ξ = 1/4.
This value is the natural limit of the Catalan series on the positive real axis. In the cosmological reparametrization,
ξ_X = Ω_X(1 − Ω_X).
As we have seen,
ξ_X ≤ 1/4.
The maximum value occurs when
Ω_X = 1/2.
But
Ω_X = H_X²/H².
If
Ω_X = 1/2,
then
H_X² = H²/2.
Since
H² = H_B² + H_X²,
it follows that
H_X² = H_B².
Therefore,
ξ_X = 1/4 ⇔ H_X² = H_B².
The branch point ξ = 1/4 of the Catalan function coincides with cosmological equal dominance.
This is the main structural observation.
The dominance transition is encoded by the approach to the Catalan boundary; equal dominance occurs exactly at ξ_X = 1/4.
At the equality point,
C(1/4) = 2,
therefore
H²/H_B² = 2,
as expected, since
H² = H_B² + H_X² = 2H_B².

9. The Branch Structure
The equation
y = 1 + ξy²
can be written as
ξy² − y + 1 = 0.
This is a quadratic equation in y. For 0 < ξ ≤ 1/4, its solutions are
y_±(ξ) = [1 ± √(1 − 4ξ)]/(2ξ).
The regular branch is the branch with the minus sign:
y₋(ξ) = [1 − √(1 − 4ξ)]/(2ξ).
This branch satisfies
lim ξ→0 y₋(ξ) = 1.
Therefore,
y₋(ξ) = C(ξ).
The branch with the plus sign satisfies
y₊(ξ) = [1 + √(1 − 4ξ)]/(2ξ),
and diverges when ξ → 0:
y₊(ξ) ∼ 1/ξ.
Therefore, only y₋ is regular in the limit where X disappears.
To understand the physical meaning of these branches, define
r = H_X²/H_B².
Then
y = 1 + r.
Moreover,
Ω_X = r/(1 + r).
If
0 ≤ r ≤ 1,
then
0 ≤ Ω_X ≤ 1/2.
This means that X is subdominant. In this case, the correct branch is the regular Catalan branch:
H²/H_B² = C(ξ_X).
If, on the other hand,
H_X² > H_B²,
then X has become dominant. The same variable ξ_X corresponds to the other branch of the quadratic solution, unless the dominant sector is redefined as the new background B.
Practical rule: choose B as the dominant sector and X as the subdominant sector. In this way, each cosmological era remains on the regular Catalan branch.

10. Example: The Matter–Λ Transition
Let us now consider the simplest application.
In the late Universe, radiation may be neglected, and we may write
H²(z) = H₀² [Ω_m0(1 + z)³ + Ω_Λ0].
Before Λ domination, we choose
B = m,
X = Λ.
Then
H_B²(z) = H₀² Ω_m0(1 + z)³,
and
H_X² = H₀² Ω_Λ0.
The Λ fraction is
Ω_Λ(z) = Ω_Λ0/[Ω_m0(1 + z)³ + Ω_Λ0].
The Catalan variable of the transition is
ξ_Λ(z) = Ω_Λ(z)[1 − Ω_Λ(z)].
The maximum occurs when
Ω_Λ(z) = 1/2.
This is equivalent to
Ω_Λ0 = Ω_m0(1 + z)³.
Hence,
1 + z_eq,Λ = (Ω_Λ0/Ω_m0)¹ᐟ³.
This is the redshift at which matter and Λ have equal weight in the Friedmann equation.
Matter–Λ equality is the point at which ξ_Λ(z) reaches its maximum value 1/4.
At this point,
H_Λ² = H_m²,
H² = 2H_m²,
and therefore
H²/H_m² = 2 = C(1/4).

11. Equality Is Not the Same as Acceleration
It is important to distinguish two ideas:
matter–Λ equality
and
the onset of cosmic acceleration.
Matter–Λ equality occurs when
ρ_m = ρ_Λ.
The onset of acceleration, however, is determined by the equation
ä/a = −(4πG/3)(ρ + 3p).
For pressureless matter,
p_m = 0.
For the cosmological constant,
p_Λ = −ρ_Λ.
Then
ρ + 3p = ρ_m − 2ρ_Λ.
Acceleration begins when
ρ_m < 2ρ_Λ.
The boundary is
ρ_m = 2ρ_Λ.
Therefore,
acceleration begins when ρ_m = 2ρ_Λ,
whereas
the Catalan boundary occurs when ρ_m = ρ_Λ.
In terms of the Λ fraction, at the acceleration boundary we have
Ω_Λ = ρ_Λ/(ρ_m + ρ_Λ).
Since ρ_m = 2ρ_Λ,
Ω_Λ = ρ_Λ/(2ρ_Λ + ρ_Λ) = 1/3.
Hence,
ξ_Λ = Ω_Λ(1 − Ω_Λ) = (1/3)(2/3) = 2/9.
Thus,
acceleration: ξ_Λ = 2/9,
equal dominance: ξ_Λ = 1/4.
The reparametrization helps separate these two events.
Acceleration is a dynamical condition, because it depends on ρ + 3p.
The Catalan boundary is an algebraic condition of equality between sectors in the Friedmann equation.

12. The Signed Dominance Variable
The variable
ξ_X = Ω_X(1 − Ω_X)
is symmetric under
Ω_X ↔ 1 − Ω_X.
This means that it measures how close we are to equality, but it does not say which sector dominates.
For this purpose, we introduce the signed variable
s_X = 1 − 2Ω_X.
Then:
s_X > 0 ⇔ B dominates,
s_X = 0 ⇔ B and X are equal,
s_X < 0 ⇔ X dominates.
Moreover,
ξ_X = ¼(1 − s_X²).
Thus, the pair
(ξ_X, s_X)
provides a complete chart of the transition:
ξ_X = proximity to equality,
s_X = orientation of dominance.
On the B-dominant branch, we have s_X ≥ 0. Since
1 − Ω_X = (1 + s_X)/2,
it follows that
y = 1/(1 − Ω_X) = 2/(1 + s_X).
On the other hand,
√(1 − 4ξ_X) = s_X
on the branch s_X ≥ 0. Hence,
C(ξ_X) = 2/[1 + √(1 − 4ξ_X)] = 2/(1 + s_X).
This form shows that the signed variable regularizes the passage through equality.
Near equal dominance, ξ_X reaches a maximum and C′(ξ_X) becomes singular as a function of ξ_X. This singularity is a singularity of the unsigned chart ξ_X. The variable s_X distinguishes the two sides of the transition and removes the branch ambiguity.

13. Dynamical Dark Energy in Catalan Coordinates
In more general models, dark energy may have a redshift-dependent equation of state,
w = w(z).
In that case, its density evolves as
ρ_de(z) = ρ_de,0 exp{3 ∫₀ᶻ [(1 + w(z̃))/(1 + z̃)] dz̃}.
The dark-energy fraction is
Ω_de(z) = ρ_de(z)/ρ_tot(z).
The corresponding Catalan coordinate is
ξ_de(z) = Ω_de(z)[1 − Ω_de(z)].
This function answers a simple question:
When does dark energy compete most strongly with the material background?
In ΛCDM, the trajectory ξ_Λ(z) is fixed by the parameters Ω_m0 and Ω_Λ0.
In models with dynamical dark energy, the trajectory ξ_de(z) changes.
Thus, different models can be compared through the shape of the function
ξ_de(z).
This comparison does not replace observables such as H(z), luminosity distances, BAO, or structure growth. It only provides a clean algebraic coordinate for visualizing when and how the competition between sectors intensifies.

14. Structure Growth
The expansion of the Universe also affects the growth of structure.
In general relativity, in the linear regime and on subhorizon scales, the matter density contrast approximately obeys
δ_m″ + [2 + H′/H]δ_m′ − (3/2)Ω_m(a)δ_m = 0,
where
N = ln a,
and the prime denotes differentiation with respect to N.
Using
H² = H_B² C(ξ_X),
we obtain
ln H = ½ ln H_B² + ½ ln C(ξ_X).
Hence,
H′/H = ½[(H_B²)′/H_B²] + ½[C′(ξ_X)/C(ξ_X)] ξ_X′.
This expression shows that the Hubble friction has two parts:
one part associated with the evolution of the background H_B,
and one part associated with the evolution of the competition ξ_X.
Therefore, the reparametrization makes explicit how the dominance transition enters the dynamics of growth.
However, near equality, ξ_X is a degenerate coordinate: it does not distinguish which sector dominates. To cross equality regularly, it is better to use the signed variable s_X.
On the B-dominant branch,
C(ξ_X) = 2/(1 + s_X).
Then
ln C = ln 2 − ln(1 + s_X),
and
d ln C/dN = −s_X′/(1 + s_X).
This form is regular once the dominance chart is specified. The singularity of C′(ξ) at ξ = 1/4 reflects the branching of the coordinate ξ, not a physical divergence in the expansion H.

15. Geometric Interpretation: Dyck Paths
The Catalan numbers also count paths called Dyck paths.
Intuitively, a Dyck path is a trajectory that goes up and down but never crosses a forbidden boundary.
This image is useful for visualizing cosmology in Catalan coordinates.
On the regular branch, we choose B as the dominant sector. The condition
H_X² ≤ H_B²
acts as a barrier: as long as this condition holds, we remain on the same branch.
The equality
H_X² = H_B²
is the contact with the boundary.
If X becomes dominant, we change charts: the former complementary sector becomes the new dominant background.
Dyck paths provide a geometric image of the Catalan form: allowed trajectories remain on the chosen dominance branch until they touch the boundary of equal dominance.
This interpretation is pedagogical. It does not claim that the Universe literally counts discrete paths. It shows that the algebraic structure of the transition is the same as that of a classical family of constrained paths.

16. Summary of Advantages
The Catalan reparametrization offers several conceptual and technical advantages.
Exactness.
The identity
y = 1 + ξ_X y²
is exact.
Natural transition variable.
The quantity
ξ_X = Ω_X(1 − Ω_X)
measures competition, not merely abundance.
Universal boundary.
Every positive binary transition has the bound
0 ≤ ξ_X ≤ 1/4.
Equality as a branch point.
The condition
ξ_X = 1/4
is equivalent to
H_X² = H_B².
Rigid series.
The expansion in ξ_X has Catalan coefficients:
1, 1, 2, 5, 14, 42, …
Separation between proximity and direction.
The pair
(ξ_X, s_X)
separates the intensity of competition from the orientation of dominance.
Comparison of models.
Different dark-energy models can be compared through the trajectory
ξ_de(z).
Branch clarity.
The regular Catalan branch describes the complementary sector while it is subdominant. When it becomes dominant, one must change the background chart or explicitly use the signed variable s_X.

17. Limits of the Construction
The construction presented in this article must be understood precisely.
It does not introduce a new energy component.
It does not alter the Friedmann equation.
It does not replace the dynamical analysis of acceleration, perturbations, or structure growth.
It does not claim that the Catalan numbers microscopically govern the Universe.
What it does is simpler and cleaner:
it shows that every positive binary decomposition of the flat Friedmann equation admits a universal transition coordinate,
ξ_X = Ω_X(1 − Ω_X),
in which the ratio
H²/H_B²
satisfies exactly the functional equation of the Catalan generating function.
Therefore, the main contribution is an exact and pedagogically useful algebraic reparametrization.

18. Conclusion
The flat Friedmann equation, when written as a sum of two positive sectors,
H² = H_B² + H_X²,
admits an exact reparametrization by the generating function of the Catalan numbers.
The core of the construction is
H²/H_B² = C[Ω_X(1 − Ω_X)],
on the regular branch in which X is subdominant.
The variable
ξ_X = Ω_X(1 − Ω_X)
is the natural coordinate of the dominance transition. It is small when there is clear dominance, maximal when the sectors are equal, and bounded by
0 ≤ ξ_X ≤ 1/4.
The point
ξ_X = 1/4
is simultaneously the branch point of the Catalan function and the condition of cosmological equality
H_X² = H_B².
Thus, the reparametrization turns the Friedmann sum into a language of transition.
FLRW history may be organized as a sequence of Catalan charts:
radiation → matter → Λ.
Each passage is described by the competition between a dominant sector and a complementary sector.
Cosmic acceleration, in turn, is a distinct dynamical event. In the matter–Λ case, acceleration begins at
ξ_Λ = 2/9,
whereas matter–Λ equality occurs at
ξ_Λ = 1/4.
Therefore, the Catalan form does not replace cosmological dynamics; it clarifies its algebraic structure.
The Friedmann equation sums densities; the Catalan reparametrization reveals the algebraic grammar of its transitions.

References
[1] S. Weinberg, Cosmology, Oxford University Press, 2008.
[2] S. Dodelson, Modern Cosmology, Academic Press, 2003.
[3] B. Ryden, Introduction to Cosmology, Cambridge University Press, 2017.
[4] Planck Collaboration, “Planck 2018 results. VI. Cosmological parameters,” Astronomy & Astrophysics 641, A6, 2020.
[5] A. G. Riess et al., “Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant,” Astronomical Journal 116, 1009, 1998.
[6] S. Perlmutter et al., “Measurements of Omega and Lambda from 42 High-Redshift Supernovae,” Astrophysical Journal 517, 565–586, 1999.
[7] R. P. Stanley, Enumerative Combinatorics, Vol. 2, Cambridge University Press, 1999.
[8] P. Flajolet and R. Sedgewick, Analytic Combinatorics, Cambridge University Press, 2009.

[supprimé]

Writing into the void of a forum sometimes feels like shouting inside an anechoic chamber: either you hear only your own pulse, or the white noise of people who never entered the room in the first place.

So this is the fourth and final act.

This is not a “theory of everything,” but a rearrangement of standard pieces. If the pieces are canonical but the mosaic annoys you, the problem may be your expectation of order, not the ingredients.

Just Landauer, horizons, and high-school algebra pushed further than people seem willing to tolerate.

The Landauer–Horizon Bridgestone

Landauer’s principle says that irreversibly erasing one bit of information in contact with a heat bath at temperature T costs at least

Eₘᵢₙ ≥ k_B T ln 2.

Near a causal horizon, the temperature is not decorative. It is fixed by geometry.

For a horizon with surface gravity κ,
T_H = ℏκ / 2πc k_B.

Insert this into Landauer:

E_bit ≥ k_B · ℏκ / 2πc k_B · ln 2.

The k_B cancels:

E_bit ≥ ℏκ ln 2 / 2πc.

So the cost of a bit is set not by a laboratory machine, but by the causal geometry of the horizon.

The issue is operational:
causal inaccessibility ≟ operational erasure.

More precisely, let

𝒩_R : ρ ↦ ρ_R

be the CPTP restriction map associated with the observer’s causal domain R. For a full-rank local reference state σ, define

𝒟_DPI(R)
= 1 / ln 2 [D(ρ ∥ σ) − D(𝒩_Rρ ∥ 𝒩_Rσ)].

By data processing,

D(𝒩_Rρ ∥ 𝒩_Rσ) ≤ D(ρ ∥ σ),

so

𝒟_DPI(R) ≥ 0.

Equivalently, one may use the recoverability deficit
𝒟_rec(R)

= 1 / ln 2 inf_{ℛ_R} D(ρ ∥ (ℛ_R ∘ 𝒩_R)(ρ)).

Below, 𝒟 denotes the hydrodynamic, finite, non-area-extensive part of this operational opacity.

If information crosses a horizon and becomes inaccessible to a given observer, it is not globally destroyed. The precise claim is operational: the observer is restricted to a causal algebra, and the inaccessible degrees of freedom are traced out or coarse-grained. If the lost correlations are not recoverable by an admissible recovery map, the loss is an operational erasure for that observer.

Rejecting this is not a minor objection; it strikes directly at the logic behind black-hole thermodynamics and horizon entropy.

Jacobson with the entropy bookkeeping exposed

Now use natural units:

ℏ = c = k_B = 1.

For a local Rindler horizon,

T_U = κ / 2π.

If δI_erased bits become operationally inaccessible across that horizon, Landauer gives

δQ = T_U ln 2 · δI_erased.

Define the operational entropy variation

δS_op ≡ ln 2 · δI_erased.

Then

δQ = T_U δS_op.

That is Clausius, but with the entropy bookkeeping made explicit.

In Einstein gravity,

S_H = A / 4G,

so

δS_H = δA / 4G.

If we identify

δS_op = δS_H,

then

δQ = T_U δA / 4G.

This is precisely the local Clausius relation used by Jacobson.

Jacobson’s argument then does the rest: impose

δQ = TδS

for all local Rindler horizons, use Raychaudhuri, and obtain

G_μν + Λg_μν = 8πG T_μν.

The Landauer–Horizon bridge does not replace Jacobson. It sharpens the informational interpretation of the Clausius step.

The nontrivial assumption is exactly this:

δS_op = δA / 4G.

This is not derived from Landauer alone. It is the geometric closure assumption: operational entropy is identified with horizon entropy in the local-equilibrium Jacobson limit.

Why H⁴?

Move from a local Rindler horizon to a cosmological horizon.

In de Sitter, and approximately in quasi-de Sitter or adiabatic FLRW regimes,

T_H ≈ H / 2π.

Therefore

T_H ∼ H.

Landauer gives the cost per bit:

E_bit ≈ H ln 2 / 2π ∼ H.

In spatially flat adiabatic FLRW, the apparent-horizon/Hubble radius scales as

R_H ∼ H⁻¹,

so the corresponding causal volume scales as

V_H ∼ H⁻³.

Hence the causal-volume density scales as

V_H⁻¹ ∼ H³.

Now assume that the effective operational information sector per Hubble volume is finite and non-area-extensive:

N_eff = η𝒟.

Then

ρ_I ∼ N_eff · E_bit · V_H⁻¹.

Substitute the scalings:

ρ_I ∼ constant · H · H³.

Therefore

ρ_I ∼ H⁴.

One H comes from horizon temperature.
Three H’s come from causal volume density.
That is the whole mechanism.

More explicitly, using

V_H = 4π/3 H⁻³

and

E_bit = H ln 2 / 2π,

ρ_I ∼ N_eff (H ln 2 / 2π)(3H³ / 4π),

so

ρ_I ∼ 3 ln 2 / 8π² N_eff H⁴.

The coefficient is model-dependent. The H⁴ power follows once the participating sector is finite per causal volume rather than area-extensive.

There is one important caveat. If N_eff scaled like the full Bekenstein–Hawking horizon entropy,

N_eff ∼ M_Pl² / H²,

then the result would become

ρ ∼ M_Pl²H²,

which is the usual area-law/holographic scaling.

So the distinction is clean:

H² = area-extensive holographic sector.
H⁴ = finite local operational sector per causal volume.

A convenient parametrization is

ρ_I = 3M_Pl² αη𝒟H⁴,

with

α = ℓ_P² ln 2 / π.

Here M_Pl is the reduced Planck mass,

M_Pl⁻² = 8πG,

ℓ_P² = G.

Hence

3M_Pl²α = 3 ln 2 / 8π²,

matching the explicit causal-volume coefficient when

N_eff = η𝒟.

The dimensions are correct:

[ρ_I] = [H⁴].

The algebraic self-regulation at ξ = 1/4

This is the part people should not be able to casually shrug off.

Insert the H⁴ term into the Friedmann equation:

H² = H_bg² + αη𝒟H⁴.

Define

A ≡ αη𝒟 ≥ 0.

Then

H² = H_bg² + AH⁴.

Let

X ≡ H².

Then

X = H_bg² + AX².

Rearrange:

AX² − X + H_bg² = 0.

The H⁴ correction has turned the Friedmann equation into a quadratic equation in

X = H².

Now define the dimensionless variables

y ≡ H² / H_bg²,

ξ ≡ AH_bg² = αη𝒟H_bg².

Then

y = 1 + ξy²,

or

ξy² − y + 1 = 0.

Solve:

y_±(ξ) = [1 ± √(1 − 4ξ)] / 2ξ.

For y to be real,

1 − 4ξ ≥ 0.

Therefore

ξ ≤ 1/4.

That is the origin of the bound. The number 1/4 is not fitted.
It is the discriminant.

Now choose the physical branch. The branch must recover ordinary GR when the correction disappears:

ξ → 0 ⇒ y → 1.

That branch is

y_−(ξ) = [1 − √(1 − 4ξ)] / 2ξ.

Equivalently, in a form regular at ξ = 0,

y_−(ξ) = 2 / [1 + √(1 − 4ξ)].

Expanding around ξ = 0,

y_−(ξ) = 1 + ξ + 2ξ² + 5ξ³ + ⋯,

so

y_− → 1

as

ξ → 0.

The other branch behaves as

y_+ ∼ 1/ξ

and diverges in the GR limit.

So the GR-continuous branch is unique.
At the critical point,

ξ = 1/4,

we get

y_− = 2.

Since

y = H² / H_bg²,

the physical branch obeys

1 ≤ H² / H_bg² ≤ 2.

Therefore

H² ≤ 2H_bg²,

or

H ≤ √2 H_bg.

So the H⁴ correction can increase the expansion, but not without limit, as long as we remain on the real GR-continuous branch.

The discriminant regulates the correction.

Strictly speaking, this is an algebraic branch-regulation, not yet a full dynamical self-regulation theorem. A dynamical proof further requires the evolution of 𝒟 to preserve

ξ(N) ≤ 1/4.

The informational fraction is bounded too.

Since

ρ_I = 3M_Pl²AH⁴,

the physical fractional density is

Ω_I = ρ_I / 3M_Pl²H² = AH².

Using

ξ = AH_bg²

and

y = H² / H_bg²,

we get

Ω_I = ξy.

From the quadratic,

y = 1 + ξy²,

so

ξy² = y − 1.

Divide by y:

ξy = 1 − 1/y.

Therefore

Ω_I = 1 − 1/y.

Since the physical branch satisfies

1 ≤ y ≤ 2,

we obtain

0 ≤ Ω_I ≤ 1/2.

So the informational sector cannot exceed half of the physical critical density on the GR-continuous branch.

The whole chain

Landauer:
E_bit = T ln 2.

Horizon temperature:
T_H ∼ H.

Therefore:
E_bit ∼ H.

Causal volume:
V_H⁻¹ ∼ H³.

Finite non-area-extensive operational sector:
N_eff = η𝒟,

finite and non-area-extensive per causal volume.
Therefore:
ρ_I ∼ H · H³ ∼ H⁴.

Insert into Friedmann:
H² = H_bg² + AH⁴.

Let
y = H² / H_bg²,
ξ = AH_bg².

Then

ξy² − y + 1 = 0.

Discriminant:
1 − 4ξ ≥ 0.

Therefore:
ξ ≤ 1/4.

On the GR-continuous branch:
H² ≤ 2H_bg²,
Ω_I ≤ 1/2.

If horizon-limited information loss is treated operationally through a CPTP restriction and recoverability deficit, if its irreversible component carries a Landauer cost at the Unruh/Hawking/Gibbons–Hawking or apparent-horizon temperature, and if the participating degrees of freedom form a finite non-area-extensive sector per causal volume,
N_eff = η𝒟,
then the associated effective density scales as
ρ_I ∝ 𝒟H⁴.
Once this term enters Friedmann as
H² = H_bg² + αη𝒟H⁴,
the GR-continuous branch is algebraically regulated by the discriminant
ξ = αη𝒟H_bg² ≤ 1/4.
On that branch,
H² ≤ 2H_bg²,
Ω_I ≤ 1/2.

So the chain is:
information
→ Landauer cost
→ horizon temperature
→ operational entropy
→ Clausius
→ Einstein/Jacobson
→ H⁴
→ ξ = 1/4.
Where exactly does the chain fail?

Minimal assumption
Consider an effective dimensionless number of participating operational bits per Hubble volume:
N_eff = η𝒟
The crucial assumption is:
N_eff does not scale like A/G.
In other words, we are not counting the full Gibbons–Hawking area entropy of the horizon. We are counting only the effective operational degrees of freedom that actually participate.

1. Horizon temperature
For de Sitter, or quasi-de Sitter adiabatic expansion,
T_H = H / 2π
Therefore,
T_H ∼ H
The horizon temperature is set by the expansion rate.

2. Landauer cost per bit
Landauer’s principle gives the minimum energy cost per erased bit:
E_bit = T_H ln 2
Substituting the horizon temperature:
E_bit = (H ln 2) / 2π
Therefore,
E_bit ∼ H
Each erased operational bit carries an energy scale proportional to H.

3. Causal Hubble volume
The Hubble radius is
R_H = H⁻¹
So the Hubble volume is
V_H = (4π/3) H⁻³
Hence,
V_H⁻¹ = (3/4π) H³
The causal volume density scales as
V_H⁻¹ ∼ H³

4. Energy density
The total information energy in one Hubble volume is
E_I = N_eff E_bit
Therefore, the associated energy density is
ρ_I = E_I / V_H
or
ρ_I = N_eff E_bit V_H⁻¹
Substitute the expressions above:
ρ_I = N_eff · (H ln 2 / 2π) · (3H³ / 4π)
Therefore,
ρ_I = (3 ln 2 / 8π²) N_eff H⁴
Since
N_eff = η𝒟
we obtain
ρ_I = (3 ln 2 / 8π²) η𝒟 H⁴
Thus,
ρ_I ∝ H⁴

Equivalent Planck-normalized form
Using the reduced Planck mass convention
M_Pl⁻² = 8πG
and
ℓ**_P² = G**
define
α = ℓ**_P² ln 2 / π**
Then
3M_Pl²α = 3(1/8πG)(G ln 2 / π)
so
3M_Pl²α = 3 ln 2 / 8π²
Therefore the same density can be written as
ρ_I = 3M_Pl² α η𝒟 H⁴

One-line derivation
T_H ∼ H
→ Landauer:
E_bit ∼ T_H ln 2 ∼ H
→ causal density:
V_H⁻¹ ∼ H³
→ energy density:
ρ_I ∼ E_bit V_H⁻¹ ∼ H · H³ ∼ H⁴
Therefore,
ρ_I ∼ N_eff H⁴
and for dimensionless, non-area-scaling N_eff,
ρ_I ∝ H⁴

Final statement
The H⁴ scaling follows cleanly if the horizon-information sector contains a finite effective number of operational degrees of freedom per Hubble volume.
The logic is:
horizon temperature:
T_H ∼ H
Landauer cost per bit:
E_bit ∼ H
causal volume density:
V_H⁻¹ ∼ H³
therefore
ρ_I ∼ H × H³ = H⁴
The coefficient η𝒟 encodes how many operational information units participate and how efficiently they gravitate.

We work in natural units,
ℏ = c = k_B = 1.

The derivation is built from four canonical ingredients:

□ Landauer + Unruh + Bekenstein–Hawking + Jacobson

Each ingredient plays a distinct role:

□ Landauer gives the cost; Unruh gives the temperature; Bekenstein–Hawking gives the entropy; Jacobson gives Einstein.

1. Local horizon and Unruh temperature
Consider an arbitrary point p in spacetime and a locally accelerated observer. By the equivalence principle, in a sufficiently small neighborhood of p, this observer possesses a local causal Rindler horizon.
The temperature associated with this horizon is the Unruh temperature,
□ T_U = κ / 2π
where κ is the proper acceleration, or local surface gravity.
Thus, the local horizon functions as the natural thermal reservoir associated with degrees of freedom that are causally inaccessible to the observer.

2. Erasure cost: Landauer
If degrees of freedom cross the horizon, the corresponding microscopic information becomes inaccessible to the local observer. Operationally, this loss of access can be represented as irreversible information erasure.
By Landauer’s principle, erasing δI_erased bits of information costs, at minimum,
δQ ≥ T_U ln 2 · δI_erased.
In the reversible limit, the inequality is saturated:
□ δQ = T_U ln 2 · δI_erased.
We define the erased operational entropy as
□ δS_op ≡ ln 2 · δI_erased.
Therefore,
□ δQ = T_U δS_op.
This is the Landauer–Unruh bridge:
□ causal informational erasure ⇒ thermal flux through the local horizon.

3. Identification with horizon entropy
To connect this bridge to gravity, the relevant operational entropy is identified with the variation of horizon entropy:
□ δS_op = δS_H.
For Einstein gravity, the Bekenstein–Hawking entropy is
S_H = A / 4G.
Therefore,
□ δS_H = δA / 4G.
Thus,
□ δS_op = δA / 4G
and the Landauer–Unruh relation becomes
□ δQ = T_U δA / 4G.
This is precisely the local Clausius relation,
□ δQ = T δS,
with
T = T_U, δS = δS_H.

4. Energy flux through the horizon
Let k^μ be the null tangent vector to the generators of the local horizon, and let λ be an affine parameter chosen so that λ = 0 at the point p.
The approximate Killing vector generating local boosts is
χ^μ = −κλ k^μ.
The horizon surface element is
dΣ^ν = k^ν dλ dA.
The energy flux through the horizon is
δQ = ∫_H T_μν χ^μ dΣ^ν.
Substituting the expressions above,
□ δQ = −κ ∫_H λ T_μν k^μ k^ν dλ dA.
This is the matter side of the derivation.

5. Area variation via Raychaudhuri
The area variation of the horizon is controlled by the expansion θ of the null generators:
δA = ∫_H θ dλ dA.
The Raychaudhuri equation for a null congruence is
dθ/dλ = −(1/2)θ² − σ_μν σ^μν − R_μν k^μ k^ν.
We choose the local horizon to be in instantaneous equilibrium at the point p, namely,
θ(p) = 0, σ_μν(p) = 0.
To linear order in λ,
θ = −λ R_μν k^μ k^ν.
Hence,
□ δA = − ∫_H λ R_μν k^μ k^ν dλ dA.
This is the geometric side of the derivation.

6. Local Clausius implies Einstein
The local Clausius relation is
δQ = T_U δS_H = (κ/2π)(δA/4G).
Therefore,
□ δQ = κ δA / 8πG.
Substituting the expressions for δQ and δA,
−κ ∫_H λ T_μν k^μ k^ν dλ dA
= (κ/8πG) [− ∫_H λ R_μν k^μ k^ν dλ dA].
Canceling the common factors,
T_μν k^μ k^ν = (1/8πG) R_μν k^μ k^ν.
Equivalently,
□ R_μν k^μ k^ν = 8πG T_μν k^μ k^ν
for every null vector k^μ.
Since this holds for every null direction,
R_μν + Φ g_μν = 8πG T_μν
for some scalar function Φ.
Taking the covariant divergence,
∇^μ(R_μν + Φ g_μν) = 8πG ∇^μ T_μν.
Using local energy-momentum conservation,
∇^μ T_μν = 0,
and the Bianchi identity,
∇^μ G_μν = 0,
one obtains
Φ = −(1/2)R + Λ,
where Λ appears as an integration constant.
Therefore,
□ G_μν + Λg_μν = 8πG T_μν.
This is Einstein’s equation.

One-line version

□ δI_erased →[Landauer] δS_op = ln 2 · δI_erased →[Unruh] δQ = T_U δS_op →[δS_op = δA/4G] δQ = T_U δA/4G →[Jacobson] G_μν + Λg_μν = 8πG T_μν.

Conceptual synthesis
The Landauer–Unruh bridge converts operational information loss into thermal horizon flux:
□ δQ = T_U ln 2 · δI_erased.
The Bekenstein–Hawking entropy identifies this information loss with area variation:
□ ln 2 · δI_erased = δA / 4G.
Thus,
□ δQ = T_U δA / 4G.
This is the local Clausius relation. Requiring it to hold for all local Rindler horizons, together with the Raychaudhuri equation, forces the geometry to satisfy
□ G_μν + Λg_μν = 8πG T_μν.

Step 1: The Information Pillar — Landauer’s Principle

Landauer’s Principle states that the irreversible erasure — or the loss of operational access — of 1 logical bit of information in a system in equilibrium with a thermal bath of temperature T requires a minimum energy dissipation given by:

Eₘᵢₙ = kᴮT ln 2

Where:

• kᴮ is the Boltzmann constant — the bridge between the microscopic world and macroscopic thermodynamics.

• T is the temperature of the thermal bath.

• ln 2 arises from the change in Shannon entropy for the fundamental binary choice: 0 or 1.

Step 2: The Quantum-Relativistic Pillar — Horizon Temperature

In Quantum Field Theory in curved spacetime, any observer limited by a causal horizon — whether the event horizon of a black hole, the cosmological horizon of an expanding universe, or the Rindler horizon for constant acceleration — perceives a thermal bath.

The unified temperature for a causal horizon is determined by its surface gravity κ. The general Unruh-Hawking form is:

Tₕ = ℏκ / 2πckᴮ

Where:

• ℏ is the reduced Planck constant — the quantum of action.

• c is the speed of light — the causal limit of spacetime.

• κ is the surface gravity — the geometric intensity associated with the horizon.

Step 3: The Fusion — Causal Erasure

If we assume that the horizon acts as the thermal reservoir that “absorbs” the information that has become causally inaccessible to the observer, the temperature T in Landauer’s Principle is replaced by the horizon temperature Tₕ.

We substitute the quantum-relativistic expression into the information-theoretic one:

Eₘᵢₙ = kᴮ(ℏκ / 2πckᴮ) ln 2

Step 4: The Thermodynamic Cancellation and the Fundamental Result

The Boltzmann constant kᴮ appears both outside the temperature expression and inside its denominator. It cancels exactly:

Eₘᵢₙ = ℏκ ln 2 / 2πc

or, equivalently:

Eₘᵢₙ = (ℏκ / 2πc) ln 2

Dissection of the “DNA Equation”

The beauty of this final equation lies in the elimination of kᴮ. This does not mean thermodynamics has disappeared. Rather, it means that once the temperature is supplied by a horizon, the thermal scale is already encoded in quantum and geometric quantities.

Look at the structure of the final equation:

1. ℏ — Quantum Mechanics It sets the quantum scale of the process. The cost per bit is not purely classical; it carries the quantum grain of action.
2. κ and c — General Relativity κ encodes the surface gravity of the horizon, while c encodes the causal speed limit of spacetime. Together, they show that the cost is dictated by causal geometry.

Therefore, the most precise statement is:

when a horizon renders 1 bit operationally inaccessible to a given observer, the minimum Landauer cost associated with that loss of access is

Eₘᵢₙ = (ℏκ / 2πc) ln 2.

Thus, the equation acts as a compact bridge between information, quantum mechanics, and spacetime geometry: the price of a bit is not set by a material machine, but by the surface gravity of the horizon that limits what the observer can access.