The Theory of Persistence
Essay · Standard · 16 min

Quasi-Perelman: PT as Asymptotic Expanding Soliton

PT's Fisher–Bianchi metric is a quasi-soliton expanding Ricci soliton whose constant $\lambda(\mu) = -1/\mu^4 + O(1/\mu^5)$ vanishes asymptotically (coefficient $-1$ exact, pair-independent). This nonzero residue co-occurs with the entropic arrow of PT-time — two independent signatures of the same dynamical non-triviality of the sieve.

Go deeper: T5 , T6 , BT13 , BT17 , BT18

The question

Hamilton (1982) introduced the Ricci flow: an equation that deforms a geometry to gradually bring it toward an equilibrium shape. Perelman (2002) showed that this machinery resolves the Poincaré conjecture — a century-old open problem, which earned him the Fields medal (declined). At the heart of the proof: the notion of soliton, a state of absolute equilibrium where the geometry no longer changes.

PT (Persistence Theory) also describes a geometry that evolves with a parameter μ\mu and approaches a stable configuration at μ=15\mu^* = 15. The natural question:

Is PT a Perelman soliton?

The short answer: PT is an asymptotic expanding quasi-soliton of Ricci. The Fisher–Bianchi metric satisfies the soliton equation with constant λ(μ)<0\lambda(\mu) < 0 that vanishes as 1/μ4-1/\mu^4 at large scales, but is never exactly zero at finite μ\mu. A steady soliton in the strict sense would require λ0\lambda \equiv 0 — PT is not one. This nonzero residue co-occurs with the entropic arrow of PT-time — two independent signatures of the same dynamical non-triviality of the sieve.

The picture

Imagine a ball of clay deformed slowly according to a well-defined rule. Perelman proved that under certain conditions, the evolution makes the clay perfectly spherical: it reaches its soliton, its absolute equilibrium form where nothing moves anymore.

PT is not a flow toward a soliton; it is a family of geometries parameterized by μ\mu, with a fixed point at μ=15\mu^* = 15. At that point, the PT metric satisfies the Ricci-soliton equation up to a residue. This residue is measurable, structured, and nonzero. It crushes as 1/μ41/\mu^4 at large scales, but at μ=15\mu^* = 15 it remains small yet significant (105\sim 10^{-5}).

This difference is not an error. It co-occurs with the entropic arrow of PT-time (the qq^- branch). A perfect soliton would correspond to a geometry without this dynamical signature. PT avoids this fate by not quite closing its soliton equation.

A perfect soliton is a frozen ball. PT is a ball that never quite manages to freeze. This residual gap is time.

The slow-roll analogy

In cosmology, we know a similar situation: inflation. The primordial universe was nearly in equilibrium — nearly, but not exactly. This “nearly” — a tiny parameter called slow-roll — is precisely what allowed the universe to expand. Exact equilibrium would have been sterile: nothing would have happened.

PT works the same way, but the slow-roll parameter is forced by arithmetic: it comes directly from the structure of the prime number sieve.

Two branches of the sieve

PT explicitly distinguishes two readings at μ=15\mu^* = 15:

  • q+q^+: the static branch, which measures forces (vertex, couplings, αEM\alpha_{\rm EM}).
  • qq^-: the dynamic branch, which measures flow (propagation, dissipation, arrow of time).

If PT were a steady soliton (λ0\lambda \equiv 0), the anisotropy between the three Hubble rates H3,H5,H7H_3, H_5, H_7 would have to vanish. The entropic arrow on qq^- (BT13) remains independently true, but it co-occurs with the residue: the two say, from two different angles, that PT is dynamically non-trivial.

Why this matters

If the PT–Perelman correspondence were strictly exact, PT would be a particular case of a theory that has already earned a Fields medal. That would be encouraging but not distinctive. The fact that there is a residual gap — with a precise, computable mathematical form — makes PT a theory that dialogues with Perelman while retaining its arithmetic specificity.

Why almost, and why it matters

What an exact soliton would do

If PT were a Ricci soliton in the strict sense, its geometry at μ=15\mu^* = 15 would be in absolute equilibrium. No evolution, no fluctuation, no possible trajectory. Mathematically attractive; physically sterile.

Yet PT explicitly possesses two natural branches at μ=15\mu^* = 15:

  • q+=12/μq^+ = 1 - 2/\mu: the static branch (couplings, αEM\alpha_{\rm EM}, PMNS). The frozen.
  • q=e1/μq^- = e^{-1/\mu}: the dynamic branch (propagator, dissipation, arrow of time). The flow.

Theorem BT13 (EML closure) states strict monotonicity of the Shannon entropy HH on qq^-:

dHdμ>0strictly.\frac{dH}{d\mu} > 0 \quad \text{strictly}.

This property is independently true on the qq^- branch (it does not depend on the Perelman analysis). But it co-occurs with the non-triviality of the Perelman residue: both properties are jointly true in PT, and together they express the dynamical non-triviality of the sieve.

The computation

We write Perelman’s equation for PT:

Ric[gPT]+Hess(f)=λ(μ)gPT\mathrm{Ric}[g_{\rm PT}] + \mathrm{Hess}(f) = \lambda(\mu)\, g_{\rm PT}

with gPTg_{\rm PT} the Fisher-Bianchi metric and f=lnαf = -\ln\alpha the potential. For a perfect soliton, λ0\lambda \equiv 0. For PT, numerically:

λ(15)2.5×105,\lambda(15) \approx -2.5 \times 10^{-5},

a tiny but nonzero value. And crucially, when μ\mu grows large:

λ(μ)=1μ4+O ⁣(1μ5).\lambda(\mu) = -\frac{1}{\mu^4} + O\!\left(\frac{1}{\mu^5}\right).

The coefficient 1-1 is exact. Verified at μ=30,000\mu = 30{,}000 in 60-digit arithmetic: λμ4=1.00011\lambda \cdot \mu^4 = -1.00011.

The conceptual inversion

The initial reading would be: “PT misses the perfect soliton, too bad.” The correct reading is opposite:

The fact that λ0\lambda \neq 0 IS the geometric condition that makes PT time possible. An exact soliton would correspond to a frozen universe, with no temporal arrow, no dynamics.

Proposition (proved in Appendix U and in the transversal chapter on spectral geometry): for the PT Fisher–Bianchi metric, the following two properties are equivalent:

  1. Effective anisotropy of Hubble rates: there exist two active primes p1p2p_1 \neq p_2 with Hp1(μ)Hp2(μ)H_{p_1}(\mu) \neq H_{p_2}(\mu).
  2. Unique determination of λ\lambda: the system of (pp)(pp) equations for p{3,5,7}p \in \{3,5,7\} admits a unique solution λ\lambda.

Proof intuition: with two independent equations (p1p1)(p_1 p_1) and (p2p2)(p_2 p_2) in two unknowns (f˙,λ)(\dot f, \lambda), Cramer’s method gives a unique solution if and only if the determinant Hp2Hp1H_{p_2} - H_{p_1} is nonzero. If all HpH_p are equal to a common HH, the system degenerates into a single equation in two unknowns: λ\lambda becomes a free degree of freedom, undetermined.

Asymptotic corollary: as μ\mu \to \infty, the anisotropy Hp1Hp2=(p2p1)/μ3+O(1/μ4)0H_{p_1} - H_{p_2} = (p_2 - p_1)/\mu^3 + O(1/\mu^4) \neq 0 is automatic (active primes are distinct in Z\Z), so the unique value of λ\lambda satisfies λ(μ)=1/μ4+O(1/μ5)\lambda(\mu) = -1/\mu^4 + O(1/\mu^5), nonzero.

Co-occurrence with BT13. A third property is jointly verified in PT, but it is independent of the previous two:

  • (BT13) dH/dμ>0dH/d\mu > 0 strictly on the qq^- branch: the entropic temporal arrow is carried by qq^- independently of the anisotropy of the γp\gamma_p.

The three properties (1)–(2)–(BT13) are thus jointly true in PT, but their strict equivalence is not proved: (1) ⟺ (2) holds, while (BT13) holds for its own reason. It is this co-occurrence that legitimizes the reading “the residue is the geometric signature of the qq^- branch”, without making it a formal equivalence.

The basin of stability

We work inside the reduced stable basin of the odd sieve: (μ7,μ11)(11.63,17.98)(\mu_7, \mu_{11}) \approx (11.63,\, 17.98) which contains the fixed point μ=15\mu^* = 15. In this interval, the three active primes are exactly {3,5,7}\{3, 5, 7\}, and the reduced map sends all integers {12,13,14,15,16,17}\{12, 13, 14, 15, 16, 17\} to 1515 (theorems BT17 and BT18 of the bridge program). All numerical results are systematically verified on this basin.

Three complementary readings

The deviation from a steady soliton can be read in three complementary ways. Only one (the Perelman reading) is a strict mathematical description; the other two are analogies (qualified below).

Perelman reading (strict mathematical). Perelman studies flows that converge toward a soliton. PT is not a flow: it is a family of geometries parameterized by μ\mu, with a fixed point at μ=15\mu^* = 15. The metric gPT(μ)g_{\rm PT}(\mu^*) satisfies the soliton equation with a constant λ0\lambda \neq 0 that vanishes as 1/μ41/\mu^4 at large scales. PT is therefore an expanding quasi-soliton (λ<0\lambda < 0 in the Perelman convention) whose constant tends to zero asymptotically.

Cosmological analogy (slow-roll). In inflation, an exact de Sitter universe would be sterile: one needs a small parameter ϵsr0\epsilon_{\rm sr} \neq 0 that generates the dynamics. PT’s residue λ1/μ4|\lambda| \sim 1/\mu^4 plays qualitatively the same structural role as a slow-roll parameter: without residue, no dynamics. The analogy is purely structural, not a formal identification (the two parameters do not match term-for-term).

Speculative analogy (Connes, NCG). In noncommutative geometry, time emerges from the modular flow of a KMS state. No modular flow \Leftrightarrow no time. PT’s residue plays a qualitatively comparable role: without residue, no proper dynamics. This analogy remains speculative and is not formalized in the present work.

Closure of the temporal triptych

Three independent results, one mechanism:

ResultTime signature
Chapter 13 (g00<0g_{00} < 0 for μ>μc\mu > \mu_c)Time exists (Lorentzian signature)
Appendix R, BT13Time is oriented (qq^- carries the arrow)
Appendix U, λ0\lambda \neq 0Time is never frozen (arithmetic slow-roll)

Technical proof

Metric, gauge, and method

We work in the Fisher–Bianchi metric in Lorentzian signature (,+,+,+)(-,+,+,+) — the component gμμg_{\mu\mu} becomes negative for μ>μc6.97\mu > \mu_c \approx 6.97, the Lorentz frontier determined by the vanishing of d2(lnα)/dμ2=0d^2(\ln\alpha)/d\mu^2 = 0 (signature theorem 12.4 of Chapter 13):

gPT(μ)=N(μ)2dμ2+p{3,5,7}ap(μ)2dxp2,g_{\rm PT}(\mu) = -N(\mu)^2\, d\mu^2 + \sum_{p \in \{3,5,7\}} a_p(\mu)^2\, dx_p^2,

where ap(μ):=γp(μ)/μa_p(\mu) := \gamma_p(\mu) / \mu and γp\gamma_p is the anomalous dimension at prime pp. We choose the gauge N(μ)=μN(\mu) = \mu (technically favorable), giving proper time τ=μ2/2\tau = \mu^2/2. The basin (μ7,μ11)(11.63,17.98)(\mu_7, \mu_{11}) \approx (11.63,\, 17.98) is strictly above μc\mu_c, hence in the Lorentzian regime.

Pair-compatibility method: the soliton equation Ric+Hess(f)=λg\mathrm{Ric} + \mathrm{Hess}(f) = \lambda g has 4 diagonal components (1 temporal (ττ)(\tau\tau) + 3 spatial (pp)(pp) for p{3,5,7}p \in \{3,5,7\}), for 2 unknowns (λ,f˙)(\lambda, \dot f). The system is overdetermined. The (pp)(pp) components read

ApHpf˙=λ(minus sign in Lorentzian signature),A_p - H_p\,\dot f = \lambda \quad \text{(minus sign in Lorentzian signature)},

with Ap:=H˙p+HpkHkA_p := \dot H_p + H_p \sum_k H_k and Hp:=μ1dlnap/dμH_p := \mu^{-1} d\ln a_p/d\mu. We extract λ\lambda via the 3 pairs (3,5),(5,7),(3,7)(3,5), (5,7), (3,7) by elimination of f˙\dot f. If all three yield the same value, we have a soliton; otherwise, the spread measures the defect.

Asymptotic expansion of γp\gamma_p

For large μ\mu, q+=12/μ1q^+ = 1 - 2/\mu \to 1. The expansion proceeds in four steps:

  1. Binomial: qp=12p/μ+2p(p1)/μ2+O(1/μ3)q^p = 1 - 2p/\mu + 2p(p-1)/\mu^2 + O(1/\mu^3).
  2. Deficit: δp:=(1qp)/p=2/μ2(p1)/μ2+O(1/μ3)\delta_p := (1-q^p)/p = 2/\mu - 2(p-1)/\mu^2 + O(1/\mu^3).
  3. Holonomy T6: sin2θp=δp(2δp)=4/μ4p/μ2+O(1/μ3)\sin^2\theta_p = \delta_p(2-\delta_p) = 4/\mu - 4p/\mu^2 + O(1/\mu^3).
  4. Logarithm and derivation: lnsin2θp=ln4lnμp/μ+(p24)/(6μ2)+O(1/μ3)\ln\sin^2\theta_p = \ln 4 - \ln\mu - p/\mu + (p^2-4)/(6\mu^2) + O(1/\mu^3), then γp=μd(lnsin2θp)/dμ\gamma_p = -\mu\, d(\ln\sin^2\theta_p)/d\mu:
γp(μ)=1pμ+p243μ2+O ⁣(1μ3).\gamma_p(\mu) = 1 - \frac{p}{\mu} + \frac{p^2 - 4}{3\mu^2} + O\!\left(\frac{1}{\mu^3}\right).

Asymptotically, γp1\gamma_p \to 1 for every active prime: the pp-dependence disappears in the limit.

Kinematics in gauge N=μN = \mu

From lnap=lnμ+lnγp\ln a_p = -\ln\mu + \ln\gamma_p we derive (gauge τ\tau, d/dτ=(1/μ)d/dμd/d\tau = (1/\mu)\,d/d\mu):

Hp(μ):=d(lnap)dτ=μ1d(lnap)dμ=1μ2+pμ3+O ⁣(1μ4),H_p(\mu) := \frac{d(\ln a_p)}{d\tau} = \mu^{-1}\, \frac{d(\ln a_p)}{d\mu} = -\frac{1}{\mu^2} + \frac{p}{\mu^3} + O\!\left(\frac{1}{\mu^4}\right),H˙p:=dHpdτ=μ1dHpdμ=2μ43pμ5+O ⁣(1μ6).\dot H_p := \frac{dH_p}{d\tau} = \mu^{-1}\, \frac{dH_p}{d\mu} = \frac{2}{\mu^4} - \frac{3p}{\mu^5} + O\!\left(\frac{1}{\mu^6}\right).

And the useful combination:

Ap(μ):=H˙p+HpkHk=5μ415+6pμ5+O ⁣(1μ6).A_p(\mu) := \dot H_p + H_p\, \sum_k H_k = \frac{5}{\mu^4} - \frac{15 + 6p}{\mu^5} + O\!\left(\frac{1}{\mu^6}\right).

Key step: the coefficient -1

For a pair (p1,p2)(p_1, p_2), the compatibility of the (pp)(pp) equation gives:

λp1p2=Hp2Ap1Hp1Ap2Hp2Hp1.\lambda_{p_1 p_2} = \frac{H_{p_2} A_{p_1} - H_{p_1} A_{p_2}}{H_{p_2} - H_{p_1}}.

We compute the numerator, keeping all terms up to 1/μ71/\mu^7:

Hp2Ap1Hp1Ap2=(1μ2+p2μ3)(5μ415+6p1μ5)(1μ2+p1μ3)(5μ415+6p2μ5)+O(1/μ8).\begin{aligned} H_{p_2} A_{p_1} - H_{p_1} A_{p_2} &= \Big(-\tfrac{1}{\mu^2} + \tfrac{p_2}{\mu^3}\Big)\Big(\tfrac{5}{\mu^4} - \tfrac{15+6p_1}{\mu^5}\Big) \\ &\quad - \Big(-\tfrac{1}{\mu^2} + \tfrac{p_1}{\mu^3}\Big)\Big(\tfrac{5}{\mu^4} - \tfrac{15+6p_2}{\mu^5}\Big) + O(1/\mu^8). \end{aligned}

The 1/μ61/\mu^6 terms cancel (the factor 5/μ6-5/\mu^6 appears in both). Remaining terms in 1/μ71/\mu^7:

15+6p1+5p2μ715+6p2+5p1μ7=p1p2μ7+O ⁣(1μ8).\frac{15+6p_1 + 5p_2}{\mu^7} - \frac{15+6p_2 + 5p_1}{\mu^7} = \frac{p_1 - p_2}{\mu^7} + O\!\left(\frac{1}{\mu^8}\right).

And the denominator:

Hp2Hp1=p2p1μ3+O ⁣(1μ4).H_{p_2} - H_{p_1} = \frac{p_2 - p_1}{\mu^3} + O\!\left(\frac{1}{\mu^4}\right).

The quotient thus gives:

  λp1p2(μ)=1μ4+O ⁣(1μ5).  \boxed{\;\lambda_{p_1 p_2}(\mu) = -\frac{1}{\mu^4} + O\!\left(\frac{1}{\mu^5}\right).\;}

The coefficient 1-1 is independent of the choice of pair among the three active primes {3,5,7}\{3, 5, 7\}. This is what makes the value consistent — without this independence, the pair decomposition would not be well-defined.

The correction term (numerical conjecture)

Numerically (mpmath 60 digits),

λ(μ)=1μ4103μ5+O ⁣(1μ6),\lambda(\mu) = -\frac{1}{\mu^4} - \frac{10}{3\mu^5} + O\!\left(\frac{1}{\mu^6}\right),

with the factor 10/3-10/3 observed at μ=10,000\mu = 10{,}000 and 30,00030{,}000. This is a numerical conjecture not yet analytically proved: the analytic derivation requires extending the γp\gamma_p expansion to order 1/μ31/\mu^3, a heavy but elementary calculation.

At the fixed point μ=15\mu^* = 15

At μ=15\mu = 15 we are far from the asymptotic regime. The approximation 1/μ4-1/\mu^4 gives 1.975×105-1.975 \times 10^{-5}, whereas the exact value is 2.5215×105-2.5215 \times 10^{-5}. The gap is 22%22\% relative to the exact value (or 28%28\% relative to the approximation 1/μ4-1/\mu^4). The 1/μk51/\mu^{k\geq 5} terms exactly account for this gap.

Temptation: a local fit gives λ(12/5)μ127/30\lambda \approx -(12/5)\, \mu^{-127/30} with rel_RMS of 0.07 % on the basin. Rejected as an exact identity: the mpmath calculation shows that the logarithmic derivative dlnλ/dlnμd \ln|\lambda| / d \ln\mu varies from 4.244-4.244 to 4.211-4.211 across the basin (not constant). The fit is a local artifact of the expansion λ=μ4(1+10/(3μ)+)\lambda = -\mu^{-4}(1 + 10/(3\mu) + \ldots), which looks like an effective power law in a limited μ\mu range but is not one.

What remains open

  • Exact closed form of λ(μ)\lambda(\mu) at finite μ\mu: the function is rational (all PT operations are) but its symbolic expressions explode in size. CRT decomposition modulo 3573 \cdot 5 \cdot 7 or factorization in sin2θp\sin^2\theta_p factors.

  • Analytic confirmation of 10/3-10/3 in the correction term.

  • Global bound λ(μ)C/μ4|\lambda(\mu)| \leq C/\mu^4 on the entire basin with explicit CC.

  • Mirzakhani link: the non-power-law form of λ\lambda suggests a multi-parameter recursion, which is exactly the territory of Eynard-Orantin topological recursion.

To close

The gap between PT and a steady Ricci soliton is not a margin of error, nor an artifact of our formulation. It is a geometric invariant, measurable and computable, whose dominant coefficient 1-1 is proved exact (and pair-independent). Its structure — the pp-dependence that vanishes only asymptotically — is the direct trace of the fact that PT operates on a finite number of active directions, {3,5,7}\{3, 5, 7\}.

The residual gap is thus not the thickness of an error margin: it is the thickness of time itself, measured in units of the arithmetic sieve. Where Perelman proves convergence toward a soliton, PT proves the controlled non-realization of that same soliton (in the sense of an expanding soliton whose constant tends to zero), and finds in it the condition for a temporally living geometry.

References


← All essays