arXiv:2009.14070v3 [math.NT] 15 Feb 2021 ei ucin . functions metic yFet[6,Sgl[6 n eag 1] h rgnmti eisi qu in series trigonometric The (1.1) [14]. and Delange [17], ot and Littlewood the [26] and on Segal and Hardy [16], Beurling, by Flett an and by studied functional Nymann a first to hand series due one trigonometric hypothesis the Riemann On the function: to zeta Riemann’s of facets tdffr rmtefiiesum finite the from differs It ad n iteodpoe 1]ta,as that, [17] proved Littlewood and Hardy h anproeo hswr st rn olgtanwrltosi betw relationship new a light to bring to is work this of purpose main The e od n phrases. and words Key 2010 RTMTCADAAYI FTESERIES THE OF ANALYSIS AND ARITHMETIC ahmtc ujc Classification. Subject Mathematics iteod(e 34) ehglgtsm fiscrosaayia and analytical curious its of some properties. highlight We (3.4)). (see Littlewood n nteeuvlnebtenteReanhptei n h dens the and hypothesis Riemann in the space between functional equivalence the on ing Abstract. nti ae ecnetaclbae hoe fNmnadBeurl- and Nyman of theorem celebrated a connect we paper this In L ad-iteodfnto,Fae nerl erigsterm A theorem, Beurling’s integral, Franel function, Hardy-Littlewood X n> HE EBRADRGRGAY ROGER AND SEBBAR AHMED 2 f (0 ( 1 x , )t rgnmti eiscniee rtb ad and Hardy by first considered series trigonometric a to 1) n 1 = ) X n sin ≤ x O f 1. ( n 1 n x x introduction = ) sin = (log otemmr forfin alsBerenstein Carlos friend our of memory the To 13,1M8 17,11K65. 11K70, 11M38, 11M32, O n x X n x ∞ =1 as ,
1 ) X n>x 3 4 x n 1 lglog (log x sin ed to tends n ed to tends x 2 n x ! . x = ) 3 4 O ∞ + (1) ∞ ǫ , by , . arithmetical t fsome of ity X n lssapproach alysis ∞ =1 hnfollowed then n 1 e ada hand her sinis estion sin e two een n x rith- 2 AHMED SEBBAR AND ROGER GAY and that 3 f(x) = Ω (log log x) 4 from the fact that for x 5, the number of n x whose prime divisors are equivalent x≥ ≤ to 1 modulo 4 is C 1 , where C is a constant. Delange [14] showed that f(x) (log x) 2 is not bounded on the real line only from the following result on the reciprocals of primes in arithmetic progressions 1 = p ∞ p prime, p 1(modX 4) ≡ and obtained the Ω-result of Hardy and Littlewood just because 1 1 = log log x + c + o(1). p 2 p prime x, p 1(modX≤ 4) ≡ This trigonometric series, despite its simplicity, has many similarities with the Riemann zeta function [16] and deep relation to the divisor functions through the sawtooth function t t 1 if t = t 1 ∞ sin 2mπt 2 (1.2) t = = − ⌊ ⌋ − 6 ⌊ ⌋ { } −π m m=1 0 if t = t . X ⌊ ⌋ For s C we define ∈ s s s σ (n)= d , σs(n)= d− d n d n | | s s X X so that n σs(n)= σ (n). For example if we define 1 1 S1(x)= σ1(n), S (x)= σ (n) n x n x X≤ X≤ and 1 x 1 x x 1 ρ(x)= = n{n} n n − n − 2 n x n x X≤ X≤ j k then the divisors and the fractional parts functions are related by 1 x 1 S (x)= = x ρ(x) 1 n n n2 − n x n x X≤ j k X≤ π2 1 = x log x ρ(x)+ O(1). 2 − 2 − Similarly [33] (p.70): π2 S1(x)= x2 xρ(x)+ O(x). 12 − ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 3 X n n We will see ((3.5) with f(2πx) = sin x) an integral representationn=1 of the partial sums of f(x), using the sawtooth function.
2. Nyman-Beurling criterion for the Riemann hypothesis
2.0.1. Nyman-Beurling theorem. For x > 0, let ρ(x) be the fractional part of x so θ that ρ(x) = x x . To each 0 < θ 1 we associate the function ρ (x) = ρ( ). − ⌊ ⌋ ≤ θ x θ Then 0 ρ (x) 1 and ρ (x) = if θ < x. We introduce, as in [7], [15], [4], [5], ≤ θ ≤ θ x [6], [23], [30] and the more recent book [24] N θ N = f, f(x)= a ρ( n ), a R, θ (0, 1], a θ =0, N 1 . M n x n ∈ n ∈ n n ≥ ( n=1 n=1 ) X X Each function in has at most a countable set of points of discontinuity, and is identically zero forMx> 0. Theorem 2.1 (Nyman-Beurling). Let 1 < p . The subspace is dense in the Banach space Lp(0, 1) if and only if the Riemann≤∞ zeta function ζM(s) has no zero in 1 the right half plane Res> . p The fundamental relations in the proof of this theorem are 1 θ s 1 θ s ζ(s) (2.1) ρ( )x − dx = θ , Res> 1, x −1 s − s Z0 − which is just a variant of the classical representation
s ∞ u u (2.2) ζ(s)= s − ⌊ ⌋ du. s 1 − (u + 1)s+1 − Z0 It follows from (2.1) that for f(x) ∈M 1 N s 1 ζ(s) s f(x)x − dx = akθk. 0 − s Z Xk=1 The study of the function f(x) is intimately linked to that of following function 1 t t 2 if t = t t = − ⌊ ⌋ − 6 ⌊ ⌋ { } 0 if t = t . ⌊ ⌋ We have the formal Fourier series expansion [12], [13] ∞ a 1 ∞ A (2.3) n nθ = n sin 2πnθ n { } −π n { } n=1 n=1 X X 4 AHMED SEBBAR AND ROGER GAY where
An = ad. d n X| Davenport considered the cases of 2 a = µ(n); a = λ(n); a = Λ(n); a 2 = µ(n), a =0, n = m . n n n n n 6 These arithmetical functions have their usual number-theoretic meanings. For exam- ple if ω(n) is the number of distinct prime factors of n or, in other terms, ω(n)= 1 b n X| and ω(1) = 0, then the M¨obius function µ(n) is defined by 0 if n is divisible by a perfect square > 1 µ(n)= ( 1)ω(n) otherwise − and the Von Mangoldt function Λ(n) is defined by log p if n = pα for a prime p and some α N Λ(n)= ∈ 0 otherwise In the case of the M¨obius function an = µ(n), Davenport uses Vinogradov’s method, a refinement of Weyl’s method on estimating trigonometric sums, to prove that for any fixed h
2iπnx h (2.4) µ(n)e = O y(log y)− n y X≤ uniformly in x R/Z. The implied constants are not effective. There have been ∈ several results justifying (2.3) for other particular sequences (an). The most general problem is considered in [18]. It should be noted that the Davenport or Hardy Littlewood estimates admit a common analysis. For the convenience of the reader we gather together a few classical results on exponential sums. Let I be an interval of length at most N 1 and let f : I R be a smooth function satisfying the estimates x I, 2 N≥ T, j 1 → ∈ ≤ ≪ ≥ T f (j)(x) = exp O(j2) | | N j then with f(x)= ex (1) Van der Corput estimate: For any natural number k 2, we have ≥ 1 1 T 2k−2 1 (2.5) e(f(n)) = O log 2 (2 + T ) N N k n I ! X∈ ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 5 X n n (2) Vinogradov estimate: For some absolute constant c>n=1 0.
1 c (2.6) e(f(n)) N − k2 . N ≪ n I X∈ 2.0.2. The functions x , ρ(x) and x . The Hardy-Littlewood-Flett function f(x) is related, in many ways,⌊ ⌋ to the three{ functions} x , ρ(x) and x . The floor function x is related to the divisor function d(n)=1⌊ ⌋⋆ 1(n), the{ multiplicative} square ⌊convolution⌋ product of the constant function 1, through the Dirichlet hyperbola method. More generally if g, h are two multiplicative functions and f = g ⋆ h. The Dirichlet hyperbola method is just the evaluation of a sum in two different ways:
f(n) = g(a)h(b)= g(a)h(b)+ g(a)h(b) n x n x ab=n a √x b x b √x a x X≤ X≤ X X≤ X≤ a X≤ X≤ a g(a)h(b). − a √x b √x X≤ X≤ If g = h, then 2 f(n)=2 g(a)h(b) g(a) . − n x a √x b x a √x X≤ X≤ X≤ a X≤ As an application we have the estimate [29] (p. 262) for the divisor function d =1⋆1:
1 d(x)= x log x + (2γ 1)x + O(x 2 ). − The importance of the functions x and ρ(x) lies in the integral representations of the Riemann zeta-function: { } 1 ∞ x ∞ ρ(x) ζ(s)= s { } − 2 dx = s dx − xs+1 − xs+1 Z0 Z0 1 valid for 1 < Res< 0. Making the change of variable x = and applying Mellin − u inversion formula gives c+i 1 1 ∞ ζ(s) s ρ( )= u− ds. u −2iπ c i s Z − ∞ For later use, we give some details on the case considered by Davenport in (2.3). From (2.4) we obtain for 1 By the functional equation of the Riemann ζ-function and the functional equation of the Γ-function we obtain for 0 ∞ a+i µ(n) 1 ∞ 1 s 1 ρ(nx)= Γ(s) sin( πs)(2πx)− ds = sin(2πx). n −2iπ2 2 −π n=1 a i X Z − ∞ Using the classical equivalent formulation of the Prime Number Theorem that ∞ µ(n) = 0 we obtain Davenport’s relation n n=1 X ∞ µ(n) 1 (2.7) nx = sin(2πx) n { } −π n=1 X where the convergence is uniform by Davenport estimate (2.4). We will need two important properties of the function x : Kubert identity: { } l (2.8) x + = mx m { } l modXm Franel formula: 1 lcm(a, b) (2.9) ax bx dx = . { }{ } 12ab Z0 Kubert identity and Franel’s formula are interesting features shared by many func- tions. Let Br(x) be the Bernoulli polynomial defined by tetx ∞ = B (x)tr, t < 2π, et 1 r | | r=0 − X so that 1 1 B (x)= x , 2!B (x)= x2 x + , 1 − 2 2 − 6 ··· If r 2 is even then for 0 x 1 ≥ ≤ ≤ ∞ r 1 1 2 cos(2lπx) (2π) B (x)=( 1) 2 − r − lr Xl=1 with absolute convergence of the series. The Hurwitz zeta function ζ(s, x) is defined for Res> 1 by ∞ 1 ζ(s, x)= . (x + n)s n=0 X Then Br(x) and ζ(s, x) both satisfy the functional equation [22] 1 k 1 f(x)+ f(x + )+ + f(x + − )= f (k)f(kx), k ··· k ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 7 X n n (k) 1 n (k) s n=1 where f = k − if f(x) = Bn(x) and f = k if f(x) = ζ(s, x). Furthermore, if a, b denote arbitrary positive integers and (a, b) = gcd(a, b), [a, b] = lcm(a, b) the greatest common divisor and least common multiple respectively of a and b, then [22]: 1 r r 1 B2r (a, b) B (ax)B (bx) dx =( 1) − r r − (2r)! [a, b] Z0 1 and for Res> 2 1 2Γ2(s)ζ(2s) (a, b) s ζ(1 s, ax)ζ(1 s, bx) dx = . − − (2π)2s [a, b] Z0 Similarly to (2.2) we have 1 1 ∞ u u ζ(s,w)= + s − ⌊ ⌋ du, (s 1)ws 1 ws − (u + w)s+1 − − Z0 1 and the function ζ(s,w) s 1 is analytic in Res> 0 . In the next section − (s 1)w − { } we use two summation formulas.− If F is an antiderivative of f, then, formally [1] 1 θ 1 f(t) ∞ θ θ (2.10) ρ( )f(t)dt = θ dt n F ( ) F ( ) t t − n − n +1 0 0 n=1 Z Z X and if µ is the M¨obius function and if 0 <θ, x 1 , we have, pointwise [2] ≤ ∞ θ 1 θ µ(n) ρ ρ = χ (x). nx − n x − ]0,θ] n=1 X 3. From Beurling’s theorem to Hardy-Littlewood-Flett function f(x) 3.1. The emergence of Franel integral type. To show that the constant func- tion 1 one has, as in [2], to minimize the norms in L2([0, 1]) ∈ M N α (3.1) 1+ a ρ j k j k j=1 X · which brings back to the evaluation of integrals of Franel type, computed in [2]: 1 1 β J(β)= ρ ρ dx, β [0, 1]. x x ∈ Z0 8 AHMED SEBBAR AND ROGER GAY To show that the function sin x one has, this time, to minimize the norms ∈ M N α (3.2) sin x + a ρ j k j k j=1 X · Using (2.7) the minimization problem reduces to evaluation of the scalar products in L2 (0, 1) giving the Fourier sine series of the function θ , that is { x } θ 1 θ µ 1 θ jnx a =( √2 sin(nπx)) = √2 sin(nπx)dx = π√2 j dx n {x} {x} − j {x}{ 2 } Z0 j 1 Z0 X≥ 1 a and then to the evaluation of bx dx, another kind of integrals of Franel 0 {x}{ } type. We compute these integralsZ in the case a = m, b = n, m and n being integers. 3.2. The second kind of Franel type integrals I = 1 nx m dx, n, m n,m 0 { }{ x } ∈ N∗. The values of the integrals In,m are given by the following R Theorem 3.1. For positive integers m, n, the modified Franel integrals are given by n I = + m log m + m(n 1) log(mn) m(log((n 1)!)) n,m m − − − n(n 1) nm2 m mn − (ζ(2) (1 )) + (1 ). − 2 − 2 − − j − jk 1 j m 1 k n 1,mn jk ≤X≤ ≤ ≤ X− ≥ Let us first give few examples: 5 25 3 I = log(2) ζ(2); I = + log(2) 2 log(3) ζ(2) (2,1) 2 − − (3,1) 6 − − 2 35 35 I = 5 log(2) + log(3) 2ζ(2); I = 5 log(2) + log(3) 2ζ(2) (4,1) 6 − − (5,1) 6 − − 7 61 9 I = 2ζ(2); I = ζ(2) (1,2) 2 − (1,3) 8 − 2 5989 25 49 I = ζ(2); I = 2 log(2) 4ζ(2) (1,4) 288 − 2 (2,2) 6 − − 171 18469 I = 3 log(2) 9ζ(2); I = 4 log(2) 16ζ(2) (2,3) 10 − − (2,4) 630 − − 15059 196 I = 5 log(2) 25ζ(2); I = + 2 log(2) 4 log(3) 6ζ(2) (2,5) 336 − − (3,2) 15 − − π2 We observe that in all these examples the factor ζ(2) = is present. 6 For the proof we consider the two functions defined on ]0, + [ ∞ f(x)= fn(x)= xχ[0,1](x) nx , g(x)= x χ[1,+ ](x) { } { } ∞ ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 9 X n n and their multiplicative convolution n=1 + ∞ a dx (f⋆g)(a)= f(x)g( ) , (f⋆g)(m)= I . x x n,m Z0 We split the computations in several steps. A natural method is to use first the Mellin transform with its property (f⋆g)(s)= (f)(s) (g)(s), followed by an M M M 1 f(x) inversion. The main idea is the decomposition formula (2.10), valid if | | dx 0 x is finite: Z 1 θ 1 ∞ n θ θ ρθ(x)f(x)dx = ( n)f(x)dx + f(x)dx, θ x − x 0 n=1 n+1 θ Z X Z Z or in an generalized function form, ∞ θ θ ρθ(x)= ( n)χ[ θ , θ ](x)+ χ[θ,1] x − n+1 n x n=1 X where χB denotes the characteristic function of the set B. 3.2.1. Computations of different integrals. (1) Computation of F (s)= M(f)(s) For σ = Re s> 2 we have − 1 F (s)= nx xsdx 0 { } Z k+1 n = (nx k)xsdx k − 0 k n 1 n ≤ ≤ − Z X k+1 1 n = n xs+1dx k xsdx − k 0 1 k n 1 n Z ≤X≤ − Z n 1 = k((k + 1)s+1 ks+1) s +2 − (s + 1)ns+1 − 1 k n 1 ≤X≤ − n 1 = ns+2 (1+2s+1 +3s+1 + ns+1) s +2 − (s + 1)ns+1 { − ··· } 10 AHMED SEBBAR AND ROGER GAY (2) Computation of G(s)= M(g)(s) For 2 < σ = Res< 1 we have − − + ∞ s 1 G(s)= x − dx 1 { } Z k+1 s 1 = (x k)x − dx − k 1 k X≥ Z + k+1 ∞ s s 1 = x dx k x − dx − 1 k 1 k Z ≥ Z + X ∞ 1 = xsdx k((k + 1)s ks) − s − 1 k 1 Z X≥ 1 ζ( s) = − σ < 1 s +1 − s − Hence for 2 c+i s 1 ∞ n.m In,m = ds 2iπ c i (1 s)(2 s) Z − ∞ c+ − s − 1 ∞ m 1 s 1 s 1 s 1 s − − − − 2 1 s (n (1+2 +3 + (n 1) ))ds − 2iπ c i (1 s) n − − ··· − Z − ∞ c+i − 1 ∞ ζ(s) n 1 1 s 1 s 1 s s − − − + ( 1 s (n (1+2 + +(n 1) ))m ds 2iπ c i s 2 s − (1 s)n − − ··· − Z − ∞ − − We treat the last integral by expanding the ζ function in Dirichlet series. We will treat each type of integrals appearing separately. Then we proceed to the necessary groupings in order to conclude. c+i ∞ In the following we write instead of , with 1 By putting together these partial results we end the proof of Theorem (3.1). θ 3.3. Second approach x . The most interesting approach for the evaluation of 1 θ the integral sin(nπt )dt is to use (2.10): t Z0 1 θ θ θ 1 1 sin(nπt)dt = sin(nπt)dt + θ sin(nπt)dt. t t t Z0 Z0 Zθ 12 AHMED SEBBAR AND ROGER GAY Moreover θ θ θ p θ sin(nπt)dt = sin(nπt)( p)dt t θ t − 0 p 1 p+1 Z ≥ Z X θ θ p sin(nπt) p = θ dt p sin(nπt)dt θ t − θ p 1 p+1 p 1 p+1 X≥ Z X≥ Z θ sin(nπt) 1 nπθ nπθ = θ dt + p cos cos t nπ p − p +1 0 p 1 Z X≥ Hence the n-th Fourier coefficient θ 1 θ a =( √2 sin(nπx)) = √2 sin(nπx)dx n {x} {x} Z0 is also 1 sin(nπt) 1 nπθ nπθ a = √2 θ dt + p cos cos . n t nπ p − p +1 0 p 1 ! Z X≥ θν Seeking for the coefficient corresponding to f(x)= 1 ν N cν x the first integral ≤ ≤ { } does not matter since 1 ν N cνθν = 0. It remains to compute ≤ ≤ P P nπθ nπθ A = p cos cos . p − p +1 p 1 X≥ A depends on n and θ. We first consider the finite sum nπθ nπθ A = p cos cos N p − p +1 1 p N ≤X≤ and set x = nπθ. We have by a partial summation x x x A = cos + + cos N cos N 1 ··· N − N +1 x x x = (cos 1) + + (cos 1) + N(1 cos ) 1 − ··· N − − N +1 N x x = 2 sin2 +2N sin2 . − k N +1 Xk=1 Hence ∞ 2 x (3.3) lim AN = 2 sin . N + → ∞ − k Xk=1 ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 13 X n n n=1 We thus obtain one of our main results: the n-th Fourier coefficient an of the funda- mental function θ is related to the value at n of the antiderivative of the function f(x) given in (1.1){ • } 1 ∞ sin(nπt) 1 2 nπθ (3.4) an = √2 θ dt sin , 0 t − nπ k ! Z Xk=1 ∞ u bearing in mind that the derivative of sin2 is f(2u). k Xk=1 To give some useful integral representations we adapt an interesting method, due to Delange [14], and use a result of Saffari and Vaughan [25]. First we introduce for 0 <α 1 ≤ 1 if u u = ρ(u) <α c (u)= − ⌊ ⌋ α 0 otherwise Furthermore for x> 0,y> 1 let 1 1 x ϑ (u)= c ( ). x,y log y n α n n y X≤ According to [25] we have Lemma 3.1. We have the estimate 2 1 ϑx,y(u)= u + O (log x) 3 (log y)− , the O is uniform in u. If f is continuously differentiable function on [0, 1] x 1 1 x f(2π )= 2π f ′(2πu) du = 2π cu( )f ′(2πu) du. n − x − 0 n Z{ n } Z Hence 1 x 1 (3.5) f(2π )= 2π(log x) ϑ (u)f ′(2πu), n n − x,x n x 0 X≤ Z since 1 1 1 ϑ (u)f ′(2πu) du = f ′(2πu) du + (ϑ (u) u) f ′(2πu) du x,x x,x − Z0 Z0 Z0 1 = (ϑ (u) u) f ′(2πu) du. x,x − Z0 14 AHMED SEBBAR AND ROGER GAY From the Lemma (3.1) we get, since f ′ is bounded on (0, 1) 1 x 2 f(2π )= O (log x) 3 . n n n x X≤ A natural example is to consider a Dirichlet character modulo N, χ. In this case χ(n) x 2 sin(2π )= O (log x) 3 . n n n x X≤ We shall not try to give sufficient conditions to justify the process here. The main interest of the remark is that it suggests a method of dealing with various other sums than f(x). 4. Almost periodicity The goal of this section is to show, by elementary methods, that the Hardy- ∞ 1 x Littelwood-Flett function f(x) = sin is not bounded on the real line. First n n n=1 we recall two fundamental resultsX on Bohr-almost periodic functions [9] (p.39, 44, and 58). Theorem 4.1 (The Mean value theorem). For every almost periodic function f(x), there exists a mean value 1 T lim f(x)dx = M f(x) T + T { } → ∞ Z0 and 1 a+T lim f(x)dx = M f(x) T + T { } → ∞ Za uniformly with respect to a. In particular if f is an odd almost periodic function, then its mean M f(x) is zero. { } x Theorem 4.2 (The antiderivative theorem). The integral F (x) = f(t)dt of an 0 almost-periodic function f(x) is almost-periodic if and only if it is bounded.Z Let ∞ 1 x F(x)= cos . n2 n n=1 X The series defining F(x) is uniformly convergent on the real line. The partial sums n 1 x F (x)= cos n p2 p p=1 X ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 15 X n n are almost periodic [9] (Corollary, p.38), and then F(x) isn=1 also almost periodic [9] (Theorem IV, p.38). It is interesting to note that Fn is periodic of period pn = lcm(1, 2, , n) = eψ(n), with ψ(x) is the Chebyshev function, given by ψ(x) = ··· Λ(p), where Λ(n) is the Mangoldt function. p x X≤ n(1+o(1)) The prime number theorem asserts that pn = e as n [27] (p.261). Actually p 3n. → ∞ n ≤ Lemma 4.1. We have 1 ∞ x π lim sin2 = . x + x n 2 → ∞ n=1 X 2x 1 Let x> 0 and n = . The function h : x sin2 , being bounded on [0, π ] x π → x 2 π and continuous on each [α, ], is Riemann-integrable on [0, π ], so by considering 2 2 Riemann sums nx 2 2 1 x π 1 ∞ sin u (4.1) lim sin2 = sin2 dt = du. x + x n t π u2 → ∞ n=1 0 2 X Z Z 2x For x> 0 the function of g(t) = sin2 x is decreasing on ( , + ) and thus t π ∞ ∞ x ∞ x (4.2) sin2 sin2 dt 1. n − 2x t ≤ n=n +1 π Xx Z π 2 ∞ x 2 sin u 2 Since sin dt = x 2 du we deduce the lemma from (4.8) (4.9) and the 2x t 0 u Z π Z relations 2 ∞ sin u ∞ sin u π du = du = . u2 u 2 Z0 Z0 ∞ 1 x Corollary 4.2. The function f(x)= sin is not bounded on the real line. n n n=1 X Proof. Assume that f(x) is bounded on R then it would be almost periodic by the antiderivative theorem (4.2) and the remark that f′(x) = F(x). Since f(x) is odd π its mean is zero. This is in contradiction with the limit given by the Lemma 2 (4.1). ∞ χ(n) x Remark 4.3. The same analysis applies to the series sin( ), χ being a n n n=1 Dirichlet character modulo N. X 16 AHMED SEBBAR AND ROGER GAY We will need to consider some Bessel functions. We recall that for Res > 0 the Γ-function is ∞ s 1 u Γ(s)= u − e− du. Z0 By Fubini’s theorem 2 ∞ ∞ s 1 (u+v) ∞ s 1 Γ (s)= (uv) − e− dudv = u − ξ0(u) du, Z0 Z0 Z0 where 2t√u ∞ 2e ξ0(u)= dt. 2 1 √t 1 Z − More generally the iterated integrals [31], [32] ∞ ∞ ξ1(x)= ξ0(t) dt, ,ξk(x)= ξk 1(t) dt. ··· − Zx Zx satisfy the differential equation of Bessel type: d2ξ (x) dξ (x) x k + (1 k) k ξ (x)=0. dx2 − dx − k The ordinary Bessel function of order ν is ∞ m 2m+ν ( 1) (z/2) ν J (z)= − , I (z)= i− J (iz), z < . ν m!Γ(m +1+ ν) ν ν | | ∞ m=0 X The K-Bessel function of order ν, for ν not an integer, is π I ν(z) Iν (z) K (z)= − − . ν 2 sin πν When ν is an integer we take the limiting value. It could be also defined by 1 ∞ ν 1 z/2(t+1/t) (4.3) K (z)= t − e− dt, Reν 0. ν 2 ≥ Z0 The Mellin transform of the J0-Bessel function is: ∞ s 1 s Γ(s) J (√x)x − dx =4 . 0 Γ(1 s) Z0 − We will need two Mellin transforms, due essentially to Voronoi ∞ s 1 1 2s 2 x − K (4π√x) dx = (2π)− Γ (s), 0 2 Z0 ∞ s 1 1 2s 2 x − Y (4π√x) dx = (2π)− cos(πs)Γ (s). 0 −π Z0 ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 17 X n n 4.0.1. Summations formulas and beyond. Various classicaln summation=1 formulas, as Poisson summation formula, Voronoi summation formula or Hardy-Ramanujan sum- mation formula can all be given a unified formulation. The following Generalized Poisson summation formula is proved in [10] Theorem 4.3. Let a = a(n) be an arithmetic function with moderate growth. We define the Dirichlet series ∞ s L(a, s)= a(n)n− , Res> 1 n=1 X and we suppose that L(a, s) has an analytic continuation to C with only a possible pole at s = 1. We suppose also that there are positive constants A, a1, , ag such that with the Γ-factors ··· g s γ(s)= A Γ(ajs) j=1 Y L(a, s) satisfies the functional equation γ(s)L(a, s)= γ(1 s)L(a, 1 s). − − Furthermore for f (R), the Schwartz space, we define a very special Hankel’ s transform: ∈ S ∞ γ(s) s g(x)= f(y)K(xy)dy, with K(x)= x− ds. 3 γ(1 s) Z0 ZRes= 2 − Then, ∞ ∞ a(n)f(n)= f(0)L(a, 0) + Res (f)(s)L(a, s)+ a(n)g(n) s=1M n=1 n=1 X X where Ress=1 is the evaluation of the residue at s = 1. 4.0.2. 2 classical choices. (1) For a(n) = 1 we have L(a, s)= ζ(s) and s s γ(s)= π− 2 Γ( ), K(x) = 2cos(2πx). 2 We recover Poisson summation formula for even functions in f(x) (R): ∈ S ∞ 1 ∞ ∞ ∞ (4.4) f(n)= f(0) + f(x)dx +2 f(y) cos(2πny)dy. −2 n=1 0 n=1 0 X Z X Z (2) If a(n)= d(n) we have L(d,s)= ζ2(s) and s s 2 γ(s) 2s 2 γ(s)= π− Γ( ) , = (2π)− (2+2cos πs)Γ(s) 2 γ(1 s) − 18 AHMED SEBBAR AND ROGER GAY and K(x)=4K (4π√x) 4Y (4π√x). 0 − 0 We recover Voronoi summation formula ∞ 1 ∞ f(n)d(n)= f(0) + f(x) (2γ + log x) dx+ 4 n=1 0 X Z ∞ ∞ 1 1 d(n) f(y) 4K 4π(ny) 2 2πY 4π(ny) 2 dy. 0 − 0 n=1 0 X Z As a consequence we have Koshliakov’s formula valid for a> 0: 4π ∞ √a γ log +4 d(n)K (2πan) − a 0 n=1 ! X 1 ∞ 2πn = γ log(4πa)+4 d(n)K . √a − 0 a n=1 ! X This formula was proved by Ramanujan about ten years earlier (He did not appeal to Voronoi’s formula) and by many authors later. 4.1. Another function of Hardy and Littlewood. Hardy and Littlewood gave in [17] (p.269) the following relation ∞ 1 z/n (4.5) F(z)= 1 e− = 2 log z +2γ n − n=1 X ∞ 2 K √2nπiz + K √ 2nπiz − 0 0 − n=1 X n o where Rez > 0,γ is Euler’s constant. For z < 1: | | n ∞ 1 z ∞ ( z) (4.6) 1 e− n = ζ(n + 1) − . n − − n! n=1 n=1 X X An immediate consequence of this expansion is obtained by taking real and imagi- nary parts with z = ix, x R, x < 1: ∈ | | ∞ 1 x ∞ 1 x x4k x4k+2 1 cos =2 sin2 = ζ(4k + 1) + ζ(4k + 3) n − n n 2n − (4k)! (4k + 2)! n=1 n=1 k 0 k 0 X X X≥ X≥ 2k k 1 x = ( 1) − ζ(2k + 1) − (2k)! k 0 X≥ ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 19 X n n n=1 ∞ 1 x x4k+1 x4k+3 sin = ζ(4k + 2) ζ(4k + 4) . n n (4k + 1)! − (4k + 3)! n=1 k 0 k 0 X X≥ X≥ More generally we define the series ( z)n G (z)= ζ(n ν) − ν − n! n>Reν+1 X which has a Mellin-Barnes type integral representation when x > 0,c is fixed with Reν +1 1 K 1 lim G(2iπnk)= . K K d2 →∞ k=1 d n X X| 4.2. Laplace transform of a√tJ1(a√t), a> 0, t> 0, another approach to Segal’s formula. In [26] (formula (12)) Segal proves the following result z Theorem 4.4. If g(z) := (1 cos ) then − k k 1 X≥ πz 1 1 2√2kπ√z g(z)= + J (2√2kπ√z). 2 − 2 4 k 1 k 1 X≥ This formula is interesting compared to (4.5), as we have for real z, g(z)=ReF(iz). The proof given in [26] uses a rather elaborated tools such the three Bessel functions J1,J2,J3, the functional equation of the Riemann ζ-function etc. We give here a proof which we think is simpler. 20 AHMED SEBBAR AND ROGER GAY Let πz 1 1 2√2kπ√z g (z)= + J (2√2kπ√z). 1 2 − 2 4 k 1 k 1 X≥ In the Laplace transform + ∞ tp (a√tJ (a√t))(p)= a √tJ (a√t)e− dt, Re p> 0 L 1 1 Z0 we set u2 = t and obtain + ∞ pu2 2 (a√tJ (a√t))(p)=2a J (au)e− u du, Re p> 0. L 1 1 Z0 According to [34] (page 394, formula(4)) we have for with Arg p < π | | 4 + ν 2 2 2 a ∞ p u ν+1 a 2 J (au)e− u du = e− 4p . ν (2p2)ν+1 Z0 Replacing p by √p with Arg p < π and taking ν = 1 we obtain | | 2 + 2 ∞ pu2 2 a a J (au)e− u du = e− 4p . 1 4p2 Z0 Hence 2 2 a −a (a√tJ (a√t))(p)= e 4p Rep> 0. L 1 2p2 Note that a√tJ (a√t) is not in L2([0, + [) since its Laplace transform is not 1 ∞ bounded in the L2-norm on the lines Rep = c. With a =2√2kπ we get 4kπ 2kπ (2√2kπtJ (2√2kπt))(p)= e− p . L 1 p2 As we have πt 1 π 1 ( − )(p)= L 2 2p2 − 2p and, by continuity, the Laplace transform of the sum in πt 1 1 2√2kπt g (t)= − + J (2√2kπt) 1 2 4 k 1 k 1 X≥ is 2 π π k (e− p ) p2 k 1 X≥ which converges in Rep> 0, the Laplace transform of g1(t) is 2π π 1 π e− p π 1 π 1 2 + 2 2π = 2 + 2 2π . 2p − p p 1 e− p 2p − p p e p 1 − − ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 21 X n n On the other hand n=1 t 1 p 1 (1 cos )(p)= 2 2 = 2 2 L − k p − p + k− p(p k + 1) and ∞ 1 (g)(p)= . L p(p2k2 + 1) Xk=1 The equality (g)(p)= (g1)(p) is obtained by using the well known partial fractions decompositionL L z z 2z2 =1 + , z C 2iπZ, ez 1 − 2 z2 +4k2π2 ∈ \ k 1 − X≥ 2π where we have to set z = . Hence g = g by injectivity of Laplace Transform. p 1 4.3. Some Mellin transforms and the cube of theta functions. It has been noticed in [16] (p.14) that the function 1 it R(t)= e n n n t X≤ is very similar to ζ(1 + it) in its asymptotic behaviour as t + . This could → ∞ n2 suggest a link between this function and the theta series ϑ3(q) = q . In this n Z section, following a suggestion of Crandall [11], we would like toX briefly∈ show by considering Mellin transforms an unexpected link to the third power of the (fourth) 2 Jacobi theta function ϑ (q)= ( 1)nqn , q < 1. We define 4 − | | n Z X∈ t ∞ e n (4.8) χ˜(s, t) = − , ns n=1 X t ∞ e n (4.9) χ(s, t) = ( 1)n − . − ns n=1 X These two functions are defined for s C and Res > 1 forχ ˜(s, t), Res > 0 for χ(s, t). They are related by ∈ 1 s χ(s, t)= s 1 χ˜( , t) χ˜(s, t). 2 − 2 − We have ∞ ∞ s 1 2 1 t − χ˜ (s, t) dt = Γ(s) = Γ(s)(ζ(s 1) ζ(s)) (n + m)s − − 0 n,m=1 Z X 22 AHMED SEBBAR AND ROGER GAY ∞ n+m ∞ s 1 2 ( 1) t − χ (s, t) dt = Γ(s) − = Γ(s)(η(s 1) η(s)) , (n + m)s − − 0 n,m=1 Z X where ∞ ( 1)n η(s)= − (2n + 1)s n=0 X is the Dirichlet η-function. Furthermore ∞ 1 ∞ s 1 3 1 t − χ˜ (s, t) dt = Γ(s) (pq + qr + rp)s 0 p,q,r=1 Z X and for χ(s, t) we have a more interesting result ∞ p+q+r 1 ∞ s 1 3 ( 1) (4.10) t − χ (s, t) dt = − . Γ(s) (pq + qr + rp)s 0 p,q,r=1 Z X We have the following lemma due to Andrews [3] (p.124) Lemma 4.4. For q < 1 | | n n n2 j2 n j ∞ ( 1) q ∞ q − (1 q )( 1) (4.11) ϑ3(q)= ( 1)nr (n)qn =1+4 − 2 − − 4 − 3 1+ qn − 1+ qn n Z n=1 n=1 X∈ X Xj n | |≤ where r3(n) is the number of representations of n as sum of three squares. According to a result of Fermat an integer is a sum of three squares if and only if it is not of form 4n(8m + 7). There are some gaps in the expansion in power series of the left hand side of (4.11). Similarly to (4.10) we have ∞ p+q+r 1 ∞ s 1 3 ′ ( 1) (4.12) t − ϑ (q) 1 dt = − . Γ(s) 4 − (p2 + q2 + r2)s Z0 p,q,r Z X∈ A link between (4.10) and (4.12) is given by Crandall [11] (p.372) as a relation between two Epstein zeta functions associated with the two not equivalent ternary forms 2 2 2 q1(u,v,w)= u + v + w , q2(u,v,w)= uv + vw + wu in the form ∞ p+q+r ∞ p+q+r ′ ( 1) 1 s 2 2 ( 1) (4.13) − = 6(1 2 − ) ζ (s) 4 − . (p2 + q2 + r2)s − − − (pq + qr + rp)s p,q,r Z p,q,r=1 X∈ X Next we establish a functional equation ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 23 X n n 1 n=1 Theorem 4.5. For t> 0 the function χ( , t) satisfies the following functional equa- 2 tion √ ∞ n γ 2πdt 1 ( 1) t/n e− (4.14) χ( , t)= − e− = √i , 2 √n √ n=1 d X XO with γ =1 i and is the set of odd integers. − O We could also seek for a result similar to (4.14) for ∞ 1 t/n (4.15)χ ˜(s, t)= e− , Res> 1. ns n=1 X The relevance of this function lies in its relation to a Hardy-Littlewood-Flett like function: ∞ 1 ∞ 1 t Im˜χ(1, it) = Im eit/n = sin( ). − n n n n=1 n=1 X X In order to prove (4.14) we consider, for a fixed t> 0, the function e t/x f(x)= eiπx − , x> 0 √x extended to the origin by f(0) = 0 and to R as an even function. The obtained function is ∞ on the real line to which we apply the Poisson summation formula (4.4) to obtainC ∞ n t/x ∞ t/x ( 1) t/n ∞ iπx e− ∞ iπx e− (4.16) − e− = e dx +2 e cos(2πnx)dx. √n √x √x n=1 0 n=1 0 X Z X Z −t e y Remark 4.1. The function y is continuous on [0, [ and decreases to 0 at 7→ √y ∞ ∞ infinity, hence the proper integrals f(y) cos(2πny)dy,n 0 are convergent. ≥ Z0 We compute (f)(n) as follows F −t ∞ e y (f)(n)= eiπy cos(2nπy) dy F √y Z0 1 ∞ (iπ+2iπn)y t/y 1/2 ∞ (iπ 2iπn)y t/y 1/2 (4.17) = e − y− dy + e − − y− dy . 2 Z0 Z0 We recall the modified Bessel function (4.3), written in the form ν/2 ∞ ν 1 w a/w 1 w − e− − dw =2 K 2√a a ν Z0 24 AHMED SEBBAR AND ROGER GAY that we use in the form ν/2 ∞ ν 1 bw a/w a (4.18) w − e− − dw =2 Kν 2√ab . 0 b Z Actually for (4.17) we need only the simplest case of π z (4.19) K 1 (z)= e− . 2 2z r The first integral in the left hand side of (4.17), with a = t, b = iπ +2iπn = iπ(2n + 1) − is equal to t 1/4 (4.20) 2 K 2 t ( iπ(2n + 1)) . iπ(2n + 1) 1/2 − − The second integral with p a = t, b = iπ 2iπn = iπ( 2n + 1) − − − is equal to t 1/4 (4.21) 2 K 2 t ( iπ( 2n + 1)) . iπ( 2n + 1) 1/2 − − − − By using (4.19) we see that (4.17) is the sump of 1/4 t π 2√t( iπ(2n+1)) e− − iπ(2n + 1) 4 t ( iπ(2n + 1)) − s − and p 1/4 t π 2√t( iπ( 2n+1)) e− − − . iπ( 2n + 1) 4 t ( iπ( 2n + 1)) − − s − − As in [11] we denote by γ =1 i, dp= 2n +1, n N⋆ with d = i d . Then d describes 1 and − ± ∈ −| | | | O\{ } p p 1/4 γ√tπd t π 2√ itπd 1 e− e− − = √i . iπd 4 t ( iπd) 2 √d − s − Hence p γ√tπd ∞ ∞ e (4.22) 2 f(y) cos(2πny)dy = √i − . √ n=1 0 d , d=1 d X Z ∈OX6 For the remaining term in (4.16) we use (4.20), with n = 0, to obtain ∞ γ√2πt f(x)dx = (f)(0) = √i e− F Z0 ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 25 X n n which together with (4.22) gives finally (4.14): n=1 1 e γ√2πdt χ( , t)= √i − . 2 √d XO In close analogy to Jacobi’s transformation of Theta functions (4.14) appears as a convergence acceleration of a slowly convergent series. 1 t2 Incidentally χ( ; ) is a Fourier transform of a function of the Schwartz class. 2 4 1 Indeed let g(x)= , the reciprocity formulas are 1+ ex2 1 1 itx itx − gˆ(t)= x2 e dx, g(x)= gˆ(t)e dt 2π R 1+ e R Z Z with 1 n itx nx2 gˆ(t)= ( 1) e e− dx −2π − R n>0 Z Xn 2 1 ( 1) t2 1 1 t = − e− 4n = χ( ; ). −2√π √n −2√π 2 4 n>0 X Remark 4.2. The convolution of three functions f,g,h (R) is, as well known, ∈ S (f ⋆ (g ⋆ h))(x) = f(y)(g ⋆ h)(x y)dy R − Z = f(y) g(z)h(x y z) dy R R − − Z Z = f(y)g(z)h(x y z)dzdy. R R − − Z Z × 1 With f(x)= g(x)= h(x)= we have 1+ ex2 dtdu (f ⋆ (g ⋆ h))(x)= y2 z2 (x y z)2 R2 (1 + e )(1 + e )(1 + e ) Z Z − − so that for the Fourier transform 1 1 t2 √π 1 t2 f⋆g⋆h\ (t)=(2π)2gˆ3(t)=(2π)2 χ3( ; )= χ3( ; ) 8π√π 2 4 2 2 4 or 2 √π 3 1 t itx f⋆g⋆h(x)= χ ( ; )e− dt. 2 R 2 4 Z Evaluating at x = 0 we obtain dy dz dy dz y2 z2 ( y z)2 = y2 z2 (y z)2 R2 (1 + e )(1 + e )(1 + e ) R2 (1 + e )(1 + e )(1 + e ) Z Z − − Z Z − 26 AHMED SEBBAR AND ROGER GAY 2 √π 1 t ∞ 1 1 = χ3( ; ) dt = √π χ3( ; u) du. 2 R 2 4 0 √u 2 1 Z Z From (4.10), with s = 2 , we have (Compare with [11]) ∞ ( 1)p+q+r 1 dy dz (4.23) − 1 = y2 z2 (y z)2 . 2 −π R2 (1 + e )(1 + e )(1 + e ) p,q,r=1 (pq + qr + rp) − X Z Z We end this study by using an interesting integral representation due to Mellin [21] (p. 22-23): κ1+i κq+i Γ(s) 1 ∞ ∞ Γ(s z z ) (4.24) = 1 q s q −s z1 ···zq (w0 + w1 + + wq) (2iπ) κ1 i ··· κq i w − ··· ··· Z − ∞ Z − ∞ 0 Γ(z1) Γ(zq) z1 ··· zq dz1 dzq, w wq ··· 1 ··· κ > 0, ν =1, , q; Res > κ + + κ > 0. ν ··· 1 ··· q In the case of q = 2 we obtain at once, as in (4.23) p+q+r κ1+i κ2+i ( 1) 1 ∞ ∞ − = Γ(s u u )Γ(u )Γ(u ) (pq + qr + rp)s − 4π2 Γ(s) − 1 2 1 2 p,q,r 1 κ1 i κ2 i X≥ Z − ∞ Z − ∞ ( 1)p ( 1)q ( 1)r − − − du1du2 ps+u1 u1u2 ps+u2 u1u2 ru1+u2 p 1 − q 1 − r 1 X≥ X≥ X≥ κ1+i κ2+i 1 ∞ ∞ = 2 Γ(s u1u2)Γ(u1)Γ(u2)K(s; u1,u2) du1du2, −4π κ1 i κ2 i − Z − ∞ Z − ∞ where K(s; u ,u )= η(s + u u u )η(s + u u u )η(u + u ). 1 2 1 − 1 2 2 − 1 2 1 2 This representation of the Epstein zeta function of the ternary form q2(u,v,w) = uv + vw + wu in terms of the Dirichlet η-function and similar other representations can shed some light on their analytic continuation. References [1] J. Alcantara-Bode, An integral equation formulation of the Riemann Hypothesis. Integr. Equations Oper, Theory, 17 (2) (1993), 151-168. [2] J. Alcantara-Bode, Some properties of the Beurling function, Pro-Mathematica, XIV, No. 27-28 (2000), 1-11. [3] G. E. Andrew, The fifth and seventh order mock theta functions, Trans. Amer. Math. Soc.293 (1986), no. 1, 113-134. [4] L. B´aez-Duarte, A strengthening of the Nyman–Beurling criterion for the Riemann hypoth- esis, Atti Accad. Naz. Lincei Cl. Sci. Fis. Mat. Natur. Rend. Lincei (9) Mat. Appl. 14 (1) (2003) 5-11. [5] M. Balazard, E. Saias, Notes sur la fonction ζ de Riemann, 1, Adv. Math. 139 (1998) 310-321. ∞ 1 x ARITHMETIC AND ANALYSIS OF THE SERIES sin 27 X n n n=1 [6] M. Balazard, E. Saias, Notes sur la fonction ζ de Riemann, 4, Adv. Math. 188 (2004) 69-86. [7] A. Beurling, A closure problem related to the Riemann zeta-function, Proc. Natl. Acad. Sci.41 (1955) 312-314. [8] J. Boersma, J., R.J. Chapman, Solution to problem 94-12(by M. L. Glasser): A nonharmonic trigonometric series. SIAM Review, 37(3), (1995) 443-448. [9] H. Bohr, Almost Periodic Functions. Julius Springer, Berlin (1933) (Chelsea Publishing Com- pany, N.Y., 1947). [10] H. Cohen, Number theory, volume II: Analytic and modern tools, Graduate Texts in Mathe- matics, vol. 240, Springer, 2007. [11] R. E. Crandall, New representations for the Madelung constant, Experimental Mathematics, 8 (1999) 367-379. [12] H. Davenport, On some infinite series involving arithmetical functions, Quarterly J. Math (2) 8 (1937), 8-13. [13] H. Davenport, On some infinite series involving arithmetical functions (II), Quarterly J. Math(2)8(1937), 313-320. [14] H. Delange, Sur la fonction f(x) = 1/n sin(x/n). Th´eorie analytique et ´el´ementaire des nombres, Caen, 29-30, Septembre 1980,X Journ´ees math´ematiques SMF-CNRS, 1980. [15] W. F. Donoghue, Distributions and Fourier Transforms, Academic, N.Y., 1969. [16] T.M. Flett, On the function (1/n) sin(t/n). J. London Math. Soc., 25:5-9, 1950. [17] G.H. Hardy and J.E. Littlewood,X Notes on the theory of series (XX): on Lambert series. Proc. London Math. Soc., (2) 41 (1936), 257-270. [18] P. Hartman, On a class of arithmetical Fourier series, Amer. J. Math. 60 (1938), 66-74. [19] M. Katsurada, On Mellin-Barnes Type of Integrals and Sums associated With the Riemann Zeta function. Publ. Inst. Math.(Beograd) (Nouvelle Ser.) 62 (76), (1997) 13-25. [20] M. Katsurada, Power series with the Riemann zeta-function in the coefficients, Proc. Japan Acad. 72 (1996), 61-63. [21] H. Mellin, Eine Formel fur den Logarithmus transcendenter Funktionen von endlichen Geschlecht, Acta Sot. Sci. Fennicae 29 (1900) 3-49. [22] L. J. Mordell, Integral formulae of arithmetical character, Journal London Math. Soc, 33 (1958), 371-375. [23] N. Nikolski, Distance formulae and invariant subspaces, with an application to localization of zeros of the Riemann ζ-function. Ann. Inst. Fourier (Grenoble), 45 (1)(1995), 143-159. [24] N. Nikolski, Hardy Spaces (Cambridge Studies in Advanced Mathematics). Cambridge: Cam- bridge University Press.(2019). x [25] B. Safari, R.C. Vaughan, On the fractional parts of and related sequences II, Ann. Inst. n Fourier (Grenoble) 27, 2 (1977), 1-30. [26] S.L. Segal, On f(x)= 1/n sin(x/n), J. London Math. Soc., (2)4, 385-393(1972). [27] G. Tenenbaum, IntroductionX to analytic and probabilistic number theory, volume 163, Grad- uate Studies in Mathematics. American Mathematical Society, Providence, RI, third edition, 2015. Translated from the 2008 French edition by Patrick D. F. Ion ∞ ζ(n + 2) n [28] J. Tennenbaum, On the function x , Math. Scand. 41 (1977), 242-248. n=0 n! [29] E.C. Titchmarsh, The theory of theX Riemann zeta-function, Clarendon Press, Oxford, 1951 [30] V.I. Vasyunin, On a biorthogonal system related with the Riemann hypothesis, St. Petersburg Math. J. 7 (1996) 405-419. 28 AHMED SEBBAR AND ROGER GAY [31] G. Voronoi, Sur une fonction transcendante et ses applications `ala sommation de quelques s´eries, Annales Scientifiques de l’Ecole´ Normale Sup´erieure, troisi`eme s´erie 21 (1904), 207- 267. [32] G. Voronoi, Sur une fonction transcendante et ses applications `ala sommation de quelques s´eries, seconde partie, Annales Scientifiques de l’Ecole´ Normale Sup´erieure, troisi`eme s´erie 21 (1904), 459-533. [33] A. Walfisz, Teilerprobleme, Math. Z. 26 (1927) 66-88. [34] G.N. Watson, A Treatise on the Theory of Bessel Functions, Cambridge University Press (1922) Chapman University, One University Drive,, Orange CA 92866, USA Email address: [email protected] Universite´ de Bordeaux, IMB, UMR 5251, F-33405 Talence, France Email address: [email protected] Universite´ de Bordeaux, IMB, UMR 5251, F-33405 Talence, France Email address: [email protected]