A discrete mean-value theorem for the higher derivatives of the Riemann zeta function
Abstract.
We show that the th derivative of the Riemann zeta function, when summed over the non-trivial zeros of zeta, is real and positive/negative in the mean for odd/even, respectively. We show this by giving a full asymptotic expansion of these sums.
1. Shanks’ Conjecture and statement of results
Shanks’ Conjecture from Shanks [12] says for a zero of , that
This conjecture was first proven by evaluating the discrete first moment of to leading order in Conrey, Ghosh, Gonek [2] and for all other terms in Fujii [3], with a minor correction to its lowest order terms in Fujii [4]. We state the precise asymptotic formula shortly, in Theorem 1. Other proofs of Shanks’ Conjecture have since been given (either explicitly or implied), for example, in Trudgian in [15] and in Stopple [13].
In this paper we will generalise this to all higher derivatives, by giving explicit asymptotic formulae (including the lower order terms) for the sums of the th derivatives of evaluated at the zeros of . From the results of Theorem 3 it will follow that
Throughout we assume is sufficiently large and that , where is the imaginary part of any zero . We write where for a general complex number . We let and be the coefficients in the Laurent expansion of about , which is given by
It is clear that Shanks’ conjecture follows from the following result.
Theorem 1 (Fujii).
With the notation and assumptions given above, we have
where, unconditionally,
where is a positive constant. Conditionally upon the Riemann Hypothesis we have
The aim of this paper is to generalise the earlier result on Shanks’ Conjecture to all th derivatives of zeta. Kaptan, Karabulut, Yıldırım [9] found the main term for the mean of .
Theorem 2 (Kaptan, Karabulut and Yıldırım).
With the notation and assumptions given above, we have
Similar to us, and all other previous papers on this and other calculations, they start by rewriting the sum as a Cauchy residue integral, but they only proceed to calculate the main term in their subsequent calculations, while we also calculate the full asymptotic and error terms. (In fact, we only became aware of their paper while writing up our results).
To give the lower order terms in the expansion, in addition to the coefficients defined previously, let be the coefficients in the Laurent expansion for about , given by
Note that the are related to the by the following recursive formula, as shown in Israilov [7], with equivalent forms in Maslanka [10] and Bomberi [1]:
Remark.
Our constants are related to the Stieltjes constants which are normally used in the Laurent expansion for about by
We have decided to use for simplicity in our formulae, and to remain consistent with Fujii’s papers on Shanks’ Conjecture. Clearly our coefficients are then also related to .
Theorem 3.
With the notation and assumptions given above, we have
where
with is a positive constant. If we assume the Riemann Hypothesis, then
In the case , we are able to show that under RH,
2. Overview of the Proof
In this section we give a brief outline of the paper, including the method we will follow to prove Theorem 3.
In Section 3, we state some preliminary lemmas that we will use throughout the paper.
In Section 4, we begin by considering the integral given by
where is the positively oriented rectangular contour with vertices with . The non-trivial zeros of up to a height are contained within and so by Cauchy’s Theorem the integral represents the summation
in question. We then show that the contribution to the integral from the bottom, top and right-hand side of the contour is contained within the error term of the theorem, so the main contribution comes from the integral along the left-hand side of the integral.
In Section 5, we evaluate this part of the contour integral. Through this, in Lemma 10 we link the integral to a summation
which is then the main object in question.
In Section 6, Lemma 11 we use Perron’s formula to evaluate this sum as a complex integral up a vertical line just to the right of the critical strip.
In Section 7, we will evaluate that complex integral and show that
where and is given explicitly to optimise the error term , which we can describe both unconditionally and under the assumption of RH. We will also show, via a different argument, a slight improvement on the error term for .
Finally, in Section 8, we will evaluate the residue at to find the asymptotic expansion described in Theorem 3 by considering the Laurent expansions around the pole of the terms in the above integrand.
Combining all of these steps together gives the result.
We finish the paper by giving some examples in Section 9. Firstly, we can recover the results in Theorem 1 by specialising to the case , and in particular recover Shanks’ Conjecture. We then give the first new case where , as well as providing some good numerical evidence for the strength of the result.
3. Preliminary Lemmas
We start by stating the strong form of the convexity bound for the zeta function and the th derivatives of the zeta function. The case can be found in Ivić [8, Ch. I, Sect. 1.5], while the cases for all other can be derived from this using Cauchy’s theorem on derivatives of analytic functions on in a small disc of radius centred at .
Lemma 1.
For uniformly in ,
The following result follows from Gonek [6, Sect. 2, p. 126].
Lemma 2.
If is such that for any ordinate , uniformly for we have
(3.1) |
The functional equation for can be given by
(3.2) |
where
and is the Gamma function. Taking the logarithmic derivative gives
(3.3) |
To obtain the form of the functional equation for that we will need, we will use the following result from Gonek [6], which is his Lemma 6.
Lemma 3.
For fixed, and we have
Lemma 4.
For and
Proof.
4. Beginning the proof
Throughout we assume is sufficiently large and satisfies , where is the imaginary part of any zero . (This constraint simplifies the arguments, but has no effect on the resulting expressions.) Consider the integral
(4.1) |
where the contour is the positively oriented rectangular contour with vertices with .
By Cauchy’s Residue Theorem,
We need to evaluate in another way to determine the behaviour of . To do this, we begin by splitting the integral along each part of the contour, so
say. We will first bound trivially within error terms as follows. Our main aim will then be to evaluate which we will do in the next section.
Lemma 5.
The integral along the bottom of the contour is .
Proof.
This follows as the integral is of finite length with a bounded integrand. ∎
Lemma 6.
The integral along the top of the contour is .
Proof.
Lemma 7.
The integral along the RHS of the contour is .
Proof.
Since are all within a remainder term of , only the integral over the LHS of the contour will contribute in any meaningful way. Observe that
where
(4.3) |
5. Evaluating
Our main aim then for this section is to evaluate . To do this we want to write in such a way that we can apply Gonek’s Lemma 5 from [6] which we will state below when we use it later in this section.
Lemma 8.
We can write the integral given in (4.3) as
Proof.
We will now require the use of Lemma 5 from Gonek [6, Sect. 4, p. 131]. As it is such an important result for our proof, we state it here, with the phrasing adapted to suit our needs.
Lemma 9 (Gonek).
Let be a sequence of complex numbers such that for any , . Let be as before and let be a non-negative integer. Then for sufficiently large,
We can finally simplify to get a single sum that we will work on evaluating in the next section.
Lemma 10.
The integral can be written as
where the sum is taken over all integers and such that .
Proof.
Before we apply Lemma 9, we will split the integral given in (5.1) substituting the Dirichlet series to simplify the argument slightly. To do this, we write
say.
For , we have
since by the Binomial Theorem,
For , we have
where the last three sums are over all integers and such that . Our results follow from combining and . ∎
Therefore
Since we have shown that the terms from are harmless within the error term of , we have that the integral given in equation (4.1) is equal to the sum
(5.2) |
6. Evaluating the sum
As we have discussed above, all that remains to do is to evaluate the sum in (5.2). To do this, we first note that by Perron’s formula [11, Ch. 5, Sect. 5.1], we have
where we set and as before.
Since we want to be able to evaluate the integral on the RHS, we will modify this slightly and instead use a truncated Perron formula.
Lemma 11.
For , as ,
where
We will use a specific later in this section. The choice of will depend on whether we assume RH or not.
Proof.
If we let denote the coefficients in Dirichlet series for , namely
then by the truncated Perron formula [11, Ch. 5, Sect. 5.1] we have that the integral in the Lemma equals
with
Writing we see equals the sum in the Lemma. To evaluate the error, note that (with equality only if is prime). Therefore
For , since , the Dirichlet series converges, and so
where we use the fact that .
Combining and gives the required bound on . ∎
7. The Error Term
We will show that the error term can be described explicitly as follows, depending on whether we assume RH or not.
Lemma 12.
For , as , we have
where is an error term given by one of the following two cases:
-
(1)
Unconditionally, by setting , we obtain .
-
(2)
Assuming RH, setting , we obtain .
where a positive constant that is not necessarily the same in each instance.
We will calculate the residue in Section 8. First, we will show how we can obtain the different expressions for in the following subsections. We will also show in the last subsection that for , we are able to obtain a slightly better error term than that given above or in Fujii [3].
7.1. The Unconditional Case
From Titchmarsh [14, Sect. 3.8, p. 54], we know there is a positive constant such that for , all the zeros of are away from the line running from to . By Cauchy’s Residue Theorem, the integral is
By Lemmas 1 and 2, if , the integral on the horizontal lines can be estimated as
where we use the fact that and since .
For the integral on the vertical line, we have
Since , balancing the two errors comes from taking for some positive constant , and so we have
7.2. The Conditional Case
Throughout this subsection we assume the Riemann Hypothesis, RH.
One approach would be to choose to guarantee a zero-free region, and set . In that case, an application of Cauchy’s Residue Theorem yields
for , where the horizontal pieces of the contour are estimated in a manner similar to that below, and the vertical piece of the contour uses the bound for .
However, we can get a better bound, one that depends explicitly upon by choosing (that is, just to the left of the critical strip) and . By Cauchy’s Residue Theorem,
where the sum runs through , the zeros of , lying inside the contour.
To estimate the integral on the horizontal lines, choose so that all the zeros of zeta are bounded by away from the horizontal line. Therefore we may use Lemma 2 to bound the logarithmic derivative by along the line. Using the convexity result from Lemma 1, we obtain, unconditionally,
For the integral on the vertical line, since and we may again use Lemma 2 as the vertical line is bounded away from any zeros of zeta. Using Lemma 1 to bound zeta just to the left of the critical strip, we obtain (again unconditionally),
We now consider the poles at for each with , where is a zero of .
First we note that for and , we have
Being careful with the error term, we may differentiate this times to give
For each , we need to consider the coefficient of in the expansion of to find the residue at . For this, note that by the triple product rule, we may write
where is the multinomial coefficent given by
Therefore, at each zero ,
As we are just bounding these terms, we do not worry about the coefficient such as the , and the multinomial coefficients. Assuming RH, so and summing over all zeros with , we have
(7.1) |
Now assuming RH, by Gonek [6] we have that for all positive integers ,
We will use this result together with the Cauchy-Schwarz inequality and partial summation to bound each term in the our summation. We see that the terms in the right hand side of (7.1) can be bounded by
Clearly the dominant error term is when (which forces ), and so the sum over the residues at the zeros in (7.1) is bounded by
Balancing this error term with the error term from the vertical line comes from taking , and so we have
7.3. An improvement on the error term in the conditional case for
In this section we assume RH. When we have an error bound of which agrees with Fujii’s result in [3].
From the previous section, we see that
Garaev [5] has shown that under RH, we have
Remark.
Garaev’s paper assumes additionally that all the zeros are simple. However, a close examination of the proof shows that this particular result does not need that additional assumption.
By partial summation
so setting to balance this error term with the error coming from the vertical line in the contour yields
an improvement of in our error term.
8. Finding the Leading Asymptotic Terms
We now evaluate
from Lemma 12. We expand each of the terms in this residue calculation in their Laurent expansions about . Since
we have
-
(1)
The Laurent expansion for about :
-
(2)
The Laurent expansion for about :
-
(3)
The Laurent expansion for about :
-
(4)
The Laurent expansion for about :
We now work out the residue at by considering the terms in powers of , as runs from down to . We consider a combinatorial style argument to make sure we consider all possible terms. The following tables show this for various powers of .
To calculate the resulting contribution to the residue for the leading order behaviour, that is, the coefficient of , there is only one way to obtain a factor of for the residue calculation. We have to have the term from , the term from , the term from and the term from . Combined, the coefficient for is .
To calculate the sub-leading term there are two distinct possibilities.
-
(1)
The first possibility is to have the term from , the term from , the term from and the term from . The resulting contribution to the residue for this term is .
-
(2)
The other possibility is to have the term from , the term from , the term from and the term from . The resulting contribution to the residue for this term is .
Combined, the coefficient for the term is
The other terms in all lower order cases follow in a similar manner — Table 1 shows all possibilities, with one extra caveat for the lowest order term which has an additional term.
Contribution to the residue | ||||
⋮ | ⋮ | ⋮ | ⋮ | ⋮ |
In the lowest-leading order case, the coefficient of the term, there are all the usual terms for the sub-leading behaviour, with one additional term. The additional term comes from the term from , the term from , the term from , and the term from . Combined, the extra term for has coefficient .
Combining all of the terms in the above tables, and putting we have shown that
9. Examples of first moments of specific th derivatives of the Riemann zeta function
The first test of a good general theorem is if it can recover any previously known results. We begin this section by doing exactly that in the case where we recover Shanks’ Conjecture and Fujii’s asymptotics (Theorem 1) from Chapter 1. After that we give another example, as well as showing some numerical data showing that our result is in excellent agreement with the true values. We have chosen to show the most simple new case in full detail. The calculations for higher level derivatives become unwieldy to do by hand so we could employ a computer package to manage this for us to obtain any specific th order derivative that we choose.
9.1. The case (Shanks’ Conjecture)
We can use Theorem 3 to recover Fujii’s asymptotic formula Theorem 1 (with an improved error term) and hence Shanks’ Conjecture from Section 1. To do this, set in the theorem to obtain
Note that
since (mentioned in the preamble before Theorem 3). Substituting this into the expression above and simplifying gives
in perfect agreement with Theorem 1 (with the improved error term from Theorem 3 under RH).
9.2. The case
We set in Theorem 3 to obtain
Note that
Substituting this into the expression above and simplifying gives
The error term is unconditionally, and assuming RH can be improved to .
We show some numerical evidence of the result for . The 100,000th zero has imaginary part roughly equal to , and the true value of is . Clearly the imaginary part is small as we expect, and we can tell by the graph that the result for the real part of the sum is extremely accurate even for this relatively small number of zeros, since the error never exceeds 2,400.

Acknowledgements
The authors are grateful to the anonymous referee for their comments on the paper. We are particularly thankful for their suggestions with regard to our calculation for , which shortened and neatened this argument considerably.
This work will form part of the second author’s PhD thesis at the University of York.
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
References
- Bombieri [1999] Bombieri, E., 1999. Complements to Li’s Criterion of the Riemann Hypothesis. J. Number Theory, 77, pp.274–287.
- Conrey, Ghosh and Gonek [1988] Conrey, J.B., Ghosh, A. and Gonek, S.M., 1988. Simple Zeros of Zeta functions. Colloque de Théorie Analytique des Nombres “Jean Coquet” (Marseille, 1985). Orsay: Univ. Paris XI, Publ. Math. Orsay, vol. 88, pp.77–83.
- Fujii [1994] Fujii, A., 1994. On a Conjecture of Shanks. Proc. Japan Acad., 70(4), pp.109–114.
- Fujii [2012] Fujii, A., 2012. On the Distribution of Values of the Derivative of the Riemann Zeta Function at Its Zeros. I. Proc. Steklov Inst. Math., 276, pp.51–76.
- Garaev [2003] Garaev, M.Z., 2003. One inequality involving simple zeros of . Hardy-Ramanujan J., 26, pp.18–22.
- Gonek [1984] Gonek, S.M., 1984. Mean values of the Riemann zeta-function and its derivatives. Invent. math., 75, pp.123–141.
- Israilov [1981] Israilov, M.I., 1981. The Laurent expansion of the Riemann zeta function. Trudy Mat. Inst. Steklov, 158, pp.98–104.
- Ivić [1985] Ivić, A., 1985. The Riemann Zeta-Function: Theory and Applications. Dover Publications, Inc.
- Kaptan, Karabulut and Yıldırım [2011] Kaptan, D., Karabulut, Y. and Yıldırım, C., 2011. Some mean value theorems for the Riemann zeta-function and Dirichlet -functions. Comment. Math. Univ. St. Pauli, 60(1-2), pp.83–87.
- Maslanka [2004] Maslanka, K., 2004. Effective method of computing Li’s Coefficients and their properties. Unpublished. arXiv:math/0402168v5.
- Montgomery and Vaughan [2006] Montgomery, H.L. and Vaughan, R.C., 2006. Multiplicative number theory I. Cambridge University Press.
- Shanks [1961] Shanks, D., 1961. Review of ‘Tables of the Riemann zeta function’ by C.B. Haselgrove in collaboration with J.C.P. Miller. Math. Comp., pp.84–86.
- Stopple [2020] Stopple, J., 2020. Notes on the phase statistics of the Riemann zeros. Unpublished. arXiv:2007.08008.
- Titchmarsh [1986] Titchmarsh, E.C., 1986. The Theory of the Riemann Zeta-Function. Clarendon Press.
- Trudgian [2010] Trudgian, T.S., 2010. On a conjecture of Shanks. J. Number Theory, 130(12), pp.2635–2638.