Current task: —

← Session 21 — Tier ii′ planning + launch: close SfBoundProp in Lean (D-064)

#901 assistant 2026-08-21 18:40:51.648 UTC · 11030→5538 tok open ⤢
thinking (13864 chars)
σ₂ from 4624 zeros = −0.0456736, and (4.10) predicts −σ₂ = 0.0461543, i.e. σ₂ = −0.0461543. The computed partial sum is −0.0456736 with a tail bound of 5.3e-4 — the difference −0.00048 is within the tail bound (barely, and the tail Σ_{γ>5096} 2Re(1/ρ²) is negative-ish... tail ≈ −Σ 2γ²/γ⁴ = −2Σ1/γ² ~ −2·log(5096)/(2π·5096)·... hmm the tail of the REAL part: 2Re(1/ρ²) = 2(1/4−γ²)/|ρ|⁴ ≈ −2/γ²: tail ≈ −2Σ_{γ>5096}1/γ² ≈ −2·log(5096)/(2π·5096) ≈ −0.00053. So partial + tail ≈ −0.0456736 − 0.00053 ≈ −0.04620 ≈ −0.0461543 within 5e-5 ✓ consistent!). So (4.10) at j=1 is CONFIRMED (σ₂ negative, matching magnitude within tail error). So the crux identity is established with certainty: **λ_n = S_∞(n) + 1 − Res₁(n), equivalently Sf(n) = Res₁(n) = Σ_{j=1}^n C(n,j)c_{j−1}** where c_k are the Laurent coefficients of −ζ′/ζ at s=1. This follows from Lagarias's (4.10) [τ_j = η_j + (−1)^jσ_{j+1} − (−1)^jδ] summed against binomial coefficients — itself the coefficient-wise form of the Hadamard product identity for ξ. AND the bound assembly becomes clean and tail-free: - Residue theorem (exact): sides_k(T) = λ_n(T) + triv + Res₁(n) ⟹ Res₁ = sides − λ_n(T) − triv. - |Sf| = |Res₁| ≤ |λ_n(T)| + |sides(T)| + |triv| — the D-037 B(n) form: |λ_n(T)| ≤ 2N(T) [bridge via hRH], |sides| ≤ [I_I..I_IV pieces], |triv| ≤ 1. NO tail term. NO prime-side η machinery needed for the BOUND (the η's/BL Thm 2 are needed only for the IDENTITY Sf = Res₁ — which is exactly the "BL Thm 2" content in our setup: it identifies the Lean-defined Sf := S_∞+1−λ_n with the contour-computable Res₁). Wait — actually, hold on. Let me reconsider what "Sf = Res₁" requires as INPUT. The derivation: Sf := S_∞+1−λ_n; the claim S_∞+1−λ_n = Res₁ ⟺ S_∞+1 = λ_n + Res₁. And the proof I found: (4.10) [τ vs η vs σ relation] summed binomially. The (4.10) relation comes from the Hadamard product: ξ′/ξ(s+1) = Σ(−1)^n σ_{n+1} s^n [zero side] = [Γ-part + L′/L-part − e/s] [archimedean+finite]. So the identity "Sf = Res₁" is EXACTLY the content of "λ_n = S_∞ − S_f + 1 with S_f = ΣC(n,j)c_{j−1}" — i.e., the BL decomposition. In Lean terms: we must prove S_∞(n) + 1 − liCoeff n = Σ_j C(n,j)c_{j−1} — an exact identity connecting: - the archimedean sum S_∞ (defined via odd-m tsum), - the zero sum liCoeff (paired), - the Laurent coefficients c_j of −ζ′/ζ. The cleanest proof route: the GENERATING FUNCTION identity. Consider F(z) := Σ_{n≥0} [S_∞(n)+1−Res₁(n)−λ_n]... hmm, or the direct coefficient route: both sides are finite binomial sums; the identity reduces (via λ_n = Σ_j (−1)^{j+1}C(n,j)σ_j and S_∞ = Σ_j C(n,j)τ_{j−1} — itself needing proof!) to the coefficient identity (4.10): τ_j = c_j + (−1)^j(σ_{j+1}−1). And (4.10) is the Laurent-coefficient reading of: ξ′/ξ(s+1) = Σ_j (−1)^j σ_{j+1} s^j [Hadamard product, zero side] ξ′/ξ(s+1) = [τ-series] + [L′/L series] + [pole terms] [definition of ξ] with −ζ′/ζ(s+1) = 1/s + Σ c_j s^j and ξ′/ξ = ζ′/ζ + 1/s + 1/(s−1) − (1/2)log π + (1/2)ψ(s/2)... (the ξ-ζ log-deriv relation from D-036 §0). So (4.10) follows from the HADAMARD PRODUCT of ξ — which brings us back to needing ξ's Hadamard product formalized (the L2-adjacent big lift) — OR the identity can be proved via the EXPLICIT FORMULA machinery that zeta-23-lean ALREADY HAS (EF_lit_zeta — the full Weil explicit formula, proved!). The power sums σ_j are moments of the zero measure; the τ/c coefficients are moments of the archimedean measure; their relation is a Weil-EF instance!! That might be the formalization shortcut: derive (4.10) for j = 0..n−1 from zeta-23-lean's EF with suitable test functions. n is arbitrary (up to N₁!) — so we need (4.10) for ALL j — the EF route gives each j as an instance; the UNIFORM-in-j derivation is the Laurent-coefficient comparison (one complex-analysis theorem), which is the Hadamard route. Hmm. OK — I now have a complete, numerically-verified picture of the mathematical content. Time to write D-064 with all of this. Let me also quickly re-examine what the B(n) side pieces require given the Sf = Res₁ form: FINAL ASSEMBLY (all exact, no tail!): 1. Sf(n) := S_∞(n) + 1 − (liCoeff n).re [existing Lean def, with liCoeff FIXED to the paired form] 2. THEOREM (crux, "BL2-identity"): Sf(n) = Res₁(n) := Σ_{j=1}^n C(n,j)·c_{j−1}, c_k = Laurent coefficients of −ζ′/ζ at 1. 3. THEOREM (residue theorem, L0): for T zero-avoiding in [√n, √n+1]: Res₁(n) = sides_k(T) − λ_n(T) − triv(n), where - sides_k(T) = (1/2πi)∮_{∂R_T} k_n(w−1)(−ζ′/ζ(w)) dw - λ_n(T) = Σ_{|γ|<T} mult·[1−(ρ/(ρ−1))ⁿ] = Σ_{|γ|<T} mult·[1−(1−1/ρ)ⁿ] (finite reindexing) - triv(n) = 1 − (2/3)ⁿ. 4. |Sf| ≤ |λ_n(T)| + |sides| + |triv| ≤ 2N(T) [bridge, hRH] + [side pieces] + 1 = B(n). And the side pieces need: kernel bounds (L1), zero-avoiding T with H(T) (L3), R-bound (L2), M₀ (L6), Z₂ ≤ 0.57 (certified), I_I decay. Now — the components and their formalization routes: A. **liCoeff fix (paired summable form)**: NEW def; prove Summable via pair bound |pair| ≤ C(n)mult/|ρ|² + zero_sum_inv_sq. The old def replaced. Downstream (Sf, SfBoundProp, LiVerifiedUpTo, lambda_pos_finite_height(_verified)) recompile unchanged (they only use liCoeff abstractly). ALSO: `LiVerifiedUpTo` and `SfBoundProp` become TRUE statements with the fix. B. **Crux identity (2)**: Sf = Res₁. Content: S_∞(n)+1−λ_n = Σ_j C(n,j)c_{j−1}. Route: define c_j as Laurent coeffs (via h(u) = u·(−ζ′/ζ(1+u)) analytic, c_j = h^{(j+1)}(0)/(j+1)!). Prove the Laurent/Hadamard identity: for each j, (−1)^j σ_{j+1} = τ_j − c_j + (−1)^j where σ_j = zero power sums (paired, absolutely convergent for j ≥ 2; j=1 conditionally), τ_j = archimedean coefficients — THIS is the Hadamard-product Laurent comparison. Formalization: this is the single hardest "new analysis" piece along with L2. Options: (i) via ξ's Hadamard product (needs entire-function product theory — heavy); (ii) via zeta-23-lean's EF_lit (Weil explicit formula already proved!) applied to test functions picking out Laurent coefficients — investigate in D-064; (iii) direct: prove both sides of "ξ′/ξ(s+1) = Σ(−1)^jσ_{j+1}s^j" and "ξ′/ξ(s+1) = τ-side + c-side" as power series identities near s=0 — the zero-side series converges for |s| < γ₁-ish; the identity of two analytic functions on a disc... Actually the cleanest classical proof: log ξ(s) = A+Bs − Σ_{n≥2} σ_n s^n/n (from the product), differentiate and compare with the ξ = Γ·ζ decomposition. The PRODUCT is the input. zeta-23-lean might have Hadamard-for-ξ?? They proved the full Weil EF — which is morally equivalent; their EF proof route (WeilEF/Contour.lean "rectangle identity / weighted argument principle") likely contains the argument-principle infrastructure we can reuse. INVENTORY ITEM for D-064. ALSO: the τ_j side: S_∞(n) = Σ_j C(n,j)τ_{j−1} with τ from digamma/ζ* — but our Lean SInfty is DEFINED as the odd-m tsum! So we ALSO need: SInfty(n) = Σ_j C(n,j)τ_{j−1} — the identity between the odd-m series and the binomial-τ form (this is Lagarias (4.11) ↔ (4.7) — the digamma partial-fraction expansion (4.19)-(4.22): τ_{j−1} = (−1)^j(1−2^{−j})ζ(j) via Hurwitz — Euler-Maclaurin-class identities; the odd-m tsum form ↔ binomial form is exactly the computation D-037 §3 does (absolute convergence swap). Formalization: finite binomial algebra + the double-sum swap ✓ tractable (tier-(i)-class work). Then λ_n = Σ_j (−1)^{j+1}C(n,j)σ_j: expand 1−(1−1/ρ)ⁿ = Σ_j(−1)^{j+1}C(n,j)ρ^{−j} and swap the finite j-sum with the zero sum (needs Σ mult/|ρ|^j summable — j=1 conditionally: the swap at j=1 needs care — pair first, then swap ✓). Assemble: S_∞+1−λ_n = Σ_j C(n,j)[τ_{j−1} − (−1)^{j+1}σ_j] + 1 = Σ_j C(n,j)[τ_{j−1} + (−1)^jσ_j] + 1; with (4.10) τ_{j−1} + (−1)^jσ_j = c_{j−1} + (−1)^j·... (4.10): τ_j = η_j + (−1)^jσ_{j+1} − (−1)^j ⟹ τ_{j−1} + (−1)^j σ_j·... set j→j−1: τ_{j−1} = c_{j−1} + (−1)^{j−1}σ_j − (−1)^{j−1} = c_{j−1} − (−1)^jσ_j + (−1)^j. So τ_{j−1}+(−1)^jσ_j = c_{j−1} + (−1)^j. Sum: Σ_j C(n,j)[c_{j−1} + (−1)^j] + 1 = Res₁(n) + Σ_j(−1)^j C(n,j) + 1 = Res₁(n) + (−1) + 1 = Res₁(n) ✓✓✓. PERFECT — the crux identity reduces EXACTLY to (4.10) + swap justifications. C. **Residue theorem (3)**: rectangle residue for J = k_n(w−1)(−ζ′/ζ) with poles: zeros (simple, residue via logDerivResidue ✓ vendored), w=1 (order n+1 — the higher-order machinery lift + the residue VALUE = Σ_j C(n,j)c_{j−1} — computable via the general residue formula: Res = (1/n!)·(d/dw)^n[(w−1)^{n+1}J] at w=1 — with (w−1)^{n+1}J = (w−1)^{n+1}k_n(w−1)·h(w)/w·... = [(w/(w−1))ⁿ−1](w−1)^{n+1}(−ζ′/ζ) = [wⁿ − (w−1)ⁿ](w−1)·(−ζ′/ζ(w))·... hmm: k_n(w−1)(w−1)^{n+1} = [(1+1/(w−1))ⁿ−1](w−1)^{n+1} = (w−1)[wⁿ − (w−1)ⁿ]·... (1+1/(w−1))ⁿ(w−1)^{n+1} = wⁿ(w−1). So (w−1)^{n+1}k_n(w−1) = wⁿ(w−1) − (w−1)^{n+1}. So J·(w−1)^{n+1} = [wⁿ(w−1) − (w−1)^{n+1}]·(−ζ′/ζ(w)) — and −ζ′/ζ = h(w)/(w−1) with h analytic, h(1) = 1: = [wⁿ − (w−1)ⁿ]·h(w). So Res_{w=1} J = (1/n!)(d/dw)ⁿ{[wⁿ − (w−1)ⁿ]h(w)}|_{w=1} = (1/n!)(d/dw)ⁿ[wⁿh(w)]|_1 − (1/n!)(d/dw)ⁿ[(w−1)ⁿh(w)]|_1. Second term: (w−1)ⁿh(w) vanishes to order n at 1; its n-th derivative at 1 = n!·h(1) = n!. So Res = (1/n!)(wⁿh(w))^{(n)}(1) − 1. And (wⁿh(w))^{(n)} = Σ_j C(n,j)(wⁿ)^{(j)}h^{(n−j)}(w): at w=1: Σ_j C(n,j) n!/(n−j)! h^{(n−j)}(1). So Res = Σ_{j=0}^n C(n,j) h^{(n−j)}(1)/(n−j)! − 1 = Σ_{j=0}^n C(n,j) c_{n−j−1}... wait h^{(k)}(1)/k! = c_{k−1} (since h(w) = (w−1)(−ζ′/ζ) = 1 + Σ_{k≥1} c_{k−1}(w−1)^k): h^{(k)}(1)/k! = c_{k−1}. So Res = Σ_{j=0}^n C(n,j)c_{(n−j)−1} − 1 = Σ_{i=0}^n C(n,n−i)c_{i−1} − 1 (i = n−j) = Σ_{i=0}^n C(n,i)c_{i−1} − 1. c_{−1} = 1: i=0 term = C(n,0)·1 = 1: so Res = 1 + Σ_{i=1}^n C(n,i)c_{i−1} − 1 = Σ_{i=1}^n C(n,i)c_{i−1} = Res₁ ✓✓✓ — the residue formula is consistent (good — no sign error), and its Lean proof is a Leibniz-rule computation. CLEAN. D. **Bridge (4)**: |λ_n(T)| ≤ 2N(T): each summand |1−(1−1/ρ)ⁿ| ≤ 2 on the line (|1−1/ρ| = 1 ⟺ Re ρ = 1/2 ✓ via hRH) — plus the WINDOW: zeros with |γ| ∈ (window around T)... T chosen zero-avoiding so λ_n(T) = Σ_{|γ|<T} and the bound |each| ≤ 2 needs |1−1/ρ| ≤ 1 for ALL zeros in range — off-line zeros have |1−1/ρ| > 1 if Re < 1/2. With hRH: all |γ| ≤ T ≤ T₀ on line ✓ so ≤ 2 ✓. Count: N(T) ≤ FAlt(T) (perimeter Trudgian). ✓. Note D-037's "4C₀ window slack" — with T zero-avoiding and λ_n(T) the strict truncation, no slack needed... their 2F(T) at T = √n+1 already absorbs. ✓ E. **Side pieces (L1, L2, L3, L6)**: as triaged. F. **Perimeter (already decided)**: Trudgian N(T)/S(T) (L4), RHUpTo (H1), LiVerifiedUpTo (H4), M₀ (L6 after E-023b). Now the full sub-structure is CLEAR and every piece has a home. The two big lifts identified: - LIFT-α: the higher-order residue machinery (rectangle residue theorem beyond simple poles + the Leibniz residue computation at w=1). Bounded engineering. - LIFT-β: the Laurent coefficient relation (4.10) — i.e., the Hadamard/EF content connecting zero power sums to the archimedean/finite Laurent coefficients. THE crux. Route decision needed (Hadamard-direct vs EF-reuse). - LIFT-γ: L2 R-bound (D-056) — machinery: BorelCaratheodory + Jensen exist; the partial-fraction route needs Hadamard — SAME obstacle as LIFT-β! Interesting: both LIFT-β and LIFT-γ classically need ξ's Hadamard product. If we formalize the Hadamard product for ξ once (as a packaged result with the zero-count convergence from perimeter), both lifts share it. ALTERNATIVE shared route: zeta-23-lean's EF machinery. Hmm — wait, actually for LIFT-γ (R-bound) there's the local Borel–Carathéodory route avoiding Hadamard; and for LIFT-β there's the EF route avoiding Hadamard. Either way, D-064 must decide. Let me check quickly whether zeta-23-lean has any Hadamard/ξ-factorization theorems, since that decides. Actually — the agent report mentioned StrongPNTPrefix has LogOfAnalyticFunction + Blaschke — enough for the local route. And for (4.10), the EF route: zeta-23-lean proved EF_lit_zeta (the full explicit formula!). The Laurent-coefficient identity (4.10) for each j is a MOMENT identity: σ_{j+1} = Σ 1/ρ^{j+1} (zero moment) = τ_j − η_j + (−1)^j (archimedean+finite moments). Such moment identities ARE explicit-formula instances with test functions designed to extract 1/ρ^{j+1} — but EF test functions need C²_c decay; extracting individual Laurent coefficients via EF is the Mellin-inversion game — nontrivial but the zeta-23 WeilEF/Contour.lean "weighted argument principle" might literally contain a version. INVENTORY for D-064 (part of ii′-0's remaining work). Alternatively — a much more elementary route to (4.10) avoiding BOTH Hadamard and EF: the identity Sf = Res₁ can be attacked DIRECTLY as a power-series identity: define Δ(n) := S_∞(n)+1−λ_n−Res₁(n). Hmm — we know Δ(1) = 0 (verified: S_f(1) = Res₁(1) = −γ... wait is that PROVABLE elementarily? Sf(1) = S_∞(1)+1−λ_1: S_∞(1) = −a₁ (exact); λ_1 = Σ1/ρ (paired); Sf(1) = 1−a₁−λ_1·... = 1−(log4π+γ)/2−λ_1. Res₁(1) = c₀ = −γ. Δ(1) = 0 ⟺ λ_1 = 1−a₁+γ = 1−(log4π+γ)/2+γ = 1 + γ/2 − (1/2)log 4π — the CLASSICAL Σ1/ρ formula (Davenport Ch. 12!). Proving λ_1 = 1+γ/2−(1/2)log 4π in Lean is itself the n=1 case of the whole thing — via the argument principle on ξ′/ξ (which zeta-23-lean's EF machinery essentially contains: EF with test function picking Σ1/ρ). So no free lunch, but the EF machinery is RIGHT THERE proved — the marginal work is the test-function instances. OK — I now have everything needed for D-064. Let me also decide the λ₂ memory-discrepancy note: with Sf = Res₁ exact, λ_2 = S_∞(2)+1−Res₁(2) = −1.874539+1+0.966885 = 0.092346. E-015's certified λ_2 (dps 1450, flint) — the archive's value: I should check E-015's CSV to confirm 0.092346 vs 0.092383. Quick check — this matters as a cross-validation of the crux identity (E-021 V1 already validated λ_n = S_∞−S_f+1 to 4.8e-25 using EXACT S_f(n) from the arithmetic side — meaning the archive's S_f values equal Res₁(n) if the identity holds!). E-021's S_f exact: λ_2 = 0.092346?? Let me check E-015's data file.
The identity web is now settled numerically. One final cross-check: the crux identity predicts λ₂ = S_∞(2)+1−Res₁(2) = 0.092346 — let me verify against the archive's certified λ₂ (E-015, dps 1450).
tool: Bash
input
{
  "command": "cd /d/Python/Quantum && head -5 research/numerics/data/E-015*.csv 2>/dev/null; ls research/numerics/data/ | grep \"E-015\"",
  "description": "Check E-015 certified lambda_2 value"
}
raw content_json (database row)
[{"type":"thinking","thinking":"σ₂ from 4624 zeros = −0.0456736, and (4.10) predicts −σ₂ = 0.0461543, i.e. σ₂ = −0.0461543. The computed partial sum is −0.0456736 with a tail bound of 5.3e-4 — the difference −0.00048 is within the tail bound (barely, and the tail Σ_{γ>5096} 2Re(1/ρ²) is negative-ish... tail ≈ −Σ 2γ²/γ⁴ = −2Σ1/γ² ~ −2·log(5096)/(2π·5096)·... hmm the tail of the REAL part: 2Re(1/ρ²) = 2(1/4−γ²)/|ρ|⁴ ≈ −2/γ²: tail ≈ −2Σ_{γ>5096}1/γ² ≈ −2·log(5096)/(2π·5096) ≈ −0.00053. So partial + tail ≈ −0.0456736 − 0.00053 ≈ −0.04620 ≈ −0.0461543 within 5e-5 ✓ consistent!). So (4.10) at j=1 is CONFIRMED (σ₂ negative, matching magnitude within tail error).\n\nSo the crux identity is established with certainty:\n\n**λ_n = S_∞(n) + 1 − Res₁(n), equivalently Sf(n) = Res₁(n) = Σ_{j=1}^n C(n,j)c_{j−1}**\n\nwhere c_k are the Laurent coefficients of −ζ′/ζ at s=1. This follows from Lagarias's (4.10) [τ_j = η_j + (−1)^jσ_{j+1} − (−1)^jδ] summed against binomial coefficients — itself the coefficient-wise form of the Hadamard product identity for ξ.\n\nAND the bound assembly becomes clean and tail-free:\n- Residue theorem (exact): sides_k(T) = λ_n(T) + triv + Res₁(n) ⟹ Res₁ = sides − λ_n(T) − triv.\n- |Sf| = |Res₁| ≤ |λ_n(T)| + |sides(T)| + |triv| — the D-037 B(n) form: |λ_n(T)| ≤ 2N(T) [bridge via hRH], |sides| ≤ [I_I..I_IV pieces], |triv| ≤ 1.\n\nNO tail term. NO prime-side η machinery needed for the BOUND (the η's/BL Thm 2 are needed only for the IDENTITY Sf = Res₁ — which is exactly the \"BL Thm 2\" content in our setup: it identifies the Lean-defined Sf := S_∞+1−λ_n with the contour-computable Res₁).\n\nWait — actually, hold on. Let me reconsider what \"Sf = Res₁\" requires as INPUT. The derivation: Sf := S_∞+1−λ_n; the claim S_∞+1−λ_n = Res₁ ⟺ S_∞+1 = λ_n + Res₁. And the proof I found: (4.10) [τ vs η vs σ relation] summed binomially. The (4.10) relation comes from the Hadamard product: ξ′/ξ(s+1) = Σ(−1)^n σ_{n+1} s^n [zero side] = [Γ-part + L′/L-part − e/s] [archimedean+finite]. So the identity \"Sf = Res₁\" is EXACTLY the content of \"λ_n = S_∞ − S_f + 1 with S_f = ΣC(n,j)c_{j−1}\" — i.e., the BL decomposition. In Lean terms: we must prove S_∞(n) + 1 − liCoeff n = Σ_j C(n,j)c_{j−1} — an exact identity connecting:\n- the archimedean sum S_∞ (defined via odd-m tsum),\n- the zero sum liCoeff (paired),\n- the Laurent coefficients c_j of −ζ′/ζ.\n\nThe cleanest proof route: the GENERATING FUNCTION identity. Consider F(z) := Σ_{n≥0} [S_∞(n)+1−Res₁(n)−λ_n]... hmm, or the direct coefficient route: both sides are finite binomial sums; the identity reduces (via λ_n = Σ_j (−1)^{j+1}C(n,j)σ_j and S_∞ = Σ_j C(n,j)τ_{j−1} — itself needing proof!) to the coefficient identity (4.10): τ_j = c_j + (−1)^j(σ_{j+1}−1). And (4.10) is the Laurent-coefficient reading of:\n\nξ′/ξ(s+1) = Σ_j (−1)^j σ_{j+1} s^j   [Hadamard product, zero side]\nξ′/ξ(s+1) = [τ-series] + [L′/L series] + [pole terms]   [definition of ξ]\n\nwith −ζ′/ζ(s+1) = 1/s + Σ c_j s^j and ξ′/ξ = ζ′/ζ + 1/s + 1/(s−1) − (1/2)log π + (1/2)ψ(s/2)... (the ξ-ζ log-deriv relation from D-036 §0). So (4.10) follows from the HADAMARD PRODUCT of ξ — which brings us back to needing ξ's Hadamard product formalized (the L2-adjacent big lift) — OR the identity can be proved via the EXPLICIT FORMULA machinery that zeta-23-lean ALREADY HAS (EF_lit_zeta — the full Weil explicit formula, proved!). The power sums σ_j are moments of the zero measure; the τ/c coefficients are moments of the archimedean measure; their relation is a Weil-EF instance!! That might be the formalization shortcut: derive (4.10) for j = 0..n−1 from zeta-23-lean's EF with suitable test functions. n is arbitrary (up to N₁!) — so we need (4.10) for ALL j — the EF route gives each j as an instance; the UNIFORM-in-j derivation is the Laurent-coefficient comparison (one complex-analysis theorem), which is the Hadamard route. Hmm.\n\nOK — I now have a complete, numerically-verified picture of the mathematical content. Time to write D-064 with all of this. Let me also quickly re-examine what the B(n) side pieces require given the Sf = Res₁ form:\n\nFINAL ASSEMBLY (all exact, no tail!):\n1. Sf(n) := S_∞(n) + 1 − (liCoeff n).re  [existing Lean def, with liCoeff FIXED to the paired form]\n2. THEOREM (crux, \"BL2-identity\"): Sf(n) = Res₁(n) := Σ_{j=1}^n C(n,j)·c_{j−1}, c_k = Laurent coefficients of −ζ′/ζ at 1.\n3. THEOREM (residue theorem, L0): for T zero-avoiding in [√n, √n+1]:\n   Res₁(n) = sides_k(T) − λ_n(T) − triv(n), where\n   - sides_k(T) = (1/2πi)∮_{∂R_T} k_n(w−1)(−ζ′/ζ(w)) dw\n   - λ_n(T) = Σ_{|γ|<T} mult·[1−(ρ/(ρ−1))ⁿ] = Σ_{|γ|<T} mult·[1−(1−1/ρ)ⁿ] (finite reindexing)\n   - triv(n) = 1 − (2/3)ⁿ.\n4. |Sf| ≤ |λ_n(T)| + |sides| + |triv| ≤ 2N(T) [bridge, hRH] + [side pieces] + 1 = B(n).\n\nAnd the side pieces need: kernel bounds (L1), zero-avoiding T with H(T) (L3), R-bound (L2), M₀ (L6), Z₂ ≤ 0.57 (certified), I_I decay.\n\nNow — the components and their formalization routes:\n\nA. **liCoeff fix (paired summable form)**: NEW def; prove Summable via pair bound |pair| ≤ C(n)mult/|ρ|² + zero_sum_inv_sq. The old def replaced. Downstream (Sf, SfBoundProp, LiVerifiedUpTo, lambda_pos_finite_height(_verified)) recompile unchanged (they only use liCoeff abstractly). ALSO: `LiVerifiedUpTo` and `SfBoundProp` become TRUE statements with the fix.\n\nB. **Crux identity (2)**: Sf = Res₁. Content: S_∞(n)+1−λ_n = Σ_j C(n,j)c_{j−1}. Route: define c_j as Laurent coeffs (via h(u) = u·(−ζ′/ζ(1+u)) analytic, c_j = h^{(j+1)}(0)/(j+1)!). Prove the Laurent/Hadamard identity: for each j, (−1)^j σ_{j+1} = τ_j − c_j + (−1)^j where σ_j = zero power sums (paired, absolutely convergent for j ≥ 2; j=1 conditionally), τ_j = archimedean coefficients — THIS is the Hadamard-product Laurent comparison. Formalization: this is the single hardest \"new analysis\" piece along with L2. Options: (i) via ξ's Hadamard product (needs entire-function product theory — heavy); (ii) via zeta-23-lean's EF_lit (Weil explicit formula already proved!) applied to test functions picking out Laurent coefficients — investigate in D-064; (iii) direct: prove both sides of \"ξ′/ξ(s+1) = Σ(−1)^jσ_{j+1}s^j\" and \"ξ′/ξ(s+1) = τ-side + c-side\" as power series identities near s=0 — the zero-side series converges for |s| < γ₁-ish; the identity of two analytic functions on a disc... Actually the cleanest classical proof: log ξ(s) = A+Bs − Σ_{n≥2} σ_n s^n/n (from the product), differentiate and compare with the ξ = Γ·ζ decomposition. The PRODUCT is the input. zeta-23-lean might have Hadamard-for-ξ?? They proved the full Weil EF — which is morally equivalent; their EF proof route (WeilEF/Contour.lean \"rectangle identity / weighted argument principle\") likely contains the argument-principle infrastructure we can reuse. INVENTORY ITEM for D-064.\n   ALSO: the τ_j side: S_∞(n) = Σ_j C(n,j)τ_{j−1} with τ from digamma/ζ* — but our Lean SInfty is DEFINED as the odd-m tsum! So we ALSO need: SInfty(n) = Σ_j C(n,j)τ_{j−1} — the identity between the odd-m series and the binomial-τ form (this is Lagarias (4.11) ↔ (4.7) — the digamma partial-fraction expansion (4.19)-(4.22): τ_{j−1} = (−1)^j(1−2^{−j})ζ(j) via Hurwitz — Euler-Maclaurin-class identities; the odd-m tsum form ↔ binomial form is exactly the computation D-037 §3 does (absolute convergence swap). Formalization: finite binomial algebra + the double-sum swap ✓ tractable (tier-(i)-class work).\n   Then λ_n = Σ_j (−1)^{j+1}C(n,j)σ_j: expand 1−(1−1/ρ)ⁿ = Σ_j(−1)^{j+1}C(n,j)ρ^{−j} and swap the finite j-sum with the zero sum (needs Σ mult/|ρ|^j summable — j=1 conditionally: the swap at j=1 needs care — pair first, then swap ✓).\n   Assemble: S_∞+1−λ_n = Σ_j C(n,j)[τ_{j−1} − (−1)^{j+1}σ_j] + 1 = Σ_j C(n,j)[τ_{j−1} + (−1)^jσ_j] + 1; with (4.10) τ_{j−1} + (−1)^jσ_j = c_{j−1} + (−1)^j·... (4.10): τ_j = η_j + (−1)^jσ_{j+1} − (−1)^j ⟹ τ_{j−1} + (−1)^j σ_j·... set j→j−1: τ_{j−1} = c_{j−1} + (−1)^{j−1}σ_j − (−1)^{j−1} = c_{j−1} − (−1)^jσ_j + (−1)^j. So τ_{j−1}+(−1)^jσ_j = c_{j−1} + (−1)^j. Sum: Σ_j C(n,j)[c_{j−1} + (−1)^j] + 1 = Res₁(n) + Σ_j(−1)^j C(n,j) + 1 = Res₁(n) + (−1) + 1 = Res₁(n) ✓✓✓. PERFECT — the crux identity reduces EXACTLY to (4.10) + swap justifications. \n\nC. **Residue theorem (3)**: rectangle residue for J = k_n(w−1)(−ζ′/ζ) with poles: zeros (simple, residue via logDerivResidue ✓ vendored), w=1 (order n+1 — the higher-order machinery lift + the residue VALUE = Σ_j C(n,j)c_{j−1} — computable via the general residue formula: Res = (1/n!)·(d/dw)^n[(w−1)^{n+1}J] at w=1 — with (w−1)^{n+1}J = (w−1)^{n+1}k_n(w−1)·h(w)/w·... = [(w/(w−1))ⁿ−1](w−1)^{n+1}(−ζ′/ζ) = [wⁿ − (w−1)ⁿ](w−1)·(−ζ′/ζ(w))·... hmm: k_n(w−1)(w−1)^{n+1} = [(1+1/(w−1))ⁿ−1](w−1)^{n+1} = (w−1)[wⁿ − (w−1)ⁿ]·... (1+1/(w−1))ⁿ(w−1)^{n+1} = wⁿ(w−1). So (w−1)^{n+1}k_n(w−1) = wⁿ(w−1) − (w−1)^{n+1}. So J·(w−1)^{n+1} = [wⁿ(w−1) − (w−1)^{n+1}]·(−ζ′/ζ(w)) — and −ζ′/ζ = h(w)/(w−1) with h analytic, h(1) = 1: = [wⁿ − (w−1)ⁿ]·h(w). So Res_{w=1} J = (1/n!)(d/dw)ⁿ{[wⁿ − (w−1)ⁿ]h(w)}|_{w=1} = (1/n!)(d/dw)ⁿ[wⁿh(w)]|_1 − (1/n!)(d/dw)ⁿ[(w−1)ⁿh(w)]|_1. Second term: (w−1)ⁿh(w) vanishes to order n at 1; its n-th derivative at 1 = n!·h(1) = n!. So Res = (1/n!)(wⁿh(w))^{(n)}(1) − 1. And (wⁿh(w))^{(n)} = Σ_j C(n,j)(wⁿ)^{(j)}h^{(n−j)}(w): at w=1: Σ_j C(n,j) n!/(n−j)! h^{(n−j)}(1). So Res = Σ_{j=0}^n C(n,j) h^{(n−j)}(1)/(n−j)! − 1 = Σ_{j=0}^n C(n,j) c_{n−j−1}... wait h^{(k)}(1)/k! = c_{k−1} (since h(w) = (w−1)(−ζ′/ζ) = 1 + Σ_{k≥1} c_{k−1}(w−1)^k): h^{(k)}(1)/k! = c_{k−1}. So Res = Σ_{j=0}^n C(n,j)c_{(n−j)−1} − 1 = Σ_{i=0}^n C(n,n−i)c_{i−1} − 1 (i = n−j) = Σ_{i=0}^n C(n,i)c_{i−1} − 1. c_{−1} = 1: i=0 term = C(n,0)·1 = 1: so Res = 1 + Σ_{i=1}^n C(n,i)c_{i−1} − 1 = Σ_{i=1}^n C(n,i)c_{i−1} = Res₁ ✓✓✓ — the residue formula is consistent (good — no sign error), and its Lean proof is a Leibniz-rule computation. CLEAN.\n\nD. **Bridge (4)**: |λ_n(T)| ≤ 2N(T): each summand |1−(1−1/ρ)ⁿ| ≤ 2 on the line (|1−1/ρ| = 1 ⟺ Re ρ = 1/2 ✓ via hRH) — plus the WINDOW: zeros with |γ| ∈ (window around T)... T chosen zero-avoiding so λ_n(T) = Σ_{|γ|<T} and the bound |each| ≤ 2 needs |1−1/ρ| ≤ 1 for ALL zeros in range — off-line zeros have |1−1/ρ| > 1 if Re < 1/2. With hRH: all |γ| ≤ T ≤ T₀ on line ✓ so ≤ 2 ✓. Count: N(T) ≤ FAlt(T) (perimeter Trudgian). ✓. Note D-037's \"4C₀ window slack\" — with T zero-avoiding and λ_n(T) the strict truncation, no slack needed... their 2F(T) at T = √n+1 already absorbs. ✓\n\nE. **Side pieces (L1, L2, L3, L6)**: as triaged.\n\nF. **Perimeter (already decided)**: Trudgian N(T)/S(T) (L4), RHUpTo (H1), LiVerifiedUpTo (H4), M₀ (L6 after E-023b).\n\nNow the full sub-structure is CLEAR and every piece has a home. The two big lifts identified:\n- LIFT-α: the higher-order residue machinery (rectangle residue theorem beyond simple poles + the Leibniz residue computation at w=1). Bounded engineering.\n- LIFT-β: the Laurent coefficient relation (4.10) — i.e., the Hadamard/EF content connecting zero power sums to the archimedean/finite Laurent coefficients. THE crux. Route decision needed (Hadamard-direct vs EF-reuse).\n- LIFT-γ: L2 R-bound (D-056) — machinery: BorelCaratheodory + Jensen exist; the partial-fraction route needs Hadamard — SAME obstacle as LIFT-β! Interesting: both LIFT-β and LIFT-γ classically need ξ's Hadamard product. If we formalize the Hadamard product for ξ once (as a packaged result with the zero-count convergence from perimeter), both lifts share it. ALTERNATIVE shared route: zeta-23-lean's EF machinery.\n\nHmm — wait, actually for LIFT-γ (R-bound) there's the local Borel–Carathéodory route avoiding Hadamard; and for LIFT-β there's the EF route avoiding Hadamard. Either way, D-064 must decide. Let me check quickly whether zeta-23-lean has any Hadamard/ξ-factorization theorems, since that decides. Actually — the agent report mentioned StrongPNTPrefix has LogOfAnalyticFunction + Blaschke — enough for the local route. And for (4.10), the EF route: zeta-23-lean proved EF_lit_zeta (the full explicit formula!). The Laurent-coefficient identity (4.10) for each j is a MOMENT identity: σ_{j+1} = Σ 1/ρ^{j+1} (zero moment) = τ_j − η_j + (−1)^j (archimedean+finite moments). Such moment identities ARE explicit-formula instances with test functions designed to extract 1/ρ^{j+1} — but EF test functions need C²_c decay; extracting individual Laurent coefficients via EF is the Mellin-inversion game — nontrivial but the zeta-23 WeilEF/Contour.lean \"weighted argument principle\" might literally contain a version. INVENTORY for D-064 (part of ii′-0's remaining work).\n\nAlternatively — a much more elementary route to (4.10) avoiding BOTH Hadamard and EF: the identity Sf = Res₁ can be attacked DIRECTLY as a power-series identity: define Δ(n) := S_∞(n)+1−λ_n−Res₁(n). Hmm — we know Δ(1) = 0 (verified: S_f(1) = Res₁(1) = −γ... wait is that PROVABLE elementarily? Sf(1) = S_∞(1)+1−λ_1: S_∞(1) = −a₁ (exact); λ_1 = Σ1/ρ (paired); Sf(1) = 1−a₁−λ_1·... = 1−(log4π+γ)/2−λ_1. Res₁(1) = c₀ = −γ. Δ(1) = 0 ⟺ λ_1 = 1−a₁+γ = 1−(log4π+γ)/2+γ = 1 + γ/2 − (1/2)log 4π — the CLASSICAL Σ1/ρ formula (Davenport Ch. 12!). Proving λ_1 = 1+γ/2−(1/2)log 4π in Lean is itself the n=1 case of the whole thing — via the argument principle on ξ′/ξ (which zeta-23-lean's EF machinery essentially contains: EF with test function picking Σ1/ρ). So no free lunch, but the EF machinery is RIGHT THERE proved — the marginal work is the test-function instances.\n\nOK — I now have everything needed for D-064. Let me also decide the λ₂ memory-discrepancy note: with Sf = Res₁ exact, λ_2 = S_∞(2)+1−Res₁(2) = −1.874539+1+0.966885 = 0.092346. E-015's certified λ_2 (dps 1450, flint) — the archive's value: I should check E-015's CSV to confirm 0.092346 vs 0.092383. Quick check — this matters as a cross-validation of the crux identity (E-021 V1 already validated λ_n = S_∞−S_f+1 to 4.8e-25 using EXACT S_f(n) from the arithmetic side — meaning the archive's S_f values equal Res₁(n) if the identity holds!). E-021's S_f exact: λ_2 = 0.092346?? Let me check E-015's data file.","signature":""},{"type":"text","text":"The identity web is now settled numerically. One final cross-check: the crux identity predicts λ₂ = S_∞(2)+1−Res₁(2) = 0.092346 — let me verify against the archive's certified λ₂ (E-015, dps 1450)."},{"type":"tool_use","id":"Bash_51","name":"Bash","input":{"command":"cd /d/Python/Quantum && head -5 research/numerics/data/E-015*.csv 2>/dev/null; ls research/numerics/data/ | grep \"E-015\"","description":"Check E-015 certified lambda_2 value"},"extra_content":null}]