From sum of two squares to arithmetic Siegel–Weil formulas

Abstract

The main goal of this expository article is to survey recent progress on the arithmetic Siegel–Weil formula and its applications. We begin with the classical sum of two squares problem and put it in the context of the Siegel–Weil formula. We then motivate the geometric and arithmetic Siegel–Weil formula using the classical example of the product of modular curves. After explaining the recent result on the arithmetic Siegel–Weil formula for Shimura varieties of arbitrary dimension, we discuss some aspects of the proof and its application to the arithmetic inner product formula and the Beilinson–Bloch conjecture. Rather than being intended as a complete survey of this vast field, this article focuses more on examples and background to provide easier access to several recent works by the author with W. Zhang and Y. Liu.

1. Sum of two squares

1.1. Which prime $p$ can be written as the sum of two squares?

For the first few primes we easily find that

$$\begin{equation*} 5=1^2+2^2, \quad 13=2^2+3^2,\quad 17=1^2+4^2,\quad 29=2^2+5^2 \end{equation*}$$

are the sums of two squares, while other primes like $3, 7, 11, 19,23$ are not. The answer seems to depend on the residue class of $p$ modulo 4.

Theorem 1.1.1.

A prime $p\ne 2$ is the sum of two squares if and only if $p\equiv 1\ (\mathrm{mod}\ 4)$.

Theorem 1.1.1 is usually attributed to Fermat and appeared in his letter to Mersenne dated Dec 25, 1640 (hence the name Fermat’s Christmas Theorem), although the statement can already be found in the work of Girard in 1625. The “only if” direction is obvious, but the “if” direction is far from trivial. Fermat claimed that he had an irrefutable proof, but nobody was able to find the complete proof among his work—apparently margins were often too narrow for Fermat. The only clue (in his letters to Pascal and to Digby) is that he used a “descent argument”: if such a prime $p$ is not of the required form, then one can construct another smaller prime and so on, until a contradiction occurs when one encounters 5, the smallest such prime. More than 100 years later, Euler (1755) gave the first rigorous proof of Theorem 1.1.1 based on infinite descent. For a detailed history of Theorem 1.1.1, see Dickson Dic66, Ch. VI, pp. 227–231.

1.2. Which positive integer $n$ can be be written as the sum of two squares?

If $n=x^2+y^2$ ($x,y\in \mathbb{Z}$) and $p\mid n$, then either $p\mid \mathrm{gcd}(x,y)$ or $(x/y)^2\equiv -1\ (\mathrm{mod}\ p)$, and hence either $n/p^2$ is also the sum of two squares or $p\not \equiv 3\ (\mathrm{mod}\ 4)$ (by the quadratic reciprocity). It follows that each $p\mid n$ with $p\equiv 3\ (\mathrm{mod}\ 4)$ must appear to an even power. On the other hand, the familiar Diophantus identity

$$\begin{equation*} \left(a^2 + b^2\right)\left(c^2 + d^2\right) =\left(ac-bd\right)^2 + \left(ad+bc\right)^2 \end{equation*}$$

shows that a product of integers of the form $x^2+y^2$ is also of the same form. Combining with Theorem 1.1.1 we obtain:

Corollary 1.2.1.

A positive integer $n$ is of the form $n=x^2+y^2$ if and only if each prime factor $p\equiv 3\ (\mathrm{mod}\ 4)$ of $n$ appears to an even power.

1.3. In how many different ways can one represent $n$ as the sum of two squares?

Definition 1.3.1.

To answer this question, we naturally define the representation number

$$\begin{equation*} r(n)\coloneq \#\{(x,y)\in \mathbb{Z}^2: n=x^2+y^2\}. \end{equation*}$$

In particular, $n$ is of the form $x^2+y^2$ if and only if $r(n)>0$.

Example 1.3.2.

$$\begin{align*} 4&= 0^2+(\pm 2)^2=(\pm 2)^2+0^2, & r(4)=4,\\ 5&=(\pm 1)^2+(\pm 2)^2=(\pm 2)^2+(\pm 1)^2, &r(5)=8,\\ 25&=0^2+(\pm 5)^2=(\pm 3)^2+(\pm 4)^2=(\pm 4)^2+(\pm 3)^2, &r(25)=12. \\ \end{align*}$$

In his book Fundamenta nova theoriae functionum ellipticarum (1829), Jacobi proved the following general formula for the representation numbers.

Theorem 1.3.3 (Jacobi).

$$\begin{equation*} r(n)=4\left(\sum _{d|n\atop d\equiv 1\bmod 4}1-\sum _{d|n\atop d\equiv 3\bmod 4}1\right). \end{equation*}$$

As a byproduct, Jacobi’s formula shows that

$$\begin{align*} p\equiv 1\ (\mathrm{mod}\ 4)&\Rightarrow \ r(p)=4(2-0)=8>0, \\ p\equiv 3\ (\mathrm{mod}\ 4)&\Rightarrow \ r(p)=4(1-1)=0, \end{align*}$$

which gives an immediate (and different) proof of Theorem 1.1.1!

1.4. Jacobi’s proof

Jacobi’s proof of Theorem 1.3.3 involves Jacobi’s theta series,

$$\begin{equation*} \theta \coloneq \sum _{n\in \mathbb{Z}}q^{n^2}=1+2q+2q^4+2q^9+\cdots . \end{equation*}$$

The representation numbers $r(n)$ naturally appear as the $n$th coefficients of the square of Jacobi’s theta series

$$\begin{equation*} \theta ^2=\left(\sum _{n\in \mathbb{Z}}q^{n^2}\right)^2=\sum _{n\ge 0}r(n)q^n=1 + 4 q + 4 q^2 + 4 q^4 + 8 q^5 +4 q^8+\cdots . \end{equation*}$$

Jacobi used his theory of elliptic functions (including his famous Triple Product Identity) to derive the formula (Jac1820, p. 107)

$$\begin{equation*} \theta ^2=1+4\left(\frac{q}{1-q}-\frac{q^3}{1-q^3}+\frac{q^5}{1-q^5}-\frac{q^7}{1-q^7}+\cdots \right), \end{equation*}$$

which is easily seen to be equivalent to Theorem 1.3.3.

1.5. Another proof using modular forms

An alternative way of evaluating $\theta ^2$ is to view $q=e^{2\pi i\tau }$, and $\theta =\theta (\tau )$ as a holomorphic function on the upper half-plane

$$\begin{equation*} \mathcal{H}\coloneq \{\tau =\mathsf{x}+i\mathsf{y}\in \mathbb{C}:\mathsf{y}=\operatorname {Im}(\tau )>0\}. \end{equation*}$$

The function $\theta (\tau )$ satisfies two transformation rules (see Zag08, Proposition 9):

$$\begin{equation*} \theta (\tau +1)=\theta (\tau ),\quad \theta \left(-\frac{1}{4\tau }\right)=(-2i\tau )^{1/2}\theta (\tau ). \end{equation*}$$

The first rule is clear by the periodicity of the exponential function. The second rule can be proved using the Poisson summation formula and also plays a key role in Riemann’s proof of the functional equation of the Riemann zeta function (see DS05, §4.9). These rules amount to saying that

$$\begin{equation*} \theta (\tau )\in M_{1/2}(\Gamma _1(4)) \end{equation*}$$

is a modular form of weight $1/2$ and level $\Gamma _1(4)$. Jacobi’s theta series $\theta (\tau )$ and its variants (under the general name of theta series) form one of most important classes of modular forms.

It follows that

$$\begin{equation*} \theta ^2(\tau )\in M_1(\Gamma _1(4)) \end{equation*}$$

is a modular form of weight 1 and level $\Gamma _1(4)$. The space $M_1(\Gamma _1(4))$ is in fact one dimensional (Zag08, Proposition 3 or DS05, Theorem 3.6.1), so if one can construct another a modular form of weight 1 and level $\Gamma _1(4)$, then it has to be a scalar multiple of $\theta ^2(\tau )$. We next construct such a modular form using Eisenstein series, another of the most important classes of modular forms.

Definition 1.5.1.

Let $\chi : (\mathbb{Z}/4 \mathbb{Z})^\times {\overset{\sim }{\longrightarrow }} \{\pm 1\}$ be the (unique) nontrivial character. We define an Eisenstein series

$$\begin{equation} G_k^{\chi }(\tau )=\sum _{(0,0)\ne (c,d)\in \mathbb{Z}^2, 4|c}\frac{\chi (d)}{(c\tau +d)^k}, {\tag{1.5.1.1}} \end{equation}$$

where $\chi (d)$ is understood to be 0 when $(4,d)\ne 1$.

When $k\ge 3$, the series 1.5.1.1 is absolutely convergent and is nonzero only when $k$ is odd. When $k\ge 3$ is odd, it defines a modular form $G_{k}^\chi (\tau )\in M_k(\Gamma _1(4))$ of weight $k$, level $\Gamma _1(4)$, and character $\chi$. The constant term of the $q$-expansion of $G_k^\chi (\tau )$ is nonzero, and we let $E_k^{\chi }(\tau )$ be a scalar multiple of $G_k^\chi (\tau )$ so the constant term is normalized to be 1. This normalized Eisenstein series $E_k^\chi (\tau )$ then has the explicit $q$-expansion (see DS05, §4.5)

$$\begin{equation} E_k^{\chi }(\tau )=1+c_k^\chi \cdot \sum _{n\ge 1}\left(\sum _{d|n}\chi (d)d^{k-1}\right)q^n, {\tag{1.5.1.2}} \end{equation}$$

where $c_k^\chi =2/L(1-k,\chi )$ is related to a special value of the Dirichlet $L$-function $L(s,\chi )$.

When $k=1$, the series 1.5.1.1 is not absolutely convergent, but one can still suitably modify it to obtain a modular form

$$\begin{equation*} E_1^\chi (\tau )\in M_1(\Gamma _1(4)) \end{equation*}$$

with the same formula 1.5.1.2 for its $q$-expansion, either using the Weierstrass $\sigma$-function (see DS05, §4.8) or using the analytic continuation of

$$\begin{equation} G_1^{\chi }(\tau ,s)\coloneq \sum _{(0,0)\ne (c,d)\in \mathbb{Z}^2, 4|c}\frac{\chi (d)}{(c\tau +d)|c\tau +d|^{2s}},\quad {\mathrm{Re}}(s)>1/2 {\tag{1.5.1.3}} \end{equation}$$

to $s=0$ (see Miy89, §7.2). In particular, when $k=1$, formula 1.5.1.2 simplifies to

$$\begin{equation*} E_1^\chi (\tau )=1+c_1^\chi \cdot \sum _{n\ge 1}\left(\sum _{d|n}\chi (d)\right)q^n. \end{equation*}$$

As explained, $E_1^\chi (\tau )$ must be a scalar multiple of $\theta ^2(\tau )$. Since both of them have constant coefficient 1, we indeed have the equality

$$\begin{equation} \theta ^2(\tau )=E_1^\chi (\tau ). {\tag{1.5.1.4}} \end{equation}$$

Comparing the coefficient before $q$, we obtain $c_1^\chi =4$ and hence

$$\begin{equation*} r(n)=4\sum _{d|n}\chi (d), \end{equation*}$$

which proves Theorem 1.3.3.

Remark 1.5.2.

As a byproduct of the proof, we also obtain $L(0,\chi )=\frac{1}{2}$ from $c_1^\chi =4$. Via the functional equation of $L(s,\chi )$, this is equivalent to the famous Leibniz formula for $\pi$ (1676):

$$\begin{equation*} L(1,\chi )=1-\frac{1}{3}+\frac{1}{5}-\frac{1}{7}+\cdots =\frac{\pi }{4}. \end{equation*}$$

To summarize, Jacobi’s Theorem 1.3.3 can be proved using the identity of two modular forms 1.5.1.4, namely using a relation of the form

theta series $\longleftrightarrow$ Eisenstein series.

Notice that the Fourier coefficients of theta series encode representation numbers of quadratic forms, while the Fourier coefficients of Eisenstein series are generalized divisor sums which are more explicit.

2. Siegel–Weil formula

2.1. Siegel’s formula

Siegel Sie35 generalizes formula 1.5.1.4 from the binary quadratic form $x^2+y^2$ to more general quadratic forms in an arbitrary number of variables. Let $\Lambda$ be a positive definite quadratic lattice over $\mathbb{Z}$ of rank $m$ with quadratic form $Q: \Lambda \rightarrow \mathbb{Z}$. Denote by $(\ ,\ )$ the associated symmetric bilinear form, defined by

$$\begin{equation*} (x,y)\coloneq Q(x+y)-Q(x)-Q(y) \end{equation*}$$

(so $Q(x)=\frac{1}{2}(x,x)$). Denote by $\operatorname {Sym}_n(\mathbb{Z})$ the set of symmetric matrices whose diagonal entries are in $\mathbb{Z}$ and whose off-diagonal entries are in $\frac{1}{2}\mathbb{Z}$. Denote by $\operatorname {Sym}_n(\mathbb{Z})_{\ge 0}\subseteq \operatorname {Sym}_n(\mathbb{Z})$ the subset of positive semidefinite matrices.

Definition 2.1.1.

For $T\in \operatorname {Sym}_n(\mathbb{Z})_{\ge 0}$, define the (generalized) representation number

$$\begin{equation*} r_\Lambda (T)\coloneq \#\{(x_1,\ldots ,x_n)\in \Lambda ^n: {\frac{1}{2}}((x_i,x_j))_{1\le i,j\le n}=T\}. \end{equation*}$$

Define Siegel’s theta series

$$\begin{equation} \theta _\Lambda (\tau )\coloneq \sum _{T\in \operatorname {Sym}_n(\mathbb{Z})_{\ge 0}}r_\Lambda (T)q^T, \quad q^T\coloneq e^{2\pi i\operatorname {tr}T\tau }, {\tag{2.1.1.1}} \end{equation}$$

and a holomorphic function on Siegel’s half-space

$$\begin{equation*} \mathcal{H}_n\coloneq \{\tau =\mathsf{x}+i \mathsf{y}: \mathsf{x}\in \operatorname {Sym}_n(\mathbb{R}), \mathsf{y}\in \operatorname {Sym}_n(\mathbb{R})_{>0}\}. \end{equation*}$$

Using the Poisson summation formula, Siegel proved that $\theta _\Lambda (\tau )$ is a Siegel modular form on $\mathcal{H}_n$ of weight $m/2$.

Example 2.1.2.

Notice that when $n=1$, Siegel’s half-space $\mathcal{H}_n$ recovers the usual upper half-plane $\mathcal{H}$, and Siegel’s theta series $\theta _\Lambda (\tau )$ recovers Jacobi’s theta series

$$\begin{equation*} \theta _\Lambda (\tau )\coloneq \sum _{n\ge 0}r_\Lambda (n)q^n,\quad r_\Lambda (n)\coloneq \{x\in \Lambda : Q(x)=n\}. \end{equation*}$$

In general, the theta series $\theta _\Lambda$ for a lattice $\Lambda$ may fail to be an Eisenstein series on the nose, but Siegel’s formula shows that the weighted average of theta series within its genus class is always a Siegel Eisenstein series on $\mathcal{H}_n$:

weighted average of theta series $\longleftrightarrow$ Siegel Eisenstein series.

More precisely, recall that two quadratic lattices $\Lambda ,\Lambda '$ are in the same genus, if they are isomorphic over $\mathbb{R}$ and over $\mathbb{Z}_p$ for all primes $p$. Denote by $\mathrm{Gen}(\Lambda )$ the set of isomorphism classes of quadratic lattices in the same genus of $\Lambda$. Denote by $\operatorname {Aut}(\Lambda )$ the automorphism group of $\Lambda$ as a quadratic lattice.

Theorem 2.1.3 (Siegel).

The following identity holds:

$$\begin{equation} \frac{ \sum _{\Lambda '\in \mathrm{Gen}(\Lambda )}\frac{1}{\#\operatorname {Aut}(\Lambda ')}\cdot \theta _{\Lambda '}(\tau )}{ \sum _{\Lambda '\in \mathrm{Gen}(\Lambda )}\frac{1}{\#\operatorname {Aut}(\Lambda ')}} =E_\Lambda (\tau ). {\tag{2.1.3.1}} \end{equation}$$

Here $E_\Lambda (\tau )$ is a certain normalized Siegel Eisenstein series on $\mathcal{H}_n$ of weight $m/2$.

Example 2.1.4.

Consider the case $m=2$, $n=1$ and $\Lambda =\mathbb{Z}^2$ equipped with the quadratic form $Q=x^2+y^2$. Then

$$\begin{equation*} \theta _\Lambda (\tau )=\theta ^2(\tau ),\quad E_\Lambda (\tau )= E_1^\chi (\tau ). \end{equation*}$$

In this case $\mathrm{Gen}(\Lambda )$ is a singleton and Siegel’s formula recovers 1.5.1.4.

Example 2.1.5 (cf. Ser73, V.2.3).

Siegel’s formula is extremely useful in studying the arithmetic of quadratic forms. For example, one can deduce his famous mass formula (also known as the Smith–Minkowski–Siegel mass formula), which computes the mass of $\mathrm{Gen}(\Lambda )$, defined to be weighted size

$$\begin{equation*} \sum _{\Lambda '\in \mathrm{Gen}(\Lambda )}\frac{1}{\#\operatorname {Aut}(\Lambda ')} \end{equation*}$$

as an Euler product of local factors indexed by primes $p$.

For example, consider the simplest case when $\Lambda$ is

•

unimodular, i.e., $\det (\frac{1}{2}((x_i,x_j)_{i,j=1}^n)\in \{\pm 1\}$ for a $\mathbb{Z}$-basis $\{x_1,\ldots ,x_n\}$ of $\Lambda$, and

•

even, i.e., 2 divides $Q(x)$ for all $x\in \Lambda$.

The rank $m$ of any unimodular even lattice $\Lambda$ is necessarily a multiple of 8. Siegel’s mass formula computes the mass of $\mathrm{Gen}(\Lambda )$ explicitly as

$$\begin{align*} \sum _{\Lambda '\in \mathrm{Gen}(\Lambda )}\frac{1}{\#\operatorname {Aut}(\Lambda ')}&=\frac{2\zeta (2)\zeta (4)\cdots \zeta (m-2)\zeta (m/2)}{\operatorname {vol}(S^0)\operatorname {vol}(S^1)\cdots \operatorname {vol}(S^{m-1})} = {B_{{m/2}} \over m}\prod _{{1\leq j<m/2}}{B_{{2j}} \over 4j}, \end{align*}$$

where $\zeta (s)$ is the Riemann zeta function, $\operatorname {vol}(S^{k-1})=\frac{2\pi ^{k/2}}{\Gamma (k/2)}$ is the volume of the unit $(k-1)$-sphere, and $B_k$ is the $k$th Bernoulli number.

Example 2.1.6 (cf. Ser73, VII.6.6).

Let $\Lambda$ be the root lattice of type $E_8$, defined by

$$\begin{equation*} \Lambda \coloneq \{x=(x_1,\ldots , x_8)\in \mathbb{Z}^8\cup (\mathbb{Z}+1/2)^8: \sum _{i=1}^8x_i\equiv 0\ (\mathrm{mod}\ 2)\},\quad Q=\sum _{i=1}^8 x_i^2. \end{equation*}$$

Then $\Lambda$ is unimodular and even. Siegel’s mass formula computes that

$$\begin{equation*} \sum _{\Lambda '\in \mathrm{Gen}(\Lambda )}\frac{1}{\#\operatorname {Aut}(\Lambda ')}={B_{4} \over 8}{B_{2} \over 4}{B_{4} \over 8}{B_{6} \over 12}={-1/30 \over 8}\;{1/6 \over 4}\;{-1/30 \over 8}\;{1/42 \over 12}={1 \over 696729600}. \end{equation*}$$

The fact that $\#\operatorname {Aut}(\Lambda )=696729600$ (which is also the order of the Weyl group of type $E_8$) then implies that $\Lambda$ is the unique unimodular and even lattice of rank 8.

In this case $E_\Lambda (\tau )$ is related to the classical Eisenstein series $E_4(\tau )\in M_4({\mathrm{SL}}_2(\mathbb{Z}))$ of weight 4 and level 1 by

$$\begin{equation*} \begin{split} E_\Lambda (\tau )&=E_4(2\tau )\\ &=1+240\sum _{n\ge 1}\left(\sum _{d|n}d^3\right)q^{2n} =1+240q^2+2160q^4+6720q^6+\cdots . \end{split} \end{equation*}$$

Here, the factor $240$ agrees with the constant $c_k=2/\zeta (1-k)$ for the weight $k=4$. Siegel’s formula then implies that

$$\begin{equation*} r_\Lambda (2n)\coloneq \#\left\{x=(x_i)\in \Lambda : \sum _{i=1}^8 x_i^2=2n\right\}=240\sum _{d|n}d^3. \end{equation*}$$

In particular, we recover that there are $r_\Lambda (2)=240$ roots in the $E_8$ root system (see Figure 1).

2.2. Siegel–Weil formula

In a series of works (Sie36Sie37Sie51Sie52), Siegel further generalized his formula from definite to indefinite quadratic forms and from the base field $\mathbb{Q}$ to totally real fields. In the indefinite case the formula is more difficult to state, as the theta series $\theta _\Lambda (\tau )$ is divergent, and it is necessary to introduce extra weight functions in the definition to ensure convergence. The situation was greatly clarified by Weil Wei65 using powerful tools from representation theory, especially due to his use of the Weil representation.

Let $V$ be a quadratic space over $\mathbb{Q}$ of dimension $m$ with bilinear form $(\ ,\ )$. For simplicity we assume that $m$ is even (so the weight $m/2$ of the relevant Siegel modular forms is integral; see Remark 2.2.12). Consider the reductive dual pair $(G, H)=({\mathrm{Sp}}(2n), \mathrm{O}(V))$, where ${\mathrm{Sp}}(2n)$ is the symplectic group of the standard $2n$-dimensional symplectic space over $\mathbb{Q}$ and $\mathrm{O}(V)$ is the orthogonal group of $V$. Let $P=MN\subseteq G={\mathrm{Sp}}(2n)$ be the standard Siegel parabolic subgroup, so that under the standard basis we have

$$\begin{align*} M(\mathbb{Q})&=\left\{m(a)=\begin{pmatrix} a & 0\\ 0 &{}^ta^{-1}\end{pmatrix}: a\in \mathrm{GL}_n(\mathbb{Q})\right\},\\ N(\mathbb{Q}) &= \left\{n(b)=\begin{pmatrix} 1_n & b \\ 0 & 1_n\end{pmatrix}: b\in \operatorname {Sym}_n(\mathbb{Q})\right\}. \end{align*}$$

Definition 2.2.1.

Let $\mathbb{A}$ be the ring of adèles of $\mathbb{Q}$. We fix the standard additive character $\psi : \mathbb{A}\rightarrow \mathbb{C}^\times$ whose archimedean component is given by $\psi _\infty : \mathbb{A}_\infty =\mathbb{R}\rightarrow \mathbb{C}^\times , x\mapsto e^{2\pi i x}$. The (Schrödinger model of the) Weil representation $\omega =\omega _{V,\psi }$ is the representation of $G(\mathbb{A})\times H(\mathbb{A})$ on the space of Schwartz functions $\mathscr{S} (V(\mathbb{A})^n)$ such that for any $\varphi \in \mathscr{S} (V(\mathbb{A})^n)$ and $\mathbf{x}\in V(\mathbb{A})^n$,

$$\begin{align} \begin{aligned} \omega (m(a))\varphi (\mathbf{x})&=\chi _V(\det a)|\det a|^{m/2}\varphi (\mathbf{x}\cdot a),\quad &m(a)\in M(\mathbb{A}),\\ \omega (n(b))\varphi (\mathbf{x})&=\psi (\frac{1}{2}\operatorname {tr}b(\mathbf{x},\mathbf{x}))\varphi (\mathbf{x}),&n(b)\in N(\mathbb{A}),\\ \omega (w)\varphi (\mathbf{x})&= \widehat{\varphi }(\mathbf{x}),&w=\left(\begin{smallmatrix} 0 & 1_n\\ -1_n & 0\\ \end{smallmatrix}\right),\\ \omega (h)\varphi (\mathbf{x})&=\varphi (h^{-1}\cdot \mathbf{x}),& h\in H(\mathbb{A}). \end{aligned} {\tag{2.2.1.1}} \end{align}$$

Here

•

$\chi _V: \mathbb{A}^\times /\mathbb{Q}^\times \rightarrow \mathbb{C}^\times$ is the quadratic character that corresponds to the quadratic extension $\mathbb{Q}(\sqrt {{\mathrm{disc}}(V)})/\mathbb{Q}$, and ${\mathrm{disc}}(V)\in \mathbb{Q}^\times /(\mathbb{Q}^\times )^2$ is the discriminant of $V$ defined to be$$\begin{equation*} {\mathrm{disc}}(V)\coloneq (-1)^{m \choose 2}\det ((x_i,x_j))_{i,j=1}^m \end{equation*}$$

for any $\mathbb{Q}$-basis $\{x_1,\ldots ,x_m\}$ of $V$.

•

$|\cdot |: \mathbb{A}^\times \rightarrow \mathbb{R}_{>0}$ is the normalized absolute value,

•

$(\mathbf{x},\mathbf{y})\coloneq ((x_i,y_j))_{i,j=1}^n\in \mathrm{Mat}_n(\mathbb{A})$ for $\mathbf{x}=(x_1,\ldots ,x_n)\in V(\mathbb{A})^n$ and $\mathbf{y}=(y_1,\ldots , y_n)\in V(\mathbb{A})^n$,

•

$\widehat{\varphi }$ is the Fourier transform of $\varphi$ using the self-dual Haar measure on $V(\mathbb{A})^n$ with respect to $\psi$,$$\begin{equation*} \widehat{\varphi }(\mathbf{x})\coloneq \int _{V(\mathbb{A})^n}\varphi (\mathbf{y})\psi (\operatorname {tr}(\mathbf{x},\mathbf{y}))\mathrm{d}\mathbf{y}. \end{equation*}$$

Remark 2.2.2.

There are many nice expository articles on the Weil representation; see e.g., Kud96Pra93Pra98. See also a more general formula of the Weil representation in Kud96, Proposition 4.3).

Example 2.2.3.

When $n=1$ we have $G={\mathrm{Sp}}(2)={\mathrm{SL}}(2)$. The standard Siegel parabolic is the standard Borel subgroup of ${\mathrm{SL}}(2)$ consisting of upper triangle matrices $P=\left\{\left(\begin{smallmatrix} a & b \\ 0 & a^{-1}\end{smallmatrix}\right)\right\}$, and $M=\left\{\left(\begin{smallmatrix} a & 0 \\ 0 & a^{-1}\end{smallmatrix}\right)\right\}$, $N=\left\{\left(\begin{smallmatrix} 1 & b \\ 0 & 1\end{smallmatrix}\right)\right\}$ are the diagonal and upper unipotent matrices, respectively. In this case the first two formulas in 2.2.1.1 simplify to

$$\begin{equation} \omega \left(\begin{smallmatrix} a & 0 \\ 0 & a^{-1}\end{smallmatrix}\right)\varphi (x)=\chi _V(a)|a|^{m/2}\varphi (xa),\quad \omega \left(\begin{smallmatrix} 1 & b \\ 0 & 1\end{smallmatrix}\right)\varphi (x)=\psi (\tfrac{1}{2}b(x,x))\varphi (x), {\tag{2.2.3.1}} \end{equation}$$

for any $a\in \mathbb{A}^\times$, $b\in \mathbb{A}$, $\varphi \in \mathscr{S} (V(\mathbb{A}))$ and $x\in V(\mathbb{A})$.

Our next goal is to use the Weil representation $\omega$ to construct theta series and Siegel Eisenstein series, starting from any common choice of a Schwartz function $\varphi \in \mathscr{S} (V(\mathbb{A})^n)$.

Definition 2.2.4.

Associated to $\varphi \in \mathscr{S} (V(\mathbb{A})^n)$, define the (two-variable) theta function

$$\begin{equation*} \theta (g,h,\varphi )\coloneq \sum _{\mathbf{x}\in V^n}\omega (g)\varphi (h^{-1}\mathbf{x}),\quad g\in G(\mathbb{A}), h\in H(\mathbb{A}). \end{equation*}$$

Then $\theta (g,h,\varphi )$ is invariant under $G(\mathbb{Q})\times H(\mathbb{Q})$ (automorphic on both $G$ and $H$ in a broad sense).

Remark 2.2.5.

Using $\theta (g,h,\varphi )$ as an integral kernel allows one to lift automorphic forms on $G$ to automorphic forms on $H$ (and vice versa): for a cuspidal automorphic representation $\pi$ of $G(\mathbb{A})$ and $\phi \in \pi$, define the theta lift $\theta _\varphi (\phi )$ of $\phi$ to $H(\mathbb{A})$ by the Petersson inner product on $G(\mathbb{Q})\backslash G(\mathbb{A})$,

$$\begin{equation*} \theta _\varphi (\phi )(h)\coloneq \langle \theta (-,h,\varphi ),\phi \rangle _\mathrm{Pet}=\int _{G(\mathbb{Q})\backslash G(\mathbb{A})}\theta (g,h,\varphi )\overline{\phi (g)}\mathrm{d}g. \end{equation*}$$

Then $\theta _\varphi (\phi )$ is an automorphic form on $H(\mathbb{A})$. This may be viewed as the starting point of the modern theory of theta correspondence, which is indispensable in the study of automorphic forms and the Langlands correspondence. We refer to Gan Gan14 for an excellent recent survey on theta correspondence.

Example 2.2.6.

Assume that $V$ is positive definite, then the theta integral

$$\begin{equation} \int _{H(\mathbb{Q})\backslash H(\mathbb{A})}\theta (g,h,\varphi )\ \mathrm{d}h {\tag{2.2.6.1}} \end{equation}$$

(or in other words, the theta lift of the constant function $\mathbf{1}$ on $H(\mathbb{A})$ to $G(\mathbb{A})$) is closely related to Siegel’s theta series. More precisely, for $\Lambda \subseteq V$ a lattice over $\mathbb{Z}$, we take the Schwartz function $\varphi =(\otimes _p\varphi _p)\otimes \varphi _\infty \in \mathscr{S} (V(\mathbb{A})^n)$ such that

•

$\varphi _p\in \mathscr{S} (V(\mathbb{Q}_p)^n)$ is the characteristic function of $(\Lambda \otimes \mathbb{Z}_p)^n$,

•

$\varphi _\infty \in \mathscr{S} (V(\mathbb{R})^n)$ is the standard Gaussian function $\varphi _\infty (\mathbf{x})=e^{-\pi \operatorname {tr}(\mathbf{x},\mathbf{x})}$.

For $\tau =\mathsf{x}+i\mathsf{y}\in \mathcal{H}_n$, we consider $g_\tau =n(\mathsf{x})m(a)\in G(\mathbb{R})$, where $a\in \mathrm{GL}_n(\mathbb{R})$ such that $\mathsf{y}=a{}^{t} a$. Then $g_\tau \cdot i 1_n=\tau$. By 2.2.1.1 we have

$$\begin{equation} \omega _\infty (g_\tau )\varphi _\infty (\mathbf{x})=\chi _\infty (\det a) |\det a|^{m/2}\cdot q^{\frac{1}{2}(\mathbf{x},\mathbf{x})}. {\tag{2.2.6.2}} \end{equation}$$

Define the classical theta integral

$$\begin{equation} \chi _\infty (\det a)^{-1}|\det a|^{-m/2}\cdot \int _{H(\mathbb{Q})\backslash H(\mathbb{A})}\theta (g_{\tau },h,\varphi )\ \mathrm{d}h. {\tag{2.2.6.3}} \end{equation}$$

Then it recovers the weighted average of theta series in 2.1.3.1 (see e.g., Han13, §4.6, KR14, §7).

In fact, let $K\subseteq H(\mathbb{A})$ be the stabilizer of $\Lambda$, then we have a bijection

$$\begin{equation*} H(\mathbb{Q})\backslash H(\mathbb{A})/K {\overset{\sim }{\longrightarrow }} \mathrm{Gen}(\Lambda ),\quad h \mapsto \Lambda '\coloneq h(\Lambda \otimes \widehat{\mathbb{Z}})\cap V. \end{equation*}$$

Let $\{h_i\}$ be a complete set of representatives of $H(\mathbb{Q})\backslash H(\mathbb{A})/K$ and let $\{\Lambda _i\}$ be the corresponding representatives of $\mathrm{Gen}(\Lambda )$ under this bijection. Then

$$\begin{equation*} \int _{H(\mathbb{Q})\backslash H(\mathbb{A})}\theta (g_{\tau },h,\varphi )\ \mathrm{d}h=\sum _i\int _{H(\mathbb{Q})\backslash H(\mathbb{Q})h_i K}\theta (g_{\tau },h,\varphi )\mathrm{d}h. \end{equation*}$$

Using $H(\mathbb{Q})\cap h_i Kh_i^{-1}=\operatorname {Aut}(\Lambda _i)$, each summand evaluates to

$$\begin{equation*} \int _{H(\mathbb{Q})\backslash H(\mathbb{Q})h_i Kh_i^{-1}} \theta (g_{\tau }, h h_i,\varphi )\mathrm{d}h=\frac{1}{\#\operatorname {Aut}(\Lambda _i)}\int _{h_iKh_i^{-1}}\theta (g_{\tau },h h_i, \varphi )\mathrm{d}h. \end{equation*}$$

Unfolding the definition, the second integral equals

$$\begin{equation*} \int _{K}\theta (g_{\tau }, h_i h, \varphi )\mathrm{d}h=\int _K\sum _{\mathbf{x}\in V^n}\omega _\infty (g_{\tau })\varphi (h^{-1}h_i^{-1}\mathbf{x})\mathrm{d}h, \end{equation*}$$

which by our choice of $\{\varphi _p\}$ evaluates to

$$\begin{equation*} \operatorname {vol}(K)\sum _{\mathbf{x}\in \Lambda _i^n}\omega _\infty (g_{\tau })\varphi _\infty (\mathbf{x}). \end{equation*}$$

Thus, combining with 2.2.6.2, we know that the classical theta integral 2.2.6.3 evaluates to

$$\begin{equation*} \operatorname {vol}(K)\sum _i \frac{1}{\#\operatorname {Aut}(\Lambda _i)}\cdot \sum _{\mathbf{x}\in \Lambda _i^n}q^{\frac{1}{2}(\mathbf{x},\mathbf{x})}=\operatorname {vol}(K)\sum _i \frac{1}{\#\operatorname {Aut}(\Lambda _i)}\cdot \theta _{\Lambda _i}(\tau ). \end{equation*}$$

Finally notice that

$$\begin{equation*} \operatorname {vol}(K)\sum _i \frac{1}{\#\operatorname {Aut}(\Lambda _i)}=\operatorname {vol}(H(\mathbb{Q})\backslash H(\mathbb{A})), \end{equation*}$$

thus if we normalize the Haar measure $\mathrm{d}h$ such that $\operatorname {vol}(H(\mathbb{Q})\backslash H(\mathbb{A}))=1$, then the classical theta integral 2.2.6.3 recovers the weighted average of theta series in 2.1.3.1.

Definition 2.2.7.

Associated to $\varphi \in \mathscr{S} (V(\mathbb{A})^n)$, also define the Siegel Eisenstein series

$$\begin{equation*} E(g, s,\varphi )\coloneq \sum _{\gamma \in P(\mathbb{Q})\backslash G(\mathbb{Q})}\Phi _\varphi (\gamma g,s),\quad g\in G(\mathbb{A}), s\in \mathbb{C}, \end{equation*}$$

where

$$\begin{equation} \mathscr{S} (V(\mathbb{A})^n)\rightarrow {\mathrm{Ind}}_{P(\mathbb{A})}^{G(\mathbb{A})}(\chi _V|\cdot |^{s+\frac{n+1}{2}}),\quad \varphi \mapsto \Phi _\varphi (g,s)\coloneq \omega (g)\varphi (0)\cdot |\det a(g)|^{s-s_0} {\tag{2.2.7.1}} \end{equation}$$

is the standard Siegel–Weil section of the degenerate principal series representation of $G(\mathbb{A})$ and

$$\begin{equation*} s_0\coloneq \frac{m-(n+1)}{2}. \end{equation*}$$

Here we write $g=nm(a)k$ under the Iwasawa decomposition $G(\mathbb{A})=N(\mathbb{A})M(\mathbb{A})K$ for $K$ the standard maximal open compact subgroup of $G(\mathbb{A})$, and the quantity $|\det a(g)|\coloneq |\det a|$ is well defined.

Example 2.2.8.

Similarly, define the classical Siegel Eisenstein series

$$\begin{equation*} E(\tau ,s,\varphi )\coloneq \chi _\infty (\det a)^{-1}|\det a|^{-m/2}\cdot E(g_{\tau },s, \varphi ). \end{equation*}$$

When $V$ is positive definite and $\varphi$ is chosen as in Example 2.2.6, the special value $E(\tau , s_0,\varphi )$ at $s=s_0$ essentially recovers $E_\Lambda (\tau )$ in 2.1.3.1.

For example, consider the case $m=2$, $n=1$ and $\Lambda =\mathbb{Z}^2$ equipped with the quadratic form $Q=x^2+y^2$. We have ${\mathrm{disc}}(V)=-1$ and the quadratic character $\chi _V$ corresponds to the quadratic extension $\mathbb{Q}(i)/\mathbb{Q}$, and hence corresponds to the Dirichlet character $\chi : (\mathbb{Z} /4 \mathbb{Z})^\times \simeq \{\pm 1\}$. For $\tau =\mathsf{x}+i\mathsf{y}\in \mathcal{H}$, we have

$$\begin{equation*} g_\tau =\left(\begin{smallmatrix} 1 & \mathsf{x}\\ 0 & 1\end{smallmatrix}\right)\left(\begin{smallmatrix} \mathsf{y}^{1/2} & 0 \\ 0 & \mathsf{y}^{-1/2}\end{smallmatrix}\right)\in {\mathrm{SL}}_2(\mathbb{R}). \end{equation*}$$

For $\gamma = \left(\begin{smallmatrix} a & b\\ c & d\end{smallmatrix}\right)\in {\mathrm{SL}}_2(\mathbb{Q})$, one can compute that the Siegel–Weil section evaluates to

$$\begin{equation*} \Phi _\varphi (\gamma g_\tau ,s)=\frac{\chi (d)\cdot \mathsf{y}^{1/2}}{c\tau +d}\cdot \operatorname {Im}(\gamma \tau )^{s/2}=\frac{\chi (d)}{c\tau +d}\cdot \frac{\mathsf{y}^{(1+s)/2}}{|c\tau +d|^s}, \end{equation*}$$

using 2.2.6.2 together with 2.2.3.1 (or the more general Kud96, Proposition 4.3). Comparing with 1.5.1.3, we see that the classical Eisenstein series

$$\begin{equation*} E(\tau , 0, \varphi )=\mathsf{y}^{-1/2}\cdot E(g_\tau , 0,\varphi ) \end{equation*}$$

at $s_0=0$ essentially recovers the Eisenstein series $E_1^\chi (\tau )$ of weight 1 (up to a nonzero constant).

The Siegel Eisenstein series $E(g,s,\varphi )$ converges absolutely when ${\mathrm{Re}}(s)>\frac{n+1}{2}$. It has a meromorphic continuation to $s\in \mathbb{C}$ and satisfies a functional equation relating $s\leftrightarrow -s$ (i.e., centered at $s=0$). The Siegel–Weil formula gives a precise identity of the form

theta integral $\longleftrightarrow$ value of Siegel Eisenstein series at $s=s_0$.

Notice that $s=s_0$ is the unique point such that the map $\varphi \mapsto \Phi _\varphi (g,s)$ in 2.2.7.1 is $G(\mathbb{A})$-equivariant, so that both sides of the identity at least have the same transformation behavior with respect to the Weil representation.

Theorem 2.2.9 (Siegel–Weil formula Wei65KR88aKR88b).

Let $\alpha$ be the dimension of a maximal isotropic subspace of $V$. If $\alpha =0$ (i.e., $V$ is anisotropic) or $\alpha >0$ and $m-\alpha >n+1$, then $E(g,s,\varphi )$ is holomorphic at $s_0$ and

$$\begin{equation*} \kappa \cdot \int _{H(\mathbb{Q})\backslash H(\mathbb{A})}\theta (g,h,\varphi )\ \mathrm{d}h= E(g,s_0,\varphi ) , \end{equation*}$$

where $\kappa =1$ if $m>n+1$ or $\kappa =2$ otherwise. Here the Haar measure $\mathrm{d}h$ is normalized so that $\operatorname {vol}(H(\mathbb{Q})\backslash H(\mathbb{A}))=1$.

Example 2.2.10.

When $V$ is positive definite (so $\alpha =0$), the Siegel–Weil formula recovers Siegel’s formula (see Examples 2.2.8 and 2.2.6), and in particular recovers Jacobi’s formula as a special case of $m=2$, $n=1$, $(G,H)=({\mathrm{Sp}}(2),\mathrm{O}(2,0))$ (see Example 2.1.4).

Remark 2.2.11.

The condition in Theorem 2.2.9 is known as Weil’s convergence condition, which ensures the convergence of the theta integral. It has been a long effort, starting with the work of Kudla–Rallis KR94, to generalize the Siegel–Weil formula outside the convergence range and for all reductive dual pairs of classical groups. We refer to Gan–Qiu–Takeda GQT14 for the most general Siegel–Weil formula and a nice summary of its literature and history.

Remark 2.2.12.

When $m$ is odd, the theta series and Eisenstein series have half integral weights $m/2$, and are automorphic forms on the metaplectic cover $\tilde{G}={\mathrm{Mp}}(2n)$ of ${\mathrm{Sp}}(2n)$. In this case the Weil representation needs to be modified to be a representation of $\tilde{G}(\mathbb{A})\times H(\mathbb{A})$ and the Siegel–Weil formula still holds after modification.

3. Geometric Siegel–Weil formula

In this section we discuss an example of the Siegel–Weil formula in the indefinite case $(G,H)=({\mathrm{Sp}}(4),\mathrm{O}(2,2))$ originating from the classical work of Hurwitz, and we use it to motivate the more general geometric Siegel–Weil formula.

3.1. Hurwitz class number relation

Definition 3.1.1.

For any positive integer $D$, the Hurwitz class number $H(D)$ is defined to be the weighted size of ${\mathrm{SL}}_2(\mathbb{Z})$-equivalence classes of positive definite binary quadratic forms

$$\begin{equation*} ax^2+bxy+cy^2 \text{ with discriminant } b^2-4ac=-D, \quad a,b,c\in \mathbb{Z}. \end{equation*}$$

Here the forms equivalent to $a(x^2+y^2)$ ($D=4a^2$) and $a(x^2+xy+y^2)$ ($D=3a^2$) are counted with multiplicities $1/2$ and $1/3$, respectively, due to extra symmetry.

Example 3.1.2.

$$\begin{align*} H(3) =\frac{1}{3}& \longleftrightarrow \{x^2+xy+y^2\},\\ H(8) = 1 & \longleftrightarrow \{x^2+2y^2\},\\ H(11) = 1& \longleftrightarrow \{x^2+xy+3y^2\},\\ H(12) =\frac{4}{3}& \longleftrightarrow \{2(x^2+xy+y^2), x^2+3y^2\},\\ H(20)= 2 &\longleftrightarrow \{x^2+5y^2, 2x^2+2xy+3y^2\}.\\ \end{align*}$$

Example 3.1.3.

Table 1 lists the first few Hurwitz class numbers.

Example 3.1.4.

When $-D$ is a fundamental discriminant and $D>4$, the Hurwitz class number $H(D)$ is equal to the class number of the imaginary quadratic field $\mathbb{Q}(\sqrt {-D})$.

Understanding these class numbers $H(D)$ remains a central subject in algebraic number theory. The following remarkable formula, which we call the Hurwitz class number relation or the Hurwitz formula, gives an elementary expression for a certain sum of Hurwitz class numbers.

Theorem 3.1.5 (Kronecker Kro60, Gierster Gie83, Hurwitz Hur85).

If $m$ is not a perfect square, then

$$\begin{equation} \sum _{d d'=m} \max \{ d, d'\}=\sum _{t\in \mathbb{Z}, 4m-t^2>0}H(4m-t^2). {\tag{3.1.5.1}} \end{equation}$$

Example 3.1.6.

When $m=3$, the Hurwitz class number relation says

$$\begin{align*} 3+3 & = H(3)+H(8)+H(11)+H(12)+H(11)+H(8)+H(3) \\ & =\frac{1}{3}+1+1+\frac{4}{3}+1+1+\frac{1}{3}. \end{align*}$$

When $m=5$, the Hurwitz class number relation says

$$\begin{align*} 5 + 5 & = H(4)+H(11)+H(16)+H(19)+H(20)+H(19)+H(16)+H(11)+H(4) \\ & =\frac{1}{2}+1+\frac{3}{2}+1+2+1+\frac{3}{2}+1+\frac{1}{2}, \end{align*}$$

A quite nontrivial way to decompose the integers 6 and 10, respectively!

Hurwitz Hur85 proved this formula using the modular $j$-invariant and the modular polynomial $\Phi _m(x, y)\in \mathbb{Z}[x, y]$ of level $m$ (which defines the modular curve $Y_0(m)$ over $\mathbb{Z}$). He recognized the left-hand side (LHS) of 3.1.5.1 as $\deg \Phi _m(x,x)$ and computed this degree in a different way involving Hurwitz class numbers to arrive at formula 3.1.5.1.

3.2. A geometric proof

(cf. GK93). From the modern point of view, we have a nice geometric interpretation of Hurwitz’s proof, in terms of the geometry of the modular curve

$$\begin{equation*} Y(\mathbb{C})={\mathrm{SL}}_2(\mathbb{Z})\backslash \mathcal{H}. \end{equation*}$$

The modular curve $Y(\mathbb{C})$ is the moduli space of elliptic curves:

$$\begin{align*} Y&=\{E: \text{elliptic curve over }\mathbb{C} \text{ (up to isomorphism)}\}\\ \tau &\mapsto E_\tau =\mathbb{C}/(\mathbb{Z}+\mathbb{Z}\tau ), \end{align*}$$

which allows one to define a canonical model of $Y$ as an algebraic curve over $\mathbb{Q}$. Each elliptic curve $E_\tau$ has a Weierstrass equation

$$\begin{equation*} E_\tau : y^2=x^3+A_\tau x+B_\tau . \end{equation*}$$

The $j$-invariant

$$\begin{equation*} j(E_\tau )\coloneq 1728\cdot \frac{4A_\tau ^3}{4A_\tau ^3+27B_\tau ^2} \end{equation*}$$

only depends on the isomorphism class of $E_\tau$ and gives rise to an isomorphism

$$\begin{equation} j: Y\xrightarrow {\sim } \mathbb{A}^1, \quad \tau \mapsto j(E_\tau ). {\tag{3.2.0.1}} \end{equation}$$

Definition 3.2.1.

Define the surface $X$ to be the product of two modular curves,

$$\begin{equation*} X\coloneq Y\times Y=\{(E,E')\}, \end{equation*}$$

which is the moduli space of pairs of elliptic curves $(E,E')$. For each positive integer $m$, we define the modular correspondence $Z(m)$ over the surface $X$ by

$$\begin{equation*} Z(m)\coloneq \{ (E, E',\varphi ): \phi : E\xrightarrow {\deg m} E'\}\rightarrow X, \end{equation*}$$

parameterizing a pair of elliptic curves $(E,E')$ together with a degree $m$ isogeny $\phi : E\rightarrow E'$.

The isogeny $\phi$ imposes one nontrivial condition, and thus $Z(m)$ defines a divisor on $X$. For example, when $m=1$ the modular correspondence is nothing but the diagonally embedded modular curve

$$\begin{equation*} Z(1)=\Delta (Y)\subseteq X=Y\times Y. \end{equation*}$$

Given two divisors $Z(m)$ and $Z(1)$ on the surface $X$, one expects that the intersection $Z(m)\cap Z(1)$ should be zero dimensional. When this is the case (i.e., when $Z(m)$ and $Z(1)$ intersect properly), we obtain a geometric intersection number $\langle Z(m),Z(1)\rangle _X$ by counting the number of intersection points weighted by intersection multiplicities (see Figure 2).

A curious observation comes: the geometric intersection number $\langle Z(m),Z(1)\rangle _X$ is equal to the LHS of 3.1.5.1,

$$\begin{equation} \langle Z(m),Z(1)\rangle _X=\sum _{dd'=m}\max \{d,d'\}. {\tag{3.2.1.1}} \end{equation}$$

In fact, under the isomorphism 3.2.0.1 we know that $X$ has a natural compactification

$$\begin{equation*} \overline{X}=\mathbb{P}^1\times \mathbb{P}^1. \end{equation*}$$

Since $\overline{Z(1)}=\Delta (\mathbb{P}^1)$, the well-known cohomological equivalence on $\mathbb{P}^1\times \mathbb{P}^1$,

$$\begin{equation*} \Delta (\mathbb{P}^1)\sim \mathbb{P}^1\times \{\text{pt}\}+ \{\text{pt}\}\times \mathbb{P}^1 \end{equation*}$$

then implies that

$$\begin{equation*} \langle \overline{Z(m)},\overline{Z(1)}\rangle _{\overline{X}}=\langle \overline{Z(m)}, \mathbb{P}^1\times \{\text{pt}\}\rangle _{\overline{X}}+\langle \overline{Z(m)}, \{\text{pt}\}\times \mathbb{P}^1\rangle _{\overline{X}}=2\sum _{d|m}d, \end{equation*}$$

where the last equality comes from counting the number of degree $m$ isogenies with a fixed source (resp., target) elliptic curve. The desired identity 3.2.1.1 then follows from subtracting the contribution at $\{\infty \}=\mathbb{P}^1\setminus \mathbb{A}^1$.

On the other hand, via the moduli interpretation we have

$$\begin{equation*} Z(m)\cap Z(1)=\{(E,E): E\xrightarrow {\deg m} E\}. \end{equation*}$$

When $m$ is a perfect square, we know that $Z(m)$ contains $Z(1)$ (by considering the multiplication-by-$\sqrt {m}$ isogeny $E\xrightarrow {[\sqrt {m}]} E$), and thus $Z(m)$ and $Z(1)$ do not intersect properly. However, when $m$ is not a perfect square (as assumed in Theorem 3.1.5), $E$ has to be an elliptic curve with complex multiplication by an imaginary quadratic order $\mathbb{Z}[\alpha ]$, where $\alpha \overline{\alpha }=m$, and thus $Z(m)$ and $Z(1)$ do intersect properly. Using the theory of complex multiplication, counting the weighted number of such elliptic curves with complex multiplication exactly gives the sum of Hurwitz class numbers as the right-hand side (RHS) of 3.1.5.1,

$$\begin{equation} \langle Z(m),Z(1)\rangle _X=\sum _{t\in \mathbb{Z}, 4m-t^2>0}H(4m-t^2). {\tag{3.2.1.2}} \end{equation}$$

Combining 3.2.1.1 and 3.2.1.2 completes our sketch of the geometric proof of the Hurwitz formula (see Gross–Keating GK93 for complete details).

3.3. Hurwitz formula as a geometric Siegel–Weil formula

The Hurwitz class numbers appearing in the Hurwitz formula 3.1.5.1 also naturally appear as Fourier coefficients of Siegel Eisenstein series. More precisely, consider the Siegel Eisenstein series $E(\tau ,s)$ on ${\mathrm{Sp}}(4)$ of weight 2,

$$\begin{align*} E(\tau ,s)& \coloneq \sum _{\left(\begin{smallmatrix} A & B\\ C & D\end{smallmatrix}\right)\in P(\mathbb{Z})\backslash {\mathrm{Sp}}_4(\mathbb{Z})}\det (C\tau +D)^{-2}\frac{\det \operatorname {Im}(\tau )^{s-1/2}}{|\det (C\tau +D)|^{2s-1}}\\ & =\sum _{T\in \operatorname {Sym}_2(\mathbb{Z})_{>0}}E_T(\tau ,s)q^T+\cdots . \end{align*}$$

Then we have (up to a normalizing constant)

$$\begin{equation*} H(4m-t^2)=E_T(\tau ,s_0=1/2), \quad T=\left(\begin{smallmatrix} m & t/2 \\ t/2 & 1\end{smallmatrix}\right). \end{equation*}$$

Notice that the condition $4m-t^2>0$ neatly translates to the condition $T>0$.

We summarize our discussion with the following diagram.

$$\begin{equation} \vcenter{\img[][317pt][125pt][{$\xymatrix@R=.5em{\displaystyle\sum_{d d'=m} \max\{ d, d'\} \ar@{=}[d] &amp; = &amp; \displaystyle\sum_{t\in\mathbb{Z}, 4m-t^2&gt;0}H(4m-t^2) \ar@{=}[d]\\ \langle Z(m), Z(1)\rangle_{X} &amp; = &amp; \displaystyle\sum_{T=\left(\begin{smallmatrix}m &amp; t/2 \\ t/2 &amp; 1\end{smallmatrix}\right)&gt;0}E_T(\sz,1/2) \\ &amp;&amp;\\ &amp;&amp;\\&amp;&amp;\\ {\begin{array}[h]{c} \text{intersection number of} \\ \text{two divisors on } X \end{array}} \ar@{&lt;~&gt;}[uuuu] &amp; = &amp; {\begin{array}[h]{c} \text{sum of Fourier coefficients of}\\ \text{a Siegel Eisenstein series on } \Sp(4). \end{array}} \ar@{&lt;~&gt;}[uuuu]}$}]{Images/img3b21a4212a99a8d381b6514a22090bb4.svg}} {\tag{3.3.0.1}} \end{equation}$$

In this way the Hurwitz formula can be viewed as a geometric Siegel–Weil formula for the pair $(G, H)=({\mathrm{Sp}}(4), \mathrm{O}(2,2))$, where one replaces the theta integral on $\mathrm{O}(2,2)$ by a geometric theta series, i.e., the generating series of geometric intersection numbers of modular correspondences for the surface $X=Y\times Y$,

geometric theta series on $X\longleftrightarrow$ value of Siegel Eisenstein series on ${\mathrm{Sp}}(4)$ at $s_0=1/2$.

Notice here the natural appearance of the product of modular curves due to the exceptional isomorphism ${\mathrm{SO}}(2,2)\simeq {\mathrm{SL}}_2\times {\mathrm{SL}}_2/\{\pm 1\}$. This geometric Siegel–Weil formula further computes a more general geometric intersection number $\langle Z(m), Z(n)\rangle _{X}$ as the sum of $E_T(\tau ,1/2)$ with $T=\left(\begin{smallmatrix} m & t/2 \\ t/2 & n\end{smallmatrix}\right)>0$.

Remark 3.3.1.

The remarkable discovery that generating series involving intersection numbers of cycles are modular originates from the work of Hirzebruch–Zagier HZ76 on Hilbert modular surfaces. Historically HZ76 was the primary motivation in Kudla’s work discussed below (cf. the introduction of KM90Kud97a) and also in the work of Gross–Kohnen–Zagier GKZ87. See Example 3.5.1 and Remark 3.5.5(i), (ii).

3.4. Orthogonal Shimura varieties

Kudla proved a more general geometric Siegel–Weil formula by replacing the surface $X=Y\times Y$ by an orthogonal Shimura variety of arbitrary dimension. Our next goal is to discuss Kudla’s formula.

Let $F$ be a totally real number field. Pick a real place $w$ of $F$. Let $V$ be a quadratic space over $F$ of dimension $m$ such that for any place $v|\infty$ of $F$,

$$\begin{equation*} \text{ the $F_v$-quadratic space } V_v\coloneq V \otimes _F F_v \text{ has signature } \begin{cases} (m-2,2), & \text{if $v=w$}, \\ (m,0), & \text{if $v\ne w$}. \end{cases} \end{equation*}$$

Let $G=\mathrm{GSpin}(V)$, which sits in an exact sequence

$$\begin{equation*} 1\rightarrow \mathbb{G}_m\rightarrow G\rightarrow {\mathrm{SO}}(V)\rightarrow 1. \end{equation*}$$

Associated to any open compact subgroup $K\subseteq G(\mathbb{A}_F^\infty )$, we have a GSpin Shimura variety $X=\mathrm{Sh}_G$, which has a smooth canonical model of dimension $m-2$ over the reflex field $F$ (viewed as a subfield of $\mathbb{C}$ via the embedding induced by the place $w$) and admits complex uniformization

$$\begin{equation*} X(\mathbb{C})=G(F)\backslash [ \mathbb{D}\times G(\mathbb{A}_F^\infty )/K]. \end{equation*}$$

Here $\mathbb{D}$ is the hermitian symmetric domain of oriented negative 2-planes in $V_w$. Unfolding the definition we may rewrite $X(\mathbb{C})$ as a disjoint union of quotients of $\mathbb{D}$ by congruence subgroups of $G(F)$ (see Kud04, §1). The Shimura variety $X$ is quasi-projective, and it is projective when $V$ is anisotropic (e.g., when $F\ne \mathbb{Q}$, by the signature condition).

Remark 3.4.1.

One technical reason that one would prefer to work with the Shimura variety $X$ associated to $G$ (instead of ${\mathrm{SO}}(V)$) is that it is of Hodge type (instead of abelian type), and it admits an embedding into a Siegel modular variety (the moduli space of polarized abelian varieties) of larger dimension.

Example 3.4.2 (cf. Kud04HP14FH00).

Consider $F=\mathbb{Q}$. Via accidental isomorphisms between $G=\mathrm{GSpin}(m-2,2)$ and classical groups of symplectic type in low ranks, the Shimura varieties $X$ recover many classical modular varieties in low dimensions (see Table 2, where $D$ is a quaternion algebra over $\mathbb{Q}$ and $E=\mathbb{Q}(\sqrt {d})$ is a real quadratic field). When $\dim X\le 19$, $X$ is also closely related to the moduli space of polarized K3 surfaces and has proved to be useful for studying the arithmetic of K3 surfaces.

3.5. Kudla’s generating series of special cycles and the modularity conjecture

The Shimura variety $X$ is equipped with special divisors $Z(m)\rightarrow X$ generalizing in the case $X=Y\times Y$. Via an embedding into a Siegel modular variety, $Z(m)$ parameterizes certain polarized abelian varieties together with a special endomorphism of degree $m$ (see MP16, §5).

Example 3.5.1.

The special divisor $Z(m)$ recovers

•

Heegner points on modular curves and Shimura curves when $\dim X=1$ (cf. Kud04, Appendix),

•

modular correspondences $Z(m)$ on $X=Y\times Y$ considered in §3.2 when $\dim X=2$ and $G$ is split,

•

Hirzebruch–Zagier cycles HZ76 on Hilbert modular surfaces when $\dim X=2$ and $G$ is nonsplit.

More generally, for any $y\in V$ with $(y,y)>0$, its orthogonal complement $V_y\subseteq V$ has rank $n-1$. The embedding $G_y\coloneq \mathrm{GSpin}(V_y)\hookrightarrow G=\mathrm{GSpin}(V)$ defines a Shimura subvariety of codimension 1,

$$\begin{equation*} \mathrm{Sh}_{G_y}\rightarrow X=\mathrm{Sh}_{G}. \end{equation*}$$

For any $x\in V(\mathbb{A}_F^\infty )$ with $(x,x)\in F_{>0}$, there exists $y\in V$ and $g\in G(\mathbb{A}_F^\infty )$ such that $y=gx$. Define the special divisor

$$\begin{equation*} Z(x)\rightarrow X \end{equation*}$$

to be the $g$-translate of $\mathrm{Sh}_{G_y}$. For any $\mathbf{x}=(x_1,\ldots , x_n)\in V(\mathbb{A}_F^\infty )^n$ with $(\mathbf{x},\mathbf{x})\in \operatorname {Sym}_{n}(F)_{>0}$, define the special cycle (of codimension $n$)

$$\begin{equation*} Z(\mathbf{x})=Z(x_1)\cap \cdots \cap Z(x_n)\rightarrow X. \end{equation*}$$

Here $\cap$ denotes the fiber product over $X$. More generally, for a $K$-invariant Schwartz function $\varphi \in \mathscr{S} (V(\mathbb{A}_F^\infty )^n)^K$ and $T\in \operatorname {Sym}_{n}(F)_{>0}$, define the weighted special cycle

$$\begin{equation*} Z_\varphi (T)=\sum _{\mathbf{x}\in K\backslash V(\mathbb{A}_F^\infty )^n\atop (\mathbf{x},\mathbf{x})=T}\varphi (\mathbf{x}) Z(\mathbf{x})\in {\mathrm{CH}}^n(X)_\mathbb{C}\coloneq {\mathrm{CH}}^n(X)\otimes \mathbb{C}. \end{equation*}$$

Here ${\mathrm{CH}}^n(X)$ is the Chow group of algebraic cycles of codimension $n$ on $X$ (up to rational equivalence). With extra care, we can also define $Z_\varphi (T)\in {\mathrm{CH}}^n(X)_\mathbb{C}$ for any $T\in \operatorname {Sym}_{n}(F)_{\ge 0}$ (see Kud04YZZ09).

Definition 3.5.2.

Define Kudla’s generating series of special cycles,

$$\begin{equation} Z_\varphi (\tau )=\sum _{T\in \operatorname {Sym}_n(F)_{\ge 0}}Z_\varphi (T)q^T, {\tag{3.5.2.1}} \end{equation}$$

as a formal sum valued in ${\mathrm{CH}}^n(X)_\mathbb{C}$, where

$$\begin{equation*} \tau \in \mathcal{H}_n=\{\mathsf{x}+i\mathsf{y}: \mathsf{x}\in \operatorname {Sym}_n(F_\infty ),\ \mathsf{y}\in \operatorname {Sym}_n(F_\infty )_{>0}\},\quad q^T\coloneq \prod _{v|\infty }e^{2\pi i \operatorname {tr}T\tau _v}. \end{equation*}$$

Remark 3.5.3.

Analogous constructions of special divisors and Kudla’s generating series also apply to Shimura varieties of unitary type associated to hermitian spaces with signature $(m-1,1)$ at one archimedean place and signature $(m,0)$ at all other archimedean places (see Liu11a). These Shimura varieties of orthogonal/unitary type can be naturally viewed as Shimura varieties associated to totally definite incoherent quadratic/hermitian spaces (Z19Gro20); see §6.5. In the unitary case, we obtain a generating series of the form

$$\begin{equation} Z_\varphi (\tau )=\sum _{T\in \mathrm{Herm}_n(F)_{\ge 0}}Z_\varphi (T)q^T, {\tag{3.5.3.1}} \end{equation}$$

where we replace positive semidefinite symmetric matrices $\operatorname {Sym}_n(F)_{\ge 0}$ by positive semidefinite hermitian matrices $\mathrm{Herm}_n(F)_{\ge 0}$, and we replace Siegel’s half-space by the hermitian half-space

$$\begin{equation*} \tau \in \mathcal{H}_n\coloneq \{\mathsf{x}+i\mathsf{y}: \mathsf{x}\in \mathrm{Herm}_n(F_\infty ),\ \mathsf{y}\in \mathrm{Herm}_n(F_\infty )_{>0}\}. \end{equation*}$$

We may view $Z_\varphi (\tau )$ as a geometric theta series, now valued in Chow groups for cycles of arbitrary codimension $n$. The analogy to Siegel’s theta series 2.1.1.1 and the theta integral 2.2.6.1 leads to Kudla’s modularity conjecture:

Conjecture 3.5.4 (Kudla’s modularity).

The formal generating series $Z_\varphi (\tau )$ converges absolutely and defines a modular form on $\mathcal{H}_n$ of weight $m/2$ valued in ${\mathrm{CH}}^n(X)_\mathbb{C}$.

Remark 3.5.5.

(i)

The analogous modularity in Betti cohomology, i.e., the modularity of the generating series valued in $H^{2n}(X(\mathbb{C}), \mathbb{C})$ defined by the image of $Z_\varphi (T)$ under the cycle class map$$\begin{equation*} {\mathrm{CH}}^n(X)\rightarrow H^{2n}(X(\mathbb{C}), \mathbb{Z}), \end{equation*}$$

is known by the classical work of Kudla–Millson KM90. The special case of special divisors on Hilbert modular surfaces dates back to Hirzebruch–Zagier HZ76 (see also Funke–Millson FM14).

(ii)

Kudla’s modularity conjecture was originally formulated for orthogonal Shimura varieties over $\mathbb{Q}$ (Kud97aKud04). In this case, Borcherds Bor99 proved the conjecture for the divisor case $m=1$ (the special case of Heegner points on modular curves dates back to the classical work of Gross–Kohnen–Zagier GKZ87). Zhang Zha09 proved the modularity for general $n$ assuming the absolute convergence of the series. Bruinier–Westerholt-Raum BWR15 proved the desired convergence and hence established Kudla’s modularity conjecture for orthogonal Shimura varieties over $\mathbb{Q}$. More recently, Bruinier–Zemel BZ22 have extended the modularity to toroidal compactifications of orthogonal Shimura varieties when $n=1$.

(iii)

For orthogonal Shimura varieties over totally real fields, Yuan–Zhang–Zhang YZZ09 proved the modularity for $n=1$ (see also Bruinier Bru12 for a different proof) and reduced the $n>1$ case to the convergence.

(iv)

Conjecture 3.5.4 in the unitary case was formulated by Liu Liu11a, who also proved the case $n=1$ and reduced the $n>1$ case to the convergence. Recently Xia Xia21 proved the desired convergence when $E=\mathbb{Q}(\sqrt {-d})$ for $d=1,2,3,7,11$ (in the notation of §6.5), and thus established Conjecture 3.5.4 in these cases.

(v)

Kudla Kud04, Problem 4 also proposed the modularity problem in the arithmetic Chow group ${\widehat{\mathrm{CH}}}^n(\mathcal{X})$ of a suitable (compactified) integral model $\mathcal{X}$ of $X$ (see GS90BGKK07 and also Sou92). The problem seeks to define canonically an explicit arithmetic generating series $\widehat{\mathcal{Z}}_\varphi (\tau )$ valued in ${\widehat{\mathrm{CH}}}^n(\mathcal{X})_\mathbb{C}$ which lifts $Z_\varphi (\tau )$ under the restriction map$$\begin{equation*} {\widehat{\mathrm{CH}}}^n(\mathcal{X})\rightarrow {\mathrm{CH}}^n(X), \end{equation*}$$

and such that $\widehat{\mathcal{Z}}_\varphi (\tau )$ is modular. When $n=1$, this arithmetic modularity was proved by Howard–Madapusi Pera HMP20 (orthogonal groups over $\mathbb{Q}$) and Bruinier–Howard–Kudla–Rapoport–Yang BHK$^{+}$20 (unitary groups over $\mathbb{Q}$). Several low dimensional cases were also proved:

•

Shimura/modular curves (Kudla–Rapoport–Yang KRY06, Sankaran San14, and Du–Yang DY19),

•

Hilbert modular surfaces (Bruinier–Burgos Gil–Kühn BBGK07),

•

Product of modular curves (Berndt–Kühn BK12BK12a).

We also mention the arithmetic modularity of the difference of two arithmetic theta series by Ehlen–Sankaran ES18 for $n=1$ (unitary groups over $\mathbb{Q}$), the almost arithmetic modularity by Mihatsch–Zhang MZ21, Theorem 4.3 for $n=1$ (unitary groups over totally real fields $F\ne \mathbb{Q}$), the arithmetic modularity of Fourier–Jacobi coefficients for general $n$ by Sankaran San20 (anisotropic orthogonal groups), and several striking recent works involving applications of arithmetic modularity AGHM18Zha21SSTT22.

(vi)

Finally, we mention several recent works on the modularity conjecture for more general classes of orthogonal and unitary Shimura varieties (indefinite at more than one archimedean place) by Rosu–Yott RY20, Kudla Kud21, and Maeda Mae21Mae22.

3.6. Kudla’s geometric Siegel–Weil formula

In the special case $n=\dim X$, the generating series $Z_\varphi (\tau )$ in 3.5.2.1 is valued in ${\mathrm{CH}}^{\dim X}(X)_\mathbb{C}$ (i.e., the Chow group of 0-cycles). When $X$ is projective, composing with the degree map

$$\begin{equation*} \deg :{\mathrm{CH}}^{\dim X}(X)_\mathbb{C}\rightarrow \mathbb{C}, \end{equation*}$$

we obtain a generating series $\deg Z_\varphi (\tau )$ valued in $\mathbb{C}$. Its terms encode geometric intersection numbers between special divisors on $X$ generalizing the case $X=Y\times Y$ considered in §3.3. Kudla proved the following remarkable geometric version of the Siegel–Weil formula (by analogy with Theorem 2.2.9 specialized to the case $n=\dim X=m-2$ and so $s_0=1/2$).

Theorem 3.6.1 (Kudla’s geometric Siegel–Weil formula Kud97a, Corollary 10.5).

Assume that $X$ is projective (i.e., $V$ is anisotropic). Take $n=\dim X$. Then for any $\varphi \in \mathscr{S} (V(\mathbb{A}_F^\infty )^n)^K$ the following identity holds (up to a nonzero constant depending only on choices of measures)

$$\begin{equation} \deg Z_\varphi (\tau )\doteq E(\tau , 1/2,\varphi \otimes \varphi _\infty ). {\tag{3.6.1.1}} \end{equation}$$

Here $\varphi _\infty \in \mathscr{S} (V(F_\infty )^n)$ is a certain Schwartz function constructed from the Kudla–Millson Schwartz form (KM86).

Thus Kudla’s geometric Siegel–Weil formula is a precise identity of the form

geometric theta series on $X \longleftrightarrow$ value of Siegel Eisenstein series on ${\mathrm{Sp}}(2n)$ at $s_0=1/2$.

Remark 3.6.2.

Kudla in fact proved a geometric Siegel–Weil formula for the generating series of special cycles of all dimensions, i.e., without assuming $n=\dim X$. Here and in the next section we focus on the case of 0-cycles and refer to Kudla’s excellent surveys Kud04Kud02bKud02a for the general case.

4. Arithmetic Siegel–Weil formula

In this section we discuss an arithmetic version of the Siegel–Weil formula. Parallel to the previous section, we will use an example in the case $(G,H)=({\mathrm{Sp}}(6),\mathrm{O}(2,2))$ considered by Gross–Keating to motivate the more general case.

4.1. Gross–Keating formula

Gross–Keating took the geometric point of view of the Hurwitz class number relation and found a remarkable generalization for arithmetic intersection numbers. As the moduli space of elliptic curves, the modular curve $Y$ has a canonical integral model $\mathcal{Y}$ over $\operatorname {Spec}\mathbb{Z}$ such that $\mathcal{Y}(\mathbb{C})=Y$. The integral model $\mathcal{Y}$ is an arithmetic surface fibered over the arithmetic curve $\operatorname {Spec}\mathbb{Z}$, and its fiber above $p$ is a smooth curve in characteristic $p$ (see Figure 3).

Analogously, the surface $X$ has a canonical integral model $\mathcal{X}$ over $\operatorname {Spec}\mathbb{Z}$, which is an arithmetic threefold. The modular correspondence $Z(m)\rightarrow X$ naturally extends to a divisor $\mathcal{Z}(m)\rightarrow \mathcal{X}$. Now on the arithmetic threefold $\mathcal{X}$, we need three (instead of two) divisors so that the intersection has the expected dimension 0. Define the arithmetic intersection number by

$$\begin{equation*} \langle \mathcal{Z}(m_1), \mathcal{Z}(m_2), \mathcal{Z}(m_3)\rangle _{\mathcal{X}}\coloneq \sum _p \langle \mathcal{Z}(m_1),\mathcal{Z}(m_2), \mathcal{Z}(m_3)\rangle _{\mathcal{X}_p}\cdot \log p, \end{equation*}$$

where $\langle \mathcal{Z}(m_1),\mathcal{Z}(m_2), \mathcal{Z}(m_3)\rangle _{\mathcal{X}_p}$ encodes the intersection number supported in the fiber $\mathcal{X}_p$. This definition makes sense when these three divisors intersect properly (or more generally when their intersection is supported in finitely many fibers $\mathcal{X}_p$).

Is there a formula for this arithmetic intersection number analogous to 3.3.0.1? The index set for $T$ should be $\operatorname {Sym}_3(\mathbb{Z})$ to put the three diagonal entries $m_1,m_2,m_3$ in $T$, and thus we need to look at the Siegel Eisenstein series $E(\tau ,0)$ on ${\mathrm{Sp}}(6)$ of weight 2 (instead of on ${\mathrm{Sp}}(4)$). Then the relevant special point in the Siegel–Weil formula is $s_0=0$, the central point. However, it is too naive to expect $\langle \mathcal{Z}(m_1), \mathcal{Z}(m_2), \mathcal{Z}(m_3)\rangle _{\mathcal{X}}$ to be equal to $E(\tau ,0)$: the arithmetic intersection number involves a $\mathbb{Q}$-linear combination of $\log p$’s and hence is no longer a rational number like $E_T(\tau ,1/2)$ in 3.3.0.1. Moreover, the Eisenstein series $E(\tau ,s)$ turns out to have an odd functional equation at the center $s=0$, and hence $E(\tau ,0)=0$ automatically! This automatic vanishing suggests that it would be interesting to look at its first derivative $E'(\tau ,0)$ at $s=0$.

After these two appropriate modifications—replacing ${\mathrm{Sp}}(4)$ by ${\mathrm{Sp}}(6)$ and replacing the value at $s_0=1/2$ by the central derivative at $s_0=0$—it turns out that we do have the following remarkable formula relating the arithmetic intersection numbers to $E'(\tau ,0)$ (see a nice exposition of the proof in VGW+07).

Theorem 4.1.1 (Gross–Keating GK93, Gross–Kudla–Zagier).

Assume there is no positive definite binary quadratic form representing $m_1,m_2,m_3$ simultaneously. Then (up to an explicit constant)

$$\begin{equation} \vcenter{\img[][336pt][105pt][{$\xymatrix@R=4em{\langle\mathcal{Z}(m_1), \mathcal{Z}(m_2), \mathcal{Z}(m_3)\rangle_{\mathcal{X}} &amp; \!\!\!\!\doteq\!\!\!\! &amp; \displaystyle\sum_{T=\left(\begin{smallmatrix}m_1 &amp; * &amp; *\\ * &amp; m_2 &amp; * \\ * &amp; * &amp;m_3 \end{smallmatrix}\right)&gt;0} E_T'(\sz,0)\\ {\begin{array}[h]{c} \text{arithmetic intersection number}\\ \text{of three divisors on } \mathcal{X} \end{array}} \ar@{&lt;~&gt;}[u] &amp; &amp; {\begin{array}[h]{c} \text{sum of Fourier coefficients}\\ \text{of the derivative of a}\\ \text{Siegel Eisenstein series on } \Sp(6) \end{array}} \ar@{&lt;~&gt;}[u] }$}]{Images/img350d8579bb8828efd0c6e453176b2622.svg}} {\tag{4.1.1.1}} \end{equation}$$

Remark 4.1.2.

The assumption on $m_1,m_2,m_3$ is analogous to the assumption that $m$ is not a perfect square in Theorem 3.1.5, which guarantees that the three divisors intersect properly.

4.2. Arithmetic Siegel–Weil formula

Parallel to 3.3.0.1, the Gross–Keating formula can be viewed as an arithmetic Siegel–Weil formula for the pair $(G, H)=({\mathrm{Sp}}(6), \mathrm{O}(2,2))$, where one replaces the theta integral on $\mathrm{O}(2,2)$ by a generating series of arithmetic intersection numbers on the arithmetic threefold $\mathcal{X}$,

arithmetic theta series on $\mathcal{X}\longleftrightarrow$ central derivative of Siegel Eisenstein series on ${\mathrm{Sp}}(6)$.

Of course there is nothing stopping us from considering the higher dimensional case. In fact Kudla (Kud97b) and Kudla–Rapoport (KR99KR00aKR14) proposed vast conjectural generalizations of the Gross–Keating formula by

(1)

taking $\mathcal{X}$ to be a suitable regular integral model $\mathcal{X}$ of the Shimura variety $X$ of orthogonal/unitary type considered in §3.4–3.5 associated to quadratic/hermitian spaces.

(2)

defining suitable integral models $\mathcal{Z}(m)\rightarrow \mathcal{X}$ of the special divisors $Z(m)\rightarrow X$, and more generally integral models $\mathcal{Z}(T)$ for special divisors $Z(T)$.

Remark 4.2.1.

Such a regular integral model $\mathcal{X}$ is known to exist for certain level subgroups $K$ defining $X$ (cf. LZ22a, §14 in the unitary case). To avoid technicalities, we will be vague about the choice of level structures and integral models, but see Remark 5.3.2(iii), (iv).

Now take $n=\dim \mathcal{X}$ (which equals $\dim X+1$) so that the intersection of $n$ special divisors has expected dimension 0. We have a natural decomposition

$$\begin{equation*} \mathcal{Z}(m_1)\cap \cdots \cap \mathcal{Z}(m_n)=\bigsqcup _{T=\left(\begin{smallmatrix} m_1 & * & *\\ * & \cdots & * \\ * & * &m_n \end{smallmatrix}\right)} \mathcal{Z}(T) \end{equation*}$$

index by symmetric/hermitian matrices $T$ with diagonal entries $m_1,\ldots , m_n$. Here $\cap$ denotes the fiber product over $\mathcal{X}$. When $T>0$, it turns out that $\mathcal{Z}(T)$ is supported in finitely many fibers $\mathcal{Z}(T)_p$, and we have a well-defined $T$-part of the arithmetic intersection number

$$\begin{equation*} \langle \mathcal{Z}(m_1), \dots , \mathcal{Z}(m_n)\rangle _T\coloneq \sum _p \langle \mathcal{Z}(m_1), \dots , \mathcal{Z}(m_n)\rangle _{\mathcal{Z}(T)_p}\cdot \log p. \end{equation*}$$

Now we are ready to state the conjecture on an arithmetic version of the Siegel–Weil formula (by analogy with Theorem 2.2.9 in the orthogonal case specialized to $n=m+1$ and so $s_0=0$ is the central point), which is known as the Kudla–Rapoport conjecture in the unitary case (KR14, Conjecture 11.10).

Conjecture 4.2.2 (Arithmetic Sigel–Weil formula, nonsingular part).

Take $n=\dim \mathcal{X}$. Then for any $T\in \operatorname {Sym}_n(F)_{>0}$ (resp., $T\in \mathrm{Herm}_n(F)_{>0}$) in the orthogonal (resp., unitary) case with diagonal entries $m_1,\ldots , m_n$, the following identity holds (up to a nonzero constant depending only on choices of measures).

$$\begin{equation} \vcenter{\img[][313pt][75pt][{$\xymatrix@R=1em{\langle\mathcal{Z}(m_1),\cdots, \mathcal{Z}(m_n)\rangle_T &amp; \questeq&amp; E_T'(\sz,0)\\ {\begin{array}[h]{c} \text{arithmetic intersection number}\\ \text{of $n$ divisors on } \mathcal{X} \end{array}} \ar@{&lt;~&gt;}[u] &amp; &amp; {\begin{array}[h]{c} \text{Fourier coefficient of the}\\ \text{central derivative of a}\\ \text{ Siegel Eisenstein series} \\ \text{on } \Sp(2n) \text{\textup{(}resp., $\UU(n,n)$\textup{)}} \end{array}} \ar@{&lt;~&gt;}[u] }$}]{Images/imgfa154e673d6fbe76dd82d4845c5d56ca.svg}} {\tag{4.2.2.1}} \end{equation}$$

Remark 4.2.3.

In general the $n$ special divisors do not intersect properly, and a more sophisticated definition of the arithmetic intersection numbers is needed (cf. Definition 5.1.3). In particular, with the correct definition the conjecture works even for improper intersections.

Thus the arithmetic Siegel–Weil formula is a precise conjectural identity of the form

arithmetic theta series on $\mathcal{X}\longleftrightarrow$ central derivative of Siegel Eisenstein series on ${\mathrm{Sp}}(2n)$ ($\mathrm{U}(n,n)$).

Now we can state one of the main results of LZ22a (see LZ22a, Theorem 1.3.1 for more precise technical assumptions).

Theorem 4.2.4 (Li and Zhang LZ22a).

Conjecture 4.2.2 holds for arbitrary $n$ in the unitary case.

Our recent work with Zhang LZ22b has also established a slightly weaker semi-global (at a good odd prime $p$) version of Conjecture 4.2.2 in the orthogonal case. We will discuss some key ideas of the proof in §5.

Remark 4.2.5.

(i)

When $n=3$ and $G$ is split (orthogonal case), the arithmetic Siegel–Weil formula recovers the Gross–Keating formula, and also treats the case of improper intersections (see Remarks 4.1.2, 4.2.3).

(ii)

Theorem 4.2.4 was previously proved when $n=3$ by Terstiege Ter11Ter13 and when the intersection is zero dimensional (GK93KR99KR00aKR11BY21).

(iii)

There is also an archimedean part of the arithmetic Siegel–Weil formula, relating archimedean arithmetic intersection numbers with the nonsingular but indefinite Fourier coefficients of $E'(\tau ,0)$. These Fourier coefficients are nonholomorphic, unlike the positive definite Fourier coefficients in Conjecture 4.2.2. This archimedean arithmetic Siegel–Weil formula was proved by Liu Liu11a (unitary case), and Garcia–Sankaran GS19 in full generality (see also Bruinier–Yang BY21 for an alternative proof in the orthogonal case).

(iv)

Kudla conjectured that there should also be a singular part of the arithmetic Siegel–Weil formula, relating the singular Fourier coefficients of $E'(\tau ,0)$ to certain arithmetic intersection numbers. However the singular part is more difficult to prove, or even to formulate precisely, cf. Kud04, Problem 6. As a special case, the constant term of the arithmetic Siegel–Weil formula should roughly relate the arithmetic volume of $\mathcal{X}_n$ to logarithmic derivatives of Dirichlet $L$-functions. Such an explicit arithmetic volume formula was proved by Hörmann Hör14 (orthogonal case) and Bruinier–Howard BH21 (unitary case), though a precise comparison with the constant term of $E'(\tau ,0)$ is yet to be formulated and established.

(v)

Ideally, putting all singular/nonsingular and archimedean/nonarchimedean parts together, one should arrive at a full arithmetic Siegel–Weil formula of the form in complete analogy to 3.6.1.1,$$\begin{equation*} \widehat{\deg }\ \widehat{\mathcal{Z}}_\varphi (\tau )\questeq E'(\tau ,0,\varphi \otimes \varphi _\infty ). \end{equation*}$$

Here $\widehat{\mathcal{Z}}_\varphi (\tau )$ is the conjectural arithmetic theta series in Remark 3.5.5(v), $\widehat{\deg }$ is the arithmetic degree map$$\begin{equation*} \widehat{\deg }: {\widehat{\mathrm{CH}}}^{\dim \mathcal{X}}(\mathcal{X})_\mathbb{C}\rightarrow \mathbb{C}, \end{equation*}$$

and $\varphi _\infty$ is the standard Gaussian function on the totally positive definite space over $F_\infty$. The full arithmetic Siegel–Weil formula was established by Kudla, Rapoport, and Yang (KRY99Kud97bKR00bKRY06) for $n=1,2$ (orthogonal case) in great generality. However, it remains an open problem to formulate such a precise full arithmetic Siegel–Weil formula in higher dimensions.

(vi)

Recently Feng–Yun–Zhang FYZ21 proved a higher Siegel–Weil formula over function fields for unitary groups, which relates nonsingular coefficients of the $r$th derivative of Siegel Eisenstein series and intersection numbers of special cycles on moduli spaces of Drinfeld shtukas with $r$ legs. The case $r=0$ (resp., $r=1$) can be viewed as an analogue of the Siegel–Weil formula (resp., the arithmetic Siegel–Weil formula). Over function fields, the possibility of relating higher derivatives of analytic objects to intersection numbers was first discovered by Yun–Zhang YZ17YZ19 in the context of the higher Gross–Zagier formula. Over number fields, however, no analogue of such a higher Siegel–Weil formula (resp., higher Gross–Zagier formula) is currently known when $r>1$. Feng–Yun–Zhang FYZ21a also defined higher theta series over function fields (including all singular terms) and conjectured their modularity. Previously, an arithmetic Siegel–Weil formula over function fields was proved by Wei Wei19 for special cycles on moduli spaces of Drinfeld modules of rank 2 with complex multiplication (analogue of the special case $n=r=1$).

5. Local arithmetic Siegel–Weil formula

In order to prove the arithmetic Siegel–Weil formula (Conjecture 4.2.2), one first notices that it can be reduced to a local identity:

(1)

Geometric side (LHS): The arithmetic intersection numbers corresponding to nonsingular matrices $T$ can be computed as a sum indexed by primes of $\mathbb{Z}$ (or the ring of integers $O_F$ in general). The local term at a finite prime $p$ can be further reduced to an arithmetic intersection on a Rapoport–Zink space $\mathcal{N}_n$, which is a local analogue of Shimura varieties over $\mathbb{Z}_p$ (or a completion of $O_F$ in general), via the theory of $p$-adic uniformization of Shimura varieties (RZ96).

(2)

Analytic side (RHS): The nonsingular Fourier coefficients $E_T(\tau ,s)$ has a product expansion indexed by primes of $\mathbb{Z}$, and thus the derivative $E_T'(\tau ,0)$ can also be written as a sum indexed by the primes of $\mathbb{Z}$. The term indexed by a finite prime $p$ can be further reduced to the derivative of the local representation density of quadratic/hermitian forms over $\mathbb{Z}_p$.

This reduction step is illustrated in the following diagram.

$$\begin{equation*} \vcenter{\img[][327pt][137pt][{$\xymatrix@R=4em{ {\begin{array}[h]{c} \text{arithmetic intersection number}\\ \text{of $n$ special divisors on } \mathcal{X} \end{array}} \ar@{&lt;~&gt;}[d] &amp; \questeq&amp; {\begin{array}[h]{c} \text{sum of Fourier coefficients}\\ \text{of a central derivative}\\ \text{of a Siegel Eisenstein series}\\ \text{on } \Sp(2n) \text{ or } \UU(n,n) \end{array}} \ar@{&lt;~&gt;}[d]\\ {\begin{array}[h]{c} \text{arithmetic intersection number}\\ \text{of $n$ special divisors on } \mathcal{N}_n \end{array}} &amp; \questeq&amp; {\begin{array}[h]{c} \text{central derivative}\\ \text{of the local density of a}\\ \text{hermitian or quadratic form}\\ \text{in $n$-variables} \end{array}} }$}]{Images/img47dddbf785087a2886a1e5d263f34d4d.svg}} \end{equation*}$$

The conjectural local identity on the bottom is known as the local arithmetic Siegel–Weil formula. The local arithmetic Siegel–Weil formula has been recently proved in our work with Zhang LZ22a (resp., LZ22b) in the unitary (resp., orthogonal) case. Next we will make this local conjecture more precise in the unitary case. In the unitary case, this reduction step was made precise by Kudla–Rapoport (KR11, KR14), and the local conjecture is also known as the local Kudla–Rapoport conjecture.

5.1. Geometric side

Let $p$ be an odd prime. Let $F_0$ be a finite extension of $\mathbb{Q}_p$ with residue field $k=\mathbb{F}_q$ and uniformer $\varpi$. Let $F$ be the unramified quadratic extension of $F_0$ (e.g., $F/F_0=\mathbb{Q}_{p^2}/\mathbb{Q}_p$). Associated to this datum, we have the unitary Rapoport–Zink space $\mathcal{N}_n$:

•

It is a formal scheme over $\operatorname {Spf}_{{O_{\breve{F}}}}$ of relative dimension $n-1$, parameterizing deformations (up to quasi-isogenies) of a fixed $O_F$-hermitian formal $p$-divisible group $\mathbb{X}_n/\bar{k}$ of relative height $2n$, dimension $n$, and signature $(1,n-1)$. Here ${O_{\breve{F}}}$ is the completion of the maximal unramified extension of $O_F$.

•

The space of special homomorphisms$$\begin{equation*} \mathbb{V}_n=\operatorname {Hom}_{O_F}(\mathbb{X}_1, \mathbb{X}_n)\otimes _{O_F} F \end{equation*}$$

has a structure of a (nonsplit) $F$-hermitian space of dimension $n$, coming from the principal polarization on $\mathbb{X}_n$.

•

The unitary group $\mathrm{U}(\mathbb{V}_n)$ naturally acts on $\mathcal{N}_n$ via the action on $\mathbb{X}_n$.

•

Each vector $x\in \mathbb{V}_n$ gives rise to a special divisor or Kudla–Rapoport (KR) divisor $\mathcal{Z}(x)\subseteq \mathcal{N}_n$, defined to be the locus where the homomorphism $x$ deforms. This is the local analogue of the special divisor considered in Conjecture 4.2.2.

The Rapoport–Zink space $\mathcal{N}_n$ is formally smooth over $\operatorname {Spf}{O_{\breve{F}}}$, but its geometric structure is rather complicated. For example, $\mathcal{N}_n$ is highly nonreduced: the reduced subscheme $\mathcal{N}_n^\mathrm{red}$ has dimension $\lfloor \frac{n-1}{2}\rfloor$, near the middle dimension of $\mathcal{N}_n$. The structure of $\mathcal{N}^\mathrm{red}_n$ was studied by Vollaard–Wedhorn (VW11), and they showed that $\mathcal{N}_n^\mathrm{red}$ has a nice stratification into smooth varieties, known as the Bruhat–Tits stratification. Each closed stratum of the Bruhat–Tits stratification is isomorphic to a smooth projective variety $\mathrm{DL}_i$ of dimension $i$ over $\bar{k}$ ($0\le i\le \lfloor \frac{n-1}{2}\rfloor$), and the incidence relation between the closed strata resembles the combinatorial structure of the Bruhat–Tits building for quasi-split unitary groups over $F_0$. Here each $\mathrm{DL}_i$ is a generalized Deligne–Lusztig variety associated to the unitary group $\mathrm{U}(2i+1)$ over $k$.

Example 5.1.1.

Take $n=1$. Then $\mathcal{N}_1\cong \operatorname {Spf}{O_{\breve{F}}}$, and $\mathcal{N}_1^\mathrm{red}=\{\text{pt}\}$.

Example 5.1.2.

Take $n=3$. Then $\mathcal{N}_3$ has relative dimension 2 over $\operatorname {Spf}{O_{\breve{F}}}$, while the reduced subscheme $\mathcal{N}_3^\mathrm{red}$ has dimension 1 (see Figure 4). In this case only two types of Deligne–Lusztig variety show up:

(1)

$\mathrm{DL}_0=\{\text{pt}\}$, a single point.

(2)

$\mathrm{DL}_1=\{x^{q+1}+y^{q+1}+z^{q+1}=0\}\subseteq \mathbb{P}^2$, the Fermat curve of degree $q+1$.

And $\mathcal{N}_3^\mathrm{red}$ is an infinite tree, where

(1)

the number of $\mathrm{DL}_1$ containing a given $\mathrm{DL}_0$ is exactly $q+1$.

(2)

the number of $\mathrm{DL}_0$ contained in a given $\mathrm{DL}_1$ is exactly $q^3+1$.

Definition 5.1.3.

Let $L\subseteq \mathbb{V}_n$ be an $O_F$-lattice of rank $n$. Let $x_1,\ldots , x_n$ be an $O_F$-basis of $L$. Define the special cycle or Kudla–Rapoport (KR) cycle

$$\begin{equation*} \mathcal{Z}(L)\coloneq \mathcal{Z}(x_1)\cap \cdots \cap \mathcal{Z}(x_n)\subseteq \mathcal{N}_n. \end{equation*}$$

Define the arithmetic intersection number

$$\begin{equation*} \operatorname {\mathrm{Int}}(L) \coloneq \chi (\mathcal{O}_{\mathcal{Z}(x_1)} \otimes ^\mathbb{L}\cdots \otimes ^\mathbb{L}\mathcal{O}_{\mathcal{Z}(x_n)} ), \end{equation*}$$

where $\chi$ denotes the Euler–Poincaré characteristic, $\mathcal{O}_{\mathcal{Z}(x_i)}$ denotes the structure sheaf of the Kudla–Rapoport divisor $\mathcal{Z}(x_i)$, and $\otimes ^\mathbb{L}$ denotes the derived tensor product of coherent sheaves on $\mathcal{N}_n$. It is known (by Terstiege Ter13, as extended in LZ22a, Corollary 2.8.2 or by Howard How19) that $\operatorname {\mathrm{Int}}(L)$ is independent of the choice of the basis $x_1,\ldots , x_n$ and hence is a well-defined invariant of $L$ itself, justifying the notation.

Remark 5.1.4.

When the intersection $\mathcal{Z}(L)$ is zero dimensional, we have

$$\begin{equation*} \mathcal{O}_{\mathcal{Z}(x_1)}\otimes ^\mathbb{L}\cdots \otimes ^\mathbb{L}\mathcal{O}_{\mathcal{Z}(x_n)}=\mathcal{O}_{\mathcal{Z}(x_1)}\otimes \cdots \otimes \mathcal{O}_{\mathcal{Z}(x_n)}=\mathcal{O}_{\mathcal{Z}(L)}, \end{equation*}$$

and thus $\operatorname {\mathrm{Int}}(L)$ is nothing but the ${O_{\breve{F}}}$-length of $\mathcal{O}_{\mathcal{Z}(L)}$ (which equals the sum of intersection multiplicities at all points). Even though $\mathcal{Z}(L)$ is the intersection of $n$ divisors in a $n$-dimensional formal scheme, in general $\mathcal{Z}(T)$ may be fail to have the expected dimension 0 due to improper intersection. In this case, the derived intersection is needed so that the intersection number is well behaved.

Example 5.1.5.

Take $n=3$ and $L\cong \langle \varpi \rangle ^3$ (the hermitian form with respect to an $O_F$-basis of $L$ is $\operatorname {diag}\{\varpi ,\varpi ,\varpi \}$). In this case the intersection $\mathcal{Z}(L)$ is not zero dimensional, and in fact

$$\begin{equation*} \mathcal{Z}(L)\cong \mathrm{DL}_1=\{x^{q+1}+y^{q+1}+z^{q+1}=0\}\subseteq \mathbb{P}^2. \end{equation*}$$

The arithmetic intersection number turns out to be its topological Euler characteristic

$$\begin{equation*} \operatorname {\mathrm{Int}}(L)=\chi (\mathrm{DL}_1)=2-2g=2-q(q-1)=2+q-q^2. \end{equation*}$$

5.2. Analytic side

Definition 5.2.1.

Let $L, M$ be two hermitian $O_F$-lattices of rank $n$, $m$, respectively. Let $S=\mathrm{Rep}(M, L)$ be the $O_{F_0}$-scheme such that for any $O_{F_0}$-algebra $R$,

$$\begin{equation*} S(R)=\mathrm{Herm}(L \otimes _{F_0}R, M \otimes _{F_0}R), \end{equation*}$$

where $\mathrm{Herm}$ denotes the set of hermitian module homomorphisms. The local density of representations of $M$ by $L$ is defined to be

$$\begin{equation*} \mathrm{Den}(M,L)\coloneq \lim _{N\rightarrow +\infty }\frac{\# S(O_{F_0}/\varpi ^N)}{q^{N\cdot \dim S_{F_0}}}. \end{equation*}$$

It gives a quantitative measure of “how many different ways” one can embed $L$ into $M$ as a hermitian submodule.

Example 5.2.2.

Consider $L=M=\langle 1\rangle ^n$, the rank $n$ lattice with hermitian form given by the identity matrix. Then one can compute that

$$\begin{equation*} \mathrm{Den}(\langle 1\rangle ^n, \langle 1\rangle ^n)=\prod _{i=1}^{n}(1-(-q)^{-i}). \end{equation*}$$

One can recognize it as the number $\frac{G(k)}{q^{\dim G}}$, where $G$ is the unitary group $\mathrm{U}(\langle 1\rangle ^n)$ defined over $O_{F_0}$, or in fancier language, as the local $L$-factor of the Gross motive Gro97 of the quasi-split unitary group in $n$ variables.

The local density $\mathrm{Den}(M, L)$ has nice compatibility when replacing $M$ by $M \oplus \langle 1\rangle ^k$ for $k\ge 0$. More precisely, it is known (Hir98, Theorem II) that $\mathrm{Den}(\langle 1\rangle ^{n+k}, L)$ is a polynomial in $(-q)^{-k}$ with $\mathbb{Q}$-coefficients.

Example 5.2.3 (KR11, p.677).

$$\begin{equation} \mathrm{Den}(\langle 1\rangle ^{n+k}, \langle 1\rangle ^n)=\prod _{i=1}^n(1-(-q)^{-i}X)\bigg |_{X= (-q)^{-k}}. {\tag{5.2.3.1}} \end{equation}$$

Definition 5.2.4.

Define the normalized Siegel series $\mathrm{Den}(X,L)\in \mathbb{Z}[X]$ such that

$$\begin{equation*} \mathrm{Den}((-q)^{-k},L)=\frac{\mathrm{Den}(\langle 1\rangle ^{n+k}, L)}{\mathrm{Den}(\langle 1\rangle ^{n+k}, \langle 1\rangle ^n)}. \end{equation*}$$

These polynomials $\mathrm{Den}(X,L)$ play an important role in computing the Fourier coefficients of Siegel Eisenstein series, and many works are devoted to proving more explicit formulas for them (see, e.g., Kit83Kat99Hir98Hir12IK22CY20). The local Siegel series satisfies a functional equation (Hir12, Theorem 5.3)

$$\begin{equation} \mathrm{Den}(X,L)=(-X)^{{\mathrm{val}}(L)}\cdot \mathrm{Den}\left(\frac{1}{X},L\right). {\tag{5.2.4.1}} \end{equation}$$

Here ${\mathrm{val}}(L)\coloneq {\mathrm{val}}(\det (L))\in \mathbb{Z}$ is the valuation of $L$. This may be viewed as a local analogue of the functional equation of Eisenstein series.

Definition 5.2.5.

If ${\mathrm{val}}(L)$ is odd (equivalently, $L \otimes _{O_F}F$ is a nonsplit hermitian space), then $\mathrm{Den}(1,L)=0$ by the functional equation 5.2.4.1. In this case, define the central derivative of the local density by

$$\begin{equation*} \partial \mathrm{Den}(L)\coloneq -\frac{\mathrm{d}}{\mathrm{d}X}\bigg |_{X=1}\mathrm{Den}(X,L). \end{equation*}$$

5.3. Local arithmetic Siegel–Weil formula

Now we are ready to state the local arithmetic Siegel–Weil formula in the unitary case, originally conjectured by Kudla–Rapoport KR11, Conjecture 1.3.

Theorem 5.3.1 (Local arithmetic Siegel–Weil formula, Li and Zhang LZ22a, Theorem 1.2.1).

Let $L\subseteq \mathbb{V}_n$ be an $O_F$-lattice of full rank $n$. Then

$$\begin{equation*} \operatorname {\mathrm{Int}}(L)=\partial \mathrm{Den}(L). \end{equation*}$$

Remark 5.3.2.

(i)

The theorem was proved by Kudla–Rapoport KR11 for $n=2$ and Terstiege Ter13 for $n=3$.

(ii)

The local arithmetic Siegel–Weil formula is proved for the orthogonal case in our work with Zhang LZ22b for arbitrary $n$. The case $n=3$ was previously proved by Gross–Keating GK93 and Terstiege Ter11.

(iii)

The local arithmetic Siegel–Weil formula is proved when the quadratic extension $F/F_0$ is ramified for exotic smooth models in our work with Liu LL21 for arbitrary even $n$, and for Krämer models⁠² by Shi and by He–Shi–Yang Shi20HSY20 for $n=2$.

²

During the refereeing process of this article, He–Shi–Yang HSY21 formulated a conjectural local arithmetic Siegel–Weil formula for Krämer models for arbitrary $n$ and proved it for the case $n=3$. The case for arbitrary $n$ has been proved in our work with He–Shi–Yang HLSY22.

✖

(iv)

It is more difficult to prove or formulate the local arithmetic Siegel–Weil formula in the presence of more general level structures (even when the quadratic extension $F/F_0$ is unramified). In the unitary case LZ22a formulates and proves a local arithmetic Siegel–Weil formula when the level is the parahoric subgroup given by the stabilizer of an almost self-dual lattice (the case $n=2$ was previously proved by Sankaran San17). Recently Cho Cho20 proposed a general formulation for all minuscule parahoric levels in the unitary case.

Example 5.3.3.

Take $n=3$ and $L\cong \langle \varpi \rangle ^3$. Specializing the formula of Cho–Yamauchi CY20 (extended to the unitary case in LZ22a, Theorem 3.5.1) gives

$$\begin{align*} \mathrm{Den}(X,L)&=(1-X)(1+qX)(1-q^2X)+(q^3+1)(1-X)X^2\\ & = 1-(1-q+q^2)X+(1-q+q^2)X^2-X^3. \end{align*}$$

It satisfies the functional equation

$$\begin{equation*} \mathrm{Den}(X,L)=-X^3\cdot \mathrm{Den}\left(\frac{1}{X},L\right). \end{equation*}$$

It is easy to compute

$$\begin{align*} \partial \mathrm{Den}(L)&=-\frac{\mathrm{d}}{\mathrm{d}X}\bigg |_{X=1}\mathrm{Den}(X,L)\\ &= (1-q+q^2)-2(1-q+q^2)+3\\ &= 2+q-q^2. \end{align*}$$

So combining Example 5.1.5 we obtain $\operatorname {\mathrm{Int}}(L)=\partial \mathrm{Den}(L)$ in this case! It is miraculous that the purely analytic quantity $\partial \mathrm{Den}(L)$ secretly knows about the Euler characteristic $\chi (\mathrm{DL}_1)=\operatorname {\mathrm{Int}}(L)$ of the Deligne–Lusztig curve $\mathrm{DL}_1$.

5.4. Strategy of the proof: uncertainty principle

The previously known special cases of Theorem 5.3.1 were proved via explicit computation of both the geometric and analytic sides. Explicit computation seems infeasible for the general case. The proof in LZ22a instead proceeds via induction on $n$ using the uncertainty principle, a standard tool from local harmonic analysis. Even for $n=2,3$, this proof is different from the previous proofs.

More precisely, for a fixed $O_F$-lattice $L^\flat \subseteq \mathbb{V}\coloneq \mathbb{V}_n$ of rank $n-1$, consider functions on $x\in \mathbb{V}\setminus L^\flat _F$,

$$\begin{equation*} \mathrm{Int}_{L^\flat }(x)\coloneq \mathrm{Int}(L^\flat +\langle x\rangle ),\quad \partial \mathrm{Den}_{L^\flat }(x)\coloneq \partial \mathrm{Den}(L^\flat +\langle x\rangle ). \end{equation*}$$

Then it remains to show the equality of the two functions

$$\begin{equation*} \mathrm{Int}_{L^\flat }=\partial \mathrm{Den}_{L^\flat }. \end{equation*}$$

By the definition of $\operatorname {\mathrm{Int}}(L)$ and $\partial \mathrm{Den}(L)$, both functions are easily seen to vanish when $x$ is nonintegral, i.e., $\mathrm{val}(x)<0$. Here $\mathrm{val}(x)$ denotes the valuation of the norm of $x$. By utilizing the inductive structure of Rapoport–Zink spaces and local densities, it is not hard to see that if $x\perp L^\flat$ with $\mathrm{val}(x)=0$, then

$$\begin{equation*} \mathrm{Int}_{L^\flat }(x)=\mathrm{Int}(L^\flat ),\quad \partial \mathrm{Den}_{L^\flat }(x)=\partial \mathrm{Den}(L^\flat ) \end{equation*}$$

for the lattice $L^\flat \subseteq \mathbb{V}_{n-1}\cong \langle x\rangle ^\perp _F$ of full rank $n-1$. Thus by induction on $n$, the difference function $\phi =\mathrm{Int}_{L^\flat }-\partial \mathrm{Den}_{L^\flat }$ already vanishes on a large subset

$$\begin{equation*} \{x \in \mathbb{V}: x\perp L^\flat , \mathrm{val}(x)\le 0\}. \end{equation*}$$

We would like to deduce that $\phi$ indeed vanishes identically. To this end, we apply the following uncertainty principle.

Proposition 5.4.1 (Uncertainty Principle, LZ22a, Proposition 8.1.6).

Let $\phi \in \mathscr{S} (\mathbb{V})$ be a Schwartz function on $\mathbb{V}$. If both $\phi$ and its Fourier transform $\hat{\phi }$ vanish on $\{x\in \mathbb{V}: \mathrm{val}(x)\le 0\}$, then $\phi =0$.

In other words, $\phi ,\hat{\phi }$ cannot simultaneously have small support unless $\phi =0$. Applying the uncertainty principle to the difference function $\phi$, we can then finish the proof as long as we get a good control over the support of $\hat{\phi }$. However, both functions $\mathrm{Int}_{L^\flat },\partial \mathrm{Den}_{L^\flat }$ have singularities along the hyperplane $L^\flat _F\coloneq L^\flat \otimes _{O_F}F\subseteq \mathbb{V}$. Intuitively, if $x$ is closer to $L^\flat$, then $\mathcal{Z}(x)$ and $\mathcal{Z}(L^\flat )$ intersect more improperly, which results in the blowup of $\operatorname {\mathrm{Int}}_{L^\flat }(x)$ along $L^\flat _F$. These singularities cause trouble in computing the Fourier transforms or even in showing that $\phi \in \mathscr{S} (\mathbb{V})$.

5.5. Strategy of the proof: decomposition and local modularity

To overcome this difficulty, we isolate the singularities by decomposing

$$\begin{equation*} \mathrm{Int}_{L^\flat }=\mathrm{Int}_{L^\flat ,\mathscr{H}}+\mathrm{Int}_{L^\flat ,\mathscr{V}},\quad \partial \mathrm{Den}_{L^\flat }=\partial \mathrm{Den}_{L^\flat ,\mathscr{H}}+\partial \mathrm{Den}_{L^\flat ,\mathscr{V}} \end{equation*}$$

into horizontal and vertical parts. Here on the geometric side $\mathrm{Int}_{L^\flat ,\mathscr{H}}$ (resp., $\mathrm{Int}_{L^\flat ,\mathscr{V}}$) is the contribution from the horizontal (resp., vertical part) of the KR cycles, illustrated in Figure 5 (by the red (resp., blue) part in the online color version of this article). One can hope to understand the horizontal part explicitly using deformation theory, and the vertical part using algebraic geometry over the residue field $\bar{k}$.

By the uncertainty principle and the induction on $n$, it remains to prove the following.

Theorem 5.5.1 (Key Theorem).

(KR1)

$\operatorname {\mathrm{Int}}_{L^\flat ,\mathscr{H}}=\partial \mathrm{Den}_{L^\flat ,\mathscr{H}}$.

(KR2)

$\operatorname {\mathrm{Int}}_{L^\flat ,\mathscr{V}}\in \mathscr{S} (\mathbb{V})$ and $\widehat{\operatorname {\mathrm{Int}}_{L^\flat ,\mathscr{V}}}=- \operatorname {\mathrm{Int}}_{L^\flat ,\mathscr{V}}$.

(KR3)

$\partial \mathrm{Den}_{L^\flat ,\mathscr{V}}\in \mathscr{S} (\mathbb{V})$ and  $\widehat{\partial \mathrm{Den}_{L^\flat ,\mathscr{V}}}=- \partial \mathrm{Den}_{L^\flat ,\mathscr{V}}$.

In other words, the horizontal part $\partial \mathrm{Den}_{L^\flat ,\mathscr{H}}$ matches $\mathrm{Int}_{L^\flat ,\mathscr{H}}$, and subtracting the horizontal parts removes the singularities along $L^\flat _F$ so that vertical parts are indeed in $\mathscr{S} (\mathbb{V})$. The extra invariance under Fourier transform,

$$\begin{equation} \widehat{\mathrm{Int}}_{L^\flat ,\mathscr{V}}=- \mathrm{Int}_{L^\flat ,\mathscr{V}},\quad \widehat{\partial \mathrm{Den}}_{L^\flat ,\mathscr{V}}=-\partial \mathrm{Den}_{L^\flat ,\mathscr{V}}, {\tag{5.5.1.1}} \end{equation}$$

can be thought of as a local modularity, by analogy with the global modularity of arithmetic generating series (such as in Bruinier–Howard–Kudla–Rapoport–Yang BHK$^{+}$20 discussed in Remark 3.5.5(v)) encoding an extra global ${\mathrm{SL}}_2$-symmetry.

5.6. Ingredients of the Key Theorem

Some key ingredients of the proof of the Key Theorem 5.5.1 include the following.

(KR1)

We describe explicitly the horizontal part of KR cycles in terms of Gross’s quasi-canonical liftings Gro86, using the work of Tate, Grothendieck–Messing and Breuil on the deformation theory of $p$-divisible groups.

(KR2)

On the geometric side we show 5.5.1.1 by reducing to the case of intersection with Deligne–Lusztig curves $\mathrm{DL}_1$. This reduction requires the Bruhat–Tits stratification of $\mathcal{N}^\mathrm{red}$ into the Deligne–Lusztig varieties $\mathrm{DL}_i$, as discussed in §5.1, and the Tate conjecture Tat94 for these Deligne–Lusztig varieties. We prove the latter by reducing to a cohomological computation of Lusztig Lus76.

(KR3)

On the analytic side we show 5.5.1.1 using Cho–Yamauchi’s explicit formula CY20 for $\partial \mathrm{Den}(L)$ in terms of weighted lattice counting, and reduce to a (rather subtle) lattice theoretic problem. (In fact we only show directly something weaker than 5.5.1.1, which is enough to imply Theorem 5.3.1, and we then deduce 5.5.1.1 a posteriori.)

6. Arithmetic inner product formula

In this last section we discuss an application of the arithmetic Siegel–Weil formula to the Beilinson–Bloch conjecture and the arithmetic inner product formula.

6.1. Birch and Swinnerton-Dyer conjecture

One long-standing problem in number theory is the determination of the rational points $E(\mathbb{Q})$ for an elliptic curve $E: y^2=x^3 +a x+b$ defined over $\mathbb{Q}$. The celebrated Birch and Swinnerton-Dyer (BSD) conjecture predicts a deep link between $E(\mathbb{Q})$ and its $L$-function $L(E,s)$. Define the algebraic rank

$$\begin{equation*} r_\mathrm{alg}(E)\coloneq \operatorname {rank}E(\mathbb{Q}) \end{equation*}$$

to be the rank of the finitely generated abelian group $E(\mathbb{Q})$. Define the analytic rank

$$\begin{equation*} r_\mathrm{an}(E)\coloneq \operatorname {ord}_{s=1} L(E,s) \end{equation*}$$

to be the order of vanishing of $L(E,s)$ at the central point $s=1$. The BSD conjecture predicts the rank equality between these two notions of ranks of seemingly different nature,

$$\begin{equation} r_\mathrm{an}(E)\questeq r_\mathrm{alg}(E), {\tag{6.1.0.1}} \end{equation}$$

and a refined BSD formula

$$\begin{equation} \frac{L^{(r)}(E,1)}{r!}\stackrel{?}{=}\frac{\Omega (E) R(E) \prod _p c_p(E)\cdot |\mathrm{Ш}(E)|}{|E(\mathbb{Q})_\mathrm{tor}|^2} {\tag{6.1.0.2}} \end{equation}$$

for the Taylor expansion of $L(E/\mathbb{Q},s)$ at $s=1$ (here $r=r_\mathrm{an}(E)$) in terms of various important arithmetic invariants of $E$. Among these invariants are the order of the (mysterious) Tate–Shafarevich group $\mathrm{Ш}(E)$, and the regulator $R(E)\coloneq \det (\langle P_i, P_j\rangle _{\text{NT}})_{r\times r}$, where

$$\begin{equation*} \langle \ , \ \rangle _\text{NT}: E(\mathbb{Q})\times E(\mathbb{Q})\rightarrow \mathbb{R} \end{equation*}$$

is the Néron–Tate height pairing and $\{P_i\}$ is a basis of the free part of $E(\mathbb{Q})$.

The BSD conjecture is still widely open in general, but much progress has been made in the rank 0 or 1 case. The seminal work of Gross–Zagier GZ86, Kolyvagin Kol90 proved the implications

$$\begin{equation} r_\mathrm{an}(E)=0\Rightarrow r_\mathrm{alg}(E)=0,\quad r_\mathrm{an}(E)=1\Rightarrow r_\mathrm{alg}(E)=1, {\tag{6.1.0.3}} \end{equation}$$

confirming 6.1.0.1 when $r_\mathrm{an}(E)\le 1$. Due to the work of many people, many cases of 6.1.0.2 are also known when $r_\mathrm{an}(E)\le 1$.

The key to relate $r_\mathrm{an}(E)$ and $r_\mathrm{alg}(E)$ is the Gross–Zagier formula

$$\begin{equation} L'(E,1)\doteq \langle P, P\rangle _\mathrm{NT}, {\tag{6.1.0.4}} \end{equation}$$

(up to an explicit nonzero constant, including the period of $E$ and other rational factors) relating the first derivative of the $L$-function at $s=1$ and the Néron–Tate height of certain rational points $P$ on $E$ known as Heegner points (cf. Example 6.4.1). It gives one crucial implication in 6.1.0.3,

$$\begin{equation} r_\mathrm{an}(E)=1\Longrightarrow r_\mathrm{alg}(E)\ge 1. {\tag{6.1.0.5}} \end{equation}$$

The tools of Heegner points and $L$-functions, linked via the Gross–Zagier formula, are indispensable in studying the arithmetic of elliptic curves. We refer to Zhang Zha14a for an excellent recent survey on Heegner points and the Birch–Swinnerton-Dyer conjecture (see also Gross Gro04 and Darmon Dar04).

6.2. Beilinson–Bloch conjecture

In the 1980s, Beilinson (Bem̆87, Conjecture 5.9) and Bloch (Blo84a, Recurring Fantasy and Blo84b, Conjecture) proposed vast generalizations of the BSD conjecture to higher dimensional varieties.

Let $X$ be a smooth projective variety of dimension $n$ over a number field $F$. For $1\le m\le n$, denote by ${\mathrm{CH}}^m(X)$ the Chow group of algebraic cycles of codimension $m$ on $X$ defined over $F$ (up to rational equivalence), and ${\mathrm{CH}}^m(X)^0\subseteq {\mathrm{CH}}^m(X)$ the subgroup of geometrically cohomologically trivial cycles. Denote by $L(H^{2m-1}(X), s)$ the motivic $L$-function associated to the $(2m-1)$-th étale cohomology $H^{2m-1}(X_{\overline{F}}, \mathbb{Q}_\ell )$. Then the Beilinson–Bloch (BB) conjecture predicts the equality between analytic and algebraic ranks

$$\begin{equation} \operatorname {ord}_{s=m}L(H^{2m-1}(X), s)= \operatorname {rank}{\mathrm{CH}}^m(X)^0, {\tag{6.2.0.1}} \end{equation}$$

and a refined formula

$$\begin{equation*} L^{(r)}(H^{2m-1}(X),m)\doteq \det \langle \ , \ \rangle _{\mathrm{BB}} \end{equation*}$$

for the leading coefficient at $s=m$ in terms of the determinant of the Beilinson–Bloch height pairing⁠³

³

When $m>1$, the Beilinson–Bloch height pairing is only defined assuming certain conjectures on algebraic cycles on $X$ (see Bem̆87, Conjectures 2.2.1 and 2.2.3). This important technical issue is addressed in LL21LL22 so that the RHS of 6.7.1.1 in Theorem 6.7.1 can be defined unconditionally, but we will intentionally ignore it for the purpose of this article.

✖

$$\begin{equation*} \langle \ ,\ \rangle _\mathrm{BB}: {\mathrm{CH}}^m(X)^0\times {\mathrm{CH}}^{n+1-m}(X)^0\rightarrow \mathbb{R}. \end{equation*}$$

Example 6.2.1.

When $X/F=E/\mathbb{Q}$ and $m=1$, we recover the BSD conjecture as

$$\begin{equation*} {\mathrm{CH}}^1(E)^0\simeq E(\mathbb{Q}), L(H^1(E), s)=L(E,s) \text{ and } \langle \ ,\ \rangle _\mathrm{BB}=-\langle \ ,\ \rangle _\mathrm{NT}. \end{equation*}$$

The BB conjecture is even more elusive than the BSD conjecture: in general we do not know that ${\mathrm{CH}}^m(X)^0$ is finitely generated, nor do we know that $L(H^{2m-1}(X), s)$ can be analytically continued to the central point $s=m$, so neither side of 6.2.0.1 is well defined! This may be more an exciting challenge than disappointment for mathematicians—after all we were in a similar circumstance when the BSD conjecture was formulated in the 1960s: we knew neither the analytic continuation of $L(E,s)$ (except when $E$ has complex multiplication) nor the finiteness of $\mathrm{Ш}(E)$ in order to make sense of either side of 6.1.0.2.

A good testing ground for the BB conjecture is by taking $X$ to be Shimura varieties. By Langland’s philosophy, their $L$-functions can be computed in terms of automorphic $L$-functions, so the analytic rank in 6.2.0.1 becomes accessible. Even though we do not known if ${\mathrm{CH}}^m(X)^0$ is finitely generated, it still makes sense to test if it is nontrivial. It is thus tempting to relate special cycles on Shimura varieties to automorphic $L$-functions, and in particular, to generalize the Gross–Zagier formula 6.1.0.4 to higher dimensional Shimura varieties and prove the analogue of 6.1.0.5 towards the BB conjecture

$$\begin{equation} \operatorname {ord}_{s=m}L(H^{2m-1}(X), s)=1\Longrightarrow \operatorname {rank}{\mathrm{CH}}^m(X)^0\ge 1. {\tag{6.2.1.1}} \end{equation}$$

Here we use the notation $\operatorname {rank}{\mathrm{CH}}^m(X)^0\ge 1$ to stand for the nontriviality of ${\mathrm{CH}}^m(X)^0_\mathbb{Q}$.

6.3. Arithmetic inner product formula

As discussed in §3.5, Shimura varieties $X$ of type $\mathrm{GSpin}(n-1,2)$ and $\mathrm{U}(n-1,1)$ admit a rich supply of special cycles, recovering Heegner points in the case of modular curves. As explained in Kud04, III, when $n$ is even, the arithmetic Siegel–Weil formula together with the doubling method of Piatetski-Shapiro–Rallis PSR84 has important applications to the arithmetic inner product formula of the form

$$\begin{equation*} L'(1/2, \pi )\doteq \langle \Theta _\varphi (\phi ), \Theta _\varphi (\phi )\rangle _\mathrm{BB}. \end{equation*}$$

Here

•

$\pi$ is a certain cuspidal automorphic representation on ${\mathrm{Mp}}(n)$ or $\mathrm{U}(n)$.

•

$L'(1/2,\pi )$ is the central derivative the standard $L$-function $L(s,\pi )$ of $\pi$ (cf. Yam14) with global root number $\varepsilon (\pi )=-1$.

•

$\Theta _\varphi (\phi )$ is an algebraic cycle on $X$ of codimension $n/2$, constructed from holomorphic forms $\phi \in \pi$ using the method of arithmetic theta lifting (see Definition 6.6.3).

This arithmetic inner product formula can be viewed as an arithmetic analogue of the Rallis inner product formula (see GQT14) of the form

$$\begin{equation} L(1/2, \pi )\doteq \langle \theta _\varphi (\phi ),\theta _\varphi (\phi )\rangle _\mathrm{Pet}, {\tag{6.3.0.1}} \end{equation}$$

which relates the central value $L(1/2,\pi )$ when $\varepsilon (\pi )=+1$ and the Petersson inner product of the (usual) theta lift $\theta _\varphi (\phi )$ discussed in Remark 2.2.5. It can also be viewed as a higher dimensional generalization of the Gross–Zagier formula (the case $n=2$). Thus the arithmetic Siegel–Weil formula is intimately linked to the Birch–Swinnerton-Dyer conjecture, and more generally the Beilinson–Bloch conjecture for higher dimensional Shimura varieties.

The conjectural arithmetic inner product formula was formulated by Kudla Kud97b using the Gillet–Soulé height and in more generality by Liu Liu11a using the Beilinson–Bloch height (in the unitary case). In the unitary case the arithmetic inner product formula has been recently proved under mild local assumptions and Kudla’s modularity conjecture in our works with Liu LL21LL22, which has been applied to prove the first unconditional theorem for the Beilinson–Bloch conjecture for higher dimensional Shimura varieties. Our remaining goal is to explain some details about the main results of LL21LL22.

6.4. Arithmetic theta lifting on modular curves

To motivate the construction of arithmetic theta lifting, let us first consider an example of Heegner points on elliptic curves.

Example 6.4.1 (cf. Zag85).

Consider the curve $37a1$ in Cremona’s table,

$$\begin{equation*} E=37a1: y^2 + y = x^3 - x, \end{equation*}$$

It is the rank one optimal curve over $\mathbb{Q}$ of smallest conductor ($N=37$). It is isomorphic to the modular curve $X_0^+(N)=X_0(N)/\langle w_{N}\rangle$ for $N=37$ ($N$ is the smallest such that $X_0^+(N)$ has positive genus). The Mordell-Weil group $E(\mathbb{Q})\cong \mathbb{Z}$, generated by $P=(0,0)$. Let

$$\begin{equation*} \begin{split} f_E&=\sum _{n\ge 1}a_nq^n\\ &=q - 2 q^2- 3 q^3+ 2 q^4- 2 q^5+ 6 q^6- q^7+ 6 q^9+ 4 q^{10}- 5 q^{11}- 6 q^{12} +\cdots \in S_2(N) \end{split} \end{equation*}$$

be the weight 2 newform of level $N=37$ associated to $E$. It gives a modular parametrization

$$\begin{gather*} \varphi _E: \mathcal{H}\rightarrow E=\mathbb{C}/\Lambda ,\\ \varphi _E(\tau )=\sum _{n\ge 1}\frac{a_nq^n}{n}=q-q^2-q^3+\frac{1}{2}q^4-\frac{2}{5}q^5+q^6-\frac{1}{7}q^7+\cdots \end{gather*}$$

such that $\varphi _E^*(\omega _E)=2\pi if_E(\tau )d\tau$. The Shimura–Waldspurger correspondence, which can be viewed as an instance of the theta correspondence for the pair $(G,H)=(\mathrm{Mp}(2),\mathrm{O}(3))$,

$$\begin{equation*} \theta : S_{3/2}^+(4N) \rightarrow S_2(N), \end{equation*}$$

gives a weight 3/2 newform in Kohnen’s plus space $S_{3/2}^+(4N)$,

$$\begin{equation*} g_E=\sum _{d\ge 1}c_dq^d=-q^3-q^4+q^7-q^{11}+q^{12}+2q^{16}+3q^{27}+\cdots -6q^{67}+\cdots \end{equation*}$$

such that $\theta (g_E)=f_E$.

Let $d$ be a positive integer such that $-d\equiv 0,1\ (\mathrm{mod}\ 4)$. If $N=37$ splits in $K=\mathbb{Q}(\sqrt {-d})$, then one can construct a rational point $P_d\in E(\mathbb{Q})$ using Heegner points on $X_0(N)$ associated to quadratic orders of discriminant $-d$. For example, when $-d$ is a fundamental discriminant, we have a Heegner point $y_K\in X_0(H_K)$ defined over the Hilbert class field $H_K$ of $K$, and $P_d=\operatorname {tr}_{H_K/K}\varphi _E(y_K-\infty )$. The point $P_d$ may depend on the choice of $d$, even when $E(\mathbb{Q})\cong \mathbb{Z}$. In Table 3, we compute a list of Heegner points $P_d\in E(\mathbb{Q})$ for small $d$’s, and also the integer multiples $n_d$ such that $P_d=n_d\cdot P$ for the generator $P=(0,0)$.

From Table 3, we observe the miraculous coincidence that the integer $n_d$ exactly matches the coefficient $c_d$ of $g_E$! In other words, the generating series

$$\begin{equation*} Z(\tau )=\sum _{d\ge 1}P_d q^d=g_E\cdot P\in S_{3/2}^+(4\cdot 37) \otimes E(\mathbb{Q}) \end{equation*}$$

is a modular form valued in $E(\mathbb{Q}) \otimes \mathbb{C}$. This allows us to define the arithmetic theta lift by taking the Petersson inner product of $Z(\tau )$ with $g_E$,

$$\begin{equation*} \Theta (g_E)\coloneq \langle Z(\tau ), g_E\rangle _{\mathrm{Pet}}=\langle g_E,g_E\rangle _\mathrm{Pet }\cdot P \in E(\mathbb{Q})\otimes \mathbb{C}, \end{equation*}$$

which is now a canonical element in $E(\mathbb{Q}) \otimes \mathbb{C}$ (i.e., no need to choose any particular $d$). The arithmetic inner product formula in this case asserts the identity

$$\begin{equation} L'(E,1)\doteq \frac{\langle \Theta (g_E), \Theta (g_E)\rangle _\mathrm{NT}}{\langle g_E, g_E\rangle _\mathrm{Pet}}. {\tag{6.4.1.1}} \end{equation}$$

In fact, in this case we have $\frac{\langle \Theta (g_E), \Theta (g_E)\rangle _\mathrm{NT}}{\langle g_E, g_E\rangle _\mathrm{Pet}}=\langle g_E, g_E\rangle _\mathrm{Pet}\cdot \langle P, P\rangle _\mathrm{NT}$ and we can compute explicitly

$$\begin{align*} \langle g_E, g_E\rangle _\mathrm{Pet}&=0.7146356107\cdots =\frac{3\omega ^+(E)}{4\pi },\\ \langle P, P\rangle _\mathrm{NT}&= 0.0511114082\cdots ,\\ L'(E,1)&= 0.3059997738 \cdots . \end{align*}$$

Here $\omega ^+(E)$ is the real period of $E$. Therefore

$$\begin{equation*} \frac{L'(E,1)}{\langle g_E, g_E\rangle _\mathrm{Pet}\cdot \langle P, P\rangle _\mathrm{NT}}=8.37758040956\dots =\frac{8\pi }{3} \end{equation*}$$

and the equality 6.4.1.1 indeed holds up to an elementary constant $\frac{8\pi }{3}$.

6.5. Unitary Shimura varieties

Now let us come to the setting of unitary Shimura varieties. Let $E/F$ be a CM extension of a totally real number field. Let $\mathbb{V}$ be a totally definite incoherent hermitian space over $\mathbb{A}_E$ of rank $n$. Here “incoherent” means that $\mathbb{V}$ is not the base change of a global $E/F$-hermitian space, or equivalently the product of the local Hasse invariants of the local hermitian spaces $\mathbb{V}_v\coloneq \mathbb{V} \otimes _{\mathbb{A}_F}F_v$ is

$$\begin{equation*} \prod _{v} \varepsilon (\mathbb{V}_v)=-1. \end{equation*}$$

Picking any place $w|\infty$ of $F$ gives a nearby global $E/F$-hermitian space $V$ such that

$$\begin{equation*} V_v\cong \mathbb{V}_v,\text{ if } v\ne w, \quad \text{but $V_{w}$ has signature $(n-1,1)$.} \end{equation*}$$

Associated to any open compact subgroup $K\subseteq \mathrm{U}(V)(\mathbb{A}_F^\infty )$, we have a unitary Shimura variety $X=\mathrm{Sh}_{\mathrm{U}(V)}$ (cf. Z19Gro20), which has a smooth canonical model of dimension $n-1$ over $E$ (viewed as a subfield of $\mathbb{C}$ via the embedding induced by a place above $w$) and admits complex uniformization

$$\begin{equation*} X(\mathbb{C})=\mathrm{U}(V)(F)\backslash [ \mathbb{D}\times \mathrm{U}(V)(\mathbb{A}_F^\infty )/K], \end{equation*}$$

where