Estimates for 1728·J and s₂ (Milla, arXiv:1809.00533v6, Chapter 5) #
This file states the explicit estimates and q-series approximations of Chapter 5
("Estimates") of Milla's proof of the Chudnovsky formula (arXiv:1809.00533v6):
Chudnovsky.Jtilde: the truncationJtildeof Klein'sJ-invariant (TheoremtheonaeherJof the paper),Jtilde(τ) = (1 + 240(q + 9q²))³ / (1728·q·(1 - q - q²)²⁴);Chudnovsky.s₂tilde: the truncationstilde₂of Ramanujan'ss₂(Theoremtheonaehers2of the paper),stilde₂(τ) = (1 + 240(q + 9q²))/(1 - 504(q + 33q²)) · (1 - 24(q + 3q²) - 3/(π·Im τ));Chudnovsky.E₄trunc,Chudnovsky.E₆trunc,Chudnovsky.E₂starTrunc: the quadratic truncationsX,Y,Zof Lemmalemxyof the paper;Chudnovsky.eisensteinTail: the Lambert-series tailsR_w⁽ˡ⁾ = ∑_{n≥l} σ_{w-1}(n)·qⁿ(Lemmalemrestchaetz), witheisensteinTail k lcorresponding to weightw = k + 1;Chudnovsky.sigma_le_pow_succ:σₖ(n) ≤ n^(k+1)(Lemmasigmaschaetz);Chudnovsky.norm_q_lt_of_mem_Region:|q| < e^(-7.852)forIm τ > 5/4(Lemmaarchim-lem, with Mathlib'sπ > 3.1415replacing Archimedes);Chudnovsky.norm_eisensteinTail_leand the concrete boundsnorm_eisensteinTail_sigma₁/₃/₅(Lemmaslemrestchaetz,lemrestkonkr);Chudnovsky.lemE6:|E₆(τ)| > 0.8on the region (LemmalemE6);Chudnovsky.theonaeherJ:|1728·J - 1728·Jtilde| < 0.2together with the bracketing0.737/|q| < |1728·J| < 1.321/|q|and|J| > 1.096(TheoremtheonaeherJ);Chudnovsky.theonaehers2:|s₂ - stilde₂| < 222000·|q|³(Theoremtheonaehers2).
All estimates hold on Chudnovsky.Region = {τ | Im τ > 5/4}.
Explicit q-truncations #
The quadratic truncations X, Y, Z of the Eisenstein series (paper Lemma lemxy),
and the resulting approximations Jtilde (paper Thm. theonaeherJ) and stilde₂
(paper Thm. theonaehers2). All are honest rational functions of the nome q τ
(plus, for Z, the real quantity 3/(π·Im τ)), so the numerical phase can evaluate
them directly.
The quadratic truncation X := E₄⁽²⁾ = 1 + 240(q + 9q²) of E₄ (paper Lemma lemxy).
Equations
- Chudnovsky.E₄trunc τ = 1 + 240 * (Chudnovsky.q τ + 9 * Chudnovsky.q τ ^ 2)
Instances For
The quadratic truncation Y := E₆⁽²⁾ = 1 - 504(q + 33q²) of E₆ (paper Lemma lemxy).
Equations
- Chudnovsky.E₆trunc τ = 1 - 504 * (Chudnovsky.q τ + 33 * Chudnovsky.q τ ^ 2)
Instances For
The quadratic truncation
Z := E₂⁽²⁾ - 3/(π·Im τ) = 1 - 24(q + 3q²) - 3/(π·Im τ) of E₂* (paper Lemma lemxy).
Equations
- Chudnovsky.E₂starTrunc τ = 1 - 24 * (Chudnovsky.q τ + 3 * Chudnovsky.q τ ^ 2) - 3 / (↑Real.pi * ↑τ.im)
Instances For
The approximation Jtilde of Klein's J-invariant from Theorem theonaeherJ of the paper:
Jtilde(τ) = (1 + 240(q + 9q²))³ / (1728·q·(1 - q - q²)²⁴).
The denominator (1 - q - q²)²⁴ is the truncation of ∏ (1-qⁿ)²⁴ = Δ/q coming from
Euler's pentagonal number theorem (paper Remark penta).
Equations
- Chudnovsky.Jtilde τ = (1 + 240 * (Chudnovsky.q τ + 9 * Chudnovsky.q τ ^ 2)) ^ 3 / (1728 * Chudnovsky.q τ * (1 - Chudnovsky.q τ - Chudnovsky.q τ ^ 2) ^ 24)
Instances For
The approximation stilde₂ of Ramanujan's s₂ from Theorem theonaehers2 of the paper:
stilde₂(τ) = (1 + 240(q + 9q²))/(1 - 504(q + 33q²)) · (1 - 24(q + 3q²) - 3/(π·Im τ)).
Equations
- One or more equations did not get rendered due to their size.
Instances For
The σₖ bound (paper Lemma sigmaschaetz) #
Paper Lemma sigmaschaetz: σₖ(n) ≤ n^(k+1).
(The paper states this for n ≥ 1; with Mathlib's convention σ k 0 = 0 it holds for
all n.)
The Archimedes lemma (paper Lemma archim-lem) #
Paper Lemma archim-lem: if Im τ > 5/4 then |q| < e^(-7.852).
Archimedes' bound π > 3 + 10/71 is replaced by Mathlib's Real.pi_gt_d4.
q-expansions of the Eisenstein series (paper Prop. sigmaeisen) #
We derive the divisor-sum q-expansions of E₄, E₆ and E₂ from Mathlib's
EisensteinSeries.q_expansion_bernoulli and EisensteinSeries.E2_eq_tsum_cexp,
in a ℕ-indexed form suitable for splitting off the leading terms.
Summability of the divisor-sum q-series ∑ σₖ(n)·qⁿ (paper Prop. sigmaeisen).
Summability of the (real, absolute) divisor-sum q-series ∑ σₖ(n)·|q|ⁿ.
Paper Prop. sigmaeisen: E₄(τ) = 1 + 240·∑_{n≥1} σ₃(n)·qⁿ.
Paper Prop. sigmaeisen: E₆(τ) = 1 - 504·∑_{n≥1} σ₅(n)·qⁿ.
Paper Prop. sigmaeisen: E₂(τ) = 1 - 24·∑_{n≥1} σ₁(n)·qⁿ.
Split off the first two terms of the divisor-sum q-series, leaving the tail
∑_{n≥3} σₖ(n)·qⁿ (the paper's R_{k+1}⁽³⁾).
Numerical bound e^(-7.852) < 0.000389, replacing the Archimedes step of the paper.
Proved from Real.exp_one_gt_d9 and an 8-term Taylor lower bound for exp 0.852.
Geometric tail bounds for the Lambert-type series #
(paper Lemmas lemrestchaetz and lemrestkonkr)
The tail of the divisor-sum (Lambert-type) q-series of the Eisenstein series of
weight k + 1 (paper Lemma lemrestchaetz):
eisensteinTail k l τ = ∑_{n ≥ l} σₖ(n)·qⁿ, i.e. the paper's R_w⁽ˡ⁾ with w = k + 1.
Equations
- Chudnovsky.eisensteinTail k l τ = ∑' (n : ℕ), ↑((ArithmeticFunction.sigma k) (n + l)) * Chudnovsky.q τ ^ (n + l)
Instances For
The tail R⁽ˡ⁾ peels off its leading term: R⁽ˡ⁾ = σₖ(l)·qˡ + R⁽ˡ⁺¹⁾.
|q| > 0 on the region (q τ = e^{2πiτ} ≠ 0).
Paper Lemma lemrestkonkr, first bound: |R₂⁽³⁾| ≤ 4.007·|q|³ for Im τ > 5/4.
Paper Lemma lemrestkonkr, second bound: |R₄⁽³⁾| ≤ 28.1·|q|³ for Im τ > 5/4.
Paper Lemma lemrestkonkr, third bound: |R₆⁽³⁾| ≤ 245.6·|q|³ for Im τ > 5/4.
Truncation errors (paper δX, δY, δZ of Lemma lemxy) #
The Eisenstein series equal their quadratic truncations plus 240/-504/-24 times the
corresponding tail R⁽³⁾.
E₆ = Y + δY with δY = -504·R₆⁽³⁾ (paper Lemma lemxy).
E₂* = Z + δZ with δZ = -24·R₂⁽³⁾ (paper Lemma lemxy).
Paper Lemma lemE6: |E₆(τ)| > 0.8 for Im τ > 5/4.
Paper Lemma lemE6, "in particular": E₆(τ) ≠ 0 for Im τ > 5/4.
Truncation norm bounds (paper Lemma lemxy) #
Paper Lemma lemxy: 0.8014 ≤ |Y| where Y = E₆⁽²⁾.
Paper Lemma lemxy: |X| ≤ 1.0937 where X = E₄⁽²⁾.
Paper Lemma lemxy: |Z| ≤ 1.0094 where Z = E₂⁽²⁾ - 3/(π Im τ).
Paper: |stilde₂| ≤ 1.3776.
The analytic function k(τ) = (E₄³ - E₆²)/(1728·q) = Δ/q (paper Lemma lemk).
Equations
- Chudnovsky.kfun τ = (ModularForm.E₄ τ ^ 3 - ModularForm.E₆ τ ^ 2) / (1728 * Chudnovsky.q τ)
Instances For
The truncation ktilde(τ) = (1 - q - q²)²⁴ (paper Lemma lemk).
Equations
- Chudnovsky.ktilde τ = (1 - Chudnovsky.q τ - Chudnovsky.q τ ^ 2) ^ 24
Instances For
1728·q·J = E₄³ / k.
Paper Lemma lemk: 0.9907 ≤ |ktilde|.
Paper Lemma lemk: |ktilde| ≤ 1.0094.
Paper Lemma lemxy: |Y| ≤ 1.1987 where Y = E₆⁽²⁾.
Horner-form norm bound: for ‖p‖ ≤ ρ (0 ≤ ρ), the norm of a polynomial
∑ⱼ cⱼ·pʲ (written in Horner form via List.foldr) is bounded by the same fold
with |cⱼ| and ρ. Used to bound the degree-46 polynomial remainder in lemk.
Coefficients (ascending) of the degree-46 polynomial G in the exact identity
X³ - Y² - 1728·q·(1-q-q²)²⁴ = 1728·q³·G(q), used in lemk.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Paper Lemma lemk (the one analytic input still admitted): the difference between
Klein's k = Δ/q and its truncation ktilde = (1-q-q²)²⁴ is O(|q|²):
|k - ktilde| ≤ 365.6·|q|².
Milla proves this by Taylor-expanding (1-q-q²)²⁴ to second order (bounding the third
derivative) and comparing with the explicit expansion of (X³-Y²)/(1728q); the sharper
|k - ktilde| < 25|q|⁵ needs Euler's pentagonal number theorem. Here we instead use the exact
polynomial identity X³ - Y² - 1728·q·(1-q-q²)²⁴ = 1728·q³·G(q) (G of degree 46, with
|G(q)| ≤ 179.8 bounded via horner_norm_bound), together with the truncation-error
estimates |E₄³-X³| ≤ 24202|q|³ and |E₆²-Y²| ≤ 296780|q|³.
Paper Lemma lemxy: 0.9063 ≤ |X| where X = E₄⁽²⁾.
Theorem theonaeherJ #
Paper Theorem theonaeherJ, headline approximation bound:
|1728·J(τ) - 1728·Jtilde(τ)| < 0.2 for Im τ > 5/4.
Paper Theorem theonaeherJ, lower bracketing:
0.737/|q| < |1728·J(τ)| for Im τ > 5/4.
Paper Theorem theonaeherJ, upper bracketing:
|1728·J(τ)| < 1.321/|q| for Im τ > 5/4.
Paper Theorem theonaeherJ: |J(τ)| > 1.096 for Im τ > 5/4.
|J(τ)| > 1 for Im τ > 5/4; in particular J(τ) ≠ 0 and J(τ) ≠ 1 there.