Arithmetic statistics of modular symbols

TL;DR

This paper proves conjectures about the statistical behavior of modular symbols, including their moments and distribution, using analytic techniques involving Eisenstein series, advancing understanding of their asymptotic properties.

Contribution

It establishes average growth of moments and the Gaussian distribution of modular symbols, refining previous conjectures with new analytic proofs.

Findings

Proved asymptotic growth of first and second moments of modular symbols.

Established Gaussian distribution of normalized modular symbols with cusp restrictions.

Used analytic properties of Eisenstein series twisted by modular symbols.

Abstract

Mazur, Rubin, and Stein have recently formulated a series of conjectures about statistical properties of modular symbols in order to understand central values of twists of elliptic curve $L$-functions. Two of these conjectures relate to the asymptotic growth of the first and second moments of the modular symbols. We prove these on average by using analytic properties of Eisenstein series twisted by modular symbols. Another of their conjectures predicts the Gaussian distribution of normalized modular symbols ordered according to the size of the denominator of the cusps. We prove this conjecture in a refined version that also allows restrictions on the location of the cusps.

Tables

Table 1. Table 1. The variance shifts for 15.a1

$d$	$D_{f,d}$	Experimental Shift
1	-0.440048	-0.440
3	-0.244592	-0.246
5	-0.153710	-0.153
15	-0.041745	-0.040

Equations

$$\{\infty ,\mathfrak{a}\}\mapsto 2\pi i{\int }_{i\infty }^{\mathfrak{a}}f(z)dz.$$

$$⟨r⟩=2\pi i{\int }_{i\infty }^r\alpha .$$

$$⟨\gamma ,\alpha ⟩:=⟨\gamma (\infty )⟩=2\pi i{\int }_{z_0}^{\gamma z_0}\alpha ,$$

$$f(z)=\sum _{n\ge 1}a(n)e(nz).$$

$$⟨r{⟩}^{\pm }=\pi i{\int }_{i\infty }^{r}f(z)dz\pm \pi i{\int }_{i\infty }^{-r}f(z)dz.$$

$$\tau (\chi )L(E,\bar{\chi },1)=\sum _{a\in (\mathbb{Z}/m\mathbb{Z}{)}^{\ast }}\chi (a)⟨a/m{⟩}^{\pm }$$

$$G_c(x)=\frac{1}{c}\sum _{0\le a\le cx}⟨a/c⟩.$$

$$g(x)$$

$$G_c(x)\to g(x).$$

$$c^{-1}\sum _{0\le a\le cx}{\int }_{a/c+i\delta }^{i\infty }\alpha \to {\int }_{[0,x]\times [\delta ,\infty ]}\alpha ,$$

$$E(f,c)=\frac{1}{\varphi (c)}\sum _{\begin{array}{c}a\mathrm{mod}c\\ (a,c)=1\end{array}}⟨a/c⟩,\mathrm{Var}(f,c)=\frac{1}{\varphi (c)}\sum _{\begin{array}{c}a\mathrm{mod}c\\ (a,c)=1\end{array}}{\left(⟨a/c⟩-E(f,c)\right)}^2,$$

$$\underset{\stackrel{c\to \infty }{(c,q)=d}}{\lim }(\mathrm{Var}(f,c)-C_f\log c)=D_{f,d}.$$

$$C_f=-\frac{6}{{\pi }^2}\prod _{p|q}(1+p^{-1}{)}^{-1}L({\mathrm{sym}}^{2}f,1).$$

$$\frac{⟨a/c⟩}{(C_f\log c+D_{f,d}{)}^{1/2}},c\in \mathbb{N}\text{ with }(c,q)=d,a\in (\mathbb{Z}/c\mathbb{Z}{)}^{\ast }$$

$$\frac{1}{M}\sum _{1\le c\le M}{G}_c(x)\to g(x),\text{ as }M\to \infty .$$

$$\frac{1}{\sum _{\stackrel{c\le M}{(c,q)=d}}\varphi (c)}\sum _{\begin{array}{c}c\le M\\ (c,q)=d\end{array}}\varphi (c)(\mathrm{Var}(f,c)-C_f\log c)\to D_{f,d},\text{ as }M\to \infty ,$$

$$D_{f,d}=A_{d,q}L({\mathrm{sym}}^{2}f,1)+B_q{L}^{\prime }({\mathrm{sym}}^{2}f,1).$$

$$\begin{array}[]{ccc}Q_{d}\cap I&\to&\displaystyle{\mathbb{R}}\\ \displaystyle\frac{a}{c}&\mapsto&\frac{\displaystyle\langle{a}/{c}\rangle}{\displaystyle(C_{f}\log c)^{1/2}},\end{array}$$

$$T_{\mathfrak{a}\mathfrak{b}}=\left\{r=\frac{a}{c}\mathrm{mod}1,\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)\in {\Gamma }_{\infty }\backslash {\sigma }_{\mathfrak{a}}^{-1}\Gamma {\sigma }_{\mathfrak{b}}/{\Gamma }_{\infty }\text{ and }c>0\right\}\subseteq \mathbb{R}/\mathbb{Z},$$

$$r=\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)\infty \mathrm{mod}1.$$

$$T_{\mathfrak{a}\mathfrak{b}}(M)=\{r\in T_{\mathfrak{a}\mathfrak{b}},c(r)\le M\}.$$

$$⟨r{⟩}_{\mathfrak{a}\mathfrak{b}}=2\pi i{\int }_{\mathfrak{b}}^{{\sigma }_{\mathfrak{a}}r}\alpha .$$

$$\sum _{r\in T_{\mathfrak{a}\mathfrak{b}}(M)}⟨r{⟩}_{\mathfrak{a}\mathfrak{b}}{1}_{[0,x]}(r)=$$

$$+O_f(M^{2-\delta }).$$

$$\frac{\sum _{r\in T_{\mathfrak{a}\mathfrak{b}}(M)}⟨r{⟩}_{\mathfrak{a}\mathfrak{b}}^2}{\sum _{r\in T_{\mathfrak{a}\mathfrak{b}}(M)}1}=C_f\log M+D_{f,\mathfrak{a}\mathfrak{b}}+o_f(1),asM\to \infty ,$$

$$C_f=\frac{-16{\pi }^2{\Vert f\Vert }^2}{\mathrm{vol}(\Gamma \backslash \mathbb{H})}.$$

$$g(r)=\frac{⟨r{⟩}_{\mathfrak{a}\mathfrak{b}}}{\sqrt{C_f\log c(r)}}$$

$$\frac{\#\left\{r\in T_{\mathfrak{a}\mathfrak{b}}(M)\cap I,\frac{⟨r{⟩}_{\mathfrak{a}\mathfrak{b}}}{\sqrt{C_f\log c(r)}}\in [a,b]\right\}}{\#(T_{\mathfrak{a}\mathfrak{b}}(M)\cap I)}\to \frac{1}{\sqrt{2\pi }}{\int }_a^b\exp \left(-\frac{t^2}{2}\right)dt,$$

$$⟨\gamma ,\alpha ⟩:=⟨\gamma (\infty )⟩=2\pi i{\int }_{z_0}^{\gamma z_0}\alpha$$

$${\chi }_{\epsilon }(\gamma )=\exp \left(2\pi i\epsilon {\int }_{z_0}^{\gamma z_0}\alpha \right).$$

$$A_{\mathfrak{a}}(z)=2\pi i{\int }_{\mathfrak{a}}^z\alpha .$$

Full text

The second author was supported by a Sapere Aude grant from The Danish Council for Independent Research (Grant-id:0602-02161B). The first author would like to acknowledge the hospitality of the University of Copenhagen, while this paper was written.

1. Introduction

Modular symbols are fundamental tools in number theory. By the work of Birch, Manin, Cremona and others they can be used to compute modular forms, the homology of modular curves, and to gain information about elliptic curves and special values of $L$-functions. In this paper we study the arithmetical properties of the modular symbol map

[TABLE]

Here $f\in S_{2}(\Gamma)$ is a holomorphic cusp form of weight $2$ for the group ${\Gamma}={\Gamma}_{0}(q)$, and $\{\infty,{\mathfrak{a}}\}$ denotes the homology class of curves between the cusps $\infty$ and ${\mathfrak{a}}$.

For our purposes it is convenient to work with the real-valued, cuspidal one-form ${\alpha}=\Re(f(z)dz)$. The finite cusps ${\mathfrak{a}}$ are parametrized by ${\mathbb{Q}}$, so for $r\in{\mathbb{Q}}$ we write

[TABLE]

Note that the path can be taken to be the vertical line connecting $r\in\mathbb{Q}$ to $\infty$.

If ${\mathfrak{a}}$ is equivalent to $\infty$ under the ${\Gamma}$-action such that $\{\infty,{\mathfrak{a}}\}=\{\infty,{\gamma}(\infty)\}$ for some ${\gamma}\in{\Gamma}$ we write the map (1.1) as

[TABLE]

where in the last expression we have replaced $\infty$ by any $z_{0}\in{\mathbb{H}}^{*}$.

Mazur, Rubin, and Stein [MR16, Ste15] have recently formulated a series of conjectures about the value distribution of $\langle r\rangle$. We now describe these conjectures. Let $E$ be an elliptic curve over ${\mathbb{Q}}$ of conductor $q$ with associated holomorphic weight 2 cusp form $f(z)$. We write the Fourier expansion of $f$ at $\infty$ as

[TABLE]

It is a fundamental question in number theory to understand how often the central value of $L(E,\chi,s)$, vanishes, when $\chi$ runs over the characters of $\hbox{Gal}(\bar{\mathbb{Q}}/{\mathbb{Q}})$.

Mazur, Rubin, and Stein used raw modular symbols. For $r\in{\mathbb{Q}}$ these are defined by

[TABLE]

This corresponds roughly to taking $\alpha$ in (1.2) to be the real or imaginary part of the 1-form $f(z)dz$. See Remark 2.1 for the precise statement. Modular symbols and the central value of twists of the corresponding $L$-function are related by the Birch–Stevens formula e.g. [Pol14, Eq. 2.2], [MTT86, Eq 8.6] that is

[TABLE]

for a primitive character of conductor $m$ (here the choice of $\pm$ corresponds to the sign of $\chi$). To understand the vanishing of $L(E,\bar{\chi},s)$ at $s=1$ Mazur, Rubin, and Stein were led to investigate the distribution of modular symbols and theta constants. In this paper we investigate modular symbols but not theta constants.

Mazur, Rubin, and Stein studied computationally the statistics of (raw) modular symbols. Since $f$ has period $1$, the same is true for the modular symbols: $\langle r+1\rangle=\langle r\rangle.$ They observed the behavior of contiguous sums of modular symbols, defined for each $x\in[0,1]$ by

[TABLE]

Based on their computations they defined

[TABLE]

and arrived at the following conjecture.

Conjecture 1.1* (Mazur–Rubin–Stein).*

As $c\to\infty$ we have

[TABLE]

They added credence to this conjecture with the following heuristics. If we cut off the paths in (1.2) for modular symbols at height $\delta>0$, then

[TABLE]

because the left-hand side is a Riemann sum for the integral. The heuristics involves interchanging this limit with the limit as $\delta\to 0$.

In another direction Mazur and Rubin investigated the distribution of $\langle a/c\rangle$ for $(a,c)=1$ as $c\to\infty$. Define the usual mean and variance by

[TABLE]

where $\phi$ is Euler’s totient function. They conjectured the following asymptotic behavior of the variance.

Conjecture 1.2* (Mazur–Rubin).*

There exist a constant $C_{f}$ and constants $D_{f,d}$ for each divisor $d$ of $q$, such that

[TABLE]

Moreover,

[TABLE]

The constant $C_{f}$ is called the variance slope and the constant $D_{f,d}$ the variance shift. In Conjecture 1.2 the symmetric square $L$-function $L(\textrm{sym}^{2}f,s)$ is normalized such that 1 is at the edge of the critical strip.

Moreover, the numerics suggest that the normalized raw modular symbols obey a Gaussian distribution law.

Conjecture 1.3* (Mazur–Rubin).*

Let $d|q$. The data

[TABLE]

has limit the standard normal distribution.

In this paper we prove average versions of Conjectures 1.1 and 1.2 when we average over $c$. We only work with $q$ squarefree. The restriction to $q$ squarefree may not be necessary. On the other hand we can only work with averages over $c$ and not individual $c$.

Moreover, we prove a refined version of Conjecture 1.3. We can restrict the rational number $a/c$ to lie in any prescribed interval and we can restrict to rational numbers $a/c$ with $(c,q)$ a fixed number.

To prove these results we specialize to $\Gamma_{0}(q)$ more general results on modular symbols for cofinite Fuchsian groups with cusps. Here is a statement of our results for $\Gamma_{0}(q)$, when $q$ is squarefree.

Theorem 1.4**.**

Conjecture 1.1 holds on average. More precisely we have

[TABLE]

Remark 1.5*.*

In fact we can restrict the summation to $(c,q)=d$ for $d$ a fixed divisor of $q$. See Corollary 8.1 and the discussion following it.

Theorem 1.6**.**

Conjecture 1.2 holds on average. More precisely let $d|q$, and let $C_{f}$ be given by (1.4). Then

[TABLE]

where

[TABLE]

Here $A_{d,q}$, $B_{q}$ are explicitly computable constants given by (8.12).

Theorem 1.7**.**

Let $I\subseteq{\mathbb{R}}/{\mathbb{Z}}$ be an interval of positive length, and consider for $d|q$ the set $Q_{d}=\{{a}/{c}\in{\mathbb{Q}},(a,c)=1,(c,q)=d\}$. Then the values of the map

[TABLE]

ordered according to $c$, have asymptotically a standard normal distribution.

Remark 1.8*.*

Putting $I={\mathbb{R}}/{\mathbb{Z}}$ in Theorem 1.7 we prove Conjecture 1.3. The difference in normalization, i.e. the appearance of $D_{f,d}$ in the denominator, is irrelevant as explained in Remark 7.9.

Theorems 1.4, 1.6, and 1.7 will follow rather easily from the following general theorems for finite covolume Fuchsian groups $\Gamma$ with cusps. Let ${\mathfrak{a}}$, ${\mathfrak{b}}$ be cusps of $\Gamma$, not necessarily distinct. Let $\alpha=\Re(f(z)dz)$, where $f\in S_{2}(\Gamma)$ is a holomorphic cusp form of weight $2$. We do not assume that $f$ has real coefficients. Fix scaling matrices $\sigma_{{\mathfrak{a}}}$ and $\sigma_{\mathfrak{b}}$ for the two cusps. Define

[TABLE]

see Proposition 2.2. Note that

[TABLE]

We denote the denominator of $r$ by $c(r)$. We order the elements of $T_{{\mathfrak{a}}{\mathfrak{b}}}$ according to the size of $c(r)$ and define

[TABLE]

We define general modular symbols as

[TABLE]

Theorem 1.9**.**

Let $x\in[0,1]$. Let $a_{{\mathfrak{a}}}(n)$ be the Fourier coefficients of $f$ at the cusp ${\mathfrak{a}}$. There exists a $\delta>0$ such that, as $M\to\infty$,

[TABLE]

Theorem 1.10**.**

Let $\left\lVert f\right\rVert$ be the Petersson norm of $f$. There exists an explicit constant $D_{f,{\mathfrak{a}}{\mathfrak{b}}}$, depending on ${\Gamma}$, $f$, ${\mathfrak{a}},{\mathfrak{b}}$, which we call the variance shift, such that

[TABLE]

where

[TABLE]

The constant $C_{f}$ is the variance slope.

The formula for $D_{f,{\mathfrak{a}}{\mathfrak{b}}}$ is explicit but complicated, see (7.4). It depends on $\left\lVert f\right\rVert$, the period $\int_{{\mathfrak{a}}}^{{\mathfrak{b}}}\alpha$ and data from the Kronecker limit formula for the Eisenstein series for the cusps ${\mathfrak{a}}$ and ${\mathfrak{b}}$.

Theorem 1.11**.**

Let $I\subseteq{\mathbb{R}}/{\mathbb{Z}}$ be an interval of positive length. The values of the map $g:T_{{\mathfrak{a}}{\mathfrak{b}}}\cap I\to{\mathbb{R}}$ with

[TABLE]

ordered according to $c(r)$ have asymptotically a standard normal distribution, that is, for every $a,b\in[-\infty,\infty]$ with $a\leq b$ we have

[TABLE]

as $M\to\infty$.

Remark 1.12*.*

We have obtained different but related normal distribution results for modular symbols in [PR04, Ris04, PR05, PR09]. One difference between these papers and the current one is in the ordering and normalization of the values of $\left\langle\gamma,\alpha\right\rangle$. The orderings in these papers were more geometric (action of $\Gamma$ on ${\mathbb{H}}$) and less arithmetic, and the modular symbols used closed paths in ${\Gamma\backslash{\mathbb{H}}}$. However, in Theorem 1.7 we need to combine statistics from various cusps, since not all rational cusps are equivalent to $\infty$ for $\Gamma_{0}(q)$. Moreover, we allow to restrict $\gamma(\infty)=a/c$ to a general $I\subseteq{\mathbb{R}}/{\mathbb{Z}}$. This is a new feature, making our current results significantly more refined.

The expression for the variance shift in Theorems 1.6 and 1.10, see also (7.4), is very explicit. This is an unexpected facet. The analogue in [PR04, Theorem 2.19] involves also the reduced resolvent of the Laplace operator, which is much harder to understand.

Remark 1.13*.*

To prove Theorem 1.11 we study the asymptotic $k$th moments of the modular symbols $\langle r\rangle_{{\mathfrak{a}}{\mathfrak{b}}}$ for all $k$, see Theorem 7.5. We make no effort to optimize the error bounds for the moments in Theorems 1.9, 1.10, and 7.5, but they can all be made explicit in terms of spectral gaps.

Remark 1.14*.*

An important tool in this paper and in [PR04] is non-holomorphic Eisenstein series twisted with modular symbols $E^{m,n}(z,s)$. These were introduced by Goldfeld [Gol99b, Gol99a] and studied extensively by many authors, see e.g. [O’S00, DO00, JO08, JOS16, BD16] and the references therein. For their definition see (2.11). In [PR04] we used the Eisenstein series twisted with the $k$th power of modular symbols as a generating series itself to study the $k$th moment of modular symbols. In this paper we need to understand the $n$th Fourier coefficient of $E_{{\mathfrak{a}}}^{(k)}(\sigma_{\mathfrak{b}}z,s)$, which involves twists by the $k$th power of $\langle r\rangle_{{\mathfrak{a}}{\mathfrak{b}}}$. The reason why the results here are more arithmetic is because the Fourier coefficients of Eisenstein series and of twisted Eisenstein series encode arithmetic data and modular symbols. To isolate the $n$th coefficient we use inner products with Poincaré series.

Remark 1.15*.*

The structure of the paper is as follows. In Section 2 we introduce the generating functions (Dirichlet series) $L_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s,m,n)$ for the powers of modular symbols. We also introduce Poincaré series and Eisenstein series twisted by powers of modular symbols.

In Section 3 we analyze $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ for $\Re(s)>1/2$, and conclude with the statement that $T_{{\mathfrak{a}}{\mathfrak{b}}}$ is equidistributed modulo $1$.

In Section 4 we study Eisenstein series twisted by modular symbols, or rather a related series $D_{{\mathfrak{a}}}^{(k)}(z,s)$. We prove the analytic continuation for $\Re(s)>1/2$ and study the order of the poles and leading singularity at $s=1$. The crucial identity is Eq. (4.4) that allows us to understand $D_{{\mathfrak{a}}}^{(k)}(z,s)$ inductively using the resolvent of the Laplace operator $R(s)$.

In Section 5 we find explicit expressions for the functional equations of $D^{(k)}_{\mathfrak{a}}(z,s)$ and $E^{(k)}_{{\mathfrak{a}}}(z,s)$, see Theorems 5.2 and 5.4.

In Section 6 we study the analytic properties of the derivatives $L_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s,0,n)$ for $k\geq 1$. When $k=1$ we find the residue at $s=1$, and the whole singular part when $k=2$. Finally, we identify the order of the pole and the leading singularity for all $k$.

In Section 7 we prove Theorems 1.9, 1.10, and 1.11 for general finite covolume Fuchsian groups with cusps. We use the method of contour integration.

In Section 8 we specialize the general results to $\Gamma_{0}(q)$ for squarefree $q$. This leads to the proofs of Theorems 1.4, 1.6, and 1.7.

Remark 1.16*.*

In a recent preprint [KS17] the authors prove that, when $c\to\infty$ through primes, the limiting behavior in Conjecture 1.1 holds. Their method is different from ours.

Acknowledgments

We would like to thank Peter Sarnak for alerting us to [MR16]. Also we would like to thank Barry Mazur and Karl Rubin for useful comments, and for providing us with numerical data.

2. Generating series for powers of modular symbols

From now on we allow ${\Gamma}$ to be any cofinite Fuchsian group with cusps. All implied constants in our estimates depend on $\Gamma$ and $f$. In this section we define a generating series for modular symbols, and explain how it can be understood in terms of derivatives of Eisenstein series with characters.

2.1. Modular symbols

The modular symbols defined by

[TABLE]

are independent of $z_{0}\in{\mathbb{H}}^{*}$, and independent of the path between $z_{0}$ and $\gamma z_{0}$. Fix a set of inequivalent cusps for $\Gamma$. For such a cusp ${\mathfrak{a}}$ we fix a scaling matrix $\sigma_{{\mathfrak{a}}}$, i.e. a real matrix of determinant 1 mapping $\infty$ to ${\mathfrak{a}}$, and satisfying ${\Gamma}_{\mathfrak{a}}=\sigma_{{\mathfrak{a}}}\Gamma_{\infty}\sigma_{{\mathfrak{a}}}^{-1}$. Here $\Gamma_{\mathfrak{a}}$ is the stabilizer of ${\mathfrak{a}}$ in ${\Gamma}$, and ${\Gamma}_{\infty}$ is the standard parabolic subgroup. We have $\left\langle\gamma,{\alpha}\right\rangle=0$ for $\gamma$ parabolic since ${\alpha}$ is cuspidal.

For any real number $\varepsilon$ and every cuspidal real differential 1-form ${\alpha}$ we have a family of unitary characters $\chi_{\epsilon}:{\Gamma}\to S^{1}$ defined by

[TABLE]

Note that this character is the conjugate of the one we considered in [PR04]. We also need the antiderivative of ${\alpha}$. We define

[TABLE]

We expand ${\alpha}$ at a cusp ${\mathfrak{b}}$. Let us assume that ${\alpha}=\Re(f(z)dz)$, where $f(z)\in S_{2}(\Gamma)$, and

[TABLE]

where $j(\gamma,z)=cz+d$. Then

[TABLE]

and

[TABLE]

We consider the line integral

[TABLE]

By [JO08, Eq. (3.3), (3.5)] we have the uniform bound

[TABLE]

for $z\in{\mathbb{H}}$. Consequently, we have the estimate

[TABLE]

Remark 2.1*.*

If $f(z)$ be a cusp form for ${\Gamma}_{0}(q)$ with real Fourier coefficients at infinity, then $\overline{f(z)}=f(-\bar{z})$. Since $\overline{\pi i\int_{i\infty}^{r}f(z)dz}=\pi i\int_{i\infty}^{-r}f(z)dz$ it follows that

[TABLE]

where ${\alpha}_{g}=\Re(g(z)dz)$. Consequently, taking $\alpha=\Re(f(z)dz)$ covers both cases of raw modular symbols.

2.2. The generating series

Let ${\mathfrak{a}},{\mathfrak{b}}$ be cusps and $\sigma_{\mathfrak{a}},\sigma_{\mathfrak{b}}$ fixed scaling matrices. Then we define

[TABLE]

It is easy to see that ${a}/{c}\bmod 1$ is well-defined for the double coset containing $\gamma$. Note also that $\gamma\infty=a/c\bmod 1$.

Proposition 2.2**.**

Let $r\in T_{{\mathfrak{a}}{\mathfrak{b}}}$. Then there exists a unique

[TABLE]

satisfying

[TABLE]

Proof.

We imitate [Iwa02, p. 50]. Assume

[TABLE]

and $r=\gamma\infty$, $r^{\prime}=\gamma^{\prime}\infty$. We may assume that $0\leq a<c$ and $0\leq a^{\prime}<c^{\prime}$. The matrix $\gamma^{\prime\prime}=\gamma^{\prime-1}\gamma\in\sigma_{{\mathfrak{b}}}^{-1}\Gamma\sigma_{{\mathfrak{b}}}$ has lower left entry $-ac^{\prime}+a^{\prime}c$. If this is zero then $\gamma^{\prime\prime}\in\Gamma_{\infty}$ and $r=r^{\prime}\bmod 1$. If not, then by [Iwa02, Eq. (2.30)] we have $\left\lvert-ac^{\prime}+a^{\prime}c\right\rvert\geq c_{{\mathfrak{b}}}$, which implies that

[TABLE]

Therefore $r\neq r^{\prime}$ and since $0\leq r,r^{\prime}<1$ this implies that $r\neq r^{\prime}\bmod 1$. ∎

From Proposition 2.2 we conclude the following result:

Corollary 2.3**.**

Any $r\in T_{{\mathfrak{a}}{\mathfrak{b}}}$ determines a unique number $c(r)>0$ and unique cosets $a(r)\bmod c(r)$, and $d(r)\bmod c(r)$ satisfying $a(r)d(r)=1\bmod c(r)$ and $a(r)/c(r)=r\bmod 1$.

From Corollary 2.3 it follows that there exists a unique pair $(c(r),a(r))$ such that $0\leq a(r)<c(r)$ and $a(r)/c(r)=r\bmod 1$. We can therefore put an ordering on the elements of $T_{{\mathfrak{a}}{\mathfrak{b}}}$ by putting the lexicographical ordering on the $(c(r),a(r))$ i.e.

[TABLE]

For $r\in T_{{\mathfrak{a}}{\mathfrak{b}}}$ we define $\bar{r}\bmod 1$ by $\bar{r}=d(r)/c(r)$. For $r\in T_{{\mathfrak{a}}{\mathfrak{b}}}$ we define

[TABLE]

We suppress ${\alpha}$ from the notation. We notice that if $\gamma\in\Gamma_{\infty}\backslash\sigma_{{\mathfrak{a}}}^{-1}\Gamma\sigma_{\mathfrak{b}}/\Gamma_{\infty}$ corresponds to $r$ as in Proposition 2.2 we have

[TABLE]

which we will also refer to as a modular symbol. The map $r\mapsto\langle r\rangle_{{\mathfrak{a}}{\mathfrak{b}}}$ does not grow too fast in terms of $c(r)$:

Proposition 2.4**.**

The following estimate holds

[TABLE]

Proof.

We use (2.3) with $\sigma_{{\mathfrak{a}}}\gamma\sigma_{{\mathfrak{b}}}^{-1}$ and a fixed $z$. Writing the lower row of ${\gamma}$ as $(c(r),d(r))$ we may assume that $\left\lvert d(r)\right\rvert\leq\left\lvert c(r)\right\rvert$. Writing $\sigma_{{\mathfrak{b}}}^{-1}z=w$, we use the elementary inequalities

[TABLE]

from which the result follows. ∎

We can now define the main generating series:

Definition 2.5**.**

For $\Re(s)>1$ we define

[TABLE]

By Proposition (2.4) and [Iwa02, Prop 2.8] we see that $L^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,m,n)$ is absolutely convergent for $\Re(s)>1$, and uniformly convergent on compacta of $\Re(s)>1$. It is the analytic properties of this series that will eventually allow us to prove our main results.

2.3. Relation to Eisenstein series

To explain how $L^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,m,n)$ relates to Eisenstein series we recall that the generalized Kloosterman sums are defined by

[TABLE]

For $m,n\in{\mathbb{Z}}$ we define

[TABLE]

where the sum is over $c$ for $\begin{pmatrix}*&*\\ c&*\end{pmatrix}\in\sigma_{\mathfrak{a}}^{-1}\Gamma\sigma_{\mathfrak{b}}.$ This is a version of the Selberg–Linnik zeta function. When ${\varepsilon}=0$, that is when the character is trivial, we omit it from the notation. Using [Iwa02, Prop. 2.8] we see that for $\Re(s)>1$ this is an absolutely converging Dirichlet series. Note that

[TABLE]

where $r$ and $\begin{pmatrix}a&b\\ c&d\end{pmatrix}$ are related as in Proposition 2.2. It therefore follows that

[TABLE]

Proposition 2.6**.**

For any cusps ${\mathfrak{a}},{\mathfrak{b}}$ and any $m,n\in{\mathbb{Z}}$ we have

[TABLE]

Proof.

This follows easily from (2.4), and the following basic properties of Kloosterman sums: By inspection we see that

[TABLE]

Also, we have

[TABLE]

as seen by inverting ${\gamma}$ in the definition of the Kloosterman sums. ∎

We now recover $L^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ as Fourier coefficients of Eisenstein series twisted by modular symbols. For $m\in{\mathbb{N}}\cup\{0\}$ we define Poincaré series

[TABLE]

When $m=0$, i.e. in the case of the usual Eisenstein series we will often omit the subscript $m$. When we omit to specify ${\varepsilon}$ in the notation, we have set ${\varepsilon}=0$. For $m>0$ it is known that $E_{{\mathfrak{a}},m}(z,s,{\varepsilon})\in L^{2}({\Gamma\backslash{\mathbb{H}}})$ and that it admits meromorphic continuation to $s\in{\mathbb{C}}$, see [GS83, p. 247]. The usual Eisenstein series has Fourier expansion at a cusp ${\mathfrak{b}}$ given by (see e.g [Sel89, p. 640–641])

[TABLE]

where

[TABLE]

As usual $\delta_{{\mathfrak{a}}{\mathfrak{b}}}=1$ if ${\mathfrak{a}}={\mathfrak{b}}$ and is [math] otherwise.

The derivatives of $\phi_{{\mathfrak{a}}{\mathfrak{b}}}(s,{\varepsilon})$ and $\phi_{{\mathfrak{a}}{\mathfrak{b}}}(n,s,{\varepsilon})$ in ${\varepsilon}$ are given by

[TABLE]

It follows that the series $L_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s,0,n)$ can be understood by understanding derivatives of Eisenstein series.

We consider the $k$th derivative in ${\varepsilon}$ of $E_{\mathfrak{a}}(z,s,{\varepsilon})$ at ${\varepsilon}=0$, which we denote by $E_{\mathfrak{a}}^{(k)}(z,s)$. It is easily seen that, when $\Re(s)>1$, we have

[TABLE]

This series is absolutely and uniformly convergent on compact subsets of $\Re(s)>1$, as can be seen by (2.3) and using the standard region of absolute convergence of the Eisenstein series $E_{{\mathfrak{a}}}(z,s)$. We note that for $k\in{\mathbb{N}}$

[TABLE]

Here the Fourier expansion is computed by termwise differentiation of the Fourier expansion of $E_{{\mathfrak{a}}}(z,s,{\varepsilon})$ at the cusp ${\mathfrak{b}}$, which is allowed. Hence, the generating series $L_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s,0,n)$ for the modular symbols appear as Fourier coefficients of $E_{{\mathfrak{a}}}^{(k)}(z,s)$, see (2.8) and (2.9).

2.4. Automorphic Poincaré series with modular symbols

While the generating series for the modular symbols appear as Fourier coefficients of $E_{{\mathfrak{a}}}^{(k)}(z,s)$, the series $E_{{\mathfrak{a}}}^{(k)}(z,s)$ are not automorphic when $k>0$. They are higher order modular forms. Properties of such forms has been studied extensively in many papers such as [CDO02, DS06, DD09, DS08, Sre09, IO09, Dei11, BD12].

Since $2\left\langle\gamma,{\alpha}\right\rangle=\left\langle\gamma,f(z)dz\right\rangle+\overline{\left\langle\gamma,f(z)dz\right\rangle}$, we see that our $E^{(k)}(z,s)$ is indeed a linear combination of the Eisenstein series twisted by powers of $\left\langle\gamma,f(z)dz\right\rangle$ and $\overline{\left\langle\gamma,f(z)dz\right\rangle}$ with $m+n=k$:

[TABLE]

For background on $E_{{\mathfrak{a}}}^{m,n}(z,s)$ see Remark 1.14. For our purpose it is convenient to consider a related function, which has the advantage of also being automorphic: Recall (2.1). Let $m,k\in{\mathbb{N}}\cup\{0\}$. For $\Re(s)>1$ we define

[TABLE]

Using (2.2), and by comparison with the standard Eisenstein series, we see that the function $D^{(k)}_{{\mathfrak{a}},m}(z,s)$ is absolutely and uniformly convergent on compact subsets of $\Re(s)>1$. It follows immediately that $D^{(k)}_{{\mathfrak{a}},m}(z,s)$ is ${\Gamma}$-automorphic, and holomorphic for $s$ in this half-plane.

We consider also

[TABLE]

so that $D^{(k)}_{{\mathfrak{a}},m}(z,s)$ is the $k$th derivative in ${\varepsilon}$ of $D_{{\mathfrak{a}},m}(z,s,{\varepsilon})$ at ${\varepsilon}=0$. As usual when $m=0$ we omit it from the notation.

We now explain how $E_{{\mathfrak{a}},m}^{(k)}(z,s)$ and $D_{{\mathfrak{a}},m}^{(k)}(z,s)$ are related. We arrange our Eisenstein and Poincaré series in column vectors indexed by the cusps as follows:

[TABLE]

and

[TABLE]

We define the diagonal matrix $U(z,{\varepsilon})$ with diagonal entries

[TABLE]

so that

[TABLE]

Let $A(z)$ be the diagonal matrix with diagonal entries the antiderivatives $2\pi i\int_{{\mathfrak{a}}}^{z}\alpha=A_{{\mathfrak{a}}}(z)$. It follows from (2.13) by differentiation at ${\varepsilon}=0$ that for $\Re(s)>1$ we have the vector equations

[TABLE]

Hence, we can freely translate between $D_{m}^{(k)}(z,s)$ and $E_{m}^{(k)}(z,s)$.

3. The generating series $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$

In this section we discuss the analytic properties of $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ for $\Re(s)>1/2$. In order to do so we first need some general bounds on Eisenstein series and Poincaré series with modular symbols.

3.1. Bounds on Eisenstein series

We first discuss a construction of Eisenstein series of weight $k$ for $\sigma=\Re(s)>1/2$, generalizing the construction for weight [math] in [CdV83, Section 2]. This approach is useful to estimate Eisenstein series in the cusps for $\Re(s)>1/2$.

For $\gamma\in{\hbox{SL}_{2}({\mathbb{R}})}$ we define $j_{\gamma}(z)=(cz+d)/\left\lvert cz+d\right\rvert$. Let

[TABLE]

be the Laplace operator of weight $k$ and $\tilde{\Delta}_{k}$ the closure of $\Delta_{k}$ acting on smooth, weight $k$ functions such that $f,\Delta_{k}f$ are square integrable. We fix a fundamental domain $F$ of $\Gamma$, and notice that $\sigma_{{\mathfrak{a}}}^{-1}F$ is a fundamental domain for $\sigma_{{\mathfrak{a}}}^{-1}\Gamma\sigma_{{\mathfrak{a}}}$. However, the spectral analysis of the $k$-Laplacian remains the same, since the manifold ${\Gamma\backslash{\mathbb{H}}}$ is isometric to $\sigma_{{\mathfrak{a}}}^{-1}\Gamma\sigma_{{\mathfrak{a}}}\backslash{\mathbb{H}}$. Recall the decomposition of $F$ as $F(Y)\cup_{{\mathfrak{a}}}F_{{\mathfrak{a}}}(Y)$ for $Y$ sufficiently large, see [Iwa02, p. 40].

Lemma 3.1**.**

Let $k\in\mathbb{Z}$. Let $h(y)$ be a smooth function that is identically $1$ for $y\geq Y+1$ and identically [math] for $y<Y+1/2$. Then for $\Re(s)>1/2$ and $s(1-s)\not\in{\operatorname{spec}}(-\tilde{\Delta}_{k})$ there exists a unique function $E_{{\mathfrak{a}}}(z,s,k)$ satisfying the eigenvalue equation

[TABLE]

and such that

[TABLE]

Moreover, the $L^{2}$-norm in (3.1) is $O_{k}((2\sigma-1)^{-1})$.

Proof.

If such a solution exists we denote the left-hand side of (3.1) by $g(z,s,k)$. We apply $\Delta_{k}+s(1-s)$ to deduce

[TABLE]

where

[TABLE]

is a compactly supported function, is independent of $k$, and is holomorphic in $s$.

We now define $H(z,s)$ by (3.2), and apply to it the resolvent $R(s,k)=(\tilde{\Delta}_{k}+s(1-s))^{-1}$ of $\tilde{\Delta}_{k}$, which is defined for $\Re(s)>1/2$ with $s(1-s)\not\in{\operatorname{spec}}(-\tilde{\Delta}_{k})$. This produces a unique function $g(z,s,k)\in L^{2}(\sigma_{{\mathfrak{a}}}^{-1}\Gamma\sigma_{{\mathfrak{a}}}\backslash{\mathbb{H}},k)$. The standard inequality for the operator norm of the resolvent

[TABLE]

allows to estimate $\left\lVert g(z,s,k)\right\rVert_{L^{2}}\ll(2\sigma-1)^{-1}$. ∎

Lemma 3.2**.**

Let $\Re(s)\in[1+{\varepsilon},A]$. Then for cusps ${\mathfrak{a}},{\mathfrak{b}}$ we have

[TABLE]

In particular

[TABLE]

Proof.

We use

[TABLE]

and the estimate [Iwa02, corollary 3.5], which can be made uniform for $\Re(s)\in[1+{\varepsilon},A]$, $y\geq Y_{0}$. We estimate the right-hand side of (3.4) on the compact part $F(Y)$ of $F$ and on the cuspidal zones

[TABLE]

see [Iwa02, p. 40]. For different cusps ${\mathfrak{a}}$ and ${\mathfrak{b}}$, the matrix $\sigma_{\mathfrak{a}}^{-1}\sigma_{{\mathfrak{b}}}$ has nonzero $c$ bounded from below by $c({{\mathfrak{a}},{\mathfrak{b}}})$, see [Iwa02, Section 2.6]. We have $\Im(\sigma_{{\mathfrak{a}}}^{-1}\sigma_{{\mathfrak{b}}}z)^{\sigma}\leq y^{\sigma}/(cy)^{2\sigma}\ll y^{-\sigma}$. For ${\mathfrak{a}}={\mathfrak{b}}$ the terms $y^{\sigma}$ cancel. ∎

As usual we define the Eisenstein series of weight $k$ by

[TABLE]

Using Lemma 3.2 we see that

[TABLE]

This equation and Lemma 3.1 show that $E_{{\mathfrak{a}}}(z,s,k)$ agrees with the construction of Lemma 3.1. Therefore, the conclusions of Lemma 3.1 hold for the Eisenstein series of weight $k$ in the region $\Re(s)>1/2$, $s(1-s)\not\in{\operatorname{spec}}(-\tilde{\Delta}_{k})$.

We can also estimate $D^{(k)}_{{\mathfrak{a}},m}(z,s)$ defined in (2.12) using (2.2) and Lemma 3.2. We write

[TABLE]

If $m>0$, and $1+{\varepsilon}\leq\Re(s)\leq A$ we see that the contribution of the identity term $A_{\mathfrak{a}}(z)^{k}e(m\sigma_{\mathfrak{a}}^{-1}z)\Im(\sigma_{\mathfrak{a}}^{-1}z)^{s}$ decays exponentially at the cusp ${\mathfrak{a}}$ i.e. $\ll\exp(-(2\pi-{\varepsilon})y)$ and is $O(y^{-\sigma+{\varepsilon}})$ at the other cusps. When $m=0$, $k>0$ we notice that in the cuspidal zone for ${\mathfrak{a}}$, the expansion (2.1) shows that $A_{{\mathfrak{a}}}(\sigma_{{\mathfrak{a}}}z)$ decays exponentially. We can deduce that for all cusps ${\mathfrak{b}}$ we have $D^{(k)}_{\mathfrak{a}}(\sigma_{{\mathfrak{b}}}z,s)=O(y^{1-\sigma+{\varepsilon}})$. These estimates show that $D^{(k)}_{{\mathfrak{a}},m}(z,s)$ is square integrable uniformly in the strip $1+{\varepsilon}\leq\Re(s)\leq A$ for $m>0$ or $k>0$.

Lemma 3.3**.**

Let $m\geq 1$, $1/2+{\varepsilon}<\Re(s)<1+{\varepsilon}$, $1+2{\varepsilon}<\Re(w)<A$. Moreover, assume that $s(1-s)\not\in{\operatorname{spec}}(-\tilde{\Delta})$. Then the $\Gamma$-invariant functions

[TABLE]

belong to $L^{1}({\Gamma\backslash{\mathbb{H}}})$. In fact

[TABLE]

where the implied constant depends only on $k,\Gamma,{\varepsilon}$, and $A$.

Proof.

We use the result from Lemma 3.1 and the notation in its proof. We have

[TABLE]

We need to analyze also

[TABLE]

It suffices to concentrate on the cuspidal sector for ${\mathfrak{a}}$, since $h(\Im(\sigma_{{\mathfrak{a}}}^{-1}z))$ vanished elsewhere. We have

[TABLE]

We use (3.5). When ${\mathfrak{b}}={\mathfrak{a}}$, we see that

[TABLE]

decays exponentially, otherwise it is bounded by $h(y)y^{\sigma-\Re(w)+{\varepsilon}}$. We easily see that $\int_{Y}^{\infty}h(y)y^{\sigma}y^{1-\Re(w)+{\varepsilon}}d\mu(z)$ is bounded. ∎

3.2. Meromorphic continuation and bounds on $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$

It is well known that the inner product of an automorphic function $G$ and $E_{\mathfrak{b}}(z,s,{\varepsilon})$ is directly related to the $n$th Fourier coefficient of $G$ at the cusp ${\mathfrak{b}}$. We first use this to find an integral expression for $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ for $n>0$.

Lemma 3.4**.**

Let $n>0$ and $\Re(s),\Re(w)>1$. The Dirichlet series $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n,{\varepsilon})$ has the integral representation

[TABLE]

Proof.

We unfold and, using (2.5), we find

[TABLE]

where we have used [GR07, 6.621.3, p. 700]. The result follows from the above computations and (2.7). ∎

We can now use the integral expression in Lemma 3.4 to find the analytic properties of $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$.

Lemma 3.5**.**

For any cusps ${\mathfrak{a}},{\mathfrak{b}}$, and $n\in{\mathbb{Z}}$ the Dirichlet series

[TABLE]

admits meromorphic continuation to $s\in{\mathbb{C}}$. For $n\neq 0$ the continuation is holomorphic at $s=1$, while for $n=0$ the continuation has a simple pole with residue

[TABLE]

For $1/2+{\varepsilon}<\Re(s)<1+{\varepsilon}$, and $s(1-s)$ bounded away from ${\operatorname{spec}}(-\tilde{\Delta})$, the following estimate holds:

[TABLE]

Proof.

It is well known that the Fourier coefficients of $E_{\mathfrak{a}}(z,s)$ admit meromorphic continuation to $s\in{\mathbb{C}}$ and the meromorphic continuation of the functions $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ follows from (2.6) and (2.7).

To analyze further when $\Re(s)>1/2$, we use Lemma 3.4 when $n\neq 0$. Using Proposition 2.6 we see that it suffices to consider $n>0$. We let $w=1+2{\varepsilon}+it$, where $t=\Im(s)$. With this choice of $w$ Stirling’s formula gives that the quotient of Gamma factors is $O(\left\lvert s\right\rvert^{1/2})$. The claim about growth on vertical lines follows from Lemma 3.3. For $n>0$ the residue at $s=1$ has $\left\langle 1,E_{{{\mathfrak{b}}},n}(z,\bar{w},0)\right\rangle$ as a factor. But this vanishes by unfolding.

For $n=0$ we examine (2.6). Using the well-known fact that $\phi_{{\mathfrak{a}}{\mathfrak{b}}}(s)$ has a pole of order $1$ at $s=1$ with residue $1/{\operatorname{vol}}({\Gamma\backslash{\mathbb{H}}})$, the claim for the residue of $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,0)$ follows. For the growth on vertical lines, we observe that $\left\lvert\phi_{{\mathfrak{a}}{\mathfrak{b}}}(s)\right\rvert=O(1)$ in the region $\Re(s)\geq 1/2+{\varepsilon}$, and away from the spectrum of $\Delta$, see [Sel89, Eq. (8.5)–(8.6), p. 655]. The result follows from Stirling’s formula. ∎

Remark 3.6*.*

The analytic properties of $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,m,n,{\varepsilon})$ for $mn\neq 0$ have been studied in [GS83].

Remark 3.7*.*

If ${\Gamma}$ is the full modular group and ${\mathfrak{a}}={\mathfrak{b}}=\infty$ with trivial scaling matrices then $T_{{\mathfrak{a}}{\mathfrak{b}}}={\mathbb{Q}}/{\mathbb{Z}}$ so that in this case

[TABLE]

Here $S(0,n,c)$ is the standard Ramanujan sum. We notice that in this case the bound in Lemma 3.5 on $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ for $n\in{\mathbb{Z}}$ is far from optimal.

3.3. Equidistribution of $T_{{\mathfrak{a}}{\mathfrak{b}}}$.

We can use Lemma 3.5 to observe that $T_{{\mathfrak{a}}{\mathfrak{b}}}$ is equidistributed modulo 1. The generating series of $e(nr)$ for $r\in T_{{\mathfrak{a}}{\mathfrak{b}}}$ is $L_{{\mathfrak{a}}{\mathfrak{b}}}(s,n,0)$. Proposition 2.6 allows us to understand its behavior through $L_{{\mathfrak{b}}{\mathfrak{a}}}(s,0,-n)$. We define

[TABLE]

Using contour integration and the polynomial estimates on vertical lines for $L_{{\mathfrak{b}}{\mathfrak{a}}}(s,0,-n)$ in Lemma 3.5 we deduce that, for some $\delta>0$ depending on the spectral gap for $\Delta$,

[TABLE]

These are the Weyl sums for the sequence $T_{{\mathfrak{a}}{\mathfrak{b}}}$. As usual $\delta_{m}(n)=1$ if $n=m$ and is [math] otherwise and similarly for a set $A$, $\delta_{A}(n)=1$ if $n\in A$ and is [math] otherwise. Good studied the asymptotics of such Weyl sums in [Goo83].

To state our results we need to introduce norms on $L^{1}$-functions on ${\mathbb{R}}/{\mathbb{Z}}$: For $A\geq 0$ let

[TABLE]

where $\hat{h}(n)$ denotes the $n$th Fourier coefficient of $h$. We then have the following result:

Theorem 3.8**.**

The sequence $T_{{\mathfrak{a}}{\mathfrak{b}}}$ is equidistributed modulo $1$, i.e. for any continuous function $h:{\mathbb{R}}/{\mathbb{Z}}\to{\mathbb{C}}$ we have

[TABLE]

If, moreover, $\left\lVert h\right\rVert_{H^{1/2}}<\infty,$ then

[TABLE]

Proof.

The first claim follows from Weyl’s equidistribution criterion. For the second we write the Fourier series of $h$, which is absolutely convergent, and interchange the summation over $n$ and over $r\in T_{{\mathfrak{a}}{\mathfrak{b}}}(M)$. We have

[TABLE]

∎

Remark 3.9*.*

In Theorem 3.8 the Sobolev norm $\left\lVert h\right\rVert_{H^{1/2}}$ can be replaced by $\left\lVert h\right\rVert_{H^{{\varepsilon}}}$ for any $0<{\varepsilon}<1/2$, at the expense of a worse exponent $\delta=\delta({\varepsilon})$. This can be done by making the appropriate changes in the contour integration argument leading to (3.6). Note that if we do so the exponent depends both on the spectral gap and on ${\varepsilon}$.

Remark 3.10*.*

Usually equidistribution modulo 1 of a sequence $\{x_{n}\}_{n=1}^{\infty}$ is stated as

[TABLE]

In Theorem 3.8 we are only looking at a subsequence of the left-hand side, namely only $N$ equal to $\#T_{{\mathfrak{a}}{\mathfrak{b}}}(M)$ for some $M$. It is possible to consider the whole sequence by noticing that it follows from (3.6) that

[TABLE]

This improves on the trivial bound in [Iwa02, (2.37)].

4. Eisenstein series with modular symbols

In this section we analyze further the series $D^{(k)}_{{\mathfrak{a}},m}(z,s)$ defined in (2.12), when $m=0$. We will omit $m$ from the notation and simply write $D_{{\mathfrak{a}}}^{(k)}(z,s)$. We will show that this function admits meromorphic continuation and prove $L^{2}$-bounds for it. Let

[TABLE]

Define

[TABLE]

Remark 4.1*.*

These definitions are motivated by perturbation theory. Even if we are not using perturbation theory directly, as in our previous work [PR04, PR13], it is useful to have the following in mind.

Let $\alpha$ be a closed smooth $1$-form rapidly decaying at the cusps but not necessarily harmonic. We define unitary operators

[TABLE]

We also define

[TABLE]

where the automorphic Laplacian $\tilde{L}({\varepsilon})$ is the closure of the operator $\Delta$ acting on smooth functions $f$ with $f,\Delta f\in L^{2}({\Gamma\backslash{\mathbb{H}}},\overline{\chi}(\cdot,\epsilon))$. This ensures that $L({\varepsilon})$ acts on the fixed space $L^{2}({\Gamma\backslash{\mathbb{H}}})$. It is then straightforward to verify that

[TABLE]

We notice that $L({\varepsilon})$ does not depend on the cusp ${\mathfrak{a}}$. We observe that $\delta({\alpha})=0$, if ${\alpha}$ is harmonic. From now on we assume $\alpha$ to be harmonic. We remark that $L({\varepsilon})$ as a function of ${\varepsilon}$ is a polynomial of degree 2, and that $L^{(1)}$ and $L^{(2)}$ defined in (4.1) and (4.2) are the first and second derivative of $L({\varepsilon})$ at ${\varepsilon}=0$.

The eigenvalue equation for $E_{{\mathfrak{a}}}(z,s,{\varepsilon})$ and (2.13) imply that

[TABLE]

A formal differentiation of this eigenvalue equation leads to the formula in Lemma 4.2 below. This lemma is the main ingredient in understanding the meromorphic continuation of $D^{(k)}_{\mathfrak{a}}(z,s)$:

Lemma 4.2**.**

The function $D^{(k)}_{\mathfrak{a}}(z,s)$ satisfies the relation

[TABLE]

when $k\geq 1$, where the last term should be omitted for $k=1$.

Proof.

Since $\alpha=\Re(f(z)dz)$ we get

[TABLE]

so that $dA_{\mathfrak{a}}(z)=2\pi i\alpha$. Using the product rule we have that

[TABLE]

where we have used that $(\Delta+s(1-s))\Im(\sigma^{-1}_{\mathfrak{a}}z)^{s}=0$. We recognize the first term as

[TABLE]

The other two terms give

[TABLE]

This proves that

[TABLE]

Now we automorphize this equation over ${\gamma}\in{\Gamma}_{\mathfrak{a}}\backslash{\Gamma}$, and use that $\Delta T_{\gamma}=T_{\gamma}\Delta$. We notice that for any two differential forms $\omega_{1}$, $\omega_{2}$ we have for the action of $\Gamma$ on functions

[TABLE]

Using that $\alpha$ is invariant under $\Gamma$, we arrive at the result. ∎

From Lemma 4.2 we find that, if $L^{(1)}D^{(k-1)}_{\mathfrak{a}}(z,s)$ and $L^{(2)}D^{(k-2)}_{\mathfrak{a}}(z,s)$ are square integrable, then

[TABLE]

where $R(s)$ is the resolvent of $\tilde{\Delta}$. We recall that

[TABLE]

is a bounded operator on $L^{2}({\Gamma\backslash{\mathbb{H}}},d\mu(z))$ satisfying (3.3) and

[TABLE]

where $P_{0}$ is the projection to the eigenspace of $\tilde{\Delta}$ for the eigenvalue [math], and $R_{0}(s)$ is holomorphic in a neighborhood of $s=1$.

Before we can prove the meromorphic continuation of $D^{(k)}_{\mathfrak{a}}(z,s)$ we need a small technical lemma about Eisenstein series:

Lemma 4.3**.**

For $1/2+{\varepsilon}\leq\Re(s)\leq 2$ the function $L^{(1)}E_{\mathfrak{a}}(z,s)$ is square integrable over ${\Gamma\backslash{\mathbb{H}}}$. More precisely we have

[TABLE]

Proof.

We have

[TABLE]

where

[TABLE]

are the raising and lowering operators as in [Fay77, Eq. (3)]. It follows that

[TABLE]

Applying this with $h=E_{\mathfrak{a}}(z,s)$ we see that it suffices to study $y\overline{f(z)}K_{0}E_{\mathfrak{a}}(z,s)$ and $yf(z)L_{0}E_{\mathfrak{a}}(z,s)$. We analyze the first and notice that the analysis of the latter is similar. We have $K_{0}E_{{\mathfrak{a}}}(z,s)=sE_{{\mathfrak{a}}}(z,s,2)$, where $E_{{\mathfrak{a}}}(z,s,2)$ is the Eisenstein series of weight $2$, see [Roe66, Eq. (10.8), (10.9)]. Using that $f$ decays exponentially at all cusps the claim now follows easily from Lemma 3.1.

∎

We are now ready to prove the meromorphic continuation of $D_{{\mathfrak{a}}}^{(k)}(z,s)$.

Theorem 4.4**.**

Let $k>0$. The function $D^{(k)}_{\mathfrak{a}}(z,s)$ admits meromorphic continuation to $\Re(s)\geq 1/2+{\varepsilon}$. For $1/2+{\varepsilon}\leq\Re(s)\leq 2$ and $s(1-s)\not\in{\operatorname{spec}}(-\tilde{\Delta})$ the functions $D^{(k)}_{\mathfrak{a}}(z,s)$ and $L^{(1)}D^{(k)}_{\mathfrak{a}}(z,s)$ are smooth and square integrable. Moreover, we have

[TABLE]

Proof.

We use induction on $k$. For $k=1$ Lemma 4.3, Eq. (4.4), and the mapping properties of the resolvent imply that $D_{\mathfrak{a}}^{(1)}(z,s)=-R(s)L^{(1)}E_{\mathfrak{a}}(z,s)$ is meromorphic when $\Re(s)\geq 1/2+{\varepsilon}$ and, using (3.3), we easily prove the bound $\left\lVert D_{\mathfrak{a}}^{(1)}(z,s)\right\rVert_{L^{2}}\ll 1$. Recall that for a twice differentiable function $h\in L^{2}({\Gamma\backslash{\mathbb{H}}})$ with $\Delta h\in L^{2}({\Gamma\backslash{\mathbb{H}}})$ we have $K_{0}h,L_{0}h\in L^{2}({\Gamma\backslash{\mathbb{H}}})$ and

[TABLE]

see [Roe66, Satz 3.1]. Since $L^{(1)}E_{{\mathfrak{a}}}(z,s)$ is in $C^{\infty}({\Gamma\backslash{\mathbb{H}}})$, we use Lemma 4.2 for $k=1$ and elliptic regularity to conclude that $D_{{\mathfrak{a}}}^{(1)}(z,s)\in C^{\infty}({\Gamma\backslash{\mathbb{H}}})$. By (4.7) the function $L^{(1)}D_{{\mathfrak{a}}}^{(1)}(z,s)$ is smooth as well. Using Lemma 4.2 for $k=1$ again we deduce that $\Delta D_{{\mathfrak{a}}}^{(1)}(z,s)\in L^{2}({\Gamma\backslash{\mathbb{H}}})$. We find from (4.7) and the estimate $\left\lVert yf(z)\right\rVert_{\infty}\ll 1$ that it suffices to estimate $\left\lVert K_{0}D^{(1)}_{{\mathfrak{a}}}(z,s)\right\rVert_{L^{2}}$ and $\left\lVert L_{0}D^{(1)}_{{\mathfrak{a}}}(z,s)\right\rVert_{L^{2}}$. We use (4.8) and Lemma 4.2 to conclude

[TABLE]

where we have used Cauchy-Schwarz and Lemma 4.3. This proves the case $k=1$.

Assume now that the claim has been proved for every $l<k$. Then for $\Re(s)\geq 1/2+{\varepsilon}$ the function $\binom{k}{1}L^{(1)}D^{(k-1)}_{\mathfrak{a}}(z,s)+\binom{k}{2}L^{(2)}D^{(k-2)}_{\mathfrak{a}}(z,s)$ is meromorphic, smooth, and square integrable. Hence, by (4.4), the mapping properties of the resolvent, and (3.3), we find that $D^{(k)}_{\mathfrak{a}}(z,s)$ is square integrable and satisfies the bound

[TABLE]

By elliptic regularity and Lemma 4.2 it follows that $D_{{\mathfrak{a}}}^{(k)}(z,s)\in C^{\infty}({\Gamma\backslash{\mathbb{H}}})$. Therefore, the function $L^{(1)}D_{{\mathfrak{a}}}^{(k)}(z,s)$ is also smooth. We now use (4.7), (4.8), and Lemma 4.2 to see that

[TABLE]

This completes the inductive step. ∎

Lemma 4.3 and Theorem 4.4 validate the conditions for (4.4), so we can now conclude the following fundamental recurrence relation for $D^{(k)}_{{\mathfrak{a}}}(z,s)$.

Corollary 4.5**.**

For $\Re(s)>1/2$ the following identity holds

[TABLE]

when $k\geq 1$, where the last term should be omitted for $k=1$.

Proposition 4.6**.**

For all $k\geq 0$ and $\Re(s)>1/2$ we have

[TABLE]

Proof.

Applying Stokes’ theorem, as in e.g. [PR04, p. 1026], we find that $L^{(1)}$ is a self-adjoint operator. Alternatively, we notice that $L^{(1)}$ is the infinitesimal variation of the family of self-adjoint operators $L({\varepsilon})$ given in Remark 4.1. Therefore, for $k\geq 1$ we have

[TABLE]

since $L^{(1)}$ is a differentiation operator, see (4.1). The case $k=0$ involves $E_{{\mathfrak{a}}}(z,s)$, which is not square integrable. This is easily compensated by the fact that $\alpha$ is cuspidal. ∎

From Corollary 4.5 it is evident that $D_{\mathfrak{a}}^{(k)}(z,s)$ has singularities at the spectrum of $\Delta$. We now describe the nature of the pole at $s=1$.

Proposition 4.7**.**

The functions $L^{(1)}E_{\mathfrak{a}}(z,s)$ and $D^{(1)}_{\mathfrak{a}}(z,s)$ are regular at $s=1$.

Proof.

We have

[TABLE]

We note that $L^{(1)}D_{\mathfrak{a}}^{(0)}(z,s)$ is regular since $D_{\mathfrak{a}}^{(0)}(z,s)=E_{\mathfrak{a}}(z,s)$ has a constant residue and $L^{(1)}$ is a differentiation operator. It then follows from (4.5) that $D^{(1)}_{\mathfrak{a}}(z,s)$ can have at most a simple pole at $s=1$. By (4.5) and Proposition 4.6 we see

[TABLE]

It follows that $D^{(1)}_{\mathfrak{a}}(z,s)$ is regular at $s=1$. ∎

Theorem 4.8**.**

Let $\left\lVert f\right\rVert$ denote the Petersson norm of $f$. If $k$ is even, then $D^{(k)}_{\mathfrak{a}}(z,s)$ has a pole at $s=1$ of order $1+k/2$. The leading term in the corresponding expansion around $s=1$, that is, the coefficient of $(s-1)^{-k/2-1}$ is the constant

[TABLE]

If $k$ is odd, then $D^{(k)}_{\mathfrak{a}}(z,s)$ has a pole at $s=1$ of order at most $(k-1)/2$.

Proof.

The case $k=0$ simply describes the well-known pole and residue of the standard Eisenstein series. For $k=1$ the result is Proposition 4.7.

To run an inductive argument we assume that the claim has been proved up to some $k$ odd. Then $k+1$ is even and we see from (4.5) that $-R(s)L^{(1)}D^{(k)}_{\mathfrak{a}}(z,s)$ can have at most a pole of order $(k-1)/2+1$. The leading term comes from

[TABLE]

which is zero by Proposition 4.6. So $-R(s)L^{(1)}D^{(k)}_{\mathfrak{a}}(z,s)$ has pole of order at most $(k-1)/2$.

On the other hand by inductive hypothesis $-R(s)L^{(2)}D^{(k-1)}_{\mathfrak{a}}(z,s)$ has a pole or order $1+(k-1)/2+1=1+(k+1)/2$ with leading coefficient

[TABLE]

Here we have again used (4.5).

It follows from Corollary 4.5 that $D^{(k+1)}_{{\mathfrak{a}}}(z,s)$ has a pole of order $1+(k+1)/2$ with leading coefficient

[TABLE]

We observe that $\left\langle 1,L^{(2)}1\right\rangle_{L^{2}}=-8\pi^{2}\left\lVert f\right\rVert^{2}$. The order of the pole and leading singularity of $D_{{\mathfrak{a}}}^{(k+1)}(z,s)$ agrees with the claim of the theorem.

We now prove that $D^{(k+2)}_{{\mathfrak{a}}}(z,s)$ has at most a pole of order $(k+1)/2$. We use Corollary 4.5 for $k+2$. Since $L^{(1)}$ annihilates the leading term in $D^{(k+1)}_{{\mathfrak{a}}}(z,s)$, as it is a constant, see (4.10), the function

[TABLE]

can have at most a pole of order $(k+1)/2+1$. This order of singularity at $s=1$ is attained only if $\left\langle 1,L^{(1)}D^{(k+1)}_{{\mathfrak{a}}}(z,s)\right\rangle$ is not identically zero. But it is indeed identically zero by Proposition 4.6. Hence, $-R(s)L^{(1)}D^{(k+1)}_{{\mathfrak{a}}}(z,s)$ has at most a pole of order $(k+1)/2$. By (4.5) and the inductive hypothesis on $D^{(k)}_{{\mathfrak{a}}}(z,s)$ it is straightforward that $-R(s)L^{(2)}D^{(k)}_{{\mathfrak{a}}}(z,s)$ has at most a pole of order $(k+1)/2$. This concludes the inductive step. ∎

5. Functional equations

Selberg’s theory of Eisenstein series [Iwa02, p. 84-94] gives that $E(z,s)$ satisfies the functional equation

[TABLE]

where the scattering matrix $\Phi(s)=(\phi_{{{\mathfrak{a}}}{{\mathfrak{b}}}}(s))$ is determined by (2.5). Recall also that

[TABLE]

In this section we show that $D^{(k)}(z,s)$ and $E^{(k)}(z,s)$ have analogous properties.

Recall the weighted $L^{2}$-spaces in [Mül96, p. 573]. We choose a smooth function $\rho:{\Gamma\backslash{\mathbb{H}}}\to{\mathbb{R}}$ with $\rho(z)=1$ for $z\in F(Y)$ and $\rho(z)=\Im(\sigma_{{\mathfrak{a}}}^{-1}z)$ for $z\in F_{{\mathfrak{a}}}(Y+1)$. For $\delta\in{\mathbb{R}}$ we define

[TABLE]

The resolvent of the Laplace operator $R(s)=(\tilde{\Delta}+s(1-s))^{-1}$ defined on $L^{2}({\Gamma\backslash{\mathbb{H}}})$ for $\Re(s)>1/2$ and $s(1-s)\not\in{\operatorname{spec}}(-\tilde{\Delta})$ admits meromorphic continuation to $\mathbb{C}$ if we restrict the domain to a smaller function space. Müller in [Mül96, Theorem 1] showed that

[TABLE]

can be defined as a bounded operator on weighted $L^{2}$-spaces for $s$ away from its poles. This is achieved by first continuing meromorphically the resolvent kernel (automorphic Green’s function) $r(z,z^{\prime},s)$ to $\mathbb{C}$. The analytic continuation of the resolvent kernel $r(z,z^{\prime},s)$ to $\mathbb{C}$ satisfies the limiting absorption principle:

[TABLE]

see [Lan85, p. 352].

We choose $0<\delta<2\pi$ so that $e^{\delta\rho(z)}y|f(z)|$ is decaying exponentially at the cusps.

Lemma 5.1**.**

Let $k\geq 0$. The function $D^{(k)}_{\mathfrak{a}}(z,s)$ admits meromorphic continuation to ${\mathbb{C}}$. Moreover, we have

1. (i)

$D_{{\mathfrak{a}}}^{(k)}(z,s)\in C^{\infty}({\Gamma\backslash{\mathbb{H}}})$, 2. (ii)

$D^{(k)}_{\mathfrak{a}}(z,s)\in L^{2}_{-\delta}({\Gamma\backslash{\mathbb{H}}})$, 3. (iii)

$K_{0}D^{(k)}_{{\mathfrak{a}}}(z,s)$, $L_{0}D^{(k)}_{{\mathfrak{a}}}(z,s)\in L^{2}_{-\delta}({\Gamma\backslash{\mathbb{H}}})$, and 4. (iv)

$L^{(1)}D^{(k)}_{\mathfrak{a}}(z,s)\in L^{2}_{\delta}({\Gamma\backslash{\mathbb{H}}})$.

Proof.

The Eisenstein series twisted by modular symbols $E_{\mathfrak{a}}^{m,n}(\sigma_{\mathfrak{b}}z,s)$ in (2.11) has Fourier expansion

[TABLE]

with

[TABLE]

see Jorgenson and O’Sullivan [JO08, Eq. (2.4)]. Here $\delta_{0,0}(m,n)=1$ if $m=m=0$ and is [math] otherwise. We quote their work [JO08, Thm 2.2] for the meromorphic continuation of the series $E_{{\mathfrak{a}}}^{m,n}(z,s)$ for $s\in{\mathbb{C}}$. Since $E^{(k)}_{{\mathfrak{a}}}(z,s)$ is a linear combination of $E^{m,n}_{{\mathfrak{a}}}(z,s)$ for $m+n=k$, it admits meromorphic continuation to ${\mathbb{C}}$. The function $D^{(k)}_{{\mathfrak{a}}}(z,s)$ is related to the $E^{(k-j)}_{{\mathfrak{a}}}(z,s)$ for $j\leq k$ through (2.14). Therefore, $D^{(k)}_{{\mathfrak{a}}}(z,s)$ admits meromorphic continuation to ${\mathbb{C}}$ as well.

Furthermore, we need bounds for the Fourier coefficients of $E^{m,n}_{{\mathfrak{a}}}(z,s)$ for all cusps. Jorgenson and O’Sullivan [JO08, Thm 2.3] proved the following: for $s$ in a compact set $S$, there exists a holomorphic function $\xi^{m,n}(s)$ such that for all $k\neq 0$ we have

[TABLE]

It is easy to see that the derivatives of $W_{s}(y)$ decay exponentially in $y$: we can use the integral representation of the $K$-Bessel function [GR07, 8.432.1, p. 917] to see that

[TABLE]

or, alternatively, use repeatedly [GR07, 8.486.11, p. 929]. Combining with (5.3) we see that, for every $a,b\in{\mathbb{N}}\cup\{0\}$, the function $E^{m,n}_{{\mathfrak{a}}}(z,s)$ is smooth and

[TABLE]

Since $E^{(k)}_{{\mathfrak{a}}}(z,s)$ is a linear combination of $E^{m,n}_{{\mathfrak{a}}}(z,s)$ for $m+n=k$, its derivatives satisfy the same upper bounds. By (2.14) the functions $D^{(k)}_{{\mathfrak{a}}}(z,s)$ are also smooth. This proves claim (i).

For the claims (ii), (iii) we need also to control the derivatives $\partial^{a+b}/\partial z^{a}\partial\bar{z}^{b}$ of $A_{{\mathfrak{a}}}(\sigma_{{\mathfrak{b}}}z)^{j}$. For $a+b=0$, we use (2.2). For $a+b\geq 1$ we have exponential decay by (4.3). We conclude that

[TABLE]

The claim (ii) follows by taking $a=b=0$, and the claim (iii) by noticing that $K_{0}$ and $L_{0}$ in (4.6) are expressed in terms of $\partial/\partial x$ and $\partial/\partial y$. Finally, claim (iv) follows from (4.7) using the cuspidality of $f(z)$ at all cusps. ∎

5.1. Functional Equation for $D_{{\mathfrak{a}}}^{(k)}(z,s)$

Theorem 5.2**.**

The meromorphic continuation of the vector-valued automorphic function $D^{(k)}(z,s)$ satisfies the functional equation

[TABLE]

with $\Psi^{(j)}(s)$ meromorphic matrices given by

[TABLE]

if $k>0$ and $\Psi^{(0)}(s)=\Phi(s)$ is the standard scattering matrix.

Proof.

For notational purposes we define $D^{(-1)}(z,s)$ to be [math]. Using Corollary 4.5 we get

[TABLE]

for $k\geq 1$, and $\Re(s)>1/2$. We can extend the validity of this equation to ${\mathbb{C}}$, since the resolvent is applied to a function belonging to $L^{2}_{\delta}({\Gamma\backslash{\mathbb{H}}})$, which follows from Lemma 5.1. Furthermore we see that $D^{(k)}_{\mathfrak{a}}(z,s)$ is holomorphic outside the poles of $R(s)$.

The proof of (5.4) is a relatively obvious generalization of [PR13, Prop. 3.1]. To justify the arguments below we quote Lemma 5.1. We proceed by induction. For $k=0$ the claim of the theorem is the standard functional equation for the vector of Eisenstein series.

Assume the result for $m<k$. We define the matrix $\Psi^{(k)}(s)$ indexed by the cusps by (5.4). Then

[TABLE]

where we have used (5.2) and (5.5). For the first term above we now use the inductive hypothesis and finally (5.5) to see that it equals

[TABLE]

This completes the inductive step. ∎

Selberg proved [Sel89, Eq. (8.5)–(8.6), p. 655] that for $\Re(s)>1/2$ with $s(1-s)$ bounded away from the spectrum the function $\Phi_{{\mathfrak{a}}{\mathfrak{b}}}(s)$ is bounded. We now show how this generalizes to $\Psi^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s)$.

Lemma 5.3**.**

Fix $k\in{\mathbb{N}}$. For $1/2+{\varepsilon}\leq\Re(s)\leq 2$ with $s(1-s)$ bounded away from ${\operatorname{spec}}(-\tilde{\Delta})$ the function $\Psi^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s)$ is bounded, i.e.

[TABLE]

Proof.

Considering (5.4) it suffices to show that for every $k$ we have

[TABLE]

For (5.6) we use Stokes’ theorem, as in the proof of Proposition 4.6, and bound

[TABLE]

By Theorem 4.4, Lemma 4.3, and Cauchy–Schwarz this is bounded by a constant times $\left\lvert s\right\rvert$.

For (5.7) we note that we can move $L^{(2)}$ in front of $E_{{\mathfrak{b}}}(z,s)$, since $L^{(2)}$ is a multiplication operator. Since $\left\lVert L^{(2)}E_{\mathfrak{b}}(z,s)\right\rVert\ll 1$ by Lemma 3.1, and $\left\lVert D^{(k)}_{\mathfrak{b}}(z,s)\right\rVert\ll 1$, see Theorem 4.4, the result follows using Cauchy–Schwarz. ∎

5.2. Functional equation for $E_{{\mathfrak{a}}}^{(k)}(z,s)$

Let $k$ be a natural number. Using (2.15) and (2.14) we can find the functional equation for $E^{(k)}(z,s)$ as follows:

[TABLE]

Setting $r=j+h+l$ we see that

[TABLE]

If we set

[TABLE]

then

[TABLE]

We can rewrite this to see that

[TABLE]

We emphasize that $\Phi^{(r)}_{{\mathfrak{a}}{\mathfrak{b}}}(s)$ does not depend on $z$. We also note that if there is only one cusp we have $\Psi^{(r)}(s)=\Phi^{(r)}(s)$.

Looking at the ${\mathfrak{a}}$-entry of (5.8) and its Fourier expansion (2.10) at the cusp ${\mathfrak{b}}$ we get for the zero Fourier coefficients:

[TABLE]

This gives $\phi_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s)=\Phi_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s)$. We summarize the results for $E^{(k)}(z,s)$.

Theorem 5.4**.**

The vector of Eisenstein series twisted by modular symbols $E^{(k)}(z,s)$ satisfies the functional equation

[TABLE]

where

[TABLE]

and the $\Psi_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s)$ are given by (5.4). Moreover, $\Phi^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s)=\phi^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s)$, where $\phi^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s)y^{1-s}$ is the zero Fourier coefficient of $E^{(k)}_{{\mathfrak{a}}}(\sigma_{{\mathfrak{b}}}z,s)$, see (2.10).

Remark 5.5*.*

The functional equations for $E^{(k)}(z,s)$ can be deduced also from the functional equations for $E^{m,n}_{{\mathfrak{a}}}(z,s)$, see [JO08, Thm. 7.1]. However, the explicit expressions for $\Phi_{{\mathfrak{a}}{\mathfrak{b}}}^{(r)}(s)$ in Theorem 5.4 are new. In this work we need the integral representation of $\Psi_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s)$ in Section 6 below.

The matrices $\Phi^{(k)}(s)$ satisfy the functional equation

[TABLE]

cf. [JO08, Th. 2.2]. For $k=0$ this is due to Selberg and for $k\geq 1$ it follows by comparing the coefficient of $y^{s}$ in (5.9). We notice also that for $1/2+{\varepsilon}\leq\Re(s)\leq A$ with $s$ bounded away from the spectrum we have

[TABLE]

as follows from Lemma 5.3.

6. Analytic properties of the generating series $L_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s,0,n)$

In this section we use the results from Sections 3, 4, 5 to study the analytic continuation of $L_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s,0,n)$ for $k\geq 0$ and $\Re(s)>1/2$.

6.1. Meromorphic continuation

If $k=0$ we obtained the meromorphic continuation of $L_{{\mathfrak{a}}{\mathfrak{b}}}^{(k)}(s,0,n)$ in Lemma 3.5. For $k\geq 1$ we consider first the case $n=0$. From (2.8), and Theorem 5.4 we find that

[TABLE]

The right-hand side is meromorphic by Theorem 5.2. If $n\geq 1$ we deduce from Lemma 3.4 and (2.13) that, for $\Re(s)>1$, $\Re(w)>1$,

[TABLE]

where

[TABLE]

We differentiate to get

[TABLE]

The differentiation is allowed and the right-hand side is meromorphic for $\Re(s)>1/2+{\varepsilon}$ by Theorem 4.4, Lemma 3.3 and the fact that $D_{{\mathfrak{b}},n}^{(k_{3})}(z,w)$ is bounded for $\Re(w)>1$. Using Proposition 2.6 we can deal also with $n\leq-1$. To summarize we have proved the following result:

Theorem 6.1**.**

For any cusps ${\mathfrak{a}},{\mathfrak{b}}$ and any integers $k\geq 0$, $n\in{\mathbb{Z}}$ the function $L^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ admits meromorphic continuation to $\Re(s)>1/2+{\varepsilon}$.

6.2. The first derivative.

We now study in more detail the analytic properties of $L^{(1)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$.

Theorem 6.2**.**

The function $L^{(1)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ has a simple pole at $s=1$ with residue

[TABLE]

For $s(1-s)$ bounded away from ${\operatorname{spec}}(-\tilde{\Delta})$, and $1/2+{\varepsilon}\leq\Re(s)\leq 1+\varepsilon$ we have

[TABLE]

Proof.

Using Proposition 2.6 the claim for $n\leq-1$ follows from the case $n\geq 1$. So we can assume that $n\geq 1$. Consider (6.2) when $k=1$. For $\Re(w)\geq 1+2{\varepsilon}$ fixed, the functions $F(s,w,n)$ is holomorphic as long as $\Re(s)>0$, so we must analyze the three expressions

[TABLE]

To analyze (6.3) we get by Lemma 3.3

[TABLE]

There is a pole of the Eisenstein series $D_{{\mathfrak{a}}}(z,s)$ at $s=1$, which gives rise to a residue for (6.3)

[TABLE]

To see that this vanishes we unfold the integral as in the Rankin method. The integrand contains the factor $e(n\sigma^{-1}_{{\mathfrak{b}}}z)$, with $n\neq 0$, and we notice that, as $\gamma$ varies over the cosets $\Gamma_{{\mathfrak{b}}}\backslash\Gamma$, the sets $\sigma_{{\mathfrak{b}}}^{-1}\gamma F$ cover the strip $\{z\in{\mathbb{H}}:\Re(z)\in[0,1]\}.$

To analyze (6.4) we note that by Theorem 4.4, Proposition 4.7 and (3.5) the term is holomorphic at $s=1$ and, for $s$ bounded away from the spectrum, satisfies

[TABLE]

Finally, we analyze (6.5). Since the Eisenstein series $D_{{\mathfrak{a}}}(z,s)$ has a pole at $s=1$ we find that the integral has a simple pole at $s=1$ with residue

[TABLE]

where we have unfolded using (2.12) and (2.1). We notice also that by Lemma 3.3 we have

[TABLE]

Using the above analysis of the three integrals we can finish the proof for $n>0$ as follows. Observing that $F(1,w,n)={(4\pi n)^{w-1}}/({\Gamma(w-1)\pi})$ we get the residue at $s=1$ for $n>0$. For the growth on vertical lines we choose $\Im(w)=\Im(s)$ and use Stirling’s asymptotics to get

[TABLE]

The bound on vertical lines now is obvious when we choose $\Re(w)=1+2{\varepsilon}$.

As far as $L^{(1)}(s,0,0)$ is concerned we use (6.1) with $k=1$, which leads to analyze $\Psi_{{\mathfrak{a}}{\mathfrak{b}}}^{(0)}(s)$ and $\Psi_{{\mathfrak{a}}{\mathfrak{b}}}^{(1)}(s)$. We start by noticing that by Lemma 5.3 they are both bounded for $1/2+{\varepsilon}\leq\Re(s)\leq 1+{\varepsilon}$. With the help of the Stirling asymptotics on the quotient of Gamma factors we easily prove the bound on vertical lines for $L^{(1)}(s,0,0)$.

Since $\Psi^{(0)}(s)=\Phi(s)$ is the standard scattering matrix it is well known that $\Psi_{{\mathfrak{a}}{\mathfrak{b}}}^{(0)}(s)$ has a simple pole with residue ${\operatorname{vol}}({\Gamma\backslash{\mathbb{H}}})^{-1}$. Using Theorem 5.2 and Proposition 4.7 we see that, if $\Psi_{{\mathfrak{a}}{\mathfrak{b}}}^{(1)}$ has a pole, it must be a simple pole with residue a constant times $\int L^{(1)}D_{\mathfrak{a}}(z,1)d\mu(z)$. This vanishes by Proposition 4.6 so $\Psi_{{\mathfrak{a}}{\mathfrak{b}}}^{(1)}$ is regular at $s=1$. The conclusion follows. ∎

6.3. The second derivative

We will now describe the full singular part of $L^{(2)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,0)$ at $s=1$. We denote the constant term in the Laurent expansion of $E_{\mathfrak{a}}(z,s)$ by $B_{{\mathfrak{a}}}(z)$, i.e.

[TABLE]

as $s\to 1$. For ${\Gamma}={\hbox{PSL}_{2}({\mathbb{Z}})}$, the function $B_{\infty}(z)$ can be described in terms of the Dedekind eta function (Kronecker’s limit formula). For general groups $\Gamma$ the function $B_{{\mathfrak{a}}}(z)$ is given in terms of generalized Dedekind sums, see e.g. [Gol73].

Theorem 6.3**.**

The function $L^{(2)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,0)$ has a pole of order 2 at $s=1$. The full singular part of the Laurent expansion at $s=1$ equals

[TABLE]

where

[TABLE]

For $s(1-s)$ bounded away from ${\operatorname{spec}}(-\tilde{\Delta})$, and $1/2+{\varepsilon}\leq\Re(s)\leq 1+{\varepsilon}$ we have

[TABLE]

Proof.

Using (6.1) we see that $L^{(2)}_{{{\mathfrak{a}}}{{\mathfrak{b}}}}(s,0,0)$ equals

[TABLE]

We consider each of the three terms separately:

We start by noting that since $\Psi^{(0)}(s)=\Phi(s)$, the first term has singular part

[TABLE]

To analyze the second term we note that by Theorem 5.2 we have

[TABLE]

which is regular by Proposition 4.7 and Proposition 4.6.

To analyze the third term we note that by Theorem 5.2

[TABLE]

and analyze the contribution of the two summands. Since $D^{(1)}_{{\mathfrak{a}}}(z,s)$ is regular at $s=1$ (Proposition 4.7), the first summand has at most a first order pole. The corresponding residue is zero by Proposition 4.6, so the first summand is regular.

It follows that the singular part of $\Psi^{(2)}_{{{\mathfrak{a}}}{{\mathfrak{b}}}}(s)$ equals the singular part of ${(2s-1)}^{-1}\int_{{\Gamma\backslash{\mathbb{H}}}}E_{{\mathfrak{b}}}(z,s)L^{(2)}E_{{\mathfrak{a}}}(z,s)d\mu(z)$. It follows that the singular part of the third term of (6.8) equals the singular part of

[TABLE]

The result follows using (6.7), (4.2), and standard values of $\Gamma^{\prime}(z)/\Gamma(z)$ [GR07, 8.366].

The bound on vertical lines follow from Theorem 5.3 and Stirling’s asymptotics. ∎

Remark 6.4*.*

We remark that Theorem 6.3 allows us to write the singular expansion of $L^{(2)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,0)$ exclusively in terms of data of Rankin–Selberg integrals and periods. Indeed, writing

[TABLE]

as $s\to 1$, we have

[TABLE]

so that the singular expansion of $L^{(2)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,0)$ equals

[TABLE]

6.4. Higher derivatives

Theorem 6.5**.**

If $k$ is even the function $L^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,0)$ has a pole at $s=1$ of order $k/2+1$. The leading term in the singular expansion around $s=1$ equals

[TABLE]

If $k$ is odd the function $L^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,0)$ has a pole at $s=1$ of order less than or equal to $(k-1)/2+1$. For $s(1-s)$ bounded away from ${\operatorname{spec}}(-\tilde{\Delta})$, and $1/2+{\varepsilon}\leq\Re(s)\leq 1+{\varepsilon}$ we have

[TABLE]

Proof.

By (6.1) we must understand the leading expansion of each $\Psi^{(h)}_{{\mathfrak{a}}{\mathfrak{b}}}(s)$ for $h\leq k$. The claim about the order of the pole for all $k$, and the leading singularity for $k$ even follows from (5.4), Theorem 4.8, and Proposition 4.6.

The claim on bounds on vertical lines follow from (6.1), Stirling’s formula, and Lemma 5.3. ∎

Theorem 6.6**.**

Let $n\neq 0$. Then we have:

1. (i)

The function $L^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ has a pole at $s=1$ of order strictly less than $[k/2]+1$. 2. (ii)

For $s(1-s)$ bounded away from ${\operatorname{spec}}(-\tilde{\Delta})$, and $1/2+{\varepsilon}\leq\Re(s)\leq 1+{\varepsilon}$ we have

[TABLE] 3. (iii)

All coefficients in the singular expansion of $L^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ are bounded independently of $n$.

Proof.

Using Proposition 2.6 it suffices to treat the case $n\geq 1$. Considering (6.2) we see that claim (i) about the orders of the pole follows from Theorem 4.8 and the fact that $\int_{{\Gamma\backslash{\mathbb{H}}}}\overline{D_{{{\mathfrak{b}}},n}(z,\bar{w})}\,d\mu(z)=0$ as is seen by unfolding.

Claim (ii) follows from the bound (6.6) (with $\Re(w)=1+2{\varepsilon}$) valid when $\Im(w)=\Im(s)$ combined with Lemma 3.3, Theorem 4.4 and the fact that $D_{{\mathfrak{b}},n}^{(k_{3})}(z,w)$ is bounded independently of $w$ and $n$ for $1+2{\varepsilon}\leq\Re(w)\leq A$, see the discussion after (3.5).

For claim (iii) we note that the constants in all singular expansions are linear combinations of

[TABLE]

where $g(z)$ is one of the coefficients in the singular expansion of $D_{{\mathfrak{a}}}^{(k_{2})}(z,s)$. For $k_{2}=0$ the function $g(z)$ is constant so, in particular, is square integrable. For $k_{2}>0$ we note that

[TABLE]

for some $j\geq 0$ and sufficiently small $r$. Here $C(r)$ is the circle centered at $1$ with radius $r$. The radius $r$ is chosen so that there are no other singularities of $D^{(k_{2})}(z,s)$ inside $C(r)$ apart from $s=1$. It follows from Theorem 4.4 that $g(z)$ is square integrable. By using Cauchy–Schwarz we see that (6.9) is bounded independently of $n$. See again the discussion after (3.5). ∎

7. Distribution results

We are now ready to prove Theorems 1.9, 1.10 and 1.11. Since we have identified the behavior of the generating functions $L^{(k)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,0,n)$ at $s=1$ and on vertical lines in Section 6, we can use the well-known method of contour integration to deduce the asymptotics of $\langle r\rangle_{{\mathfrak{a}}{\mathfrak{b}}}$.

7.1. First moment with restrictions.

In this subsection we study sums of the form

[TABLE]

for smooth functions $h$ or indicator functions $h=1_{[0,x]}$. Hence we are studying a (partial) first moment of the modular symbol but with restrictions on $r$ imposed by $h$.

Remark 7.1*.*

We present a variant of the Mazur–Rubin–Stein heuristics: By Theorem 3.8 $T_{{\mathfrak{a}}{\mathfrak{b}}}$ is equidistributed on ${{\mathbb{R}}/{\mathbb{Z}}}$. If it had been possible to extend the function $h(r)=\langle r\rangle_{{\mathfrak{a}}{\mathfrak{b}}}1_{[0,x]}(r)$ to a continuous function of $r$, this would give the asymptotics of (7.1) immediately. Using (2.1) it would be tempting to define the modular symbol for all $r\in{\mathbb{R}}/{\mathbb{Z}}$ by

[TABLE]

We cannot do so as the series is not convergent, even if Wilton’s classical estimate, see [Wil29], [Iwa97, Thm 5.3], shows that it just barely fails to converge conditionally.

By termwise integration against $1_{[0,x]}$ we would get the result

[TABLE]

This series converges to a continuous function as is easily seen from Hecke’s average bound, [Iwa97, Thm 5.1].

If instead we consider, for a fixed $\delta>0$,

[TABLE]

then this function does indeed define a continuous function on ${\mathbb{R}}/{\mathbb{Z}}$, and we can use equidistribution with$\langle r\rangle_{{\mathfrak{a}}{\mathfrak{b}},\delta}1_{[0,x]}(r)$ as a test function. If we do so, and then let $\delta\to 0$ we arrive again at (7.2). However, it is not easy to justify that one can interchange the limits $M\to\infty$ and $\delta\to 0$. On the other hand Mazur, Rubin and Stein have numerics suggesting that (7.2) is indeed the correct limit.

The above heuristics gives the correct answer. This is the content of Theorem 7.2 below. For a formal series

[TABLE]

we have a linear functional (distribution) $h\mapsto\langle h,F\rangle_{{{\mathbb{R}}/{\mathbb{Z}}}}$ from the set of smooth functions on ${{\mathbb{R}}/{\mathbb{Z}}}$ given by

[TABLE]

where $\hat{h}(n)$ denotes the $n$th Fourier coefficient of $h$.

Recall the norm (3.7). We are now ready to prove the main result of this section:

Theorem 7.2**.**

Let $h$ be a smooth function on ${{\mathbb{R}}/{\mathbb{Z}}}$ with $\left\lVert h\right\rVert_{H^{1/2}}<\infty$. Then there exists a $\delta>0$ such that

[TABLE]

where $F_{{\mathfrak{a}}{\mathfrak{b}}}$ is the formal series given by

[TABLE]

Proof.

The generating series of $\langle r\rangle_{{\mathfrak{a}}{\mathfrak{b}}}e(nr)$ is $L^{(1)}_{{\mathfrak{a}}{\mathfrak{b}}}(s,n,0)$. Writing $h(t)=\sum_{n\in{\mathbb{Z}}}\hat{h}(n)e(nt),$ we have by Proposition 2.6, Theorem 6.2, and a complex integration argument that

[TABLE]

from which the result follows. ∎

Let $0\leq x\leq 1$. Approximating $1_{[0,x]}$ by smooth periodic functions we can conclude the following result, which makes rigorous the heuristic conclusions in Remark 7.1.

Corollary 7.3**.**

Let $x\in[0,1]$. There exists a $\delta>0$ such that

[TABLE]

7.2. The variance

In this subsection we study the second moment of the modular symbols, i.e. the variance. Following Mazur and Rubin we denote the variance slope by

[TABLE]

Recall (6.7). We also define the variance shift by

[TABLE]

The following theorem follows directly from Theorem 6.3 and a complex integration argument.

Theorem 7.4**.**

There exists a $\delta>0$ such that

[TABLE]

We deduce from Theorem 7.4 and Theorem 3.8 that

[TABLE]

7.3. Normal distribution

In this subsection we show that the value distribution of modular symbols (appropriately normalized) obeys a standard normal distribution, even if we restrict $r$ to any interval.

Theorem 7.5**.**

Let $h$ be a function on ${{\mathbb{R}}/{\mathbb{Z}}}$ satisfying $\left\lVert h\right\rVert_{H^{\varepsilon}}<\infty$, and let $k\in{\mathbb{N}}$. Then there exist $\delta,B>0$ such that

[TABLE]

Proof.

We use Proposition 2.6, Theorem 6.5, and Theorem 6.6. We apply a complex integration argument in a strip of width ${\varepsilon}$ around $\Re(s)=1$ to deduce that

[TABLE]

We insert this in

[TABLE]

and use $\hat{h}(0)=\int_{{\mathbb{R}}/{\mathbb{Z}}}h(t)dt$ to get the result. ∎

Remark 7.6*.*

The result (7.5) can be strengthened to the following: There exists a polynomial $P_{k,m}$ such that

[TABLE]

The degree of $P_{k,m}$ is strictly less than $k/2$ if either $m\neq 0$ or $k$ is odd, and exactly $k/2$ for $k$ even and $m=0$.

In Theorem 7.4 we identify this polynomial when $k=2$ and $m=0$.

Using a standard approximation argument based in Theorem 7.5 we arrive at the following corollary:

Corollary 7.7**.**

Let $I\subseteq{\mathbb{R}}/{\mathbb{Z}}$ be an interval and let $k\in{\mathbb{N}}$. Then

[TABLE]

The above corollary allows us to renormalize the modular symbol map and determine the distribution of the renormalized map using the method of moments:

Corollary 7.8**.**

Let $I\subseteq{\mathbb{R}}/{\mathbb{Z}}$ be an interval of positive length. Then the values of the map

[TABLE]

ordered according to $c(r)$ have asymptotically a standard normal distribution, i.e. for every $-\infty\leq a\leq b\leq\infty$ we have

[TABLE]

as $M\to\infty$.

Proof.

Using summation by parts we find from Corollary 7.7 and Theorem 3.8 that

[TABLE]

which is the $k$th moment of the standard normal distribution. The result follows from a classical result due to Fréchet and Shohat [Loè77, 11.4.C].

∎

Remark 7.9*.*

From the above proof we see that the asymptotic moments do not change if you replace $C_{f}\log c(r)$ with $C_{f}\log c(r)+D$ for any constant $D$. Therefore Corollary 7.8 also holds if we normalize accordingly.

8. Results for Hecke congruence groups

In this section we translate the distribution results of Section 7 to the case of Hecke congruence groups ${\Gamma}={\Gamma}_{0}(q)$, where $q$ is a squarefree integer. In this case the cusps of ${\Gamma}$ and their scaling matrices can be described as follows, see [DI82, Section 2.2]): a complete set of inequivalent cusps of ${\Gamma}_{0}(q)$ are given by ${\mathfrak{a}}_{d}=1/d$ with $d|q$. Notice that if $d=q$ then $1/d$ is equivalent to the cusp at infinity. Write $q=dv$. We may take

[TABLE]

for the corresponding scaling matrix. It follows that

[TABLE]

Using this and definition (1.5) we easily see that

[TABLE]

and for $r={a}/{c}\in T_{\infty\frac{1}{d}}$ we have

[TABLE]

Therefore,

[TABLE]

where $\langle r\rangle$ is as in (1.2).

8.1. First moment with restrictions

From Theorem 7.2 and Theorem 3.8 we now deduce the following corollary.

Corollary 8.1**.**

Let $d|q$, and let $h$ be a smooth function on ${{\mathbb{R}}/{\mathbb{Z}}}$ with $\left\lVert h\right\rVert_{H^{1/2}}<\infty$. Then there exists $\delta>0$ such that

[TABLE]

where $F$ is the formal series

[TABLE]

with $a(n)$ the Fourier coefficients at infinity of $f(z)$.

Summing Corollary 8.1 over all positive divisors $d|q$ and using that

[TABLE]

we may remove the divisibility condition on $c$ and conclude that

[TABLE]

We can also remove the condition $(a,c)=1$, by summing according to $(a,c)=k$. Using that $\zeta(2)=\pi^{2}/6$ we find that

[TABLE]

Using partial summation we find that

[TABLE]

By an approximation argument, where we approximate $h(t)=1_{[0,x]}(t)$ by appropriate smooth functions, we find

[TABLE]

as $M\to\infty$. This completes the proof of Theorem 1.4.

8.2. The variance

Similarly to the analysis above we can use Theorem 7.4, Corollary 7.3, and Theorem 3.8 to conclude the following corollary.

Corollary 8.2**.**

There exists a $\delta>0$ such that

[TABLE]

Recall now the mean and variance from (1.3).

Lemma 8.3**.**

Assume $f$ is a Hecke eigenform. Then ${\operatorname{E}(f,c)}\ll c^{-1/2+{\varepsilon}}$.

Proof.

Since $f$ is an eigenfunction of all Hecke operator $T_{n}$ with eigenvalue $a(n)$ it is easy to see that for any rational $r$ we have

[TABLE]

If $r$ is a fixed integer, then we deduce that

[TABLE]

Using Möbius inversion and the Eichler bound for weight $2$ holomorphic cusp forms on congruence groups [Eic54], i.e. $a_{n}\ll n^{1/2+{\varepsilon}}$, we find that

[TABLE]

Another application of Möbius inversion, and the well-known lower bound $\phi(c)^{-1}\ll c^{-1+{\varepsilon}}$ [HW79, Thm. 329] give ${\operatorname{E}(f,c)}\ll c^{-1/2+{\varepsilon}}$. ∎

We define the variance shift by

[TABLE]

Using Lemma 8.3 we see that

[TABLE]

Using this and Corollary 8.2 we deduce that

[TABLE]

as $M\to\infty$. Here we have used (3.6) for the asymptotics of the last sum.

8.2.1. Relation with the symmetric square $L$-function

We explain how to relate $C_{f}$, and $D_{f,d}$ to the symmetric square $L$-function. We recall the definitions but refer to [ILS00, Sec. 2-3] for additional details. Assume that $f$ is a Hecke eigenform normalized with first Fourier coefficient equal to 1, and let $\lambda_{f}(n)$ be its $n$th Hecke eigenvalues. Let

[TABLE]

be the Rankin–Selberg $L$-function. It is known that $L(f\otimes f,s)$ admits meromorphic continuation to $s\in{\mathbb{C}}$ with a simple pole at $s=1$ with corresponding residue ${(4\pi)^{2}\left\lVert f\right\rVert^{2}}/{{\operatorname{vol}}({\Gamma}_{0}(q)\backslash{\mathbb{H}})}$. We have

[TABLE]

where

[TABLE]

and $\zeta_{(q)}(s)$ is the Riemann zeta function with the Euler factors at $p|q$ removed. The symmetric square $L$-function of $f$ is defined by

[TABLE]

Using these definitions, the formula ${\operatorname{vol}}({\Gamma}_{0}(q)\backslash{\mathbb{H}})=({\pi}/{3})q\prod_{p|q}(1+p^{-1})$, and that $\zeta(s)$ has a simple pole at $s=1$ with residue $1$, we find that

[TABLE]

This verifies the variance slope of Conjecture 1.2.

To express in more arithmetic terms the constant $D_{f,d}$ we notice that

[TABLE]

is the constant term in Laurent expansion at $s=1$ of

[TABLE]

by definition of $B_{\infty}(z)$. By unfolding we see that the last integral equals

[TABLE]

Hence we conclude that

[TABLE]

where

[TABLE]

To understand the part of $D_{d}$ involving $B_{1/d}$ we recall the Atkin–Lehner involutions [AL70]. We also recall that we assume that $q$ is squarefree. For every $d|q$ there exists an integer matrix of determinant $d$ of the form

[TABLE]

with $y,w\in{\mathbb{Z}}$. It is straightforward to verify that $W_{d}$ normalizes ${\Gamma}_{0}(q)$, that is, $W_{d}^{-1}{\Gamma}_{0}(q)W_{d}={\Gamma}_{0}(q)$. Since $f$ is assumed to be a Hecke eigenform it follows from Atkin–Lehner theory that

[TABLE]

with $e_{f,d}=\pm 1$ the Atkin–Lehner eigenvalues. It follows easily that

[TABLE]

Lemma 8.4**.**

The Atkin–Lehner involutions permute the Eisenstein series. More precisely, for every $d|q$ we have

[TABLE]

Proof.

Let $\sigma_{1/d}^{\prime}=\frac{1}{\sqrt{q/d}}W_{q/d}$. Then a direct computation shows that

[TABLE]

so $\sigma^{\prime}_{1/d}$ is an admissible scaling matrix for the cusp $1/d$. Since the Eisenstein series $E_{1/d}(z,s)$ is independent of the choice of scaling matrix we find that

[TABLE]

where we have used that $W_{q/d}$ normalize ${\Gamma}_{0}(q)$.

∎

Combining Lemma 8.4 and (8.10) we find that

[TABLE]

It follows that

[TABLE]

Using the definition of $G(s)$ i.e. (8.9), we see that

[TABLE]

Combining (8.5), (8.6), (8.8), and (8.11) we find the following expression for the variance shift. We have

[TABLE]

where

[TABLE]

This completes the proof of Theorem 1.6.

8.3. Numerical investigations

As an example we consider the elliptic curve 15.a1. We computed $L({\operatorname{sym}}^{2}f,1)=0.9364885435$ and $L^{\prime}({\operatorname{sym}}^{2}f,1)=0.03534541$ using lcalc in Sage [Sag17] and the data from [LMF17]111See lmfdb.org/L/SymmetricPower/2/EllipticCurve/Q/15.a/. These numbers should be accurate to at least $6$ decimal places. This allows to estimate the value of the variance shift $D_{f,d}$ and compare with the experimental values in Table 1. We would like to thank Karl Rubin for providing us with the experimental values of the variance shift. Notice that these are the opposite of what appears in [MR16], because our modular symbols are purely imaginary.

8.4. Normal distribution

With the description of the modular symbols and $T_{\infty\frac{1}{d}}$ from (8.2), (8.3), and (8.4), we can apply Corollary 7.7 to compute the moments of $\langle r\rangle$. Since $r=\frac{a\sqrt{q/d}}{c\sqrt{q/d}}\in T_{\infty\frac{1}{d}}$ has $c(r)=c\sqrt{q/d}$, we may use $\log c(r)=\log c+2^{-1}\log(q/d)$ and Remark 7.9 to conclude the following corollary.

Corollary 8.5**.**

Let $I\subseteq{\mathbb{R}}/{\mathbb{Z}}$ be any interval of positive length, and consider for $d|q$ the set $Q_{d}=\{{a}/{c}\in{\mathbb{Q}},(a,c)=1,(c,q)=d\}$. Then the values of the map

[TABLE]

ordered according to $c$ have asymptotically a standard normal distribution.

This completes the proof of Theorem 1.7.