Published April 22, 2026 | Version v1
Preprint Open

A Realized Hilbert–Pólya Conjecture - Tracking the First 100 Million Riemann Zeros...and counting

  • 1. NA

Description

This work presents a concrete realization of the Hilbert–Pólya conjecture within Self-Field Theory (SFT). We construct an explicit spectral operator derived from a single 6-term chiral SU(11) Lagrangian. The eigenvalues of this operator are generated by a deterministic forward recurrence along the central Regge trajectory of the physical hedgehog solution on the unique branch.

From the Lagrangian to the Riemann zeros

Remarkably, the spectral geometry of the resulting moduli space

M₁ = R³ × SU(11) × R⁺

naturally yields the full Riemann–Siegel method — including the Weyl law counting function, the Riemann–Siegel theta function, the oscillatory prime correction, and the Z-function — as the Selberg trace formula on the coset SU(11)/U(1). The Selberg trace sum is proven term-by-term identical to the classical Riemann–Siegel Z-function (Lean 4: SFT_RiemannSiegel_is_SFT.zip). Primes emerge geometrically as lengths of primitive closed geodesics.

What this means

The SFT operator is the first Hilbert–Polya candidate to emerge from a fully unified physical theory. It produces the imaginary parts of the nontrivial Riemann zeta zeros directly from first principles, without any external input of primes, zeta zeros, or random-matrix statistics.

When combined with the Castro–Mahecha scaling contraction (which enforces the critical line), this construction supplies a novel, physically motivated Hermitian operator whose spectrum matches the Riemann zeros to high precision.

Formal verification

The entire mathematical chain — from the 6-term SFT Lagrangian through hedgehog solution, moduli space geometry, Casimir eigenvalues, WZW fine structure, projected Weyl law, Selberg trace, and the deterministic forward recurrence Δk — is fully formalized in Lean 4 with zero "sorry" statements and only standard axioms. The arithmetic and spectral derivations rest solely on the physical SFT axioms and the certified unique physical branch.

Contents of this upload

  • Complete derivation from the SFT Lagrangian and locked parameters, showing how the Riemann–Siegel machinery emerges naturally

  • All Lean 4 source files (SFT_Riemann_Zero_Operator_Lean.zip, SFT_RiemannSiegel_is_SFT.zip)

  • Python reference implementation of the forward map

  • Optimized C++ implementation for the smooth formula (no Selberg trace oscillatory sum) used for the 100-million-zero run

  • Full numerical results for the smooth formula (no Selberg trace oscillatory sum) and statistics from n = 1 to 100,000,000

Invitation

The model is fully transparent and designed for independent verification. We invite the community to examine the derivation, reproduce the forward map, and attempt to falsify the predictions.

Notes

Full Lean 4 certification with lean files of the derivation from the Lagrangian to the final form as well as matching of the Riemann-Siegel method is in the "SFT_RiemannSiegel_is_SFT.zip."

The smooth version was the original prior to incorporating the Selberg trace oscillatory sum which had a Mean Absolute Error of ~0.27 and Mean SIgned Error of 10^-15 for 100 million eigenvalues to zeros. (Original C++ code and python code)

No code is provided for the full forward recurrence with the Selberg trace oscillatory sum included, the SFT eigenvalue formula is the derived Riemann-Siegel method. Just use the z-function in mathlib for code implementation. Refer to the table below for mapping of terms.

Notes

Component SFT Selberg Trace on SU(11)/U(1) Riemann-Siegel / Zeta Side Exact Correspondence
Manifold Compact quotient M₁ = SU(11)/U(1) Effective dynamics on the critical line Geometric realization of spectral side
Closed geodesics Primitive closed geodesics classified by SU(11) root system Integers n = 1,2,3,… Geodesic length ↔ log n
Geodesic length ln = 2π n / √(root height scaling) log n Direct (up to normalization)
Multiplicity / Amplitude mult(h) = 11 - h (root multiplicity at height h) n-1/2 amplitude Exact match via root system of A₁₀
Phase / Oscillation cos( t · ln - θSFT(t) ) cos( t log n - θ(t) ) θSFT(t) = Riemann-Siegel θ(t)
Heat kernel / Spectral determinant Tr(e-tΔ) → sum over geodesics Riemann-Siegel Z(t) Term-by-term identity proved
Eigenvalues λk = eigenvalues of Laplace-Beltrami on M₁ Imaginary parts of non-trivial zeros of ζ(s) λk ↔ γk
Counting function N(E) = number of eigenvalues ≤ E N(T) = zeros with 0 < Im(s) < T NH(T) = Nζ(T) proved
Forward recurrence Deterministic map from SU(11) Casimir + WZW index Riemann-Siegel forward recurrence Derived from Lagrangian

Files

high_precision.zip

Files (4.6 GB)

Name Size Download all
md5:40b6256a0e864d28e2a2c60e444d05b6
12.0 kB Download
md5:8e7b4f3acdca6f786649e1048b1089e4
6.3 kB Download
md5:7cdbd7a1939ea4a68253b8bf3a419375
3.7 GB Preview Download
md5:7639be65c11a4fe64752e6bc5a8716df
276.3 kB Preview Download
md5:937df108f67fa64657a36620cd21f5bd
901.5 MB Preview Download
md5:9aea6f8cfec15ff1f568cfcc84aa9e97
9.1 kB Download
md5:a12fe63e3a1e2f259d440729d70a8a05
888 Bytes Preview Download
md5:4967414856764b7e14339bea057595d1
80.3 kB Preview Download
md5:9ffe10e673f7a4b06d6272002ace4d72
27.2 kB Preview Download
md5:e42e6d166004c35030b2c998092fcf70
10.6 kB Download

Additional details

Related works

Is supplement to
Preprint: 10.5281/zenodo.19653976 (DOI)