How the Riemann Hypothesis emerges as the unitarity condition for holographic RG flow
1. The Two Physical Anchors
The argument rests on two facts that are not assumptions but theorems of physics.
The first is the Planck luminosity bound. In classical General Relativity, the maximum rate at which energy can be transported across any surface is:
$$P_{\max} = \frac{c^5}{2G}$$
In a de Sitter universe of radius $R$, the Hawking-Gibbons temperature is $T_{\text{dS}} = \hbar c / 2\pi R$. The entropy production rate cannot exceed $P_{\max}/T_{\text{dS}}$, which means the information scrambling rate — measured by the Lyapunov exponent $\lambda_L$ governing the growth of out-of-time-order correlators — is strictly bounded. This is the Maldacena-Shenker-Stanford theorem (2016):
$$\lambda_L \leq \frac{2\pi T_{\text{dS}}}{\hbar} = \frac{c}{R}$$
The vacuum, in this framework, is a maximally chaotic system operating at precisely this ceiling.
The second anchor defines the resolution of the vacuum. For any mass $M$, two characteristic scales exist: the Schwarzschild radius $r_s = 2GM/c^2$ and the Compton wavelength $\lambda_C = \hbar/Mc$. At the Planck mass $M_p = \sqrt{\hbar c / G}$, these intersect:
$$r_s(M_p) = 2L_p, \qquad \lambda_C(M_p) = L_p$$
This intersection is not a coincidence of units. It identifies the unique scale at which gravitational collapse and quantum uncertainty become simultaneously operative. Any attempt to localise energy exceeding $M_p c^2$ into a region smaller than $2L_p$ produces a black hole — the mode ceases to be a particle and becomes a geometric object. The UV cutoff $y_{\min} = 2L_p$ is not a regulator inserted by hand. It is the minimum support radius for any quantum field mode, forced by the simultaneous validity of the uncertainty principle and the formation condition for a black hole.
These two anchors together define the physical arena: a quantum system saturating the MSS bound, defined on a radial domain $y \in [2L_p, R]$, requiring an effective theory that takes both boundaries seriously.
2. JT Gravity as the Exact Effective Theory
A near-extremal horizon saturating the MSS bound has a uniquely determined effective description. The near-AdS$_2$ region of any such horizon — whether a de Sitter static patch, a near-extremal Reissner-Nordström black hole, or an SYK fixed point — flows in the infrared to Jackiw-Teitelboim gravity with a Schwarzian boundary action. This is a theorem, not a model choice (Almheiri-Polchinski 2015; Maldacena-Stanford-Yang 2016). The JT action is the unique relevant deformation of AdS$_2$ consistent with near-extremality, and it is exactly the theory produced by the spontaneous breaking of reparametrisation symmetry:
$$\text{Diff}(S^1) \to \text{SL}(2,\mathbb{R})$$
The residual unbroken symmetry is the one-dimensional conformal group, generated by three operators: time translations $H$, the dilation operator $D$, and special conformal transformations $K$. The dilation operator takes the canonical form:
$$D = \frac{1}{2}(xp + px)$$
In the holographic dictionary, $D$ generates the renormalization group flow from the UV boundary at $y = 2L_p$ to the IR horizon at $y = R$. The bulk Euclidean AdS$_2$ geometry is precisely the hyperbolic plane $\mathbb{H}^2 \cong \text{SL}(2,\mathbb{R})/\text{SO}(2)$, so the bulk geometry already encodes the representation theory governing the boundary spectrum.
The spectral problem is now precisely defined: find the spectrum of $D$ on $L^2([2L_p, R], , dy/y)$ with a reflecting boundary condition at the UV wall $y = 2L_p$. Because the domain is bounded, the spectrum is discrete. The operator $D$ on the half-line has von Neumann deficiency indices $(1,1)$, admitting a one-parameter family of self-adjoint extensions parametrized by a boundary phase $\theta \in [0, 2\pi)$. The physical content of the remainder of this argument is to identify which $\theta$ the vacuum selects, and why.
3. The Non-Perturbative Path Integral and Arithmetic Modularity
Pure JT gravity with a fixed topology has a continuous $\text{SL}(2,\mathbb{R})$ symmetry and cannot generate the discrete integer structure required for spectral rigidity. The transition from continuous to arithmetic symmetry occurs when the gravitational path integral is taken seriously as a sum over all topologies.
Saad, Shenker, and Stanford (2019) proved that the non-perturbative completion of JT gravity — summing over all two-dimensional geometries consistent with given boundary data — produces a partition function whose spectral statistics are those of the Gaussian Unitary Ensemble. The spectral density is the Plancherel measure of $\text{SL}(2,\mathbb{R})$:
$$\rho(E) = \frac{1}{4\pi^2}\sinh(2\pi\sqrt{E})$$
and the spectral form factor develops the linear ramp $K(t) \sim t/2\pi$ at late times, the universal GUE signature.
Crucially, the partition function of the full non-perturbative theory exhibits modular covariance under $\beta \to 4\pi^2/\beta$, reflecting an underlying symmetry of the space of boundary conditions. The work of Mertens and Turiaci (2022) sharpens this: in the Virasoro minimal string completion of JT gravity, the partition function transforms as a modular object under a discrete subgroup of $\text{SL}(2,\mathbb{R})$.
The precise claim required to complete the argument is that the full modular symmetry group acting on the boundary moduli space is the arithmetic group $\text{SL}(2,\mathbb{Z})$, making the effective moduli space the modular surface:
$$\mathcal{M} = \text{SL}(2,\mathbb{Z}) \backslash \mathbb{H}^2$$
This is an open conjecture, supported by the SSS and Mertens-Turiaci results but not yet established in full generality. It is the single unproven step in what follows. It should be understood as such. Every subsequent step in the chain is either a proved theorem or an exact equivalence.
Note carefully what this quotient means: $\text{SL}(2,\mathbb{Z})$ acts on the boundary moduli space — the space of boundary temperatures and couplings — not on the bulk AdS$_2$ spacetime. The bulk geometry remains the hyperbolic plane. It is the space of boundary conditions that acquires arithmetic structure, which is the physically correct statement.
4. The Lax-Phillips Theorem and the Riemann Scattering Matrix
Once the boundary moduli space is identified as $\mathcal{M} = \text{SL}(2,\mathbb{Z})\backslash\mathbb{H}^2$, the spectral problem for $D$ becomes the automorphic wave equation on $\mathcal{M}$, and its scattering theory is governed by a classical result.
Theorem (Faddeev-Pavlov 1967; Lax-Phillips 1976): For the automorphic wave equation on $\text{SL}(2,\mathbb{Z})\backslash\mathbb{H}^2$, the scattering matrix at spectral parameter $s = 1/2 + iE$ is:
$$S\left(\frac{1}{2} + iE\right) = \frac{\xi(1/2 + iE)}{\xi(1/2 - iE)} = e^{2i\theta_{\text{RS}}(E)}$$
where $\xi(s) = \pi^{-s/2}\Gamma(s/2)\zeta(s)$ is the completed Riemann zeta function and:
$$\theta_{\text{RS}}(E) = \operatorname{Im}\ln\Gamma!\left(\frac{1}{4} + \frac{iE}{2}\right) - \frac{E}{2}\ln\pi$$
is the Riemann-Siegel theta function. This is a proved theorem requiring no additional input.
The identification $s = 1/2 + iE$ with $E$ a real physical energy is not an assumption. It arises directly from the representation theory of $\text{SL}(2,\mathbb{R})$: the principal continuous series representations — those corresponding to unitary, non-decaying modes — have scaling dimensions $\Delta = 1/2 + i\nu$ with $\nu \in \mathbb{R}$. The critical line $\text{Re}(s) = 1/2$ is not chosen for convenience; it is the exact location of unitary representations of the symmetry group. Off-critical values of $\text{Re}(s)$ correspond to non-unitary, exponentially growing or decaying representations.
The reflecting boundary condition at $y = 2L_p$ requires the total phase accumulated by a wave travelling from the UV wall to the IR horizon and back to be an integer multiple of $2\pi$. This quantisation condition, combined with the Lax-Phillips scattering phase, gives the Riemann-von Mangoldt formula (proved 1905):
$$N(E) = \frac{E}{2\pi}\ln\frac{E}{2\pi} - \frac{E}{2\pi} + \frac{7}{8} + \mathcal{O}(1/E)$$
The discrete eigenvalues satisfying this condition are exactly the values $E_n = \gamma_n$ for which $\zeta(1/2 + i\gamma_n) = 0$.
5. Unitarity, Self-Adjointness, and the Riemann Hypothesis
We now have a precisely defined operator on a precisely defined domain, with a precisely defined boundary condition. The physical requirement of unitarity — that probability is conserved, that no information leaks through the UV wall — translates into an exact mathematical condition.
The boundary phase selected by the Lax-Phillips scattering matrix is $\theta = \theta_{\text{RS}}(E)$. This is the unique phase simultaneously consistent with: the $\text{SL}(2,\mathbb{Z})$ arithmetic structure of $\mathcal{M}$, the functional equation $\xi(s) = \xi(1-s)$, and the requirement that $S(E)$ is a pure phase (unitary) for all real $E$. The self-adjointness of $D$ under this extension is not assumed — it is the content of the claim to be examined.
$D$ is self-adjoint under the Lax-Phillips extension if and only if all its eigenvalues $E_n$ are real. The eigenvalues are the Riemann zeros $\gamma_n$. Therefore:
$$D \text{ is self-adjoint} \iff \gamma_n \in \mathbb{R} \text{ for all } n \iff \text{Re}!\left(\rho_n\right) = \frac{1}{2} \text{ for all non-trivial zeros } \rho_n$$
This is the Riemann Hypothesis.
To be precise about what would happen if RH is false: a zero $\rho = \sigma + i\gamma$ with $\sigma \neq 1/2$ contributes a pair of modes with scaling dimensions $\Delta = \sigma \pm i\gamma$. Since zeros come in conjugate pairs $\rho, \bar{\rho}, 1-\rho, 1-\bar{\rho}$, an off-critical zero generates both exponentially decaying and exponentially growing modes simultaneously. The growing mode violates the MSS bound, since its Lyapunov exponent $\text{Im}(\Delta) \cdot (c/R) > c/R$ would exceed the maximum scrambling rate. The decaying mode breaks time-reversal symmetry. Together, they destroy the unitarity of the holographic RG flow: probability current leaks through the UV boundary, and the vacuum is no longer a stable fixed point of the renormalization group.
The Riemann Hypothesis is therefore the exact statement that the Planck-scale reflective boundary at $2L_p$, combined with the arithmetic structure of the gravitational path integral, produces a self-adjoint dilation operator — that holographic RG flow is exactly unitary — that the vacuum neither creates nor destroys information at any energy scale.
6. Independent Confirmation: GUE Statistics
The argument above is conditional on the Mertens-Turiaci conjecture. The following provides independent evidence that the physical framework is correct, requiring no conjectures.
Saad, Shenker, and Stanford (2019) proved that the JT gravity spectral form factor satisfies $K(t) \sim t/2\pi$ — the universal ramp of the Gaussian Unitary Ensemble. Montgomery (1973) proved, and Odlyzko confirmed numerically to $10^{22}$ zeros, that the pair correlation function of Riemann zeros is:
$$R_2(r) = 1 - \left(\frac{\sin\pi r}{\pi r}\right)^2$$
This is also the universal GUE pair correlation function.
Two completely different computations — one from the gravitational path integral over two-dimensional topologies, one from the analytic number theory of the Riemann zeta function — produce identical spectral statistics. This is not a numerical coincidence. It reflects the fact that both objects are governed by the same underlying structure: the $\text{SL}(2,\mathbb{R})$ representation theory of the modular surface. The GUE universality class is the fingerprint.
The Weil explicit formula makes this structural identity precise. For any suitable test function $h$:
$$\sum_n h(\gamma_n) = h!\left(\tfrac{i}{2}\right) + h!\left(-\tfrac{i}{2}\right) - \sum_p \sum_{k=1}^\infty \frac{\ln p}{p^{k/2}} \hat{h}(k \ln p) + \cdots$$
This is simultaneously the gravitational trace formula of the holographic boundary theory: the left side is the spectrum of quantum scrambling modes; the right side is the prime geodesic data of the bulk geometry. Prime numbers are the lengths of closed geodesics in the bulk. Riemann zeros are the quasinormal frequencies of the boundary. The Weil formula is the holographic dictionary in spectral form.
7. The Complete Logical Chain
The argument proceeds as follows, with the epistemic status of each step made explicit.
The Planck luminosity bound and MSS theorem together establish that the vacuum is a maximally chaotic system with Lyapunov exponent $\lambda_L = c/R$. This is a theorem.
The Schwarzschild-Compton intersection at $M_p$ establishes that no quantum mode can be localized below $2L_p$, giving a UV cutoff that is derived, not chosen. This is established physics.
The effective theory of a system saturating the MSS bound on a domain $[2L_p, R]$ is uniquely JT gravity with the Schwarzian action, and its primary operator is the dilation operator $D$. This is a theorem.
The non-perturbative sum over topologies in JT gravity produces a partition function with GUE statistics and modular covariance. The SSS result is proved. The full $\text{SL}(2,\mathbb{Z})$ modular invariance of the boundary moduli space — making $\mathcal{M} = \text{SL}(2,\mathbb{Z})\backslash\mathbb{H}^2$ — is an open conjecture supported by Mertens-Turiaci. This is the single unproven step.
Conditional on this conjecture, the Faddeev-Pavlov and Lax-Phillips theorem gives the scattering matrix of $D$ as $\xi(1/2+iE)/\xi(1/2-iE)$. This is a proved theorem.
The quantisation condition between $y = 2L_p$ and $y = R$ produces eigenvalues $E_n = \gamma_n$, the Riemann zeros. This is a proved theorem.
The self-adjointness of $D$ under the Lax-Phillips boundary phase — equivalently, the unitarity of holographic RG flow, equivalently, the stability of the vacuum as a maximally chaotic fixed point — is exactly equivalent to all $\gamma_n$ being real.
The Riemann Hypothesis is the unitarity condition for the renormalisation group flow of the vacuum.
Comments
Post a Comment