High-Precision Approximation of Riemann Zeros via the Truncated Weil Form

The Connes–van Suijlekom truncated Weil quadratic form, indexed by a cutoff parameter c that controls the primes p ≤ c entering the operator, produces a ground state whose Fourier–Mellin zeros provably lie on the critical line; whether they converge to the Riemann zeros as c → ∞ is open (Connes 2026; Connes–Consani–Moscovici 2025). We present, to our knowledge, the first independent public implementation of the Connes–van Suijlekom Galerkin matrix at sixteen cutoffs (c = 13 through 67, plus c = 100). Across the in-sample window c = 13 through c = 67 at N = 100, the first-zero absolute error |γ1 − γ1Riemann| shrinks monotonically from ∼2×10−55 to ∼1.5×10−168, a 113-OOM convergence across fifteen cutoffs. The smallest-positive even-sector eigenvalue λmineven separately reaches ∼10−334 at c = 100, N = 250 (275-OOM span from c = 13). Out-of-sample test at c = 100. On the four-point N-sweep N ∈ {100, 150, 200, 250} at dps = 500, consecutive first-difference ratios 0.837 and 0.836 match to two decimal places. Aitken-Δ2 on the two overlapping triples yields log10|λ∞even| ≈ −536.8 and ≈ −533.7, approaching the Connes 2026 §6.4 heuristic prediction (≈ −530.4) monotonically with N (6.4 and 3.3 OOM gaps out of |x∞| ∼ 530). The same eigenvector recovers γ1, …, γ10 to 307–329 matching digits at N = 250, dps = 500. Under the unitary equivalence with Connes–Consani–Moscovici Lemma 5.1, this is the deepest such Galerkin-truncation recovery in the public Connes–van Suijlekom / Connes–Consani–Moscovici literature, subject to a hypothesis-status caveat. The raw finite-N matrix carries a small block of negative-sign eigenvalues at the finite archimedean cutoff T = 800; these are an artifact of that cutoff and are absent once T is increased, so the smallest-positive even-sector eigenvalue is the genuine smallest one (continuum positivity of QWλ is RH-equivalent and is not assumed at λ = √100). The fit |log10 λmin| ≈ 13.24 c0.634 on c ≤ 67 at N = 100 is shown to be a finite-N rate, falsified at c = 100, N = 200 by 49 OOM in the direction of faster decay. Structural observations include approximate eigenvector c-invariance (overlap ≥ 0.9498 on all 105 cutoff pairs despite eigenvalues differing by 113 OOM), multi-zero convergence universality (all ten detectable zeros within 3.8% of each other), an empirical Galerkin-convergence exponent s(c) ≈ 55 log c − 128, un-rescaled Galerkin bulk-spectrum Poisson statistics (β < 0.05; this is a structural diagnostic of the truncated operator, not a test of Montgomery's conjecture, which applies to locally-rescaled zero spacings), and tight bulk invariants log|det Qc| ≈ −65.6 c + 542 (R2 = 0.997). We make no claim of proof; the contribution is reproducible numerical data and its careful interpretation under the existing CvS / CCM framework. All code, data, and ancillary files are publicly available. Version 3.3 (2026-06-26) correction. The negative-sign eigenvalue blocks reported at c = 100 and for L(s, χ3) at c = 23, 29 are a finite archimedean-cutoff (T) artifact, not a feature of the operator: they are stable in working precision but vanish once T is increased, so cutoff-free the relevant even sectors are non-negative and the smallest-positive branch is the genuine smallest eigenvalue. No quantitative result changes. See ERRATA.md and the paper's note added in revision. The cutoff sensitivity was independently identified by B. W. A. Silva, consistent with the naturally even, positive ground state reported by R. Andrews; the investigation was prompted by A. Connes.