/teal-sea
teal-sea / zeta-labstate of record · compiled 02 Nov 2026 · revision b09fa4dd · source

Library · hunts/flow_repair/NOTES.md

NOTES: what the probe measured

3,305 words · 351 lines · source

Run of 2026-08-07, probe.py, raw numbers in results.json. Probe language throughout: measured, observed. Everything here is the accurate regime (mpmath floats with measured cross-route defects); no enclosure claims are made anywhere in this hunt.

Correction (2026-10-19): an absolute floor in the DH weight

Added after the original record. Every section after this one is the original text, unchanged; where a number in it is affected, this section says so.

The defect. omega_dh dropped every term n·e^(−πn²x/5) below an absolute floor 10^-(P+8), P the ambient precision, so Φ_DH was exactly 0 wherever e^(−πx/5) fell below that floor: for u > ½·log(5(P+8)·ln 10/π), x = e^(2u). The same construction as hunts/jensen_clock's, where hunt #135's Davenport-Heilbronn gate found it (hunts/jensen_wedge/DH_GATE.md section 1) and its own correction fixed it (PR #315).

The fix. omega_dh now factors out the leading term and stops at 2^-(prec+20) of it, jensen_clock's corrected logic, reimplemented. The original is kept as omega_dh_absolute_floor, for the record and as a planted fault. test_flow_repair.py checks the fixed weight against a cutoff-free 30-term sum past the cut at the pair-1 flow's precision (P = 64, cut at u = 2.788, where the original returns 0) and in the moment ∫u⁸⁰⁰Φ_DH du at dps 75 (the original is off by 1.7e-4 in its log); with the original floor put back, 4 of its 5 tests go red.

What the floor reached here: no quadrature node. DHFlow evaluates the weight at dps + 10 digits, so its floor sits at 10^-(dps+18), and ends its integral where the integrand of its design box (t ≤ 0.12, |Im z| ≤ 1.6) falls to 10^-(dps+12). In every rule a published stage built, that limit stops short of the cut, by 0.0125 in u at dps 40 (validate) closing to 0.0013 at dps 165 (survey pair 9); the default rule would first cross it between dps 250 and 300. So no node ever carried weight 0. The floor did cut the n = 2 term early in a band of u, so the weights at the nodes move with the fix, by at most 5.7e-45 relative (dps 40) down to 1.7e-138 (dps 165), and H and H′ at sampled contour points by at most 9.2e-44 relative (measured with both cutoffs on every published rule).

Re-verification, 2026-10-19. Every stage that evaluates the DH weight was re-run with the fix on GitHub Actions (run 39355357518, RUNS.md) and every field compared with results.json by compare_rerun.py; the report is rerun_2026-10-19.json. 12 stage blocks, 370 fields, 364 identical as printed (only seconds is skipped), 6 changed, all in one figure.

Readers of these numbers: hunts/jensen_clock reads the pair1 and precision blocks; its three stages that were waiting on this were re-run on the re-run blocks and are unchanged (its RESULTS.md, "Correction"). hunts/li_dh_onset reads results.json only. hunts/lambda_dh_bounds builds DHFlow at dps up to 205, also short of the cut, so its flows lose no node either; its own published figures were not re-run here.

0. The instrument is telling the truth (stage validate, pair1 t=0)

1. The headline: the repair times, measured

The backward-heat flow H_t = exp(−t∂²)H₀ pulls each off-line quadruple of the Davenport–Heilbronn Ξ onto the real axis at a landing time t\, read off as the root of the pair discriminant Δ(t) = 2q₂ − q₁² (contour moments), which is analytic through the collision. Every t\ is a lower bound for Λ_DH := inf{t : H_t^{DH} real-rooted}, the DH analogue of the de Bruijn–Newman constant.

#height γy₀ = β−1/2naive y₀²/2t\* measuredshave
185.69930.30851720.04759140.044126344507.28%
2114.16330.15082700.01137490.011129587942.16%
3166.47930.07435050.00276440.002747848490.60%
4176.70250.22427630.02514580.023664731725.89%
5240.40460.36952610.06827640.0576518403515.56%
6320.87650.31954960.05105600.0446814689312.49%
7331.05030.26822310.03597180.0321781939710.55%
8366.64090.12850810.00825720.008030419202.75%
9411.79670.31587370.04988810.0426532833214.50%

(Pair 1 is the repo's pinned zero; the other seeds are Spira 1994 and Balanzario–Sánchez-Ortiz 2007, each re-polished in-tree by mp.findroot on dh_f with a unit winding check before use.)

Measured lower bound: Λ_DH ≥ 0.0576518, attained by pair 5, not the famous pair 1, which comes third. Prediction P4 (the crown goes to the height-240 quadruple, the deepest β in the surveyed list) held; the other deep pairs, 6 and 9, land at 0.0447 and 0.0427, well behind.

Scope, stated plainly: nine quadruples out of infinitely many. Nothing here bounds Λ_DH from above; the measured maximum is a floor, not the constant. And the strip confinement −1 < Re s < 2 (pinned by the |aₙ| argument in zeta/epstein.py) means every DH pair sits at y₀ < 1.5, so no single pair can land later than y₀²/2 < 1.125 under dynamics in which every zeta-like neighbour accelerates the landing, an observation about the surveyed mechanism, not a theorem about the sup.

2. Predictions, settled

3. The null control: the repair clock reads geometry, not arithmetic

For each pair, the probe censused the t = 0 zero configuration in a ±40 window (line zeros by phase-refined sign scan; total strip count by argument principle; the accounting closed exactly in all five windows: line + 2·quadruples = strip, e.g. 49 + 4 = 53 at pair 1), then integrated the bare N-body dynamics ż_k = 2Σ 1/(z_k − z_j), no Dirichlet series, no character, no conductor, just the measured starting positions, in the collision-safe variable Q = Δ/4 with the universal identity dQ/dt = 2 − 4Q Σ_a 1/((x−a)² − Q).

pairt\*_ODE (geometry only)t\*_PDE (measured flow)difference
10.044110470.04412634−0.036%
20.011128500.01112959−0.010%
30.002747810.00274785−0.002%
40.023663270.02366473−0.006%
50.057665650.05765184+0.024%

Sub-0.04% across the board, sign scattering with the truncation knobs (window ±40, density-model tail, RK4 step), the ~1% target of P3 beaten by a factor of ~25. The flow-repair time contains no information about the Davenport–Heilbronn function beyond where its zeros start. Any entire function with the same zero layout would repair on the same clock. Λ-style quantities measure a configuration's geometry; what would be special about ζ (if RH holds) is that its configuration needs no repair, a restatement of RH, not an explanation. This is docs/18's position-sensitivity lesson meeting docs/09's gate #3 on the flow axis.

4. Lesions: the refusals fire, and the newborn pair hides

5. Precision response: the standing rule, passed flat

t\* for pair 1 re-measured at dps 44 / 54 / 70 and contour node counts 96 / 192, bracket tolerance 1e-10: all four runs agree to the last digit (0.0441263445516239; spread 0.0). Δ(0) + 4y₀² = 4.7e-29 against the 50-digit pinned zero. The quantity is pinned; nothing wanders.

6. Literature position (searched 2026-08-07, two passes)

No tabulated de Bruijn–Newman-type constant for the Davenport–Heilbronn function was found; adjacent work exists (generalized Newman constants for other families; the de Bruijn/Newman/Ki–Kim–Lee flow theory; the Rodgers–Tao lower bound and the Polymath 15 upper bound for ζ; Spira's and Balanzario–Sánchez-Ortiz's DH zero computations, which supplied the survey seeds). The measurement itself: Λ_DH ≥ 0.0577, repair times for nine off-line quadruples, the null-control agreement, appears unpublished as of this search. One search is a check, not a survey; the numbers above are measurements made for this tree, with no novelty claim attached.

Calogero–Moser connection. The ODE ż_k = 2Σ 1/(z_k − z_j) that the null control integrates (§3) is formalized as a Calogero–Moser particle system for polynomial zeros under heat flow in Cuenca & McSwiggen, "The Rectangular Finite Free Heat Flow" (arXiv:2606.06859, 2026), and explored for random polynomials in Hall & Ho, "Zeros of random polynomials undergoing the heat flow" (arXiv:2308.11685, 2023). The 0.04% PDE-vs-ODE agreement measured here is an empirical instance of their universality on a non-polynomial (DH) case. Both added to references/papers.md §6.

7. The v6 census sweep: every off-line pair below height 20,000 (2026-10-30)

Added after the original record; nothing above changes. census_sweep.py lands all 1,268 off-line pairs of the Davenport-Heilbronn v6 census (Nathan Berzai, Zenodo 10.5281/zenodo.23585903, data CC-BY-4.0) with §3's null control, raw numbers in census_sweep.jsonl, pins in test_census_sweep.py. Grade: measured, a float ODE model of zero motion. §3 compared it with the measured flow only up to height 412; above that it is the model alone.

Inputs, no zero search. The window zeros come from the replayed census receipts in data/dh_v6_replay/receipts.tar.gz (25,078 on-line ordinates; window closure exact in all 400 windows, docs/46), the pairs from the census table in the same tarball, both read after the manifest hash check. Completeness is the census's own: docs/46 replayed it by one route, the author's verifier, which stops short of hardened (its section 3). The dynamics are a numpy transcription of probe.ode_null_control, same window ±40, 800 RK4 steps, mirrors, tail and t_end. The full census takes about three minutes on 16 cores (RUNS.md).

Checks. The transcription equals ode_null_control on identical inputs to 1e-13. Pairs 1 to 5 from the receipts reproduce §3's results.json landings to 5e-8 relative with identical zero counts. hunts/lambda_dh_exact landed the same model on windows it built itself (sign scan, located off-line zeros, closed strip counts): its 21 pairs below 10^4 agree to 1.1e-6 and its 19 pairs in [10^4, 2·10^4] to 4.3e-6. Deleting one on-line zero next to pair 5 moves its landing by more than 1e-4, so the reproduction check can fail.

Results. 1,268 of 1,268 pairs land, every one before its isolated y0²/2. Three pairs (v6 ranks 1266 to 1268) are flagged edge: their window runs past the census top, which lacks the zeros above 20,000. All three are shallow, landing near 0.01. The latest landings:

rankv6 nheight γy0t\*_ODEy0²/2shave
15240.400.369530.057665660.068276415.54%
21072442.530.397750.055864570.079102629.38%
324857.300.356110.050651160.063407020.12%
412531.280.346950.050339760.060188516.36%
521783.650.353000.050001680.062305019.75%

In the ODE model, the pair at height 240.4 lands last of all 1,265 non-edge pairs, so P4's winner stands across the census, extending hunts/lambda_dh_exact's census of the 21 pairs with β > 0.85 below 10^4 (its RESULTS.md section 8.6) to every depth below 20,000. The deepest pair (height 2442.53) has the largest isolated estimate and is shaved 29.4% into second place.

Added at landing, 2026-10-31. This sweep is also the antecedent that hunts/lambda_dh_exact's own P4 ("a complete census to 10^4 raises the best landing") lacked: its RESULTS.md section 8.7 grades it "unresolved, evidence against" because pairs with y0 at most 0.35 were never landed. With every v6 pair below 10^4 landed by the same ODE, the best landing is not raised. That verdict belongs to that hunt and is not changed here.

The deep predictions file. deep_predictions_v1.json (SHA-256 2616116f…) holds a crowding-model refit's predictions for the ten pairs with y0 > 0.37. Its estimand is the flow landing; it is compared here with the ODE landing, which §3 tied to the flow only below height 412.

Corrected at landing, 2026-10-31. This paragraph first called the file frozen and said it was written before any of the ten was landed. The tree supports neither. The file first appears in the commit that compares it, with no timestamp, and neither its source (source_sha256 18010823…) nor the refit's code is committed, so the predictions cannot be reproduced from here. The ODE landings of all ten pairs had been in the tree since 2026-10-26 (hunts/lambda_dh_exact/census_1e4.json and census_runs/full_r2/census_1e5.json), and this sweep reproduces them to 1e-5 (test_census_sweep.py). Read the table as a measured disagreement between a committed file and the ODE, not as a blind test of the model.

v6 nheight γy0t\*_ODEpredictedover
1072442.530.397750.05586460.05663781.38%
3025833.280.374090.04759980.04888782.71%
4428235.590.375340.04661940.04807193.12%
4658561.590.383640.04346820.049487213.85%
67011535.290.379600.04148440.047823115.28%
72512313.800.386630.04668230.04887294.69%
84814063.600.372780.04275760.04603087.66%
102816620.380.374190.04449430.04579422.92%
108117363.410.373770.03508590.045596829.96%
121119239.530.390630.04242460.048157913.51%

All ten predictions are too high and none lands inside its parameter-only 95% band: mean overprediction 9.51%, maximum 29.96%. The bands carry parameter uncertainty only, not model error, so the miss points at the model's form rather than its fitted constant (a reading, not a test).

Scope: one census to height 20,000, a measured maximum over 1,265 pairs, a floor for Λ_DH as in §1 and never its value.

Standing-checklist accounting

Disposition

Instrument (probe.py: the generic-Φ flow evaluator, contour-moment pair tracker, collision-safe ODE) retained. No claim promoted; no ledger entry. The headline numbers are measurements about a rival function's zero geometry under a classical flow; nothing here is evidence for or against RH (Littlewood, docs/08), and a measured max over nine quadruples is a floor for Λ_DH, never the value. Candidate for the spine, if anyone wants it: zeta/heatflow.py could grow a Phi-parametric entry point (the DH weight differs only in (q, gamma-shift, coefficients)), but that is a zeta/ change and belongs to a session that wants it, not to this hunt.