1 Introduction
1.1 Definitions
Let
be an orthonormal basis of Hecke–Maass cusp forms on the modular surface
$ \mathrm{SL}_2 (\mathbf{Z}) \backslash \mathbf{H} $
. Assume that each
is either even or odd in the sense that
$ f (-\overline {z}) = (-1)^{\delta _f} f (z)$
for
$\delta _f = 0$
or
$ 1$
. Let
or
be the subset of even or odd forms in
respectively. For
, let
$\lambda _{f} = s_{f} (1-s_{f})$
be its Laplace eigenvalue, with
$s_f = 1/2+ i t_f$
(
$t_f> 0$
). The Fourier expansion of
$f (z)$
reads
where, as usual,
$K_{\nu } (x)$
is the K-Bessel function and
$e (x) = \exp (2\pi i x)$
. Let
$\lambda _f (n)$
(
$n \geqslant 1$
) be the Hecke eigenvalues of
$f (z)$
. It is well known that
$\rho _f (\pm n) = \rho _f (\pm 1) \lambda _f (n) $
, while
$\rho _f (-1) = (-1)^{\delta _f} \rho _f (1)$
. Define the harmonic weight
The Hecke–Maass L-function of
$f (z)$
is defined by
for
$\mathrm {Re} (s)> 1$
, and by analytic continuation on the whole complex plane.
1.2 Effective non-vanishing results
As the motivation of this article, the author [Reference QiQi] recently established an effective refinement (recorded below) of the non-vanishing result of Luo [Reference LuoLuo3] for the Rankin–Selberg special L-values
$L (s_f, Q \otimes f)$
. The non-vanishing condition for
$L (s_f, Q \otimes f)$
arises in the Phillips–Sarnak deformation theory of Maass cusp forms [Reference Phillips and SarnakPS] and was studied by Deshouillers, Iwaniec, and Luo [Reference Deshouillers and IwaniecDI, Reference LuoLuo1, Reference LuoLuo3, Reference LuoLuo4].
Theorem [Reference QiQi]
Let
$Q (z)$
be a fixed holomorphic cusp newform of even weight
$2k$
and square-free level q. Let
$3/4 < \mu < 1$
. We have


where
$\gamma (Q )$
is the Euler product:
with
$a (p)$
the p-th Hecke eigenvalue of
$Q (z)$
.
It is natural to consider the more fundamental family of Hecke–Maass L-functions
$L (s, f)$
and to see how large a non-vanishing proportion at
$s = s_f$
one might attain. It is also interesting to compare our non-vanishing results for
$L (s_f, f)$
with those for central L-values
$L(1/2, f)$
, with or without the Riemann hypothesis.
To be explicit, our main theorems for the non-vanishing statistic of the special L-values
$L (s_f, f)$
are as follows.
On the long interval
$t_f \leqslant T$
, Theorem 1 manifests that at least
$2 / \pi ^2 \approx 20 \%$
of the special L-values
$L (s_f, f) $
do not vanish for
, while it is proven in [Reference Balkanova, Huang and SödergrenBHS, Reference Liu and QiLQ2] that 25% of the central L-values
$L(1/2, f)$
do not vanish for
, although, in the latter case,
$L (1/2, f) = 0$
trivially on the whole odd sub-basis
as the root number
$\epsilon _f = - 1$
.
Theorem 2 Let
$ 0 < \mu < 1$
. We have

This is an effective refinement of the result of Xu [Reference XuXu] in which he proved that (on the even sub-basis
) the ratio in (1.5) is bounded below by an unexplicit constant as
$T \rightarrow \infty $
.
Theorem 3 Assume the Riemann hypothesis for every
$L (s, f)$
with
and for all Dirichlet L-functions (including the Riemann zeta function). Then,

Thus, on the Riemann hypothesis, for
the non-vanishing percentage 20% for
$ L (s_f, f) $
is raised drastically to 50%, but 25% for
$L(1/2, f)$
is raised only to 28% (see Corollary A.1).
Theorem 4 Assume the Riemann hypothesis as in Theorem 3. Then,

These conditional results will be deduced from the extended density theorem for
$L (s, f)$
at
$s = s_f$
.

Figure 1. Long description
A line graph with a horizontal axis labeled mu and a vertical axis with fractional values. The horizontal axis has tick marks at 0, 1 over 3, 1 over 2, and 1. The vertical axis has tick marks at 2 over pi-squared, 1 over 3, and 1 over 2.
Two lines are plotted:
* A black line starts at the origin (0,0) and follows a concave curve upward until it reaches the point (1 over 2, 2 over pi-squared). From that point, it continues as a horizontal flat line to the right until mu equals 1.
* A grey line starts on the horizontal axis at (1 over 3, 0). It rises steeply in a linear fashion, intersecting the black line before reaching a vertex at (1 over 2, 1 over 3). From there, it continues to rise with a shallower linear slope until it terminates at the point (1, 1 over 2).
Horizontal and vertical dashed grid lines connect the key vertices of both plots to their respective values on the axes.
1.3 Twisted spectral moments
Our proof of Theorems 1 and 2 will rely on the asymptotic formulas for the twisted first and second moments of
$ L (s_f , f ) $
.
Let
$T, \varPi $
be large parameters such that
$T^{\varepsilon } \leqslant \varPi \leqslant T^{1-\varepsilon }$
. Define the smoothly weighted twisted moments on the (even or odd) cuspidal spectrum:


Theorem 5 Let the notation be as above. We have
where
$\delta (m ,1)$
is the Kronecker
$\delta $
that detects
$m = 1$
.
Theorem 6 Let the notation be as above. Assume that
$m_1, m_2 < \varPi ^{1-\varepsilon }$
. Set
Then,

where
with
$\gamma $
the Euler constant, and
Corollary 7 Let
$T^{\varepsilon } \leqslant \varPi \leqslant T^{1-\varepsilon }$
. We have

Corollary 8 We have
where the superscript
$\delta $
indicates summation over
.
1.4 Density theorems
Our proof of Theorems 3 and 4 will rely on the extended density theorem under the Riemann hypothesis for Dirichlet L-functions.
Let us denote the non-trivial zeros of
$L (s, f)$
by
if the generalized Riemann hypothesis for
$L (s, f)$
is assumed, then
$\gamma _f$
are all real.
Definition 1.1 Define the
$1$
-level density of
$L (s, f)$
at
$s_f = 1/2+ it_f$
to be
where T is a large scaling parameter and
$\unicode{x3c6} : \mathbf {R} \rightarrow \mathbf {C}$
is an even Schwartz function such that its Fourier transform
has compact support.
For
$T, \varPi $
as above, define the smoothly averaged
$1$
-level density

Theorem 9 Let the notation be as above. Assume the Riemann hypothesis (for the Riemann zeta function). Then, for
$\hat {\unicode{x3c6} }$
of support in
$(-1, 1)$
, we have
for any fixed
$0 < \mu < 1$
.
Recall that
$W (\mathrm {U}) (x) \equiv 1$
according to Katz and Sarnak [Reference Katz and SarnakKS].
Theorem 9 manifests that the symmetry types of
$ L (s_f , f) $
on
and
are both unitary, while it is proven in [Reference Liu and QiLQ1, Appendix] that the symmetry types of
$ L (1/2 , f) $
on
and
are orthogonal of different parity:
$\mathrm {SO} (\mathrm {even})$
and
$\mathrm {SO} (\mathrm {odd})$
, respectively. Note that for
$\hat {\unicode{x3c6} }$
of support in
$(-1, 1),$
one is unable to differentiate
$\mathrm {SO} (\mathrm {even})$
and
$\mathrm {SO} (\mathrm {odd})$
.
However, in order to obtain a non-trivial non-vanishing result for
$ L (s_f , f) $
(as in Theorems 3 and 4), we must let the support of
$ \hat {\unicode{x3c6} } $
go beyond
$[-1, 1]$
as in the following extended density theorem.
Theorem 10 Define
$v (\mu ) = \min \{ 3\mu , 1 + \mu \}$
. Then, for
$1/3 < \mu < 1$
, we have
for
$\hat {\unicode{x3c6} }$
of support in
$(-v (\mu ) , v(\mu ) )$
if we assume, in addition, the Riemann hypothesis for Dirichlet L-functions.
Remark 1.1 We have the smoothly weighted harmonic Weyl law (see (2.13)):

Notation
By
or
$X = O (Y),$
we mean that
$|X| \leqslant c Y$
for some constant
${c> 0}$
, and by
$X \asymp Y,$
we mean that
and
. We write
or
$ {X = O_{\alpha , \beta , \ldots } (Y)} $
if the implied constant c depends on
$\alpha $
,
$\beta $
,
$\ldots $
.
By “negligibly small,” we mean
$ O_A ( T^{-A} )$
for arbitrarily large but fixed
$A $
.
The notation
$x \sim X$
stands for
$ X \leqslant x \leqslant 2 X $
.
Throughout the article,
$\varepsilon $
is arbitrarily small and its value may differ from one occurrence to another.
2 Kuznetsov trace formulas
Let
$h (t)$
be an even function satisfying the conditions:
-
(i) $h (t)$
is holomorphic in
$|\mathrm {Im}(t)|\leqslant {1}/{2}+\varepsilon $
; -
(ii)
in the above strip.
For
$m, n \geqslant 1$
, define the cuspidal term
and the Eisenstein term
where
$\omega _f$
and
$\lambda _{f} (n)$
are defined before in Section 1.1, and
On the other hand, define the Plancherel integral
and the Kloosterman–Bessel sum
where
$S (m, \pm n; c)$
is the Kloosterman sum and
$H^{\pm } (x)$
are the Bessel integrals
By [Reference WatsonWat, Formula 3.7 (6)],
where
$I_{\nu } (z)$
is the I-Bessel function. Hence, we may rewrite (2.7) in a form which is similar to (2.6):
The Kuznetsov trace formula for the even or odd cusp forms reads (see [Reference Conrey and IwaniecCI, Section 3]):
where
$\delta (m, n)$
is the Kronecker
$\delta $
-symbol.
2.1 The Kloosterman sum
Recall that
where the
$\star $
indicates the condition
$(\alpha , c) = 1$
and
$ \overline {\alpha }$
is given by
$\alpha \overline {\alpha } \equiv 1 (\mathrm {mod} \, c)$
. We have the Weil bound:
where, as usual,
$\tau (c) = \tau _0 (c) $
is the number of divisors of c.
2.2 The Hecke relation
We have
2.3 The Weyl law
Finally, we record here the harmonic-weighted Weyl law for
:

By the method in [Reference LiLi2], this follows from the Kuznetsov trace formula for
.
3 Hecke–Maass L-functions
For
, define
for
$\mathrm {Re} (s)> 1$
. It is known that
$L (s, f)$
has analytic continuation to the whole complex plane and satisfies the functional equation
with
where
Note that, in parallel to (3.1), if we replace
$\lambda _f (n)$
by
$ \tau _{ it} (n) $
as defined in (2.3), then
Moreover, we have the Euler factorization
for
$\mathrm {Re} (s)> 1$
, and the Kim–Sarnak bound [Reference KimKim]
4 Analysis for the Bessel integrals
Let
$T, \varPi $
be large parameters with
$ T^{\varepsilon } \leqslant \varPi \leqslant T^{1-\varepsilon } $
. In this section, we consider the Bessel transforms of test functions in the form
for
$\unicode{x3b3} (t) $
in the function space
$ \mathscr {Q}_{\kappa } $
or
$ \mathscr {T}_{\tau } $
defined below in Definitions 4.1 or 4.2, and
In the case that
$\unicode{x3b3} (t)$
is even, we also need to consider
Let
$H^{\pm } (x, y)$
and
$H^{\pm }_2 (x, y)$
denote the Bessel transforms of
$h (t; y)$
and
$ h_2 (t; y) $
as in (2.6) and (2.7) or (2.8). Note that
so the results below for
$H^{\pm } (x, y) $
will be valid for
$H^{\pm }_2 (x, y) $
as well.
Definition 4.1 For fixed real
$\kappa $
, let
$ \mathscr {Q }_{\kappa } $
denote the space of holomorphic functions
$\unicode{x3b3} (t) $
on the strip
$ |\mathrm {Im} (t)| \leqslant 1/2 + \varepsilon $
such that
Definition 4.2 For fixed real
$\tau $
, let
$ \mathscr {T }_{\tau } \subset C^{\infty } (\mathbf {R})$
denote the space of smooth functions on
$ \mathbf {R} $
such that

for any
$n \geqslant 0$
.
The Bessel integral
$H^{+} (x, y)$
has been studied thoroughly in [Reference QiQi, Section 4]. The results therein may be extended to
$ H^{-} (x, y) $
by literally the same proofs.
Lemma 4.1 Let
$\unicode{x3b3} (t) $
be in the space
$ \mathscr {Q}_{\kappa }$
. Set
$u = xy +x/y$
. Then, for
, we have
$H^{\pm } (x, y) = O_{\kappa } ( \varPi T^{ \kappa } u )$
.
Proof This lemma follows from shifting the integral contour in (2.6) and (2.8) up to
$\mathrm {Re} (it) = 1/2 + \varepsilon $
. For more details, see the proof of Lemma 4.1 in [Reference QiQi]. Note that for
$\mathrm {Re}(\nu )> -1/2$
, we have the Poisson integral representation (see [Reference WatsonWat, Formula 3.71 (9)]):
For the rest of this section, let us assume that
$ \unicode{x3b3} (t) $
lies in the space
$ \mathscr {T}_{\tau }$
.
Lemma 4.2 Let
$\unicode{x3b3} (t) $
be in the space
$ \mathscr {T}_{\tau }$
. Let
$A \geqslant 0$
. Then, for
, we have
where
and
$g (r)$
is a Schwartz function (namely,
$ g^{(n)} (r) = O ( (|r|+1)^{-m})$
for any
$m, n \geqslant 0$
).
Later in our analysis, we shall always split
so that we may extract the exponential factor
$ \exp (\pm i v) $
or
$\exp (\pm i w)$
, and we shall usually drop
$g ( \varPi r) \exp ( 2i T r)$
as it does not involve
$x, y$
or
$v, w$
and is trivially
$O (1)$
.
Proof This lemma follows essentially from inserting the Mehler–Sonine integrals (see [Reference WatsonWat, Formulas 6.21 (12) and 6.22 (13)])
into the Bessel integrals as in (2.6) and (2.7). For more details, see the proof of Lemma 4.3 in [Reference QiQi].
Moreover, Lemma 4.5 in [Reference QiQi] may be extended easily by applying stationary phase analysis to the integral in (4.7).
Lemma 4.3 Let
$T^{\varepsilon } \leqslant \varPi \leqslant T^{1-\varepsilon }$
. Then,
$H^{\pm } (x, y) = O (T^{-A} )$
if
.
By Lemmas 4.1–4.3, we conclude that if
$\unicode{x3b3} (t)$
is in
$ \mathscr {Q}_{\tau } \cap \mathscr {T}_{\tau } $
, then the Bessel integrals
$H^{\pm } (x, y)$
are negligibly small for
.
Corollary 4.1 Let
$\unicode{x3b3} (t)$
be in the space
$ \mathscr {Q}_{\tau } \cap \mathscr {T}_{\tau } $
. Set
$u = xy +x/y$
as before. Then, for
, we have
for any
$A \geqslant 0$
.
Remark 4.1 Later, in practice,
$ \unicode{x3b3} (t)$
will be those
as defined in (6.8) and (6.9) for
$|\mathrm {Im}(v)| \leqslant \log T$
or simply
$\unicode{x3b3} (t) \equiv 1$
in Section 13. Note that Definition 4.1 requires
$\unicode{x3b3} (t) $
to be holomorphic on the strip
$ |\mathrm {Im} (t)| \leqslant 1/2 + \varepsilon $
. For the odd case
$\delta = 1$
, it follows from Lemma 6.4 that
$ \unicode{x3b3} _{\sigma}^{1} (v, t) $
lies in
$ \mathscr {Q}_{ \unicode{x3c3} \varepsilon /2} \cap \mathscr {T}_{ \unicode{x3c3} \varepsilon /2} $
if we let
$\mathrm {Re} (v) = \varepsilon $
in (6.14), so we may just choose
$\kappa = \tau = \unicode{x3c3} \varepsilon /2$
. For the even case
$\delta = 0$
, however, it is slightly more subtle: In view of Lemma 6.4, we may apply Corollary 4.1 with
$\mathrm {Re} (v) = 1$
and
$ \tau = \unicode{x3c3} /2 $
say for
, and use the integral representation (4.7) in Lemma 4.2 with
$ \mathrm {Re} (v) = \varepsilon $
and
$ \tau = \unicode{x3c3} \varepsilon /2$
for
.
Remark 4.2 In [Reference QiQi, Section 4], we also studied the Bessel transform of a certain
$h_{^ -}^{} (t; y)$
that contains the factor
$\unicode{x3b5} ^{\delta } (t)$
defined similarly to (6.10), but here we shall not deal with this case (see the discussion at the end of Section 6.3).
Remark 4.3 In the study of related Maass central L-values, the
$H^{+}$
and
$H^{-}$
Bessel integrals have phases
$x \cosh r$
and
$x \sinh r$
, and hence behave quite differently in [Reference LiLi1] or [Reference YoungYou, Section 7]; but here, as we have just seen, the phases are
$f^{_{\phantom {i}}}_{^ \pm } (r; v, w) = v \exp (r) \pm w \exp (- r)$
and may be analyzed more or less in the same way.
5 A useful unsmoothing lemma
As in (4.2), for large parameters
$T, \varPi $
with
$T^{\varepsilon } \leqslant \varPi \leqslant T^{1-\varepsilon }$
, define
In this section, we introduce a useful unsmoothing process as in [Reference Ivić and JutilaIJ, Section 3] by an average of the weight
$\unicode{x3D5} _{T, \varPi } $
in the T-parameter.
Definition 5.1 Let
$ H \leqslant T / 3 $
and
$T^{\varepsilon } \leqslant \varPi \leqslant H^{1-\varepsilon } $
. Define
The next lemma is essentially Lemma 5.3 in [Reference Liu and QiLQ2], adapted from the arguments in [Reference Ivić and JutilaIJ, Section 3].
Lemma 5.1 Let
$\lambda $
be a real constant. Suppose that
$a_f $
is a sequence such that
for any
$\varPi $
with
$ T^{\varepsilon } \leqslant \varPi \leqslant T^{1-\varepsilon } $
. Then, for
$ \varPi ^{1+\varepsilon } \leqslant H \leqslant T / 3, $
we have
Later, we shall choose
$a_f$
to be
$\delta (L (s_f, f ) \neq 0)$
(the Kronecker
$\delta $
that detects
$L (s_f, f ) \neq 0$
),
$L (s_f, f ) $
, or
$|L (s_f, f )|^2$
(see Sections 10.4, 11, and 16). Note that in all these cases, we may let
$\lambda = 1$
or
$ 1+\varepsilon $
in (5.3) in light of the (harmonic-weighted) Weyl law or the mean Lindelöf bounds.
Part I
Twisted and mollified moments
6 Preliminaries
6.1 Poisson summation formula
The following Poisson summation formula is essentially [Reference Iwaniec and KowalskiIK, Formula (4.25)].
Lemma 6.1 Let
$\alpha , c $
be integers with
$c \geqslant 1$
. Then, for
$F (x) \in C_c^{\infty } (-\infty , \infty ),$
we have
where
$ \hat {F} (y)$
is the Fourier transform of
$F (x)$
defined by
6.2 Mean-value bounds for the Riemann
$\zeta $
function
The best bound for the mean square of
$ \zeta (1/2+it) $
on short intervals is due to Bourgain and Watt [Reference Bourgain and WattBW, Theorem 3]:
From this and the Cauchy–Schwarz inequality, we deduce easily the next lemma.
Lemma 6.2 Let
$1 \leqslant H \leqslant T/3$
. We have
6.3 Approximate functional equations
Let
$X> 1$
be a large parameter. By applying [Reference Iwaniec and KowalskiIK, Theorem 5.3] to
$L ( s , f )$
at
$s = s_f$
and to
$ L (s+it_f , f ) L (s-it_f, f ) $
at
$s = 1/2$
, we deduce, from (3.1) and (3.2), the approximate functional equations:
with
where
Moreover, in parallel, we have
up to certain residual error terms of bounds
$ O \big ( (|t| + 1)^{1/4+\varepsilon } X^{1/2} \big ) $
and
$O \big (\exp (- t^2/2) \big ) $
, respectively.
The next two lemmas are standard consequences of the Stirling formula. Compare [Reference Iwaniec and KowalskiIK, Proposition 5.4], [Reference BlomerBlo, Lemma 1], and [Reference QiQi, Lemmas 3.6 and 3.7].
Lemma 6.3 Let
and
$\theta> 0$
. Then, for
$\unicode{x3c3} = 1, 2$
, we have

for any
$A \geqslant 0$
, and
Lemma 6.4 Let
$\mathrm {Re}(v)> - 1/2$
be fixed. Then, for
$|\mathrm {Re} (it)| < \mathrm {Re}(v/2) + 1/4 + \delta /2 $
, we have
and, for t real, we have

for any
$n \geqslant 0$
.
The reader should compare Lemma 6.4 with Definitions 4.1 and 4.2.
Later, the analysis will be restricted to the range
$|t_f-T|, \, |t - T| \leqslant \varPi \log T$
, so it follows from (6.13) that the second sum in (6.5) or (6.11) (after
$\unicode{x3b5} ^{\delta _f} (t_f)$
or
$\unicode{x3b5} ^{0} (t) $
) will be negligibly small if we choose
while the first sum may be truncated effectively at
$ n = T^{1+\varepsilon } $
. Moreover, the residual error term for (6.11) has bound
$O (T^{1/2+\varepsilon })$
. Note that the error term in (6.14) is negligibly small if we choose
6.4 Bi-variable notation
Subsequently, particularly in our study of the second twisted moment in Section 9, we shall use bold
$\boldsymbol {m}$
,
$\boldsymbol {n}$
,
$\ldots $
for the (positive integral) pairs
$(m_1, m_2)$
,
$(n_1, n_2)$
,
$\ldots $
. Define the norm
$ \|\boldsymbol {n}\| = n_1 n_2 $
, the fraction
$\langle \boldsymbol {n} \rangle = n_1/ n_2$
, the dual
$\widetilde {\boldsymbol {n}} = (n_2, n_1)$
, the reduction
$\boldsymbol {n}^{\star } = \boldsymbol {n} / \mathrm {gcd}(\boldsymbol {n}) $
, the product
$ \boldsymbol {m n} = (m_1 n_1, m_2 n_2) $
, and the quotient
$ {\boldsymbol {m}/ \boldsymbol {n} = (m_1 / n_1, m_2 / n_2)} $
. Let
$\delta (\boldsymbol {n})$
be the Kronecker
$\delta $
that detects
$n_1 = n_2$
.
Define
For an arithmetic function
$a (n)$
, define
For instance, (1.9) now reads
Occasionally, we shall use
$\boldsymbol {x} $
for the positive real pair
$(x_1 , x_2)$
, and, similar to (6.18), for a smooth function
$\textit {w} (x)$
write
7 Setup
For large parameters
$T, \varPi $
with
$T^{\varepsilon } \leqslant \varPi \leqslant T^{1-\varepsilon }$
, let us define
and consider
Note that
or
respectively, differs from
or
as in (1.8) or (1.9) (see also (6.19)) merely by an exponentially small error term. In parallel, the Eisenstein contributions read
where we have used the t-evenness of
$\omega (t)$
,
$ \tau _{it} (m) $
, and
$ \mathrm {c}_{ t } (\boldsymbol {m}) $
.
Lemma 7.1 We have
8 The twisted first moment
8.1 Application of the Kuznetsov formula
By the approximate functional equations (6.5) and (6.11), with
we have

where the error terms are due to the discussions at the end of Section 6.3, and
It follows from the Kuznetsov trace formula in (2.9) that
with the diagonal
where
and the off-diagonal

where
8.2 Asymptotic for the diagonal term
In view of (8.1) and (8.4), we have
Keep in mind that
$\unicode{x3D5} (t)$
as in (7.1) is negligibly small unless
$|t - T | \leqslant \varPi \log T$
. Thus, up to an exponentially small error,
$\tanh (\pi t) $
is removable as
$\tanh (\pi t) = 1 + O (\exp (-\pi T))$
on this range.
Let us first evaluate
. By shifting the integral contour in (6.7) from
$\mathrm {Re} (v) = 3$
down to
$\mathrm {Re} (v) = \varepsilon -1/2 $
, we have
by the Stirling formula. Therefore,
as
$X = T^{1/2+\varepsilon }$
, whereas
Consequently,
Note that in the odd case of
the error term may be improved into
$O (\varPi / T^{1/2-\varepsilon })$
by shifting the integral contour further down to
$\mathrm {Re} (v) = \varepsilon -3/2 $
.
Next, we consider the case
$m> 1$
. By inserting the integral in (6.14), with
$\theta = \varepsilon $
and
$U = \log T$
, we arrive at
For
$t $
and
$ v$
on the integral domains, in view of (6.8), it follows from the Stirling formula that
Therefore,
For
$m> 1$
, the inner t-integral is Fourier of phase
$ - 2 t \log m$
; so, by repeated partial integration, we infer that the integral is negligibly small. Consequently, we have
$ H_1^{\delta } (m) = O (\varPi T^{\varepsilon }) $
, and hence

To summarize, we conclude with the following lemma for
.
Lemma 8.1 We have

8.3 Estimation of the off-diagonal sums
Let us consider the off-diagonal sums
given by (8.1) and (8.5)–(8.7).
By (6.13) in Lemma 6.3 (see also the discussion at the end of Section 6.3) along with a smooth dyadic partition, at the cost of a negligible error, we may confine the n-sum in (8.5) on the range
$ n \sim N $
for dyadic
(namely, N is of the form
$2^{j/2}$
). Next, we insert the expression of
$ V_{1}^{\delta } ( n / X; t ) $
as in (6.14) in Lemma 6.3 into (8.1), (8.6), and (8.7). After this, there arise the Bessel transforms
$H_1^{\delta \, \pm } (x, y; v)$
of the function
as in Section 4, for
Since
$\delta $
and v are not essential in our analysis, they will usually be omitted in the sequel. By Corollary 4.1 and Remark 4.1, along with the Weil bound (2.11), we may further restrict the c-sum in (8.5) to the range
. It follows that up to a negligible error, the off-diagonal sum
has bound

where, if we omit
$\delta $
and v from the notation,

with
for a certain fixed smooth weight
$\textit {v} \in C_c^{\infty }[1,2]$
. Note that v will not play a role since we have uniformly
Next, in (8.10), we open the Kloosterman sum
$S (m, \pm n;c)$
(see (2.10)) and insert the formula of the Bessel integrals
$H_1^{ \pm } (x, y)$
as in Lemma 4.2. For the latter, we have
then, we split
$\cos ( f^{_{\phantom {i}}}_{^{\pm } } (r; v, w) )$
as in (4.9) and extract the exponential factor
$e (mn/ c)$
or
$ e (- mn/ c) $
. It follows that

where

with
$\textit {w}^{_{\phantom {i}}}_{^ +} (x) = \textit {w} (x) $
, and
$\textit {w}^{_{\phantom {i}}}_{^ -} (x) = \overline {\textit {w} (x) }$
. On applying to the n-sum the Poisson summation formula in Lemma 6.1, we have

Note that the Fourier transform
$ \hat {\textit {w}}^{_{\phantom {i}}}_{^ \pm } ( N (n-m \unicode{x3c1} (r))/ c ) $
is negligibly small unless
As
we may restrict the n-sum to the range
. Consequently, up to a negligible error,

Now the r-integral may be confined further to the interval defined by (8.13) of length
$ O (cT^{\varepsilon } / m N) $
, and hence

On inserting (8.14) into (8.11), we have
Note that the implicit constant in (8.15) depends on
$\varepsilon $
but is uniform for all
$ v \in [\varepsilon - i \log T, \varepsilon + i \log T]$
. Consequently, from (8.9) and (8.15), we deduce the following bound for
.
Lemma 8.2 We have
9 The twisted second moment
As the bi-variable notation will be used frequently in what follows, the reader is advised to read Section 6.4 before proceeding with this section.
9.1 Application of the Kuznetsov formula
First, in the bi-variable notation, by the approximate functional equations (6.6) and (6.12) (here, we need to expand
$2 \mathrm {c}_t ( \boldsymbol {m} ) \mathrm {c}_t ( \boldsymbol {n}) = \mathrm {c}_t ( \boldsymbol {m} \boldsymbol {n}) + \mathrm {c}_t ( \boldsymbol {m} / \boldsymbol {n}) $
(see (6.17)) and use the symmetry in the
$\boldsymbol {n}$
-sums), we obtain

where
$ \boldsymbol {n} = (n_1, n_2)$
and
$ \| \boldsymbol {n}\| = n_1 n_2 $
as in Section 6.4, and
Note that
$ {V}_{2}^{\delta } (\boldsymbol {m}, \boldsymbol {n}; t)$
is even in t and real-valued. By the Hecke relation (2.12),

where, according to the bi-variable notation in Section 6.4, the first summations are over
$ \boldsymbol {d} = (d_1, d_2) $
and
$ \boldsymbol {h} = (h_1, h_2)$
such that
$ (d_1 h_1, d_2 h_2) = (m_1, m_2)$
. Now an application of the Kuznetsov trace formula in (2.9) yields
with the diagonal

where
$ \delta (\boldsymbol {hn}) $
is the Kronecker
$\delta $
that detects
$ h_1 n_1 = h_2 n_2 $
, and since
$ \mathrm {c}_{t} ( \boldsymbol {d^2 h n} ) = \mathrm {c}_{2t} ( \boldsymbol {d} ) $
if
$\delta (\boldsymbol {hn}) = 1$
,
and the off-diagonal

where
$\boldsymbol {1}_{\pm } = (1, \pm 1)$
; hence,
$S ( \boldsymbol {1}_{\pm } \boldsymbol {h n}; c) = S (h_1 n_1, \pm h_2 n_2; c)$
, and
9.2 Asymptotic for the diagonal sum
Set
then,
$ \delta ({\boldsymbol {h n} }) = 1 $
if and only if there is n such that
$ \boldsymbol {n} = n \widetilde {\boldsymbol {h^{\star } }}. $
Recall from Section 6.4 that
$ \widetilde {\boldsymbol {h^{\star }}} $
is dual to
$\boldsymbol {h^{\star }}$
. Therefore, (9.2) is turned into

In view of (6.13) and (9.3), the summation in (9.7) may be restricted effectively to
$ n \sqrt {\| \boldsymbol {m} \|} / h \leqslant T^{1/2+\varepsilon }$
. It follows from (6.13) and (7.1) that
$ H_2^{\delta } (\boldsymbol {d}, \boldsymbol {n}) = O (\varPi T) $
; so trivial estimation yields

where
$\boldsymbol {m^{\star }} = \boldsymbol {m} / \mathrm {gcd} (\boldsymbol {m})$
. However, our aim here is to derive an asymptotic formula for
.
Lemma 9.1 We have
if
Proof By (6.14), with
$\theta = \varepsilon $
and
$U = \log T$
, we have
For t and v on the integral domains, in view of (6.9), we have by the Stirling formula that
Consequently,
The inner t-integral is a cosine Fourier integral of phase
$ 2 t \log \, \langle \boldsymbol {d} \rangle $
(see (6.17)); so, by repeated partial integration, we infer that it is negligibly small since
by (9.9).
By Lemma 9.1, the off-diagonal terms with
$ \delta (\boldsymbol {d}) = 0 $
contribute
$ O (\varPi T^{\varepsilon } / \sqrt {\|\boldsymbol {m^{\star }}\| }) $
since our assumption
$ \max \{ \boldsymbol {m} \} < \varPi ^{1-\varepsilon }$
yields
$ |\langle \boldsymbol {d} \rangle - 1|> \varPi ^{\varepsilon }/ \varPi $
.
It is left to consider those diagonal terms with
$ \delta (\boldsymbol {d}) = 1 $
. Note that, in this case,
$ \boldsymbol {h^{\star }} = \boldsymbol {m^{\star }} $
. More explicitly, (9.7) is turned into

where, after the removal of
$\tanh (\pi t)$
,
As, in (1.11), we set
It follows from (6.7), (9.10), and (9.11) that

Now we shift the integral contour from
$\mathrm {Re}(v) = 3$
down to
$\mathrm {Re} (v) = \varepsilon -1/2$
and collect the residue at the double pole
$v = 0$
. In view of (6.9) and (6.15), the resulting integral has bound
$O ((|t|+1)^{\varepsilon -1/2} )$
and it yields the error
$O (\varPi T^{1/2+\varepsilon })$
, while we have Laurent expansions (at
$v = 0$
):
where
$\gamma $
is the Euler constant and
$\psi (s) = \Gamma '(s)/\Gamma (s)$
is the di-gamma function, and it follows that the residue is equal to
By the Stirling formula,
Next, the t-integrals may be evaluated easily by
It is clear that the error after this process is
$O (\varPi ^3T^{\varepsilon }/ T {\displaystyle \sqrt {r }})$
, while the main term reads
Moreover, by the Gauss di-gamma formula in [Reference Magnus, Oberhettinger and SoniMOS, Section 1.2],
Finally, we conclude with the following asymptotic formula for
.
Lemma 9.2 Set
Define
Then,

9.3 Estimation of the off-diagonal sums
Similar to the initial arguments for
in Section 8.3, by applying Lemma 6.3, Corollary 4.1, and Remark 4.1 to (9.4)–(9.6), we have

where
$\boldsymbol {N} $
are dyadic, and

in which
the weight
$\textit {w} (\boldsymbol {x})$
is of the product form as in (6.20) for the weight function
$\textit {w} \in C_c^{\infty } [1, 2]$
as in Section 8.3, and
$H_{2}^{ \pm } (x, y ) $
are the Bessel transforms of
as in (4.3); here
$\delta $
and v have been suppressed again from the notation.
By the formula of the Bessel integrals
$H^{ \pm } (x, y ) $
as in Lemma 4.2, along with (4.4) and (4.9), it follows that
splits into four similar parts, one of which is

with

where
for
Note that the other
are defined similarly with
$ \textit {w} $
changed into
$\textit {w}^{_{\phantom {i}}}_{^\pm } $
(
$\textit {w}^{_{\phantom {i}}}_{^+} = \textit {w} $
and
$\textit {w}^{_{\phantom {i}}}_{^-} = \overline {\textit {w} }$
as in Section 8.3) or with
$ \unicode{x3c1} _{1} $
and
$ \unicode{x3c1} _{2} $
interchanged.
Now let us assume
$ m_1 N_1 \geqslant m_2 N_2 $
, so that
Our strategy is to apply Poisson summation to the
$n_1$
-sum but trivial estimation to the
$n_2$
-sum. Note that in the worst-case scenario, the
$n_2$
-sum could have bounded length, rendering trivial estimation optimal.
More explicitly, after opening the Kloosterman sum
$S (n_1, \pm n_2; c)$
in (9.22) into an
$\alpha $
-sum modulo c (see (2.10)), the
$\alpha $
- and
$n_1$
-sums read
By applying the Poisson summation formula modulo
$c d_2$
as in Lemma 6.1, this is turned into
Note that the Fourier transform is negligibly small unless
so the
$n_1$
-sum may be restricted effectively to the range
. However, our assumption
$ m_1 < \varPi /\varPi ^{\varepsilon } $
implies that only the zero frequency survives if we adjust
$\varepsilon $
by our
$\varepsilon $
-convention.
Given that
$n_1 = 0$
, from (9.25) and (9.26), we obtain the congruence equation
and the shrunk integral domain
Note that the former has at most
$ (c, h_1) $
many solutions of
$\alpha $
modulo c. Moreover, a trivial estimation of the
$n_2$
-sum yields a factor
$O (N_2/d_2)$
. Consequently,
Next, we insert this into (9.20) and bound the r-integral by
$O (cd_1 d_2 T^{\varepsilon } / m_1 N_1)$
, which is the length of the new integral domain defined by (9.27). At any rate, we have

Of course, the bound in the case
$ m_1 N_1 < m_2 N_2 $
is the same, and hence

Thus, by (9.17) and (9.29), we conclude with the following estimate for
.
Lemma 9.3 We have
10 The mollified moments
Historically, the mollifier method was first used on zeros of the Riemann zeta function by Bohr, Landau, and Selberg [Reference Bohr and LandauBL, Reference SelbergSel]. In this section, we shall mainly follow the ideas and maneuvers in [Reference Iwaniec and SarnakIS, Reference Kowalski and MichelKM2].
10.1 Setup
We introduce the mollifier
Assume that the coefficients
$ x_m $
are bounded, real-valued, and normalized so that
supported on square-free
$m $
with
$m \leqslant M$
. Note that
For
$\unicode{x3D5} (t)$
given in (7.1), consider the mollified moments
By the Cauchy inequality,
where, of course, the
$\delta $
indicates summation over
. Our task is to calculate the mollified moments
and
and then carefully choose the coefficients
$x_m$
to maximize the ratio
.
10.2 The mollified first moment
By the definitions in (1.8), (10.1), and (10.4), we have
Thus, if we insert the asymptotic formula of
in Theorem 5, then we obtain the main term
$\varPi T /\pi \sqrt {\pi }$
, while, as
$x_m$
are bounded, we estimate the error term by
Lemma 10.1 Let
$\varPi = T^{ \nu }$
and
$M = T^{\varDelta }$
. Then, for
$\nu < 1$
and
$\varDelta < (2\nu + 1)/3,$
there exists
$\varepsilon = \varepsilon (\nu , \varDelta )> 0$
such that
10.3 The mollified second moment
In view of the definitions in (1.9), (10.3), and (10.5), we have
Assume that
$ M < \varPi ^{1-\varepsilon }$
. Now we insert the asymptotic formula of
in Theorem 6 and bound the error term by the sum of
and

Then, we arrive at the asymptotic formula for
in the next lemma.
Lemma 10.2 Let
$\varPi = T^{ \nu }$
and
$M = T^{\varDelta }$
. Then, for
$\nu < 1$
and
$\varDelta < \min \{ 1/2, \nu \} $
, there exists
$\varepsilon = \varepsilon (\nu , \varDelta )> 0$
such that
with

where
$T_{\delta }$
is a multiple of T defined by
By removing the co-primality condition
$(m_1^{\star }, m_2^{\star }) = 1$
in (10.9) by the Möbius function, we have

By the definitions of
$\sigma _{-1} (m)$
and
$\sigma _{-1}' (m)$
in (1.14), some simple calculations yield
Consequently,

Thus, if we set
then (10.11) reads

By Möbius inversion, it follows from (10.12) that
where
$\varLambda (n) $
is the von Mangoldt function. Note that the assumption
$x_1 = 1$
in (10.2) now becomes the linear constraint
By Cauchy, the optimal choice of
$ y_{h} $
which minimizes the first quadratic form
in (10.13) with respect to (10.15) is given by
We have
Moreover, by (10.14) and (10.16), it is clear that our
$x_{m}$
are indeed bounded.
By definition,
By the prime number theorem (see [Reference TitchmarshTit, Formula (3.14.3)]),
and, for h square-free,
It follows from (10.14), (10.16), (10.19), and (10.20) that
Consequently,

and a partial summation using (10.17) yields

By (10.8), (10.13), (10.17), (10.18), and (10.21), we conclude with the following asymptotic formula for
.
Lemma 10.3 Let
$\varPi = T^{ \nu }$
and
$M = T^{\varDelta }$
. Then, for
$\varDelta < \min \{ 1/2, \nu \}, $
we have

10.4 Conclusion
For
$0 < \nu < 1,$
let
$\varPi = T^{ \nu }$
. By Lemmas 10.1 and 10.3, it follows from Cauchy in (10.6) that
as
$T \rightarrow \infty $
for any
$\varepsilon> 0$
. Next, we apply Lemma 5.1 with
$a_n = \delta (L (s_f, f ) \neq 0)$
(namely, the Kronecker
$\delta $
symbol) and
$H = \varPi ^{1+\varepsilon }$
, obtaining
for any
$ T^{ \varepsilon } \leqslant H \leqslant T/3$
. Recall the Weyl law (see (2.13)):
as a consequence,
By (10.24)–(10.26), we deduce the following harmonic-weighted non-vanishing results.
Theorem 10.4 Let
$\varepsilon> 0$
be arbitrarily small. As
$T \rightarrow \infty $
, we have
and, for
$0 < \mu < 1$
,
where the sums are over
.
Note that (10.27) follows from (10.24) and (10.25) with
$ H = T/3 $
along with a dyadic summation, while (10.28) follows directly from (10.24) and (10.26).
Since the harmonic weight
one may remove
$\omega _f$
without compromising the strength of results by the method of Kowalski and Michel [Reference Kowalski and MichelKM1], adapted for the Maass-form case in [Reference Balkanova, Huang and SödergrenBHS] (see also the paragraph below (2.9) in [Reference Iwaniec and SarnakIS]). Thus, after the removal of
$\omega _f$
, we obtain the non-vanishing results in Theorems 1 and 2.
11 Moments with no twist and smooth weight
Theorem 11.1 Let
$ T^{ \varepsilon } \leqslant H \leqslant T/3$
. Then,
Proof By applying Lemma 5.1 to the first asymptotic formula in Corollary 7, we have
for any
$T^{\varepsilon } \leqslant \varPi ^{1+\varepsilon } \leqslant H \leqslant T / 3 $
, and we obtain (11.1) on the choice
$ \varPi = H / \sqrt {T} + T^{\varepsilon } $
.
Theorem 11.2 Let
$ T^{ \varepsilon } \leqslant H \leqslant T/3$
. Then,
Proof By applying Lemma 5.1 to the second asymptotic formula in Corollary 7, we infer that the spectral second moment on the left of (11.2) equals
for any
$\varPi ^{1+\varepsilon } \leqslant H \leqslant T/3$
, and we obtain (11.2) on the choice
$\varPi = \sqrt {H}$
.
Corollary 8 is a direct consequence of the above asymptotics with
$H = T/3$
along with a dyadic summation.
Part II
Density theorems
12 Sums of variant Kloosterman sums
Let us introduce the variant Kloosterman sum
The definition of
$V_{q} (m, n; c)$
is essentially from Iwaniec and Li [Reference Iwaniec and LiIL, Formula (2.17)] and arose in the identity of Luo [Reference LuoLuo2, Section 3]:
The purpose of this section is to prove the following analog of [Reference Iwaniec, Luo and SarnakILS, Lemma 6.1].
Lemma 12.1 Assuming the Riemann hypothesis for Dirichlet L-functions, for any
$x \geqslant 2$
, we have
where
$\varphi (c) $
is the Euler function,
$R (m; c) = S (m, 0; c)$
, and
$R_q (n; c) = V_q (0, n; c) $
.
Proof As in [Reference Iwaniec, Luo and SarnakILS, Section 6], we shall use the Riemann hypothesis for Dirichlet L-functions in the form
where
$ \chi $
is any character to modulus c and
$\delta ({\chi })$
is the indicator of the principal character. By definition,
and, on the change
$ \alpha \rightarrow \alpha \beta , $
it splits into the product
$ G (m; \chi ) G_{ q} (n; \chi )$
with
Note that for the principal character, these (Gauss) sums become the (Ramanujan) sums
$R (m; c)$
and
$R_q (n; c)$
, and that by the orthogonality of characters,
Finally, from these formulas, we derive
and hence (12.3).
Later in our application,
$m = 1$
and
$n = \pm 1$
. It is well known that the Ramanujan sum
By the Chinese remainder theorem, it is easy to prove
13 Further analysis for the Bessel integrals
Set
$\unicode{x3b3} (t) \equiv 1$
in (4.3) so that
According to (4.4) and (4.7), up to a negligible error, we split
$H_2^{\pm } (x, y) $
into four similar parts, one of which reads
with
for
$ \unicode{x3c1} (r) = \exp (r) - 1. $
More explicitly, the pair of
$\pm $
in the subscript of
$I^{_\pm }_{^{\pm \pm } } (v, w) $
is determined by that in
$ \pm i \unicode{x3c1} (\pm r) v $
contained in the integral. Later, we shall have
and
so, in practice,
$w $
is small and inessential and will be suppressed from the notation.
Lemma 13.1 Write
Then, we have bounds
Moreover, if
$ w < \varPi $
, then
$I^{_\pm }_{^{-+} } (v, w), I^{_{\pm \natural }}_{^{-+} } (v, w) = O (T^{-A})$
unless
.
Proof Note that
$I^{_{\pm \natural }}_{^{-+} } (v, w) $
is very similar to
$I^{_{\pm }}_{^{-+} } (v, w) $
:
Clearly, (13.4) follows from (13.2) and (13.5) by trivial estimation. As for the second statement, we use the method of stationary phase. By
$w < \varPi $
, we may as well absorb the factor
$ \exp (\mp i \unicode{x3c1} (-r)w)$
into the weight
$g (\varPi r)$
so that the phase function is
$ 2 T r - \unicode{x3c1} (r) v $
, while its derivative equals
Thus, in the case
, by an application of Lemma B.1 in [Reference QiQi] with
$S = 1$
,
$ P = 1/\varPi $
,
$Q = 1$
,
$Z = T $
, and
$ R = \varPi ^{\varepsilon } ( \varPi +T / \varPi ) $
, we infer that the integrals in (13.2) and (13.5) are negligibly small. Note here that
$R $
is larger than both
$ \varPi ^{1+\varepsilon }$
and
$ T^{1/2+\varepsilon }$
. However, for
, the integrals are always negligibly small if we apply Lemma B.1 in [Reference QiQi] with
$S = 1$
,
$ P = 1/\varPi $
,
$Q = 1$
,
$Z = v$
, and
$ R = v + T $
.
Remark 13.1 As for the other three integrals
$I^{_\pm }_{^{+-} } (v, w) $
,
$I^{_\pm }_{^{+ +} } (v, w) $
, and
$I^{_\pm }_{^{- -} } (v, w) $
, the first one is similar, while the latter two are always negligible.
14 The explicit formula
By applying [Reference Rudnick and SarnakRS, Proposition 2.1] with
we obtain the explicit formula
with
as defined before in (1.19), and
where
(see (3.3) and (3.4)), as usual,
$\Gamma _{\mathbf {R}} (s) = \pi ^{-s/2} \Gamma (s/2) $
,
$\Lambda (n) $
is the von Mangoldt function (
$\Lambda (n) = \log p$
if
$n = p^{\nu }$
and
$0$
otherwise), and
for
$\alpha _f (p)$
and
$ \beta _f (p)$
defined as in the Euler product (3.6).
Lemma 14.1 Let
$\mathrm {Re} (s) \geqslant \delta> 0$
. Then, for
, we have
Proof By definition,
$\Gamma _{\mathbf {R}} (s) = \pi ^{-s/2} \Gamma (s/2) $
, and by the change
$r \rightarrow 2\pi r/ \log T$
, the integral in (14.5) is turned into
It follows from the Stirling formula
and the crude uniform bound
that the integral above equals
as desired.
Remark 14.1 By modifying the proof of Lemma 14.1, one may improve the error term
$O_{\delta } (1/\log T)$
into
$ O_{\delta , \varepsilon } (1 / |s| \log T) $
if
.
By applying Lemma 14.1 with
$s = \kappa _f$
or
$s = \kappa _f + it_f$
to (14.3), we have
Lemma 14.2 Assume the Riemann hypothesis (for the Riemann zeta function). Let
$ s = 1 + it $
and
. Then,
for
$ x < t$
, where the implied constant is absolute.
Proof As usual, write
$s = \sigma + it$
. By [Reference TitchmarshTit, Sections 14.5 and 14.33], under the Riemann hypothesis, we have
uniformly for
$1/2 < \sigma _0 \leqslant \sigma \leqslant \sigma _1 < 1$
, while
for
$ \sigma = 1$
. By (14.8), it follows from the Perron formula in [Reference TitchmarshTit, Lemma 3.12] and the contour-shift argument in [Reference TitchmarshTit, Sections 3.14 and 3.15] with
$c = 1/ \log x$
and
$\delta = 1 - \sigma _0 $
that
in the case that
$\sigma = 1$
and
$T < t/2$
; here, one may let x be half-integral. For
$x < t$
, on choosing
$\sqrt {T} = \log x + \log ^{\delta } t$
, we obtain
Note here that
Also, trivially, we have
Accordingly, as
$x> \log ^2 t$
or not, the bound (14.7) follows from (14.9), (14.10), or (14.11).
Now we turn to the sum
$P ( f; \unicode{x3c6}; T)$
as defined in (14.4): unlike the cases of central L-values in [Reference Alpoge and MillerAM, Reference Iwaniec, Luo and SarnakILS, Reference Liu and QiLQ1], this will not contribute a main term. Suppose that
$\hat {\unicode{x3c6} }$
has support in
$(-2, 2)$
, say, and
$ t_f $
is close to T.
Firstly, we may discard those terms with
$n = p^{\nu }$
for
$\nu \geqslant 3$
since
$|a_f (p^{\nu }) | \leqslant 2 p^{\frac 7 {64} \nu }$
by the Kim–Sarnak bound (3.7), and hence their contribution is
$O (1/\log T)$
. Next, by partial summation, it follows from Lemma 14.2 that
Since
$a_f (p) = \lambda _{f} (p)$
and
$ a_f (p^2) = \lambda _{f} (p^2) - 1 $
, we conclude that
with
Finally, we obtain the following asymptotic formula by combining (14.1), (14.6), and (14.13).
Lemma 14.3 We have
if
$\hat {\unicode{x3c6} }$
has support in
$(-2, 2)$
.
15 Proof of the density theorems
As in (4.2), let us define
and consider the
-averaged
$1$
-level density
Clearly,
differs from
as in (1.20) only by an exponentially small error term. By inserting (14.15), we have

with
Note here that the main term
$ \varPi T \hat {\unicode{x3c6} } (0) / \pi \sqrt {\pi } $
is due to the Weyl law as in (2.13). Therefore, Theorems 9 and 10 are deducible easily from the following bounds for
and
.
Lemma 15.1 We have
if
$\hat {\unicode{x3c6} }$
has support in
$(-1, 1)$
.
Lemma 15.2 Let
$\varPi = T^{\mu }$
for
$0 < \mu < 1$
. Assume the Riemann hypothesis for Dirichlet L-functions. Define
Then, we have
if
$\hat {\unicode{x3c6} }$
has support in
$(- v(\mu ), v (\mu ))$
.
Subsequently, we shall focus on
as the treatment of
is much easier and yields a better bound.
Notation
Suppose
$\hat {\unicode{x3c6} }$
has support in the closed interval
$ [-v, v] $
, say for
$v < 2$
. Let
$\varPi = T^{\mu }$
and
$P = T^{v}$
.
For brevity, let us write
$ \unicode{x3D5} (t) = \unicode{x3D5} _{T, \varPi } (t)$
as before and suppress
$T, \varPi $
from the notation like
.
15.1 Application of the Kuznetsov formula
By (14.14) and (15.4), we may write

where, as in (4.3),
By applying the Kuznetsov trace formula in (2.9), we have
where


First of all, the Eisenstein term
is indeed small. Since the integral in (15.9) is trivially
$O (\varPi \log ^2 T)$
(by [Reference TitchmarshTit, (3.11.10)]), we have
Note that the similar Eisenstein term
is even smaller:
$ O (\varPi \log ^2 T) $
.
It is now left to consider the off-diagonal sums
as in (15.10).
15.2 Proof of Lemma 15.1
It follows from Corollary 4.1 that the Bessel integral
$ H_2^{\pm } ( {4\pi \sqrt {p}} / {c} , \sqrt {p} )$
in (15.10) is negligibly small for
. Therefore,
is negligible for
$v < 1$
so that
$ p / T = o (1) $
for all
$p \leqslant P = T^{v}$
. Similarly, the off-diagonal contribution to
is also negligible for
$v < 1$
.
Note that, in the odd case,
and
are actually negligibly small as there is no Eisenstein contribution.
15.3 Proof of Lemma 15.2
In light of (15.8) and (15.11), it suffices to prove
for any
$v < v (\mu )$
. As in Section 13, let us split
into four similar parts, one of which contains
$ I^{_\pm }_{^{- +} } (2\pi p/c, 2\pi / c) $
as in (13.1) and (13.2) and reads

For brevity, let us write
$ I^{_\pm }_{^{- +} } ( 2\pi p/c ) = I^{_\pm }_{^{- +} } (2\pi p/c, 2\pi / c) $
. By the Luo identity (12.2),

Next, we change the order of summations and truncate the c- and q-sums at
(again by Corollary 4.1), obtaining

where

In the asymptotic formula (12.3), we drop the condition
$p \nmid c$
and change
$\log p$
into
$ \log p /\sqrt {p}$
by partial summation, and it follows that
By (15.13) and (15.14), we get

On partial integration for the main-term integral and the change of variable
${v = 2\pi x/ cq}$
, in view of (12.4) and (12.5), we deduce that

where
$ I^{_{\pm \natural }}_{^{- +} } ( v ) $
is the (partial) derivative of
$ I^{_\pm }_{^{- +} } ( v ) $
defined as in (13.3). By Lemma 13.1, the domain of integration may be restricted to
, and it yields the bound:

Now we insert (15.15) into (15.12),

Finally, we verify
in various cases. Keep in mind that
$\varPi = T^{\mu }$
and
${P = T^{v}}$
. For
$1/2 < \mu < 1$
,
and this is
$o (\varPi T)$
as long as
${v < 1 + \mu }$
. For
$1/3 < \mu \leqslant 1/ 2$
,
and this is satisfactory for
$v < 3 \mu $
. For
$0 < \mu \leqslant 1/3$
,
is negligible for
$v < 1$
as the c-sum is void (of course, in this case, (15.7) also follows from (15.5) in Lemma 15.1). Moreover, in view of Remark 13.1, one may treat the other
in a similar or easier manner.
16 Proof of the non-vanishing results on the Riemann hypothesis
Let us modify the argument in [Reference Iwaniec, Luo and SarnakILS] as follows. Put
where the
$\delta $
in the sum is the Kronecker
$\delta $
-symbol (not to be confused with the parity
$\delta $
in
). Clearly,
On the other hand, we choose, as in [Reference Iwaniec, Luo and SarnakILS, Formula (1.42)], the Fourier pair
for every
$v < v (\mu )$
(see (15.6)). By the definition in (1.19), it follows from
$ \unicode{x3D5} (0) = 1 $
and the non-negativity of
$\unicode{x3c6} $
that
Also note that
Thus, for
$ 1/3 < \mu < 1$
, we derive from (1.20), (1.22), and (1.23) that
for any
$\varepsilon> 0$
, provided that T is sufficiently large. By subtracting (16.4) from (16.2), we obtain the lower bound
Note that (16.5) is non-trivial only if
$ v (\mu )> 1 $
; so, we do need
$ 1/3 < \mu < 1$
as above. In view of (15.6) and (16.1), let us reinterpret (16.5) as

Similar to the argument in Section 10.4, we may remove the smooth and harmonic weights by Lemma 5.1 and the technique of Kowalski and Michel [Reference Kowalski and MichelKM1]. Consequently, we arrive at (1.6) and (1.7) in Theorems 3 and 4.
A Addendum: Non-vanishing of central L-values on the Riemann hypothesis
In this addendum, we take the opportunity to improve Theorem A.3 in [Reference Liu and QiLQ1] as follows, with the restriction
$\mu> 1/3$
removed, and also to establish conditional non-vanishing results by (1.40)–(1.45) in [Reference Iwaniec, Luo and SarnakILS] for the central L-values
$L (1/2 , f)$
or
$L' (1/2, f)$
in the even or odd case.
Theorem A.1 Let the notation be as in the Appendix in [Reference Liu and QiLQ1]. Assume the Riemann hypothesis for Dirichlet L-functions. Let
$ T, M> 1$
be such that
$M = T^{\mu }$
with
$0 < \mu < 1$
. Define
Let
$\unicode{x3c6} $
be an even Schwartz function with the support of
$\hat {\unicode{x3c6} }$
in
$(-v (\mu ), v (\mu ) )$
. Then,
Proof It only requires minor modifications to the proof of Theorem A.3 (indeed, Lemmas A.8 and A.9 in [Reference Liu and QiLQ1]).
As in (A.21) and (A.23) in [Reference Liu and QiLQ1], we need to analyze integrals of the form:
where
$ \widehat {k} (- r/\pi ) $
is a Schwartz function similar to the
$g ( r )$
as in Lemma 4.2. As in the second statement of Lemma A.8 in [Reference Liu and QiLQ1], under the condition
$ \mu> 1 / 3 $
, it is proven that
$H_{T, M}^{- \sharp } (x) $
and
$ H_{T, M}^{- \natural } (x) $
are negligibly small unless
$ |T-2\pi x| \leqslant M^{1+\varepsilon } $
. The reader should compare these integrals with those in (13.2) and (13.5). Note that the derivatives of the phase functions here satisfy
For this theorem, it suffices to remove the condition
$\mu> 1/3$
but weaken the range
$ |T-2\pi x| \leqslant M^{1+\varepsilon } $
into
$ |T-2\pi x| \leqslant M^{\varepsilon } (M + T/M^2) $
. To this end, we may apply Lemma B.1 in [Reference QiQi] with
$S = 1$
,
$ P = 1/ M $
,
$Q = M^{\varepsilon } / M$
,
$Z = T / M^3 $
, and
$ R = M^{\varepsilon } ( M +T / M^2 ) $
. This is essentially in the spirit of our Lemma 13.1. Now, for
$P = T^{2v}$
, it is easy to verify that the error bound
$O (\sqrt {P} T^{\varepsilon })$
becomes
$ O (\sqrt {P} (1+T/M^3) T^{\varepsilon }) $
(for
; otherwise, the c-sum is void) in the first part of the proof of Lemma A.9 in [Reference Liu and QiLQ1].
Note that in the notation of this article,
,
, and
$\varPi = M$
.
Next, we assume the Riemann hypothesis for every
$L (s, f)$
with
$ f \in H^{\pm } (1) $
and for all Dirichlet L-functions. The following results may be derived from Theorem A.1 by adapting the arguments from Sections 10.4 and 16 and also [Reference Iwaniec, Luo and SarnakILS]. The reader may compare Corollary A.1 with Corollary 1.7 in [Reference Iwaniec, Luo and SarnakILS].
Corollary A.1 We have
Corollary A.2 We have
if
$ 0 < \mu \leqslant 1/4$
, and
if
$1/4 < \mu < 1$
, for
Acknowledgments
The author gratefully acknowledges the referees for identifying several errors in the manuscript and offering numerous valuable suggestions that substantially improved the exposition.











