Sphere Packing in R^8

8. Fourier Eigenfunctions with Double Zeroes🔗

First we define the function a. To this end, we consider the following functions.

Definition8.1
Group: Auxiliary modular expressions and integral formulas defining a. (6)
Group member previews
Statement uses 2
Statement dependency previews
Preview
Definition 7.13
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
Used by 5
Reverse dependency previews
markupTeXXL∃∀N

Define \phi_{-4} := \frac{E_4^2}{\Delta} \phi_{-2} := \frac{E_4(E_2 E_4 - E_6)}{\Delta} \phi_{0} := \frac{(E_2 E_4 - E_6)^2}{\Delta}.

The function \phi_0(z) is not modular; however, it satisfies the following transformation rules.

Lemma8.2
Group: Auxiliary modular expressions and integral formulas defining a. (6)
Group member previews
Statement uses 4
Statement dependency previews
Preview
Lemma 7.14
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
Used by 2
Reverse dependency previews
Preview
Corollary 8.13
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXXL∃∀N

We have \phi_0(z + 1) = \phi_0(z) \phi_0\left(-\frac{1}{z}\right) = \phi_0(z)-\frac{12i}{\pi}\,\frac{1}{z}\,\phi_{-2}(z)-\frac{36}{\pi^2}\,\frac{1}{z^2}\,\phi_{-4}(z).

Proof for Lemma 8.2
uses 0

Equation \phi_0(z + 1) = \phi_0(z) follows easily from periodicity of the Eisenstein series and of \Delta(z). For the second transformation law, \phi_{0}\left(-\frac{1}{z}\right) = \frac{(E_2(-1/z) E_4(-1/z) - E_6(-1/z))^{2}}{\Delta(-1/z)} = \frac{((z^2 E_2(z) - 6iz / \pi) \cdot z^4 E_4(z) - z^6 E_6(z))^{2}}{z^{12} \Delta(z)} = \frac{\left(E_2(z) E_4(z) - E_6(z) - \frac{6i}{\pi z} E_4(z)\right)^2}{\Delta(z)} = \frac{(E_2(z) E_4(z) - E_6(z))^{2} - \frac{12i}{\pi z}(E_2(z) E_4(z) - E_6(z)) E_4(z) - \frac{36}{\pi^2 z^2} E_4(z)^{2}}{\Delta(z)} = \phi_0(z) - \frac{12 i}{\pi z} \phi_{-2}(z) - \frac{36}{\pi^2 z^2} \phi_{-4}(z).

For our formalization, we choose rectangular contours for the first and second integrals, and write the definition as follows.

Definition8.3
Group: Auxiliary modular expressions and integral formulas defining a. (6)
Group member previews
uses 1
Used by 4
Reverse dependency previews
Preview
Lemma 8.11
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXL∃∀N

Define a_{\rad} : \R \to \C by a_\rad(r) := I_1(r) + I_2(r) + I_3(r) + I_4(r) + I_5(r) + I_6(r) where I_1(r) := \int_{-1}^{-1 + i} \phi_0 \left(\frac{-1}{z+1}\right) (z + 1)^2 e^{\pi i r z} \dd z I_2(r) := \int_{-1 + i}^{i} \phi_0 \left(\frac{-1}{z+1}\right) (z + 1)^2 e^{\pi i r z} \dd z I_3(r) := \int_{1}^{1 + i} \phi_0 \left(\frac{-1}{z-1}\right) (z - 1)^2 e^{\pi i r z} \dd z I_4(r) := \int_{1 + i}^{i} \phi_0 \left(\frac{-1}{z-1}\right) (z - 1)^2 e^{\pi i r z} \dd z I_5(r) := -2 \int_{0}^{i} \phi_0 \left(\frac{-1}{z}\right) z^2 e^{\pi i r z} \dd z I_6(r) := 2 \int_{i}^{i\infty} \phi_0(z) e^{\pi i r z} \dd z. Here all contours are chosen to be straight line segments. Define a(x) := a_{\rad}(\|x\|^2) for x \in \R^8.

Lean code for Definition8.32 definitions

In the original paper, the integrals I_1 and I_2 are combined, as are I_3 and I_4. We write them separately to simplify the formalization.

Now we prove that a satisfies condition eqn:a-fourier. The following lemma will be used to prove Schwartzness of both a and b.

Lemma8.4
Group: Auxiliary modular expressions and integral formulas defining a. (6)
Group member previews
uses 0
Used by 4
Reverse dependency previews
Preview
Corollary 8.5
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXL∃∀N

Let f(z) be a holomorphic function with Fourier expansion f(z) = \sum_{n \ge n_0} c_f(n) e^{\pi i n z} with c_f(n_0) \ne 0. Assume that c_f(n) has polynomial growth, that is, |c_f(n)| = O(n^k) for some k \in \N. Then there exists a constant C_f > 0 such that \left|\frac{f(z)}{\Delta(z)}\right| \le C_f e^{-\pi (n_0 - 2) \Im z} for all z with \Im z > 1/2.

Lean code for Lemma8.41 theorem
  • theorem MagicFunction.PolyFourierCoeffBound.DivDiscBoundOfPolyFourierCoeff
      (z : UpperHalfPlane) (hz : 1 / 2 < z.im) (c :   ) (n₀ : )
      (hcsum :
        Summable fun i 
          MagicFunction.PolyFourierCoeffBound.fouterm c (↑z) (i + n₀))
      (k : ) (hpoly : c =O[Filter.atTop] fun n  n ^ k)
      (f : UpperHalfPlane  )
      (hf :
         (x : UpperHalfPlane),
          f x =
            ∑' (n : ),
              MagicFunction.PolyFourierCoeffBound.fouterm c (↑x)
                (n + n₀)) :
      f z / Δ z 
        MagicFunction.PolyFourierCoeffBound.DivDiscBound c n₀ *
          Real.exp (-Real.pi * (n₀ - 2) * z.im)
    theorem MagicFunction.PolyFourierCoeffBound.DivDiscBoundOfPolyFourierCoeff
      (z : UpperHalfPlane) (hz : 1 / 2 < z.im)
      (c :   ) (n₀ : )
      (hcsum :
        Summable fun i 
          MagicFunction.PolyFourierCoeffBound.fouterm
            c (↑z) (i + n₀))
      (k : )
      (hpoly :
        c =O[Filter.atTop] fun n  n ^ k)
      (f : UpperHalfPlane  )
      (hf :
         (x : UpperHalfPlane),
          f x =
            ∑' (n : ),
              MagicFunction.PolyFourierCoeffBound.fouterm
                c (↑x) (n + n₀)) :
      f z / Δ z 
        MagicFunction.PolyFourierCoeffBound.DivDiscBound
            c n₀ *
          Real.exp
            (-Real.pi * (n₀ - 2) * z.im)

When f(z) is a quasi-modular form, we can take k to be the weight of f.

Proof for Lemma 8.4

By the product formula for \Delta, \left|\frac{f(z)}{\Delta(z)}\right| = \left|\frac{\sum_{n \ge n_0} c_f(n) e^{\pi i n z}}{e^{2 \pi i z}\prod_{n \ge 1} (1 - e^{2\pi i n z})^{24}}\right| = |e^{\pi i (n_0 - 2)z}| \cdot \frac{|\sum_{n \ge n_0} c_f(n) e^{\pi i (n - n_0) z}|}{\prod_{n \ge 1} |1 - e^{2\pi i n z}|^{24}} \le e^{-\pi (n_0 - 2) \Im z} \cdot \frac{\sum_{n \ge n_0} |c_f(n)| e^{-\pi (n - n_0) \Im z}}{\prod_{n \ge 1} (1 - e^{- 2\pi n \Im z})^{24}} \le e^{-\pi (n_0 - 2) \Im z} \cdot \frac{\sum_{n \ge n_0} |c_f(n)| e^{-\pi (n - n_0) / 2}}{\prod_{n \ge 1} (1 - e^{-\pi n})^{24}} = C_f \cdot e^{-\pi (n_0 - 2) \Im z}. Here C_f = \frac{\sum_{n \ge n_0} |c_f(n)| e^{-\pi (n - n_0) / 2}}{\prod_{n \ge 1} (1 - e^{-\pi n})^{24}}. The numerator converges absolutely because of polynomial growth, and the denominator also converges; indeed it is simply e^{\pi}\cdot \Delta(i/2).

As corollaries, we obtain the following bounds for \phi_0, \phi_{-2}, and \phi_{-4}.

Corollary8.5
Group: Auxiliary modular expressions and integral formulas defining a. (6)
Group member previews
Statement uses 4
Statement dependency previews
Preview
Lemma 7.15
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
Used by 5
Reverse dependency previews
markupTeXXL∃∀N

There exists a constant C_0 > 0 such that |\phi_0(z)| \le C_0 e^{-2 \pi \Im z} for all z with \Im z > 1/2.

Proof for Corollary 8.5

By Ramanujan's formula, E_2 E_4 - E_6 = 3E_4' = 720 \sum_{n \ge 1} n \sigma_3(n) e^{2 \pi i n z}, and (E_2(z) E_4(z) - E_6(z))^{2} = 720^{2} e^{4 \pi i z} + O(e^{5 \pi i z}). The result then follows from Lemma Lemma 8.4 with f(z) = (E_2 E_4 - E_6)^2 and n_0 = 4.

Corollary8.6
Group: Auxiliary modular expressions and integral formulas defining a. (6)
Group member previews
Statement uses 3
Statement dependency previews
Preview
Lemma 7.16
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
used by 1markupTeXXL∃∀N

There exists a constant C_{-2} > 0 such that |\phi_{-2}(z)| \le C_{-2} for all z with \Im z > 1/2.

Corollary8.7
Group: Auxiliary modular expressions and integral formulas defining a. (6)
Group member previews
Statement uses 3
Statement dependency previews
Preview
Lemma 7.16
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
used by 1markupTeXXL∃∀N

There exists a constant C_{-4} > 0 such that |\phi_{-4}(z)| \le C_{-4} e^{2 \pi \Im z} for all z with \Im z > 1/2.

Note that we can take the constants C_0, C_{-2}, and C_{-4} as C_0 = 9 \cdot 240^2 \cdot e^{\pi} \cdot \frac{E_4'(i/2)^{2}}{\Delta(i/2)} C_{-2} = 3 \cdot \frac{E_4(i/2) E_4'(i/2)}{\Delta(i/2)} C_{-4} = e^{-\pi} \cdot \frac{E_4(i/2)^{2}}{\Delta(i/2)}.

Lemma8.8
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of a. (8)
Group member previews
uses 0used by 1markupTeXXL∃∀N

For all n \in \N, there exists a constant C' such that for all r \geq 0, r^n \cdot \int_{1}^{\infty} e^{-2\pi s} \, e^{-\pi r /s} \, \dd s \leq C'.

Proof for Lemma 8.8
uses 0

Fix n \in \N. There exists a constant C such that for all x \ge 0, |x|^n \cdot |e^{-\pi x}| \le C. In particular, for all r \ge 0 and s \ge 1, r^n \cdot e^{-\pi r/s} \le C s^n. Hence r^n \cdot \int_{1}^{\infty} e^{-2\pi s} \, e^{-\pi r /s} \, \dd{s} = \int_{1}^{\infty} e^{-2\pi s} \, \left(|r|^n \cdot e^{-\pi r /s}\right) \, \dd{s} \leq C \int_{1}^{\infty} e^{-2\pi s} \, s^n \, \dd{s}. The Gamma function is given by \Gamma(x) = \int_{0}^{\infty} e^{-u} \, u^{x-1} \, \dd{u} for all x > 0. Writing u = 2\pi s and using \Gamma(n+1) = n!, we obtain C \int_{1}^{\infty} e^{-2\pi s} \, s^n \, \dd{s} \leq C \int_{0}^{\infty} e^{-2\pi s} \, s^n \, \dd{s} = C \int_{0}^{\infty} \frac{1}{(2\pi)^{n+1}} e^{-u} \, u^n \, \dd{u} = \frac{C \cdot n!}{(2\pi)^n}. Defining C' := \frac{C \cdot n!}{(2\pi)^n} finishes the proof.

Lemma8.9
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of a. (8)
Group member previews
uses 0used by 1markupTeXXL∃∀N

There exists a constant C > 0 such that for all r \geq 0, |I_1(r)|, |I_3(r)|, |I_5(r)| \leq C \int_1^{\infty} e^{-2\pi s} \, e^{-\pi r / s} \, \dd s.

Proof for Lemma 8.9

We prove the bound only for I_1(r), as the others are similar. With the change of variable z = -1 + i t for t \in [0,1], we have I_1(r) = -i \int_0^1 \phi_0\left(\frac{-1}{i t}\right) t^2 e^{-\pi i r} \cdot e^{-\pi r t} \, \dd t. With the change of variable s = 1 / t, this becomes I_1(r) = i \int_1^{\infty} \phi_0(i s) \cdot s^{-4} \cdot e^{-\pi i r} \cdot e^{-\pi r / s} \, \dd s, so |I_1(r)| \leq \int_1^{\infty} |\phi_0(i s)| \cdot s^{-4} \cdot |e^{-\pi i r}| \cdot e^{-\pi r / s} \, \dd s \le \int_1^{\infty} |\phi_0(is)| \cdot e^{-\pi r / s} \, \dd s. By Corollary Corollary 8.5, we conclude that |I_1(r)| \leq C_0 \int_1^{\infty} e^{-2\pi s} \, e^{-\pi r / s} \, \dd s.

Lemma8.10
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of a. (8)
Group member previews
uses 0used by 1markupTeXXL∃∀N

There exist C_1, C_2 > 0 such that for all r \geq 0, |I_2(r)|, |I_4(r)| \leq C_1 e^{-\pi r} and |I_6(r)| \leq C_2 \frac{e^{-\pi(r + 2)}}{r + 2}.

Proof for Lemma 8.10

For I_2(r), parametrize z by z = t + i for t \in [-1,0]. Then I_2(r) = \int_{-1}^0 \phi_0\left(\frac{-1}{t + 1 + i}\right) (t + 1 + i)^2 e^{\pi i r t} e^{-\pi r} \, \dd t. Since the function z \mapsto \phi_0\left(\frac{-1}{z+1}\right) (z+1)^2 is holomorphic and bounded on the contour, there exists C > 0 such that |\phi_0\left(\frac{-1}{z+1}\right) (z + 1)^2| \leq C, and therefore |I_2(r)| \le \int_{-1}^{0} C e^{-\pi r} \, \dd t = C e^{-\pi r}. The bound for I_4(r) is similar. /- source paragraph break -/ For I_6(r), parametrize z = i t for t \in [1, \infty), giving I_6(r) = 2 i \int_1^{\infty} \phi_0(i t) e^{-\pi r t} \, \dd t. Using Corollary Corollary 8.5, the absolute value is bounded by |I_6(r)| \leq 2 \int_1^{\infty} |\phi_0(i t)| e^{-\pi r t} \, \dd t \leq \frac{2C_0}{\pi} \int_1^{\infty} e^{-2\pi t} e^{-\pi r t} \, \dd t = \frac{2C_0}{\pi} \frac{e^{-\pi (r + 2)}}{r + 2}.

Combining these bounds, one can show that a is a Schwartz function.

Lemma8.11
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of a. (8)
Group member previews
Statement uses 2
Statement dependency previews
Preview
Definition 8.3
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
Used by 2
Reverse dependency previews
Preview
Lemma 8.12
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXL∃∀N

a(x) is a Schwartz function.

Lean code for Lemma8.111 definition
  • def MagicFunction.FourierEigenfunctions.a :
      SchwartzMap (EuclideanSpace  (Fin 8)) 
    def MagicFunction.FourierEigenfunctions.a :
      SchwartzMap (EuclideanSpace  (Fin 8)) 
    The +1-Fourier Eigenfunction of Viazovska's Magic Function. 
Proof for Lemma 8.11
Proof uses 4
Proof dependency previews
Preview
Theorem 5.8
Loading preview
Proof dependency preview content is loaded from the rendered-fragment cache.

By Theorem Theorem 5.8, it suffices to show that the function is smooth and decays faster than any polynomial. Smoothness follows from differentiation under the integral, already formalized in mathlib. The decay follows from Lemmas Lemma 8.9 and Lemma 8.10, together with Lemma 8.8.

Lemma8.12
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of a. (8)
Group member previews
Statement uses 5
Statement dependency previews
Used by 2
Reverse dependency previews
Preview
Lemma 8.26
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXL∃∀N

a(x) satisfies equation eqn:a-fourier.

Lean code for Lemma8.121 theorem
  • theorem MagicFunction.a.Fourier.eig_a :
      (FourierTransform.fourierCLE 
            (SchwartzMap (EuclideanSpace  (Fin 8)) ))
          MagicFunction.FourierEigenfunctions.a =
        MagicFunction.FourierEigenfunctions.a
    theorem MagicFunction.a.Fourier.eig_a :
      (FourierTransform.fourierCLE 
            (SchwartzMap
              (EuclideanSpace  (Fin 8)) ))
          MagicFunction.FourierEigenfunctions.a =
        MagicFunction.FourierEigenfunctions.a
Proof for Lemma 8.12
uses 0

We recall that the Fourier transform of a Gaussian function is \mathcal{F}(e^{\pi i \|x\|^2 z})(y)=z^{-4}\,e^{\pi i \|y\|^2 \,(\frac{-1}{z}) }. Next, exchange contour integration in the z variable with Fourier transform in the x variable in the definition of a. This is justified because the corresponding double integral converges absolutely. In this way we obtain \widehat{a}(y)=\int\limits_{-1}^i\phi_0\Big(\frac{-1}{z+1}\Big)\,(z+1)^2\,z^{-4}\,e^{\pi i \|y\|^2 \,(\frac{-1}{z})}\,dz +\int\limits_{1}^i\phi_0\Big(\frac{-1}{z-1}\Big)\,(z-1)^2\,z^{-4}\,e^{\pi i \|y\|^2 \,(\frac{-1}{z})}\,dz -2\int\limits_{0}^i\phi_0\Big(\frac{-1}{z}\Big)\,z^2\,z^{-4}\,e^{\pi i \|y\|^2 \,(\frac{-1}{z})}\,dz +2\int\limits_{i}^{i\infty}\phi_0(z)\,z^{-4}\,e^{\pi i \|y\|^2 \,(\frac{-1}{z})}\,dz. Making the change of variables w=\frac{-1}{z}, we arrive at \widehat{a}(y)=\int\limits_{1}^i\phi_0\Big(1-\frac{1}{w-1}\Big)\,(\frac{-1}{w}+1)^2\,w^{2}\,e^{\pi i \|y\|^2 \,w}\,dw +\int\limits_{-1}^i\phi_0\Big(1-\frac{1}{w+1}\Big)\,(\frac{-1}{w}-1)^2\,w^2\,e^{\pi i \|y\|^2 \,w}\,dw -2\int\limits_{i \infty}^i\phi_0(w)\,e^{\pi i \|y\|^2 \,w}\,dw +2\int\limits_{i}^{0}\phi_0\Big(\frac{-1}{w}\Big)\,w^{2}\,e^{\pi i \|y\|^2 \,w}\,dw. Since \phi_0 is 1-periodic we have \begin{aligned} \widehat{a}(y)\,=\,&\int\limits_{1}^i\phi_0\Big(\frac{-1}{z-1}\Big)\,(z-1)^2\,e^{\pi i \|y\|^2 \,z}\,dz +\int\limits_{-1}^i\phi_0\Big(\frac{-1}{z+1}\Big)\,(z+1)^2\,e^{\pi i \|y\|^2 \,z}\,dz \\ +&2\int\limits_{i}^{i\infty}\phi_0(z)\,e^{\pi i \|y\|^2 \,z}\,dz -2\int\limits_{0}^{i}\phi_0\Big(\frac{-1}{z}\Big)\,z^{2}\,e^{\pi i \|y\|^2 \,z}\,dz \\ \,=\,&a(y). \end{aligned} This finishes the proof of the proposition.

Next, we check that a has double zeroes at all \Lambda_8-lattice points of length greater then \sqrt{2}. Using eqn:phi0-bound, eqn:phi2-bound, and eqn:phi4-bound, we can control the behavior of \phi_0 near 0 and i\infty.

Corollary8.13
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of a. (8)
Group member previews
Statement uses 4
Statement dependency previews
Preview
Lemma 8.2
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
used by 1markupTeXXL∃∀N

We have \phi_0\left(\frac{i}{t}\right) = O(e^{-2 \pi / t}) \quad \text{as } t \to 0 \phi_0\left(\frac{i}{t}\right) = O(t^{-2}e^{2 \pi t}) \quad \text{as } t \to \infty.

Proof for Corollary 8.13
uses 0

The first estimate follows from eqn:phi0-bound with z = i/t. For the second estimate, by eqn:phi0-trans-S, eqn:phi2-bound, and eqn:phi4-bound, \left|\phi_0\left(\frac{i}{t}\right)\right| \le |\phi_0(it)| + \frac{12}{\pi t} |\phi_{-2}(it)| + \frac{36}{\pi^2 t^2} |\phi_{-4}(it)| \le C_0 e^{-2 \pi t} + \frac{12}{\pi t} \cdot C_{-2} + \frac{36}{\pi^2 t^2} \cdot C_{-4} e^{2 \pi t} = O(t^{-2}e^{2 \pi t}).

Lemma8.14
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of a. (8)
Group member previews
Statement uses 3
Statement dependency previews
Preview
Corollary 7.26
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
Used by 2
Reverse dependency previews
Preview
Lemma 8.15
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXXL∃∀N

For r>\sqrt{2} we can express a(r) in the form a(r)=-4\sin(\pi r^2/2)^2\,\int\limits_{0}^{i\infty}\phi_0\Big(\frac{-1}{z}\Big)\,z^2\,e^{\pi i r^2 \,z}\,dz.

Proof for Lemma 8.14

Denote the right-hand side by d(r). Convergence of the integral for r > \sqrt{2} follows from Corollary Corollary 8.13. We can write d(r)=\int\limits_{-1}^{i\infty-1}\phi_0\Big(\frac{-1}{z+1}\Big)\,(z+1)^2\,e^{\pi i r^2 \,z}\,dz- 2\int\limits_{0}^{i\infty}\phi_0\Big(\frac{-1}{z}\Big)\,z^2\,e^{\pi i r^2 \,z}\,dz +\int\limits_{1}^{i\infty+1}\phi_0\Big(\frac{-1}{z-1}\Big)\,(z-1)^2\,e^{\pi i r^2 \,z}\,dz. From eqn:phi0-trans-S we deduce that if r>\sqrt{2}, then \phi_0\Big(\frac{-1}{z}\Big)\,z^2\,e^{\pi i r^2 \,z}\to 0 as \Im(z)\to\infty. Therefore we can deform the paths of integration and rewrite d(r)=\int\limits_{-1}^{i}\phi_0\Big(\frac{-1}{z+1}\Big)\,(z+1)^2\,e^{\pi i r^2 \,z}\,dz +\int\limits_{i}^{i\infty}\phi_0\Big(\frac{-1}{z+1}\Big)\,(z+1)^2\,e^{\pi i r^2 \,z}\,dz -2\int\limits_{0}^{i}\phi_0\Big(\frac{-1}{z}\Big)\,z^2\,e^{\pi i r^2 \,z}\,dz -2\int\limits_{i}^{i\infty}\phi_0\Big(\frac{-1}{z}\Big)\,z^2\,e^{\pi i r^2 \,z}\,dz +\int\limits_{1}^{i}\phi_0\Big(\frac{-1}{z-1}\Big)\,(z-1)^2\,e^{\pi i r^2 \,z}\,dz +\int\limits_{i}^{i\infty}\phi_0\Big(\frac{-1}{z-1}\Big)\,(z-1)^2\,e^{\pi i r^2 \,z}\,dz. Now eqn:phi0-trans-S implies \phi_0\Big(\frac{-1}{z+1}\Big)\,(z+1)^2-2\phi_0\Big(\frac{-1}{z}\Big)\,z^2+ \phi_0\Big(\frac{-1}{z-1}\Big)\,(z-1)^2=2\phi_0(z), and hence d(r)=\int\limits_{-1}^{i}\phi_0\Big(\frac{-1}{z+1}\Big)\,(z+1)^2\,e^{\pi i r^2 \,z}\,dz -2\int\limits_{0}^{i}\phi_0\Big(\frac{-1}{z}\Big)\,z^2\,e^{\pi i r^2 \,z}\,dz +\int\limits_{1}^{i}\phi_0\Big(\frac{-1}{z-1}\Big)\,(z-1)^2\,e^{\pi i r^2 \,z}\,dz +2\int\limits_{i}^{i\infty}\phi_0(z)\,e^{\pi i r^2 \,z}\,dz=a(r).

Finally, we find another convenient integral representation for a and compute values of a(r) at r=0 and r=\sqrt{2}.

Lemma8.15
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of a. (8)
Group member previews
Statement uses 4
Statement dependency previews
Preview
Definition 8.1
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
Used by 3
Reverse dependency previews
Preview
Lemma 8.16
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXXL∃∀N

For r\geq0 we have a(r)=4i\,\sin(\pi r^2/2)^2\,\Bigg(\frac{36}{\pi^3\,(r^2-2)}-\frac{8640}{\pi^3\,r^4}+\frac{18144}{\pi^3\,r^2} +\int\limits_0^\infty\,\left(t^2\,\phi_0\Big(\frac{i}{t}\Big)-\frac{36}{\pi^2}\,e^{2\pi t}+\frac{8640}{\pi}\,t-\frac{18144}{\pi^2}\right)\,e^{-\pi r^2 t}\,dt \Bigg). The integral converges absolutely for all r\in\R_{\geq 0}.

Proof for Lemma 8.15

Suppose that r>\sqrt{2}. Then by Proposition prop:a-double-zeros, a(r)=4i\,\sin(\pi r^2/2)^2\,\int\limits_{0}^{\infty}\phi_0(i/t)\,t^2\,e^{-\pi r^2 t}\,dt. From eqn:phi0-trans-S we obtain the asymptotic expansion \phi_0(i/t)\,t^2=\frac{36}{\pi^2}\,e^{2 \pi t}-\frac{8640}{\pi}\,t+\frac{18144}{\pi^2}+O(t^2\,e^{-2\pi t}) as t\to\infty. For r>\sqrt{2}, \int\limits_0^\infty \left(\frac{36}{\pi^2}\,e^{2 \pi t}+\frac{8640}{\pi}\,t+\frac{18144}{\pi^2}\right)\,e^{-\pi r^2 t}\,dt =\frac{36}{\pi^3\,(r^2-2)}-\frac{8640}{\pi^3\,r^4}+\frac{18144}{\pi^3\,r^2}. Therefore the claimed identity holds for r>\sqrt{2}. /- source paragraph break -/ On the other hand, the definition of a shows that a(r) is analytic in a neighborhood of [0,\infty), and the asymptotic expansion above implies that the right-hand side is also analytic there. Hence the identity holds on the whole interval [0,\infty).

From eqn:a-another-integral we see that the values a(r) lie in i\R for all r\in\R_{\geq0}.

Lemma8.16
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of a. (8)
Group member previews
uses 1used by 1markupTeXL∃∀N

We have a(0) = -\frac{i}{8640}.

Lean code for Lemma8.161 theorem, incomplete
  • contains sorry
    theorem MagicFunction.a.SpecialValues.a_zero :
      MagicFunction.FourierEigenfunctions.a 0 = -8640 * Complex.I / Real.pi
    theorem MagicFunction.a.SpecialValues.a_zero :
      MagicFunction.FourierEigenfunctions.a
          0 =
        -8640 * Complex.I / Real.pi
Proof for Lemma 8.16
uses 0

These identities follow immediately from the previous proposition.

Now we construct function b. To this end we consider the function.

Definition8.17
Group: Auxiliary forms and integral formulas defining b. (4)
Group member previews
Preview
Definition 8.18
Loading preview
Group member preview content is loaded from the rendered-fragment cache.
uses 1used by 1markupTeXXL∃∀N

Define h(z) := 128 \frac{H_3(z) + H_4(z)}{H_2(z)^2}.

It is easy to see that h\in M^!_{-2}(\Gamma_0(2)). Indeed, first check that h|_{-2}\gamma=h for all \gamma\in\Gamma_0(2). Since \Gamma_0(2) is generated by \left(\begin{smallmatrix}1&0\\2&1\end{smallmatrix}\right) and \left(\begin{smallmatrix}1&1\\0&1\end{smallmatrix}\right), it suffices to check invariance under their action, which follows immediately from the transformation laws of the theta-null forms and eqn: h define. Next analyze the poles of h. It is known that \theta_{10} has no zeros in the upper half-plane, so h has poles only at the cusps. At the cusp i\infty, this modular form has Fourier expansion h(z)\,=\,q^{-1} + 16 - 132 q + 640 q^2 - 2550 q^3+O(q^4).

Let I=\left(\begin{smallmatrix}1&0\\0&1\end{smallmatrix}\right), T=\left(\begin{smallmatrix}1&1\\0&1\end{smallmatrix}\right), and S=\left(\begin{smallmatrix}0&-1\\1&0\end{smallmatrix}\right) be elements of \Gamma_1.

Definition8.18
Group: Auxiliary forms and integral formulas defining b. (4)
Group member previews
Preview
Definition 8.17
Loading preview
Group member preview content is loaded from the rendered-fragment cache.
uses 1
Used by 5
Reverse dependency previews
markupTeXXL∃∀N

Define \psi_I\,:=\,h-h|_{-2}ST \psi_T\,:=\,\psi_I|_{-2}T \psi_S\,:=\,\psi_I|_{-2}S.

Lemma8.19
Group: Auxiliary forms and integral formulas defining b. (4)
Group member previews
Preview
Definition 8.17
Loading preview
Group member preview content is loaded from the rendered-fragment cache.
uses 0
Used by 2
Reverse dependency previews
Preview
Lemma 8.20
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXXL∃∀N

\psi_I(z), \psi_S(z), and \psi_T(z) can be written as \psi_I(z) = \frac{H_4^3 (5 H_2^2 + 5 H_2 H_4 + 2 H_4^2)}{2\Delta} \psi_S(z) = -\frac{H_2^3 (2 H_2^2 + 5 H_2 H_4 + 5 H_4^2)}{2 \Delta} \psi_T(z) = \frac{H_3^3 (5 H_2^2 - 5 H_2 H_3 + 2 H_3^2)}{2 \Delta}.

Proof for Lemma 8.19
Proof uses 2
Proof dependency previews
Preview
Lemma 7.33
Loading preview
Proof dependency preview content is loaded from the rendered-fragment cache.

By the transformation laws of the theta-null functions, H_2|_{-2}ST = (-H_4)|_{-2}T = -H_3 H_3|_{-2}ST = (-H_3)|_{-2}T = -H_4 H_4|_{-2}ST = (-H_2)|_{-2}T = H_2. Using these identities and the level-one/level-two identities, we can rewrite \psi_I(z) as \psi_I(z) = h(z) - h|_{-2}ST(z) = 128 \frac{H_3 + H_4}{H_2^2} - 128 \frac{-H_4 + H_2}{H_3^2} = 128 \frac{H_3^2 (H_3 + H_4) - H_2^2 (H_2 - H_4)}{H_2^2 H_3^2} = 128 \frac{(H_2 + H_4)^2 (H_2 + 2 H_4) - H_2^2 (H_2 - H_4)}{H_2^2 H_3^2} = 128 \frac{H_4(5 H_2^2 + 5 H_2 H_4 + 2 H_4^2)}{ H_2^2 H_3^2} = \frac{H_4^3 (5 H_2^2 + 5 H_2 H_4 + 2 H_4^2)}{2\Delta}. Applying |_{-2}S and |_{-2}T to the expression for \psi_I proves the formulas for \psi_S and \psi_T.

Lemma8.20
Group: Auxiliary forms and integral formulas defining b. (4)
Group member previews
Preview
Definition 8.17
Loading preview
Group member preview content is loaded from the rendered-fragment cache.
uses 1
Used by 2
Reverse dependency previews
Preview
Lemma 8.28
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXXL∃∀N

The Fourier expansions of these functions are \psi_I(z)\,=\,q^{-1} + 144 + O(q^{1/2}) \psi_T(z)\,=\,q^{-1} + 144 + O(q^{1/2}).

Now we are ready to define b. Again, the definition here is slightly different from that in Viazovska (2017), where all contours are chosen to be straight lines.

Definition8.21
Group: Auxiliary forms and integral formulas defining b. (4)
Group member previews
Preview
Definition 8.17
Loading preview
Group member preview content is loaded from the rendered-fragment cache.
uses 1
Used by 2
Reverse dependency previews
Preview
Lemma 8.26
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXXL∃∀N

Define b_\rad : \R \to \C by b_\rad(r) := J_1(r) + J_2(r) + J_3(r) + J_4(r) + J_5(r) + J_6(r) where J_1(r) := \int_{-1}^{-1 + i} \psi_T(z) e^{\pi i r z} \, \dd z J_2(r) := \int_{-1 + i}^{i} \psi_T(z) e^{\pi i r z} \, \dd z J_3(r) := \int_{1}^{1 + i} \psi_T(z) e^{\pi i r z} \, \dd z J_4(r) := \int_{1 + i}^{i} \psi_T(z) e^{\pi i r z} \, \dd z J_5(r) := -2 \int_{0}^{i} \psi_I(z) e^{\pi i r z} \, \dd z J_6(r) := -2 \int_{i}^{i \infty} \psi_S(z) e^{\pi i r z} \, \dd z. Here all contours are straight line segments. Define b : \R^8 \to \C by b(x) := b_\rad(\|x\|^2).

Now we prove that b is a Schwartz function and satisfies eqn:b-fourier. As in the case of a(x), we start with an exponential bound for \psi_S(z) and \psi_I(z).

Lemma8.22
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of b. (8)
Group member previews
uses 0
Used by 2
Reverse dependency previews
Preview
Lemma 8.25
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXXL∃∀N

There exist constants C_I, C_S, C_T > 0 such that |\psi_I(z)| \le C_I e^{2\pi \Im z} |\psi_T(z)| \le C_T e^{2\pi \Im z} |\psi_S(z)| \le C_S e^{- \pi \Im z}.

Proof for Lemma 8.22
Proof uses 2
Proof dependency previews
Preview
Lemma 8.4
Loading preview
Proof dependency preview content is loaded from the rendered-fragment cache.

The proof is similar to that of Lemma cor:phi0-bound and follows from Lemma Lemma 8.4 together with the fact that the vanishing orders of the numerators of \psi_I, \psi_T, and \psi_S at infinity are respectively 0, 0, and \frac{3}{2}.

Lemma8.23
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of b. (8)
Group member previews
uses 0used by 1markupTeXXL∃∀N

There exists a constant C > 0 such that |J_1(r)|, |J_3(r)|, |J_5(r)| \le C \int_1^{\infty} e^{-\pi s} e^{\pi r / s}\, \dd s.

Lemma8.24
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of b. (8)
Group member previews
uses 0used by 1markupTeXXL∃∀N

There exist C_1, C_2 > 0 such that |J_2(r)|, |J_4(r)| \le C_1 e^{-\pi r} |J_6(r)| \le C_2 \frac{e^{\pi (r + 1)}}{r + 1}.

Lemma8.25
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of b. (8)
Group member previews
uses 1used by 1markupTeXL∃∀N

b(x) is a Schwartz function.

Lean code for Lemma8.251 definition
  • def MagicFunction.FourierEigenfunctions.b :
      SchwartzMap (EuclideanSpace  (Fin 8)) 
    def MagicFunction.FourierEigenfunctions.b :
      SchwartzMap (EuclideanSpace  (Fin 8)) 
    The -1-Fourier Eigenfunction of Viazovska's Magic Function. 
Proof for Lemma 8.25
Proof uses 4
Proof dependency previews
Preview
Theorem 5.8
Loading preview
Proof dependency preview content is loaded from the rendered-fragment cache.

Similar to the proof of Lemma 8.11.

Lemma8.26
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of b. (8)
Group member previews
Statement uses 4
Statement dependency previews
Preview
Lemma 5.2
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
used by 1markupTeXL∃∀N

b(x) satisfies equation eqn:b-fourier.

Lean code for Lemma8.261 theorem
  • theorem MagicFunction.b.Fourier.eig_b :
      (FourierTransform.fourierCLE 
            (SchwartzMap (EuclideanSpace  (Fin 8)) ))
          MagicFunction.FourierEigenfunctions.b =
        -MagicFunction.FourierEigenfunctions.b
    theorem MagicFunction.b.Fourier.eig_b :
      (FourierTransform.fourierCLE 
            (SchwartzMap
              (EuclideanSpace  (Fin 8)) ))
          MagicFunction.FourierEigenfunctions.b =
        -MagicFunction.FourierEigenfunctions.b
Proof for Lemma 8.26

We repeat the argument used in the proof of Proposition prop:a-fourier. Using the Gaussian Fourier identity and exchanging the contour integration in z with the Fourier transform in x, we get \mathcal{F}(b)(x)= \int\limits_{-1}^{i}\psi_T(z)\,z^{-4}\,e^{\pi i \|x\|^2 (\frac{-1}{z})}\,dz + \int\limits_{1}^{i}\psi_T(z)\,z^{-4}\,e^{\pi i \|x\|^2 (\frac{-1}{z})}\,dz - 2\,\int\limits_{0}^{i}\psi_I(z)\,z^{-4}\,e^{\pi i \|x\|^2 (\frac{-1}{z})}\,dz - 2\,\int\limits_{i}^{i\infty}\psi_S(z)\,z^{-4}\,e^{\pi i \|x\|^2 (\frac{-1}{z})}\,dz. With the change of variables w=\frac{-1}{z}, we arrive at \mathcal{F}(b)(x)= \int\limits_{1}^{i}\psi_T\Big(\frac{-1}{w}\Big)\,w^{2}\,e^{\pi i \|x\|^2 w}\,dw + \int\limits_{-1}^{i}\psi_T\Big(\frac{-1}{w}\Big)\,w^{2}\,e^{\pi i \|x\|^2 w}\,dw - 2\,\int\limits_{i\infty}^{i}\psi_I\Big(\frac{-1}{w}\Big)\,w^{2}\,e^{\pi i \|x\|^2 w}\,dw - 2\,\int\limits_{i}^{0}\psi_S\Big(\frac{-1}{w}\Big)\,w^{2}\,e^{\pi i \|x\|^2 w}\,dw. Now the definitions imply \psi_T|_{-2}S=-\psi_T \psi_I|_{-2}S=\psi_S \psi_S|_{-2}S=\psi_I, so this becomes \mathcal{F}(b)(x)= \int\limits_{1}^{i}-\psi_T(z)\,e^{\pi i \|x\|^2 z}\,dz + \int\limits_{-1}^{i}-\psi_T(z)\,e^{\pi i \|x\|^2 z}\,dz + 2\,\int\limits_{i}^{i\infty}\psi_S(z)\,e^{\pi i \|x\|^2 z}\,dz + 2\,\int\limits_{0}^{i}\psi_I(z)\,e^{\pi i \|x\|^2 z}\,dw. From the definition of b, we conclude that \mathcal{F}(b)(x)=-b(x).

Now we regard the radial function b as a function on \R_{\geq0} and check that it has double roots at the \Lambda_8 points.

Corollary8.27
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of b. (8)
Group member previews
Statement uses 2
Statement dependency previews
Preview
Definition 8.18
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
used by 1markupTeXXL∃∀N

We have \psi_I(it) = O(t^2 e^{\pi/t}) \quad \text{as } t \to 0 \psi_I(it) = O(e^{2 \pi t}) \quad \text{as } t \to \infty.

Proof for Corollary 8.27

By eqn:psiS-define, \psi_I(it) = (it)^{-2} \psi_S\left(\frac{-1}{it}\right) = -t^{-2} \psi_S\left(\frac{i}{t}\right), and combined with eqn:psiS-bound this gives eqn:psiI-near-0. Equation eqn:psiI-near-infty follows from lemma:psi-bound.

Lemma8.28
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of b. (8)
Group member previews
Statement uses 4
Statement dependency previews
Used by 2
Reverse dependency previews
Preview
Lemma 8.29
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXXL∃∀N

For r>\sqrt{2}, the function b(r) can be expressed as b(r)=-4\sin(\pi r^2/2)^2\,\int\limits_{0}^{i\infty}\psi_I(z)\,e^{\pi i r^2 \,z}\,dz.

Proof for Lemma 8.28

Denote the right-hand side by c(r). By Corollary cor:psiI-near-0-infty, the integral converges for r>\sqrt{2}. Rewrite it as c(r)=\int\limits_{-1}^{i\infty-1}\psi_I(z+1)\,e^{\pi i r^2 \,z}\,dz-2\int\limits_{0}^{i\infty}\psi_I(z)\,e^{\pi i r^2 \,z}\,dz+ \int\limits_{1}^{i\infty+1}\psi_I(z-1)\,e^{\pi i r^2 \,z}\,dz. From the Fourier expansion of \psi_I, we know that \psi_I(z)=e^{-2\pi i z}+O(1) as \Im(z)\to\infty. Because r^2>2, we can deform the paths of integration and write \int\limits_{-1}^{i\infty-1}\psi_I(z+1)\,e^{\pi i r^2 \,z}\,dz= \int\limits_{-1}^{i}\psi_T(z)\,e^{\pi i r^2 \,z}\,dz+\int\limits_{i}^{i\infty}\psi_T(z)\,e^{\pi i r^2 \,z}\,dz \int\limits_{1}^{i\infty+1}\psi_I(z-1)\,e^{\pi i r^2 \,z}\,dz= \int\limits_{-1}^{i}\psi_T(z)\,e^{\pi i r^2 \,z}\,dz+\int\limits_{i}^{i\infty}\psi_T(z)\,e^{\pi i r^2 \,z}\,dz. Hence c(r)=\int\limits_{-1}^{i}\psi_T(z)\,e^{\pi i r^2 \,z}\,dz+\int\limits_{1}^{i}\psi_T(z)\,e^{\pi i r^2 \,z}\,dz -2\int\limits_{0}^{i}\psi_I(z)\,e^{\pi i r^2 \,z}\,dz +2\int\limits_{i}^{i\infty}(\psi_T(z)-\psi_I(z))\,e^{\pi i r^2 \,z}\,dz. Next, the functions \psi_I, \psi_T, and \psi_S satisfy \psi_T+\psi_S=\psi_I. Using this identity, we find c(r)=\int\limits_{-1}^{i}\psi_T(z)\,e^{\pi i r^2 \,z}\,dz+\int\limits_{1}^{i}\psi_T(z)\,e^{\pi i r^2 \,z}\,dz -2\int\limits_{0}^{i}\psi_I(z)\,e^{\pi i r^2 \,z}\,dz -2\int\limits_{i}^{i\infty}\psi_S(z)\,e^{\pi i r^2 \,z}\,dz=b(r).

At the end of this section we find another integral representation of b(r) for r\in\R_{\geq0} and compute special values of b.

Lemma8.29
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of b. (8)
Group member previews
Statement uses 3
Statement dependency previews
Preview
Lemma 8.20
Loading preview
Statement dependency preview content is loaded from the rendered-fragment cache.
Used by 2
Reverse dependency previews
Preview
Lemma 8.30
Loading preview
Reverse dependency preview content is loaded from the rendered-fragment cache.
markupTeXXL∃∀N

For r\geq0 we have b(r)=4i\,\sin(\pi r^2/2)^2\,\left(\frac{144}{\pi\,r^2}+\frac{1}{\pi\,(r^2-2)}+\int\limits_0^\infty\,\left(\psi_I(it)-144-e^{2\pi t}\right)\,e^{-\pi r^2 t}\,dt\right). The integral converges absolutely for all r\in\R_{\geq 0}.

Proof for Lemma 8.29
Proof uses 2
Proof dependency previews
Preview
Lemma 8.15
Loading preview
Proof dependency preview content is loaded from the rendered-fragment cache.

The proof is analogous to the proof of Proposition prop:a-another-integral. First suppose that r>\sqrt{2}. Then by Proposition prop:b-double-zeros, b(r)=4i\,\sin(\pi r^2/2)^2\,\int\limits_{0}^{\infty}\psi_I(it)\,e^{-\pi r^2 t}\,dt. From the Fourier expansion of \psi_I we obtain \psi_I(it)=e^{2\pi t}+144+O(e^{-\pi t}) as t\to\infty. For r>\sqrt{2}, \int\limits_0^\infty \left(e^{2\pi t}+144\right)\,e^{-\pi r^2 t}\,dt =\frac{1}{\pi\,(r^2-2)}+\frac{144}{\pi\,r^2}. Therefore the identity holds for r>\sqrt{2}. /- source paragraph break -/ On the other hand, the definition of b shows that b(r) is analytic in a neighborhood of [0,\infty), and the asymptotic expansion above implies that the right-hand side is also analytic there. Hence the identity holds on the whole interval [0,\infty).

From eqn:b-another-integral we see that b(r)\in i\R for all r\in\R_{\geq 0}.

Lemma8.30
Group: Schwartz, Fourier-eigenfunction, and double-zero properties of b. (8)
Group member previews
uses 1used by 1markupTeXL∃∀N

We have b(0) = 0.

Lean code for Lemma8.301 theorem, incomplete
  • contains sorry
    theorem MagicFunction.b.SpecialValues.b_zero :
      MagicFunction.FourierEigenfunctions.b 0 = 0
    theorem MagicFunction.b.SpecialValues.b_zero :
      MagicFunction.FourierEigenfunctions.b
          0 =
        0
Proof for Lemma 8.30
uses 0

These identities follow immediately from the previous proposition.