Cotangent sums, quantum modular forms, and the generalized Riemann hypothesis

TL;DR

This paper links the asymptotic behavior of cotangent sum matrices to the generalized Riemann hypothesis (GRH) for odd Dirichlet characters and explores their connection to quantum modular forms associated with Maass Eisenstein series.

Contribution

It establishes an equivalence between cotangent sum matrix asymptotics and the GRH, and connects this to quantum modular forms related to Maass Eisenstein series.

Findings

Asymptotic properties of cotangent sum matrices are equivalent to GRH for odd Dirichlet characters.

Connection between cotangent sums and quantum modular forms is demonstrated.

Provides a new perspective on GRH through quantum modular forms.

Abstract

We show that an asymptotic property of the determinants of certain matrices whose entries are finite sums of cotangents with rational arguments is equivalent to the GRH for odd Dirichlet characters. This is then connected to the existence of certain quantum modular forms related to Maass Eisenstein series.

Equations

$$\chi (n)=\{\begin{array}{ll}(-1{)}^{(n-1)/2}& \text{if n is odd,}\\ 0& \text{if n is even,}\end{array}L(s)=1-\frac{1}{3^s}+\frac{1}{5^s}-\frac{1}{7^s}+\cdots .$$

$$c_{m,n}:=\frac{4}{\pi }\sum _{j,k>0}\frac{\chi (j)\chi (k)}{\max (j/m,k/n)},$$

$$c_{3,4}\;=\;\cot\Bigl{(}\frac{\pi}{12}\Bigr{)}\,+\,\cot\Bigl{(}\frac{5\pi}{12}\Bigr{)}\,+\,\cot\Bigl{(}\frac{9\pi}{16}\Bigr{)}+\cot\Bigl{(}\frac{13\pi}{16}\Bigr{)}\;=\;6\,-\,2\,\sqrt{2+\sqrt{2}}\;.$$

$${\sf C}_{N}\;:=\;\begin{pmatrix}c_{1,1}&\cdots&c_{1,N}\\ \vdots&\ddots&\vdots\\ c_{N,1}&\cdots&c_{N,N}\end{pmatrix}\;,\qquad\widehat{\sf C}_{N}\;:=\;\left(\begin{array}[]{ccc|c}&&&1\\[-4.0pt] &{\sf C}_{N}&&\vdots\\ &&&N\\ \hline\cr 1&\cdots&N&4N/\pi\\[1.0pt] \end{array}\right)\,,$$

$$R(N):=N\frac{\det {\mathsf{C}}_N}{\det {\hat{\mathsf{C}}}_N}.$$

$$c_{m,n}\;=\;\frac{4n}{\pi}\,\,C\Bigl{(}\frac{m}{n}\Bigr{)}\,,\qquad C(x)\;:=\;\sum_{j,\,k>0}\chi(j)\,\chi(k)\,\min\Bigl{(}\frac{x}{j},\,\frac{1}{k}\Bigr{)}\qquad(x\in\mathbb{Q}_{>0})\,.$$

$$c_{m,n}\;=\;h_{m,n}\,+\,h_{n,m}\,,\qquad C(x)=H(x)\,+\,x\,H\Bigl{(}\frac{1}{x}\Bigr{)}$$

$$H(x)\;:=\;\sideset{}{{}^{*}}{\sum}_{0<j\leq kx}\frac{\chi(k)\,\chi(j)}{k}\,,\qquad h_{m,n}\;:=\;\frac{4n}{\pi}\,H\Bigl{(}\frac{m}{n}\Bigr{)}\,.$$

$$H(x)+|x|H(1/x)=C(x)(x\ne 0).$$

$$H(x)+H(x+2)=\frac{\pi }{4},$$

$$\varepsilon(\gamma)\,|cx+d|\,H\Bigl{(}\frac{ax+b}{cx+d}\Bigr{)}\;=\;\,H(x)\;+\;\text{\rm(continuous function)}$$

$$H(x)=\sum _{k=1}^{\infty }\frac{\chi (k)}{k}S(kx)(x\in \mathbb{Q}),$$

$$S(x):=\text{sideset}{}^{\ast }{\sum }_{0<j\le x}\chi (j).$$

$$S(x)=\{\begin{array}{ll}1& \text{if 4k+1<x<4k+3 for some k∈Z,}\\ 0& \text{if 4k−1<x<4k+1 for some k∈Z,}\\ \frac{1}{2}& \text{if x is an odd integer,}\end{array}$$

$$c_{m,n}=\frac{4}{\pi }{\int }_0^{\infty }S(mt)S(nt)\frac{dt}{t^2}(m,n\in \mathbb{N}).$$

$$\int_{0}^{N}\!S(mt)\,S(nt)\,\frac{dt}{t^{2}}\;=\;\int_{0}^{N}\!\sideset{}{{}^{*}}{\sum}_{0<j\leq mt\atop 0<k\leq nt}\chi(j)\chi(k)\,\frac{dt}{t^{2}}\;=\;\sum_{0<j\leq mN\atop 0<k\leq nN}\biggl{(}\frac{\chi(j)\,\chi(k)}{\max(j/m,k/n)}\,-\,\frac{\chi(j)\,\chi(k)}{N}\biggr{)}\,,$$

$$c_{m,n}\,=\sum_{0\leq j\leq m/2\atop 0\leq k\leq n/2}\begin{cases}\cot(\pi a)-\cot(\pi b)&\text{if $a:=\max\bigl{(}\frac{4j+1}{4m},\frac{4k+1}{4n}\bigr{)}<b:=\min\bigl{(}\frac{4j+3}{4m},\frac{4k+3}{4n},\frac{1}{2}\bigr{)}\,,$}\\ 0&\text{if $a\geq b$\,.}\end{cases}$$

$${\mathsf{C}}_4=\left(\begin{array}{cccc}1& 2-\sqrt{2}& 2-\sqrt{3}& 2-2\sqrt{2-\sqrt{2}}\\ 2-\sqrt{2}& 2& \sqrt{2}& 4-2\sqrt{2}\\ 2-\sqrt{3}& \sqrt{2}& 3& 6-2\sqrt{2+\sqrt{2}}\\ 2-2\sqrt{2-\sqrt{2}}& 4-2\sqrt{2}& 6-2\sqrt{2+\sqrt{2}}& 4\end{array}\right).$$

$$c_{m,n}=\pi {\int }_0^{1/2}S(4mx)S(4nx)\frac{dx}{{\sin }^2(\pi x)}.$$

$$h_{m,n}\;=\;\sum_{k=1}^{2n-1}\chi(k)\,S\Bigl{(}\frac{km}{n}\Bigr{)}\,\cot\Bigl{(}\frac{\pi k}{4n}\Bigr{)}\;.$$

$$h_{m,n}$$

$$\displaystyle\;=\;\frac{1}{2}\,\sum_{0<k<4n}\chi(k)\,S\Bigl{(}\frac{km}{n}\Bigr{)}\,\cot\Bigl{(}\frac{\pi k}{4n}\Bigr{)}\;=\;\sum_{0<k<2n}\chi(k)\,S\Bigl{(}\frac{km}{n}\Bigr{)}\,\cot\Bigl{(}\frac{\pi k}{4n}\Bigr{)}$$

$$R(N)\le \frac{\pi }{8}\frac{|\rho {|}^2}{\Re (\rho )-\frac{1}{2}}(\forall N\ge 1).$$

$$G(x)=S(1/x)=\{\begin{array}{ll}1& \text{if 4k+31<x<4k+11 for some k∈Z≥0,}\\ 0& \text{otherwise}\end{array}$$

$${\tau }_{a}f(x)=f(x/a)(0<a\le 1,x\in {\mathbb{R}}_{>0})$$

$${\mathcal{V}}_{\mathcal{D}}=⟨G_a|a\in \mathcal{D}\cap (0,1]⟩,{\mathcal{V}}_N=⟨G_{n/N}|1\le n\le N⟩,$$

$${\sf d}(\xi,\mathcal{V}_{N})\;=\;\bigl{\|}\,\xi\,-\,\mathcal{V}_{N}\bigr{\|}\;=\;\inf_{x\in\mathbb{R}^{N}}\,\Bigl{\|}\,\xi\,-\,\sum_{n=1}^{N}x_{n}\,\tau_{n/N}G\,\Bigr{\|}\;.$$

$$\mathsf{d}(\xi ,{\mathcal{V}}_N{)}^2=\frac{\pi /4}{R(N)}(N\in \mathbb{N}).$$

$$\langle\xi,\xi\rangle\,=\,\int_{0}^{1}dx\,=\,1,\qquad\langle\xi,G_{n/N}\rangle\,=\,\int_{0}^{1}S\bigl{(}\frac{n}{Nx}\Bigr{)}\,dx\,=\,\frac{n}{N}\,L(1)\,=\,\frac{\pi n}{4N}$$

$$\langle G_{m/N},\,G_{n/N}\rangle\,=\,\int_{0}^{1}S\bigl{(}\frac{m}{Nx}\Bigr{)}\,S\bigl{(}\frac{n}{Nx}\Bigr{)}\,dx\;=\;\frac{m}{N}\,C\Bigl{(}\frac{m}{n}\Bigr{)}\;=\;\frac{\pi}{4N}\,c_{m,n}$$

$${\mathsf{p}}_s(t)=t^{s-1}\xi (t)=t^{s-1}(0<t\le 1,s\in \mathbb{C}).$$

Full text

Cotangent sums, quantum modular forms,

and the generalized Riemann hypothesis

John Lewis

 and 

Don Zagier

1. Introduction and main results

To answer immediately a question that the reader may be asking, we should say from the outset that the generalized Riemann hypothesis (henceforth, just GRH) of the title refers only to the $L$-series associated to odd real Dirichlet characters and that we do not aim to prove it, but merely to show its equivalence to an asymptotic statement about the determinants of certain matrices whose entries $c_{m,n}$ are given by finite sums of cotangents of rational multiples of $\pi$. To answer a second natural question, since it may seem unusual for an author to contribute to the proceedings volume of a conference partly intended to celebrate his own birthday, we should also say a few words about the genesis of this joint paper. Both the definition of the cotangent sums considered and their relation to GRH were found by the first author and presented in his talk at the conference in question. In the course of that talk he mentioned the homogeneity property $c_{\ell m,\ell n}=\ell c_{m,n}$, leading the second author to ask whether the function $C:\mathbb{Q}\to\mathbb{R}$ defined by $C(m/n)=c_{m,n}/n$ might be related to a quantum modular form in the sense of [16]. The answer to this question turned out to be positive. Furthermore, and unexpectedly, this quantum modular aspect turned out also to be related to the earlier joint work of the authors [13] on the analogues of classical period polynomials for Maass wave forms. A glance at the title of the final section of this paper will show why it seemed natural, and even irresistible, to publish both parts of the story in these proceedings.

We now describe the main results of the paper in more detail. We denote by

[TABLE]

the primitive Dirichlet character of conductor 4 and its associated $L$-series. (In Section 5 we will describe the generalization to other odd primitive real Dirichlet characters.) For positive integers $m$ and $n$ we define a real number $c_{m,n}$ by

[TABLE]

where the conditionally convergent sum is to be interpreted as the limit for $N\to\infty$ of the sum over the rectangle $1\leq j\leq mN$, $1\leq k\leq nN$. These numbers are studied in more detail in §2, where we will see that they are algebraic numbers (in fact, algebraic integers) that can be expressed as finite sums of cotangents of odd multiples of $\pi/4m$ or $\pi/4n$, a typical example being

[TABLE]

For each $N\in\mathbb{N}$ we define the two symmetric matrices

[TABLE]

which we will show later are positive definite, and define a function $R:\mathbb{N}\to\mathbb{R}$ by

[TABLE]

The first main result of this paper is then the following.

Theorem 1**.**

The following two statements are equivalent.

• (a)

The function $R(N)$ is unbounded.

• (b)

The $L$-series $L(s)$ has no zeros in the half-plane $\sigma:=\Re(s)>\frac{1}{2}$ (GRH).

The proof of this theorem will be given in §3. The direction (a)$\,\Rightarrow\,$(b) is fairly elementary, while the reverse direction uses functional analysis and techniques coming from old work of Beurling. The result that we actually prove (Theorem 3 in §3) is in fact stronger than Theorem 1 in both directions: if GRH holds, then $R(N)$ is not merely unbounded, but actually tends to infinity as $N\to\infty$, and if $L(s)$ has a zero $\rho$ with $\Re(\rho)>\frac{1}{2}$, then the function $R(N)$ is bounded by a number depending explicitly on $\rho$. Of course we believe that the two statements in Theorem 1 are both true rather than both false, the evidence being on the one hand the generally held belief in the validity of (b) and on the other hand the following plot of the function $R(N)$ for $N\leq 100000$, which suggests that $R(N)$ tends to infinity and perhaps even grows roughly linearly with $\log N$. (The straight-line fit shown in the picture is the graph of the function $\,5.18\,\log N\,$.)

Figure 1. Graph of the function $R(N)$

The second main result concerns the quantum modular nature of certain functions related to the coefficients $c_{m,n}$. On the one hand, the numbers $c_{m,n}$ defined by (2) have the obvious homogeneity property $c_{\ell m,\ell n}=\ell\,c_{m,n}$ and can be rewritten in the form

[TABLE]

On the other hand, by splitting the sums in (2) or (5) into two pieces depending on which of its two arguments realizes the “max” or “min”, we obtain the decompositions

[TABLE]

where the function $H:\mathbb{Q}_{>0}\to\mathbb{R}$ and numbers $h_{m,n}\in\mathbb{R}$ are defined by

[TABLE]

(Here and later an asterisk on a summation sign means that terms with equality—in this case, those with $j=kx$—are to be counted with multiplicity 1/2.) We will show in Section 2 that $h_{m,n}$ (or its sum with 1/2 if $mn$ is odd) is a sum of cotangents of odd multiples of $\pi/4n$ and hence is an algebraic integer, refining the corresponding properties for $c_{m,n}$ given above.

We can now come to the quantum modular form property mentioned in the opening paragraph. The graph of the function $H$ on $\mathbb{Q}$ looks as follows:

Figure 2. Graph of the function $H(x)$

This graph suggests that the function $H$ is only naturally defined on $\mathbb{Q}$ and is not well-behaved as a function on $\mathbb{R}$, having a possibly dense set of discontinuities and certainly not differentiable at generic points. The graph of the function $C$ shown in Figure 3, on the other hand, suggest that it is everywhere continuous and possibly even differentiable except at rational points with odd numerator and denominator, where it seems to be left and right differentiable but with a discontinuous derivative. Our second main result, proved in Section 4, refutes the last of these statements and confirms the others.

Figure 3. Graph of the function $C(x)$

Theorem 2**.**

The function $H(x)$ is discontinuous at all rational points $x$ with odd numerator and denominator, and is continuous but not left or right differentiable at other rational points. The function $C(x)$ is continuous everywhere, is not left or right differentiable at rational points with odd numerator and denominator, but is differentiable at all other rational points.

The similarity between the continuity and differentiability statements for the two functions $H(x)$ and $C(x)$ in this theorem is not coincidental, since we will see in Section 4 that the derivative of $C(x)$ is equal to $H(1/x)$.

We now discuss the “quantum modular form” aspect of Theorem 2. Note first that the function $H$ is periodic of period 4 and hence extends to a periodic function on all of $\mathbb{Q}$. This function is even, so we could also extend $H$ from $\mathbb{Q}_{>0}$ to $\mathbb{Q}$ by $H(x):=H(|x|)$. One might then think of extending $C$ to negative values of the arguments by setting $C(x)=H(x)+xH(1/x)$ for all $x\neq 0$, but if we did that then $C(x)$ for $x<0$ would exhibit oscillatory behavior of the same sort as $H(x)$. Instead, we extend $C$ to all of $\mathbb{Q}$ as an even function, i.e., we set $C(x)=C(|x|)$ for $x\neq 0$ and extend this by continuity to $C(0)=0$. Then the relationship between $C$ and $H$ becomes

[TABLE]

On the other hand, the 4-periodicity property of $H(x)$ can be strengthened to the statement

[TABLE]

saying that up to sign and up to a constant $H$ actually has period 2. Also, since $H$ is even, we can replace the term $H(1/x)$ in (8) by $H(-1/x)$. But the two matrices S=\bigl{(}\begin{smallmatrix}0&-1\\ 1&0\end{smallmatrix}\bigr{)} and T^{2}=\bigl{(}\begin{smallmatrix}1&2\\ 0&1\end{smallmatrix}\bigr{)} generate the subgroup $\Gamma_{\vartheta}$ of index 3 in the full modular group $\,\text{SL}(2,\mathbb{Z})=\langle S,T\rangle$ consisting of matrices congruent to \bigl{(}\begin{smallmatrix}1&0\\ 0&1\end{smallmatrix}\bigr{)} or \bigl{(}\begin{smallmatrix}0&1\\ 1&0\end{smallmatrix}\bigr{)} modulo 2 (this is the so-called “theta group,” under which the Jacobi theta function $\vartheta(z)=\sum_{n}e^{\pi in^{2}z}$ transforms like a modular form of weight 1/2), so equations (8) and (9) can be combined to the following statement:

Corollary**.**

The function $H:\mathbb{Q}\to\mathbb{R}$ satisfies the transformation property

[TABLE]

for every matrix \gamma=\bigl{(}\begin{smallmatrix}a&b\\ c&d\end{smallmatrix}\bigr{)}\in\Gamma_{\vartheta}, where $\varepsilon:\Gamma_{\vartheta}\to\{\pm 1\}$ is the homomorphism sending both generators \bigl{(}\begin{smallmatrix}1&2\\ 0&1\end{smallmatrix}\bigr{)} and \bigl{(}\begin{smallmatrix}0&-1\\ 1&0\end{smallmatrix}\bigr{)} of $\Gamma_{\vartheta}$ to $-1$.

This corollary tells us precisely that the function $H$ is a quantum modular form in the sense of [16] (except that there is an absolute value sign here that was not present there). Recall that a “quantum modular form” is by definition a function on $\mathbb{Q}$ that does not extend nicely to $\mathbb{R}$ but which modulo (the restriction from $\mathbb{R}$ to $\mathbb{Q}$ of) better-behaved (here, continuous) functions on $\mathbb{R}$ transforms like a modular form with respect to some subgroup of finite index of $\text{SL}(2,\mathbb{Z})$ (here the theta group $\Gamma_{\vartheta}$) with some weight (here $-1$) and character (here $\varepsilon$). In Section 6 of this paper we will show that the quantum modularity property of $H$ is related to an actual modular form (specifically, a Maass Eisenstein series of eigenvalue $1/4$) on the group $\Gamma_{\vartheta}$.

The authors would like to thank Steven Johnson (MIT) and Alexander Weisse (MPI) for invaluable help with the numerical aspects of this paper.

2. Properties of the numbers $c_{m,n}$ and $h_{m,n}$

In this section we give an integral representation for the numbers $c_{m,n}$ that will be used later both to extend the function $C(x)$ continuously from $\mathbb{Q}$ to $\mathbb{R}$ and as the key to the proof in Section 3 of the relation between the numbers $c_{m,n}$ and the zeros of the Dirichlet $L$-series $L(s)$. We will also prove the formulas for $c_{m,n}$ and $h_{m,n}$ as algebraic numbers (finite sums of cotangents).

The first step is to rewrite the definition of $H$ given in (7) more compactly in the form

[TABLE]

where $S(x)$ is the function defined on $\mathbb{R}_{>0}$ by

[TABLE]

The function $S$, like $H$ itself, is periodic of period 4, and we can use this property to extend it to all of $\mathbb{R}$. (Equivalently, we can extend it by $S(-x)=S(x)$ and $S(0)=0$, or by observing that the summation condition in (11) can be replaced by $J<j\leq x$ for any integer $J<x$ divisible by 4, which works for all real $x$.) The function $S$ is then the step function given by

[TABLE]

as pictured in Figure 4. Note that the sum in (10) converges (conditionally) for $x\in\mathbb{Q}$ by summation by parts, because the function $\chi(k)S(kx)$ is then an odd periodic function of $k$ and hence has average value zero, so that the sums $\sum_{k=1}^{K}\chi(k)S(kx)$ are bounded.

$-1$[math]1$$2$$3$$4$$5$$\tfrac{1}{2}$$1

Figure 4. Graph of the function $S(x)$

We use the bounded periodic function $S$ to prove the following formula for the $c_{m,n}\,$:

Proposition 1**.**

The numbers $c_{m,n}$ defined by (2) have the integral representation

[TABLE]

Proof. From the definition of $S$ we find

[TABLE]

and since the second sum is identically zero for $N$ divisible by 4, equation (13) follows by letting $N$ tend to infinity in $4\mathbb{Z}$. $\quad\square$

Our next result expresses $c_{m,n}$ as a finite sum of algebraic cotangent values indexed by lattice points in the rectangle $[0,m/2]\times[0,n/2]$ satisfying a certain “selection rule.”

Proposition 2**.**

Let $m$ and $n$ be positive integers. Then the number defined by (2) is given as the following integral linear combination of cotangents of rational multiples of $\pi$:

[TABLE]

Corollary**.**

The numbers $c_{m,n}$ are algebraic integers for all $m,\,n>0$.

Using this proposition we can calculate any $c_{m,n}$ in closed form as an algebraic number. For example, the first $4\times 4$ matrix of values (cf. equation (3)) is given by

[TABLE]

Proof. Since $S(x)$ is even and periodic of period 4, we can use Euler’s formula $\sum_{k\in\mathbb{Z}}(k+t)^{-2}=\pi^{2}/\sin^{2}(\pi t)$ to rewrite the integral representation (13) as

[TABLE]

But from the definition of $S(4mx)$ and $S(4nx)$ as characteristic functions, it follows that the restriction of $S(4mx)S(4nx)$ to $[0,\frac{1}{2}]$ is the characteristic function of the union of all intervals $[a,b]$ with $a<b$ as given in the proposition. This proves the proposition because $\pi\int_{a}^{b}\csc^{2}(\pi x)\,dx=\cot(\pi a)-\cot(\pi b)$. The corollary then follows because by elementary algebraic number theory the number $\cot(\pi\alpha)$ is an algebraic integer for every rational number $\alpha$ having an even denominator. $\quad\square$

We next give the corresponding formula for the numbers $h_{m,n}$, which is in fact simpler.

Proposition 3**.**

Let $m$ and $n$ be positive integers. Then the number defined by (7) is given as the following half-integral linear combination of cotangents of odd multiples of $\pi/4n$:

[TABLE]

Proof. Since the summand in (10) is an even function of $k$, we can rewrite the sum as a sum over all non-zero integral values of $k$, and then use the periodicity (with period at most $4n$) of the factor $\chi(k)S(km/n)$ together with Euler’s formula $\sum_{n\in\mathbb{Z}}(x+n)^{-1}=\pi\cot(\pi x)$ (with the sum being interpreted as a Cauchy principal value) to get

[TABLE]

as claimed. $\quad\square$

Remark. One can also deduce the formula for $c_{m,n}$ in Proposition 2 from (14); this is left as an exercise to the reader. We also observe that the terms on the right-hand side of (14) (assuming, as we can, that $m$ and $n$ are coprime) are integral multiples of cotangents of odd multiples of $\pi/4n$, and hence are algebraic integers (in some cyclotomic field), except for the $k=n$ term if $m$ and $n$ are odd, which is half-integral because then $S(km/n)=\frac{1}{2}$ and $\cot(\pi k/4n)=1$. Hence $h_{m,n}$ is an algebraic integer unless both $m$ and $n$ are odd, in which case it equals an algebraic integer plus $1/2$, and the sum $c_{m,n}=h_{m,n}+h_{n,m}$ is always an algebraic integer.

3. The connection with the generalized Riemann hypothesis

In the first subsection of this section we outline the basic strategy for proving Theorem 1. The easier direction (a)$\,\Rightarrow\,$(b) of this theorem is then proved in Subsection 3.2 and the more difficult converse direction in Subsection 3.3, using ideas from two classical papers by A. Beurling. As already stated in the introduction, we will actually prove the following stronger result:

Theorem 3**.**

Let $L(s)$ be the $L$-series defined in (1) and $R(N)$ be the function defined in (4).

• (a)

If $L(s)$ has no zeros with $\Re(s)>\frac{1}{2}$, then $R(N)$ tends to infinity as $N$ tends to infinity.

• (b)

If $L(\rho)=0$ for some $\rho\in\mathbb{C}$ with $\Re(\rho)>\frac{1}{2}$, then $R(N)$ satisfies the bound

[TABLE]

3.1. The square-wave function $G$ and its multiplicative shifts

Set $L^{2}=L^{2}(0,1)$ and let $\xi$ be the characteristic function of the interval $(0,1]$ and $G:\mathbb{R}_{>0}\to\mathbb{R}$ the function

[TABLE]

(and $G(x)=1/2$ if $1/x$ is an odd integer, but we are working in $L^{2}$). Clearly both $\xi$ and $G$ belong to $L^{2}$, since they are bounded with support in the interval $(0,1]$. On the other hand, if we identify $L^{2}$ with the space of square-integrable functions on $\mathbb{R}_{>0}$ that vanish on $(1,\infty)$, then we have bounded operators $\tau_{a}$ on $L^{2}$ defined by

[TABLE]

We set $G_{a}=\tau_{a}G$ and for any $\mathcal{D}\subseteq\mathbb{R}$ and $N\in\mathbb{N}$ define subspaces $\mathcal{V}_{\mathcal{D}}$ and $\mathcal{V}_{N}$ of $L^{2}$ by

[TABLE]

where $\langle\,S\,\rangle$ denotes the linear span of a subset $S$ of a linear space. Finally, we denote by ${\sf d}(\xi,\mathcal{V}_{N})$ the $L^{2}$ distance from $\xi\in L^{2}$ to the subspace $\mathcal{V}_{N}$, viz.

[TABLE]

The key connection is then given by the following proposition.

Proposition 4**.**

The distance ${\sf d}(\xi,\mathcal{V}_{N})$ is related to the number $R(N)$ defined in (4) by

[TABLE]

The proof of this, which is a straightforward computation, will be given in §3.2 below. On the other hand, a quite elementary argument, also given in §3.2, proves the following

Proposition 5**.**

If $L(\rho)=0$ with $\Re(\rho)>\frac{1}{2}$ then ${\sf d}(\xi,\mathcal{V}_{N})^{2}\,\geq\,\frac{2\Re(\rho)-1}{|\rho|^{2}}$ for all $N\in\mathbb{N}$ .

Together these two results prove statement (b) of Theorem 3. The harder statement (a) follows from the two following results, proved in §3.2 and §3.3 respectively.

Proposition 6**.**

The space $\mathcal{V}_{\mathcal{D}}$ is dense in $\mathcal{V}_{\mathbb{R}}$ for any set $\mathcal{D}$ that is dense in $(0,1)$.

Proposition 7**.**

If $L(s)$ has no zeros in $\sigma>\frac{1}{2}$, then $\mathcal{V}_{\mathbb{R}}$ is dense in $L^{2}$.

The whole argument is summarized by the following diagram, where the flow of the double arrows shows that the six boxed statements are mutually equivalent.

\boxed{\sf{GRH}}$$\boxed{\xi\in\overline{\mathcal{V}_{\mathbb{Q}}}}$$\boxed{R(N)\text{ is unbounded}}$$\boxed{\overline{\mathcal{V}_{\mathbb{R}}}=L^{2}}$$\boxed{\overline{\mathcal{V}_{\mathcal{D}}}=L^{2}\text{ (\mathcal{D} dense)}}$$\boxed{\lim\nolimits_{N\to\infty}R(N)=\infty}Prop. 5Prop. 4Prop. 7Prop. 4Prop. 6trivialtrivial

(The implications labelled “Prop. 4” follow in one direction because each $\mathcal{V}_{N}$ is contained in $\mathcal{V}_{\mathbb{Q}}$ and in the other because if $R(N)$ were bounded for some infinite sequence of integers $\{N_{i}\}$ then the set $\cup_{i}\mathcal{V}_{N_{i}}$ would be dense in $\mathcal{V}_{\mathbb{R}}$ but have a non-zero distance from $\xi$.)

3.2. Proof of Propositions 4, 5 and 6

Recall from linear algebra that the distance $d(v,V)$ from a vector $v$ in a (real) Hilbert space to a finite-dimensional vector space $V$ with basis $\{v_{1},\dots,v_{N}\}$ is given by the formula $d(v,V)^{2}=\det(\widehat{G})/\det(G)$, where $G$ is the $N\times N$ Gram matrix with $(m,n)$-entry $\langle v_{m},v_{n}\rangle$ (scalar product of the basis elements $v_{m}$ and $v_{n}$) and $\widehat{G}$ the augmented $(N+1)\times(N+1)$ matrix obtained from $G$ by adding a row and column with entries $\langle v,v_{n}\rangle$ and a diagonal entry $\langle v,v\rangle$. We apply this to the vector $v=\xi\in L^{2}$ and the space $\mathcal{V}_{N}$ defined as above as the span of the vectors $G_{n/N}\in L^{2}$ ($n=1,\dots,N$). (Here we mean the real rather than the complex span, i.e., we work in the real Hilbert space of real-valued square-integrable functions on $(0,1)$.) We also denote by $\widehat{\mathcal{V}}_{N}$ the span of $\mathcal{V}_{N}$ and $\xi$. We will show in a moment that the vectors $G_{n/N}$ and $\xi$ are linearly independent, so $\dim\widehat{\mathcal{V}}_{N}=\dim\mathcal{V}_{N}+1=N+1$. We can then compute the entries of the Gram matrices as

[TABLE]

(the first is trivial and the second a special case of equation(21) below) and

[TABLE]

(combine equations (25) and (5)). Proposition 4 follows.

To find the connection with the $L$-series $L(s)$ and its complex zeros, we next consider the “complex monomial functions” ${\sf p}_{s}$ on $(0,1]$ defined by

[TABLE]

Note the key fact that ${\sf p}_{s}\in{L}^{2}$ if and only if $\sigma>\frac{1}{2}$, and also that for $f\in L^{2}$ and $\sigma>\frac{1}{2}$, we have $\langle{\sf p}_{s},f\rangle$ = $\widetilde{f}(s)$, where $\widetilde{f}(s)$ denotes the Mellin transform

[TABLE]

(Here we are again identifying $L^{2}=L^{2}(0,1)$ with the space of functions $L^{2}(\mathbb{R}_{>0})$ supported on $(0,1]$.) Then obviously

[TABLE]

and the Mellin transform of $G_{a}$ is almost equally easy to calculate using (12):

[TABLE]

the calculation being valid for all $s$ in $\Re(s)>0$ since the integral then converges absolutely. One consequence of this is that the functions $G_{a}$ ($0<a\leq 1$) are all linearly independent, since if some (finite) linear combination of them vanished then the vanishing of its Mellin transform would imply that the product of $L(s)$ with a finite Dirichlet series vanished identically (first for $\Re(s)>0$, and then for all $s$ by analytic continuation). In particular, the space $\mathcal{V}_{N}$ has dimension exactly $N$. In fact, almost the same argument also shows that the functions $G_{a}$ and $\xi$ are linearly independent, since a linear relation among them would imply that the product of $L(s)$ with a finite Dirichlet series was constant, which is obviously impossible. It follows that the extended space $\widehat{\mathcal{V}}_{N}$ spanned by the functions $G_{n/N}$ ($1\leq n\leq N$) and $\xi$ has dimension $N+1$, as claimed above.

The proof of the direction (a)$\,\Rightarrow\,$(b) of Theorem 1 now follows easily. If there exists $\rho$ with $\mathrm{Re}\,(\rho)>\frac{1}{2}$ and $L(\rho)=0$ then by (21) we have ${\sf p}_{\rho}\perp\mathcal{V}_{N}$ for all $N$. Since $\mathcal{V}_{N}$ is finite-dimensional, the infimum in equation (18) is attained, so that $R(N)$ is unbounded (and hence ${\sf d}(\xi,\mathcal{V}_{N})$ has lim inf equal to 0) if and only if there is a sequence of functions $h_{N}\in{\mathcal{V}}_{N}$ converging in $L^{2}$ to $\xi$; and this would contradict ${\sf p}_{\rho}\perp\mathcal{V}_{N}$ for all $N$ since $\langle{\sf p}_{\rho},\xi\rangle=\rho^{-1}\neq 0$.

This argument actually proves the stronger statement given in Proposition 5 and in part (b) of Theorem 3. Indeed, if $\rho$ is a zero of $L(s)$ with $\Re(\rho)>\frac{1}{2}$, then on the one hand ${\sf p}_{\rho}$ belongs to $L^{2}$ as already noted, with $L^{2}$ norm

[TABLE]

and on the other hand $\langle{\sf p}_{\rho},f\rangle=0$ for any $f\in\mathcal{V}_{N}$ by virtue of (21). We apply this to the function $f=f_{N}\in\mathcal{V}_{N}$ attaining the minimum distance to $\xi$. Then applying the Cauchy-Schwarz inequality to the scalar product $\langle{\sf p}_{\rho},\xi-f_{N}\rangle$ we get (using (20))

[TABLE]

and (15) now follows immediately from (19).

Finally, we prove the Proposition 6 (“density lemma”). We first note that the $L^{2}$ norm $\|G_{a}-G_{b}\|$ tends to 0 as $a$ tends to $b$, since by the calculation $\langle G_{a},G_{b}\rangle=aC(a/b)$ given above (for rational numbers, but true in general by the same argument) we have

[TABLE]

which tends to 0 as $b\to a$ because $C(x)$ is a continuous function. It follows immediately that if $\mathcal{D}$ is dense in $(0,1)$ then the vector $G_{a}\in L^{2}$ for any $a<1$ lies in the closure of $\mathcal{V}_{\mathcal{D}}$, and since $\mathcal{V}_{\mathbb{R}}$ is spanned by such vectors, this shows that $\overline{\mathcal{V}_{\mathbb{R}}}=\overline{\mathcal{V}_{D}}$ as claimed.

3.3. Proof of Proposition 7

It remains to show that $\sf{GRH}$ implies that $\overline{\mathcal{V}_{\mathbb{R}}}=L^{2}$.

Let $\mathrm{H}^{2}=\mathrm{H}^{2}(\mathcal{R})$ denote the Hardy space of the right half-plane $\mathcal{R}=\{s\mid\Re(s)>\frac{1}{2}\}$, that is, the space of all holomorphic functions $F(s)$ on $\mathcal{R}$ that are square-integrable on the lines $\Re(s)=c$ for all $c>\frac{1}{2}$ and with these $L^{2}$ norms $\|\ \|_{c}$ bounded in $c$, and then equipped with the Hilbert space norm $\|F\|=\sup\{\|F\|_{c}\mid c>\frac{1}{2}\}$. We quote some results from the standard analysis of these spaces, citing as we go the relevant sections in the books by Hoffman [11] and Garnett [8] for the details.

The Paley-Wiener theorem ([11] pp. 131–132) gives the fact that the Mellin transform is (up to scalar multiple) an isometry between the spaces $L^{2}$ and $\mathrm{H}^{2}$. Applying this to the function $G$ implies that its Mellin transform $\,\widetilde{G}(s)=L(s)/s\,$ lies in the space $\mathrm{H}^{2}$, a fact which is also easy to see directly from the definition of $\mathrm{H}^{2}$. Next, we note that if we assume that $\sf{GRH}$ holds, so that $\widetilde{G}(s)$ has no zeros in $\mathcal{R}$, then the factorization theorem for functions in the space $\mathrm{H}^{2}$ ([11] pp. 132–133) applied to the function $\widetilde{G}(s)$ gives the implication

[TABLE]

Indeed, the factorization theorem represents any element $F$ of $\mathrm{H}^{2}$ as the product of a “Blaschke product,” a “singular function,” and an “outer function.” Here we can omit the definitions of all three of these since the Blaschke product is defined as a product over the zeros of $F$ in $\mathcal{R}$ and hence is equal to 1 if $F$ has no such zeros, and the singular function factor is also constant because $\widetilde{G}(s)$ both continues across the line $\Re(s)=\frac{1}{2}$ and has slow rate of convergence to zero as $s\to\infty$ on the positive reals (cf. 3.14 in [1] and Chapter II, Theorem 6.3, of [8]). This leaves only the outer function factor of $\widetilde{G}(s)$, and here we can apply a corollary of a theorem of Beurling on shift-invariant subspaces in $\mathrm{H}^{2}$ ([4], [11] pp. 99–101), which tells us that a function $F\in\mathrm{H}^{2}$ is outer if and only if the space

[TABLE]

is dense in $\mathrm{H}^{2}$. (Notice that $\mathcal{W}(F)$ is contained in $\mathrm{H}^{2}$ because the definition of the $\mathrm{H}^{2}$ norm implies that $\|s^{-1}F(s)\|_{\mathrm{H}^{2}}\leq 2\|F\|_{\mathrm{H}^{2}}<\infty$ for $F\in\mathrm{H}^{2}$, so $s^{-1}\mathrm{H}^{2}\subset\mathrm{H}^{2}\,$.)

It follows that GRH implies that $\mathcal{W}(\widetilde{G})$ is dense in $\mathrm{H}^{2}$. We can now apply the inverse Mellin transform from $\mathrm{H}^{2}$ to $L^{2}$ to obtain the isometrically equivalent picture in $L^{2}$. An easy calculation shows that $s^{-1}\widetilde{f}(s)$ is the Mellin transform of $If$ for any $f\in L^{2}$, where

[TABLE]

and hence that the space $\mathcal{W}(\widetilde{f})$ is the Mellin transform of the space

[TABLE]

It is not actually obvious that $I$ preserves the space $L^{2}$, but this follows from the Paley-Wiener theorem together with the fact just shown that multiplication by $s^{-1}$ preserves $\mathrm{H}^{2}$. A more direct argument is that the operator $I$ is the adjoint ${\sf A}^{\ast}$ of the averaging operator

[TABLE]

as is easily checked, and since ${\sf A}$ is a bounded operator on $L^{2}(0,\infty)$ (cf. [10], Example 5.4, p. 23), its adjoint $I$ is also bounded on $L^{2}(0,\infty)$. But $I$ also preserves the condition of having support in (0,1), so it is a bounded operator on $L^{2}$.

As a consequence the implication (22) translates on the $L^{2}$ side to the statement

[TABLE]

Proposition 7 is then obtained by applying to $f=G$ the following lemma, in which $\mathcal{V}(f)$ for any $f\in L^{2}$ denotes the span of the multiplicative translates $\tau_{a}f$ ($0<a\leq 1$).

Lemma. For $f\in L^{2}$ we have $W(f)\subset\overline{\mathcal{V}(f)}\,$.

Proof. We first note that it suffices just to show that $If\in\overline{\mathcal{V}(f)}$ for all $f\in L^{2}$, because then $I^{n}f\in\overline{\mathcal{V}(f)}$ for all $n\geq 1$ by induction and the remark that $g\in\overline{\mathcal{V}(f)}\implies\overline{\mathcal{V}(g)}\subseteq\overline{\mathcal{V}(f)}$. Since $(W^{\perp})^{\perp}=\overline{W}$ for all subspaces $W$, we need to show that

[TABLE]

But since $a\tau_{a^{-1}}$ is the adjoint of $\tau_{a}$ and ${\sf A}$ is the adjoint of $I$, this statement just says that any function that is orthogonal to all $h(ax)$ with $0<a\leq 1$ is also orthogonal to ${\sf A}h$, and this is obvious because ${\sf A}h(x)=\int_{0}^{1}h(ax)\,da\,$. $\quad\square$

4. Analytic properties of $C(x)$ and $H(x)$

In this section we study the analytic properties of $C(x)$ and $H(x)$, discussing in particular the question of their extendability from $\mathbb{Q}$ to $\mathbb{R}$ (which we answer completely only in the case of $C(x)$) and describing the asymptotic properties in neighborhoods of rational points needed to establish Theorem 2 and its corollary (quantum modularity). We also discuss the relation of these two functions to the Dirichlet series $L(s)$. It should be mentioned in passing that much of the material in this section is similar in spirit to the papers [2] and [3] of Bettin and Conrey, with $\chi$ replaced by the trivial character and $L(s)$ by the Riemann zeta function.

4.1. Sum and integral representations of $C(x)$ and $H(x)$

Proposition 8 below gives a number of formulas for these two functions. However, the word “function” has a slighly different meaning in the two cases. We have already defined them both as functions on $\mathbb{Q}$. Part (i) of the proposition says that $C$ extends continuously to $\mathbb{R}$, and parts (iv) and (v) give further properties of that function. But in the case of $H$ we do not know to what extent it can be defined as a function on $\mathbb{R}$ (this question is discussed in Subsection 4.2, though not in great detail since this is not our main interest). However, it can certainly be defined as a distribution on $\mathbb{R}$, in fact even as the derivative of a continuous function, since by integrating (10) we get the absolutely convergent sum representation

[TABLE]

for its integral, where $A$ is the periodic continuous function on $\mathbb{R}$ defined by $A(n+\epsilon)=(-1)^{n}\epsilon$ for $n\in\mathbb{Z}$ and $|\epsilon|\leq\frac{1}{2}$. Parts (ii) and (iii) of the proposition are then to be interepreted in the sense of distributions, one giving the Fourier expansion of the even periodic distribution $H$ and the other its relation, already mentioned in the introduction, with the derivative of $C$.

Proposition 8**.**

$($i$)$ The function $C(x)$ extends continuously from $\mathbb{Q}$ to $\mathbb{R}$ via the formula

[TABLE]

$($ii$)$ The distribution $H$ has the Fourier expansion

[TABLE]

where $d(n)$ denotes the number of positive divisors of $n$.

$($iii$)$ The function $C(x)$ and the distribution $H(x)$ are related by the equation

[TABLE]

$($iv$)$ The function $C(x)$ is also given by the sum

[TABLE]

where $J(x)$ denotes the modified cosine integral

[TABLE]

$($v$)$ The Mellin transform of $C(x)$ is given by

[TABLE]

where $L(s)$ is the analytic continuation of the $L$-series defined in (1).

Proof. Formula (25) for $x\in\mathbb{Q}_{>0}$ is precisely the statement of Proposition 1 after an obvious change of variables, and the convergence and continuity of the integral for all $x\in\mathbb{R}_{>0}$ are easy consequences of the facts that $S$ is supported on $[1,\infty)$ and is locally constant and bounded. To prove (ii), we insert the (standard) Fourier expansion

[TABLE]

of the even periodic function $S(x)$ into the sum (10), combine the double sum into a single one, and use Leibniz’s formula $L(1)=\pi/4$. For (iii), we insert (11) into (25) to get

[TABLE]

from which equation (27) follows immediately by differentiating and using (10). (The easy justification of these steps in the sense of distributions by integrating against smooth test functions is left to the reader.) Finally, to prove (iv) we combine (27) and (6) to get

[TABLE]

and then insert the Fourier development (26) and integrate term-by-term. This gives (28) up to an additive term $\lambda x$ that can be eliminated by noting that $C$ is bounded. Note that the sum in (28) is uniformly and absolutely convergent since a simple integration-by-parts argument shows that $J(x)=\text{O}(x^{-2})$ as $x\to\infty$. Finally, for (v) we note first that the Mellin transform $\widetilde{S}(s):=\int_{0}^{\infty}S(x)\,x^{s-1}dx$ of $S$ is given according to equation (21) by $\widetilde{S}(s)=-L(-s)/s$ for $\sigma=\Re(s)<0$, and then use (25) to get

[TABLE]

for $-1<\sigma<0$, as claimed. (The Mellin transform exists in this strip because $C(x)=\text{O}(\min(|x|,1))\,$ as an easy consequence of (25).) In particular, $\widetilde{C}(s)$ extends meromorphically to all of $\mathbb{C}$ and is invariant under $s\mapsto-s-1$. We also note that by the well-known functional equation of $L(s)$, equation (30) can be written in the alternative form

[TABLE]

We can then use this to give a second derivation of (28), not making use of the distribution $H$, by first applying the Mellin inversion formula to write $C(x)$ as an integral over a vertical line $\sigma=c$ with $-1<c<0$ and then shifting the path of integration to the right, picking up the term $\pi/8$ in (28) from the residue at $s=0$ and permitting us to use the convergent Dirichlet series representation $L(s)^{2}=\sum_{n=1}^{\infty}\chi(n)d(n)n^{-s}$ for $\sigma>1$. (We omit the details of this calculation.) $\quad\square$

4.2. The function $H$ on the real line

In this subsection we give two one-parameter generalizations of $H(x)$ and discuss possible ways to define this function at irrational arguments.

We first generalize the definitions of $H(x)$ and $C(x)$ and the assertions of Proposition 8 to families depending on a complex parameter $s$. Specifically, we can generalize (10) to

[TABLE]

which now converges absolutely. The same calculation as for $H(x)$ then gives the Fourier expansion

[TABLE]

where $L(s)$ is the Dirichlet $L$-function defined in (1) and $\sigma_{\nu}$ is the sum-of-divisors function $\sigma_{\nu}(n)=\sum_{k|n,k>0}k^{\nu}$. Similarly, we define $C_{s}(x)=x^{s}C_{s}(1/x)$ for all $s\in\mathbb{C}$ with $\sigma>0$ by

[TABLE]

where the equality of the two expressions is proved by the same calculation as for the proof of (25). Then $C_{s}(x)$ and $H_{s}(x)$ are related by

[TABLE]

by the same calculations as in the special case $s=1$ and one has an expansion like (28) but with $J(x)$ replaced by $-\int_{x}^{\infty}t^{-s-1}\cos(t)\,dt$. As long as $\sigma>1$ all of the sums and integrals involved are absolutely convergent and all steps are justified.

A different possible way to regularize $H(x)$ for $x\in\mathbb{R}$ is to define

[TABLE]

Using the identity \sum\limits_{n>0}\dfrac{d(n)}{n}\,x^{n}=\sum\limits_{k>0}\dfrac{1}{k}\,\log\Bigl{(}\dfrac{1}{1-x^{k}}\Bigr{)} and standard trigonometric identities, we find after a short calculation that

[TABLE]

Both this series and the original one defining $T$ converge exponentially fast for any positive $\varepsilon$, giving us another possible approach to the analytic properties of the limiting function $H$.

Summarizing this discussion, we have at least five potential definitions of $H(x)$ for $x\in\mathbb{R}$:

Definition 1: Define $H(x)$ by the series (10), if this sum converges.

Definition 2: Define $H(x)$ by the Fourier series (26), if this series converges.

Definition 3: Define $H(x)$ as the limit for $s\searrow 1$ of the series (33), if this limit exists.

Definition 4: Define $H(x)$ as the limit in (36), if this limit exists.

Definition 5: Define $H(x)$ as the limit of $H(x^{\prime})$ as $x^{\prime}\in\mathbb{Q}$ tends to $x$, if this limit exists.

We can then ask—but have not been able to answer—the question whether any or all of these definitions converge for irrational values of $x$ or, if they do, whether they give the same value. The last definition is in the sense the strongest one, since any definition of $H$ on $\mathbb{R}$ or a subset of $\mathbb{R}$ that agrees with the original definition on $\mathbb{Q}$ must coincide with the value in Definition 5 at any argument $x$ at which this function is continuous. In any case, we pose the explicit question:

Question: Are the five definitions given above convergent and equal to one another for all irrational values of $x$, or for all $x$ belonging to some explicit set of measure 1 ?

It seems reasonable to expect the answer to the second question to be affirmative with the set of measure 1 being the complement of some set of irrational numbers having extremely good rational approximations, like the well-known “Brjuno numbers” in dynamics.

4.3. Asymptotic behavior of $H(x)$ and $C(x)$ near rational points

The results of Section 4.1 prove the continuity of $C(x)$, which is part of the statement of Theorem 2 in the Introduction. We now discuss the remaining statements there, concerning the behavior of $C$ and $H$ near rational points.

We begin by looking numerically at the asymptotic properties of $C(x)$ near $x=1$. Computing the values of C\bigl{(}1\pm\frac{1}{n}\bigr{)} for $1\leq n\leq 1000$ and using a numerical interpolation technique that is explained elsewhere (see e.g. [9]), we find empirically an expansion of the form

[TABLE]

with the first few coefficients having numerical values given by

[TABLE]

These numerics become clearer if we work with $H(x)$ instead, where we find the simpler expansion

[TABLE]

with no odd powers of $1/n$ and with coefficients given numerically by

[TABLE]

which then gives the expansion (37) with $c_{1}^{\pm}=-h_{0}^{\mp}-\dfrac{1}{4}$, $c_{3}=h_{2}+\dfrac{1}{24}$, and more generally $c_{r}=(-1)r\sum_{0<i<r/2}\binom{r-1}{2i-1}h_{2i}$ for $r>1$. Moreover, on calculating the next few values of $h_{2i}$ numerically to high precision we are able to guess the closed formulas

[TABLE]

for the coefficients in (38), where $\gamma$ is Euler’s constant and $B_{n}$ the $n$th Bernoulli number. In fact, it is not difficult to prove (38) (and hence also (37)) with these values of $h_{i}$ by applying a twisted Euler-Maclaurin formula (giving the asymptotics as $\varepsilon\to 0$ of sums over an interval of $\chi(n)f(n\varepsilon)$ for a smooth function $f$) to the closed formula

[TABLE]

which follows easily from (14), though we will not carry this out here.

The “$\log n$” terms in equations (37) and  (38) already show that $C(x)$ is not differentiable at $x=1$ and $H(x)$ is not continuous at $x=1$, but from these formulas one might imagine that the functions $C(x)-|x-1|\log(|x-1|)/4$ and $H(x)-\log(|x-1|)/4$ are $C^{\infty}$ to both the right and the left of this point. However, this is not the case, because both equations (37) and (38) are valid only when $n$ is an integer and change in other cases. For instance, if $n$ tends to infinity in $\mathbb{N}+\frac{1}{2}$ rather than $\mathbb{N}$, then $H(1\pm 1/n)$ has an expansion of the same form as (38) and with the same constants $h_{0}^{\pm}$, but with $h_{2}=\pi^{2}/576$ replaced by $-7\pi^{2}/1152$, $h_{4}=-49\pi^{4}/1382400$ replaced by $127\pi^{4}/2764800$, etc.) This statement, which again can be proved using an appropriate twisted version of the Euler-Maclaurin formula, is a typical phenomenon of quantum modular forms.

If we look at the asymptotics of $C(x)$ and $H(x)$ as $x$ tends to any rational number $\alpha$ with odd numerator and denominator, then we find a similar behavior, which is not surprising since any such $\alpha$ is $\Gamma_{\vartheta}$-equivalent to 1 and both $H$ and $C$ have transformation properties, modulo functions with better smoothness properties, under the action of $\Gamma_{\vartheta}$. More precisely, if we write $\alpha=a/c$ with $a$ and $c>0$ odd and coprime and complete \bigl{(}\begin{smallmatrix}a\\ c\end{smallmatrix}\bigr{)} to a matrix \bigl{(}\begin{smallmatrix}a&b\\ c&d\end{smallmatrix}\bigr{)}\in\text{SL}(2,\mathbb{Z}), then we find asymptotic expansions

[TABLE]

as $n\to\pm\infty$ with $n\in\mathbb{Z}$, and similar expansions with other coefficients $h^{\pm}_{i,\alpha,\beta}$ and $c^{\pm}_{i,\alpha,\beta}$ if $n\to\pm\infty$ with $n\in\mathbb{Z}+\beta$ for some fixed rational number $\beta$.

Similar statements hold for $x$ tending to rational numbers $\alpha$ with an even numerator or denominator, but now without the log term, the simplest case being

[TABLE]

where $A_{k}$ is the number of “up-down” permutations of $\{1,\dots,k\}$ (permutations $\pi$ where $\pi(i)-\pi(i-1)$ has sign $(-1)^{i}$ for all $i$), which is also equal to the coefficient of $x^{k}/k!$ in $\tan x+\sec x$. (Note that $A_{k-1}=(4^{k}-2^{k})|B_{k}|/k$ for $k$ even, so that the coefficients here are closely related to those in (39).) The difference between the two cases is due to the fact that the action of $\Gamma_{\vartheta}$ on $\mathbb{P}^{1}(\mathbb{Q})$ has two orbits (cusps), and comes from the relation that we will see in the final section between the quantum modularity properties of $C$ and $H$ and a specific Maass modular form on $\Gamma_{\vartheta}$ that is exponentially small at the cusp $\infty$ but has logarithmic growth at the cusp 1 .

5. Generalization to other odd Dirichlet characters

In this section we discuss the changes that are needed when the character defined in (1) is replaced by an arbitrary odd primitive Dirichlet character. Suppose that $D<0$ is the discriminant of an imaginary quadratic field $K$ (or equivalently, that $D$ is either square-free and congruent to 1 mod 4 or else equal to 4 times a square-free number not congruent to 1 mod 4), and let \chi=\chi_{D}=\bigl{(}\frac{D}{\cdot}\bigr{)} be the associated Dirichlet character, the case studied up to now corresponding to $K=\mathbb{Q}(i)$ and $D=-4$. The Dirichlet $L$-series $L(s,\chi)$ is defined as $\sum_{n>0}\chi(n)n^{-s}$ for $\sigma>1$ and by analytic continuation otherwise, and is equal to the quotient of the Dedekind zeta function of $K$ by the Riemann zeta function. Its value at $s=1$ is well-known to be non-zero and given by

[TABLE]

where $h^{\prime}(D)$ is $1/3$ or $1/2$ if $D=-3$ or $D=-4$ and otherwise is the class number of $K$.

We now define $S(x)$ by the same formula (11) as before, with $\chi=\chi_{D}$. It is still periodic (with period $|D|$) and even and hence bounded, but now has average value $h^{\prime}(D)$ rather than 1/2 and also no longer takes on only the values 0 and 1, as one sees in the two following pictures of the graphs of this function for $D=-3$ and $D=-7$.

$-1$[math]1$$2$$3$$4$$5$$1

$-1$[math]1$$2$$3$$4$$5$$6$$7$$8$$1$$2

Figure 5. Graphs of the functions $S(x)$ for $D=-3$ and $D=-7$

We define $H(x)=H_{D}(x)$ and $C(x)=C_{D}(x)$ for $x\in\mathbb{Q}$ by the same formulas (10) and (8) as before, but we define $h_{m,n}$ and $c_{m,n}$ by (7) and (5) with a factor $|D|/\pi$ instead of $4/\pi$ in order to get integral linear combinations of cotangents. The Fourier expansions (31) and (26) are then replaced by

[TABLE]

and the formula corresponding to (14) now reads

[TABLE]

by the same calculations as before. Finally we also define the matrices ${\sf C}_{N}$ and $\widehat{\sf C}_{N}$ exactly as we did in the special case $D=-4$ (except for replacing the lower right-hand entry in $\widehat{\sf C}_{N}$ in (3) by $|D|N/\pi$), and then from the scalar product calculations

[TABLE]

which are proved exactly as before, we deduce the same connection as in Theorems 1 and 3 between the unboundedness of the function $R(N)$ and the Riemann hypothesis for $L(s,\chi)$.

The only real difference with the case $D=-4$ is in the argument for quantum modularity. The function $H(x)$ on $\mathbb{Q}$ is still even and periodic of period $|D|$ (and also anti-periodic up to a constant with period $|D|/2$ if $D$ is even), it again has discontinuities at infinitely many rational points by an argument similar to the one given in 4.3, and by the analog of Proposition 8 the function $C:\mathbb{Q}\to\mathbb{R}$ again has a continuous extension to $\mathbb{R}$ and is therefore much better behaved analytically than $H(x)$, as illustrated by the following graphs of these functions for $D=-3$, which look qualitatively much like their $D=-4$ counterparts in Figures 2 and 3.

Figure 6. Graph of the function $H(x)$ for $D=-3$

Figure 7. Graph of the function $C(x)$ for $D=-3$

The difference is that this longer suffices to prove the quantum modularity of $H$ because the subgroup of $\text{SL}(2,\mathbb{Z})$ generated by the matrices $T^{|D|}$ (or $T^{|D|/2}$ if $D$ is even) and $S$ now no longer has finite index. Instead, we need a statement like the corollary to Theorem 2 with $\Gamma_{\vartheta}$ replaced by the congruence subgroup $\Gamma_{(D)}=\Gamma_{0}(D)\cup S\Gamma_{0}(D)$ of $\text{SL}(2,\mathbb{Z})$ (or by $\Gamma_{(D/2)}$ if $D$ is even, which is indeed $\Gamma_{\vartheta}$ if $D=-4$). This statement is given in the following theorem.

Theorem 4**.**

Let $D<0$ and $H(x)=\sum_{k>0}\chi_{D}(k)S_{D}(kx)/k$ be as above. Then the function

[TABLE]

extends continuously to $\mathbb{R}$ for all matrices \gamma=\bigl{(}\begin{smallmatrix}a&b\\ c&d\end{smallmatrix}\bigr{)}\in\Gamma_{(D)}, where $\varepsilon:\Gamma_{(D)}\to\{\pm 1\}$ is the homomorphism mapping $\Gamma(D)$ to $1$ and $S$ to $-1$.

This theorem is illustrated for $D=-3$ and a typical element of $\Gamma_{(-3)}$ in the following figure.

Figure 8. Graph of $C_{\gamma}(x)$ for $D=-3$, \gamma=\bigl{(}\begin{smallmatrix}4&-3\\ 3&-2\end{smallmatrix}\bigr{)}

We will deduce Theorem 4 as a consequence of the following proposition, which is a generalization to arbitrary elements $\gamma\in\Gamma_{(D)}$ of equation (32) for $D=-4$ and $\gamma=S$.

Proposition 9**.**

Let $\gamma=\begin{pmatrix}a&b\\ c&d\end{pmatrix}\in\Gamma_{(D)}$ with $c\neq 0$. Then

[TABLE]

To see that this implies the continuity, we note that we can rewrite the right-hand side by interchanging the summation and integration in the form (generalizing (25))

[TABLE]

for $x\gtrless-d/c$. The right-hand side of this formula is continuous because $S(t)$ is piecewise continuous and bounded and has support in $|t|\geq 1$, while the second factor of the integrand is piecewise continuous and bounded by a constant times $t^{-2}$, as one sees easily by partial summation. Another argument is that (44) is equivalent to the equality of distributions

[TABLE]

(generalizing  (27)) and we already know that the distribution $H$ is locally the derivative of a continuous function.

For the proof of (44) we will restrict to the case when the matrix $\gamma$ belongs to $\Gamma(D)$, which is sufficient for the proof of Theorem 4 because $\gamma\mapsto C_{\gamma}$ is a cocycle and we already know the continuity of $C_{\gamma}$ ($=C$) for the special case $\gamma=S$ representing the non-trivial coset of $\Gamma(D)$ in $\Gamma_{(D)}$. To avoid distracting case distinctions, we consider only the case $x+d/c<0$ and $\gamma x>0$. The argument is the same in principle in all other cases, but we are actually working with ordered tuples of points on the 1-manifold $\mathbb{P}^{1}(\mathbb{R})$, on which the group $\Gamma_{(D)}$ acts in an orientation-preserving way, and it is notationally simpler to fix the positions of the points occurring with respect to the point at infinity, so that we can work on $\mathbb{R}$ instead. We may also assume (by replacing $\gamma$ by $-\gamma$ if necessary, but then remembering that $\gamma$ may be congruent to $-\mathbb{1}_{2}$ rather than $\mathbb{1}_{2}$ modulo $D$) that $c>0$. Finally, to be able to work with absolutely convergent sums and to reorder the terms freely, we will replace $H(x)$ by the function $H_{s}(x)$ ($s>1$) introduced in Section 4.2 and $C_{\gamma}$ by the corresponding cocycle

[TABLE]

and only set $s=1$ at the end. The equation we have to prove then becomes

[TABLE]

We denote the last term (without the minus sign) by $A$. Since the variable $t$ in the integral is always positive (because we are considering the case $\gamma x>0$), we can replace $S(t)$ by its definition (11) and interchange the summation and integration to getf

[TABLE]

In the second sum we replace the vector $\begin{pmatrix}j\\ k\end{pmatrix}$ by the vector $-\gamma\begin{pmatrix}j\\ k\end{pmatrix}=\begin{pmatrix}-aj-bk\\ -cj-dk\end{pmatrix}$. This does not change the product $\chi(j)\chi(k)$ because $\gamma$ is congruent to plus or minus the identity modulo the period $|D|$ of $\chi$, changes the expression $cj-ak$ in the denominator to $k$, and changes the inequalities $k>0,\,j\geq k\gamma x$ to $kx\leq j<-kd/c$ (which imply $k>0$). Hence

[TABLE]

completing the proof of equation (48) and hence of the proposition and theorem.

6. Modular forms are everywhere

As one can already see from the examples in the original article [16] where this notion was introduced, quantum modular forms are sometimes related to actual modular forms. These modular forms may be either holomorphic or Maass forms, with the quantum modular form in the latter case being related to the “periods” of Maass forms in the sense developed in [13] and [7]. The quantum modular form $H(x)$ that we have been studying in this paper turns out to be of this latter type, with the associated Maass form being an Eisenstein series with eigenvalue 1/4 for the hyperbolic Laplace operator. In this final section we explain how this works, first for the case $D=-4$ studied in the first four sections of this paper. In that case the relevant modular group $\Gamma_{\vartheta}$ is generated by a translation and an inversion, so that for both the period theory of the Maass form $u$ and the quantum modularity of $H$ one needs only the functional equation of the associated $L$-series. Then at the end we indicate how the quantum modularity of $H$ for general $D$ follows from the full theory of periods of Maass forms.

We start with the case $D=-4$, so that $\chi$ is the character given by (1). The relevant modular form here is the Maass Eisenstein series

[TABLE]

where $\mathfrak{H}$ denotes the upper half-plane and $K_{0}$ the usual $K$-Bessel function of order 0. This is an eigenfunction with eigenvalue 1/4 with respect to the hyperbolic Laplace operator \Delta=-y^{2}\,\bigl{(}\frac{\partial^{2}}{\partial x^{2}}\,+\,\frac{\partial^{2}}{\partial y^{2}}\bigr{)} and is a modular function with character $\varepsilon$ for the group $\Gamma_{\vartheta}$ defined in Section 1, meaning that $u(\gamma z)=\varepsilon(\gamma)u(z)$ for all $\gamma\in\Gamma_{\vartheta}$ or, more explicitly, that

[TABLE]

To see this, we observe that $u(z)$ is proportional to E\bigl{(}\frac{z+1}{4},\frac{1}{2}\bigr{)}-E\bigl{(}\frac{z-1}{4},\frac{1}{2}\bigr{)}, where $E(z,s)$ is the usual non-holomorphic Eisenstein series of weight 0 and eigenvalue $s(1-s)$ for the Laplace operator with respect to the full modular group $\Gamma_{1}=\text{SL}(2,\mathbb{Z})$. (This follows easily from the well-known Fourier expansion of $E(z,\frac{1}{2})$ as a linear combination of the three functions $\sqrt{y}$, $\sqrt{y}\log y$, and $\sqrt{y}\sum_{n\neq 0}d(n)K_{0}(2\pi|n|y)e^{2\pi inx}$.) The transformation equations (50) then follow from the $\text{SL}(2,\mathbb{Z})$-invariance of $E(z,\frac{1}{2})$, the first one trivially since $E(z+1,s)=E(z,s)$ and the second by using the invariance of $E(z,s)$ under \bigl{(}\begin{smallmatrix}1&0\\ \mp 4&1\end{smallmatrix}\bigr{)} to get

[TABLE]

We now associate to $u(z)$ the periodic holomorphic function $f$ on $\mathbb{C}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,\mathbb{R}=\mathfrak{H}^{+}\cup\mathfrak{H}^{-}$ (where $\mathfrak{H}^{\pm}=\{z\in\mathbb{C}\mid\pm\Im(z)>0\}$) having the same Fourier coefficients as $u$, i.e.,

[TABLE]

Proposition 10**.**

The period function $\psi(z)$ defined by

[TABLE]

extends holomorphically from $\mathbb{C}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,\mathbb{R}$ to $\mathbb{C}^{\prime}=\mathbb{C}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,(-\infty,0]\,$.

Proof. We follow the proof of the corresponding result given in Chapter 1 of [13] for Maass cusp forms on the full modular group $\Gamma_{1}$. (See Theorem on p. 202 of [13], which also gives a converse statement characterizing cusp forms in terms of holomorphic functions $\psi$ on $\mathbb{C}^{\prime}$ satisfying a certain functional equation.) That proof required only the functional equation of the $L$-series associated to the Maass cusp form, which worked because the group $\Gamma_{1}$ is generated by the translation $T$ and the inversion $S$, and can be applied here because $\Gamma_{\vartheta}=\langle S,T^{2}\rangle$ has a similar structure. Our situation is also a little different because our function $u(z)$ is an Eisenstein series rather than a cusp form, but since the Fourier expansion (49) of $u(z)$ at infinity has no constant term this has no effect on the proof.

The argument in [13] was first to write the $L$-series of the Maass form $u(z)$ as the Mellin transform of the restriction of $u$ (or of its normal derivative in the case of an odd cusp form) to the imaginary axis, multiplied by a suitable gamma factor, and to deduce from this relationship and the $S$-invariance of $u$ a functional equation for the $L$-series. One then observed that the Mellin transforms of the restrictions to the positive or negative imaginary axis of both the associated periodic holomorphic function $f$ and the associated period function $\psi$ were also equal, up to different gamma factors, to the same $L$-series, and the functional equation of this $L$-series combined with Mellin inversion then led to a formula for $\psi$ that applied in all $\mathbb{C}^{\prime}$ rather than just on $\mathbb{C}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,\mathbb{R}$. Here the first step, which would require the normal derivative since the restriction of the function (49) to $i\mathbb{R}_{+}$ vanishes, can be skipped since the $L$-series $L(u;s)=\sum\chi(n)d(n)n^{-s}$ is simply the square of the Dirichlet $L$-series (1) and hence has a known functional equation. The Mellin transforms of the restrictions to $f$ and $\psi$ to the positive or negative imaginary axis are then given by

[TABLE]

[TABLE]

where in the last line we have used the functional equation of $L(s)$ and standard identities for gamma functions. By the Mellin inversion formula we deduce that

[TABLE]

for any $c\in(0,1)$, and analytic continuation from $i\mathbb{R}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,\{0\}$ to $\mathbb{C}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,\mathbb{R}$ then gives the formula

[TABLE]

for all $z\in\mathbb{C}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,\mathbb{R}\,$. But the factor in square brackets is bounded by a power of $s$ times $e^{-\pi|s|}$ for $|s|\to\infty$ on the vertical line $c+i\mathbb{R}$, so the integral on the right-hand side of this equation is absolutely convergent for all $z\in\mathbb{C}$ with $|\arg(z)|<\pi$, i.e., for all $z\in\mathbb{C}^{\prime}\,$. $\quad\square$

We now show how to connect the functions $f(z)$ and $\psi(z)$ to the functions $H(x)$ and $C(x)$, respectively, and to deduce the continuity of $C$ – and hence the quantum modularity of $H$ – from Proposition 10 for the period function $\psi$. The first step is easy, since from the Fourier expansions (51) and (26) and standard integral formulae we find the relation

[TABLE]

between the periodic function $f(z)$ and the periodic distribution $H(x)$. (Here and in what follows the symbol $\,\doteq\,$ denotes equality up to easily computed scalar factors whose values are irrelevant for the argument.) Replacing $t$ by $1/t$ in the integral, we obtain

[TABLE]

where we have used (27) to get the first equality and integration by parts for the second. But now replacing $t$ by $1/t$ and using the functional equation $C(t)=|t|C(1/t)$ we find

[TABLE]

with the same proportionality constant as in (53), and adding these two equations gives the relation

[TABLE]

between the period function of the Maass wave form $u$ and the function $C(x)$. This establishes the desired connection between the continuity of $C(x)=H(x)+|x|H(1/x)$, which expresses the quantum modularity of $H$, and Proposition 10, which expresses the modularity of $u$: in one direction, if we know that $C(x)$ is continuous (and bounded by $\,\min(1,|x|)$), then (55) immediately gives the analytic continuation of $\psi(z)$ from $\mathbb{C}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,\mathbb{R}$ to $\mathbb{C}^{\prime}$; and conversely, the inversion formula for the Stieltjes transform given in [14] lets us invert (55) to get

[TABLE]

where ${\mathcal{C}}$ is any contour with endpoints at $z=-1$ which encloses the origin, so that the continuity of $C$ is a direct consequence of the holomorphy of $\psi$ in $\mathbb{C}^{\prime}$.

We observe that the entire discussion given here applies in an essentially unchanged form to the more general functions $C_{s}$ and $H_{s}$ discussed in §4.2, with the Maass form (49) replaced by the form E\bigl{(}\frac{z+1}{4},\frac{s}{2}\bigr{)}-E\bigl{(}\frac{z-1}{4},\frac{s}{2}\bigr{)} with spectral parameter $s/2$ instead of $1/2$.

Finally, we consider the case when the character $\chi$ in (1) is replaced by an arbitrary primitive odd Dirichlet character $\chi_{D}$. Then all of the above calculations still go through: the function $u$ is defined by equation (49) with $\chi$ replaced by $\chi_{D}$ and $\pi/2$ by $2\pi/|D|$, which is again a Maass form (with spectral parameter $\frac{1}{2}$ and character $\varepsilon$ on the group $\Gamma_{(D)}$) because it is proportional to \sum_{r\!\!\pmod{D}}\chi(r)E\bigl{(}\frac{z+r}{|D|},\frac{1}{2}\bigr{)}; the associated periodic function $f$ and period function $\psi$ are defined by (51) (with the new $\chi$ and with $q^{n/4}$ replaced by $q^{n/|D|}$) and (52) (with no change at all); and Proposition 10 remains true with the same proof. The difference, however, is that this proposition is no longer equivalent to the modularity of $u$, but only to its invariance (up to sign) under the transformations $S$ and $T^{D}$ (or $T^{D/2}$ if $D$ is even), which in general generate a subgroup of $\Gamma_{(D)}$ of infinite order, as already discussed in Section 5. To get the full modularity (of this $u$ or any other potential Maass form on $\Gamma_{(D)}$), we need to generalize Proposition 10 to the statement that for any matrix \gamma=\bigl{(}\begin{smallmatrix}a&b\\ c&d\end{smallmatrix}\bigr{)}\in\Gamma_{(D)} the function

[TABLE]

extends holomorphically from $\mathbb{C}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,\mathbb{R}$ to $\mathbb{C}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,(-\infty,-d/c]\,$ if $c>0$ or $\mathbb{C}\,\hbox{ \leavevmode\hbox to3.6pt{\vbox to6.6pt{\pgfpicture\makeatletter\hbox{\hskip 0.3pt\lower-0.3pt\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }\definecolor{pgfstrokecolor}{rgb}{0,0,0}\pgfsys@color@rgb@stroke{0}{0}{0}\pgfsys@invoke{ }\pgfsys@color@rgb@fill{0}{0}{0}\pgfsys@invoke{ }\pgfsys@setlinewidth{0.4pt}\pgfsys@invoke{ }\nullfont\hbox to0.0pt{\pgfsys@beginscope\pgfsys@invoke{ }{{{}{}}{{}}{} {{}{}}{}\pgfsys@beginscope\pgfsys@invoke{ }\pgfsys@setlinewidth{0.6pt}\pgfsys@invoke{ }\pgfsys@roundcap\pgfsys@invoke{ }{}\pgfsys@moveto{3.0pt}{0.0pt}\pgfsys@lineto{0.0pt}{6.0pt}\pgfsys@stroke\pgfsys@invoke{ } \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope} \pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope{}{}{}\hss}\pgfsys@discardpath\pgfsys@invoke{\lxSVG@closescope }\pgfsys@endscope\hss}}\lxSVG@closescope\endpgfpicture}}}\,[-d/c,\infty)$ if $c<0$. This statement can be proved in several ways and can be linked by a discussion similar to the one above to the continuity property of $C_{\gamma}(x)$ stated in Theorem 4. But the simplest approach is to relate the functions $C_{\gamma}(x)$ directly to the invariant distribution associated to $u$ in the sense developed in [12], [13] and [7]. Specifically, this theory says that an eigenfunction $u$ of $\Delta$ with spectral parameter $s$ and Fourier expansion $u(x+iy)=\sqrt{y}\,\sum_{n\neq 0}A_{n}\,K_{s-1/2}(\lambda|n|y)e^{i\lambda nx}$ is invariant (possibly with character) under the action of a Fuchsian group $\Gamma$ if and only if the associated distribution $U(x)=\sum_{n\neq 0}|n|^{s-1/2}\,A_{n}\,e^{i\lambda nx}$ on $\mathbb{P}^{1}(\mathbb{R})$ is invariant (with the same character) with respect to the group action $(U|\gamma)(x)=|cx+d|^{-2s}U(\gamma x)$ for \gamma=\bigl{(}\begin{smallmatrix}a&b\\ c&d\end{smallmatrix}\bigr{)}\in\Gamma. (Here the word “distribution” must be interpreted correctly, namely, as a functional on the space of test functions $\phi(x)$ that are smooth on $\mathbb{R}$ and for which $|x|^{-2s}\phi(1/x)$ is smooth near $x=0$.) In our case the distribution $U$ associated to $u$ is given by $U(x)=\sum_{n=1}^{\infty}\chi(n)d(n)\sin(2\pi nx/|D|)$ and is related to the distribution (41) by $U(x)\doteq H^{\prime}(x)$. But that means that if we differentiate the definition (43) of $C_{\gamma}(x)$ to get

[TABLE]

then the first two terms on the right cancel and we recover equation (46).