Shimura's exact-form identity on the truncated tile: the Petersson-type area integral as a boundary integral
ProvedMTT.PeterssonTile.petersson_tile_stokes_fdLet , let be the truncated standard fundamental-domain tile, and let and be holomorphic on the upper half-plane . Write for and . Then
This is Shimura's exact-form identity (8.2.17) on the tile: since is antiholomorphic, , i.e. , so the area integral of is half the divergence integral of the planar field , and Green's theorem on the tile (green_tile_truncated_fd) turns it into the boundary integral written out over the four boundary pieces (arc, right edge, top, left edge). Applied to a primitive of the period form (coefficientwise) and a second cusp form (so that is the Petersson-type integrand), and summed over the -translates tiling a fundamental domain of (pulled back to with integral_moebius_image_hyperbolic), this is the analytic core of the Eichler–Shimura pairing identity (8.2.22).
Formalization Note The area integral is over Mathlib's regionBetween with volume.prod volume; deriv Fp is the complex derivative. Holomorphy is DifferentiableOn ℂ on the open half-plane; continuity of the derivatives on the tile follows from Cauchy's theory (DifferentiableOn.deriv).
import Mathlib.MeasureTheory.Integral.Prod import Mathlib.MeasureTheory.Integral.IntervalIntegral.FundThmCalculus import Mathlib.MeasureTheory.Integral.IntervalIntegral.IntegrationByParts import Mathlib.MeasureTheory.Measure.Lebesgue.Basic import Mathlib.Analysis.SpecialFunctions.Sqrt import Mathlib.Analysis.Complex.CauchyIntegral set_option autoImplicit false noncomputable section open MeasureTheory Set
theorem MTT.PeterssonTile.petersson_tile_stokes_fd
(H : ℝ) (hH : 1 < H) {Fp a : ℂ → ℂ}
(hFp : DifferentiableOn ℂ Fp {z : ℂ | 0 < z.im})
(ha : DifferentiableOn ℂ a {z : ℂ | 0 < z.im}) :
(∫ p in regionBetween (fun x : ℝ ↦ Real.sqrt (1 - x ^ 2)) (fun _ ↦ H)
(Set.Ioo (-(1 / 2)) (1 / 2)),
deriv Fp ((p.1 : ℂ) + p.2 * Complex.I) *
(starRingEnd ℂ) (a ((p.1 : ℂ) + p.2 * Complex.I)) ∂(volume.prod volume)) =
(1 / 2 : ℂ) *
((∫ x in (-(1 / 2) : ℝ)..(1 / 2),
(-Complex.I * (Fp ((x : ℂ) + H * Complex.I) *
(starRingEnd ℂ) (a ((x : ℂ) + H * Complex.I))) -
-Complex.I * (Fp ((x : ℂ) + Real.sqrt (1 - x ^ 2) * Complex.I) *
(starRingEnd ℂ) (a ((x : ℂ) + Real.sqrt (1 - x ^ 2) * Complex.I))))) +
(∫ x in (-(1 / 2) : ℝ)..(1 / 2),
(Fp ((x : ℂ) + Real.sqrt (1 - x ^ 2) * Complex.I) *
(starRingEnd ℂ) (a ((x : ℂ) + Real.sqrt (1 - x ^ 2) * Complex.I))) *
((-x / Real.sqrt (1 - x ^ 2) : ℝ) : ℂ)) +
(∫ y in (Real.sqrt 3 / 2)..H,
(Fp (((1 / 2 : ℝ) : ℂ) + y * Complex.I) *
(starRingEnd ℂ) (a (((1 / 2 : ℝ) : ℂ) + y * Complex.I)) -
Fp (((-(1 / 2) : ℝ) : ℂ) + y * Complex.I) *
(starRingEnd ℂ) (a (((-(1 / 2) : ℝ) : ℂ) + y * Complex.I))))) := by sorry