Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Square-torus inverse Laplacian: theta formula and Fourier normalization

Open
QAlgorithms.Papers.Independent.Intakethird100.P31.square_torus_green

by Gohan · Nov 1, 2026 · Mathlib 0df444a (Lean v4.33.1)

prove2quantum

Let GTG_TGT​ be the inverse Laplacian on the square torus of side 2π2\pi2π, with its constant mode removed. Construct it independently by integrating the torus heat kernel KsK_sKs​:

GT(x,y)=∫0∞(Ks(x,y)−14π2) ds,Ks(x,y)=14π2∑(k,l)∈Z2e−s(k2+l2)cos⁡(kx+ly).G_T(x,y)=\int_0^\infty\left(K_s(x,y)-\frac1{4\pi^2}\right)\,ds,\qquad K_s(x,y)=\frac1{4\pi^2}\sum_{(k,l)\in\mathbb Z^2}e^{-s(k^2+l^2)}\cos(kx+ly).GT​(x,y)=∫0∞​(Ks​(x,y)−4π21​)ds,Ks​(x,y)=4π21​(k,l)∈Z2∑​e−s(k2+l2)cos(kx+ly).

For (x,y)∉(2πZ)2(x,y)\notin(2\pi\mathbb Z)^2(x,y)∈/(2πZ)2, this integral converges and

GT(x,y)=−12π[log⁡∣θ1(x+iy2π,i)∣−y24π−log⁡η(i)].G_T(x,y)=-\frac1{2\pi}\left[\log\left|\theta_1\left(\frac{x+iy}{2\pi},i\right)\right|-\frac{y^2}{4\pi}-\log\eta(i)\right].GT​(x,y)=−2π1​[log​θ1​(2πx+iy​,i)​−4πy2​−logη(i)].

Here θ1\theta_1θ1​ is the odd Jacobi theta function with first argument of period one, and η(i)>0\eta(i)>0η(i)>0 is the Dedekind eta value. The kernel is integrable on a torus cell. Its normalized Fourier coefficient is zero at (k,l)=(0,0)(k,l)=(0,0)(k,l)=(0,0) and [4π2(k2+l2)]−1[4\pi^2(k^2+l^2)]^{-1}[4π2(k2+l2)]−1 otherwise. The heat and theta series are asserted to converge on their stated domains.

This identifies the Green function used in the geometric coefficients without defining it by the identity being asserted.

Preamble
import Definitions.Def_QAlgorithms_Papers_Independent_Intakethird100_P31
Formal statement
noncomputable section

open Filter Set MeasureTheory

open scoped Topology BigOperators

namespace QAlgorithms.Papers.Independent.Intakethird100.P31

/-- PDF p.27, (5.13)--(5.14): the independently constructed torus inverse
Laplacian has the theta expression and zero constant Fourier mode. The
nonzero Fourier coefficients are `1/(4*pi^2*|k|^2)` with normalized torus
measure. Summability and integrability conclusions prevent default `tsum`
and Bochner-integral values from replacing the kernels. -/
theorem square_torus_green :
    (∀ s : ℝ, 0 < s → ∀ x y : ℝ,
      Summable (fun k : ℤ × ℤ =>
        Real.exp (-s * ((k.1 : ℝ) ^ 2 + (k.2 : ℝ) ^ 2)) *
          Real.cos ((k.1 : ℝ) * x + (k.2 : ℝ) * y))) ∧
    (∀ z : ℂ, Summable (fun j : ℕ => (-1 : ℂ) ^ j *
      (Real.exp (-Real.pi * ((j : ℝ) + 1 / 2) ^ 2) : ℂ) *
        Complex.sin (((2 * j + 1 : ℕ) : ℂ) * (Real.pi : ℂ) * z))) ∧
    0 < etaAtI ∧
    (∀ x y : ℝ, OffSquareLattice x y →
      IntegrableOn (fun s : ℝ => squareHeat s x y - 1 / (4 * Real.pi ^ 2))
        (Ioi 0) ∧
      oddThetaAtI (((x : ℂ) + Complex.I * (y : ℂ)) / (2 * (Real.pi : ℂ))) ≠ 0 ∧
      squareGreen x y = -(1 / (2 * Real.pi)) *
        (Real.log ‖oddThetaAtI (((x : ℂ) + Complex.I * (y : ℂ)) /
          (2 * (Real.pi : ℂ)))‖ - y ^ 2 / (4 * Real.pi) - Real.log etaAtI)) ∧
    IntegrableOn (fun xy : ℝ × ℝ => squareGreen xy.1 xy.2)
      ((Icc (0 : ℝ) (2 * Real.pi)) ×ˢ (Icc (0 : ℝ) (2 * Real.pi))) ∧
    (∀ k l : ℤ,
      (1 / (4 * (Real.pi : ℂ) ^ 2)) *
        (∫ x in Icc (0 : ℝ) (2 * Real.pi),
          ∫ y in Icc (0 : ℝ) (2 * Real.pi),
            (squareGreen x y : ℂ) * Complex.exp
              (-Complex.I * (((k : ℝ) * x + (l : ℝ) * y : ℝ) : ℂ))) =
        if k = 0 ∧ l = 0 then 0 else
          ((1 / (4 * Real.pi ^ 2 * ((k : ℝ) ^ 2 + (l : ℝ) ^ 2)) : ℝ) : ℂ)) := by sorry

end QAlgorithms.Papers.Independent.Intakethird100.P31

end
Source
Entanglement Entropy in the O(N) model, arXiv:0904.4477v1, PDF p. 27, §V A (5.13)–(5.14): zero-mean square-torus Green function
Read-back

What the Lean code literally says, in plain math · gpt-6.1-sol

Define, for real s,x,ys,x,ys,x,y and complex zzz,

K(s,x,y)=14π2∑(m,n)∈Z2e−s(m2+n2)cos⁡(mx+ny),G(x,y)=∫(0,∞)(K(s,x,y)−14π2) ds,Θ(z)=2∑j=0∞(−1)je−π(j+12)2sin⁡ ⁣((2j+1)πz),E=∣η(i)∣,\begin{aligned} K(s,x,y)&=\frac{1}{4\pi^2} \sum_{(m,n)\in\mathbb Z^2}e^{-s(m^2+n^2)}\cos(mx+ny),\\ G(x,y)&=\int_{(0,\infty)} \left(K(s,x,y)-\frac{1}{4\pi^2}\right)\,ds,\\ \Theta(z)&=2\sum_{j=0}^{\infty}(-1)^j e^{-\pi(j+\frac12)^2}\sin\!\bigl((2j+1)\pi z\bigr), \qquad E=|\eta(i)|, \end{aligned}K(s,x,y)G(x,y)Θ(z)​=4π21​(m,n)∈Z2∑​e−s(m2+n2)cos(mx+ny),=∫(0,∞)​(K(s,x,y)−4π21​)ds,=2j=0∑∞​(−1)je−π(j+21​)2sin((2j+1)πz),E=∣η(i)∣,​

where η\etaη is the Dedekind eta function, sums are assigned 000 when their families are not summable, and integrals are assigned 000 when their integrands are not integrable. The statement asserts that for every real s>0s>0s>0 and all real x,yx,yx,y, the family (e−s(m2+n2)cos⁡(mx+ny))(m,n)∈Z2\bigl(e^{-s(m^2+n^2)}\cos(mx+ny)\bigr)_{(m,n)\in\mathbb Z^2}(e−s(m2+n2)cos(mx+ny))(m,n)∈Z2​ is summable; for every complex zzz, the family ((−1)je−π(j+12)2sin⁡((2j+1)πz))j∈N\bigl((-1)^j e^{-\pi(j+\frac12)^2}\sin((2j+1)\pi z)\bigr)_{j\in\mathbb N}((−1)je−π(j+21​)2sin((2j+1)πz))j∈N​ is summable; and E>0E>0E>0. It further asserts that for every real pair (x,y)(x,y)(x,y) for which no integers m,nm,nm,n satisfy x=2πmx=2\pi mx=2πm and y=2πny=2\pi ny=2πn simultaneously, the function s↦K(s,x,y)−1/(4π2)s\mapsto K(s,x,y)-1/(4\pi^2)s↦K(s,x,y)−1/(4π2) is Lebesgue integrable on (0,∞)(0,\infty)(0,∞), Θ((x+iy)/(2π))≠0\Theta((x+iy)/(2\pi))\neq0Θ((x+iy)/(2π))=0, and

G(x,y)=−12π(log⁡∣Θ ⁣(x+iy2π)∣−y24π−log⁡E).G(x,y)=-\frac{1}{2\pi} \left( \log\left|\Theta\!\left(\frac{x+iy}{2\pi}\right)\right| -\frac{y^2}{4\pi}-\log E \right).G(x,y)=−2π1​(log​Θ(2πx+iy​)​−4πy2​−logE).

Also, (x,y)↦G(x,y)(x,y)\mapsto G(x,y)(x,y)↦G(x,y) is Lebesgue integrable on the entire closed square [0,2π]×[0,2π][0,2\pi]\times[0,2\pi][0,2π]×[0,2π], including its lattice corners, and for every pair of integers k,ℓk,\ellk,ℓ,

14π2∫[0,2π](∫[0,2π]G(x,y)e−i(kx+ℓy) dy) dx={0,k=0 and ℓ=0,14π2(k2+ℓ2),otherwise.\frac{1}{4\pi^2} \int_{[0,2\pi]} \left(\int_{[0,2\pi]}G(x,y)e^{-i(kx+\ell y)}\,dy\right)\,dx = \begin{cases} 0,&k=0\ \text{and}\ \ell=0,\\[2pt] \dfrac{1}{4\pi^2(k^2+\ell^2)},&\text{otherwise}. \end{cases}4π21​∫[0,2π]​(∫[0,2π]​G(x,y)e−i(kx+ℓy)dy)dx=⎩⎨⎧​0,4π2(k2+ℓ2)1​,​k=0 and ℓ=0,otherwise.​

This last equality is an equality of complex numbers, with the real values on the right viewed as complex numbers; the second case includes pairs with exactly one coordinate zero.

Human review
  • Endorsed by Shuze Chen · Nov 1, 2026

    Confirmed by the moderator at approval.

  • Endorsed by Gohan · Nov 1, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me