← Session 21 — Tier ii′ planning + launch: close SfBoundProp in Lean (D-064)
thinking (13901 chars)
This is the full BL Thm 2. The decomposition λ_n = [prime part] + [archimedean] + [trivial-zeros part]. In the BL/S_f notation:
λ_n = S_∞(n) − S_f(n) + 1? Let me map: from the formula:
λ_n = −[prime part with signs] + 1 − (log 4π+γ)n/2 − Σ_{j=2}^n (−1)^{j−1}C(n,j)(1−2^{−j})ζ(j)
Hmm, and D-037's S_∞: S_∞(n) = Σ_{j=1}^n (−1)^j C(n,j) a_j with a₁ = (log 4π+γ)/2, a_j = (1−2^{−j})ζ(j). Let me expand: Σ_{j=1}^n (−1)^j C(n,j) a_j = −n·a₁ + Σ_{j=2}^n (−1)^j C(n,j)(1−2^{−j})ζ(j) = −n(log4π+γ)/2 + Σ_{j=2}(−1)^j C(n,j)(1−2^{−j})ζ(j).
BL's λ_n has: 1 − (log4π+γ)n/2 − Σ_{j=2}^n (−1)^{j−1}C(n,j)(1−2^{−j})ζ(j) + [prime part] = 1 + S_∞(n) + [prime part with the −(−1)^{j−1} = (−1)^j sign flip ✓].
So: λ_n = S_∞(n) + 1 + [prime part], i.e. **S_f(n) = −[prime part]** where [prime part] = −Σ_{j=1}^n (−1)^{j−1}C(n,j)η_{j-1} = Σ_j (−1)^j C(n,j)η_{j−1}... with η the regularized prime sums (Stieltjes-analogues). So:
**S_f(n) = Σ_{j=1}^n (−1)^{j−1}C(n,j) η_{j−1}** (up to sign convention) — the "finite part" is the PRIME-side sum, NOT a zero sum!! That's the whole point of the BL decomposition: S_∞ = archimedean, S_f = finite primes. And |S_f(n)| = 13.52 max measured (E-021, arithmetic side) ✓ consistent with Maślanka.
And how does Lagarias's Thm 6.1 bound S_f? Via the contour: the ZERO-side representation of S_f: since λ_n = Σ_ρ[1−(1−1/ρ)ⁿ] (zero side) and λ_n = S_∞ − S_f + 1 (arithmetic side), we get S_f = S_∞ + 1 − Σ_ρ[1−(1−1/ρ)ⁿ]. The contour argument computes the truncated zero sum:
λ_n(T) := Σ_{|γ|≤T}[1−(1−1/ρ)ⁿ]... and Thm 6.1 says S_f(n) = λ_n(√n) + O(√n log n)?? Hmm wait — that means S_∞ + 1 − λ_n = λ_n(√n) + O — i.e. S_∞ + 1 − [λ_n − λ_n(√n)] = O... i.e. the identity is:
S_f(n) − λ_n(√n) = S_∞(n) + 1 − λ_n + λ_n(√n) = S_∞(n) + 1 − [tail beyond √n].
And Lagarias's contour bounds THIS: S_f(n) − λ_n(T) for T ~ √n. Let me recompute with my residue bookkeeping:
S_f − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T)... no wait: S_f = S_∞ + 1 − λ_n so S_f − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T). Hmm that's not what I want either. Let me recompute from D-032's statement: "S_f(n) = λ_n(√n) + O(√n log n)". Since |S_f| is small (13.52) and |λ_n(√n)| ≤ 2N(√n) ~ √n log n — the theorem S_f = λ_n(√n) + O(√n log n) is consistent with both being O(√n log n).
So the contour's job: bound S_f(n) − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T)... hmm, no: = S_∞ + 1 − (λ_n − λ_n(T))·... S_f − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T). With λ_n = λ_n(T) + tail: S_f − λ_n(T) = S_∞ + 1 − 2λ_n(T) − tail?? That doesn't look right either. Hmm, unless λ_n(T) in Lagarias's notation means the truncated sum of the OTHER form: λ_n(T) = Σ_{|γ|≤T}[1−(ρ/(ρ−1))ⁿ]?
Let me test at n=1, T=2: λ_1(2) = 0 (no zeros below 2). S_f(1) = −0.577. So "S_f = λ_n(√n) + O" would read −0.577 = 0 + O — the O-term must BE −0.577 at n=1. The O-terms at n=1: sides + triv + Res₁ + ... : from my earlier bookkeeping: S_f = bracket_A − sides − tail with bracket_A = S_∞+1+Res₁+triv. S_f(1) = −0.577, bracket_A(1) = −0.798, sides+tail = −0.2208 (computed). OK whatever — the identity is exact but opaque.
STOP. The numerology is settled enough for the DESIGN: the exact identity structure is:
(ID) Sf(n) = λ_n(T) + [bracket_A(n) − λ_n(T) − sides(T) − tail(T)]...
No. Let me just carefully redo the algebra ONE more time, cleanly, with all signs, because this determines the skeleton.
Definitions:
- λ_n := full paired zero sum Σ_ρ[1−(1−1/ρ)ⁿ] (this is liCoeffP).
- Sf := S_∞ + 1 − λ_n (Lean def, matches BL's λ_n = S_∞ − S_f + 1).
- λ_n(T) := Σ_{|γ|≤T}′ mult·[1−(1−1/ρ)ⁿ] (truncated).
- tail(T) := λ_n − λ_n(T) (paired tail).
- Contour: integrand J(w) = k_n(w−1)·(−ζ′/ζ(w)), rectangle σ ∈ [−5/2, 2√n], |Im| ≤ T, T ∈ [√n, √n+1] zero-avoiding.
- Residue theorem: sides(T) := (1/2πi)∮ J = Σ_{|γ|<T} mult·Res_ρ + Res_{−2} + Res₁(n).
- Res_ρ J = k_n(ρ−1)·(−mult... for a zero of order m: −ζ′/ζ has residue −m). So Res = −m·k_n(ρ−1) = −m[(ρ/(ρ−1))ⁿ − 1] = m[1 − (ρ/(ρ−1))ⁿ].
- So Σ_{|γ|<T} mult·[1−(ρ/(ρ−1))ⁿ] = sides − Res_{−2} − Res₁. ... (∗)
- KEY REINDEXING: for the TRUNCATED sum (finite!): Σ_{|γ|<T} mult[1−(ρ/(ρ−1))ⁿ] = Σ_{|γ|<T} mult[1−(1−1/ρ)ⁿ] = λ_n(T) (the finite sum reindexed termwise via ρ ↦ 1−ρ, which preserves the truncated set and mult). EXACT, no convergence issue. ✓
So (∗) becomes: **λ_n(T) = sides(T) − Res_{−2} − Res₁(n)**.
Then: Sf = S_∞ + 1 − λ_n = S_∞ + 1 − λ_n(T) − tail(T) = S_∞ + 1 − [sides − Res_{−2} − Res₁] − tail.
|Sf| ≤ |λ_n(T)| + |S_∞ + 1 + Res_{−2} + Res₁ − sides − tail|?? No — the clean split for the BOUND:
|Sf| ≤ |λ_n(T)| + |Sf − λ_n(T)|, and Sf − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T) = S_∞ + 1 − 2λ_n(T) − tail — ugh, that has 2λ_n(T) ~ 2·2N(T) — no good.
The RIGHT split (Lagarias's): Sf = S_∞ + 1 − λ_n = [S_∞ + 1 + Res₁ + Res_{−2} − sides] + [sides − Res_{−2} − Res₁] − tail... circular. OK let me just substitute once:
Sf = S_∞ + 1 + Res₁ + Res_{−2} − sides − tail.
Now |Sf| ≤ |sides| + |tail| + |S_∞ + 1 + Res₁ + Res_{−2}|.
For B(n) to work: need |S_∞ + 1 + Res₁ + Res_{−2}| bounded by... 2N(T)?? At n=1: |−1.554120+1+(−0.577216)+(1/3)| = 0.798. 2N(2) = 0. Fails at n=1 but the range is n ≥ 10⁵. Hmm, but structurally "2N(T) ~ √n log n" vs "S_∞+1+Res₁+triv ~ ?". What IS the size of bracket_A(n) = S_∞+1+Res₁+triv at large n? If S_∞ + Res₁ cancels the ½n log n (which it must, since...), bracket_A could be O(n)-ish or smaller. Numerically: n=1: −0.798, n=2: −1.213. Let me estimate n=3 quickly? Tedious. The point for the DESIGN: the identity
(★★) S_∞(n) + 1 + Res₁(n) + Res_{−2}(n) = 0 for all n
would give the clean form |Sf| ≤ |sides| + |tail|, but it's numerically FALSE (−0.798 at n=1). So that's not the identity. INSTEAD, the correct identity must be:
S_∞(n) + 1 + Res₁(n) + Res_{−2}(n) = −(something of size ~λ_n(T))?? At n=1: −0.798 = ? λ_1(2)=0. Hmm. Or maybe my Res₁ is wrong by a sign or the trivial residue is different. Let me recompute Res_{−2}: integrand J(w) = k_n(w−1)(−ζ′/ζ(w)) at w=−2: −ζ′/ζ has residue −1 (simple zero), k_n(−3) = (1−1/3)ⁿ−1 = (2/3)ⁿ−1. Res_{−2} = k_n(−3)·(−1) = 1−(2/3)ⁿ ✓ (matches triv). At n=1: 1/3 ✓.
Hmm what if the w=1 pole's order and residue: k_n(w−1) = (1+1/(w−1))ⁿ − 1 = (w/(w−1))ⁿ − 1. At w=1: pole order n ✓. −ζ′/ζ ~ +1/(w−1) — WAIT: ζ has a simple pole at w=1 with residue +1, so ζ(w) ~ 1/(w−1), ζ′(w) ~ −1/(w−1)², ζ′/ζ ~ −1/(w−1). So −ζ′/ζ ~ +1/(w−1) ✓ c_{−1} = +1 ✓.
So Res₁(n) = Σ_{j=1}^n C(n,j)c_{j−1}, c_{−1}=1, c_0=−γ, c_1 = γ₁+γ². At n=1: c_0 = −γ ✓.
Let me test (★★) variant: maybe S_∞ + Res₁ + Res_{−2} = −1 + λ_n-related? Compute S_f(1) directly via prime side: S_f(1) = η_0·(−1)^{?}: S_f(1) = Σ_{j=1}^1 (−1)^{j−1}C(1,j)η_{j−1} = η_0. And −ζ′/ζ(s+1) = 1/s + Σ η_n s^n (from Theorem-2.txt line 59): so η_0 = c_0 evaluated via s: −ζ′/ζ(s+1) at s=0: 1/s + η_0 + ...; and c_0 (constant term of −ζ′/ζ at s=1) = η_0 ✓!! So S_f(1) = η_0 = c_0 = Res₁(1) ✓✓ (consistent with my earlier finding S_f(1) = Res₁(1) = −γ).
n=2: S_f(2) = Σ_{j=1}^2 (−1)^{j−1}C(2,j)η_{j−1} = 2η_0 − η_1 = −2γ − η_1. Res₁(2) = 2c_0 + c_1 = −2γ + (γ₁+γ²). And η_1 vs c_1: from −ζ′/ζ(s+1) = 1/s + η_0 + η_1 s + ...: with u = s: −ζ′/ζ(1+u) = u⁻¹ + c_0 + c_1 u + ...: so η_1 = c_1 = γ₁+γ². So S_f(2) = −2γ − (γ₁+γ²) = −1.154432 − 0.260361 = −1.414793?? But earlier I computed S_f(2) = S_∞(2)+1−λ_2 = −0.966922. CONTRADICTION! So my prime-side formula sign is wrong. From Theorem-2.txt: λ_n = −Σ_j (−1)^{j−1}C(n,j)η_{j−1} + 1 − (log4π+γ)n/2 − Σ_{j=2}(−1)^{j−1}C(n,j)(1−2^{−j})ζ(j).
n=2: λ_2 = −[2η_0 − η_1] + 1 − (log4π+γ) − [−(−1)¹C(2,2)(1−1/4)ζ(2)] = −2η_0 + η_1 + 1 − log4π − γ + (3/4)ζ(2)·(−1)^{...}: careful: −Σ_{j=2}^2 (−1)^{j−1}C(n,j)a_j = −(−1)¹·1·a₂ = +a₂ = (3/4)ζ(2) = 1.233701. So λ_2 = −2(−0.577216) + 0.260361 + 1 − 3.108240 + 1.233701 = 1.154432 + 0.260361 + 1 − 3.108240 + 1.233701 = 0.540254?? But λ_2 = 0.092383. WRONG — so η_1 ≠ c_1. Indeed η_1 is the PRIME-side regularized sum (with Λ(m)), NOT the Stieltjes polynomial! Theorem-2.txt line 53: η_n defined via Λ(m) sums (prime side); line 61: each η_n is a polynomial in Stieltjes constants (γ₀s + (γ₁−½γ₀²)s²...) — those are the coefficients of log(sζ(s+1))'s expansion — a DIFFERENT relation. The expansion at line 59: −ζ′/ζ(s+1) = 1/s + Σ η_n s^n. Hmm — but −ζ′/ζ(s+1) near s=0: 1/s + c_0 + c_1 s... where c_k are the Laurent coeffs of −ζ′/ζ at 1. So η_n = c_n EXACTLY?? But then λ_2 computation gives 0.540 ≠ 0.0924. So one of my c's is wrong. Let me recompute c_1: −ζ′/ζ(1+u) = u⁻¹ − γ + (γ₁+γ²)u + ...? From log ζ(1+u) = −log u + γu − (γ₁+γ²)u²/2 + O(u³) — let me redo: ζ(1+u) = u⁻¹ + Σ_{k≥0} (−1)^k γ_k/k! u^k = u⁻¹(1 + γu − γ₁u²/2 + γ₂u³/6 − ...). log ζ = −log u + log(1 + v) where v = γu − γ₁u²/2 + ...: log(1+v) = v − v²/2 + ... = γu − γ₁u²/2 − γ²u²/2 + O(u³) = γu − (γ₁+γ²)u²/2 + O(u³). d/du: ζ′/ζ(1+u) = −u⁻¹ + γ − (γ₁+γ²)u + O(u²). −ζ′/ζ(1+u) = u⁻¹ − γ + (γ₁+γ²)u + O(u²). So c_0 = −γ, c_1 = γ₁ + γ² = −0.072816 + 0.333177 = 0.260361. Then λ_2 from the formula: −[2c_0 − c_1] + 1 − (log4π+γ) + a₂ = −2(−0.577216) + 0.260361 + 1 − 3.108240 + 1.233701 = 0.540254 ≠ 0.092383. So EITHER the formula's η's aren't the c's, or the archimedean sign conventions differ. Theorem-2.txt line 57-61 says η_n are the Λ-sum regularized limits AND equal the Laurent coefficients of −ζ′/ζ(s+1). Both can be true (PNT connection). Then λ_2 = 0.540?? But the true λ_2 = 0.0924. So the formula as transcribed in Theorem-2.txt must have a sign convention I'm misreading (e.g. the archimedean term is 1 − (log 4π + γ)n/2 with a DIFFERENT sign structure, or the a_j signs). At n=1: λ_1 = −η_0 + 1 − (log4π+γ)/2 = 0.577216 + 1 − 1.554120 = 0.023096 ✓✓ CORRECT. n=2: −[2η_0 − η_1] = −2(−0.577216) + η_1·(+1)... wait: −Σ_{j=1}^2 (−1)^{j−1}C(2,j)η_{j−1} = −[(−1)⁰C(2,1)η_0 + (−1)¹C(2,2)η_1] = −[2η_0 − η_1] = −2η_0 + η_1 = 1.154432 + 0.260361 = 1.414793. Archimedean+trivial: 1 − (log4π+γ)·2/2 − (−1)¹C(2,2)a₂ = 1 − 3.108240 + a₂ = 1 − 3.108240 + 1.233701 = −0.874539. λ_2 = 1.414793 − 0.874539 = 0.540254. True λ_2 = 0.092383. So η_1 must be ≈ 0.260361 − (0.540254−0.092383) = −0.187510. So the true η_1 ≈ −0.1875, NOT c_1 = 0.2604. So Theorem-2.txt's "η_n = Laurent coeffs of −ζ′/ζ(s+1)" is imprecise/wrong in sign or indexing — indeed line 61 gives the polynomial structure: γ₀s + (γ₁−½γ₀²)s²... suggesting the η's relate to log(sζ(s+1)) = log(1 + γs − ...) whose expansion gives γ₀ s + (−γ₁/2 − γ₀²/2)s²... hmm: log(sζ(s+1)): sζ(s+1) = 1 + γs − γ₁s²/2 + ...: log = γs − (γ₁+γ²)s²/2 + ... So the coefficients of log(sζ(s+1)) are γ, −(γ₁+γ²)/2, ... and η_1 (true) ≈ −0.1875. Hmm: −(γ₁+γ²)/2 = −0.130181. Not −0.1875 either. What IS η_1? True value: η_1 = lim{Σ_{m≤x} Λ(m)(log m)/m − (log x)²/2} ≈ −0.1875? Hmm I recall η_0 = −γ ≈ −0.5772? Hmm no — actually the standard values: η_0 = −0.5772 (= −γ ✓), η_1 = −0.18754...? Let me check: yes I believe the BL η constants are η_0 = −γ, η_1 ≈ −0.1875, η_2 ≈ −0.21..., — wait those are related to but not equal to Stieltjes. The relation: η_n = (−1)^n[...Stieltjes polynomial...]: η_1 = −(γ₁ + γ²/2)? = −(−0.072816 + 0.166589) = −0.093773. ×2 = −0.187546 ✓!! So η_1 = −2(γ₁ + γ²/2)... hmm so η_1 = −2γ₁ − γ² = 0.145631 − 0.333177 = −0.187546 ✓✓. So the relation is η_1 = −(γ² + 2γ₁), consistent with the line-61 polynomial pattern "γ₀s + (γ₁ − ½γ₀²)s²" up to signs/conventions I'm not going to nail by hand.
CONCLUSION for the design: the exact identity content is:
- **BL Thm 2 (P1)**: λ_n = S_∞(n) − S_f(n) + 1 where S_f(n) = Σ_j (−1)^{j−1}C(n,j)η_{j−1}·(−1)... (prime side) — with the η's being the regularized Λ-sums.
- The CONTOUR identity (∗): λ_n(T) = sides(T) − Res_{−2} − Res₁(n), where Res₁(n) = Σ_j C(n,j)c_{j−1} with c's = Laurent coefficients of −ζ′/ζ at 1.
- The BRIDGE between them: S_∞ + 1 − λ_n = S_f ⟹ S_f(n) = S_∞+1+Res₁+Res_{−2} − sides − tail — and the identification needed: **S_∞(n) + 1 + Res₁(n) + Res_{−2}(n) = S_f(n) + [something]** — NO wait. We have TWO expressions for S_f: prime side (definitional) and the contour one. The assembly for the BOUND: |S_f(n)| ≤ |λ_n(T)| + |S_f(n) − λ_n(T)|. And S_f(n) − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T) — hmm this has λ_n (full) which we don't control...
Hmm hmm. OK here's the cleaner classical assembly (this must be what Lagarias does): DON'T split via |Sf| ≤ |λ_n(T)| + |Sf − λ_n(T)|. Instead substitute (∗) into Sf's definition:
Sf = S_∞ + 1 − λ_n = S_∞ + 1 − λ_n(T) − tail(T)
= S_∞ + 1 − [sides − Res_{−2} − Res₁] − tail
= [S_∞ + 1 + Res₁ + Res_{−2}] − sides − tail.
Now the bracket: S_∞(n) + 1 + Res₁(n) + Res_{−2}(n). Claim: this bracket is EXACTLY 0 for all n ≥ 1?? Numerically: n=1: −1.554120 + 1 − 0.577216 + 0.333333 = −0.798003 ≠ 0. FALSE. Damn.
Alternative: the bracket equals −2·[trivial-zero sum]?? The trivial zeros beyond −2 (−4, −6, ...) — the rectangle only encloses −2, but the FULL trivial-zero contribution to the archimedean side is the a_j = (1−2^{−j})ζ(j) sum (per Theorem-2.txt's anatomy: "the (1−2^{−j})ζ(j) terms account for the trivial zeros"). S_∞ contains Σ_j(−1)^j C(n,j)a_j — the trivial zeros are IN S_∞'s definition. So the bracket S_∞+1+Res₁+Res_{−2} mixes: S_∞ (archimedean + trivial in BL form), Res₁ (the contour's w=1 output), Res_{−2} (the one enclosed trivial zero). For the bracket to be manageable, the identity must be:
S_∞(n) + Res₁(n) + 1 = −Res_{−2}(n) + bracket₀(n) with bracket₀ small?
Numerically: S_∞+Res₁+1 at n=1: −1.554120 − 0.577216 + 1 = −1.131336. −Res_{−2} = −1/3. bracket₀ = −1.131336 + 0.333333 = −0.798003. Still −0.798. What's −0.798 at n=1? Hmm: −γ + ...? −0.798 = −γ − 0.2208 = c_0 − 0.2208. And 0.2208 ≈ ? log 2 − 0.47...? (2/3)? Meh.
You know what — I've spent enough on hand-numerology. The DEFINITIVE source is Lagarias's paper §6 (in our archive as PDF). The identity structure MUST be read from there. This is actually a perfect illustration of WHY ii′-0 exists: the exact identity bookkeeping needs to be pinned from the primary source before any Lean. Let me read Lagarias Thm 6.1 from the archived PDF (pages around §6).
The exact identity bookkeeping is subtle — I need Lagarias's Thm 6.1 from the primary source before designing the skeleton. Let me read it from the archived PDF.
tool: Bash
input
{
"command": "cd /d/Python/Quantum/research/paper && python -c \"\nimport subprocess, sys\ntry:\n from pypdf import PdfReader\nexcept ImportError:\n from PyPDF2 import PdfReader\nr = PdfReader('lagarias-li-coefficients-automorphic-0404394.pdf')\nprint('pages:', len(r.pages))\nfor i, p in enumerate(r.pages):\n t = p.extract_text() or ''\n if 'Theorem 6.1' in t or '6.1' in t[:200]:\n print(i, t[:100].replace(chr(10),' '))\n\" 2>&1 | head -20",
"description": "Locate Theorem 6.1 in Lagarias PDF"
}
raw content_json (database row)
[{"type":"thinking","thinking":"This is the full BL Thm 2. The decomposition λ_n = [prime part] + [archimedean] + [trivial-zeros part]. In the BL/S_f notation:\n\nλ_n = S_∞(n) − S_f(n) + 1? Let me map: from the formula:\nλ_n = −[prime part with signs] + 1 − (log 4π+γ)n/2 − Σ_{j=2}^n (−1)^{j−1}C(n,j)(1−2^{−j})ζ(j)\n\nHmm, and D-037's S_∞: S_∞(n) = Σ_{j=1}^n (−1)^j C(n,j) a_j with a₁ = (log 4π+γ)/2, a_j = (1−2^{−j})ζ(j). Let me expand: Σ_{j=1}^n (−1)^j C(n,j) a_j = −n·a₁ + Σ_{j=2}^n (−1)^j C(n,j)(1−2^{−j})ζ(j) = −n(log4π+γ)/2 + Σ_{j=2}(−1)^j C(n,j)(1−2^{−j})ζ(j).\n\nBL's λ_n has: 1 − (log4π+γ)n/2 − Σ_{j=2}^n (−1)^{j−1}C(n,j)(1−2^{−j})ζ(j) + [prime part] = 1 + S_∞(n) + [prime part with the −(−1)^{j−1} = (−1)^j sign flip ✓].\n\nSo: λ_n = S_∞(n) + 1 + [prime part], i.e. **S_f(n) = −[prime part]** where [prime part] = −Σ_{j=1}^n (−1)^{j−1}C(n,j)η_{j-1} = Σ_j (−1)^j C(n,j)η_{j−1}... with η the regularized prime sums (Stieltjes-analogues). So:\n\n**S_f(n) = Σ_{j=1}^n (−1)^{j−1}C(n,j) η_{j−1}** (up to sign convention) — the \"finite part\" is the PRIME-side sum, NOT a zero sum!! That's the whole point of the BL decomposition: S_∞ = archimedean, S_f = finite primes. And |S_f(n)| = 13.52 max measured (E-021, arithmetic side) ✓ consistent with Maślanka.\n\nAnd how does Lagarias's Thm 6.1 bound S_f? Via the contour: the ZERO-side representation of S_f: since λ_n = Σ_ρ[1−(1−1/ρ)ⁿ] (zero side) and λ_n = S_∞ − S_f + 1 (arithmetic side), we get S_f = S_∞ + 1 − Σ_ρ[1−(1−1/ρ)ⁿ]. The contour argument computes the truncated zero sum:\nλ_n(T) := Σ_{|γ|≤T}[1−(1−1/ρ)ⁿ]... and Thm 6.1 says S_f(n) = λ_n(√n) + O(√n log n)?? Hmm wait — that means S_∞ + 1 − λ_n = λ_n(√n) + O — i.e. S_∞ + 1 − [λ_n − λ_n(√n)] = O... i.e. the identity is:\n\nS_f(n) − λ_n(√n) = S_∞(n) + 1 − λ_n + λ_n(√n) = S_∞(n) + 1 − [tail beyond √n].\n\nAnd Lagarias's contour bounds THIS: S_f(n) − λ_n(T) for T ~ √n. Let me recompute with my residue bookkeeping:\n\nS_f − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T)... no wait: S_f = S_∞ + 1 − λ_n so S_f − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T). Hmm that's not what I want either. Let me recompute from D-032's statement: \"S_f(n) = λ_n(√n) + O(√n log n)\". Since |S_f| is small (13.52) and |λ_n(√n)| ≤ 2N(√n) ~ √n log n — the theorem S_f = λ_n(√n) + O(√n log n) is consistent with both being O(√n log n). \n\nSo the contour's job: bound S_f(n) − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T)... hmm, no: = S_∞ + 1 − (λ_n − λ_n(T))·... S_f − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T). With λ_n = λ_n(T) + tail: S_f − λ_n(T) = S_∞ + 1 − 2λ_n(T) − tail?? That doesn't look right either. Hmm, unless λ_n(T) in Lagarias's notation means the truncated sum of the OTHER form: λ_n(T) = Σ_{|γ|≤T}[1−(ρ/(ρ−1))ⁿ]?\n\nLet me test at n=1, T=2: λ_1(2) = 0 (no zeros below 2). S_f(1) = −0.577. So \"S_f = λ_n(√n) + O\" would read −0.577 = 0 + O — the O-term must BE −0.577 at n=1. The O-terms at n=1: sides + triv + Res₁ + ... : from my earlier bookkeeping: S_f = bracket_A − sides − tail with bracket_A = S_∞+1+Res₁+triv. S_f(1) = −0.577, bracket_A(1) = −0.798, sides+tail = −0.2208 (computed). OK whatever — the identity is exact but opaque.\n\nSTOP. The numerology is settled enough for the DESIGN: the exact identity structure is:\n\n(ID) Sf(n) = λ_n(T) + [bracket_A(n) − λ_n(T) − sides(T) − tail(T)]... \n\nNo. Let me just carefully redo the algebra ONE more time, cleanly, with all signs, because this determines the skeleton.\n\nDefinitions:\n- λ_n := full paired zero sum Σ_ρ[1−(1−1/ρ)ⁿ] (this is liCoeffP).\n- Sf := S_∞ + 1 − λ_n (Lean def, matches BL's λ_n = S_∞ − S_f + 1).\n- λ_n(T) := Σ_{|γ|≤T}′ mult·[1−(1−1/ρ)ⁿ] (truncated).\n- tail(T) := λ_n − λ_n(T) (paired tail).\n- Contour: integrand J(w) = k_n(w−1)·(−ζ′/ζ(w)), rectangle σ ∈ [−5/2, 2√n], |Im| ≤ T, T ∈ [√n, √n+1] zero-avoiding.\n- Residue theorem: sides(T) := (1/2πi)∮ J = Σ_{|γ|<T} mult·Res_ρ + Res_{−2} + Res₁(n).\n - Res_ρ J = k_n(ρ−1)·(−mult... for a zero of order m: −ζ′/ζ has residue −m). So Res = −m·k_n(ρ−1) = −m[(ρ/(ρ−1))ⁿ − 1] = m[1 − (ρ/(ρ−1))ⁿ].\n - So Σ_{|γ|<T} mult·[1−(ρ/(ρ−1))ⁿ] = sides − Res_{−2} − Res₁. ... (∗)\n- KEY REINDEXING: for the TRUNCATED sum (finite!): Σ_{|γ|<T} mult[1−(ρ/(ρ−1))ⁿ] = Σ_{|γ|<T} mult[1−(1−1/ρ)ⁿ] = λ_n(T) (the finite sum reindexed termwise via ρ ↦ 1−ρ, which preserves the truncated set and mult). EXACT, no convergence issue. ✓\n\nSo (∗) becomes: **λ_n(T) = sides(T) − Res_{−2} − Res₁(n)**. \n\nThen: Sf = S_∞ + 1 − λ_n = S_∞ + 1 − λ_n(T) − tail(T) = S_∞ + 1 − [sides − Res_{−2} − Res₁] − tail.\n\n|Sf| ≤ |λ_n(T)| + |S_∞ + 1 + Res_{−2} + Res₁ − sides − tail|?? No — the clean split for the BOUND:\n\n|Sf| ≤ |λ_n(T)| + |Sf − λ_n(T)|, and Sf − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T) = S_∞ + 1 − 2λ_n(T) − tail — ugh, that has 2λ_n(T) ~ 2·2N(T) — no good.\n\nThe RIGHT split (Lagarias's): Sf = S_∞ + 1 − λ_n = [S_∞ + 1 + Res₁ + Res_{−2} − sides] + [sides − Res_{−2} − Res₁] − tail... circular. OK let me just substitute once:\n\nSf = S_∞ + 1 + Res₁ + Res_{−2} − sides − tail.\n\nNow |Sf| ≤ |sides| + |tail| + |S_∞ + 1 + Res₁ + Res_{−2}|.\n\nFor B(n) to work: need |S_∞ + 1 + Res₁ + Res_{−2}| bounded by... 2N(T)?? At n=1: |−1.554120+1+(−0.577216)+(1/3)| = 0.798. 2N(2) = 0. Fails at n=1 but the range is n ≥ 10⁵. Hmm, but structurally \"2N(T) ~ √n log n\" vs \"S_∞+1+Res₁+triv ~ ?\". What IS the size of bracket_A(n) = S_∞+1+Res₁+triv at large n? If S_∞ + Res₁ cancels the ½n log n (which it must, since...), bracket_A could be O(n)-ish or smaller. Numerically: n=1: −0.798, n=2: −1.213. Let me estimate n=3 quickly? Tedious. The point for the DESIGN: the identity \n\n(★★) S_∞(n) + 1 + Res₁(n) + Res_{−2}(n) = 0 for all n \n\nwould give the clean form |Sf| ≤ |sides| + |tail|, but it's numerically FALSE (−0.798 at n=1). So that's not the identity. INSTEAD, the correct identity must be:\n\nS_∞(n) + 1 + Res₁(n) + Res_{−2}(n) = −(something of size ~λ_n(T))?? At n=1: −0.798 = ? λ_1(2)=0. Hmm. Or maybe my Res₁ is wrong by a sign or the trivial residue is different. Let me recompute Res_{−2}: integrand J(w) = k_n(w−1)(−ζ′/ζ(w)) at w=−2: −ζ′/ζ has residue −1 (simple zero), k_n(−3) = (1−1/3)ⁿ−1 = (2/3)ⁿ−1. Res_{−2} = k_n(−3)·(−1) = 1−(2/3)ⁿ ✓ (matches triv). At n=1: 1/3 ✓.\n\nHmm what if the w=1 pole's order and residue: k_n(w−1) = (1+1/(w−1))ⁿ − 1 = (w/(w−1))ⁿ − 1. At w=1: pole order n ✓. −ζ′/ζ ~ +1/(w−1) — WAIT: ζ has a simple pole at w=1 with residue +1, so ζ(w) ~ 1/(w−1), ζ′(w) ~ −1/(w−1)², ζ′/ζ ~ −1/(w−1). So −ζ′/ζ ~ +1/(w−1) ✓ c_{−1} = +1 ✓.\n\nSo Res₁(n) = Σ_{j=1}^n C(n,j)c_{j−1}, c_{−1}=1, c_0=−γ, c_1 = γ₁+γ². At n=1: c_0 = −γ ✓.\n\nLet me test (★★) variant: maybe S_∞ + Res₁ + Res_{−2} = −1 + λ_n-related? Compute S_f(1) directly via prime side: S_f(1) = η_0·(−1)^{?}: S_f(1) = Σ_{j=1}^1 (−1)^{j−1}C(1,j)η_{j−1} = η_0. And −ζ′/ζ(s+1) = 1/s + Σ η_n s^n (from Theorem-2.txt line 59): so η_0 = c_0 evaluated via s: −ζ′/ζ(s+1) at s=0: 1/s + η_0 + ...; and c_0 (constant term of −ζ′/ζ at s=1) = η_0 ✓!! So S_f(1) = η_0 = c_0 = Res₁(1) ✓✓ (consistent with my earlier finding S_f(1) = Res₁(1) = −γ).\n\nn=2: S_f(2) = Σ_{j=1}^2 (−1)^{j−1}C(2,j)η_{j−1} = 2η_0 − η_1 = −2γ − η_1. Res₁(2) = 2c_0 + c_1 = −2γ + (γ₁+γ²). And η_1 vs c_1: from −ζ′/ζ(s+1) = 1/s + η_0 + η_1 s + ...: with u = s: −ζ′/ζ(1+u) = u⁻¹ + c_0 + c_1 u + ...: so η_1 = c_1 = γ₁+γ². So S_f(2) = −2γ − (γ₁+γ²) = −1.154432 − 0.260361 = −1.414793?? But earlier I computed S_f(2) = S_∞(2)+1−λ_2 = −0.966922. CONTRADICTION! So my prime-side formula sign is wrong. From Theorem-2.txt: λ_n = −Σ_j (−1)^{j−1}C(n,j)η_{j−1} + 1 − (log4π+γ)n/2 − Σ_{j=2}(−1)^{j−1}C(n,j)(1−2^{−j})ζ(j).\n\nn=2: λ_2 = −[2η_0 − η_1] + 1 − (log4π+γ) − [−(−1)¹C(2,2)(1−1/4)ζ(2)] = −2η_0 + η_1 + 1 − log4π − γ + (3/4)ζ(2)·(−1)^{...}: careful: −Σ_{j=2}^2 (−1)^{j−1}C(n,j)a_j = −(−1)¹·1·a₂ = +a₂ = (3/4)ζ(2) = 1.233701. So λ_2 = −2(−0.577216) + 0.260361 + 1 − 3.108240 + 1.233701 = 1.154432 + 0.260361 + 1 − 3.108240 + 1.233701 = 0.540254?? But λ_2 = 0.092383. WRONG — so η_1 ≠ c_1. Indeed η_1 is the PRIME-side regularized sum (with Λ(m)), NOT the Stieltjes polynomial! Theorem-2.txt line 53: η_n defined via Λ(m) sums (prime side); line 61: each η_n is a polynomial in Stieltjes constants (γ₀s + (γ₁−½γ₀²)s²...) — those are the coefficients of log(sζ(s+1))'s expansion — a DIFFERENT relation. The expansion at line 59: −ζ′/ζ(s+1) = 1/s + Σ η_n s^n. Hmm — but −ζ′/ζ(s+1) near s=0: 1/s + c_0 + c_1 s... where c_k are the Laurent coeffs of −ζ′/ζ at 1. So η_n = c_n EXACTLY?? But then λ_2 computation gives 0.540 ≠ 0.0924. So one of my c's is wrong. Let me recompute c_1: −ζ′/ζ(1+u) = u⁻¹ − γ + (γ₁+γ²)u + ...? From log ζ(1+u) = −log u + γu − (γ₁+γ²)u²/2 + O(u³) — let me redo: ζ(1+u) = u⁻¹ + Σ_{k≥0} (−1)^k γ_k/k! u^k = u⁻¹(1 + γu − γ₁u²/2 + γ₂u³/6 − ...). log ζ = −log u + log(1 + v) where v = γu − γ₁u²/2 + ...: log(1+v) = v − v²/2 + ... = γu − γ₁u²/2 − γ²u²/2 + O(u³) = γu − (γ₁+γ²)u²/2 + O(u³). d/du: ζ′/ζ(1+u) = −u⁻¹ + γ − (γ₁+γ²)u + O(u²). −ζ′/ζ(1+u) = u⁻¹ − γ + (γ₁+γ²)u + O(u²). So c_0 = −γ, c_1 = γ₁ + γ² = −0.072816 + 0.333177 = 0.260361. Then λ_2 from the formula: −[2c_0 − c_1] + 1 − (log4π+γ) + a₂ = −2(−0.577216) + 0.260361 + 1 − 3.108240 + 1.233701 = 0.540254 ≠ 0.092383. So EITHER the formula's η's aren't the c's, or the archimedean sign conventions differ. Theorem-2.txt line 57-61 says η_n are the Λ-sum regularized limits AND equal the Laurent coefficients of −ζ′/ζ(s+1). Both can be true (PNT connection). Then λ_2 = 0.540?? But the true λ_2 = 0.0924. So the formula as transcribed in Theorem-2.txt must have a sign convention I'm misreading (e.g. the archimedean term is 1 − (log 4π + γ)n/2 with a DIFFERENT sign structure, or the a_j signs). At n=1: λ_1 = −η_0 + 1 − (log4π+γ)/2 = 0.577216 + 1 − 1.554120 = 0.023096 ✓✓ CORRECT. n=2: −[2η_0 − η_1] = −2(−0.577216) + η_1·(+1)... wait: −Σ_{j=1}^2 (−1)^{j−1}C(2,j)η_{j−1} = −[(−1)⁰C(2,1)η_0 + (−1)¹C(2,2)η_1] = −[2η_0 − η_1] = −2η_0 + η_1 = 1.154432 + 0.260361 = 1.414793. Archimedean+trivial: 1 − (log4π+γ)·2/2 − (−1)¹C(2,2)a₂ = 1 − 3.108240 + a₂ = 1 − 3.108240 + 1.233701 = −0.874539. λ_2 = 1.414793 − 0.874539 = 0.540254. True λ_2 = 0.092383. So η_1 must be ≈ 0.260361 − (0.540254−0.092383) = −0.187510. So the true η_1 ≈ −0.1875, NOT c_1 = 0.2604. So Theorem-2.txt's \"η_n = Laurent coeffs of −ζ′/ζ(s+1)\" is imprecise/wrong in sign or indexing — indeed line 61 gives the polynomial structure: γ₀s + (γ₁−½γ₀²)s²... suggesting the η's relate to log(sζ(s+1)) = log(1 + γs − ...) whose expansion gives γ₀ s + (−γ₁/2 − γ₀²/2)s²... hmm: log(sζ(s+1)): sζ(s+1) = 1 + γs − γ₁s²/2 + ...: log = γs − (γ₁+γ²)s²/2 + ... So the coefficients of log(sζ(s+1)) are γ, −(γ₁+γ²)/2, ... and η_1 (true) ≈ −0.1875. Hmm: −(γ₁+γ²)/2 = −0.130181. Not −0.1875 either. What IS η_1? True value: η_1 = lim{Σ_{m≤x} Λ(m)(log m)/m − (log x)²/2} ≈ −0.1875? Hmm I recall η_0 = −γ ≈ −0.5772? Hmm no — actually the standard values: η_0 = −0.5772 (= −γ ✓), η_1 = −0.18754...? Let me check: yes I believe the BL η constants are η_0 = −γ, η_1 ≈ −0.1875, η_2 ≈ −0.21..., — wait those are related to but not equal to Stieltjes. The relation: η_n = (−1)^n[...Stieltjes polynomial...]: η_1 = −(γ₁ + γ²/2)? = −(−0.072816 + 0.166589) = −0.093773. ×2 = −0.187546 ✓!! So η_1 = −2(γ₁ + γ²/2)... hmm so η_1 = −2γ₁ − γ² = 0.145631 − 0.333177 = −0.187546 ✓✓. So the relation is η_1 = −(γ² + 2γ₁), consistent with the line-61 polynomial pattern \"γ₀s + (γ₁ − ½γ₀²)s²\" up to signs/conventions I'm not going to nail by hand.\n\nCONCLUSION for the design: the exact identity content is:\n- **BL Thm 2 (P1)**: λ_n = S_∞(n) − S_f(n) + 1 where S_f(n) = Σ_j (−1)^{j−1}C(n,j)η_{j−1}·(−1)... (prime side) — with the η's being the regularized Λ-sums.\n- The CONTOUR identity (∗): λ_n(T) = sides(T) − Res_{−2} − Res₁(n), where Res₁(n) = Σ_j C(n,j)c_{j−1} with c's = Laurent coefficients of −ζ′/ζ at 1.\n- The BRIDGE between them: S_∞ + 1 − λ_n = S_f ⟹ S_f(n) = S_∞+1+Res₁+Res_{−2} − sides − tail — and the identification needed: **S_∞(n) + 1 + Res₁(n) + Res_{−2}(n) = S_f(n) + [something]** — NO wait. We have TWO expressions for S_f: prime side (definitional) and the contour one. The assembly for the BOUND: |S_f(n)| ≤ |λ_n(T)| + |S_f(n) − λ_n(T)|. And S_f(n) − λ_n(T) = S_∞ + 1 − λ_n − λ_n(T) — hmm this has λ_n (full) which we don't control...\n\nHmm hmm. OK here's the cleaner classical assembly (this must be what Lagarias does): DON'T split via |Sf| ≤ |λ_n(T)| + |Sf − λ_n(T)|. Instead substitute (∗) into Sf's definition:\n\nSf = S_∞ + 1 − λ_n = S_∞ + 1 − λ_n(T) − tail(T)\n = S_∞ + 1 − [sides − Res_{−2} − Res₁] − tail\n = [S_∞ + 1 + Res₁ + Res_{−2}] − sides − tail.\n\nNow the bracket: S_∞(n) + 1 + Res₁(n) + Res_{−2}(n). Claim: this bracket is EXACTLY 0 for all n ≥ 1?? Numerically: n=1: −1.554120 + 1 − 0.577216 + 0.333333 = −0.798003 ≠ 0. FALSE. Damn.\n\nAlternative: the bracket equals −2·[trivial-zero sum]?? The trivial zeros beyond −2 (−4, −6, ...) — the rectangle only encloses −2, but the FULL trivial-zero contribution to the archimedean side is the a_j = (1−2^{−j})ζ(j) sum (per Theorem-2.txt's anatomy: \"the (1−2^{−j})ζ(j) terms account for the trivial zeros\"). S_∞ contains Σ_j(−1)^j C(n,j)a_j — the trivial zeros are IN S_∞'s definition. So the bracket S_∞+1+Res₁+Res_{−2} mixes: S_∞ (archimedean + trivial in BL form), Res₁ (the contour's w=1 output), Res_{−2} (the one enclosed trivial zero). For the bracket to be manageable, the identity must be:\n\nS_∞(n) + Res₁(n) + 1 = −Res_{−2}(n) + bracket₀(n) with bracket₀ small?\n\nNumerically: S_∞+Res₁+1 at n=1: −1.554120 − 0.577216 + 1 = −1.131336. −Res_{−2} = −1/3. bracket₀ = −1.131336 + 0.333333 = −0.798003. Still −0.798. What's −0.798 at n=1? Hmm: −γ + ...? −0.798 = −γ − 0.2208 = c_0 − 0.2208. And 0.2208 ≈ ? log 2 − 0.47...? (2/3)? Meh.\n\nYou know what — I've spent enough on hand-numerology. The DEFINITIVE source is Lagarias's paper §6 (in our archive as PDF). The identity structure MUST be read from there. This is actually a perfect illustration of WHY ii′-0 exists: the exact identity bookkeeping needs to be pinned from the primary source before any Lean. Let me read Lagarias Thm 6.1 from the archived PDF (pages around §6).","signature":""},{"type":"text","text":"The exact identity bookkeeping is subtle — I need Lagarias's Thm 6.1 from the primary source before designing the skeleton. Let me read it from the archived PDF."},{"type":"tool_use","id":"Bash_43","name":"Bash","input":{"command":"cd /d/Python/Quantum/research/paper && python -c \"\nimport subprocess, sys\ntry:\n from pypdf import PdfReader\nexcept ImportError:\n from PyPDF2 import PdfReader\nr = PdfReader('lagarias-li-coefficients-automorphic-0404394.pdf')\nprint('pages:', len(r.pages))\nfor i, p in enumerate(r.pages):\n t = p.extract_text() or ''\n if 'Theorem 6.1' in t or '6.1' in t[:200]:\n print(i, t[:100].replace(chr(10),' '))\n\" 2>&1 | head -20","description":"Locate Theorem 6.1 in Lagarias PDF"},"extra_content":null}]