Computational SectionNumerical verification of Littlewood's bounds for |L(1,χ)|
Introduction
Let q be an odd prime, χ be a Dirichlet character and be the associated Dirichlet L-function. The goal of this paper is to introduce a fast algorithm to compute the values of for every non-principal primitive Dirichlet character χ defined and, using such a new method, to numerically study a generalisation of the classical bounds of Littlewood [12] for , where is a quadratic Dirichlet character. Assuming the Riemann Hypothesis for holds, in 1928 Littlewood proved, for , that as d tends to infinity, where γ is the Euler-Mascheroni constant. In 1973 Shanks [16] numerically studied the behaviour of the upper and lower Littlewood indices defined as for several small discriminants d. Such computations were extended by Williams-Broere [17] in 1976 and by Jacobson-Ramachandran-Williams [5] in 2006.
Recently Lamzouri-Li-Soundararajan [8, Theorem 1.5] proved an effective form of Littlewood's inequalities: assuming the Generalised Riemann Hypothesis holds, for every integer and for every non-principal primitive character , they obtained that and
Using our method we will compute the values of for every non-principal primitive Dirichlet character χ defined , for every odd prime q up to 107. This largely extends previous results. Moreover, letting we obtain the following Theorem 1 Let , q be a prime number and be defined in (4). Then we have . Moreover, we also have where the lower bound holds just for , and where the upper bound holds just for .
We also have an analogous result on . Theorem 2 Let , q be a prime number and be defined in (4). Then we have . Moreover, we also have where the upper bound holds just for , and where the lower bound holds just for .
The paper is organised as follows: in Section 2 we will see how to compute using the values of Euler's Γ function and the Fast Fourier Transform algorithm; we will also describe the actual computation we performed and how Theorem 1, Theorem 2 are obtained. In Sections 3-4 we will see how to efficiently evaluate , and , for using precomputed values of the Riemann zeta-function at positive integers. We will also insert some tables and figures (the scatter plots were obtained using GNUPLOT, v.5.2.8).
Section snippets
Acknowledgments
I wish to thank Luca Righi (University of Padova) for his help in developing the C language implementation of the algorithms described in Sections 3-4. I would also like to thank the referee for his/her remarks.
Computation of and proofs of Theorems 1-2
Recall that q is an odd prime and let χ be a primitive non-principal Dirichlet character mod q. The values of can be computed in two different ways. Recalling eq. (3.1) of [2], we have , so that where is the digamma function and Γ is Euler's function. As we will see later, for computational purposes it is in fact more efficient to distinguish between the parity of the Dirichlet characters. If χ is an even character
An alternative algorithm to compute ,
We describe here an alternative way of computing , , which is based on the well-known Euler formula (see, e.g., Lagarias [7, eq. (2.3.3)]): where is the Riemann zeta-function.
We follow the argument used in Languasco-Righi [11] to study the Ramanujan-Deninger Gamma function . We immediately remark that the series in (12) absolutely converges for ; this fact and the well-known relation let us obtain
An alternative algorithm to compute ,
Here we apply to the digamma function , , the same argument used in Section 3. The starting point is the well-known Euler formula (see, e.g., Lagarias [7, Theorem 3.1.2]):
We immediately remark that the series in (23) absolutely converges for ; this fact and the well-known relation let us obtain , , in two different ways. Recalling and , we also remark that, letting , , for every
References (17)
Number Theory, vol. II: Analytic and Modern Tools
(2007)- et al.
Values of the Euler ϕ-function not divisible by a given odd prime, and the distribution of Euler-Kronecker constants for cyclotomic fields
Math. Comput.
(2014) - et al.
The design and implementation of FFTW3
Proc. IEEE
(2005) The accuracy of floating point summation
SIAM J. Sci. Comput.
(1993)- et al.
Numerical results on class groups of imaginary quadratic fields
Further remarks on reducing truncation errors
Commun. ACM
(1965)Euler's constant: Euler's work and modern developments
Bull. Am. Math. Soc.
(2013)- et al.
Conditional bounds for the least quadratic non-residue and related problems
Math. Comput.
(2015)Math. Comput.
(2017)
Cited by (10)
The Kummer ratio of the relative class number for prime cyclotomic fields
2024, Journal of Mathematical Analysis and ApplicationsNumerical estimates on the Landau-Siegel zero and other related quantities
2023, Journal of Number TheoryA unified strategy to compute some special functions of number-theoretic interest
2023, Journal of Number TheoryUniform effective estimates for |L(1,χ)|
2022, Journal of Number Theory