arXiv:1702.05635v2 [math.CA] 20 May 2020 napprwitnb h elkonG .Hry[0,a neetn in interesting an [10], Hardy H. (1.1) [3]) G. in well-known (corrected presented the is by formula written paper a In u iia nerl hntelf ieo 11 ooti ad’ result Hardy’s obtain to (1.1) of side left the than integrals similar but 00MteaisSbetCasfiain1M6 33C15. 11M06, Classification Subject Mathematics 2010 Keywords: 12 2 form, different (1.2) slightly a in but well, as formula this duced 4,i or-rt h etsd f(.)a nivreMli rnfr b transform 6 Mellin inverse 5, an 4, as [3, (1.1) equation of studies functional side many classical left the the the in re-write exploited to been is 14], has which trick, the Here P where ξ eew aeue h tnadntto o h imn eafun zeta Riemann the for notation standard the used have we Here Z ( 0 2 1 n ∞ ≥ + 1 NSEILREANX UCINFRUA OF FORMULAE FUNCTION XI RIEMANN SPECIAL ON Ξ( + 1 it n ψ t/ − ) ( rpriso h iam function. i digamma Hardy’s the for of bounds properties upper and to lower connection obtain a also We provides function. that formula integral general new funct a digamma the to connections have have that Koshliakov, Abstract. , x AD NOVN H IAM FUNCTION DIGAMMA THE INVOLVING HARDY s t 2) hr 1,15] [11, where , 2 := ) ore nerl;Reanx ucin iam function. Digamma function; xi Riemann Integrals; Fourier for cosh( Z 0 cos( ∞ ∂x ℜ ∂ ( Ξ( πt/ + 1 s xt log(Γ( ecnie oepoete fitgascniee yHard by considered integrals of properties some consider We ) t/ ) 1. 2) > t 2) 2 dt nrdcinadMi formulas Main and Introduction 1 x cosh( . )) ξ = cos( ohikv[2 q(4,e.2) o 3 q(.5] pro- eq.(1.15)]) [3, (or eq.(20)] eq.(14), [12, Koshliakov ( , LXNE PATKOWSKI E ALEXANDER s 4 1 := ) Γ( πt/ e xt − x ) en h am ucin[,p.] n Ξ( and pg.1], [1, function gamma the being ) x ξ 2) 1 2 ( s s dt 2 = ) ( x s = + − ξ e 1 1) (1 1 2 x γ π Z − 0 − + ∞ s 2 s 1] icmrhhsue simpler used has Titchmarsh [15]. ) 2 1 Γ( ( ψ log( ( 2 s t ) 1) + ζ π ( log(2) + ) s ) . − log( o.W establish We ion. h polygamma the tga through ntegral t )) + e − 1 2 πt e ction and y 2 x e Z 4 utilizing y x 0 htΞ( that dt. ∞ ,12, 7, , ζ ψ ( tegral s ( t := ) t := ) +1) t ) e − πt 2 e 4 x dt, 2 ALEXANDER E PATKOWSKI has infinitely many real zeros [15]. See Csordas [2] for some more recent work in this direction. A recent paper by Dixit offers a beautiful generalization of formula (1.1), which involves a confluent hypergeometric function [3, Theorem 1.3]. Other refinements can be found in [4, 5]. The purpose of this paper is to offer some further results concerning (1.1) that appear to have been overlooked. A new generalization is offered as well that takes a different route than the ones offered by [3, 4, 5], considering what is essentially the polygamma function [9, pg.904, eq. NH 37(1)]. In the following section, we offer some new inequalities that we noticed upon examining the integral given in (1.1). Recall the Mellin transform of a suitable function f(t) is given by [11, pg.90, eq.(4.105)] ∞ f(t)ts−1dt = F (s), Z0 and its inverse, 1 c+i∞ f(t)= F (s)t−sds, 2πi Zc−i∞ provided c is a real number chosen where F (s) is analytic. We will be applying Parseval’s formula for Mellin transforms throughout the proofs of our theorems [15, pg.34, eq.(2.15.10)]
∞ 1 c+i∞ f(t)g(t)dt = F (s)G(1 s)ds. 2πi − Z0 Zc−i∞ Here the needed condition for validity is to have c chosen as a real number restricted to the region where F (s)G(1 s) is analytic. −
Theorem 1.1. For m N0, we have that ∈ ∞ 2 4x ( 1)mex tmψ¯ (t)e−πt e dt − m Z0 ∞ 1 1 1 Ξ(t) Γ( 2 it + m) −2xit Γ( 2 + it + m) 2ixt = 2 1 −1 e + 1 e dt, 2 0 cosh(πt)(t + ) Γ( it) Γ( + it) Z 4 2 − 2 ∂m ψ¯ (t) := (ψ(t + 1) log(t)). m ∂tm −
Proof. This is a special application of Parseval’s theorem for Mellin transforms 2 −πt m with one function chosen as f(t) = e , and the other as g(t) = t ψ¯m(t). From ON SPECIAL RIEMANN XI FUNCTION FORMULAE OF HARDY INVOLVING THE DIGAMMA FUNCTION3
Titchmarsh [15, eq.(2.15.7)] we have for 0
∞ ′+ ∞ 2 4 c i m x m −πt e x 1 πξ(s) Γ(1 s + m) −s ( 1) e t ψ¯m(t)e dt = − x ds, − 2πi ′ sin(πs) s(s 1)Γ(1 s) Z0 Zc −i∞ − − ′ ′ 1 for 0 Koshliakov [12, eq.(36), eq.(40)] gives a similar integral to (1.1) but squaring the 2 1 −1 integrand Ξ(t)(t + 4 ) . ∞ 2 ∞ Ξ (t/2) cos(xt) x/2 x dt = e K0(2πe t)Λ(t)dt, (1 + t2)2 cosh(πt/2) Z0 Z0 where 2 1 2 1 1 Λ(t)= ζ(2) + γ 2γ1 +2γ log(t)+ log (t)+ d(n) . − 2 t + n − n ≥1 nX Here d(n) denotes the number of divisors of n, and the Stieltjes constant is log(n) log(N)2 γ1 = lim . N→∞ n − 2 1≤ ≤ Xn N Moll and Dixit state his result in [6, Theorem 4.4], and offer further interesting generalizations. However, we noticed a slightly different form of this integral based on the type of proof we used in the previous theorem. 4 ALEXANDER E PATKOWSKI Theorem 1.2. Let Ks(t) denote the modified Bessel function. For x> 0, (1.5) ∞ x/2 ∞ 2 2πe −x/2 −x Ξ (t) cos(xt) (ψ(t+1) log(t)) e 4 K0(2πnte ) dt = dt. − t − 2 1 2 cosh(πt) 0 ≥1 0 (t + 4 ) Z nX Z −1/2 −t Proof. First, since [9, pg.920, eq. WA 231, 245(9)] K0(t)= O(t e ), for t → ∞ in arg(t) < 3π/2 where as usual K (t) is the modified Bessel function [9, pg.928], | | s it follows that the series on the right side of the formula (by [9, pg. 676, eq. EH II 51(27)]) d> 0, d+i∞ 1 2 s x −s Γ ( )( ) ds =4 K0(nπx), 2πi 2 2n ≥1 d−i∞ ≥1 nX Z nX converges absolutely. Summing through, we find for d′ > 1, ′ d +i∞ 1 2 s x −s Γ ( )( ) ζ(s)ds =4 K0(nπx), 2πi ′ 2 2 d −i∞ ≥1 Z nX 2 s x −s 2π which has a simple pole at s = 1. Note that lim →1(s 1)Γ ( )( ) ζ(s) = , s − 2 2 x 1 since Γ( )= √π and lim →1(s 1)ζ(s)=1. If we calculate this residue and move 2 s − the line of integration to (s)= b, 0 Other associated integral formulae related to Fourier cosine transforms may be obtained from using known evaluations. Theorem 1.3. We have for π > (β′) , 2 |ℑ | ∞ cosh(t) 2 −2t (1.7) et/2 2e−t/2 e−πn e dt cosh(2t) + cosh(2β′) − 0 ≥1 Z nX ′ β ∞ ′ πe 2 4β = (ψ(t + 1) log(t))e−πt e dt. 4 cosh(β′) − Z0 Proof. We apply the classical [15, eq.(2.16.2)] ∞ Ξ(t) π 2 −2x cos(xt)dt = (ex/2 2e−x/2 e−πn e ), (t2 + 1 ) 2 − 0 4 ≥1 Z nX ON SPECIAL RIEMANN XI FUNCTION FORMULAE OF HARDY INVOLVING THE DIGAMMA FUNCTION5 with an integral from [9, pg.511, eq. ET I 31(16)], valid for π (α) > (¯αβ) , and ℜ |ℑ | ∞ cos(ty) cosh(tα/2)dt π cos(βy/α) = . cosh(αt) + cosh(β) 2α cosh(β/2) cosh(πy/α) Z0 We need only let β 2β′, and α =2, and then we have π > (β′) . This coupled → 2 |ℑ | with Parseval’s theorem gives the result. The integral formula in Theorem 1.3 may be useful for fast numerical approxima- tions, since the integrand of the integral on the left side is comprised entirely of exponentials. A possible avenue would be to truncate the series at the first term, since it is a double exponential. Note that the integral on the left side of (1.7) hints to us about the non-negativity of the integral on the right side. We will investigate this further in the next section. 2. Inequalities We first start with some preliminaries about some special functions which can be found in [1]. First we recall the error function [9, pg.887], x 2 2 erf(x)= e−t dt, √π Z0 and its complement, erfc(x)=1 erf(x). Also, erfi(x) := i erf(ix). We also use − − the exponential integral [9, pg.883] ∞ e−t Ei(x)= dt. − t Z−x Theorem 2.1. ∞ Define the Hardy integral for y > 0 as I(y) := 0 (ψ(t + 1) 2 − log(t))e−yt dt. Then we have, R 2 −y (2c1) √y Ei( y) e √yπ − + erfc(√y) I(y). √π erfi(√y) − 4 − 12 12 ≤ Furthermore, π Ei( y) I(y) c2 erf( 2y) − , ≤ s 2y − 4 r p where c1 =0.952894, and c2 =1.56624. Proof. We first need a result from the paper [13, Corollary 1], for x> 0, 1 1 1 (2.1) < ψ(x + 1) log(x) < . 2x − 12x2 − 2x 6 ALEXANDER E PATKOWSKI For the upper bound we first write, 1 ∞ 2 2 I(y)= (ψ(t + 1) log(t))e−yt dt + (ψ(t + 1) log(t))e−yt dt. − − Z0 Z1 Applying (2.1) we have that this is 1 2 Ei( y) (ψ(t + 1) log(t))e−yt dt − . ≤ − − 4 Z0 By Schwarz’s inequality and the fact that 1 (ψ(t + 1) log(t))2dt =1.56624 − Z0 we have the upper bound. For the lower bound, we write the Schwarz’s inequality in the form 1 2 1 1 2 2 (ψ(t + 1) log(t))1/2dt (ψ(t + 1) log(t))e−yt dt eyt dt , − ≤ − Z0 Z0 Z0 and also note, 1 (ψ(t + 1) log(t))dt =0.952894. 0 − Z p This calculation, coupled with similar computations derived from the left side of (2.1) and the definitions of the special integral functions, gives our result. Now it can be seen that the lower bound actually implies the Hardy integral I(x) is non-negative for x> 0, and also the squeeze theorem shows that I(x) 0 when → x . → ∞ 3. Other observations and further comments While studies such as [6] show there is a plethora of integral formulae like Hardy’s (1.1) to be found, it would be interesting to see if these integrals could be incor- porated into similar work like Csordas [2]. We were interested in seeing if (1.1) could be adapted to imply information beyond Hardy’s theorem, that Ξ(t) has infinitely many real zeros. The classical ideas are included in Titchmarsh’s text [15, pg. 256-260], and are based around the idea of first assuming Ξ(t) is of one sign, and then arriving at a contradiction from properties of an integral of the form ∞ 0 Ξ(t)k(t)dt, k(t) non-negative for t> 0. The idea in [15, pg. 260] considers the R ON SPECIAL RIEMANN XI FUNCTION FORMULAE OF HARDY INVOLVING THE DIGAMMA FUNCTION7 moments k(t)= t2n. Apply the operator ∂¯ = ∂2 + 1 of Glasser [8] to (1.1) (after x − x 4 replacing t by 2t) and set x = x′ + iπ to get ∞ ∞ ′ x/2 2x ′ sin(x t) sinh(πt) lim ∂¯x(e I(πe )) = Ξ(t)cos(x t)dt + i Ξ(t) dt. x→x′+iπ cosh(πt) Z0 Z0 Differentiating 2n times and setting x′ = 0 gives us the right side as the value ∞ ( 1)n t2nΞ(t)dt. Note that the integral I(x) clearly converges only for (x) > 0, − 0 ℜ and henceR I(πe2z) converges when π < (z) < π . It would be interesting if new − 4 ℑ 4 information could be obtained about Ξ(t) from alternative forms of these moments like those found here. References [1] G. Andrews, R. Askey, and R. Roy. Special Functions, volume 71 of Encyclopedia of Mathe- matics and its Applications. Cambridge University Press, New York, 1999. [2] G. Csordas, Fourier Transforms of Positive Definite Kernels and the Riemann ξ-function, Computational Methods and Function Theory, Volume 15, Issue 3, pp 373–391 (2015). [3] A. Dixit. Analogues of the general theta transformation formula, Proc. Roy. Soc. Edinburgh, Sect. A, 143:371–399, 2013. [4] A. Dixit, N. Robles, A. Roy and A. Zaharescu, Zeros of combinations of the Riemann ξ- function on bounded vertical shifts, J. Number Theory Volume 149, April 2015, Pages 404?434 [5] A. Dixit, A. Roy and A. Zaharescu, Riesz-type criteria and theta transformation analogues, J. Number Theory 160, p. 385–408 (2016). [6] A. Dixit, V. H. Moll Self-reciprocal functions, powers of the Riemann zeta function and modular-type transformations, J. Number Theory 147 (2015), 211–249. [7] W. L. Ferrar, Some solutions of the equation F (t) = F (t−1), J. London Math. Soc. 11 (1936), 99–103. [8] M.L. Glasser A note on the Fourier transform, Journal of Physics A: Mathematical and Theoretical, Volume 49, Number 45 (2016), 454004. [9] I. S. Gradshteyn and I. M. Ryzhik. Table of Integrals, Series, and Products. Edited by A.Jeffrey and D. Zwillinger. Academic Press, New York, 7th edition, 2007. [10] G.H. Hardy, Note by Mr. G.H. Hardy on the preceding paper, Quart. J. Math. 46 (1915), 260–261. [11] H. Iwaniec and E. Kowalski, Analytic number theory, American Mathematical Society Col- loquium Publications, vol. 53, American Mathematical Society, Providence, RI, 2004. [12] N. S. Koshlyakov, On a transformation of definite integrals and its application to the theory of Riemann’s function ζ(s), Comp. Rend. (Doklady) Acad. Sci. URSS 15 (1937), 3–8. [13] Feng Qi, Run-Qing Cui, Chao-Ping Chen, Bai-Ni Guo, Some completely monotonic functions involving polygamma functions and an application, J. Math. Anal. Appl. 310 (2005) 303–308. 8 ALEXANDER E PATKOWSKI [14] S. Ramanujan, New expressions for Riemann’s functions ξ(s) and Ξ(t), Quart. J. Math., 46: 253–260, 1915. [15] E. C. Titchmarsh, The theory of the Riemann zeta function, Oxford University Press, 2nd edition, 1986. 1390 Bumps River Rd. Centerville, MA 02632 USA E-mail: [email protected]