17 August 2026 · revised 4 September 2026, v1.5
We reimplemented Anthropic's compressed Weil form from scratch and reproduced their published tables to six significant figures; then ran the one control their own numerics could not run; then built the arithmetic half of their identity, for zeta and for Dirichlet L-functions; then went after the sentence that says how far the method can go — found its proof in their Lean development rather than in the paper — ran that ceiling's own programme at every period up to 128, where no value of it had ever been computed outside their repository; then solved the small-period programme exactly, every configuration enumerated, to certify our own instrument; then built their entire Lean development and passed all of it through the proof-checker's kernel; and then measured the depth-law rate across three conductors and three heights, where it comes out flat in height and split by conductor.
An unreleased Claude research model proved unconditionally that at least two thirds of the non-trivial zeros of zeta in a dyadic window are simple and on the critical line — 0.67250 with the optimal test family, against a record of 5/12 = 0.41666 that had stood since 2020. It also gives at least 5/6 distinct (0.83625 optimised) and the same for a fixed primitive Dirichlet L-function. The work took two sessions in Claude Code, about 650 candidate ideas, roughly 60 subagents and 31 million output tokens; it was validated by two mathematicians and partly formalised in Lean.
The mechanism abandons Levinson's mollifier entirely. Weil's explicit formula gives a real symmetric form; compress it onto a family of windowed exponentials spaced at the mean zero spacing; a zero on the line contributes a positive square, a mirror pair off the line contributes a block of signature (1,1). The hypothesis was classically needed only to read the zero side as a positive sum, and here it is replaced by Sylvester's law of inertia plus a rank–trace inequality. The arithmetic input is Montgomery's 1973 pair correlation at band-width one, which was always unconditional.
anthropic.com/research/riemann-zeta · paper 35 pp · informal note 5 pp · appendix 95 pp · transcripts 116 pp · github.com/anthropics/zeta-23-lean · on Baluyot–Goldston–Suriajaya–Turnage-Butterbaugh, Goldston–Suriajaya, Bombieri, Yoshida.
We implemented the compressed form from the zero side — closed-form window transform, the taper
ϱ(x) = x − sin(2πx)/2π, the lattice, the assembly — and ran it against real zeta zeros.
Nothing was tuned to match. Their published table:
| Quantity, T = 1000 | Their value | Ours |
|---|---|---|
| ‖Â‖²/tr Â, λ=1, η=0.10 | 1.387 | 1.3872 |
| rank–trace certificate, λ=1, η=0.10 | +0.49 N | +0.4867 N |
| same, η=0.05 | +0.51 N | +0.5078 N |
| same, λ=0.9 | +0.39 N | +0.3913 N |
| Cauchy–Schwarz certificate, three rows | +0.380 / +0.402 / +0.286 | +0.3786 / +0.4009 / +0.2851 |
| interior reading 4 − 2 − ratio | +0.613 | +0.6128 |
| dimension d at T = 2000, two windows | 1650 / 1834 | 1650 / 1834 exact |
| their trace-level cross-check, two windows | 1322.28 / 1436.41 | 1322.28 / 1436.41 |
Six significant figures on their own cross-check, and all five ratio rows of their second table exact to the printed three decimals. Verified Two gates were built before any of it: the closed-form window transform against two independent quadratures (worst 9.1×10⁻¹³), and their Lemma 2.2 completeness identity (worst 1.2×10⁻¹¹), with a falsifier — detune the lattice by 1% and the identity degrades to 8.1×10⁻⁷, the aliasing signature. Along the way our own reference quadrature turned out to be the thing that was wrong, not the code; it was fixed, not loosened.
Their section 8 tests the off-line mechanism on fabricated pairs, because zeta has no off-line zeros. We ran it on genuine ones: Davenport–Heilbronn, in their own window (968.377, 2031.623], 1080 on-line zeros and 58 real mirror pairs, completeness cross-checked against an independent argument-principle staircase (1195/1079 against the list's 1196/1080).
The negative index is 32, against 58 pairs. Their signature-(1,1) prediction holds on non-synthetic data for the first time. The certificate stays below the truth at every band-width tried — +464.53 against a true 1080. Measured
And the cleanest structural fact we measured: flatten every pair onto the line and the trace moves by −1.95×10⁻³ while the Frobenius norm moves by +73.763. The trace is exactly blind to off-line depth; the entire cost of being off the line lands on the second moment, 1.2718 per pair, and the certificate loses precisely that.
| Target | What it needs | Status |
|---|---|---|
| 0.67250 | the published certificate | proved, unconditional |
| any window improvement | — | closed. Montgomery–Taylor is optimal in closed form; the whole window lever is 0.005834 |
| 0.6792 on the line, 0.8396 distinct | the out-of-band positivity trick, imported | the positivity is now unconditional (below); the import is blocked by an identity, not by a gap |
| 0.68185 | their stated ceiling | kernel-checked conditional on EnclOK, and only in the Lean development — an explicit 256-periodic law with p₀ = 0.681828687… . What the kernel proves is the implication EnclOK → ceiling; EnclOK itself is checked by interval arithmetic outside Lean against a certificate file that is not published. Not reproduced here, and not refuted. The paper asserts the number with no derivation. |
| 0.6730085279 | a band-one configuration-wise certificate with a sharper certified inequality — [Ainta26] | proved, unconditional, audited on the sibling page of this site |
| 0.6731181924629976 | the same class again, with a one-step coboundary inside the certified inequality and a sharp trace credit — this programme's own nine-point result | proved, unconditional, both certificates published beside the paper |
| 0.6735128 | the ceiling of the Ψ-defect route these three all use | measured, over three routes — the route's cap, not the class's |
| 0.72222 = 13/18 | a fourth spectral moment | Hardy–Littlewood-type, conjectural |
| 0.70 / 0.80 / 0.90 | pair correlation to support 1.04 / 1.26 / 1.70 | we reproduced these thresholds: 1.0426 / 1.2578 / 1.7015 |
| 1.00 | all moments | the pair correlation conjecture |
Their Remark 1.1 says no certificate "of this kind, reading only this bandwidth-one data and holding configuration by configuration, can certify a proportion of simple zeros exceeding 0.68185". It is the sentence on which every account of the method's limits rests.
The number appears exactly once in 251 pages, and so does the phrase "explicit
extremal law"; searches over the paper, the informal note, the 95-page appendix and the 116 pages of
transcripts turn up no derivation and no cross-reference, and the certificate class is never defined
there. It is derived in the Lean repository instead, in
Zeta23/PairCeiling/, and we reconstructed the problem before reading it. Both accounts
are below, because the reconstruction is what makes the formal statement legible — and because
it is the reader of the paper who currently has no way to know the remark has a formal derivation at all.
Measured
So we reconstructed the problem the sentence must be about, and it has a clean shape. Normalise
to mean density one and let M be the limiting second moment of the multiplicities.
Integrality gives m² ≥ 2m − 1 with equality on m ∈ {1,2}, so a
configuration's true simple-on-line proportion is exactly 2 − M, and a
certificate holding configuration by configuration cannot certify more than
2 − Mmax. The ceiling is a question about how large the second
moment can be made. In configuration terms Montgomery's pair correlation at band-width one is the
single constraint ĝ(α) = M − |α| on |α| ≤ 1 for the hole
function g, with g ≤ 1 because the pair correlation is a positive
measure.
By Krein factorisation the dual collapses to
min ⟨Ah,h⟩/|⟨h,1⟩|² for the kernel
A = I + |u − v| on [−½, ½]. It returns
1.327499296320, against the paper's own optimised
1/c*₁ = 1.32749929632059 — agreement to 6×10⁻¹³,
by a route that shares no algebra with theirs. The ceiling problem's dual is the optimisation
the paper already performs. Verified
That immediately says the gap between 0.67250 and any higher ceiling is a duality gap — and the next step decides whether the gap is real. Attaining the bound would force the configuration's inter-zero gaps onto the zero set of the optimal test function. Those zeros are not an arithmetic progression; they sit at
zk = k + c/k + O(k−2), c = 0.061183 — 1.05727830, 2.03006755, 3.02024302, 4.01523565, …
Three points in a configuration force zi + zj = zk,
which needs 1/i + 1/j = 1/(i+j) — false for every positive pair. Search over all
841 triples up to index 59: the closest miss is 3.17×10⁻³, the median miss
0.0077, none below 10⁻⁶, and the closed form predicts the closest miss exactly.
Measured
So the bound 1.3274993 is attained by no configuration containing three points, hence by nothing of positive density: the ceiling of the band-one method is strictly above 0.6725007. The structural claim their remark makes is correct. What is missing is only its number.
Where is it, then? Solving the configuration side directly — a semi-infinite linear programme
in the hole function, whose value is a lower bound on Mmax by construction
— brackets it at 0.6725 < ceiling ≤ 0.694 reading band-one pair
correlation alone. Convergence is slow for a reason the argument above supplies: the extremal
configuration is atomic, and a programme over densities reaches it only in the limit.
Measured
Zeta23/PairCeiling/ carries a stability inequality (two integrations by parts)
together with its instance at an explicit 256-periodic law of marked configurations
— marks in {1,2} summing to 256, i.e. exactly the class above — giving every band-one
certificate a ceiling of 0.6818287 + 2.55×10⁻⁶(|r′(1)| + ∫|r″|).
It is a ceiling per regularity budget, and that is the precise content of “of this kind”
in the remark: the flat 0.68185 of Remark 1.1 is what the formula gives when
|r′(1)| + ∫|r″| ≲ 8.4. A certificate whose test function is rougher than that is not
capped at 0.68185 by this theorem — it is capped higher. Whether the budget is bounded for their
own C³ ramp, and for the nine-point design, is not something we can read off the statement.
The law's simple-point fraction is an exact rational,
p₀ = 10909258999421303588095230195816054408197 / 16000000000000000000000000000000000000000
= 0.681828687463831…
Checked from their own table. The 256 enclosures sit at scale
2¹⁴⁰; all 255 interior rows bracket j/256 to within one
unit in the last place, so the law's form factor is Montgomery's
F(α) = |α| sampled on the band-one grid — the constraint above, exactly.
Their row certificate max|256·S(j) − j| ≤ 3×10⁻⁴⁰ recomputes
as 1.837×10⁻⁴⁰, and their edge constant
|D(1)| ≤ 0.82395317 recomputes from the enclosures as
0.823953161. Verified
And it lands inside everything we derived independently.
M = 2 − p₀ = 1.318171313 satisfies our dual bound 1.327499296 with slack
0.009327984 — that is the duality gap, now pinned at both ends; it sits above
our continuum programme's 1.305951, which is exactly what "the extremal configuration is atomic"
predicts; and p₀ falls inside our bracket.
And here is the exact status of that number, which the earlier version of this page
got wrong. Their remark is kernel-checked conditional on EnclOK — the hypothesis
that the true form factor of their law lies in the 256 enclosures. We opened
Zeta23/PairCeiling/CeilingLaw256.lean and LawN256.lean to settle it: every ceiling
theorem there — lawN256_rows, ceiling_law256,
ceiling_law256_decimal — is universally quantified over the form factor S
and takes EnclOK LawN256.K S 0 LawN256.encl as a hypothesis. Nothing in the development
discharges it, and the file says so itself: the enclosures were produced by 70-digit interval
arithmetic outside Lean, from a certificate file cert_N256_blk_b128m.json that is not in
the repository. The law's positions and weights are not in the Lean source, only its
enclosures — so the form factor cannot be recomputed from the law here, at forty digits or at any
number of digits. Their number is right as far as anything we can check goes, the derivation is
real and formal, and what it rests on is one unpublished file. We do not call it verified, and we do
not call it wrong. Measured
Their law lives at period 256, and their files claim optimality over 256-periodic laws only — nothing anywhere says 256 is the best period. So we ran the same linear programme — same class, same grid constraint, same objective — at periods 2 through 128, by column generation over continuum positions, gated first: period 2 returns the closed form 3/4 that the problem admits by hand (primal and dual), and period 4 beats an independent brute-force grid enumeration by exactly the discretisation margin. Verified
p₀(N) = 0.750000, 0.698231, 0.686309, 0.684517, 0.686395, 0.689838, 0.692361 at N = 2, 4, 8, 16, 32, 64, 128 — against 0.681829 at their 256
The ceiling is a strong function of the period. It falls by 0.066 from period 2 to a minimum at period 16, and periods do not embed in one another — doubling a configuration zeroes every odd row of its form factor, so each refinement of the grid genuinely re-binds the adversary. The computed values rise again beyond 16; their own 256 point, sitting well below everything reachable here, is what shows that rise to be the search's power fading with dimension rather than the true curve — these are exhibited-law upper bounds, and the bias measurable at small periods grows with N. The period-128 point has since been run to convergence — 0.692361, its final pricing pass returning no improving configuration, the run checkpointed every round across repeated external machine stops — so the rise beyond 16 is not an unfinished search at 128 either. Measured
The 256 attempt, reported because this page promised it: a run at their own period ran for four days and 312 rounds of column generation, descending to 0.694872 with 13,256 columns priced, and was still descending — a few parts in 10⁵ per round, no plateau — when it died on a memory allocation the instrument cannot make at that period. We therefore read no verdict from it, by the rule this programme fixed before the run started: an unconverged value is a fact about the search, not about the curve. What it does establish is a valid upper bound, 0.694872, sitting 1.3×10⁻² above their exhibited 0.681829 — consistent with everything above, and with the measured fact that our search bias grows with the period. The column pool is checkpointed and the arithmetic that failed is a three-line fix; whether it is worth four more days is a question about the value of a number nobody needs, not about the mathematics. Measured
Two things follow. Their 0.681829 is the lowest ceiling anyone has exhibited at any period — every period we reached yields a valid but strictly weaker one. And whether periods beyond 256 descend further is open — in their files as much as here — so the number of Remark 1.1 is an upper bound on the method's true ceiling, with the limit over periods computed by nobody. One structural fact came out of the laws themselves: at every period the optimal law is a near-unit fence carrying a dilute gas of double-plus-hole defects — at period 16 a single such configuration holds 71 per cent of the mass — which is the concrete finite-period face of the atomic extremal configuration the duality argument above predicts.
A search can exhibit a law but never certify one, so the values above are upper bounds whose search bias their own 256 point exposes. We therefore solved the grid-position subclass of the same programme to proven optimality — positions restricted to a grid of spacing 1/D, pricing by exhaustive enumeration of every configuration, no seed families, no local search, nothing a basin can hide from — with both sides certified: the exhibited law's weights re-solved at 40 digits, and a weak-duality lower bound whose minimum reduced cost is taken over the entire class. At period 8 and D = 8 that pass prices 629 million configurations. Verified
period 4: 0.698593, 0.698410, 0.698255, 0.698238 at D = 16, 32, 64, 128 — the sweep says 0.698231
period 8: 0.694812, 0.690957, 0.688970 at D = 4, 6, 8 — the sweep says 0.686309
The certified values descend monotonically onto the sweep's numbers and nothing certified ever comes out below them — which is what it would have taken to show the search was fooling us. At period 4 the finest grid lands within 6×10⁻⁶ of the sweep; at period 8 the δ²-extrapolation of the three grids lands 1.1×10⁻⁴ above it. The certified optimal laws also carry the sweep's structure exactly: a unit fence plus a dilute gas of double-plus-hole defects. So the sweep instrument is validated where an exact route exists, at stated resolution — the finest certified-to-sweep gap at period 8 is 2.7×10⁻³, so a bias smaller than that remains invisible — and period 16 is beyond any exhaustive check: the mark-assignment count alone passes 10¹⁰ at the coarsest useful grid, which is why the certified ladder stops at 8 and says so.
One caveat, and it is theirs, stated in their own file. The chain runs: an exact-rational
certificate file cert_N256_blk_b128m.json → interval arithmetic performed
outside Lean → the hypothesis EnclOK that the true form factor lies in
those enclosures → everything downstream kernel-checked. That certificate file is not in
the repository — the README says it is available from the authors. Publishing it, or
the interval-arithmetic driver, closes the one link that a reader currently cannot check, and it is
the single most useful thing that could be added to the development.
Said as a reader's need rather than a complaint, and with what it would cost. The
caveat is theirs, stated in their own file, not something we found hidden. What a published
certificate buys is specific: a reader could recompute the law's form factor and discharge
EnclOK without the authors — which is the difference between a chain that is
kernel-checked and one that is checkable end to end. It is the practice this programme has adopted
for its own computer-assisted results and the one we asked of the draft audited on the sibling page:
certificates committed, every enclosure reconstructed per run, the trust base stated in the paper.
Our own two certificates are published beside the results they carry, and the verifier reconstructs
every box from them on each run rather than from an in-memory design. We are asking for exactly what
we hold ourselves to, and nothing more.
And one arithmetic coincidence, which we report as an inference and not a finding: invert their own optimal-window formula at 0.68185 and you get band-width λ = 1.0139 — 1.4% above one, the single place a band-width-one ceiling cannot come from. Exploratory
A note on scope, since the two ceilings in the paper are easy to conflate. Section 7.5(a) and
Proposition 7.4 cap the rank: n₊ ≤ d = λ₁N, giving 100% at
λ = 1. Remark 1.1 caps the extremal law over configurations. They are different
quantities and both are upper bounds; the second is the tighter, and it is the one nobody has
written down.
The route past the window is the trick of Chirre, Gonçalves and de Laat: use test functions whose transform is allowed outside the known band, provided it has the right sign there, and pay only with the sign of Montgomery's form factor rather than its value. That buys 0.6792 on the line and 0.8396 distinct — strict improvements of +0.0067 and +0.0034 over the paper's own optimum. Under the Riemann hypothesis. The paper looked at it, priced it at the bound 1.3208, and classified it "beyond the method".
The positivity it needs is unconditional, and the proof is two lines. The zero multiset is invariant under ρ ↦ 1 − ρ̄ — the functional equation, which the paper already assumes. So with A(u) = Σρ e(ρ−½)u for real u, reindexing gives Σρ e−(ρ−½)u = A(u), and for any weight W(z) = ∫ g(t)eztdt with g ≥ 0,
Σρ,ρ′ xρ−ρ′ W(ρ − ρ′) = ∫ g(t) |A(t + log x)|² dt ≥ 0, for every x > 0, unconditionally.
The convergence needs only |β − β′| < 2, which the critical strip gives. Checked independently: 1000 evaluations over 40 configurations carrying off-line pairs to depth 0.49, the minimum was +6.45 and nothing was ever negative; the closed-form identity reproduces the direct double sum to between 7 and 10 digits. Verified
So the input that made this route look conditional is not conditional. It holds for exactly the β-weighted object the certificate reads, and it comes from the functional equation, not from the hypothesis.
Establishing the positivity is not the same as running the import, and the obstruction between
them is sharper than we first thought. Written out, the certificate's second moment is
‖Â‖²F = Σρ,ρ′ Ψ(γ−γ′) with
Ψ = (Φ₂/aL)², Φ₂ the transform of
φ². So the weight it places on the form factor is the
autocorrelation of a non-negative function, and an autocorrelation is non-negative
everywhere, not merely inside the band. Measured for the paper's own C³ ramp: the weight
never goes below −1.8×10⁻¹² anywhere, and beyond the band edge it is
1.8×10⁻¹². Measured
That is fatal to the manoeuvre in one line: to drop an out-of-band piece for free you need the weight there to be negative, so that non-negativity of the form factor pays for it. With a weight that is non-negative everywhere, out-of-band mass can only be paid for, never dropped. This is a property of the square, not of the taper, so no re-choice of window reaches it, including a window of non-compact support.
There is exactly one escape that is not obviously closed, and we priced it. Replace the Frobenius
norm by a masked form ⟨Â, S∘Â⟩ with
Skl = s(τk−τl): this is still computable from
the prime side by the same diagonal method, and its weight is signed when
ŝ is. It fails, and it fails on the case that makes the whole certificate sharp. A
single simple on-line zero contributes a rank-one P = uuT with
|u| = 1, and the masked form charges it
(u∘u)TS(u∘u) instead of 1 — measured at 0.48 to 0.81 for a mild
signed mask, and negative for a stronger one, at four different zero positions. The
counterexample needs no taper at all: for P = [[1,1],[1,1]] the rank–trace lemma
demands a second moment of at least 3, the Frobenius norm gives 4, and a mask with off-diagonal
σ gives 2 + 2σ — which fails as soon as σ < ½, while any mask with a
sign change in its transform has off-diagonals below 1 by construction.
Measured
So the barrier is the conversion, not the arithmetic. A signed second-moment weight is perfectly admissible in Montgomery's classical deduction, which uses no spectrum — but that route buys "simple", not "simple on the line". The inertia step is what buys the line, and the inertia step is exactly the one that is spectral. Anyone attacking this should know that the obvious two constructions are closed for stated reasons, and that what is left is an inertia bound whose second-moment slot is not a spectral functional of the matrix at all. The prize is unchanged: +0.0067 on the line and +0.0034 distinct, unconditionally.
Reaching band-width 1.04 is where the next 2.7 points live. The paper attributes the barrier to "the Hardy–Littlewood conjectures, or equivalently Montgomery's pair correlation conjecture for α > 1". Reading section 6 against section 5, what the certificate actually consumes is a one-sided upper bound, not an asymptotic — every step that uses the second moment uses it as an inequality — on one specific signed off-diagonal prime sum, with shifts restricted to |h| ≤ T0.04, and an average over those shifts with any power saving would do. That is materially weaker than the conjectures the paper names, and the paper never draws the bound-versus-asymptotic distinction for the second moment. Exploratory We could not find a route to it from a sieve-type upper bound either, because the sum is signed and needs cancellation in the shift variable, and we record that as an open point rather than a claim.
Everything in Part I is the zero side. The other expression for their (2.20) — the integral
of φ̂φ̂ against the density built from
Γ′/Γ and the prime powers n ≤ X — is the step that makes
the whole method arithmetic, and it is the step their own internal referee calibration ranked
lowest in confidence. We built it, as one matrix product so that d² integrals cost
one multiplication, and ran it against a zero side computed by different code in a different
session.
For zeta the two expressions agree to 5.9×10⁻¹⁰ on
I = [600,1200] and 5.3×10⁻¹⁰ on [1000,2000]; and the prime side
returns tr Ĝ = 867.4206 and ‖Ĝ‖²F = 1203.2626,
which are the zeta row of the off-line control in Part I — two expressions, two code paths, the
same numbers. Verified
Then their Theorem E, which carries no table anywhere in the 251 pages. Their
section 7.3 replaces the density by νX,χ = μχ +
PX,χ, drops the pole term, and says Section 4 goes through verbatim. Run for
the real primitive characters of conductor 3, 4 and 5, on a zero census of our own built and gated
independently, the prime and zero sides agree to
1.4×10⁻¹¹ · 1.2×10⁻¹¹ · 1.5×10⁻¹¹ and the traces to ten significant figures
— better than zeta's own residual, for a reason worth having: that residual is the precision of the ordinate list, not truncation of the zero sum. Extending zeta's zero sum from 1,794 zeros to 87,702 does not move it at all, nor does quadrupling the range of the quadrature. Round our character census to nine decimals and it reads 7.2×10⁻¹⁰ against zeta's 5.9×10⁻¹⁰; the whole ladder is exactly a factor of ten per decimal. Anyone reporting an agreement figure for this identity is reporting their zero list's precision once the sum reaches about a hundred units above the top of the lattice. Measured
And the table itself shows what section 7.3 leaves invisible: conductor buys exactly
what height buys. The band-width in units of the mean spacing is
λ₁ = log(qT/2π) / ℓ1,χ, so at fixed height it rises with the
conductor, and every downstream quantity follows.
| T = 600, λ = 1, η = 0.10 | λ₁ | ‖Â‖²/tr  | predicted 1/c | certificate / N(I) |
|---|---|---|---|---|
| ζ | 0.9219 | 1.3877 | 1.3920 | +0.4461 |
| L(s,χ), conductor 3 | 0.9361 | 1.3561 | 1.3803 | +0.4788 |
| L(s,χ), conductor 4 | 0.9390 | 1.3626 | 1.3780 | +0.4732 |
| L(s,χ), conductor 5 | 0.9411 | 1.3579 | 1.3763 | +0.4800 |
The identity itself has now been rerun at heights 1500 and 1950 — two and a half and three and a quarter times the original table — and it holds at the same floor: the characters read 1.4–2.1×10⁻¹¹ at both heights with traces agreeing to twelve figures, while zeta's own floor grows with height exactly as a fixed nine-decimal ordinate list must make it. Measured
A 7.6 per cent better certificate at the same height, purely from the conductor. Extending the
censuses to height 4100 lets the law be read along both variables at once, and it splits cleanly
in two. The band-width is a function of the product qT alone, exactly: conductor 5
at height 600 and conductor 3 at height 1000 both give
λ₁ = 0.941066813, agreeing to the last bit, as do conductor 4 at 1500 and
conductor 3 at 2000. The certificate is not — at the same
λ₁ those pairs read +0.4800 against +0.5133, and +0.5328 against +0.5505.
So conductor and height are interchangeable in the band-width and only approximately in the
certificate, and the residual halves as the height rises. The measured moment ratio sits
below the asymptotic for every character, by a margin that shrinks with height (1.76 to
0.77 per cent at conductor 3 over T = 600 to 2000), while for zeta it crosses sign.
Measured
The Davenport–Heilbronn run above is a single point. We hold a one-parameter family
f(s,τ) = (1−τ)(1+√5·5−s)ζ(s) + τL(s,ψ) that satisfies one
functional equation at every τ — re-verified here to 4.1×10⁻²⁸ at 30 digits —
with Euler products at both ends and genuine off-line zeros in between. We ran their certificate
along it in their own window, I′ = (968.377, 2031.623], the window of
their table (6), so the curve extends the control above rather than sitting beside it.
Verified
The census is complete and gated three independent ways at each τ: an argument principle
round the full box, sign changes of the real Hardy-type function on the line, and every off-line
zero located and verified individually. All three agree exactly at all thirteen parameter values,
with N(I′) = 1194 or 1196 and an off-line pair count running
0, 50, 63, 84, 88, 87, 86, 80, 63, 53, 24, 14, 0. Both endpoints return exactly
zero pairs, which is the curve's own null.
Branch A. Across 39 readings — thirteen parameter values at three band-widths — the certificate never exceeds the truth, and the negative index never exceeds the true pair count. The index tracks the pair count through its whole rise and fall: at λ₁ = 0.9292 it reads 0, 27, 40, 57, 64, 64, 64, 58, 46, 35, 16, 5, 0. The registration fraction averages 0.655 there and 0.705 at λ₁ = 0.9458, against 0.101 at λ₁ = 0.7179 — strongly governed by band-width, and only weakly by how many pairs there are. Measured
Two things the single control could not show. First, the trace is blind not merely to off-line depth but to the entire configuration: across the whole family, with 88 pairs appearing and vanishing, the trace moves by 7×10⁻⁵ of itself. Second, and sharper, see below.
The two endpoints of that family are the cleanest experiment on the page. At τ = 0 and τ = 1 every zero is on the critical line — no off-line pairs at all, same window, same functional equation, 1194 against 1196 zeros. The certificate nevertheless reads +278.4 at one endpoint and +558.5 at the other, a factor of two apart; and at band-width λ₁ = 0.7179 it reads −15.2 at τ = 0, certifying nothing whatever on an object all of whose zeros are on the line, while certifying +79.0 at τ = 0.5 where there are 86 genuine off-line pairs. Measured
The cause is measured, not inferred. At τ = 0 the zero set is the union of two independent
families: the zeros of ζ, and the zeros of the factor 1+√5·5−s,
which lie on the critical line in an arithmetic ladder of spacing 2π/log 5.
Splitting the second moment into the two families and their cross term gives
1053.8 + 257.7 + 517.6: the ladder on its own is almost correlation-free, and
the cross term between two independent families is 28 per cent of the whole second
moment. Superposition destroys level repulsion, the second moment inflates, and the
certificate collapses.
So the shortfall between what the certificate reads and the truth is governed by the pair correlation of the configuration, not by how many zeros are off the line. This is a natural, non-synthetic instance of what section 1.5 of the paper says in the abstract — that the remaining third is not shown to be off the line, merely not reached — and it puts a mechanism and a number under it. At the two higher band-widths the certificate runs between 23 and 49 per cent of the truth along the family, and at the lowest it runs from below zero to 22 per cent, so the Davenport–Heilbronn run above, where the certificate reads 43 per cent of the truth and the slack is therefore 57 per cent, is not a fixed price of reading two moments; it is a functional of the configuration.
Item Three leaves the shortfall as a mechanism without a formula. It has one, and it is
parameter-free. In the paper's units the overlap of two zeros is
Φ₂(γ−γ′)/aL exactly, because the lattice samples the continuum at
critical density, so ‖Â‖²F is a two-point sum with kernel
Ψ = (Φ₂/aL)². Split a configuration into two families; if they are
independent, each zero of one sees the other at its mean density, and the cross term must take the
value
2 d₁ d₂ Lτ ∫Ψ, ∫Ψ = 2πb/a²L (Plancherel), Lτ = the length the lattice covers
with nothing fitted: a, b, L and the lattice are fixed by
the taper. Tested against zeta's ordinates superposed with an arithmetic ladder, sweeping the
ladder's spacing so the second family carries 9 to 48 per cent of the density, and averaging over
ladder phase: the law tracks the measured cross term to between 0.1 and 0.6 per cent at
every spacing tried but one, and the certificate with it.
Measured Its resolution, measured rather than assumed: a
deliberately dependent pair of families departs from the same law by 39 per cent, so the
test separates independence from dependence by a factor of about 250.
The single exception is the spacing 2π/log 5, and it is not an accident.
There the law misses by 2.9 per cent and, more tellingly, the phase spread of the cross term is
twenty times its neighbours'. Scanning spacings settles what that is: the spread is
25 to 79 at 2π/log n for n = 2, 3, 4, 5, 7, 8, 9, 11, 13, and
1.7 to 4.1 at four off-tune spacings in between. Clean separation, no overlap. A
ladder whose reciprocal spacing lands on log n is phase-locked to zeta's own
arithmetic and is not an independent family at all.
And the strength of the lock has a closed form — now complete across the
harmonics. Poisson-summing the ladder puts its m-th harmonic at frequency
m·log n carrying the second moment's own weight there, and the coefficient is
zeta's structure factor at that frequency, which by Landau is
−(T/2π)·Λ(n)/√n. The harmonics sit at distinct frequencies in
the ladder phase, so they are orthogonal over the phase average and the spread adds them in
quadrature:
sd / cross = C · √( Σm [ Ŵ(m log n)/Ŵ(0) · Λ(nm)/nm/2 ]² ), the sum over nm ≤ eL
with one constant per window and nothing else free. Measured against the nine arithmetic spacings on three windows — the original, one 2.4 times longer and one 3.9 times longer, 27 readings whose raw spreads differ by a factor of three: the ratio to that predictor is constant to 0.7–1.6 per cent within each window. The m = 1 term alone leaves the three sparsest ladders 6 to 31 per cent high — an excess that does not move when the ladders carry four times the points, so it is not finite-sample; it is the m ≥ 2 prime-power harmonics, which for n ≥ 5 sit at the band edge where the weight vanishes and for n = 2, 3, 4 are in-band with weight growing as n falls, in exactly the observed order. Adding them linearly instead of in quadrature over-corrects to a 15–17 per cent spread, which is why the naive harmonic sum looks worse than no correction at all. And the null carries its floor: at the off-tune spacings the predictor is identically zero and the measured spread is 0.0006 to 0.009, shrinking as the window grows while the arithmetic signal stays. Measured What the single constant still carries is a drift with the window, now measured at five of them: 0.2605, 0.2223, 0.2073, 0.1941, 0.1858 at T = 1000, 2500, 4000, 6000, 9000 — the law's collapse itself holds at 1.3–1.6 per cent at the two new heights, so it is the constant alone that moves. Five windows are enough to kill the flat candidate: constant (which is what the Poisson–Landau derivation says) is excluded outright — the measured C(9000)/C(1000) = 0.7134 ± 0.0043 is sixty-seven standard errors from 1 — and the least-bad two-parameter fit needs a negative limiting value, unphysical for a variance ratio. The constant carries height structure the derivation does not produce; no closed form is claimed, and that is the finding.
Earlier versions of this page also excluded a pure 1/log T, on the strength of a
predicted ratio of 0.7628. We cannot reconstruct the height convention that
gives 0.7628, and it is withdrawn. The two natural conventions give
log 1000 / log 9000 = 0.7587 (10.5σ from the measurement) and
log(1000/2π) / log(9000/2π) = 0.6976 (3.7σ). And the convention a
dyadic window actually suggests — the mid-window height 1.5T in the standard density
normalisation — gives
log(1500/2π) / log(13500/2π) = 0.7136, which sits 0.1σ
from the measured 0.7134 ± 0.0043, with no fitted parameter. So a pure 1/log of the
mid-window height is not excluded by this measurement; on the contrary it reproduces the endpoint
ratio essentially exactly. That does not restore a closed form — one ratio at one pair of
heights is a weak test, and the five-window behaviour is what the item is about — but the flat
exclusion of 1/log T printed here before was wrong, and it was wrong because of an unrecorded
height convention. Limit
Which settles item Three's configuration on its own terms. The ladder there is the zero set of
1+√5·5−s, whose spacing is exactly 2π/log 5. Its
cross term is 517.6 against the 448.0 the independent law predicts, and since the
certificate carries the second moment with a minus sign, that arithmetic lock costs
69.6 of certificate : +278.4 measured where independent families of the
same densities would give +348.0. So the slack is now a correction with a formula plus a
named residue, rather than a range — and the residue is a resonance with the prime 5.
Their synthetic run reports no negative directions at all at pair depth 10⁻⁴, which reads as a hard blind spot. It is not. We measured it three ways at fixed height and taper: with the lattice oversampled at λ₁ = 1.1, so there are more atoms than zeros, detection is not bounded by the mechanism, and the smallest eigenvalue falls off like the square of the depth — measured down to depth 0.001, which is where we stopped. At λ₁ ≤ 1 there are fewer atoms than zeros and the negative direction has nowhere to live. Measured
On genuine pairs at λ₁ = 0.9292, 37 of 58 register; the shallowest that does sits at depth 0.086137, the deepest that does not, inside the atom span, at 0.195160 — so in that band registration depends on the local configuration rather than on depth. This is the depth-resolution curve the paper leaves qualitative, and it inverts the natural reading: the blindness is not total, and it is not about the zeros.
Off-line cost saturates completely for adjacent pairs: sixty-four consecutive off-line pairs cost the second moment the same as one isolated pair. Measured A block of consecutive pairs is invisible to the second moment except at its two ends. The certificate stayed sound in every one of the twenty-four near-extremal variants we tried, so this is not a soundness threat — it is a distinct, non-arithmetic reason the two-moment reading gives up information, and it is not mentioned anywhere in the paper.
The family of item Two settles what that effect needs, exactly rather than by fitting. Write the matrix as the all-on-line control plus one perturbation per pair; then the joint cost is the sum of the isolated costs plus the cross terms, identically, and the cross terms are a three-by-three contraction that can be computed in closed form. Across 33 readings at pair counts from 14 to 88, the cross terms are at most 0.6 per cent of the joint cost: on real configurations the cost is additive. Were the saturation above operating at these counts the cross terms would be about −99 per cent, so the test resolves the effect at roughly one part in five hundred of the size it would have. The interaction is real but short-ranged — the largest cross term always sits at the closest pair of pairs, one to three mean spacings apart where the pairs are many, and is negligible beyond about four. Saturation needs adjacency, and naturally occurring off-line pairs are never adjacent.
The same decomposition prices a single pair, and shows that there is no such thing as
a cost per pair. Over 692 real pairs the isolated cost rises from 0.011 to 6.94 as the
depth runs from 0.02 to 0.40, so the figure quoted for any one object is the mean of a steep
depth law over that object's depth distribution; along the family the same mean runs from 0.28 to 1.67,
depending on the parameter and the band-width. Two populations cost essentially nothing: pairs shallower than about 0.02, and
— at any depth at all, up to 0.32 — every pair lying in the guard band outside
the span of the sampling atoms. Those zeros are subtracted in the certificate's
N(I′) term while contributing to neither moment, which is the finite-height
boundary term of the paper's (4.6) seen directly.
Everything above is at fixed height. Run zeta itself from T = 600 to T = 5000, where every zero
in range is simple and on the line so the truth is exact, and the certificate as a fraction of the
truth climbs monotonically: 0.4131, 0.4563, 0.5048, 0.5326, 0.5471, tracking
2 − 1/c(λ₁) = 0.6080 → 0.6271 as
λ₁ creeps from 0.9219 to 0.9453.
Measured
Since λ₁ → 1 gives 2 − 1/c → 2/3,
the shortfall does not close: it converges to a third of the zeros. That residual
third is Theorem A's own gap, and no amount of height removes it — which is worth stating
plainly because the finite-height numbers, rising steadily from 41 per cent towards 55, invite the
opposite reading.
Items Five and Six are conductor 5 on both objects, which is a real limitation, so we built a
Davenport–Heilbronn function at conductor 7:
f(s) = eiθL(s,χ) + e−iθL(s,χ̄)
with χ of order 6 mod 7 and θ fixed by the root number — real period-7 coefficients,
one functional equation, and no Euler product (a(2)a(3) = 0.0819 against a(6) = −1). Gates:
the functional equation at 1.6×10⁻¹⁵, the evaluator against the raw Dirichlet series at
4×10⁻¹¹. Its census on (1000, 1400] balances against the argument principle exactly
— 458 = 404 + 2×27, twenty-seven genuine off-line pairs at depths 0.09
to 0.59. And then a third, the same construction at conductor 13, χ of order
12: the same four gates (raw series 4.3×10⁻¹¹, functional equation
1.1×10⁻¹⁵), and a census on the same window balancing exactly —
497 = 455 + 2×21, twenty-one pairs at depths 0.05 to 0.42.
Verified
The additivity result carries over and is better resolved. The cross terms are 0.066 and 0.050 per cent of the joint cost, against at most 0.64 per cent at conductor 5; saturation at 27 pairs would have made them −96 per cent, so the test resolves the effect at one part in fifteen hundred. The closed-form decomposition agrees with the independent two-build route to 1.1×10⁻¹³.
And the depth law, fitted like for like, is one law across all three
conductors. Numbers of very different appearance had been on record — 10.86 at
conductor 7 against 14.1–15.1 at conductor 5 — but they came from two different
estimators on two different populations: one fitted pairs, the other fitted binned means with the
guard-band population included. Refit with one estimator, on the atom-span populations, over a
common depth window, the rates are 12.6 ± 0.1 at conductor 5 (527 pairs), 12.3–12.5
± 0.5 at conductor 7, and 13.5–13.7 ± 0.7 at conductor 13 — all
twelve pairwise comparisons agree within twice the combined errors, with the power to resolve a
difference of 1.1 to 1.75 where the apparent gap had been 3.6. What survives, uniformly: every
like-for-like rate sits 4 to 10 per cent below the 2L of their Remark 5.10,
to 5.8 standard errors where the sample is largest. The residual systematic is the window:
moving the depth window's edges moves a fitted rate by up to two units, because the law is not
exactly log-linear in depth. At that power no conductor law was resolvable either way; the
height windows two paragraphs below return to the question with three times the precision.
Measured
And the natural reading of that deficit — the asymptotic regime is not reached
yet — is now tested at height, and it fails. The conductor-13 object was censused
again on (2000, 2600] and (3400, 4000], each window balancing exactly against the argument
principle — 809 = 733 + 2×38 and 854 = 766 +
2×44 — and the same estimator run on the same depth window at all three
heights, gated by reproducing the banked numbers digit for digit first. The rate does
not move: 13.65, 13.85, 13.97 — a shift of 0.3 against a combined error of 0.7 —
while 2L grows from 15.27 to 17.72. The deficit grows from 10 to 21 per
cent of the prediction, at fourteen standard errors by the end. Were it the approach to an
asymptotic regime it would close with height; instead it grows by exactly 2L's own
growth, so at every height reachable here the rate simply does not scale with the band
length — the remark's scaling, not merely its constant, is what fails. One conductor,
three heights; no closed form for the rate's true height-dependence is claimed.
Measured
And with the height windows comes, for the first time, the power to separate the
conductors — and they separate. The same two windows were censused at conductors
5 and 7 (the same construction, χ of order 4 mod 5 and of order 6 mod 7), all four censuses
balancing exactly against the argument principle — 717 = 657 + 2×30 and
762 = 678 + 2×42 at conductor 5, 749 = 653 + 2×48 and 795 = 677 +
2×59 at conductor 7 — and the same estimator run on the same depth window,
gated once more by reproducing every banked number digit for digit first. The gate earned its
keep a second time: the window's printed edges were a rounding of the true data-defined ones,
and the rounded lower edge silently drops conductor 7's shallowest pair, moving its banked rate
by 0.4 units; the exact edges are now pinned. The rate is flat in height at every
conductor — 12.45, 12.42, 12.60 at conductor 7 and 13.65, 13.85, 13.97 at conductor 13,
while 2L grows by 2.45 — and the two flat lines are different: 1.3 to 1.4
units apart at both heights, in both band-width conventions, at about three standard errors
each time. The one-law reading above was the truth at the power then available; with
25 to 38 pairs per window the test resolves what it could not, and the rate is a
per-conductor constant — flat where the remark says it should grow, split where
a universal number would agree, and not ordered by the conductor either: 5 and 7 sit together
near 12.5 while 13 sits near 13.9, though 2L orders all three. The deficits at the
top height are 14, 24 and 21 per cent at conductors 5, 7 and 13, at up to fourteen standard
errors. Conductor 5's wobble has since been run to ground: the same-object anchor census at
T = 1000 (661 = 605 + 2×28, balanced) lands the engine's own object at 12.27 ± 0.47
— the dial family's level, so the object caveat is dead and the conductor-5 series is
same-object at every height — and a fourth window at T = 2600 (738 = 682 + 2×28,
balanced) comes out 13.12 ± 0.46, between the two levels and compatible with all of them.
What the four points support: the endpoint move (12.07 to 13.64, 2.4 standard errors) is real
but distributed across the range rather than owned by one anomalous height; a constant rate
strains against the four points without dying (p about 0.06); resolving the adjacent steps
would need roughly twice the pairs per window, so the conductor-5 series is recorded as a
direction, not a shape. What sets each conductor's constant is now the open question.
Measured
Each of these is the obvious next move from where the paper leaves off. Three are closed, and the fourth is closed as posed and reopened as a different problem — the degeneracy claim it rested on is false, and what is left is a question about fourth-moment arithmetic. A closed route is worth as much to whoever comes next as an open one, so each is stated with the thing that killed it.
The appendix's expansion of the pair-block eigenvalue near zero depth is off by a factor of two
against its own preceding line; the correct order is δ²L²/6, confirmed numerically to four digits.
And F(λ) in equation (1.3) collides with Montgomery's form factor F(α)
inside section 5 and Remark 5.10 — same symbol, two functions, adjacent lines.
Measured Neither touches a theorem.
We have now built it. The development — 329 files, 103,067 lines, Apache-2.0 — was compiled from a clean install of the pinned toolchain, Lean v4.33.0-rc2 against its pinned Mathlib, and all 9,010 compilation targets build and pass the proof-checker's kernel, with zero unproven gaps in the checked theorems. The audit below was done first, on the sources; the build turns its conclusions from a careful reading into a machine-checked fact. (The build itself needs nothing exotic: the toolchain, the Mathlib build cache, and patience with a 16 GB machine, which cannot compile the wide tail of the dependency graph and run anything else at the same time.)
| Check | Result |
|---|---|
| new axiom declarations | none. The two that a grep finds sit inside a docstring in a tactic file and are not declarations |
sorry under Zeta23/ or Solution/ | none. All 27 are in the deliberate trusted-statement files; the two hits under FromPNTPlus are comments saying the sorry'd part was not ported |
| shape of the headline statements | unconditional and correctly quantified — for every ε > 0 there is a T₀ beyond which (2/3 − ε)·N(T,2T) ≤ N₀ˢ(T,2T) |
| the definitions behind them | honest, and in the harder direction: the denominator counts with multiplicity, the numerator counts distinct on-line points; zeros and multiplicities are Mathlib's own riemannZeta and analyticOrderAt |
the optimal-window constant cMT | 1/1.327499296320588 — our own independent computation of it differs by 6×10⁻¹³, the paper's printed digits by 2×10⁻¹⁵ |
| the ceiling of Remark 1.1 | derived here and nowhere else; its data recomputes (see Part II) |
EnclOK | the one link outside the kernel: an unpublished certificate file and interval arithmetic done outside Lean |
That last row is the whole of it, and it is flagged in their own file rather than hidden. The definitions being in the harder direction is worth dwelling on: a formalisation can be sorry-free and still state something weaker than the paper, and this one does the opposite — it makes the denominator larger and the numerator smaller than a lazy reading would.
2L grows, so the deficit is not the approach to an
asymptotic regime — and the flat lines differ by conductor, near 12.5 at conductors 5 and 7
against near 13.9 at 13, without ordering by conductor size. The remark's scaling fails at
reachable heights; what sets each conductor's constant is the open question.If only one page of this is read, this is it — ordered by what we think is worth most to whoever works on the paper next.
EnclOK. The number occurs
once in 251 pages with no derivation there; the derivation is an explicit 256-periodic law in
Zeta23/PairCeiling/ with simple-point fraction p₀ = 0.681828687463831… .
Its positions and weights are not in the Lean source — only the enclosures of its form
factor are, produced outside Lean from an unpublished certificate file — so the ceiling is not
reproduced here and is not refuted here. Reconstructing the problem independently first: its dual is
the paper's own optimisation (we get 1.327499296320 against their
1/c*₁ = 1.32749929632059), that dual bound is attained by nothing — the optimal test
function's zeros sit at k + 0.0612/k and no two ever sum to a third — and the
extremal configuration must therefore be atomic. Their law is atomic, lands inside our bracket, and
its data recomputes from their own table: the 255 interior enclosures are j/256 to one
ulp, the row certificate to 1.8×10⁻⁴⁰, the edge constant to 0.823953161.
The duality gap is exactly 0.009327984, pinned at both ends. And running their
programme ourselves at every period up to 128: the ceiling moves by 0.066 with the period, their
256 value is the lowest exhibited anywhere, and the limit over periods is computed by nobody
— their number is an upper bound on the method's true ceiling. The sweep behind that
statement is itself certified now: the same programme solved exactly on grids at periods
4 and 8, every configuration enumerated and priced — 629 million at the finest — and
the certified optima descend onto the sweep's values without ever undercutting them. The one thing
missing is that a reader of the paper has no way to know any of this, and the certificate file the
chain rests on is not published.2L grows by 2.5: the deficit widens
to 14–24 per cent of the prediction, at up to fourteen standard errors, so the remark's
scaling itself fails at reachable heights. And with the height windows' power the conductors
separate for the first time: the rate is a per-conductor constant, near 12.5 at
conductors 5 and 7 and near 13.9 at 13, three standard errors apart in every reading — not
a universal number, not 2L, and not ordered by the conductor.sorry outside the deliberate trusted-statement files, statements
unconditional and correctly quantified, and the definitions pointed the harder way
— denominator with multiplicity, numerator distinct. Its optimal-window constant agrees
with our independent computation to 6×10⁻¹³. One link still sits outside the
kernel and they say so: an unpublished certificate file and interval arithmetic done outside
Lean. Publishing that file is the single most useful addition.No zero is located, excluded or constrained by anything on this page, and nothing here bears on whether the Riemann hypothesis is true. No RH claim is made or implied. Their result is a theorem about a proportion and does not touch the hypothesis either; their authors say so, and the paper's own remark is that the mechanism puts it out of reach.