arxiv: v1 [math.nt] 6 Jun 2011

Transcription

1 VARIATIONS OF THE RAMANUJAN POLYNOMIALS AND REMARKS ON ζ(j +1)/π j+1 arxiv: v1 [mathnt] 6 Jun 011 MATILDE LALÍN AND MATHEW ROGERS Abstract We observe that five polynomial families have all of their zeros on the unit circle We prove the statements explicitly for four of the polynomial families The polynomials have coefficients which involve Bernoulli numbers, Euler numbers, and the odd values of the Riemann zeta function These polynomials are closely related to the Ramanujan polynomials, which were recently introduced by Murty, Smyth and Wang [10] Our proofs rely upon theorems of Schinzel [13], and Lakatos and Losonczi [9] and some generalizations 1 Introduction In a recent paper, Murty, Smyth and Wang considered the Ramanujan polynomials [10] They were defined by Gun, Murty and Rath [7] using R k+1 (z) := B j B k+ j (j)!(k + j)! zj, (1) where B j is the jth Bernoulli number Among other fascinating results, Murty, Smyth and Wang showed that R k+1 (z) has all of its non-real zeros on the unit circle The purpose of this paper is to study some variants of R k+1 (z), which also have many zeros on the unit circle Claim 1 Let B j denote the Bernoulli numbers, and let E j denote the Euler numbers Suppose that k The following polynomials have all of their non-zero roots Date: June 6, Mathematics Subject Classification Primary ; Secondary M L is supported by NSERC Discovery Grant and a start-up grant from Université de Montréal M R is supported by NSF award DMS

2 MATILDE LALÍN AND MATHEW ROGERS on the unit circle: P k (z) := (π) (k)! ( ) k ( 1) j B j B k j z j () j +ζ() ( z +( 1) k z ), Q k (z) := ( k +1 ) P k (z) k P k (z/) P k (z), (3) Y k (z) := π ( Qk (i z)+q k k ( i z) ) (4) W k (z) := ( + ) P k (z) k P k (z/) P k (z), (5) S k (z) := E j E k j ( k j ) z j (6) We will offer a general proof of Claim 1 for Q k (z), Y k (z), W k (z), and S k (z) It seems that P k (z) is more difficult to handle In Section 4 we offer several partial results concerning P k (z) An important secondary goal of this work, is to highlight a connection with the odd values of the Riemann zeta function Recall that the Riemann zeta function is defined by 1 ζ(s) := n s While it is a classical fact that ζ(j)/π j is rational when j 1, very little is known about the arithmetic nature of ζ(j + 1) The only theorems in this direction are celebrated irrationality results For instance, Apéry showed that ζ(3) is irrational [, 11], Rivoal proved that infinitely many odd zeta values are irrational [1], and Zudilin established that at least one element of the set {ζ(5), ζ(7), ζ(9), ζ(11)} is also irrational [14] Therefore, let us briefly consider the identity which gave birth to the Ramanujan polynomials [3, pg 76] The formula can be written as (π) (k)! = ( 1) j B j B k j ( k j n=1 ) n=1 z j + ζ() 1 z n e πn/z 1 +( 1)k+1 n=1 ( ( 1) k z +z ) 1 z n e πnz 1 This identity holds whenever z iq The restriction is necessary to ensure that both infinite series converge Not surprinsigly, this formula is also mentioned in works dealing with the irrationality of ζ(3) [11] Notice that Gun, Murty and Rath used (7) toexpress oddzeta values intermsofeichler integrals[7] Nowconsider the case when k = A brief numerical calculation shows that the polynomial obtained from the left-hand side z 4 +5z +1 90ζ(3) π 3 ( z 3 +z ) = 0, (8) (7)

3 VARIATIONS OF THE RAMANUJAN POLYNOMIALS AND REMARKS ON ζ(j +1)/π j+1 3 has all of its zeros on the unit circle Notice that if we truncate the right-hand side of (7), then we can nearly obtain a closed form for ζ(3): ( ) z 4 +5z +1 π 3 ζ(3) z 3 +z 90, (9) where z 0 = e π/z 1 + z 1 e πz 1 This approximation gives six decimal places of numerical accuracy The accuracy can be increased by including higher order terms in the truncation It would be extremely interesting if this idea could be used to say something about the irrationality of ζ(3)/π 3 Unfortunately such a theorem is well beyond the scope of this paper We will settle for the more modest goal of studying the polynomial families listed in Claim 1 S k (z), Y k (z) and the theorems of Schinzel, Lakatos and Losonczi InthissectionweprovethatS k (z) andy k (z) havealloftheir non-zerorootsonthe unit circle Our proofs follow from applying the theorems of Schinzel [13], Lakatos and Losonczi [9], and Lakatos [8] Lakatos proved that any reciprocal polynomial k A jz j, with real-valued coefficients, which satisfies A k A j A k, (10) must have all of its zeros on the unit circle If the inequality is strict then the polynomial has only simple zeros Equation (10) is a very strong restriction There have been a number of recent improvements to (10) with a similar flavor (see [13] and [9]) Schinzel proved that any self-inversive polynomial which satisfies A k inf c,d C d =1 ca j d k j A k, (11) must have all of its zeros on the unit circle [13] Self-inversive polynomials have complex-valued coefficients which satisfy A j = ǫa k j, for some fixed ǫ = 1 Notice that the class of self-inversive polynomials includes both reciprocal and antireciprocal polynomials In Theorems 1 and we apply Schinzel s theorem with d = 1 Theorem 1 Suppose that k 1 Then all of the zeros of ( ) k S k (z) = E j E k j z j j lie on the unit circle Furthermore, all of the zeros are simple

4 4 MATILDE LALÍN AND MATHEW ROGERS Proof With (11) in mind, let us begin by setting ( ) k A j := E j E k j j The sign of E n alternates with respect to n This implies that all of the coefficients of S k (z) have sign ( 1) k Our proof consists of three main steps First we remove the absolute values from the sum (11) This is easily accomplished by showing that ( 1) k (ca j A k ) > 0 π for c = Next we evaluate k 4(1+3 1 k ) A j explicitly, and finally we deduce the desired upper bound In order to remove the absolute value signs from (11), we need to demonstrate that ( 1) k (ca j A k ) > 0 Using the following bound for Euler numbers [1, pg 805]: 4 n+1 (n)! π n+1 > E n > and the fact that E 0 = 1, leads to 4(1+3 1 k ) π 4 n+1 (n)! π n+1 (1+3 1 n ), A k = 4(1+3 1 k ) E k E 0 π > 4k+ (k)! π k+ ( ) k > E j E k j = A j j The absolute values can be removed because both A k and A j have sign ( 1) k : If we take ( 1) k 4(1+3 1 k ) A k > ( 1) k A j > 0 π c = π 4(1+3 1 k ), then the previous inequality implies ( 1) k A k c( 1) k A j > 0 Let E n (z) denote the classical Euler polynomials, and recall a standard convolution identity [6]: n ( ) n E j (v)e n j (w) = (1 w v)e n (v +w)+e n+1 (v +w) j Set v = w = 1, and then use E ( n = n E 1 n ), to obtain an expression for Sk (1) We have S k (1) = k+1 E k+1 (1) k+1 ( k+ 1) B k+ = k +1

5 VARIATIONS OF THE RAMANUJAN POLYNOMIALS AND REMARKS ON ζ(j +1)/π j+1 5 The evaluation of E k+1 (1) follows from [1, pg 805] Thus ca j A k =( 1) k (A k ca j ) =(k +1)( 1) k A k c( 1) k S k (1) =(k +1) A k c S k (1) =(k +1) A k π ( k+ 1) B k+ (k +1)(1+3 1 k ) To finish theverification of (11), we needto show that thelast expression is bounded from above by A k This is equivalent to showing that π ( k+ 1) B k+ (k +1)(1+3 1 k ) k A k = k E k (1) We will resort to an inequality for Bernoulli numbers [1, pg 805]: (n)! (π) n < B n < (n)! (π) n (1 1 n ) (13) Thus we find π ( k+ 1) B k+ > π ( k+ 1) (k+)! (k +1)(1+3 1 k ) (k +1)(1+3 1 k ) (π) k+ = k+1 (1 k )(k +1)! (14) (1+3 1 k ) π k+1 On the other hand, we have already used the fact that Euler numbers are bounded from above by k E k < k 4k+1 (k)! π k+1 (15) Substituting (14) and (15) into (1), reduces the inequality to k +1 > k k 1 k This final inequality is easily verified with elementary calculus for k 1 Since the inequality is strict, we conclude immediately that S k (z) has only simple zeros which all lie on the unit circle We have proved that S k (z) has all of its zeros on the unit circle Perhaps it is interesting to note that S k (z) satisfies (π/) k+1 S k ( z ) =z k χ 4 (n)sech(πn/z) (k)! n k+1 n=1 +( 1) k n=1 χ 4 (n)sech(πnz/) n k+1, where χ 4 (n) is the non-principle character mod 4 This formula appears immediately after equation (7) in Ramanujan s notebook [3, pg 76] As a result it is

6 6 MATILDE LALÍN AND MATHEW ROGERS easy to approximate the zeros of S k (z) by the zeros of exponential polynomials It remains to be seen whether or not there are any interesting applications for this observation To illustrate this method a second time, we prove that the polynomial Y k (z) has all of its non-zero roots on the unit circle Notice that Y k (z) is a close analogue of S k (z), except that it involves Bernoulli numbers rather than Euler numbers Theorem Suppose that k The polynomial Y k (z) = π k ( Qk (i z)+q k ( iz) ) = πk (k)! ( ) k B j B k j ( j 1)( k j 1) z j, j has all of its non-zero roots on the unit circle Furthermore, all of the zeros are simple Proof Observe that Y k (z) has degree k 1, since the coefficients of z k and z 0 are identically zero We prove that Y k (z)/z satisfies thehypothesis of Schinzel s theorem [13] If we eliminate the trivial factor of z, then we obtain a polynomial of the form where Y k (z) z k = A j z j, ( ) A j = (π)k k (1 j )(1 k+j )B j+ B k j (k)! j + Notice that Y k (z)/z is reciprocal, since A k j = A j By elementary properties of Bernoulli numbers, the sign of A j is ( 1) k for all j Schinzel s theorem can be applied if the following inequality holds: k A k ca j A k, (16) for some c C We prove that (16) holds when c = π (1 k ) Our proof follows 8(1 3 k ) the same three steps as in the case of S k (z) In order to remove the absolute value signs from (16), we need to demonstrate that ( 1) k (ca j A k ) > 0 We demonstrate this by comparing an upper bound on ( 1) k A k, with a lower bound on ( 1) k A j The lower bound on A j is a consequence of an inequality from [5]: In particular we find B n > (n)! (π) n (1 n ) ( 1) k A j = A j > 4 (17)

7 VARIATIONS OF THE RAMANUJAN POLYNOMIALS AND REMARKS ON ζ(j +1)/π j+1 7 By (13), we find an upper bound for A k : ( 1) k A k = A k < π 1 Comparing (18) and (17), allows us to easily conclude 1 k 3 k (18) ( 1) k (ca j A k ) > 0, (19) whenever k > Since we have proved (19), Schinzel s sum immediately reduces to k k ca j A k = (k 1)( 1) k A k +( 1) k c A j (0) Now we simplify the remaining sum Let B j (z) denote the usual Bernoulli polynomials By standard evaluations of Bernoulli polynomials [1, pg 805], we have ( )( ( ) )( ( ) ) A j = (π)k k 1 1 B j+ B j+ (0) B k j B k j (0) 4(k)! j + Next we use a well known convolution identity for Bernoulli polynomials [6]: n ( ) n B j (v)b n j (w) = n(w +v 1)B n 1 (v+w) (n 1)B n (v +w) (1) j Considering all of the cases where (v,w) { ( 1, 1 k A j = (π)k (k)! () ( B k ( 1 ) ( (, 0, 1 ), 1,0),(0,0)}, leads to ) ) B k (0) = (π)k (k)! (k 1)(1 k )B k () Substituting () into (16), leads to a closed form expression for the sum we are interested in: k ca j A k = (k 1)( 1) k A k ( 1) k c (π)k (k)! (k 1)(1 k )B k The proof can be completed by showing that this last expression is bounded from above by A k or Employing π (1 k )(π) k 8(1 3 k ) (k)! (k 1)(1 k ) B k < k A k ζ(n) = ( 1)n+1 B n (π) n (n)! reduces the desired inequality to 1 (k 1)(1 k )ζ(k) < k(1 3 k )ζ(k ) (3)

8 8 MATILDE LALÍN AND MATHEW ROGERS It is elementary to show that this inequality holds for k > 1 3 Generalizing the criteria to other families Conditions such as(11) appear to be too restrictive to apply to polynomial families such as P k (z), W k (z) and Q k (z) In this section we prove that Q k (z) and W k (z) have all their roots on the unit circle, by extending the theorems used in the previous section Let us briefly recall how to derive results such as (10) and (11) For a real-valued reciprocal polynomial V k (z) = k A jz j, the condition immediately implies that A k > ca j A k, (4) 1 > c V k (z) v k (z) A k, (5) where v k (z) = zk+1 1 Notice that if v z 1 k(z) is expanded in a geometric series, then (5) can be derived from (4) as a simple consequence of the triangle inequality Despite the fact that (5) does not imply (4), it turns out that (5) is easily strong enough to conclude that V k (z) has all of its zeros on the unit circle To demonstrate this, first restrict z to the unit circle, and write z = e iθ with θ (0,π), and ( ) sin (k+1)θ ṽ k (θ) = z (k+1)/ v k (z) = sin ( θ ) ( If j < k+ is a positive odd integer, then it is easy to show that ṽ jπ k k+1) has sign ( 1) (j 1)/, and ( ṽk jπ k+1) 1 This implies that ṽk (θ) has at least k+1 interlacing positive and negative values in the interval (0,π), and it has absolute value 1 at each of those points By (5) we can write c A k z (k+1)/ V k (z) = ṽ k (θ)+et, where the error term has absolute value less than 1 It follows that c A k z (k+1)/ V k (z) changes sign at least k times for θ (0,π) By the intermediate value theorem we conclude that V k (z) has at least k zeros on the unit circle Since the polynomial has at most k zeros, all of its zeros must lie on the unit circle Wecaneasily extend thisidea by selecting adifferent v k (z) 1 This typically entails constructing v k (z) to approximate specific polynomial families Definition 1 Let f(θ) : (a,b) R be a continuous function We call f(θ) a kth order alternating function on (a, b), if it assumes alternating positive and negative (or negative and positive) values at points p i, where a < p 1 < < p k+1 < b We say that f(θ) has oscillation distance d, if f(p i ) > d for each i {1,,k +1} Lemma 1 Suppose that f(θ) is a kth order alternating function on (a,b), with oscillation distance d Let F(θ) : (a,b) R be a continuous function such that F(θ) f(θ) < d for all θ Then F(θ) has at least k zeros 1 This principle was inspired by a careful study of the proof in [10]

9 VARIATIONS OF THE RAMANUJAN POLYNOMIALS AND REMARKS ON ζ(j +1)/π j+1 9 Proof This lemma is essentially a restatement of the intermediate value theorem The proof follows immediately from the method described in the previous discussion 31 The zeros of W k (z) lie on the unit circle Themainresultofthissubsection is the following theorem: Theorem 3 Suppose that k The polynomial W k (z) =( +)P k (z) k P k (z/) P k (z) = (π) k ( ) k ( 1) j B j B k j (1 1 j )(1 1 k+j ) z j (k)! j has all of its zeros on the unit circle Furthermore, all of the zeros are simple In order to prove Theorem 3 we first need to establish that a certain trigonometric polynomial possesses the alternating property with oscillation distance 03 Lemma Suppose that k > 10 The function w k (θ) := cos(kθ)+ π 3 cos((k )θ)+ sin((k 3)θ) (1 1 k ) sinθ is an alternating function of order k on ( π, π), with oscillation distance 03 Proof We need to demonstrate that w k (θ) > 3 for k + 1 values of θ ( π,π) We must also show that the sign of w k (θ) alternates over successive points in this set Since w k (θ) is even with respect to θ, and since w k (0) > 3, we only need to demonstrate that there are an additional k such points in (0,π) Suppose that k > 10, let α be defined by ( ) α = 1 π arccos 3 π = 4, (6) 3 and let j 0 = [(k 1)α]+1 (7) We claim that w k (θ) satisfies the necessary conditions on the following set of points: { π S = k 1,, (j } { 0 1)π (j0 1/)π,, (k j } 0 1/)π k 1 k 1 k 1 { } { } (k j0 )π (k )π (k (1 ǫ))π,,, k 1 k 1 k 1 where ǫ > 0 is sufficiently small First expand w k (θ) using trigonometric identities ( ) π w k (θ) =cos((k 1)θ)cos(θ) (( 1 k ) ) π +sin((k 1)θ) 3 4 sinθ+ cscθ 1 1 k 1 1 k

10 10 MATILDE LALÍN AND MATHEW ROGERS Notice that S (essentially) arises from cases where either cos((k 1)θ) = 0, or sin((k 1)θ) = 0 Begin by considering the cases where θ = jπ and j {1,,k } Then ( ) ( )( ) jπ jπ π w k = ( 1) j cos k 1 k (8) 1 1 k ( In order to have w jπ k ) > 03, we need to restrict j so that ( ) jπ cos 03 > k 1 π k In other words we must have where j k 1 (α k,1 α k ), α k = 1 π arccos ( 03 π k Since k > 10 we have (k 1)(α k α) 1 Therefore it is sufficient that j k 1 (α,1 α), where αisdefined in(6) Thisimplies thatj {1,,,j 0 1} {k j ( 0,,k }, withj 0 defined in(7) If j {1,,,j 0 1}, thenby (8) w jπ k ) hassign ( 1) j If j {k ( j 0,,k } then the cosine in (8) contributes an extra minus sign, and w jπ k ) has sign ( 1) j+1 Now consider the case where θ = (j 1/)π and j {j () 0,,k j 0 } By elementary properties of trigonometric functions, w k (θ) reduces to w k ( (j 1/)π k 1 ) = ( 1) j+1 (( π k + 1 ) 1 k csc ) ( ) (j 1/)π sin k 1 ( )) (j 1/)π k 1 In order to place a lower bound on this expression, first choose an interval (β,1 β), which contains the set of rational numbers { j 0 1/,, k j 0 1/ } This can be accomplished by selecting { α if j0 > α(k 1)+ 1 β =, α 1 if j () 0 < α(k 1)+ 1 Notice that one of these situations must occur, because (7) guarantees that j 0 (α(k 1),α(k 1)+1) We obtain the following lower bound from (9): ( ) ( ) (j 1/)π w k min π k 1 θ (πβ,π(1 β)) 3 4 sinθ + cscθ 1 1 k 1 1 k (9)

11 VARIATIONS OF THE RAMANUJAN POLYNOMIALS AND REMARKS ON ζ(j +1)/π j+1 11 The right-hand side is minimized at the end points of the interval (πβ,π(1 β)), so it follows that ( ) ( ) (j 1/)π w k π k sinπβ + cscπβ 1 1 k 1 1 k Consider both choices of β, and recall the assumption that k > 10 A few easy calculations are sufficient to obtain ( ) { (j 1/)π 57 if w j0 > α(k 1)+ k > 1, k 1 34 if j 0 < α(k 1)+ 1, for( all values) of j {j 0,,k j 0 } It is easy to deduce from (9) that the sign of w (j 1/)π k is precisely ( 1) j ( ) Finally consider the value of w (k (1 ǫ))π k Notice that w k (π) = ( 1) k ( π 3 ) k k Since k > 10 it follows easily ( that w) k (π) > 19, and w k (π) has sign ( 1) k If ǫ is sufficiently small then w (k (1 ǫ))π k also has sign ( 1) k, and absolute value much larger than 3 To briefly summarize the sign values of w k (θ), we have the following table: θ Sign(w k (θ)) 0 ( 1) 0 π ( 1) 1 π ( 1) (j 0 1)π ( 1) j 0 1 (j 0 1/)π ( 1) j 0 (k j 0 1/)π ( 1) k j 0 (k j 0 )π ( 1) k j 0+1 (k )π (k (1 ǫ))π ( 1) ( 1) k This table shows that w k (θ) changes sign at least k times over the interval (0,π) Now we use Lemma to prove our main result

12 1 MATILDE LALÍN AND MATHEW ROGERS Proof (Theorem 3) Let us define A j using W k (iz) = (π) k (k)! = A j z j ( ) k B j B k j (1 1 j )(1 1 k+j ) z j j By Lemma 1 it suffices to prove that z k W k (iz) A 0 w k (z) < 03 where w k (z) = (z k +z k )+ π 6 (zk +z k z k 3 z 3 k )+, (30) (1 1 k ) z z 1 and z = e iθ Thus we write z k W k (iz) w k (z) A 0 = A j z j k (z k +z k ) π A 0 6 (zk +z k ) z k 3 z 3 k (1 1 k ) z z 1 A 1 π k A A j A 0 (1 1 k ), (31) j= where the second step makes use of a geometric series and the triangle inequality If we recall that B 0 = 1/, and use both inequalities from (13), then we find A j = B ( jb k j(1 1 j )(1 1 k+j ) k ) j A 0 B k (1 1 k ) < (1 1 k ) Additionally we have, by (3), A 1 π A 0 6 = π 6 ζ(k ) ζ(k) π 6 This second expression goes to zero as k For example, it is not hard to see that the absolute value is less than 001 for k > 4

13 VARIATIONS OF THE RAMANUJAN POLYNOMIALS AND REMARKS ON ζ(j +1)/π j+1 13 Therefore we can remove the absolute value signs from (31) We find that z k W k (iz) A 0 w k (z) A 1 π A 0 3 (k 3) + 1 k 1 k j= A j A 0 =+4 A 1 π (k 3) + A A j 1 k A 0 =+ π 3 =+ π 3 3+ ζ(k ) ζ(k) ζ(k ) ζ(k) ( 1 k π 3 π k (k 3) k k 1 1 k 3 5 ) π k Notice that we evaluated j A j using the same Bernoulli convolution identity (1) In addition, we have used the inequality 1 1 n < ζ(n) < n, which are easy to deduce from the Euler product formula and the Dirichlet eta function As k this final upper bound approaches a limit of π 3 89 It 3 is easy to see that it becomes < 03 for k > 6 In summary, we have proved the theorem for k > 10 The cases for k 10 are easily checked with the numerical method outlined in section 4 3 The zeros of Q k (z) lie on the unit circle Notice that the coefficients of Q k (z) involve the odd values of the Riemann zeta function The primary goal of this subsection is to prove the following theorem: Theorem 4 Suppose that k The polynomial Q k (z) =( k +1)P k (z) k P k (z/) P k (z) ( ) = (π) k ( 1) j B j B k j ( j 1)( k j 1) z j (k)! j +ζ()( 1)(( 1) k z +z ), has all of its non-zero roots on the unit circle Furthermore, all of the zeros are simple As in the proof of Theorem 3, the first step is to construct a trigonometric polynomial which approximates Q k (z) Notice that Q k (z) has degree k 1, and that it has a trivial zero at z = 0 Therefore we need to prove that it has k zeros on the unit circle

14 14 MATILDE LALÍN AND MATHEW ROGERS Lemma 3 Suppose that k > 5 Then q k (θ) := cos((k )θ)+ 4 π sin((k 1)θ)+ 8(1 3 k ) sin((k 3)θ) π (1 k ) sinθ is an alternating function of order k on ( π,π), with oscillation distance 003 Proof We need to demonstrate that q k (θ) > 03 for k 1 values of θ ( π,π) We must also show that the sign of q k (θ) alternates over successive points in this subset The proof is similar to the proof of Lemma, so we will be brief Let α be defined by α = 1 ( ) 03 π arccos 16 = 47, (3) π and let j 0 = [(k 1)α]+1 (33) We claim that q k (θ) > 03 on the following set of k +1 points: { S ={0} { ± π k 1,,±(j 0 1)π k 1 ± (k j 0)π )π,,±(k k 1 k 1 } { ± (j 0 1/)π,,± (k j } 0 1/)π k 1 k 1 } { } (k (1 ǫ))π ± k 1 If we consider S \{ (j 0 1/)π, (k j 0 1/)π }, then we obtain a subset of k 1 points where the sign of q k (θ) alternates over successive points In order to prove this claim, first expand q k (θ) using trigonometric identities ( ) q k (θ) =cos((k 1)θ)cos(θ) 1 8(1 3 k ) π (1 k ) ( ( ) +sin((k 1)θ) sin(θ) 1 8(1 3 k ) + 4 ) π (1 k ) π + 8(1 3 k ) π (1 k ) cscθ Now consider the cases where θ = jπ, with (k ) j (k ) and j 0 We have ( ) ( )( ) jπ jπ q k = ( 1) j cos 1 8(1 3 k ) (34) k k 1 π (1 k ) ( To ensure that q jπ k ) > 003, we need to restrict j so that ( ) jπ cos 003 > ( ) k 1 1 8(1 3 k ) π (1 k ) By similar reasoning as in the proof of Lemma, it is sufficient that j k 1 ( (1 α), α) (α,1 α), where α is defined in (3) This immediately implies that j {±1,,±(j 0 1)} {±(k j 0 ),,±(k )}, where j 0 is defined in (33) A careful inspection of (34) reveals that the function has sign ( 1) j for j {±1,,±(j 0 1)}, and sign ( 1) j+1 for j {±(k j 0 ),,±(k )}

15 VARIATIONS OF THE RAMANUJAN POLYNOMIALS AND REMARKS ON ζ(j +1)/π j+1 15 Next consider the cases where θ = ± (j 1/)π and j {j () 0,,(k j 0 )} We obtain ( ) ( ( )( ) (j 1/)π (j 1/)π q k ± = ( 1) j 1 ±sin 1 8(1 3 k ) (k 1) (k 1) π (1 k ( ) ± 8(1 3 k ) (j 1/)π π (1 k ) csc + 4 ) (35) (k 1) π { } Now select β so that j0 1/,, k j 0 1/ (β,1 β) Following Lemma, this is accomplished by selecting { α if j0 > α(k 1)+ 1 β =, α 1 if j () 0 < α(k 1)+ 1 Therefore we obtain ( ) ( ) q (j 1/)π k ± min (k 1) θ (πβ,π(1 β)) 1 8(1 3 k ) sinθ+ 8(1 3 k ) π (1 k ) π (1 k ) cscθ ± 4 π ( ) 1 8(1 3 k ) sinπβ + 8(1 3 k ) π (1 k ) π (1 k ) cscπβ ± 4 π Checking both possible values of β, and both possible signs of ±, leads to a lower bound which holds for k > 5: ( ) q (j 1)π k ± > 008 (k 1) The final signs are summarized in the table below The only remaining cases are when j {0} {± (k (1 ǫ))π } These cases can be easily dispensed with by elementary properties of trigonometric functions To briefly summarize the sign values of q k (θ), we have the following table: θ Sign(q k (θ)) θ Sign(q k (θ)) (k (1 ǫ))π ( 1) k 0 ( 1) 0 (k )π ( 1) π ( 1) 1 (k j 0)π π ( 1) ( 1) k j 0+1 (k j 0 1/)π ( 1) k j 0 (j 0 1)π (j 0 1/)π (j 0 1/)π ( 1) j 0 1 ( 1) j 0 1 ( 1) j 0 (j 0 1)π ( 1) j 0 1 (k j 0 1/)π ( 1) k j 0 1 (k j 0 )π ( 1) k j 0+1 π ( 1) π ( 1) 1 (k )π (k (1 ǫ))π ( 1) ( 1) k

16 16 MATILDE LALÍN AND MATHEW ROGERS Notice that there are precisely k + 1 values of θ in this table If we exclude the cases where θ { (j 0 1/)π, (k j 0 1/)π }, then the sine of q k (θ) alternates over the remaining k 1 values of θ Next we use Lemma 3 to establish that Q k (z) has all of its non-zero roots on the unit circle for k Proof (Theorem 4) Let us define A j using Q k (iz) = (π) (k)! ( ) k B j B k j ( j 1)( k j 1) z j j +i( 1) k ζ()( 1)(z z ) = A j z j +i( 1) k ζ()( 1)(z z ) In order to simplify the following analysis, we have intentionally defined A j to only involve the even coefficients of Q k (iz) Notice that A 0 = A k = 0, and that A 1 = (π) B k ( k k(k 1) 1) (k)! =( 1) k ζ(k )( k 1) Suppose that k > 5 Then by Lemma 1 it suffices to prove that z k Q k (iz) A 1 q k (z) < 003 where q k (z) = (z k +z k ) i π (z z 1 k )+ 8 (1 3 k ) π (1 k ) z k 3 z 3 k z z 1, (36) and z = e iθ Therefore let us begin by writing z k Q z (iz) q k (z) A 1 = A j z j k +( 1) k i ζ(k 1)( 1) (z 1 k z ) A 1 A 1 (z k +z k )+ i π (z z 1 k ) 8 (1 3 k ) z k 3 z 3 k π (1 k ) z z 1 k A j 8 (1 3 k ) A 1 π (1 k ) + ( 1) kζ()( 1) π j= The second step follows from substituting a geometric series, and then applying the triangle inequality We know from equation (19) (after noting the change in definition of A j ), that A j A 1 > 8 π (1 3 k ) (1 k ) A 1

17 VARIATIONS OF THE RAMANUJAN POLYNOMIALS AND REMARKS ON ζ(j +1)/π j+1 17 In addition ( 1)kζ()( 1) A 1 π = ζ()(1 1 k ) πζ(k )(1 k ) π This limit tends to zero A simple calculation shows that the quantity is less than 0001 for k > 3 Therefore we can remove the absolute value signs from the inequality We are left with z k Q z (iz) q k (z) A 1 = k j= j=1 A j A 1 (k 3) 8 π (1 3 k ) (1 k ) + 4 π ( 1)kζ()( 1) A 1 A j (k 3) 8 (1 3 k ) A 1 π (1 k ) + 4 π 4 ζ()(1 1 k ) πζ(k )(1 k ) = (k 1) 4 π ζ(k)(1 k ) ζ(k )(1 k ) (k 3) 8(1 3 k ) π (1 k ) + 4 π 4 π ζ()(1 1 k ) ζ(k )(1 k ) As usual, we have evaluated j A j using (1) The limit of the upper bound is 0 It is easy to see that it becomes < 003 for k > 6 The cases for k < 6 are π easily proved with the numerical method described in Section 4 We conclude this section by deriving a second approximation for ζ(3)/π 3 If we begin with the expression for Q in terms of P, and then substitute a truncated version of (7) for P, we can obtain ζ(3) z π 3 1+z 14, (37) where z 0 = e 4πz 1 + 8z3 e π/z 1 17z e πz 1 17z3 e π/z 1 + 8z e πz 1 + z3 e 4π/z 1 (38) Selecting the zero given by z 9 39i, yields 4 decimal places of accuracy in (37) Notice that this approximation is slightly worse than (9) 4 Partial results on P k (z) We have made a number of unsuccessful attempts to apply the theorems of Schinzel, Lakatos and Losonczi, and their generalizations to the case of P k (z) A piece of evidence indicating that these methods may not work is given by the result The most general result that we can prove is that P k (z) has at least k 1 zeros in half of the unit circle, using the construction from [10]

18 18 MATILDE LALÍN AND MATHEW ROGERS in [10] which shows that R k+1 (z) has four zeros that do not lie in the unit circle (by comparison Y k (z) has all of its roots on the unit circle) We will briefly describe one instance were Lakatos s condition (10) fails, because it leads to an interesting formula Let us define the A j s as follows: 4 A j z j = P k (iz) = ( π ) ζ(j)ζ(k j)z j +ζ (k 1) ( z z ) Notice that P k (iz) has all of its zeros on the unit circle, if and only if P k (z) also has all of its zeros on the unit circle Computational experiments helped us to make the following observation: Observation 1 Suppose that k, then 4k(k 1) A 4k = 4 A 4k A j (39) Formula (39) can be proved with the identities for Bernoulli numbers that we used in Theorem This immediately rules out the possibility of applying (10) (condition (11) can also be ruled out by slightly different methods) It is curious to note that the right-hand side of (39) appears to involve odd zeta values, whereas the left-hand side does not It turns out that when (39) is explicitly calculated, the odd zeta values drop out Theorem 5 Suppose that k < 1000 Then all of the zeros of P k (z) lie on the unit circle Furthermore, all of the zeros are simple While we have not been able to prove a general theorem concerning P k (z), we have been able to prove Theorem 5 for k < 1000 The proof uses a standard computational method based on the intermediate value theorem Notice that the map z +z 1 u, sends the unit circle to the real interval [ 1, 1] Under this transformation, we also have z k +z k T k (u), where T k (u) is the usual Chebyshev Polynomial If we write (z k +( 1) k )P k (z) = z k (P k (z)+p k (1/z)) = z k Pk (u), then it follows that P k(z) has all of its zeros on the unit circle, if and only if Pk (u) :=(π) (k)! ( 1) j B j B k j ( k j ) T j (u) +ζ() ( T (u)+( 1) k T 1 (u) ),

19 VARIATIONS OF THE RAMANUJAN POLYNOMIALS AND REMARKS ON ζ(j +1)/π j+1 19 has all of its zeros in the interval [ 1,1] It is easy to count real zeros of real-valued polynomials The intermediate value theorem allows one to find zeros by detecting sign changes Since Pk (u) has degree k, it is only necessary to detect k sign changes in [ 1,1] (fewer sign changes are required if zeros lie at u = ±1) We have successfully carried out these calculations for k < Conclusion In conclusion, we have shown that S k (z), Y k (z), W k (z), and Q k (z) have all of their non-zero roots on the unit circle These polynomials have a strong connection to the Ramanujan polynomials We were disappointed that we were unable to to deduce a similar theorem concerning P k (z), however we are hopeful that the approach outlined in Section 3 might eventually succeed in this case An additional avenue might involve studying the zeros of a truncated version of the right-hand side of (7) Notice that the zeros of P k (z) are very well approximated by the zeros of 0 = z z 1 k e π/z 1 +( 1)k+1 e πz 1 Thus, it should be a worthwhile endeavor to study the zeros of these auxiliary functions Acknowledgements The authors thank Wadim Zudilin for the interesting discussions, encouragement and helpful comments, and M Ram Murty for calling our attention to [7] A portion of this research was carried out while Mat Rogers was visiting the University of Georgia He is grateful for their hospitality References [1] M Abramowitz and I A Stegun, Handbook of mathematical functions with formulas, graphs, and mathematical tables National Bureau of Standards Applied Mathematics Series, Washington, DC 1964 [] R Apéry, Irrationalité de ζ() et ζ(3), Astérique 61, (1979), [3] B C Berndt, Ramanujan s Notebooks, Part II (Springer-Verlag, New York, 1985) [4] JM Borwein and PB Borwein, A cubic counterpart of Jacobi s identity and the AGM, Trans Amer Math Soc, 33 (1991), [5] C D Aniello, On some inequalities for the Bernoulli numbers, Rend Circ Mat Palermo, 43, (1994), [6] K Dilcher, Sums of products of Bernoulli numbers, J Number Theory 60, no 1, (1996), 3 41 [7] S Gun, M R Murty, P Rath, Transcendental values of certain Eichler integrals, to appear in Bull London Math Soc [8] P Lakatos, On zeros of reciprocal polynomials, C R Math Rep Acad Sci Canada 4 (00), [9] P Lakatos and L Losonczi, Polynomials with all zeros on the unit circle Acta Math Hungar 15, no 4, (009), [10] M R Murty, C Smyth and R J Wang, Zeros of Ramanujan polynomials J of the Ramanujan Math Soc, 6 (011), no 1, [11] A van der Poorten, A proof that Euler missedapéry s proof of the irrationality of ζ(3) An informal report Math Intelligencer 1 (1978/79), no 4, [1] T Rivoal, La fonction zeta de Riemann prend une infinité de valeurs irrationelles aux entiers impairs, C R Acad Sci Paris Sér I Math 331, no 4, (000)

20 0 MATILDE LALÍN AND MATHEW ROGERS [13] A Schinzel, Self-inversive polynomials with all zeros on the unit circle, Ramanujan J, 9, (005), 19 3 [14] W Zudilin, Arithmetic of linear forms involving odd zeta values J Théor Nombres Bordeaux 16, no 1 (004), Département de mathématiques et de statistique, Université de Montréal Montreal, QC, Canada address: mlalin@dmsumontrealca Department of Mathematics, University of Illinois, Urbana, IL, USA address: mdrogers@illinoisedu