Contents

Irrationality of the dilogarithm at −1/2-1/2

Ruiyi Zhang Independent Researcher, Columbia, Maryland, USA
Abstract.

We prove that Li2⁡(−1/2)\operatorname{Li}_{2}(-1/2) is irrational. Irrationality of Li2⁡(1/s)\operatorname{Li}_{2}(1/s) was known for integers s≤−5s\leq-5 and s≥6s\geq 6 by work of Hata and of Rhin and Viola, and Prellberg recently proved it for s=−4,5,−3s=-4,5,-3. We define polynomials Gm⁢(T)G_{m}(T) with rational coefficients by an interpolation determinant. The rank of this determinant modulo primes pp with m/5<p<5⁢m/4m/5<p<5m/4 bounds the denominators of the coefficients. A concavity inequality for a logarithmic energy bounds Gm⁢(Li2⁡(−1/2))G_{m}(\operatorname{Li}_{2}(-1/2)). The only interval-arithmetic step is the potential certificate. The exact data in Table A.1, the coefficients of G4G_{4}, and the 9×99\times 9 determinant modulo 1111 are printed and hand-checkable. A congruence modulo 1111 shows that GmG_{m} has no rational zero for m=4⋅11tm=4\cdot 11^{t}.

Key words and phrases:
dilogarithm, irrationality, interpolation determinant, vector equilibrium, Padé approximation
2020 Mathematics Subject Classification:
Primary 11J72; Secondary 33B30

1. Introduction

For |z|≤1|z|\leq 1 let Li2⁡(z)=∑k≥1zk/k2\operatorname{Li}_{2}(z)=\sum_{k\geq 1}z^{k}/k^{2}. Hata [3, Theorem 1.1] proved that Li2⁡(1/s)\operatorname{Li}_{2}(1/s) is irrational for every integer s≤−5s\leq-5 and every s≥7s\geq 7. Rhin and Viola settled s=6s=6 by their permutation group method [7, Section 6], and later treated negative arguments [8, Theorem 2.1 and Section 4]. Viola and Zudilin [10] introduced a six-parameter family of double integrals for dilogarithmic values. In 2024 Calegari, Dimitrov and Tang [1, Remark 2.8.2 and Section 15.2] listed the integers still open: s∈{−4,−3,−2,2,3,4,5}s\in\{-4,-3,-2,2,3,4,5\}. Prellberg’s 2026 preprint [6, Theorem 1.1] proves the cases s=−4,5,−3s=-4,5,-3. We prove the case s=−2s=-2.

Theorem 1.

The number Li2⁡(−1/2)\operatorname{Li}_{2}(-1/2) is irrational.

By Landen’s identity [2, (25.12.3)], Li2⁡(1/3)+Li2⁡(−1/2)=−12⁢log2⁡(2/3)\operatorname{Li}_{2}(1/3)+\operatorname{Li}_{2}(-1/2)=-\frac{1}{2}\log^{2}(2/3), so the case s=3s=3 does not follow from Theorem 1.

Put β=Li2⁡(−1/2)\beta=\operatorname{Li}_{2}(-1/2) and m=4⁢nm=4n. Section 2 defines Gm⁢(T)∈ℚ⁢[T]G_{m}(T)\in\mathbb{Q}[T] by an interpolation determinant. Sections 3 and 4 give, for all large mm,

deg⁡Gm≤m⁢by Lemma 2.1,Dm∗⁢Gm⁢(T)∈ℤ⁢[T],Dm∗∈ℚ>0,\displaystyle\deg G_{m}\leq m\ \text{by Lemma~{}\ref{lem:moment-identity}},% \qquad D_{m}^{*}G_{m}(T)\in\mathbb{Z}[T],\qquad D_{m}^{*}\in\mathbb{Q}_{>0},
(1.1) lim supm→∞4∣mlog⁡Dm∗m2≤τ<4.33,log⁡|Gm⁢(β)|≤−4.36⁢m2+870⁢m⁢log⁡m,\displaystyle\limsup_{\begin{subarray}{c}m\to\infty\\ 4\mid m\end{subarray}}\frac{\log D_{m}^{*}}{m^{2}}\leq\tau<4.33,\qquad\log|G_{% m}(\beta)|\leq-4.36m^{2}+870m\log m,
τ=1881503369600−5348⁢log⁡2,4.36−τ>0.0347.\displaystyle\tau=\frac{1881503}{369600}-\frac{53}{48}\log 2,\qquad 4.36-\tau>% 0.0347.

If β=u/q0\beta=u/q_{0} with q0>0q_{0}>0, then q0m⁢Dm∗⁢Gm⁢(β)∈ℤq_{0}^{m}D_{m}^{*}G_{m}(\beta)\in\mathbb{Z}. Section 5 shows that it is nonzero for m=4⋅11tm=4\cdot 11^{t}, t≥1t\geq 1, and (1.1) shows that it tends to zero.

Lemma 4.2 is computer-assisted: it uses interval arithmetic (Appendix F). The other finite data are printed in Tables 2.1, 5.1 and A.1 and in (5.2). They were also rechecked in exact rational arithmetic.

2. The determinant

Put a=−2a=-2. For n≥1n\geq 1 let

m=4⁢n,σ=k=h=n,L=m+h,I={σ+k,…,σ+k+m−1},J={σ,…,σ+2⁢m+h−1},Λ={1,…,h},\begin{gathered}m=4n,\qquad\sigma=k=h=n,\qquad L=m+h,\\ I=\{\sigma+k,\ldots,\sigma+k+m-1\},\quad J=\{\sigma,\ldots,\sigma+2m+h-1\},% \quad\Lambda=\{1,\ldots,h\},\end{gathered}

all ordered increasingly. For j≥0j\geq 0 put

tj=∑u=1jaj−uu,bj=∑u=1jaj−uu2;t0=b0=0.t_{j}=\sum_{u=1}^{j}\frac{a^{j-u}}{u},\qquad b_{j}=\sum_{u=1}^{j}\frac{a^{j-u}% }{u^{2}};\qquad t_{0}=b_{0}=0.

For i∈Ii\in I, j∈Jj\in J and ℓ∈Λ\ell\in\Lambda let

Ci⁢j={(ai−aj)/(j−i),i≠j,0,i=j,Bi⁢j={(ti−tj)/(j−i),i≠j,bi,i=j,Ai⁢j=ai⁢δi⁢j,Vℓ⁢j=1j+ℓ,Mi⁢j=∫01∫01xi⁢yj2+x⁢y⁢𝑑x⁢𝑑y,\begin{gathered}C_{ij}=\begin{cases}(a^{i}-a^{j})/(j-i),&i\neq j,\\ 0,&i=j,\end{cases}\quad B_{ij}=\begin{cases}(t_{i}-t_{j})/(j-i),&i\neq j,\\ b_{i},&i=j,\end{cases}\\ A_{ij}=a^{i}\delta_{ij},\qquad V_{\ell j}=\frac{1}{j+\ell},\qquad M_{ij}=\int_% {0}^{1}\!\int_{0}^{1}\frac{x^{i}y^{j}}{2+xy}\,dx\,dy,\end{gathered}

and put

(2.1) Gm⁢(T)=det(CB−T⁢AV),G_{m}(T)=\det\begin{pmatrix}C\\ B-TA\\ V\end{pmatrix},

a determinant of size 2⁢m+h2m+h. For an indeterminate aa the same formulas define

Gm⁢(T;a)=det(C⁢(a)B⁢(a)−T⁢A⁢(a)V)∈ℚ⁢[a,T];Gm⁢(T)=Gm⁢(T;−2).G_{m}(T;a)=\det\begin{pmatrix}C(a)\\ B(a)-TA(a)\\ V\end{pmatrix}\in\mathbb{Q}[a,T];\qquad G_{m}(T)=G_{m}(T;-2).
Lemma 2.1.

One has Gm∈ℚ⁢[T]G_{m}\in\mathbb{Q}[T], deg⁡Gm≤m\deg G_{m}\leq m, and

(2.2) M=B+log⁡(3/2)⁢C−β⁢A,Gm⁢(β)=det(CMV).M=B+\log(3/2)C-\beta A,\qquad G_{m}(\beta)=\det\begin{pmatrix}C\\ M\\ V\end{pmatrix}.
Proof.

Write λ=log⁡(3/2)\lambda=\log(3/2) and

𝔪j=∫01uj2+u⁢𝑑u,𝔫j=∫01uj⁢log⁡(1/u)2+u⁢𝑑u.\mathfrak{m}_{j}=\int_{0}^{1}\frac{u^{j}}{2+u}\,du,\qquad\mathfrak{n}_{j}=\int% _{0}^{1}\frac{u^{j}\log(1/u)}{2+u}\,du.

Polynomial division and the geometric series give

𝔪0=λ,𝔫0=−β,𝔪j=a⁢𝔪j−1+1j,𝔫j=a⁢𝔫j−1+1j2(j≥1),𝔪j=aj⁢λ+tj,𝔫j=bj−aj⁢β.\begin{gathered}\mathfrak{m}_{0}=\lambda,\quad\mathfrak{n}_{0}=-\beta,\qquad% \mathfrak{m}_{j}=a\mathfrak{m}_{j-1}+\frac{1}{j},\quad\mathfrak{n}_{j}=a% \mathfrak{n}_{j-1}+\frac{1}{j^{2}}\quad(j\geq 1),\\ \mathfrak{m}_{j}=a^{j}\lambda+t_{j},\qquad\mathfrak{n}_{j}=b_{j}-a^{j}\beta.% \end{gathered}

The substitution u=x⁢yu=xy gives

Mi⁢j=𝔪i−𝔪jj−i(i≠j),Mi⁢i=𝔫i.M_{ij}=\frac{\mathfrak{m}_{i}-\mathfrak{m}_{j}}{j-i}\quad(i\neq j),\qquad M_{% ii}=\mathfrak{n}_{i}.

Hence the first identity in (2.2). Add λ\lambda times each CC row to the corresponding row of B−β⁢AB-\beta A. Only the mm middle rows depend on TT, and their coefficients are rational. ∎

Put r=3/5r=3/5 and

γ+⁢(θ)=r⁢ei⁢θ(0≤θ≤π),γ−⁢(θ)=r⁢ei⁢θ(π≤θ≤2⁢π).\gamma_{+}(\theta)=re^{i\theta}\quad(0\leq\theta\leq\pi),\qquad\gamma_{-}(% \theta)=re^{i\theta}\quad(\pi\leq\theta\leq 2\pi).

With the real intervals oriented increasingly, let

Γ=[−2,−r]+[r,1]+(cΓ−1)⁢γ++cΓ⁢γ−,cΓ=log⁡2+i⁢π2⁢π⁢i,\Gamma=[-2,-r]+[r,1]+(c_{\Gamma}-1)\gamma_{+}+c_{\Gamma}\gamma_{-},\qquad c_{% \Gamma}=\frac{\log 2+i\pi}{2\pi i},

a linear combination of paths. For every integer dd,

(2.3) ∫Γzd−1⁢𝑑z={(1−ad)/d,d≠0,0,d=0.\int_{\Gamma}z^{d-1}\,dz=\begin{cases}(1-a^{d})/d,&d\neq 0,\\ 0,&d=0.\end{cases}

For d≠0d\neq 0 integrate zd/dz^{d}/d on each path. For d=0d=0 the real intervals give −log⁡2-\log 2 and the arcs give i⁢π⁢(2⁢cΓ−1)=log⁡2i\pi(2c_{\Gamma}-1)=\log 2.

For u=(u1,…,ud)u=(u_{1},\ldots,u_{d}) put Δ⁢(u)=∏i<j(uj−ui)\Delta(u)=\prod_{i<j}(u_{j}-u_{i}). Use variables x=(x1,…,xm)x=(x_{1},\ldots,x_{m}), v=(v1,…,vL)v=(v_{1},\ldots,v_{L}) and z=(z1,…,zm)z=(z_{1},\ldots,z_{m}).

Lemma 2.2.

With product differentials d⁢x⁢d⁢v⁢d⁢zdx\,dv\,dz,

(2.4) Gm⁢(β)=23⁢m2/2−m(m!)2⁢L!⁢∫[0,1]m∫[0,1]L∫ΓmΦ⁢(x,v,z)⁢𝑑z⁢𝑑v⁢𝑑x,G_{m}(\beta)=\frac{2^{3m^{2}/2-m}}{(m!)^{2}L!}\int_{[0,1]^{m}}\int_{[0,1]^{L}}% \int_{\Gamma^{m}}\Phi(x,v,z)\,dz\,dv\,dx,

where

Φ⁢(x,v,z)=Δ⁢(x)2⁢Δ⁢(v)2⁢Δ⁢(z)2⁢∏d=1m∏b=1L(zd−vb)⁢∏i′=1mxi′σ+k+h⁢∏b=1Lvbσ⁢∏d=1mzd−m−k∏i′=1m∏b=1L(2+xi′⁢vb).\Phi(x,v,z)=\frac{\Delta(x)^{2}\Delta(v)^{2}\Delta(z)^{2}\displaystyle\prod_{d% =1}^{m}\prod_{b=1}^{L}(z_{d}-v_{b})\displaystyle\prod_{i^{\prime}=1}^{m}x_{i^{% \prime}}^{\sigma+k+h}\prod_{b=1}^{L}v_{b}^{\sigma}\prod_{d=1}^{m}z_{d}^{-m-k}}% {\displaystyle\prod_{i^{\prime}=1}^{m}\prod_{b=1}^{L}(2+x_{i^{\prime}}v_{b})}.
Proof.

By (2.3), Ci⁢j=ai⁢∫Γzj−i−1⁢𝑑zC_{ij}=a^{i}\int_{\Gamma}z^{j-i-1}\,dz. Put fi⁢(v)=∫01xi/(2+x⁢v)⁢𝑑xf_{i}(v)=\int_{0}^{1}x^{i}/(2+xv)\,dx. Then

Mi⁢j=∫01fi⁢(v)⁢vj⁢𝑑v,Vℓ⁢j=∫01vℓ−1⁢vj⁢𝑑v.M_{ij}=\int_{0}^{1}f_{i}(v)v^{j}\,dv,\qquad V_{\ell j}=\int_{0}^{1}v^{\ell-1}v% ^{j}\,dv.

Set eΔ=m⁢(m−1)/2e_{\Delta}=m(m-1)/2 and 𝐰=(z1,…,zm,v1,…,vL)\mathbf{w}=(z_{1},\ldots,z_{m},v_{1},\ldots,v_{L}). Multilinearity in the rows and antisymmetrization, applied to the LL rows [M;V][M;V], then to the mm rows CC, and then to the mm functions fif_{i}, give

Gm⁢(β)=1(m!)2⁢L!⁢∫[0,1]m×[0,1]L×ΓmDz⁢Dz⁢v⁢Dx⁢Dx⁢v⁢𝑑x⁢𝑑v⁢𝑑z,G_{m}(\beta)=\frac{1}{(m!)^{2}L!}\int_{[0,1]^{m}\times[0,1]^{L}\times\Gamma^{m% }}D_{z}D_{zv}D_{x}D_{xv}\,dx\,dv\,dz,

where

Dz\displaystyle D_{z} =det(ai⁢zd−i−1)i∈I, 1≤d≤m=(−1)eΔ⁢am⁢(σ+k)+eΔ⁢Δ⁢(z)⁢∏dzd−σ−k−m,\displaystyle=\det(a^{i}z_{d}^{-i-1})_{i\in I,\,1\leq d\leq m}=(-1)^{e_{\Delta% }}a^{m(\sigma+k)+e_{\Delta}}\Delta(z)\prod_{d}z_{d}^{-\sigma-k-m},
Dz⁢v\displaystyle D_{zv} =det(𝐰bj)1≤b≤2⁢m+h,j∈J=(−1)m⁢L⁢Δ⁢(z)⁢Δ⁢(v)⁢∏dzdσ⁢∏bvbσ⁢∏d,b(zd−vb),\displaystyle=\det(\mathbf{w}_{b}^{j})_{1\leq b\leq 2m+h,\,j\in J}=(-1)^{mL}% \Delta(z)\Delta(v)\prod_{d}z_{d}^{\sigma}\prod_{b}v_{b}^{\sigma}\prod_{d,b}(z_% {d}-v_{b}),
Dx\displaystyle D_{x} =det(xi′i)i∈I, 1≤i′≤m=Δ⁢(x)⁢∏i′xi′σ+k,\displaystyle=\det(x_{i^{\prime}}^{i})_{i\in I,\,1\leq i^{\prime}\leq m}=% \Delta(x)\prod_{i^{\prime}}x_{i^{\prime}}^{\sigma+k},
Dx⁢v\displaystyle D_{xv} =det(((2+xi′⁢vb)−1)1≤i′≤m, 1≤b≤L(vbℓ−1)1≤ℓ≤h, 1≤b≤L)=aeΔ⁢Δ⁢(x)⁢Δ⁢(v)⁢∏i′xi′h∏i′,b(2+xi′⁢vb).\displaystyle=\det\begin{pmatrix}((2+x_{i^{\prime}}v_{b})^{-1})_{1\leq i^{% \prime}\leq m,\,1\leq b\leq L}\\ (v_{b}^{\ell-1})_{1\leq\ell\leq h,\,1\leq b\leq L}\end{pmatrix}=\frac{a^{e_{% \Delta}}\Delta(x)\Delta(v)\prod_{i^{\prime}}x_{i^{\prime}}^{h}}{\prod_{i^{% \prime},b}(2+x_{i^{\prime}}v_{b})}.

The last equality follows from the Cauchy–Vandermonde identity. Explicitly, for ti′=2/xi′t_{i^{\prime}}=2/x_{i^{\prime}},

det(((2+xi′⁢vb)−1)i′,b(vbℓ−1)ℓ,b)=∏i′xi′−1⁢det(((ti′+vb)−1)i′,b(vbℓ−1)ℓ,b).\det\begin{pmatrix}((2+x_{i^{\prime}}v_{b})^{-1})_{i^{\prime},b}\\ (v_{b}^{\ell-1})_{\ell,b}\end{pmatrix}=\prod_{i^{\prime}}x_{i^{\prime}}^{-1}% \det\begin{pmatrix}((t_{i^{\prime}}+v_{b})^{-1})_{i^{\prime},b}\\ (v_{b}^{\ell-1})_{\ell,b}\end{pmatrix}.

The Cauchy–Vandermonde identity is

det(((ti′+vb)−1)i′,b(vbℓ−1)ℓ,b)=Δ⁢(t)⁢Δ⁢(v)∏i′,b(ti′+vb).\det\begin{pmatrix}((t_{i^{\prime}}+v_{b})^{-1})_{i^{\prime},b}\\ (v_{b}^{\ell-1})_{\ell,b}\end{pmatrix}=\frac{\Delta(t)\Delta(v)}{\prod_{i^{% \prime},b}(t_{i^{\prime}}+v_{b})}.

This is the Cauchy–Vandermonde identity [4]. With the preceding row-scaling identity it gives the displayed formula for Dx⁢vD_{xv}, and continuity covers xi′=0x_{i^{\prime}}=0. We verify the identity by induction on mm. At the induction step fix i′i^{\prime}, let v′=(v1,…,vL−1)v^{\prime}=(v_{1},\ldots,v_{L-1}), and let t^i′\widehat{t}_{i^{\prime}} denote tt with ti′t_{i^{\prime}} omitted. After multiplying by vL+ti′v_{L}+t_{i^{\prime}}, the residue of the quotient at vL=−ti′v_{L}=-t_{i^{\prime}} is

(−1)i′+L⁢Δ⁢(t^i′)⁢Δ⁢(v′)∏j≠i′∏b<L(tj+vb).(-1)^{i^{\prime}+L}\frac{\Delta(\widehat{t}_{i^{\prime}})\Delta(v^{\prime})}{% \prod_{j\neq i^{\prime}}\prod_{b<L}(t_{j}+v_{b})}.

By induction this is the signed complementary minor on the determinant side, so the residues agree. After subtracting them, the difference is a polynomial in vLv_{L} of degree at most h−1h-1 and vanishes at the other L−1≥hL-1\geq h variables v1,…,vL−1v_{1},\ldots,v_{L-1}. For h=0h=0 it tends to zero at infinity. The case m=0m=0 is the Vandermonde identity.

The product of the four factors is (−1)eΔ+m⁢L⁢am⁢(σ+k)+2⁢eΔ⁢Φ(-1)^{e_{\Delta}+mL}a^{m(\sigma+k)+2e_{\Delta}}\Phi. Since m=4⁢nm=4n, both signs are positive and m⁢(σ+k)+2⁢eΔ=3⁢m2/2−mm(\sigma+k)+2e_{\Delta}=3m^{2}/2-m. ∎

An exact calculation gives

(2.5) 16336404000000000⁢G4⁢(T)=∑j=04pj⁢Tj,16336404000000000\,G_{4}(T)=\sum_{j=0}^{4}p_{j}T^{j},

with the coefficients in Table 2.1.

Table 2.1. The integer coefficients in (2.5).
jj pjp_{j}
0 11105448366574821464202411110544836657482146420241
11 99027535307004682303995609902753530700468230399560
22 3311364976073005456321920033113649760730054563219200
33 4921251432479346585600000049212514324793465856000000
44 2742680795140797235200000027426807951407972352000000

3. The divisor

Put

Nden=10⁢n−1,U=11⁢n−1,dℓ=lcm⁡(1,…,ℓ),𝒟m=dNden3⁢m⁢dUn.N_{\mathrm{den}}=10n-1,\qquad U=11n-1,\qquad d_{\ell}=\operatorname{lcm}(1,% \ldots,\ell),\qquad\mathcal{D}_{m}=d_{N_{\mathrm{den}}}^{3m}d_{U}^{n}.

For F∈ℚ⁢[T]F\in\mathbb{Q}[T] let vp⁢(F)=minj⁡vp⁢([Tj]⁢F)v_{p}(F)=\min_{j}v_{p}([T^{j}]F).

Lemma 3.1.

For an indeterminate aa,

(3.1) Pm⁢(a,T):=𝒟m⁢Gm⁢(T;a)∈ℤ⁢[a,T].P_{m}(a,T):=\mathcal{D}_{m}G_{m}(T;a)\in\mathbb{Z}[a,T].
Proof.

The denominators of a CC entry divide dNdend_{N_{\mathrm{den}}}. An off-diagonal BB entry is a difference of harmonic sums divided by j−ij-i, so its denominator divides dNden2d_{N_{\mathrm{den}}}^{2}, and so does that of bib_{i}. The denominators of a VV entry divide dUd_{U}. ∎

3.1. Residue classes

For ρ>0\rho>0 and 0≤x<ρ0\leq x<\rho let

(3.2) K⁢(x,ρ)\displaystyle K(x,\rho) ={k≥0:1/2≤x+k⁢ρ<3/2},\displaystyle=\{k\geq 0:1/2\leq x+k\rho<3/2\},
ℒ⁢(x,ρ)\displaystyle\mathcal{L}(x,\rho) ={l≥0:1/4≤x+l⁢ρ<5/2},\displaystyle=\{l\geq 0:1/4\leq x+l\rho<5/2\},
W⁢(x,ρ)\displaystyle W(x,\rho) ={w≥1:0<w⁢ρ−x≤1/4},\displaystyle=\{w\geq 1:0<w\rho-x\leq 1/4\},

and

(3.3) 𝔯⁢(x,ρ)=min⁡(2⁢|K|+|W|,|ℒ|)−𝟏{K={0},ℒ={0,1},W=∅},𝔠⁢(ρ)=1+∫0ρ𝔯⁢(x,ρ)⁢𝑑x.\mathfrak{r}(x,\rho)=\min(2|K|+|W|,|\mathcal{L}|)-\mathbf{1}_{\{K=\{0\},\,% \mathcal{L}=\{0,1\},\,W=\varnothing\}},\qquad\mathfrak{c}(\rho)=1+\int_{0}^{% \rho}\mathfrak{r}(x,\rho)\,dx.
Lemma 3.2.

Let pp be an odd prime with p2>Up^{2}>U, and let a=−2a=-2. Then

(3.4) vp⁢(Gm)≥−m−∑u=0p−1𝔯⁢(u/m,p/m)=−m⁢𝔠⁢(p/m).v_{p}(G_{m})\geq-m-\sum_{u=0}^{p-1}\mathfrak{r}(u/m,p/m)=-m\mathfrak{c}(p/m).
Proof.

Multiply the C,B,VC,B,V rows by p,p2,pp,p^{2},p. Write i=u+k⁢pi=u+kp, j=u+l⁢pj=u+lp with 0≤u<p0\leq u<p. Only the harmonic terms with pp in the denominator survive modulo pp, and ap≡a(modp)a^{p}\equiv a\pmod{p}, so

p⁢tu+k⁢p≡au⁢tk,p2⁢bu+k⁢p≡au⁢bk(modp).pt_{u+kp}\equiv a^{u}t_{k},\qquad p^{2}b_{u+kp}\equiv a^{u}b_{k}\pmod{p}.

The scaled CC and BB entries vanish modulo pp between distinct residue classes modulo pp. For a cross-class BB entry, the numerator is pp-integral after the row scaling and the remaining factor p/(j−i)p/(j-i) is zero modulo pp. Within a class, for j≠ij\neq i, j−i=p⁢(l−k)j-i=p(l-k) with 0<|l−k|<p0<|l-k|<p, so the quotient denominator is a unit. Diagonal entries have their displayed separate values. Thus within a class they reduce to au⁢Ck⁢la^{u}C_{kl} and au⁢Bk⁢la^{u}B_{kl}, and −p2⁢T⁢A≡0-p^{2}TA\equiv 0. A VV row reduces to zero unless ℓ+u=p⁢w\ell+u=pw. Then its entry is 1/(l+w)1/(l+w), without the factor aua^{u}. All displayed denominators are units because p2>Up^{2}>U. So, after permutations, the reduction has blocks with 2⁢|K|+|W|2|K|+|W| rows and |ℒ||\mathcal{L}| columns.

For K={0}K=\{0\}, ℒ={0,1}\mathcal{L}=\{0,1\}, W=∅W=\varnothing the block, with its rows divided by the unit aua^{u}, is

(3.5) (01−a0−1)=(030−1),\begin{pmatrix}0&1-a\\ 0&-1\end{pmatrix}=\begin{pmatrix}0&3\\ 0&-1\end{pmatrix},

of rank one for every odd pp. Every other block has rank at most the minimum of its dimensions. So the rank is at most R=∑u𝔯⁢(u/m,p/m)R=\sum_{u}\mathfrak{r}(u/m,p/m). The matrix has 2⁢m+n2m+n rows. Row operations over ℤ(p)\mathbb{Z}_{(p)}, independent of TT, make its last 2⁢m+n−R2m+n-R rows divisible by pp. Subtracting the scaling exponent 3⁢m+n3m+n gives the first inequality in (3.4).

Every boundary in (3.2) has the form x=c0+b⁢ρx=c_{0}+b\rho with 4⁢c0,b∈ℤ4c_{0},b\in\mathbb{Z}. At ρ=p/m\rho=p/m, m=4⁢nm=4n, it lies in m−1⁢ℤm^{-1}\mathbb{Z}. The K,ℒK,\mathcal{L} conditions are closed on the left and open on the right in xx, and the WW condition is w⁢ρ−1/4≤x<w⁢ρw\rho-1/4\leq x<w\rho. So 𝔯\mathfrak{r} is constant on each [u/m,(u+1)/m)[u/m,(u+1)/m), and

(3.6) ∑u=0p−1𝔯⁢(u/m,p/m)=m⁢∫0p/m𝔯⁢(x,p/m)⁢𝑑x.∎\sum_{u=0}^{p-1}\mathfrak{r}(u/m,p/m)=m\int_{0}^{p/m}\mathfrak{r}(x,p/m)\,dx.\qed

For 1/5<ρ<5/41/5<\rho<5/4 only 0≤k≤70\leq k\leq 7, 0≤l≤120\leq l\leq 12 and 1≤w≤21\leq w\leq 2 occur. The sets in (3.2) change only on the lines

(3.7) x=0,x=ρ,x=1/2−k⁢ρ,x=3/2−k⁢ρ(0≤k≤7),x=1/4−l⁢ρ,x=5/2−l⁢ρ(0≤l≤12),x=w⁢ρ−1/4,x=w⁢ρ(1≤w≤2).\begin{gathered}x=0,\quad x=\rho,\quad x=1/2-k\rho,\quad x=3/2-k\rho\quad(0% \leq k\leq 7),\\ x=1/4-l\rho,\quad x=5/2-l\rho\quad(0\leq l\leq 12),\\ x=w\rho-1/4,\quad x=w\rho\quad(1\leq w\leq 2).\end{gathered}

Cut (1/5,5/4)(1/5,5/4) at the values of ρ\rho where two of these lines meet in 0≤x≤ρ0\leq x\leq\rho, and let 𝒞\mathcal{C} be the set of these cut points. On each resulting interval the order of the lines in 0<x<ρ0<x<\rho is fixed, so the sets (3.2) are constant on each xx-interval between consecutive lines, and 𝔠\mathfrak{c} is affine. At both endpoints of each row, its five printed lines ℓ1,…,ℓ5\ell_{1},\ldots,\ell_{5}, in their printed order, satisfy

0≤ℓ1⁢(ρ)≤ℓ2⁢(ρ)≤⋯≤ℓ5⁢(ρ)≤ρ.0\leq\ell_{1}(\rho)\leq\ell_{2}(\rho)\leq\cdots\leq\ell_{5}(\rho)\leq\rho.

Every other distinct line in (3.7) is ≤0\leq 0 at both endpoints or ≥ρ\geq\rho at both endpoints. Affineness and distinctness make the ordering strict in the interior, apart from boundary lines. Hence no intersection in 0≤x≤ρ0\leq x\leq\rho is omitted. Table A.1 lists the 4444 intervals, the five lines in 0<x<ρ0<x<\rho in increasing order, the values of 𝔯\mathfrak{r} on the six xx-intervals, 𝔠⁢(ρ)\mathfrak{c}(\rho) and ∫𝔠\int\mathfrak{c}. Each row can be checked by hand from (3.3). The table was also rechecked in exact rational arithmetic (Appendix F). Summation gives

(3.8) ∫1/55/4𝔠⁢(ρ)⁢𝑑ρ=1121513369600.\int_{1/5}^{5/4}\mathfrak{c}(\rho)\,d\rho=\frac{1121513}{369600}.

The configuration of (3.5) occupies 5/2−2⁢ρ≤x<ρ−1/45/2-2\rho\leq x<\rho-1/4 for 11/12<ρ<111/12<\rho<1, and 1/2≤x<ρ−1/41/2\leq x<\rho-1/4 for 1≤ρ<5/41\leq\rho<5/4. Its contribution to (3.8) is

−∫11/121(3⁢ρ−11/4)⁢𝑑ρ−∫15/4(ρ−3/4)⁢𝑑ρ=−548.-\int_{11/12}^{1}(3\rho-11/4)\,d\rho-\int_{1}^{5/4}(\rho-3/4)\,d\rho=-\frac{5}% {48}.

3.2. Larger primes

For p>5⁢np>5n define row potentials λC=0\lambda_{C}=0, λB⁢(i)\lambda_{B}(i), λV⁢(ℓ)\lambda_{V}(\ell) and column potentials χ⁢(j)\chi(j) by

5⁢n<p<6⁢n:\displaystyle 5n<p<6n: (λB,λV,χ)=(1+𝟏i≥p, 1, 1j≥p+2⁢n),\displaystyle\quad(\lambda_{B},\lambda_{V},\chi)=(1+\mathbf{1}_{i\geq p},\,1,% \,\mathbf{1}_{j\geq p+2n}),
6⁢n≤p<8⁢n:\displaystyle 6n\leq p<8n: (λB,λV,χ)=(𝟏i+p<10⁢n, 1, 1j≥p),\displaystyle\quad(\lambda_{B},\lambda_{V},\chi)=(\mathbf{1}_{i+p<10n},\,1,\,% \mathbf{1}_{j\geq p}),
8⁢n≤p<10⁢n:\displaystyle 8n\leq p<10n: (λB,λV,χ)=(0, 1, 1j≥p),\displaystyle\quad(\lambda_{B},\lambda_{V},\chi)=(0,\,1,\,\mathbf{1}_{j\geq p}),
10⁢n≤p≤11⁢n−1:\displaystyle 10n\leq p\leq 11n-1: (λB,λV,χ)=(0, 1p−ℓ∈J, 0),\displaystyle\quad(\lambda_{B},\lambda_{V},\chi)=(0,\,\mathbf{1}_{p-\ell\in J}% ,\,0),

and put

Epout={19⁢n−2⁢p,5⁢n<p<8⁢n,11⁢n−p,8⁢n≤p≤U.E_{p}^{\rm out}=\begin{cases}19n-2p,&5n<p<8n,\\ 11n-p,&8n\leq p\leq U.\end{cases}
Lemma 3.3.

For an odd prime 5⁢n<p≤U5n<p\leq U one has vp⁢(Gm)≥−Epoutv_{p}(G_{m})\geq-E_{p}^{\rm out}.

Proof.

A CC entry has a pole only at j=i+pj=i+p, of order at most one. A BB entry has pole order at most

𝟏max⁡(i,j)≥p+𝟏{max⁡(i,j)≥p,i≡j(modp)},\mathbf{1}_{\max(i,j)\geq p}+\mathbf{1}_{\{\max(i,j)\geq p,\ i\equiv j\pmod{p}% \}},

including the double pole of bib_{i} for i≥pi\geq p. A VV entry has a simple pole when p∣j+ℓp\mid j+\ell. As p2>Up^{2}>U, there are no higher powers.

In the first range every BB row has potential at least one, and two if i≥pi\geq p. A collision j=i+pj=i+p has j≥p+2⁢nj\geq p+2n, and its column gives the second unit. In the second range all i<pi<p. The columns j≥pj\geq p cover the simple BB poles and the CC poles. A collision needs i+p<10⁢ni+p<10n, and its row gives the other unit. In the third range i<pi<p and i+p≥10⁢ni+p\geq 10n, so only simple BB poles remain. In these three ranges every VV row has potential one, which covers j+ℓ=pj+\ell=p and j+ℓ=2⁢pj+\ell=2p. In the fourth range i,j<pi,j<p, only VV poles remain, and the rows with λV=1\lambda_{V}=1 are those containing one. So λrow+χj\lambda_{\rm row}+\chi_{j} dominates every pole order. Multiply each row and column by pp to the power of its potential. The potentials add up to (10⁢n−p)+n+(8⁢n−p)(10n-p)+n+(8n-p), (8⁢n−p)+n+(10⁢n−p)(8n-p)+n+(10n-p), n+(10⁢n−p)n+(10n-p) and 11⁢n−p11n-p in the four ranges. ∎

3.3. The factor at two

Lemma 3.4.

Let Nm=⌈53⁢m2/48⌉−mN_{m}=\lceil 53m^{2}/48\rceil-m. Then

(3.9) Pm⁢(a,T)∈aNm⁢ℤ⁢[a,T].P_{m}(a,T)\in a^{N_{m}}\mathbb{Z}[a,T].
Proof.

Work in ℚ⁢[[a]]⁢[T]\mathbb{Q}[[a]][T]. Write

Sj⁢(a)=∑u≥1auuj,vk⁢(i)={(i−k)−1,i≠k,0,i=k,zk⁢(i)=vk⁢(i)2,ek⁢(i)=𝟏i=k.S_{j}(a)=\sum_{u\geq 1}\frac{a^{u}}{u^{j}},\qquad v_{k}(i)=\begin{cases}(i-k)^% {-1},&i\neq k,\\ 0,&i=k,\end{cases}\quad z_{k}(i)=v_{k}(i)^{2},\quad e_{k}(i)=\mathbf{1}_{i=k}.

Subtract S1⁢(a)S_{1}(a) times each CC row from the paired row of B−T⁢AB-TA. Comparison of coefficients, including the diagonal, gives

(3.10) (CB−T⁢A−S1⁢C)\displaystyle\begin{pmatrix}C\\ B-TA-S_{1}C\end{pmatrix} =∑k≥0ak⁢(vk⁢(I)−zk⁢(I)−(T+S2)⁢ek⁢(I))⁢ek⁢(J)t\displaystyle=\sum_{k\geq 0}a^{k}\binom{v_{k}(I)}{-z_{k}(I)-(T+S_{2})e_{k}(I)}% e_{k}(J)^{t}
+∑k≥0ak⁢(ek⁢(I)vk⁢(I))⁢vk⁢(J)t+∑k≥0ak⁢(0−ek⁢(I))⁢zk⁢(J)t.\displaystyle\quad+\sum_{k\geq 0}a^{k}\binom{e_{k}(I)}{v_{k}(I)}v_{k}(J)^{t}+% \sum_{k\geq 0}a^{k}\binom{0}{-e_{k}(I)}z_{k}(J)^{t}.

For example, the off-diagonal lower entry is ∑k≥0ak⁢vk⁢(i)⁢vk⁢(j)−(ai+aj)/(i−j)2\sum_{k\geq 0}a^{k}v_{k}(i)v_{k}(j)-(a^{i}+a^{j})/(i-j)^{2}, and the diagonal entry is ∑k≥0ak⁢vk⁢(i)2−ai⁢(T+S2)=bi−T⁢ai\sum_{k\geq 0}a^{k}v_{k}(i)^{2}-a^{i}(T+S_{2})=b_{i}-Ta^{i}. Each coefficient of aa involves finitely many terms.

Expand along the nn rows of VV, and expand each remaining 2⁢m2m-row minor by the rank-one summands in (3.10). A nonzero term selects 2⁢m=8⁢n2m=8n distinct rank-one summands. At weight kk at most

ck=1+𝟏k∈J+𝟏k∈Ic_{k}=1+\mathbf{1}_{k\in J}+\mathbf{1}_{k\in I}

occur: the first summand needs k∈Jk\in J, the third needs k∈Ik\in I. All entries, including T+S2T+S_{2}, have nonnegative order in aa. So every term has order at least the sum of the 8⁢n8n smallest available weights. This use of the infinite expansion is coefficientwise: modulo aNm+1a^{N_{m}+1} only finitely many kk occur, so Cauchy–Binet applies to a finite truncation and then gives the asserted formal divisibility.

There is one weight at each 0≤k<n0\leq k<n, two at each n≤k<2⁢nn\leq k<2n, and three at each 2⁢n≤k<6⁢n2n\leq k<6n. Write 5⁢n=3⁢d+ϵ05n=3d+\epsilon_{0}, 0≤ϵ0<30\leq\epsilon_{0}<3. The minimum available-weight sum is

∑k=0n−1k+2⁢∑k=n2⁢n−1k+3⁢∑k=2⁢n2⁢n+d−1k+ϵ0⁢(2⁢n+d)\displaystyle\sum_{k=0}^{n-1}k+2\sum_{k=n}^{2n-1}k+3\sum_{k=2n}^{2n+d-1}k+% \epsilon_{0}(2n+d) =53⁢n23−4⁢n+ϵ0⁢(3−ϵ0)6\displaystyle=\frac{53n^{2}}{3}-4n+\frac{\epsilon_{0}(3-\epsilon_{0})}{6}
=⌈53⁢n23⌉−4⁢n=Nm.\displaystyle=\left\lceil\frac{53n^{2}}{3}\right\rceil-4n=N_{m}.

The minors of VV do not depend on a,Ta,T, so aNm∣Gm⁢(T;a)a^{N_{m}}\mid G_{m}(T;a) in ℚ⁢[a,T]\mathbb{Q}[a,T]. By (3.1), division by the monic aNma^{N_{m}} gives (3.9). ∎

3.4. Summation over primes

Lemma 3.5.

For every m=4⁢nm=4n, there is a rational number Dm∗>0D_{m}^{*}>0 with Dm∗⁢Gm⁢(T)∈ℤ⁢[T]D_{m}^{*}G_{m}(T)\in\mathbb{Z}[T] and

log⁡Dm∗≤τ⁢m2+o⁢(m2),τ=1881503369600−5348⁢log⁡2=4.325296633153372⁢…,\log D_{m}^{*}\leq\tau m^{2}+o(m^{2}),\qquad\tau=\frac{1881503}{369600}-\frac{% 53}{48}\log 2=4.325296633153372\ldots,

so that 4.36−τ>0.03474.36-\tau>0.0347.

Proof.

Put e0⁢(p)=3⁢m⁢⌊logp⁡N⌋+n⁢⌊logp⁡U⌋=vp⁢(𝒟m)e_{0}(p)=3m\lfloor\log_{p}N\rfloor+n\lfloor\log_{p}U\rfloor=v_{p}(\mathcal{D}_% {m}). For odd pp, with the first line taking precedence, let

(3.11) Ep={e0⁢(p),p≤m/5⁢or⁢p2≤U⁢or⁢p/m∈𝒞⁢or⁢p=5⁢m/4,m⁢𝔠⁢(p/m),m/5<p<5⁢m/4,Epout,5⁢m/4<p≤U,0,p>U,E_{p}=\begin{cases}e_{0}(p),&p\leq m/5\ \text{or}\ p^{2}\leq U\ \text{or}\ p/m% \in\mathcal{C}\ \text{or}\ p=5m/4,\\ m\mathfrak{c}(p/m),&m/5<p<5m/4,\\ E_{p}^{\rm out},&5m/4<p\leq U,\\ 0,&p>U,\end{cases}

and

Dm∗=2e0⁢(2)−Nm⁢∏p⁢oddpEp.D_{m}^{*}=2^{e_{0}(2)-N_{m}}\prod_{p\ {\rm odd}}p^{E_{p}}.

By (3.6), m⁢𝔠⁢(p/m)∈ℤm\mathfrak{c}(p/m)\in\mathbb{Z}. Lemmas 3.1, 3.2, 3.3 and 3.4 give vp⁢(Dm∗⁢Gm)≥0v_{p}(D_{m}^{*}G_{m})\geq 0 for every pp. The exponent at two may be negative.

For all sufficiently large mm, p>m/5p>m/5 implies p2>Up^{2}>U. Put

ϑ⁢(y)=∑p≤ylog⁡p,ηPNT⁢(X)=supy≥X|ϑ⁢(y)/y−1|.\vartheta(y)=\sum_{p\leq y}\log p,\qquad\eta_{\rm PNT}(X)=\sup_{y\geq X}|% \vartheta(y)/y-1|.

By the prime number theorem, ηPNT⁢(X)→0\eta_{\rm PNT}(X)\to 0. For an affine ff on [t0,t1]⊂[1/5,11/4][t_{0},t_{1}]\subset[1/5,11/4], Stieltjes integration by parts gives

m⁢∑t0⁢m<p<t1⁢mf⁢(p/m)⁢log⁡p=m2⁢∫t0t1f⁢(ρ)⁢𝑑ρ+O⁢(ηPNT⁢(m/5)⁢m2+m⁢log⁡m).m\sum_{t_{0}m<p<t_{1}m}f(p/m)\log p=m^{2}\int_{t_{0}}^{t_{1}}f(\rho)\,d\rho+O% \bigl{(}\eta_{\rm PNT}(m/5)m^{2}+m\log m\bigr{)}.

There are finitely many pieces, and the excluded primes at cut points contribute O⁢(m⁢log⁡m)O(m\log m). For p≤m/5p\leq m/5 the first powers contribute (13⁢m/4)⁢ϑ⁢(m/5)=(13/20)⁢m2+o⁢(m2)(13m/4)\vartheta(m/5)=(13/20)m^{2}+o(m^{2}), and the higher powers, including those of two,

O⁢(m⁢∑j=2⌊log2⁡U⌋ϑ⁢(U1/j))=O⁢(m3/2⁢log⁡(3⁢m)).O\Bigl{(}m\sum_{j=2}^{\lfloor\log_{2}U\rfloor}\vartheta(U^{1/j})\Bigr{)}=O(m^{% 3/2}\log(3m)).

Subtract Nm⁢log⁡2N_{m}\log 2 once. For p>5⁢m/4p>5m/4 the exponent is m⁢𝔠out⁢(p/m)m\mathfrak{c}_{\rm out}(p/m) with

𝔠out⁢(ρ)={19/4−2⁢ρ,5/4<ρ<2,11/4−ρ,2≤ρ≤11/4,∫5/411/4𝔠out⁢(ρ)⁢𝑑ρ=98+932=4532.\mathfrak{c}_{\rm out}(\rho)=\begin{cases}19/4-2\rho,&5/4<\rho<2,\\ 11/4-\rho,&2\leq\rho\leq 11/4,\end{cases}\qquad\int_{5/4}^{11/4}\mathfrak{c}_{% \rm out}(\rho)\,d\rho=\frac{9}{8}+\frac{9}{32}=\frac{45}{32}.

With (3.8),

1320+1121513369600+4532−5348⁢log⁡2=τ.\frac{13}{20}+\frac{1121513}{369600}+\frac{45}{32}-\frac{53}{48}\log 2=\tau.

The total error is O⁢(ηPNT⁢(m/5)⁢m2+m3/2⁢log⁡(3⁢m))=o⁢(m2)O(\eta_{\rm PNT}(m/5)m^{2}+m^{3/2}\log(3m))=o(m^{2}). ∎

4. Analytic estimate

In (2.4) there are mm, L=5⁢m/4L=5m/4 and mm variables xx, vv, zz. Let

𝒦x=𝒦v=[0,1],𝒦z=[−2,−r]∪[r,1]∪{z:|z|=r}.\mathcal{K}_{x}=\mathcal{K}_{v}=[0,1],\qquad\mathcal{K}_{z}=[-2,-r]\cup[r,1]% \cup\{z:|z|=r\}.

For i∈{x,v,z}i\in\{x,v,z\}, let

𝒦i+={q∈ℂ:dist⁡(q,𝒦i)≤3/4096}.\mathcal{K}_{i}^{+}=\{q\in\mathbb{C}:\operatorname{dist}(q,\mathcal{K}_{i})% \leq 3/4096\}.

The sets 𝒦z+\mathcal{K}_{z}^{+} avoid zero, and log⁡|2+x⁢v|\log|2+xv| is bounded on 𝒦x+×𝒦v+\mathcal{K}_{x}^{+}\times\mathcal{K}_{v}^{+}. For positive measures μ=(μx,μv,μz)\mu=(\mu_{x},\mu_{v},\mu_{z}) of masses (1,5/4,1)(1,5/4,1) put

I⁢(σ,σ′)\displaystyle I(\sigma,\sigma^{\prime}) =∬log⁡|u−v|⁢d⁢σ⁢(u)⁢𝑑σ′⁢(v),D⁢(σ,σ′)=∬log⁡|2+u⁢v|⁢d⁢σ⁢(u)⁢𝑑σ′⁢(v),\displaystyle=\iint\log|u-v|\,d\sigma(u)\,d\sigma^{\prime}(v),\qquad D(\sigma,% \sigma^{\prime})=\iint\log|2+uv|\,d\sigma(u)\,d\sigma^{\prime}(v),
𝒬⁢(μ)\displaystyle\mathcal{Q}(\mu) =I⁢(μx,μx)+I⁢(μv,μv)+I⁢(μz,μz)+I⁢(μv,μz)−D⁢(μx,μv),\displaystyle=I(\mu_{x},\mu_{x})+I(\mu_{v},\mu_{v})+I(\mu_{z},\mu_{z})+I(\mu_{% v},\mu_{z})-D(\mu_{x},\mu_{v}),
ℰ⁢(μ)\displaystyle\mathcal{E}(\mu) =32⁢log⁡2+𝒬⁢(μ)+34⁢∫log⁡|u|⁢d⁢μx+14⁢∫log⁡|u|⁢d⁢μv−54⁢∫log⁡|u|⁢d⁢μz.\displaystyle=\frac{3}{2}\log 2+\mathcal{Q}(\mu)+\frac{3}{4}\int\log|u|\,d\mu_% {x}+\frac{1}{4}\int\log|u|\,d\mu_{v}-\frac{5}{4}\int\log|u|\,d\mu_{z}.

For a fixed ν\nu let 𝒰i⁢(q)=∫log⁡|q−t|⁢d⁢νi⁢(t)\mathcal{U}_{i}(q)=\int\log|q-t|\,d\nu_{i}(t), Di⁢(q)=∫log⁡|2+q⁢t|⁢d⁢νi⁢(t)D_{i}(q)=\int\log|2+qt|\,d\nu_{i}(t), and

Fx⁢(q)\displaystyle F_{x}(q) =2⁢𝒰x⁢(q)−Dv⁢(q)+34⁢log⁡|q|,\displaystyle=2\mathcal{U}_{x}(q)-D_{v}(q)+\tfrac{3}{4}\log|q|,
Fv⁢(q)\displaystyle F_{v}(q) =2⁢𝒰v⁢(q)+𝒰z⁢(q)−Dx⁢(q)+14⁢log⁡|q|,\displaystyle=2\mathcal{U}_{v}(q)+\mathcal{U}_{z}(q)-D_{x}(q)+\tfrac{1}{4}\log% |q|,
Fz⁢(q)\displaystyle F_{z}(q) =2⁢𝒰z⁢(q)+𝒰v⁢(q)−54⁢log⁡|q|.\displaystyle=2\mathcal{U}_{z}(q)+\mathcal{U}_{v}(q)-\tfrac{5}{4}\log|q|.
Lemma 4.1.

On triples of positive measures μ=(μx,μv,μz)\mu=(\mu_{x},\mu_{v},\mu_{z}) of masses (1,5/4,1)(1,5/4,1) with finite absolute logarithmic self-integrals and supp⁡μx⊆𝒦x+\operatorname{supp}\mu_{x}\subseteq\mathcal{K}_{x}^{+}, supp⁡μv⊆𝒦v+\operatorname{supp}\mu_{v}\subseteq\mathcal{K}_{v}^{+}, and supp⁡μz⊆𝒦z+\operatorname{supp}\mu_{z}\subseteq\mathcal{K}_{z}^{+}, respectively, the functional 𝒬\mathcal{Q} is concave. The same statement applies when all supports lie in the smaller sets 𝒦x,𝒦v,𝒦z\mathcal{K}_{x},\mathcal{K}_{v},\mathcal{K}_{z}.

(4.1) ℰ⁢(μ)≤32⁢log⁡2−𝒬⁢(ν)+∫Fx⁢𝑑μx+∫Fv⁢𝑑μv+∫Fz⁢𝑑μz.\mathcal{E}(\mu)\leq\frac{3}{2}\log 2-\mathcal{Q}(\nu)+\int F_{x}\,d\mu_{x}+% \int F_{v}\,d\mu_{v}+\int F_{z}\,d\mu_{z}.
Proof.

For positive compactly supported measures η1,η2\eta_{1},\eta_{2} in a disk DD strictly containing their supports, the Green kernel of DD is −log⁡|u−v|-\log|u-v| plus a bounded term. Cauchy–Schwarz for its energy, together with the finite absolute self-integrals, gives absolute logarithmic cross-integrability. Thus signed differences of the measures used below have finite absolute logarithmic energy.

For a real signed measure σ=σ+−σ−\sigma=\sigma_{+}-\sigma_{-} of mass zero, with positive compactly supported parts of finite absolute logarithmic self-integrals, write IST=−II_{\rm ST}=-I for the sign convention in Saff–Totik. Their positivity lemma [9, Lemma I.1.8] gives

‖σ‖2=−I⁢(σ,σ)≥0.\|\sigma\|^{2}=-I(\sigma,\sigma)\geq 0.

For smooth σ\sigma, the identity

−I⁢(σ,σ)=12⁢π⁢∫ℂ|∇𝒰σ|2⁢𝑑A-I(\sigma,\sigma)=\frac{1}{2\pi}\int_{\mathbb{C}}|\nabla\mathcal{U}^{\sigma}|^% {2}\,dA

follows from Δ⁢𝒰σ=2⁢π⁢σ\Delta\mathcal{U}^{\sigma}=2\pi\sigma and integration by parts. The boundary term vanishes because the mass is zero. Approximation extends this Dirichlet-integral formula to the stated finite-energy measures. Put ⟨σ,σ′⟩=−I⁢(σ,σ′)\langle\sigma,\sigma^{\prime}\rangle=-I(\sigma,\sigma^{\prime}).

Let (σx,σv,σz)(\sigma_{x},\sigma_{v},\sigma_{z}) have mass zero in each coordinate, and first let the vv-support avoid zero. Put R⁢(v)=−2/vR(v)=-2/v. Then

log⁡|2+x⁢v|\displaystyle\log|2+xv| =log⁡|x−R⁢(v)|+log⁡|v|,\displaystyle=\log|x-R(v)|+\log|v|,
log⁡|R⁢(v)−R⁢(z)|\displaystyle\log|R(v)-R(z)| =log⁡2+log⁡|v−z|−log⁡|v|−log⁡|z|,\displaystyle=\log 2+\log|v-z|-\log|v|-\log|z|,

and the one-variable terms integrate to zero. With σx′=σx\sigma^{\prime}_{x}=\sigma_{x}, σv′=R∗⁢σv\sigma^{\prime}_{v}=R_{*}\sigma_{v} and σz′=−R∗⁢σz\sigma^{\prime}_{z}=-R_{*}\sigma_{z},

𝒬⁢(σ)\displaystyle\mathcal{Q}(\sigma) =−‖σx′‖2−‖σv′‖2−‖σz′‖2+⟨σx′,σv′⟩+⟨σv′,σz′⟩\displaystyle=-\|\sigma^{\prime}_{x}\|^{2}-\|\sigma^{\prime}_{v}\|^{2}-\|% \sigma^{\prime}_{z}\|^{2}+\langle\sigma^{\prime}_{x},\sigma^{\prime}_{v}% \rangle+\langle\sigma^{\prime}_{v},\sigma^{\prime}_{z}\rangle
=−12⁢(‖σx′−σv′‖2+‖σv′−σz′‖2+‖σx′‖2+‖σz′‖2)≤0.\displaystyle=-\tfrac{1}{2}\bigl{(}\|\sigma^{\prime}_{x}-\sigma^{\prime}_{v}\|% ^{2}+\|\sigma^{\prime}_{v}-\sigma^{\prime}_{z}\|^{2}+\|\sigma^{\prime}_{x}\|^{% 2}+\|\sigma^{\prime}_{z}\|^{2}\bigr{)}\leq 0.

The zz-supports here avoid zero. In general fix a probability λ\lambda with smooth density supported in [1/2,3/4][1/2,3/4], and put

σv,ζ=σv|{|v|>ζ}−σv⁢({|v|>ζ})⁢λ.\sigma_{v,\zeta}=\sigma_{v}|_{\{|v|>\zeta\}}-\sigma_{v}(\{|v|>\zeta\})\lambda.

It has mass zero and support away from zero, so 𝒬⁢(σx,σv,ζ,σz)≤0\mathcal{Q}(\sigma_{x},\sigma_{v,\zeta},\sigma_{z})\leq 0. Finite energy excludes an atom at zero, so the coefficient of λ\lambda tends to zero as ζ→0\zeta\to 0. The cross-integrability just proved makes the logarithmic kernels integrable against the total variations, and the potential of λ\lambda is bounded on compact sets. Dominated convergence in the original coordinates gives

𝒬⁢(σx,σv,ζ,σz)⟶𝒬⁢(σx,σv,σz)≤0.\mathcal{Q}(\sigma_{x},\sigma_{v,\zeta},\sigma_{z})\longrightarrow\mathcal{Q}(% \sigma_{x},\sigma_{v},\sigma_{z})\leq 0.

The kernel log⁡|2+x⁢v|\log|2+xv| is bounded on the sets specified above. For the reference and the smoothed measures below also ∫|log⁡|v|⁢|d|⁢σv|⁢(v)<∞\int|\log|v||\,d|\sigma_{v}|(v)<\infty: the reference has bounded interval density, and the smoothed measures stay away from zero. With

|log⁡|R⁢(v)−R⁢(w)||≤log⁡2+|log⁡|v−w||+|log⁡|v||+|log⁡|w|||\log|R(v)-R(w)||\leq\log 2+|\log|v-w||+|\log|v||+|\log|w||

and the first kernel identity, dominated convergence also applies to the inverted pairings.

Expand 𝒬⁢(ν+σ)\mathcal{Q}(\nu+\sigma) with σ=μ−ν\sigma=\mu-\nu. The quadratic part is nonpositive and d⁢𝒬ν⁢(ν)=2⁢𝒬⁢(ν)d\mathcal{Q}_{\nu}(\nu)=2\mathcal{Q}(\nu). Adding the linear terms gives (4.1). ∎

Lemma 4.2 (computer-assisted).

There is a positive ν\nu of masses (1,5/4,1)(1,5/4,1), with 5050, 5252, 6262 constant-density intervals on [0,1][0,1], [0,1][0,1], [−2,−3/5][-2,-3/5] and 157157 Poisson kernels on |z|=3/5|z|=3/5, all with rational data, such that

Fx≤Kx=−4.626,Fv≤Kv=−4.480,Fz≤Kz=−0.090F_{x}\leq K_{x}=-4.626,\qquad F_{v}\leq K_{v}=-4.480,\qquad F_{z}\leq K_{z}=-0% .090

on [0,1][0,1], [0,1][0,1] and [−2,−3/5]∪[3/5,1]∪{|z|=3/5}[-2,-3/5]\cup[3/5,1]\cup\{|z|=3/5\}, 𝒬⁢(ν)\mathcal{Q}(\nu) satisfies (B.1), and, on the disk |q|≤δ=1/250|q|\leq\delta=1/250,

(4.2) Fx⁢(q)<−6.2823007011,Fv⁢(q)<−4.4935082437.F_{x}(q)<-6.2823007011,\qquad F_{v}(q)<-4.4935082437.

The interval densities are at most 131,153,43131,153,43, and the Poisson densities (4.4) add up to at most 2424.

Proof.

The rational data of ν\nu are printed in Tables D.1 and D.2. An interval [t0,t1][t_{0},t_{1}] of mass ww carries w⁢d⁢t/(t1−t0)w\,dt/(t_{1}-t_{0}). A Poisson entry (c,w)(c,w), |c|<r|c|<r, is

w⁢(r2−|c|2)2⁢π⁢|r⁢ei⁢θ−c|2⁢d⁢θ(0≤θ≤2⁢π).\frac{w(r^{2}-|c|^{2})}{2\pi|re^{i\theta}-c|^{2}}\,d\theta\qquad(0\leq\theta% \leq 2\pi).

The file anc/analytic/reference.json repeats these tables for the verification program. The weights are positive, the three masses are exact, and the Poisson entries are invariant under conjugation. An interval has potential

wt1−t0⁢Re⁡{f⁢(q−t0)−f⁢(q−t1)},f⁢(u)=u⁢log⁡u−u,\frac{w}{t_{1}-t_{0}}\operatorname{Re}\{f(q-t_{0})-f(q-t_{1})\},\qquad f(u)=u% \log u-u,

with u⁢log⁡|u|−uu\log|u|-u on the real line. A Poisson entry has potential

(4.3) w⁢{log⁡|r2−q⁢c¯|−log⁡r,|q|≤r,log⁡|q−c|,|q|≥r.w\begin{cases}\log|r^{2}-q\overline{c}|-\log r,&|q|\leq r,\\ \log|q-c|,&|q|\geq r.\end{cases}

Let UP⁢(q)U_{P}(q) denote the Poisson contribution to the zz-potential. Put

MP=∑cwc=2646000950940000000000,WP:=MP⁢log⁡(151/250).M_{P}=\sum_{c}w_{c}=\frac{26460009509}{40000000000},\qquad W_{P}:=M_{P}\log(15% 1/250).

Since r+δ=151/250r+\delta=151/250, the 157 Poisson rows give

sup|q|≤δUP⁢(q)≤WP.\sup_{|q|\leq\delta}U_{P}(q)\leq W_{P}.

Real interval energies use g⁢(u)=12⁢u2⁢log⁡|u|−34⁢u2g(u)=\tfrac{1}{2}u^{2}\log|u|-\tfrac{3}{4}u^{2}, with g′′⁢(u)=log⁡|u|g^{\prime\prime}(u)=\log|u|. For the denominator,

log⁡(2+x⁢v)=log⁡2+∑j=1420(−1)j+1⁢(x⁢v/2)jj+R420,|R420|≤2⁢(1/2)421421.\log(2+xv)=\log 2+\sum_{j=1}^{420}\frac{(-1)^{j+1}(xv/2)^{j}}{j}+R_{420},% \qquad|R_{420}|\leq\frac{2(1/2)^{421}}{421}.

The pair decomposition and its outward enclosures are printed in Table C.1. Since the product mass of νx⊗νv\nu_{x}\otimes\nu_{v} is 1⋅(5/4)1\cdot(5/4), the integrated remainder is at most

54⁢2⁢(1/2)421421=52⋅421⋅2421<1.10×10−129.\frac{5}{4}\frac{2(1/2)^{421}}{421}=\frac{5}{2\cdot 421\cdot 2^{421}}<1.10% \times 10^{-129}.

The bounds for Fx,Fv,FzF_{x},F_{v},F_{z} are proved by subdivision in Arb ball arithmetic [5]. On a real subinterval, each interval-potential component is quasi-convex: away from t0,t1t_{0},t_{1}, its derivative has the sign of log⁡|q−t0|−log⁡|q−t1|\log|q-t_{0}|-\log|q-t_{1}|, which changes from negative to positive at (t0+t1)/2(t_{0}+t_{1})/2, while the negative denominator term is convex. On each side of |q|=r|q|=r a Poisson-entry potential is the logarithm of the distance from qq to a fixed point, hence quasi-convex, and a row containing r=3/5r=3/5 also evaluates q=rq=r. Thus the component maximum is bounded by the relevant endpoint values. On a subarc of the circle, the largest distance to a pole is at an endpoint or at the antipodal point. Interval sign tests exclude the antipodal point when possible. If it lies on the subarc, or if a sign test is inconclusive, the full-circle upper bound is used. Otherwise the endpoint values suffice. Table B.1 covers the real paths and the upper semicircle, and conjugation covers the lower one. On the subintervals meeting zero the positive term 34⁢log⁡|q|\frac{3}{4}\log|q| or 14⁢log⁡|q|\frac{1}{4}\log|q| is increasing in |q||q| and equals −∞-\infty at zero. For a real pole tt and q=r⁢ei⁢θq=re^{i\theta}, |r⁢ei⁢θ−t|2=r2+t2−2⁢r⁢t⁢cos⁡θ|re^{i\theta}-t|^{2}=r^{2}+t^{2}-2rt\cos\theta. This is monotone in θ∈[0,π]\theta\in[0,\pi] for each fixed sign of tt, which justifies the endpoint or antipodal alternatives.

On |q|≤δ|q|\leq\delta use |q−t|≤|t|+δ|q-t|\leq|t|+\delta, |2+q⁢t|≥2−δ|2+qt|\geq 2-\delta and log⁡|r2−q⁢c¯|−log⁡r≤log⁡(r+δ)\log|r^{2}-q\overline{c}|-\log r\leq\log(r+\delta) with the same primitives. This gives (4.2). An interval has density w/(t1−t0)w/(t_{1}-t_{0}). With respect to arclength, a Poisson entry has density at most

(4.4) Hc=w⁢(r+|c|)2⁢π⁢r⁢(r−|c|).H_{c}=\frac{w(r+|c|)}{2\pi r(r-|c|)}.

The endpoint majorants and disk bounds are printed in Table C.2 and the density checks in Table C.3. The density checks are rational, by an integer square-root bound for each |c|2|c|^{2} and by π>333/106\pi>333/106, which follows from Machin’s formula:

16⁢(15−13⋅53+15⋅55−17⋅57)−4239=333106+341821341563593750.∎16\left(\frac{1}{5}-\frac{1}{3\cdot 5^{3}}+\frac{1}{5\cdot 5^{5}}-\frac{1}{7% \cdot 5^{7}}\right)-\frac{4}{239}=\frac{333}{106}+\frac{3418213}{41563593750}.\qed

For this ν\nu put

α=𝒬⁢(ν)−32⁢log⁡2−Kx−54⁢Kv−Kz;\alpha=\mathcal{Q}(\nu)-\tfrac{3}{2}\log 2-K_{x}-\tfrac{5}{4}K_{v}-K_{z};

by (B.2), α>4.365\alpha>4.365.

Lemma 4.3.

Put (Hx,Hv,Hz)=(131,153,67)(H_{x},H_{v},H_{z})=(131,153,67). If 0<d=|q−q′|≤1/160<d=|q-q^{\prime}|\leq 1/16, then

|𝒰i⁢(q)−𝒰i⁢(q′)|≤8⁢Hi⁢d⁢(1+|log⁡d|).|\mathcal{U}_{i}(q)-\mathcal{U}_{i}(q^{\prime})|\leq 8H_{i}d(1+|\log d|).
Proof.

Let E⊂ℝE\subset\mathbb{R} be a measurable set of measure ℓE≤9/5\ell_{E}\leq 9/5 containing the real support. By the triangle inequality,

|log⁡|q−t|−log⁡|q′−t||≤log⁡(1+d|q−t|)+log⁡(1+d|q′−t|).|\log|q-t|-\log|q^{\prime}-t||\leq\log\Bigl{(}1+\frac{d}{|q-t|}\Bigr{)}+\log% \Bigl{(}1+\frac{d}{|q^{\prime}-t|}\Bigr{)}.

Replacing qq or q′q^{\prime} by its real part increases each integral. Symmetric decreasing rearrangement gives

∫Elog⁡(1+d|q−t|)⁢𝑑t≤2⁢∫0ℓE/2log⁡(1+d/t)⁢𝑑t≤2⁢d⁢(1+log⁡(1+ℓE2⁢d)),\int_{E}\log\Bigl{(}1+\frac{d}{|q-t|}\Bigr{)}\,dt\leq 2\int_{0}^{\ell_{E}/2}% \log(1+d/t)\,dt\leq 2d\Bigl{(}1+\log\Bigl{(}1+\frac{\ell_{E}}{2d}\Bigr{)}\Bigr% {)},

so

∫E|log⁡|q−t|−log⁡|q′−t||⁢𝑑t≤4⁢d⁢(1+log⁡(1+ℓE2⁢d))≤4⁢d⁢(1+log⁡(1+910⁢d))<8⁢d⁢(1+|log⁡d|),\int_{E}|\log|q-t|-\log|q^{\prime}-t||\,dt\leq 4d\Bigl{(}1+\log\Bigl{(}1+\frac% {\ell_{E}}{2d}\Bigr{)}\Bigr{)}\leq 4d\Bigl{(}1+\log\Bigl{(}1+\frac{9}{10d}% \Bigr{)}\Bigr{)}<8d(1+|\log d|),

by d2+9⁢d/5≤149/1280<1d^{2}+9d/5\leq 149/1280<1. The interval components have disjoint interiors, so multiply by the density bound. By (4.3), a Poisson component is Lipschitz with constant at most w/(r−|c|)w/(r-|c|), and by (4.4)

wr−|c|=2⁢π⁢rr+|c|⁢Hc≤2⁢π⁢Hc<8⁢Hc.\frac{w}{r-|c|}=\frac{2\pi r}{r+|c|}H_{c}\leq 2\pi H_{c}<8H_{c}.

Summation and 43+24=6743+24=67 give the bound for 𝒰z\mathcal{U}_{z}. ∎

Lemma 4.4.

If m=4⁢n≥64m=4n\geq 64, then

log⁡|Gm⁢(β)|≤−4.36⁢m2+870⁢m⁢log⁡m.\log|G_{m}(\beta)|\leq-4.36m^{2}+870m\log m.
Proof.

Expand Γm\Gamma^{m} in (2.4) into products of its four oriented pieces. On each piece use normalized arclength, and keep its coefficient and tangent phase until taking absolute values. Since |cΓ−1|=|cΓ||c_{\Gamma}-1|=|c_{\Gamma}|, the absolute weights add up to

ϰ=(2−r)+(1−r)+2⁢π⁢r⁢|cΓ|=95+35⁢π2+log2⁡2,\varkappa=(2-r)+(1-r)+2\pi r|c_{\Gamma}|=\frac{9}{5}+\frac{3}{5}\sqrt{\pi^{2}+% \log^{2}2},

so the paths contribute a factor at most ϰm\varkappa^{m}. Only positive measures enter Lemma 4.1.

Put ϵ=m−2\epsilon=m^{-2} and η=3⁢ϵ\eta=3\epsilon. Replace every atom uu by the uniform measure of mass 1/m1/m on the circle u+2⁢i⁢ϵ+ϵ⁢ei⁢θu+2i\epsilon+\epsilon e^{i\theta}. The resulting μx,μv,μz\mu_{x},\mu_{v},\mu_{z} have masses (1,5/4,1)(1,5/4,1) and finite energy. The common translation preserves differences, so by subharmonicity in each variable averaging increases every positive logarithmic pair. Each circle has energy log⁡ϵ\log\epsilon, so removing the diagonal adds

−(m+5⁢m/4+m)⁢log⁡ϵ=132⁢m⁢log⁡m.-(m+5m/4+m)\log\epsilon=\frac{13}{2}m\log m.

The positive endpoint terms also increase, since

(4.5) 12⁢π⁢∫02⁢πlog⁡|u+2⁢i⁢ϵ+ϵ⁢ei⁢θ|⁢d⁢θ=log⁡|u+2⁢i⁢ϵ|≥log⁡u(u>0).\frac{1}{2\pi}\int_{0}^{2\pi}\log|u+2i\epsilon+\epsilon e^{i\theta}|\,d\theta=% \log|u+2i\epsilon|\geq\log u\quad(u>0).

If an original factor vanishes, there is nothing to prove.

For m≥64m\geq 64 the negative kernels satisfy the derivative bounds

(4.6) LD=101100,Lz=167100,1+η2−2⁢η−η2<LD,1r−η<Lz,L_{D}=\frac{101}{100},\qquad L_{z}=\frac{167}{100},\qquad\frac{1+\eta}{2-2\eta% -\eta^{2}}<L_{D},\qquad\frac{1}{r-\eta}<L_{z},

since a displacement of at most η\eta in each coordinate gives |2+x⁢v|≥2−2⁢η−η2|2+xv|\geq 2-2\eta-\eta^{2} and |z|≥r−η|z|\geq r-\eta. So the negative terms contribute at most

A0=6⋅54⁢LD+3⋅54⁢Lz=110780.A_{0}=6\cdot\frac{5}{4}L_{D}+3\cdot\frac{5}{4}L_{z}=\frac{1107}{80}.

Discard the factorials and 2−m2^{-m} in (2.4). Then the logarithm of the integrand, with its prefactor, is at most

(4.7) m2⁢ℰ⁢(μ)+132⁢m⁢log⁡m+A0.m^{2}\mathcal{E}(\mu)+\frac{13}{2}m\log m+A_{0}.

We bound the three fields in (4.1) on the smoothed measures. Put ω=1+log⁡(m2/3)\omega=1+\log(m^{2}/3). Lemma 4.3 with d≤3⁢ϵd\leq 3\epsilon gives the coefficients

Cx=16⁢Hx=2096,Cv=16⁢Hv+8⁢Hz=2984,Cz=16⁢Hz+8⁢Hv=2296,\displaystyle C_{x}=16H_{x}=2096,\qquad C_{v}=16H_{v}+8H_{z}=2984,\qquad C_{z}% =16H_{z}+8H_{v}=2296,
Cx+54⁢Cv+Cz=8122,\displaystyle C_{x}+\tfrac{5}{4}C_{v}+C_{z}=8122,

and a weighted error at most 3⋅8122⁢ϵ⁢ω3\cdot 8122\,\epsilon\omega.

For an xx or vv center u<δ−3⁢ϵu<\delta-3\epsilon the whole circle lies in |q|<δ|q|<\delta, where (4.2) applies. Otherwise u≥δ−3⁢ϵ≥δ/2u\geq\delta-3\epsilon\geq\delta/2, and (4.5) gives

log⁡|u+2⁢i⁢ϵ|−log⁡u=12⁢log⁡(1+4⁢ϵ2/u2)≤8⁢ϵ2δ2.\log|u+2i\epsilon|-\log u=\frac{1}{2}\log(1+4\epsilon^{2}/u^{2})\leq\frac{8% \epsilon^{2}}{\delta^{2}}.

So the positive terms add at most (17/16)⁢8⁢ϵ2/δ2(17/16)8\epsilon^{2}/\delta^{2}. By (4.6), the negative terms in the three fields add at most A0⁢ϵA_{0}\epsilon. At m=64m=64, and hence for all larger mm,

1716⁢8⁢ϵδ2+A0≤2656252048+110780<200.\frac{17}{16}\frac{8\epsilon}{\delta^{2}}+A_{0}\leq\frac{265625}{2048}+\frac{1% 107}{80}<200.

So the total smoothing error is at most (3⋅8122+200)⁢ϵ⁢ω=24566⁢ϵ⁢ω<26000⁢ϵ⁢ω(3\cdot 8122+200)\epsilon\omega=24566\epsilon\omega<26000\epsilon\omega. Insert the lower endpoint q−q_{-} of (B.1) in (4.1). As

−α<32⁢log⁡2−q−+Kx+54⁢Kv+Kz<−4.36,-\alpha<\frac{3}{2}\log 2-q_{-}+K_{x}+\frac{5}{4}K_{v}+K_{z}<-4.36,

we get m2⁢ℰ⁢(μ)≤−4.36⁢m2+26000⁢ωm^{2}\mathcal{E}(\mu)\leq-4.36m^{2}+26000\omega. With (4.7) and the path factor,

log⁡|Gm⁢(β)|≤−4.36⁢m2+26000⁢(1+log⁡(m2/3))+132⁢m⁢log⁡m+m⁢log⁡ϰ+A0.\log|G_{m}(\beta)|\leq-4.36m^{2}+26000(1+\log(m^{2}/3))+\frac{13}{2}m\log m+m% \log\varkappa+A_{0}.

Divided by m⁢log⁡mm\log m, no positive term increases for m≥64m\geq 64. At m=64m=64 their sum is

C∗=26000⁢(1+log⁡(642/3))64⁢log⁡64+132+log⁡ϰlog⁡64+110780⋅64⁢log⁡64<810<870.∎C_{*}=\frac{26000(1+\log(64^{2}/3))}{64\log 64}+\frac{13}{2}+\frac{\log% \varkappa}{\log 64}+\frac{1107}{80\cdot 64\log 64}<810<870.\qed

5. Nonvanishing

Let ℤ(11)={u/v∈ℚ:11∤v}\mathbb{Z}_{(11)}=\{u/v\in\mathbb{Q}:11\nmid v\}, and let v11v_{11} be the 1111-adic valuation, v11⁢(0)=+∞v_{11}(0)=+\infty.

Lemma 5.1.

For m=4m=4,

(5.1) G4⁢(T)≡g⁢(T):=4+2⁢T+7⁢T2+3⁢T3+9⁢T4(mod11),G_{4}(T)\equiv g(T):=4+2T+7T^{2}+3T^{3}+9T^{4}\pmod{11},

and gg has no root in 𝔽11\mathbb{F}_{11}.

Proof.

Here I={2,…,5}I=\{2,\ldots,5\}, J={1,…,9}J=\{1,\ldots,9\}, Λ={1}\Lambda=\{1\}, and all denominators are prime to 1111. Modulo 1111 the matrix (2.1) is

(5.2) (501517023310919807559049460211403431081+7⁢T8225601883+8⁢T71040594273+6⁢T285105921029+10⁢T31936439287510).\begin{pmatrix}5&0&1&5&1&7&0&2&3\\ 3&1&0&9&1&9&8&0&7\\ 5&5&9&0&4&9&4&6&0\\ 2&1&1&4&0&3&4&3&10\\ 8&1+7T&8&2&2&5&6&0&1\\ 8&8&3+8T&7&10&4&0&5&9\\ 4&2&7&3+6T&2&8&5&10&5\\ 9&2&10&2&9+10T&3&1&9&3\\ 6&4&3&9&2&8&7&5&10\end{pmatrix}.

Its determinant is g⁢(T)g(T). Reduction of (2.5), where 16336404000000000≡4(mod11)16336404000000000\equiv 4\pmod{11}, gives the same polynomial. Table 5.1 lists the values of gg. None is zero. ∎

Table 5.1. The polynomial gg over 𝔽11\mathbb{F}_{11}.
ww 0 1 2 3 4 5 6 7 8 9 10
g⁢(w)g(w) 4 3 6 3 2 7 7 9 5 5 4
Lemma 5.2.

For every ξ∈ℚ\xi\in\mathbb{Q} and every integer t≥1t\geq 1, G4⋅11t⁢(ξ)≠0G_{4\cdot 11^{t}}(\xi)\neq 0.

Proof.

Put Qt=11tQ_{t}=11^{t}, mt=4⁢Qtm_{t}=4Q_{t} and h=σ=k=Qth=\sigma=k=Q_{t}, and

(5.3) Pt⁢(W)=Qt3⁢mt+h⁢Gmt⁢(W/Qt2).P_{t}(W)=Q_{t}^{3m_{t}+h}G_{m_{t}}(W/Q_{t}^{2}).

This is the determinant obtained after T=W/Qt2T=W/Q_{t}^{2} by multiplying the CC, BB and VV rows by QtQ_{t}, Qt2Q_{t}^{2} and QtQ_{t}. Every nonzero index, difference and Cauchy denominator has absolute value less than 11⁢Qt11Q_{t}, hence 1111-adic valuation at most tt. So Qt⁢tiQ_{t}t_{i}, Qt2⁢biQ_{t}^{2}b_{i} and all entries of the scaled matrix lie in ℤ(11)⁢[W]\mathbb{Z}_{(11)}[W]. Hence Pt∈ℤ(11)⁢[W]P_{t}\in\mathbb{Z}_{(11)}[W] and deg⁡Pt≤mt\deg P_{t}\leq m_{t}.

Let a=−2a=-2 and 0≤u<Qt0\leq u<Q_{t}. Only the summation indices divisible by QtQ_{t} survive, and aQt≡a(mod11)a^{Q_{t}}\equiv a\pmod{11}, so

Qt⁢tQt⁢j+u≡au⁢tj,Qt2⁢bQt⁢j+u≡au⁢bj(mod11).Q_{t}t_{Q_{t}j+u}\equiv a^{u}t_{j},\qquad Q_{t}^{2}b_{Q_{t}j+u}\equiv a^{u}b_{% j}\pmod{11}.

The scaled CC and BB entries vanish modulo 1111 between distinct residue classes modulo QtQ_{t}. In the class uu write i=Qt⁢i0+ui=Q_{t}i_{0}+u, j=Qt⁢j0+uj=Q_{t}j_{0}+u. Then 2≤i0≤52\leq i_{0}\leq 5 and 1≤j0≤91\leq j_{0}\leq 9. The only VV row in this class is ℓ=Qt−u\ell=Q_{t}-u, with

Qtj+ℓ=1j0+1.\frac{Q_{t}}{j+\ell}=\frac{1}{j_{0}+1}.

Cross-class CC and BB entries have 1111-adic valuation at least 11. Within a class the nonzero differences j0−i0j_{0}-i_{0} satisfy |j0−i0|≤7<11|j_{0}-i_{0}|\leq 7<11 and are units. All other scaled VV entries vanish modulo 1111. So, after permutations, the matrix has QtQ_{t} diagonal blocks

[au⁢C4au⁢(B4−W⁢A4)V4],0≤u<Qt,\begin{bmatrix}a^{u}C_{4}\\ a^{u}(B_{4}-WA_{4})\\ V_{4}\end{bmatrix},\qquad 0\leq u<Q_{t},

where the subscript 44 refers to G4G_{4}. Each block has eight rows with the factor aua^{u}. Hence

(5.4) Pt⁢(W)≡εt⁢(−2)4⁢Qt⁢(Qt−1)⁢g⁢(W)Qt(mod11),εt∈{1,−1}.P_{t}(W)\equiv\varepsilon_{t}(-2)^{4Q_{t}(Q_{t}-1)}g(W)^{Q_{t}}\pmod{11},% \qquad\varepsilon_{t}\in\{1,-1\}.

The constant coefficient and the coefficient of WmtW^{m_{t}} are 1111-adic units, and deg⁡Pt=mt\deg P_{t}=m_{t}.

Take W=Qt2⁢ξW=Q_{t}^{2}\xi. If v11⁢(W)>0v_{11}(W)>0, the constant term is the unique term of least valuation in Pt⁢(W)P_{t}(W). If v11⁢(W)=0v_{11}(W)=0, then (5.4) and Lemma 5.1 give Pt⁢(W)≢0(mod11)P_{t}(W)\not\equiv 0\pmod{11}. If v11⁢(W)<0v_{11}(W)<0, write Pt⁢(W)=∑j=0mtcj⁢WjP_{t}(W)=\sum_{j=0}^{m_{t}}c_{j}W^{j}. Here v11⁢(cmt)=0v_{11}(c_{m_{t}})=0 and v11⁢(cj)≥0v_{11}(c_{j})\geq 0, so

v11⁢(cj⁢Wj)>mt⁢v11⁢(W)=v11⁢(cmt⁢Wmt)(0≤j<mt),v_{11}(c_{j}W^{j})>m_{t}v_{11}(W)=v_{11}(c_{m_{t}}W^{m_{t}})\qquad(0\leq j<m_{% t}),

and the leading term is the unique term of least valuation. In all three cases Pt⁢(W)≠0P_{t}(W)\neq 0, and by (5.3) Gmt⁢(ξ)≠0G_{m_{t}}(\xi)\neq 0. ∎

6. Proof of Theorem 1

Suppose β=u/q0\beta=u/q_{0} with u∈ℤu\in\mathbb{Z} and q0∈ℤ>0q_{0}\in\mathbb{Z}_{>0}. By Lemmas 2.1 and 3.5,

Zm:=q0m⁢Dm∗⁢Gm⁢(u/q0)∈ℤ.Z_{m}:=q_{0}^{m}D_{m}^{*}G_{m}(u/q_{0})\in\mathbb{Z}.

By Lemma 5.2, Zm≠0Z_{m}\neq 0 for m=4⋅11tm=4\cdot 11^{t}, t≥1t\geq 1. By Lemmas 3.5 and 4.4,

log⁡|Zm|≤−(4.36−τ)⁢m2+870⁢m⁢log⁡m+o⁢(m2)+m⁢log⁡q0.\log|Z_{m}|\leq-(4.36-\tau)m^{2}+870m\log m+o(m^{2})+m\log q_{0}.

As 4.36−τ>0.03474.36-\tau>0.0347, the right-hand side tends to −∞-\infty along this sequence. So 0<|Zm|<10<|Z_{m}|<1 for large tt, which is impossible. ∎

Appendix A Arithmetic integration

In Table A.1 the second column lists the lines (3.7) in 0<x<ρ0<x<\rho, in increasing order of xx. The third column lists 𝔯\mathfrak{r} on the six xx-intervals they cut from [0,ρ)[0,\rho). Then 𝔠⁢(ρ)=1+∑𝔯⋅(length)\mathfrak{c}(\rho)=1+\sum\mathfrak{r}\cdot(\text{length}). In the last three rows the fifth value includes the correction −1-1 of (3.3). By (3.11), the values of 𝔠\mathfrak{c} at the cut points are not used. Of the lines (3.7), 2525 meet 0<x<ρ0<x<\rho for some ρ∈(1/5,5/4)\rho\in(1/5,5/4), and 2727 meet the closed interval after adding x=0x=0 and x=ρx=\rho. The table has 44⋅6=26444\cdot 6=264 regions and 6666 distinct triples (K,ℒ,W)(K,\mathcal{L},W).

Table A.1. Exact arithmetic integration.
ρ\rho lines in 0<x<ρ0<x<\rho 𝔯\mathfrak{r} 𝔠⁢(ρ)\mathfrak{c}(\rho) ∫𝔠⁢𝑑ρ\int\mathfrak{c}\,d\rho
(1/5,9/44)(1/5,9/44) 1/4−ρ, 5/2−12⁢ρ, 3/2−7⁢ρ, 1/2−2⁢ρ, 2⁢ρ−1/41/4-\rho,\,5/2-12\rho,\,3/2-7\rho,\,1/2-2\rho,\,2\rho-1/4 11,11,11,9,11,1111,11,11,9,11,11 3+ρ3+\rho 1409/968001409/96800
(9/44,5/24)(9/44,5/24) 5/2−12⁢ρ, 1/4−ρ, 3/2−7⁢ρ, 1/2−2⁢ρ, 2⁢ρ−1/45/2-12\rho,\,1/4-\rho,\,3/2-7\rho,\,1/2-2\rho,\,2\rho-1/4 11,10,11,9,11,1111,10,11,9,11,11 21/4−10⁢ρ21/4-10\rho 841/69696841/69696
(5/24,11/52)(5/24,11/52) 3/2−7⁢ρ, 1/4−ρ, 1/2−2⁢ρ, 2⁢ρ−1/4, 5/2−11⁢ρ3/2-7\rho,\,1/4-\rho,\,1/2-2\rho,\,2\rho-1/4,\,5/2-11\rho 10,9,9,11,11,1010,9,9,11,11,10 4−4⁢ρ4-4\rho 493/48672493/48672
(11/52,3/14)(11/52,3/14) 3/2−7⁢ρ, 1/4−ρ, 1/2−2⁢ρ, 5/2−11⁢ρ, 2⁢ρ−1/43/2-7\rho,\,1/4-\rho,\,1/2-2\rho,\,5/2-11\rho,\,2\rho-1/4 10,9,9,11,10,1010,9,9,11,10,10 4−4⁢ρ4-4\rho 573/66248573/66248
(3/14,7/32)(3/14,7/32) 1/4−ρ, 1/2−2⁢ρ, 5/2−11⁢ρ, 2⁢ρ−1/4, 3/2−6⁢ρ1/4-\rho,\,1/2-2\rho,\,5/2-11\rho,\,2\rho-1/4,\,3/2-6\rho 9,9,11,10,10,109,9,11,10,10,10 5/2+3⁢ρ5/2+3\rho 1411/1003521411/100352
(7/32,2/9)(7/32,2/9) 1/4−ρ, 1/2−2⁢ρ, 5/2−11⁢ρ, 3/2−6⁢ρ, 2⁢ρ−1/41/4-\rho,\,1/2-2\rho,\,5/2-11\rho,\,3/2-6\rho,\,2\rho-1/4 9,9,11,10,9,109,9,11,10,9,10 17/4−5⁢ρ17/4-5\rho 1813/1658881813/165888
(2/9,9/40)(2/9,9/40) 1/4−ρ, 5/2−11⁢ρ, 1/2−2⁢ρ, 3/2−6⁢ρ, 2⁢ρ−1/41/4-\rho,\,5/2-11\rho,\,1/2-2\rho,\,3/2-6\rho,\,2\rho-1/4 9,9,9,10,9,109,9,9,10,9,10 9/4+4⁢ρ9/4+4\rho 283/32400283/32400
(9/40,5/22)(9/40,5/22) 5/2−11⁢ρ, 1/4−ρ, 1/2−2⁢ρ, 3/2−6⁢ρ, 2⁢ρ−1/45/2-11\rho,\,1/4-\rho,\,1/2-2\rho,\,3/2-6\rho,\,2\rho-1/4 9,9,9,10,9,109,9,9,10,9,10 9/4+4⁢ρ9/4+4\rho 347/48400347/48400
(5/22,11/48)(5/22,11/48) 1/4−ρ, 1/2−2⁢ρ, 3/2−6⁢ρ, 2⁢ρ−1/4, 5/2−10⁢ρ1/4-\rho,\,1/2-2\rho,\,3/2-6\rho,\,2\rho-1/4,\,5/2-10\rho 9,9,10,9,10,99,9,10,9,10,9 19/4−7⁢ρ19/4-7\rho 3329/5575683329/557568
(11/48,1/4)(11/48,1/4) 1/4−ρ, 1/2−2⁢ρ, 3/2−6⁢ρ, 5/2−10⁢ρ, 2⁢ρ−1/41/4-\rho,\,1/2-2\rho,\,3/2-6\rho,\,5/2-10\rho,\,2\rho-1/4 9,9,10,9,9,99,9,10,9,9,9 2+5⁢ρ2+5\rho 307/4608307/4608
(1/4,11/40)(1/4,11/40) ρ−1/4, 5/2−9⁢ρ, 3/2−5⁢ρ, 1/2−ρ, 1/4\rho-1/4,\,5/2-9\rho,\,3/2-5\rho,\,1/2-\rho,\,1/4 8,9,8,7,8,98,9,8,7,8,9 9/2−5⁢ρ9/2-5\rho 51/64051/640
(11/40,5/18)(11/40,5/18) 5/2−9⁢ρ,ρ−1/4, 3/2−5⁢ρ, 1/2−ρ, 1/45/2-9\rho,\,\rho-1/4,\,3/2-5\rho,\,1/2-\rho,\,1/4 8,8,8,7,8,98,8,8,7,8,9 7/4+5⁢ρ7/4+5\rho 451/51840451/51840
(5/18,9/32)(5/18,9/32) ρ−1/4, 3/2−5⁢ρ, 1/2−ρ, 1/4, 5/2−8⁢ρ\rho-1/4,\,3/2-5\rho,\,1/2-\rho,\,1/4,\,5/2-8\rho 8,8,7,8,9,88,8,7,8,9,8 17/4−4⁢ρ17/4-4\rho 451/41472451/41472
(9/32,2/7)(9/32,2/7) ρ−1/4, 3/2−5⁢ρ, 1/2−ρ, 5/2−8⁢ρ, 1/4\rho-1/4,\,3/2-5\rho,\,1/2-\rho,\,5/2-8\rho,\,1/4 8,8,7,8,7,88,8,7,8,7,8 17/4−4⁢ρ17/4-4\rho 349/25088349/25088
(2/7,7/24)(2/7,7/24) ρ−1/4, 3/2−5⁢ρ, 5/2−8⁢ρ, 1/2−ρ, 1/4\rho-1/4,\,3/2-5\rho,\,5/2-8\rho,\,1/2-\rho,\,1/4 8,8,7,7,7,88,8,7,7,7,8 9/4+3⁢ρ9/4+3\rho 349/18816349/18816
(7/24,3/10)(7/24,3/10) 3/2−5⁢ρ,ρ−1/4, 5/2−8⁢ρ, 1/2−ρ, 1/43/2-5\rho,\,\rho-1/4,\,5/2-8\rho,\,1/2-\rho,\,1/4 8,6,7,7,7,88,6,7,7,7,8 4−3⁢ρ4-3\rho 83/320083/3200
(3/10,11/36)(3/10,11/36) ρ−1/4, 5/2−8⁢ρ, 1/2−ρ, 1/4, 3/2−4⁢ρ\rho-1/4,\,5/2-8\rho,\,1/2-\rho,\,1/4,\,3/2-4\rho 6,7,7,7,8,76,7,7,7,8,7 5/2+2⁢ρ5/2+2\rho 559/32400559/32400
(11/36,5/16)(11/36,5/16) 5/2−8⁢ρ,ρ−1/4, 1/2−ρ, 1/4, 3/2−4⁢ρ5/2-8\rho,\,\rho-1/4,\,1/2-\rho,\,1/4,\,3/2-4\rho 6,6,7,7,8,76,6,7,7,8,7 5/2+2⁢ρ5/2+2\rho 449/20736449/20736
(5/16,9/28)(5/16,9/28) ρ−1/4, 1/2−ρ, 3/2−4⁢ρ, 1/4, 5/2−7⁢ρ\rho-1/4,\,1/2-\rho,\,3/2-4\rho,\,1/4,\,5/2-7\rho 6,7,7,7,7,76,7,7,7,7,7 5/4+6⁢ρ5/4+6\rho 353/12544353/12544
(9/28,1/3)(9/28,1/3) ρ−1/4, 1/2−ρ, 3/2−4⁢ρ, 5/2−7⁢ρ, 1/4\rho-1/4,\,1/2-\rho,\,3/2-4\rho,\,5/2-7\rho,\,1/4 6,7,7,7,6,76,7,7,7,6,7 7/2−ρ7/2-\rho 533/14112533/14112
(1/3,11/32)(1/3,11/32) ρ−1/4, 5/2−7⁢ρ, 3/2−4⁢ρ, 1/2−ρ, 1/4\rho-1/4,\,5/2-7\rho,\,3/2-4\rho,\,1/2-\rho,\,1/4 6,7,6,5,6,76,7,6,5,6,7 9/2−4⁢ρ9/2-4\rho 151/4608151/4608
(11/32,7/20)(11/32,7/20) 5/2−7⁢ρ,ρ−1/4, 3/2−4⁢ρ, 1/2−ρ, 1/45/2-7\rho,\,\rho-1/4,\,3/2-4\rho,\,1/2-\rho,\,1/4 6,6,6,5,6,76,6,6,5,6,7 7/4+4⁢ρ7/4+4\rho 251/12800251/12800
(7/20,5/14)(7/20,5/14) 5/2−7⁢ρ, 3/2−4⁢ρ,ρ−1/4, 1/2−ρ, 1/45/2-7\rho,\,3/2-4\rho,\,\rho-1/4,\,1/2-\rho,\,1/4 6,6,4,5,6,76,6,4,5,6,7 7/2−ρ7/2-\rho 881/39200881/39200
(5/14,3/8)(5/14,3/8) 3/2−4⁢ρ,ρ−1/4, 1/2−ρ, 1/4, 5/2−6⁢ρ3/2-4\rho,\,\rho-1/4,\,1/2-\rho,\,1/4,\,5/2-6\rho 6,4,5,6,7,66,4,5,6,7,6 6−8⁢ρ6-8\rho 43/78443/784
(3/8,11/28)(3/8,11/28) 1/2−ρ,ρ−1/4, 5/2−6⁢ρ, 1/4, 3/2−3⁢ρ1/2-\rho,\,\rho-1/4,\,5/2-6\rho,\,1/4,\,3/2-3\rho 4,6,6,5,6,54,6,6,5,6,5 15/4−2⁢ρ15/4-2\rho 167/3136167/3136
(11/28,2/5)(11/28,2/5) 1/2−ρ, 5/2−6⁢ρ,ρ−1/4, 1/4, 3/2−3⁢ρ1/2-\rho,\,5/2-6\rho,\,\rho-1/4,\,1/4,\,3/2-3\rho 4,6,5,5,6,54,6,5,5,6,5 15/4−2⁢ρ15/4-2\rho 207/9800207/9800
(2/5,5/12)(2/5,5/12) 5/2−6⁢ρ, 1/2−ρ,ρ−1/4, 1/4, 3/2−3⁢ρ5/2-6\rho,\,1/2-\rho,\,\rho-1/4,\,1/4,\,3/2-3\rho 4,4,5,5,6,54,4,5,5,6,5 7/4+3⁢ρ7/4+3\rho 119/2400119/2400
(5/12,7/16)(5/12,7/16) 1/2−ρ,ρ−1/4, 3/2−3⁢ρ, 1/4, 5/2−5⁢ρ1/2-\rho,\,\rho-1/4,\,3/2-3\rho,\,1/4,\,5/2-5\rho 4,5,5,5,5,54,5,5,5,5,5 1/2+6⁢ρ1/2+6\rho 49/76849/768
(7/16,9/20)(7/16,9/20) 1/2−ρ, 3/2−3⁢ρ,ρ−1/4, 1/4, 5/2−5⁢ρ1/2-\rho,\,3/2-3\rho,\,\rho-1/4,\,1/4,\,5/2-5\rho 4,5,4,5,5,54,5,4,5,5,5 9/4+2⁢ρ9/4+2\rho 251/6400251/6400
(9/20,11/24)(9/20,11/24) 1/2−ρ, 3/2−3⁢ρ,ρ−1/4, 5/2−5⁢ρ, 1/41/2-\rho,\,3/2-3\rho,\,\rho-1/4,\,5/2-5\rho,\,1/4 4,5,4,5,4,54,5,4,5,4,5 9/2−3⁢ρ9/2-3\rho 251/9600251/9600
(11/24,1/2)(11/24,1/2) 1/2−ρ, 3/2−3⁢ρ, 5/2−5⁢ρ,ρ−1/4, 1/41/2-\rho,\,3/2-3\rho,\,5/2-5\rho,\,\rho-1/4,\,1/4 4,5,4,4,4,54,5,4,4,4,5 7/4+3⁢ρ7/4+3\rho 17/12817/128
(1/2,11/20)(1/2,11/20) 1/4,ρ−1/4, 5/2−4⁢ρ, 3/2−2⁢ρ, 1/21/4,\,\rho-1/4,\,5/2-4\rho,\,3/2-2\rho,\,1/2 4,4,5,4,3,44,4,5,4,3,4 19/4−3⁢ρ19/4-3\rho 127/800127/800
(11/20,9/16)(11/20,9/16) 1/4, 5/2−4⁢ρ,ρ−1/4, 3/2−2⁢ρ, 1/21/4,\,5/2-4\rho,\,\rho-1/4,\,3/2-2\rho,\,1/2 4,4,4,4,3,44,4,4,4,3,4 2+2⁢ρ2+2\rho 249/6400249/6400
(9/16,7/12)(9/16,7/12) 5/2−4⁢ρ, 1/4,ρ−1/4, 3/2−2⁢ρ, 1/25/2-4\rho,\,1/4,\,\rho-1/4,\,3/2-2\rho,\,1/2 4,3,4,4,3,44,3,4,4,3,4 17/4−2⁢ρ17/4-2\rho 149/2304149/2304
(7/12,5/8)(7/12,5/8) 5/2−4⁢ρ, 1/4, 3/2−2⁢ρ,ρ−1/4, 1/25/2-4\rho,\,1/4,\,3/2-2\rho,\,\rho-1/4,\,1/2 4,3,4,2,3,44,3,4,2,3,4 6−5⁢ρ6-5\rho 143/1152143/1152
(5/8,2/3)(5/8,2/3) 3/2−2⁢ρ, 1/4,ρ−1/4, 1/2, 5/2−3⁢ρ3/2-2\rho,\,1/4,\,\rho-1/4,\,1/2,\,5/2-3\rho 3,2,2,3,4,33,2,2,3,4,3 19/4−3⁢ρ19/4-3\rho 15/12815/128
(2/3,11/16)(2/3,11/16) 3/2−2⁢ρ, 1/4,ρ−1/4, 5/2−3⁢ρ, 1/23/2-2\rho,\,1/4,\,\rho-1/4,\,5/2-3\rho,\,1/2 3,2,2,3,3,33,2,2,3,3,3 11/411/4 11/19211/192
(11/16,3/4)(11/16,3/4) 3/2−2⁢ρ, 1/4, 5/2−3⁢ρ,ρ−1/4, 1/23/2-2\rho,\,1/4,\,5/2-3\rho,\,\rho-1/4,\,1/2 3,2,2,2,3,33,2,2,2,3,3 11/411/4 11/6411/64
(3/4,5/6)(3/4,5/6) 5/2−3⁢ρ, 1/4, 1/2,ρ−1/4, 3/2−ρ5/2-3\rho,\,1/4,\,1/2,\,\rho-1/4,\,3/2-\rho 2,2,2,3,3,32,2,2,3,3,3 1/2+3⁢ρ1/2+3\rho 23/9623/96
(5/6,7/8)(5/6,7/8) 1/4, 1/2,ρ−1/4, 3/2−ρ, 5/2−2⁢ρ1/4,\,1/2,\,\rho-1/4,\,3/2-\rho,\,5/2-2\rho 2,2,3,3,3,22,2,3,3,3,2 33 1/81/8
(7/8,11/12)(7/8,11/12) 1/4, 1/2, 3/2−ρ,ρ−1/4, 5/2−2⁢ρ1/4,\,1/2,\,3/2-\rho,\,\rho-1/4,\,5/2-2\rho 2,2,3,2,3,22,2,3,2,3,2 19/4−2⁢ρ19/4-2\rho 71/57671/576
(11/12,1)(11/12,1) 1/4, 1/2, 3/2−ρ, 5/2−2⁢ρ,ρ−1/41/4,\,1/2,\,3/2-\rho,\,5/2-2\rho,\,\rho-1/4 2,2,3,2,1,22,2,3,2,1,2 19/4−2⁢ρ19/4-2\rho 17/7217/72
(1,9/8)(1,9/8) 1/4, 5/2−2⁢ρ, 3/2−ρ, 1/2,ρ−1/41/4,\,5/2-2\rho,\,3/2-\rho,\,1/2,\,\rho-1/4 2,2,2,0,1,22,2,2,0,1,2 15/4−ρ15/4-\rho 43/12843/128
(9/8,5/4)(9/8,5/4) 5/2−2⁢ρ, 1/4, 3/2−ρ, 1/2,ρ−1/45/2-2\rho,\,1/4,\,3/2-\rho,\,1/2,\,\rho-1/4 2,1,2,0,1,22,1,2,0,1,2 6−3⁢ρ6-3\rho 39/12839/128
∫1/55/4𝔠⁢(ρ)⁢𝑑ρ\displaystyle\int_{1/5}^{5/4}\mathfrak{c}(\rho)\,d\rho =1121513/369600,\displaystyle=1121513/369600, ∫5/411/4𝔠out⁢(ρ)⁢𝑑ρ\displaystyle\int_{5/4}^{11/4}\mathfrak{c}_{\rm out}(\rho)\,d\rho =45/32,\displaystyle=45/32,
134⋅15\displaystyle\frac{13}{4}\cdot\frac{1}{5} =13/20,\displaystyle=13/20, sum =1881503/369600,\displaystyle=1881503/369600,
coefficient of ⁢log⁡2\displaystyle\text{coefficient of }\log 2 =53/48,\displaystyle=53/48, contribution of (3.5) =−5/48.\displaystyle=-5/48.

Appendix B Numerical values

The bounds in Table B.1 are rounded upward. The rational enclosures are

(B.1) q−:=−4.911220867982619\displaystyle q_{-}:=-4.911220867982619 <𝒬⁢(ν)<−4.911220867982618,\displaystyle<\mathcal{Q}(\nu)<-4.911220867982618,
(B.2) 4.365058361177462\displaystyle 4.365058361177462 <α<4.365058361177465,\displaystyle<\alpha<4.365058361177465,
3.730290307027085\displaystyle 3.730290307027085 <ϰ<3.730290307027087;\displaystyle<\varkappa<3.730290307027087;

the lower bound for α\alpha holds with q−q_{-} in place of 𝒬⁢(ν)\mathcal{Q}(\nu).

Table B.1. Upper bounds for Fx,Fv,FzF_{x},F_{v},F_{z} in Lemma 4.2.
support subintervals depth upper bound KK
x∈[0,1]x\in[0,1] 1164 15 −4.6260004267-4.6260004267 −4.626-4.626
v∈[0,1]v\in[0,1] 1272 16 −4.4800027327-4.4800027327 −4.480-4.480
z∈[−2,−3/5]z\in[-2,-3/5] 658 14 −0.0900423691-0.0900423691 −0.090-0.090
z∈[3/5,1]z\in[3/5,1] 4 2 −0.3671210997-0.3671210997 −0.090-0.090
|z|=3/5,Im⁡z≥0|z|=3/5,\ \operatorname{Im}z\geq 0 693 11 −0.0900220839-0.0900220839 −0.090-0.090
Table B.2. Constants in Lemmas 4.2–4.4.
quantity value
intervals on xx, vv, zz; Poisson kernels (50,52,62,157)(50,52,62,157)
(Kx,Kv,Kz)(K_{x},K_{v},K_{z}) (−4.626,−4.480,−0.090)(-4.626,-4.480,-0.090)
δ\delta 1/2501/250
sup|q|≤δFx⁢(q)\sup_{|q|\leq\delta}F_{x}(q) <−6.2823007011<-6.2823007011
sup|q|≤δFv⁢(q)\sup_{|q|\leq\delta}F_{v}(q) <−4.4935082437<-4.4935082437
(Hx,Hv,Hz,real,Hz,circle)(H_{x},H_{v},H_{z,\mathrm{real}},H_{z,\mathrm{circle}}) (131,153,43,24)(131,153,43,24)
subintervals in Table B.1 37913791
precision (decimal digits) 9090
series terms 420420
series remainder per pair <8.78×10−130<8.78\times 10^{-130}
HzH_{z} 6767
(Cx,Cv,Cz)(C_{x},C_{v},C_{z}) (2096,2984,2296)(2096,2984,2296)
Cx+54⁢Cv+CzC_{x}+\tfrac{5}{4}C_{v}+C_{z} 81228122
(LD,Lz)(L_{D},L_{z}) (101/100,167/100)(101/100,167/100)
A0A_{0} 1107/801107/80
1716⋅8⁢ϵ/δ2+A0\tfrac{17}{16}\cdot 8\epsilon/\delta^{2}+A_{0} at m=64m=64 <200<200
3⁢(Cx+54⁢Cv+Cz)+2003(C_{x}+\tfrac{5}{4}C_{v}+C_{z})+200 2456624566
C∗C_{*} <810<810
endpoint-disk precision (bits) 100100

Appendix C Printed interval data

Tables C.1–C.3 give the numerical enclosures used in Lemma 4.2. Every displayed interval is outward rounded. The static row checks compare each upper endpoint with the printed bound.

Table C.1. The pair decomposition and its enclosures.

Here Iz⁢rI_{zr} is the energy of the real part of νz\nu_{z}, Iz⁢pI_{zp} the cross term of its real and Poisson parts, and Ip⁢pI_{pp} the Poisson self-energy. Iv⁢z=I⁢(νv,νz)I_{vz}=I(\nu_{v},\nu_{z}) includes both the real and Poisson parts of νz\nu_{z}. term outward enclosure Ix⁢xI_{xx} [-1.546417482857900440,-1.546417482857900439] Iv⁢vI_{vv} [-2.239637943680446541,-2.239637943680446540] Iz⁢rI_{zr} [-0.162577112801730912,-0.162577112801730911] Iz⁢pI_{zp} [0.081132573652198143,0.081132573652198144] Ip⁢pI_{pp} [-0.272227665742282938,-0.272227665742282937] Iv⁢zI_{vz} [0.206849817575968853,0.206849817575968854] Dx⁢vD_{xv} [1.059475627780623120,1.059475627780623121] 𝒬⁢(ν)=Ix⁢x+Iv⁢v+Iz⁢r+2⁢Iz⁢p+Ip⁢p+Iv⁢z−Dx⁢v\mathcal{Q}(\nu)=I_{xx}+I_{vv}+I_{zr}+2I_{zp}+I_{pp}+I_{vz}-D_{xv} [-4.911220867982618809,-4.911220867982618808]

For i∈{x,v}i\in\{x,v\} put Ai⁢(δ)=∫log⁡(|t|+δ)⁢𝑑νi⁢(t)A_{i}(\delta)=\int\log(|t|+\delta)\,d\nu_{i}(t), and put Az,real⁢(δ)=∫log⁡(|t|+δ)⁢𝑑νz,real⁢(t)A_{z,\mathrm{real}}(\delta)=\int\log(|t|+\delta)\,d\nu_{z,\mathrm{real}}(t). Put

Dx,low=log⁡(2−δ),Dv,low=54⁢log⁡(2−δ).D_{x,\mathrm{low}}=\log(2-\delta),\qquad D_{v,\mathrm{low}}=\tfrac{5}{4}\log(2% -\delta).

Then Di⁢(q)≥Di,lowD_{i}(q)\geq D_{i,\mathrm{low}} for i∈{x,v}i\in\{x,v\} and |q|≤δ|q|\leq\delta. For |q|≤δ|q|\leq\delta the displayed endpoint calculation reads

Fx⁢(q)≤Bx:=2⁢Ax⁢(δ)−Dv,low+34⁢log⁡δ,Fv⁢(q)≤Bv:=2⁢Av⁢(δ)+Az,real⁢(δ)+WP−Dx,low+14⁢log⁡δ.F_{x}(q)\leq B_{x}:=2A_{x}(\delta)-D_{v,\mathrm{low}}+\tfrac{3}{4}\log\delta,% \qquad F_{v}(q)\leq B_{v}:=2A_{v}(\delta)+A_{z,\mathrm{real}}(\delta)+W_{P}-D_% {x,\mathrm{low}}+\tfrac{1}{4}\log\delta.

The denominator expansion retains 420420 terms. Since 0≤x⁢v/2≤1/20\leq xv/2\leq 1/2 and the total xx–vv mass is 1⋅(5/4)1\cdot(5/4), the integrated remainder is bounded by

Rx⁢v≤54⁢2⁢(1/2)421421=52⋅421⋅2421<1.10×10−129.R_{xv}\leq\frac{5}{4}\,\frac{2(1/2)^{421}}{421}=\frac{5}{2\cdot 421\cdot 2^{42% 1}}<1.10\times 10^{-129}.
Table C.2. Endpoint majorants and disk bounds for δ=1/250\delta=1/250.
quantity outward enclosure
Ax⁢(δ)A_{x}(\delta) [-0.638636770206138024,-0.638636770206138023]
Av⁢(δ)A_{v}(\delta) [-1.118528621477847036,-1.118528621477847035]
Az,real⁢(δ)A_{z,\mathrm{real}}(\delta) [0.148575311504214922,0.148575311504214923]
WP:=MP⁢log⁡(151/250)W_{P}:=M_{P}\log(151/250) [-0.333515904969252,-0.333515904969250]
MPM_{P} (sum of the 157 Poisson weights) 26460009509/40000000000
Dx,lowD_{x,\mathrm{low}} [0.691145177889270,0.691145177889275]
Dv,lowD_{v,\mathrm{low}} [0.863931472361588,0.863931472361593]
BxB_{x} [-6.28230070117056,-6.28230070117054]
BvB_{v} [-4.49350824377558,-4.49350824377554]
Table C.3. Rational density checks from the printed reference tables.
support exact maximum used integer upper bound
xx 1772217287/135816471772217287/13581647 131131
vv 2076994349/135816472076994349/13581647 153153
zz real 4076507975/950715294076507975/95071529 4343
zz circle ∑Hc<24\sum H_{c}<24 by the integer square-root and π>333/106\pi>333/106 checks 2424

The first three bounds in Table C.3 are rational comparisons of rows in Table D.1. The fourth is the corresponding rational sum over Table D.2.

Appendix D Reference measure data

The two tables specify the reference measure ν\nu. All entries are exact terminating decimals.

Table D.1 lists the intervals of ν\nu on each real support in increasing order of the index, and each line holds two consecutive entries (odd index left, even index right). Entry (a,w)(a,w) is the interval [a,b][a,b] carrying w⁢d⁢t/(b−a)w\,dt/(b-a), where bb is the aa of the entry with the next index in the same block. The last bb of each block is given in its heading. Outside these intervals ν\nu has no mass. For example, νx\nu_{x} vanishes on [0,0.02653493525245)[0,0.02653493525245).

Table D.2 lists the circular part of νz\nu_{z}. Every row is a measure on the circle |q|=r=3/5|q|=r=3/5, whose upper half is the support z∘z_{\circ} of Appendix E. It is the Poisson measure of the proof of Lemma 4.2 with center c=Re⁡c+i⁢Im⁡cc=\operatorname{Re}c+i\operatorname{Im}c, |c|<3/5|c|<3/5, and total mass ww. Its logarithmic potential is (4.3). The xx- and vv-parts of ν\nu have no Poisson rows.

Table D.1. Interval data for the reference measure: left endpoint aa and mass ww of each interval.
index aa ww index aa ww
xx-intervals on [0,1][0,1]: rows 1–50, last right endpoint 11.
1 0.02653493525245 0.000000000011 2 0.03806023374436 0.000000000056
3 0.05156362923366 0.000000000066 4 0.06698729810778 0.009047636046
5 0.08426519384873 0.016164978747 6 0.10332332985438 0.019034345994
7 0.12408009626051 0.021065372172 8 0.14644660940673 0.022447569912
9 0.17032709244997 0.02345382024 10 0.19561928549564 0.024200308115
11 0.2222148834902 0.02476625456 12 0.25 0.025199538621
13 0.2788556548905 0.02553284359 14 0.30865828381746 0.025790757407
15 0.33928026734842 0.02598943694 16 0.37059047744874 0.026140976998
17 0.40245483899194 0.026255972454 18 0.43473690388997 0.026341379215
19 0.46729843538493 0.026402361973 20 0.5 0.026444293432
21 0.53270156461507 0.026470882322 22 0.56526309611003 0.026484444361
23 0.59754516100806 0.026487798034 24 0.62940952255126 0.026483565514
25 0.66071973265158 0.026473037339 26 0.69134171618254 0.026457122002
27 0.7211443451095 0.026437395821 28 0.75 0.026415300265
29 0.7777851165098 0.026391297706 30 0.80438071450436 0.026365737946
31 0.82967290755003 0.026339451805 32 0.85355339059327 0.026313133347
33 0.87591990373949 0.026286966515 34 0.89667667014562 0.026260716757
35 0.91573480615127 0.02623413097 36 0.93301270189222 0.026207204488
37 0.94843637076634 0.026178707612 38 0.96193976625564 0.026146051892
39 0.97346506474755 0.026100427024 40 0.98296291314453 0.026045558902
41 0.99039264020162 0.025800471869 42 0.99572243068691 0.02600510817
43 0.9989294616193 0.00101387447 44 0.9989702065603 0.001749595019
45 0.99908623833539 0.002577120956 46 0.99925989215865 0.003342459833
47 0.99946473080965 0.004008718139 48 0.99966956946065 0.004554823645
49 0.99984322328391 0.004774398867 50 0.999959255059 0.005316651861
vv-intervals on [0,1][0,1]: rows 51–102, last right endpoint 11.
51 0.00960735979838 0.013858434309 52 0.01703708685547 0.020659018632
53 0.02653493525245 0.022756374414 54 0.03806023374436 0.024040451088
55 0.05156362923366 0.024796480607 56 0.06698729810778 0.025298366124
57 0.08426519384873 0.025649764501 58 0.10332332985438 0.025910099392
59 0.12408009626051 0.026113452847 60 0.14644660940673 0.02628122952
61 0.17032709244997 0.026427269003 62 0.19561928549564 0.026561515238
63 0.2222148834902 0.026691316341 64 0.25 0.026821861266
65 0.2788556548905 0.026957166859 66 0.30865828381746 0.027100366588
67 0.33928026734842 0.027253308259 68 0.37059047744874 0.027416681441
69 0.40245483899194 0.027590343833 70 0.43473690388997 0.027773157352
71 0.46729843538493 0.027962829423 72 0.5 0.028155382615
73 0.53270156461507 0.028348951325 74 0.56526309611003 0.028512049639
75 0.59754516100806 0.028704024914 76 0.62940952255126 0.02895179599
77 0.66071973265158 0.029147228949 78 0.69134171618254 0.029330189043
79 0.7211443451095 0.029497893354 80 0.75 0.029651492159
81 0.7777851165098 0.029791323098 82 0.80438071450436 0.029918006281
83 0.82967290755003 0.030031872164 84 0.85355339059327 0.030133175201
85 0.87591990373949 0.030222548301 86 0.89667667014562 0.030300285976
87 0.91573480615127 0.030366116749 88 0.93301270189222 0.030419889284
89 0.94843637076634 0.030460392026 90 0.96193976625564 0.030485017656
91 0.97346506474755 0.030483706899 92 0.98296291314453 0.030461131797
93 0.99039264020162 0.030205136068 94 0.99572243068691 0.030465558928
95 0.9989294616193 0.001188270853 96 0.9989702065603 0.002050223725
97 0.99908623833539 0.003019909346 98 0.99925989215865 0.003916871095
99 0.99946473080965 0.004697775042 100 0.99966956946065 0.005337916652
101 0.99984322328391 0.005595394787 102 0.999959255059 0.006230983047
zz-intervals on [−2,−3/5][-2,-3/5]: rows 103–164, last right endpoint −3/5-3/5.
103 -2 0.002445904785 104 -1.9999429570826 0.002196289941
105 -1.99978051259747 0.002095070288 106 -1.99953739724491 0.001843711398
107 -1.99925062313351 0.001537127441 108 -1.99896384902211 0.001185108178
109 -1.99872073366955 0.000804636229 110 -1.99855828918442 0.000466572663
111 -1.99850124626702 0.011949094627 112 -1.99401140296167 0.011829831253
113 -1.98654969628226 0.011904636078 114 -1.97614807840235 0.011879377704
115 -1.96285109064657 0.011837079332 116 -1.9467156727579 0.01177589427
117 -1.92781091907288 0.011699589994 118 -1.90621778264911 0.011609193707
119 -1.88202872861178 0.011505113623 120 -1.85534733820386 0.011387043787
121 -1.82628786523528 0.011255100131 122 -1.79497474683058 0.011109018301
123 -1.76154207057005 0.010947966233 124 -1.7261330003061 0.01077167805
125 -1.68889916311372 0.010579949336 126 -1.65 0.010372009648
127 -1.6096020831533 0.010147136288 128 -1.56787840265556 0.009904930406
129 -1.52500762571221 0.009644863536 130 -1.48117333157176 0.009366270662
131 -1.43656322541129 0.009068614306 132 -1.39136833455404 0.00875137721
133 -1.3457821904611 0.008414043831 134 -1.3 0.008056385332
135 -1.2542178095389 0.007678505629 136 -1.20863166544596 0.007280632319
137 -1.16343677458871 0.006863145229 138 -1.11882666842824 0.006426779891
139 -1.07499237428779 0.005972788018 140 -1.03212159734444 0.005503141876
141 -0.9903979168467 0.00502059841 142 -0.95 0.004528595244
143 -0.91110083688628 0.004031298216 144 -0.8738669996939 0.00353382249
145 -0.83845792942995 0.003042591462 146 -0.80502525316942 0.002564906577
147 -0.77371213476472 0.002108555487 148 -0.74465266179614 0.001681749548
149 -0.71797127138822 0.001293146877 150 -0.69378221735089 0.000950929581
151 -0.67218908092712 0.000661669911 152 -0.6532843272421 0.000429830492
153 -0.63714890935343 0.000257484313 154 -0.62385192159765 0.000143377053
155 -0.61345030371774 0.000084543242 156 -0.60598859703833 0.000061123507
157 -0.60149875373298 0.000001405548 158 -0.60144171081558 0.000003128248
159 -0.60127926633045 0.00000456608 160 -0.60103615097789 0.000005613821
161 -0.60074937686649 0.00000618074 162 -0.60046260275509 0.000006409107
163 -0.60021948740253 0.00000617203 164 -0.6000570429174 0.000006452761
Table D.2. Poisson data for the circular part of the reference measure: center cc and mass ww of each Poisson measure.
index Re⁡c\operatorname{Re}c Im⁡c\operatorname{Im}c ww
1 0.4913985408708 0.32834200771458 0.000000000005
2 0.48039239059364 0.34424431884887 0.000000000009
3 0.46887182411212 0.35977800454415 0.00000000001
4 0.45684917793738 0.37492643094071 0.000000000011
5 0.44433732622008 0.38967337672414 0.000638704082
6 0.43134966696442 0.40400305049552 0.002774568368
7 0.41790010768125 0.41790010768125 0.002648678352
8 0.40400305049552 0.43134966696442 0.003134952159
9 0.38967337672414 0.44433732622008 0.003371109225
10 0.37492643094071 0.45684917793738 0.00362052744
11 0.35977800454415 0.46887182411212 0.003823085588
12 0.34424431884887 0.48039239059364 0.004004693717
13 0.32834200771458 0.4913985408708 0.004164215022
14 0.31208809973437 0.50187848928221 0.004306593404
15 0.2955 0.5118210136366 0.004434028868
16 0.27859547146416 0.52121546722988 0.004548494496
17 0.26139261591943 0.53005179024582 0.004651574501
18 0.2439098546144 0.53832052052841 0.004744445721
19 0.22616590852777 0.54601280371417 0.004828370754
20 0.20817977832145 0.55312040271358 0.004904191893
21 0.18997072399417 0.55963570653161 0.004972725611
22 0.17155824425739 0.56555173841774 0.005034518792
23 0.15296205565559 0.57086216333684 0.005090055986
24 0.13420207145331 0.5755612947529 0.005139889582
25 0.11529838031153 0.57964410071831 0.005184386706
26 0.0962712247762 0.58310620926216 0.005223998617
27 0.07714097960205 0.58594391307192 0.005258940541
28 0.05792812993477 0.58815417346327 0.005289495695
29 0.03865324937501 0.58973462363401 0.005315860418
30 0.01933697794767 0.59068357119853 0.005338252913
31 0 0.591 0.005356869923
32 -0.01933697794767 0.59068357119853 0.005371839655
33 -0.03865324937501 0.58973462363401 0.00538329912
34 -0.05792812993477 0.58815417346327 0.005391294222
35 -0.07714097960205 0.58594391307192 0.005395976833
36 -0.0962712247762 0.58310620926216 0.005397361405
37 -0.11529838031153 0.57964410071831 0.005395592218
38 -0.13420207145331 0.5755612947529 0.005390634646
39 -0.15296205565559 0.57086216333684 0.005382588982
40 -0.17155824425739 0.56555173841774 0.005371424089
41 -0.18997072399417 0.55963570653161 0.005357167758
42 -0.20817977832145 0.55312040271358 0.005339817089
43 -0.22616590852777 0.54601280371417 0.005319360695
44 -0.2439098546144 0.53832052052841 0.005295824221
45 -0.26139261591943 0.53005179024582 0.005269193571
46 -0.27859547146416 0.52121546722988 0.005239457786
47 -0.2955 0.5118210136366 0.005206556882
48 -0.31208809973437 0.50187848928221 0.005170430831
49 -0.32834200771458 0.4913985408708 0.005131018794
50 -0.34424431884887 0.48039239059364 0.005088322391
51 -0.35977800454415 0.46887182411212 0.00504230892
52 -0.37492643094071 0.45684917793738 0.004992975694
53 -0.38967337672414 0.44433732622008 0.004940222961
54 -0.40400305049552 0.43134966696442 0.00488399491
55 -0.41790010768125 0.41790010768125 0.004824218204
56 -0.43134966696442 0.40400305049552 0.004760895309
57 -0.44433732622008 0.38967337672414 0.004694049607
58 -0.45684917793738 0.37492643094071 0.004623666088
59 -0.46887182411212 0.35977800454415 0.004549682725
60 -0.48039239059364 0.34424431884887 0.004472013635
61 -0.4913985408708 0.32834200771459 0.004390610978
62 -0.50187848928221 0.31208809973437 0.00430554631
63 -0.5118210136366 0.2955 0.004216885068
64 -0.52121546722988 0.27859547146416 0.004124710654
65 -0.53005179024582 0.26139261591943 0.004028997465
66 -0.53832052052841 0.2439098546144 0.003929757184
67 -0.54601280371417 0.22616590852777 0.003827028338
68 -0.55312040271358 0.20817977832145 0.003720966921
69 -0.55963570653161 0.18997072399417 0.003611737022
70 -0.56555173841774 0.17155824425739 0.00349933102
71 -0.57086216333684 0.15296205565559 0.003383865867
72 -0.5755612947529 0.13420207145331 0.003265638507
73 -0.57964410071831 0.11529838031153 0.003143389946
74 -0.58310620926216 0.0962712247762 0.003020584609
75 -0.58594391307192 0.07714097960205 0.002890212871
76 -0.58815417346327 0.05792812993477 0.002763390356
77 -0.58973462363401 0.03865324937501 0.002615447991
78 -0.59068357119853 0.01933697794767 0.002506482686
79 -0.591 0 0.002042234839
80 -0.59068357119853 -0.01933697794767 0.002506482686
81 -0.58973462363401 -0.03865324937501 0.002615447991
82 -0.58815417346327 -0.05792812993477 0.002763390356
83 -0.58594391307192 -0.07714097960205 0.002890212871
84 -0.58310620926216 -0.0962712247762 0.003020584609
85 -0.57964410071831 -0.11529838031153 0.003143389946
86 -0.5755612947529 -0.13420207145331 0.003265638507
87 -0.57086216333684 -0.15296205565559 0.003383865867
88 -0.56555173841774 -0.17155824425739 0.00349933102
89 -0.55963570653161 -0.18997072399417 0.003611737022
90 -0.55312040271358 -0.20817977832145 0.003720966921
91 -0.54601280371417 -0.22616590852777 0.003827028338
92 -0.53832052052841 -0.2439098546144 0.003929757184
93 -0.53005179024582 -0.26139261591943 0.004028997465
94 -0.52121546722988 -0.27859547146416 0.004124710654
95 -0.5118210136366 -0.2955 0.004216885068
96 -0.50187848928221 -0.31208809973437 0.00430554631
97 -0.4913985408708 -0.32834200771459 0.004390610978
98 -0.48039239059364 -0.34424431884887 0.004472013635
99 -0.46887182411212 -0.35977800454415 0.004549682725
100 -0.45684917793738 -0.37492643094071 0.004623666088
101 -0.44433732622008 -0.38967337672414 0.004694049607
102 -0.43134966696442 -0.40400305049552 0.004760895309
103 -0.41790010768125 -0.41790010768125 0.004824218204
104 -0.40400305049552 -0.43134966696442 0.00488399491
105 -0.38967337672414 -0.44433732622008 0.004940222961
106 -0.37492643094071 -0.45684917793738 0.004992975694
107 -0.35977800454415 -0.46887182411212 0.00504230892
108 -0.34424431884887 -0.48039239059364 0.005088322391
109 -0.32834200771458 -0.4913985408708 0.005131018794
110 -0.31208809973437 -0.50187848928221 0.005170430831
111 -0.2955 -0.5118210136366 0.005206556882
112 -0.27859547146416 -0.52121546722988 0.005239457786
113 -0.26139261591943 -0.53005179024582 0.005269193571
114 -0.2439098546144 -0.53832052052841 0.005295824221
115 -0.22616590852777 -0.54601280371417 0.005319360695
116 -0.20817977832145 -0.55312040271358 0.005339817089
117 -0.18997072399417 -0.55963570653161 0.005357167758
118 -0.17155824425739 -0.56555173841774 0.005371424089
119 -0.15296205565559 -0.57086216333684 0.005382588982
120 -0.13420207145331 -0.5755612947529 0.005390634646
121 -0.11529838031153 -0.57964410071831 0.005395592218
122 -0.0962712247762 -0.58310620926216 0.005397361405
123 -0.07714097960205 -0.58594391307192 0.005395976833
124 -0.05792812993477 -0.58815417346327 0.005391294222
125 -0.03865324937501 -0.58973462363401 0.00538329912
126 -0.01933697794767 -0.59068357119853 0.005371839655
127 0 -0.591 0.005356869923
128 0.01933697794767 -0.59068357119853 0.005338252913
129 0.03865324937501 -0.58973462363401 0.005315860418
130 0.05792812993477 -0.58815417346327 0.005289495695
131 0.07714097960205 -0.58594391307192 0.005258940541
132 0.0962712247762 -0.58310620926216 0.005223998617
133 0.11529838031153 -0.57964410071831 0.005184386706
134 0.13420207145331 -0.5755612947529 0.005139889582
135 0.15296205565559 -0.57086216333684 0.005090055986
136 0.17155824425739 -0.56555173841774 0.005034518792
137 0.18997072399417 -0.55963570653161 0.004972725611
138 0.20817977832145 -0.55312040271358 0.004904191893
139 0.22616590852777 -0.54601280371417 0.004828370754
140 0.2439098546144 -0.53832052052841 0.004744445721
141 0.26139261591943 -0.53005179024582 0.004651574501
142 0.27859547146416 -0.52121546722988 0.004548494496
143 0.2955 -0.5118210136366 0.004434028868
144 0.31208809973437 -0.50187848928221 0.004306593404
145 0.32834200771458 -0.4913985408708 0.004164215022
146 0.34424431884887 -0.48039239059364 0.004004693717
147 0.35977800454415 -0.46887182411212 0.003823085588
148 0.37492643094071 -0.45684917793738 0.00362052744
149 0.38967337672414 -0.44433732622008 0.003371109225
150 0.40400305049552 -0.43134966696442 0.003134952159
151 0.41790010768125 -0.41790010768125 0.002648678352
152 0.43134966696442 -0.40400305049552 0.002774568368
153 0.44433732622008 -0.38967337672414 0.000638704082
154 0.45684917793738 -0.37492643094071 0.000000000011
155 0.46887182411212 -0.35977800454415 0.00000000001
156 0.48039239059364 -0.34424431884887 0.000000000009
157 0.4913985408708 -0.32834200771458 0.000000000005

The interval and Poisson weights sum to masses 11, 5/45/4, and 11 as stated in Lemma 4.2.

Appendix E Dyadic subdivisions

The supports are x,v∈[0,1]x,v\in[0,1], z−=[−2,−3/5]z_{-}=[-2,-3/5], z+=[3/5,1]z_{+}=[3/5,1] and z∘={|z|=3/5,Im⁡z≥0}z_{\circ}=\{|z|=3/5,\ \operatorname{Im}z\geq 0\} (by conjugation symmetry the upper half suffices). On each support ss, the leaves are ordered from left to right. Their depths ds,1,…,ds,nsd_{s,1},\ldots,d_{s,n_{s}} are listed in run-length form: an entry d:rd:r means that depth dd occurs on rr consecutive leaves. Set ts,0=0t_{s,0}=0 and ts,j=∑i=1j2−ds,it_{s,j}=\sum_{i=1}^{j}2^{-d_{s,i}}. The printed sequence satisfies ∑j=1ns2−ds,j=1\sum_{j=1}^{n_{s}}2^{-d_{s,j}}=1 for each support. On a real support [As,Bs][A_{s},B_{s}] the jjth interval is [As+(Bs−As)⁢ts,j−1,As+(Bs−As)⁢ts,j][A_{s}+(B_{s}-A_{s})t_{s,j-1},A_{s}+(B_{s}-A_{s})t_{s,j}]. On z∘z_{\circ} it is q=(3/5)⁢ei⁢π⁢tq=(3/5)e^{i\pi t} for t∈[ts,j−1,ts,j]t\in[t_{s,j-1},t_{s,j}]. On each leaf, the Arb upper bound for supFs\sup F_{s} is strictly below its cap in Table B.1, and that table keeps the worst leaf for each support. The full leafwise Arb upper bounds are in anc/verification/analytic/leaves.json.gz.

Table E.1. Run-length encoded dyadic depths of the certificate partitions.
support block depth:repeat
xx 1 5:1,6:1,7:1,8:1,9:3,10:14,11:14
xx 2 10:18,11:4,10:19,11:4,10:21,11:4,10:22
xx 3 11:4,10:24,11:4,10:26,11:2,10:27,11:4
xx 4 10:28,11:2,10:29,11:4,10:30,11:2,10:31
xx 5 11:2,10:32,11:2,10:32,11:2,10:32,11:2
xx 6 10:32,11:4,10:32,11:2,10:32,11:2,10:32
xx 7 11:2,10:32,11:2,10:31,11:2,10:30,11:4
xx 8 10:29,11:2,10:28,11:4,10:26,11:4,10:25
xx 9 11:6,10:23,11:6,10:21,11:8,10:18,11:10
xx 10 10:15,11:14,10:12,11:16,10:9,11:18,10:6
xx 11 11:12,12:2,11:11,10:1,11:13,12:4,11:21
xx 12 12:4,11:15,12:10,11:9,12:10,13:2,12:3
xx 13 11:4,12:8,13:6,12:3,11:1,12:3,13:4
xx 14 14:22,15:2,14:2,15:2,14:1,15:2
vv 1 8:1,9:1,10:2,11:17,12:6,11:18,12:2
vv 2 11:11,10:1,11:20,10:6,11:14,10:9,11:12
vv 3 10:12,11:10,10:15,11:8,10:18,11:6,10:20
vv 4 11:4,10:23,11:2,10:25,11:2,10:26,11:2
vv 5 10:27,11:4,10:28,11:2,10:30,11:2,10:30
vv 6 11:2,10:31,11:2,10:65,11:2,10:32,11:2
vv 7 10:66,11:2,10:32,11:2,10:32,11:2,10:32
vv 8 11:2,10:31,11:2,10:30,11:4,10:29,11:2
vv 9 10:28,11:4,10:26,11:4,10:25,11:6,10:23
vv 10 11:6,10:21,11:8,10:18,11:10,10:15,11:14
vv 11 10:12,11:16,10:8,11:20,10:6,11:12,12:2
vv 12 11:26,12:4,11:20,12:8,11:14,12:10,11:9
vv 13 12:10,13:2,12:3,11:4,12:7,13:8,12:3
vv 14 11:1,12:3,13:4,14:22,15:2,14:2,15:2
vv 15 14:1,16:2,15:1
z−z_{-} 1 14:4,13:10,12:2,11:4,12:4,11:9,12:2
z−z_{-} 2 11:6,10:3,11:8,10:6,11:6,10:9,11:4
z−z_{-} 3 10:12,11:2,10:15,11:2,10:45,9:3,10:16
z−z_{-} 4 9:4,10:14,9:6,10:10,9:8,10:10,9:9
z−z_{-} 5 10:8,9:11,10:6,9:12,10:4,9:13,10:4
z−z_{-} 6 9:14,10:4,9:14,10:4,9:14,10:4,9:15
z−z_{-} 7 10:2,9:16,10:2,9:15,10:4,9:15,10:2
z−z_{-} 8 9:16,10:2,9:32,10:2,9:15,10:2,9:173
z+z_{+} 1 2:4
z∘z_{\circ} 1 3:1,5:1,6:2,7:1,8:3,9:4,10:6
z∘z_{\circ} 2 11:30,10:313,9:1,10:8,9:1,10:20,9:1
z∘z_{\circ} 3 10:8,9:1,10:20,9:1,10:8,9:1,10:20
z∘z_{\circ} 4 9:1,10:8,9:1,10:8,9:2,10:8,9:1
z∘z_{\circ} 5 10:8,9:1,10:8,9:2,10:8,9:1,10:8
z∘z_{\circ} 6 9:2,10:6,9:2,10:6,9:2,10:8,9:2
z∘z_{\circ} 7 10:6,9:2,10:6,9:3,10:4,9:3,10:6
z∘z_{\circ} 8 9:3,10:4,9:3,10:4,9:4,10:2,9:4
z∘z_{\circ} 9 10:2,9:93

Appendix F Computer checks

Lemma 4.2 is proved by the interval computation of check_analytic.py, in Arb ball arithmetic [5] through python-flint 0.9.0. The other scripts recheck printed data in exact rational arithmetic. Times are for Python 3.12.3 on one core of an ARM64 Linux server.

Table F.1. Scripts in anc/ and their running times in seconds.
script statements checked time
check_determinant.py (2.5) and Table 2.1; (2.3) for d≠0d\neq 0 and the four determinant evaluations in Lemma 2.2 for small sizes. <1<1
check_arith.py the cut points and every entry of Table A.1; (3.8), 5/485/48, 45/3245/32 and τ\tau; NmN_{m} in Lemma 3.4; the potentials of Lemma 3.3 for n≤12n\leq 12. <10<10
check_analytic.py Lemma 4.2: Tables B.1, B.2, D.1 and D.2, (B.1), (B.2), (4.2); the constants of Lemma 4.4. <90<90
check_mod11.py (5.2), (5.1) and Table 5.1. <1<1

Statements and declarations

Competing interests. The author has no competing interests.

Funding. No funding was received for this work.

Data availability. The scripts of Appendix F, the data of the measure ν\nu and the outputs of the scripts are included in the anc/ attachment accompanying this manuscript.

References

  • [1] F. Calegari, V. Dimitrov, and Y. Tang, The linear independence of 11, ζ⁢(2)\zeta(2), and L⁢(2,χ−3)L(2,\chi_{-3}), preprint (2024), arXiv:2408.15403v2.
  • [2] NIST Digital Library of Mathematical Functions, https://dlmf.nist.gov/, Release 1.2.8 of 2026-09-15, F. W. J. Olver et al., eds.
  • [3] M. Hata, Rational approximations to the dilogarithm, Trans. Amer. Math. Soc. 336 (1993), no. 1, 363–387, https://doi.org/10.1090/S0002-9947-1993-1147401-5. MR1147401
  • [4] C. Krattenthaler, Advanced determinant calculus, Séminaire Lotharingien de Combinatoire 42 (1999), Article B42q, 67 pp., arXiv:math/9902004.
  • [5] F. Johansson, Arb: efficient arbitrary-precision midpoint-radius interval arithmetic, IEEE Trans. Comput. 66 (2017), no. 8, 1281–1292, https://doi.org/10.1109/TC.2017.2690633.
  • [6] T. Prellberg, Three irrationality results for the dilogarithm: Rhin–Viola and Viola–Zudilin constructions at −1/4-1/4, 1/51/5, and −1/3-1/3, preprint (2026), arXiv:2609.06895v1.
  • [7] G. Rhin and C. Viola, The permutation group method for the dilogarithm, Ann. Sc. Norm. Super. Pisa Cl. Sci. (5) 4 (2005), no. 3, 389–437, https://doi.org/10.2422/2036-2145.2005.3.02. MR2185863
  • [8] G. Rhin and C. Viola, Linear independence of 11, Li1\mathrm{Li}_{1} and Li2\mathrm{Li}_{2}, Mosc. J. Comb. Number Theory 8 (2019), no. 1, 81–96, https://doi.org/10.2140/moscow.2019.8.81. MR3864310
  • [9] E. B. Saff and V. Totik, Logarithmic Potentials with External Fields, Grundlehren der Mathematischen Wissenschaften, vol. 316, 1st ed., Springer, 1997, https://doi.org/10.1007/978-3-662-03329-6.
  • [10] C. Viola and W. Zudilin, Linear independence of dilogarithmic values, J. Reine Angew. Math. 736 (2018), 193–223, https://doi.org/10.1515/crelle-2015-0030; preprint version: Max-Planck-Institut für Mathematik Preprint Series 2014 (26), https://archive.mpim-bonn.mpg.de/id/eprint/1914/. MR3769989