Square-torus inverse Laplacian: theta formula and Fourier normalization
OpenQAlgorithms.Papers.Independent.Intakethird100.P31.square_torus_greenLet be the inverse Laplacian on the square torus of side , with its constant mode removed. Construct it independently by integrating the torus heat kernel :
For , this integral converges and
Here is the odd Jacobi theta function with first argument of period one, and is the Dedekind eta value. The kernel is integrable on a torus cell. Its normalized Fourier coefficient is zero at and 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.
import Definitions.Def_QAlgorithms_Papers_Independent_Intakethird100_P31
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
Read-back
What the Lean code literally says, in plain math · gpt-6.1-sol
Define, for real and complex ,
where is the Dedekind eta function, sums are assigned when their families are not summable, and integrals are assigned when their integrands are not integrable. The statement asserts that for every real and all real , the family is summable; for every complex , the family is summable; and . It further asserts that for every real pair for which no integers satisfy and simultaneously, the function is Lebesgue integrable on , , and
Also, is Lebesgue integrable on the entire closed square , including its lattice corners, and for every pair of integers ,
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.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.