Arxiv:1702.05635V2 [Math.CA] 20 May 2020 Napprwitnb H Elkong .Hry[0,A Neetn in Interesting an [10], Hardy H
Total Page:16
File Type:pdf, Size:1020Kb
ON SPECIAL RIEMANN XI FUNCTION FORMULAE OF HARDY INVOLVING THE DIGAMMA FUNCTION ALEXANDER E PATKOWSKI Abstract. We consider some properties of integrals considered by Hardy and Koshliakov, that have have connections to the digamma function. We establish a new general integral formula that provides a connection to the polygamma function. We also obtain lower and upper bounds for Hardy’s integral through properties of the digamma function. Keywords: Fourier Integrals; Riemann xi function; Digamma function. 2010 Mathematics Subject Classification 11M06, 33C15. 1. Introduction and Main formulas In a paper written by the well-known G. H. Hardy [10], an interesting integral formula is presented (corrected in [3]) (1.1) ∞ ∞ Ξ(t/2) cos(xt) 1 1 1 1 2 4x dt = e−x 2x + γ + log(π) + log(2) + ex ψ(t+1)e−πt e dt, 1+ t2 cosh(πt/2) 4 2 2 2 Z0 Z0 ∂ where ψ(x) := ∂x log(Γ(x)), Γ(x) being the gamma function [1, pg.1], and Ξ(t) := 1 ξ( 2 + it), where [11, 15] 1 s s ξ(s) := s(s 1)π− 2 Γ( )ζ(s). 2 − 2 arXiv:1702.05635v2 [math.CA] 20 May 2020 Here we have used the standard notation for the Riemann zeta function ζ(s) := n−s, for (s) > 1. Koshliakov [12, eq.(14), eq.(20)] (or [3, eq.(1.15)]) pro- n≥1 ℜ Pduced this formula as well, but in a slightly different form, ∞ ∞ Ξ(t/2) cos(xt) 2 4x (1.2) 2 dt = ex (ψ(t + 1) log(t))e−πt e dt. 1+ t2 cosh(πt/2) − Z0 Z0 Here the trick, which has been exploited in the many studies [3, 4, 5, 6, 7, 12, 14], is to re-write the left side of (1.1) as an inverse Mellin transform by utilizing the classical functional equation ξ(s)= ξ(1 s) [15]. Titchmarsh has used simpler − but similar integrals than the left side of (1.1) to obtain Hardy’s result that Ξ(t) 1 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 <c< 1, 1 c+i∞ πζ(1 s) (1.3) ψ(x + 1) log(x)= − x−sds. − −2πi sin(πs) Zc−i∞ Differentiating (1.3) m times, gives us 1 c+i∞ πζ(1 s) Γ(s + m) (1.4) ψ¯ (x) = ( 1)m−1 − x−s−mds. m − 2πi sin(πs) Γ(s) Zc−i∞ This may also be obtained from the integral formula ∞ 1 1 ψ¯ (x) = ( 1)m tme−xtdt, m − t − et 1 Z0 − valid for (s) > m. Equation (1.4) gives us the Mellin transform for our g(t). ℜ − Now proceeding with our choices for f(t) and g(t) in the beginning of the proof, and noting sin(( 1 + it)π) = cosh(πt), the functional equation ξ(s)= ξ(1 s) with 2 − (1.4), the Mellin transform of f(t), π−s/2x−sΓ(s/2)/2, we have ∞ ′+ ∞ 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 <c < 1. Put c = 2 and we find the right side is equal to 1 ∞ Ξ(t) Γ( 1 it + m) 2 − x−itdt. 2√x (t2 + 1 )cosh πt Γ( 1 it) Z−∞ 4 2 − Splitting this bilateral integral up into two integrals and replacing x by e2x leads to the result. 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 <b< 1, we find ℜ b+i∞ 1 2 s x −s 2π (1.6) Γ ( )( ) ζ(s)ds + =4 K0(nπx). 2πi 2 2 x b−i∞ ≥1 Z nX Now noting (1.3) and (1.6) in conjunction with Parseval’s theorem for Mellin trans- forms, we obtain the result after applying the functional equation for ξ(s). 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.