SECANT ZETA FUNCTIONS MATILDE LALÍN, FRANCIS RODRIGUE, AND MATHEW ROGERS arxiv:304.3922v2 [math.nt] 2 Jul 203 Abstract. Westudytheseries ψ sz := secnπzn s, andprovethatit converges under mild restrictions on z and s. The function possesses a modular transformation property, which allows us to evaluate ψ sz explicitly at certain quadratic irrational values of z. This supports our conjecture that π k ψ k j Q whenever k and j are positive integers with k even. We conclude with some speculations on Bernoulli numbers.. Introduction Let ζs denote the Riemann zeta function. It is well known that ζ2nπ 2n Q for n. Dirichlet L-functions and Clausen functions are modified versions of the Riemann zeta function, which also have nice properties at integer points [Lew9]. Berndt studied a third interesting modification of the Riemann zeta function, namely the cotangent zeta function [Ber76]: ξ s z := cotπnz n s.. He proved that. converges under mild restrictions on z and s, and he produced many explicit formulas for ξ k z, when z is a quadratic irrational, and k 3 is an odd integer. Consider the following examples: ξ 3 + 5 2 = π3 45 5, ξ 5 2 = π5 945 2. Berndt s work implies that j ξ k jπ k Q whenever j is a positive integer that is not a perfect square, and k 3 is odd. A natural extension of that work is to replace cotz with one of the functions {tanz, cscz, secz}. We can settle the tangent and cosecant cases via elementary trigonometric identities: tanπnz n s = ξ s z 2ξ s 2z, cscπnz n s = ξ s z/2 ξ s z, but it is more challenging to understand the secant zeta function: secπnz ψ s z := n s..2 Date: April 5, 203. 200 Mathematics Subject Classification. Primary 33E20; Secondary 33B30, L03. Key words and phrases. Secant zeta function, Bernoulli Numbers, Clausen functions, Riemann zeta function. This work has been partially supported by NSERC Discovery Grant 35542-2008 and FQRNT Subvention établissement de nouveaux chercheurs 44987. The work of FR has also been supported by a Bourse d été de premier cycle du ISM-CRM.
2 MATILDE LALÍN, FRANCIS RODRIGUE, AND MATHEW ROGERS The main goal of this paper is to prove formulas for specials values of ψ s z. In Section 2 we prove that the sum converges absolutely if z is an irrational algebraic number and s 2. In Section 4 we obtain results such as ψ 2 2 = π2 3, ψ 2 6 = 2π2 3. These types of formulas exist because ψ k z obeys a modular transformation which we prove in Section 3 see equation 3.8. Furthermore, based on numerical experiments, we conjecture: Conjecture. Assume that k and j are positive integers, and that k is even. Then ψ k jπ k Q. The results of Section 4 support this conjecture, even though there are still technical hurdles to constructing a complete proof. For instance, we prove that the conjecture holds for infinite subsequences of natural numbers. The rational numbers that appear are also interesting, and we speculate on their properties in the conclusion. 2. Convergence Since secπz has poles at the half-integers, it follows that ψ s z is only well-defined if nz Z+ 2 foranyinteger n. Thus,weexcluderationalpointswithevendenominatorsfrom the domain of ψ s z. If z = p/q with q odd, then ψ s p/q reduces to linear combinations of Hurwitz zeta functions, and.2 converges for s >. Convergence questions become more complicated if z is irrational. Irrationality guarantees that secπnz, but we still have to account for how often secπnz is large compared to n s. The Thue-Siegel-Roth Theorem gives that secπnz n +ε when z is algebraic and irrational, and this proves that.2 converges for s > 2. The case when s = 2 requires a more subtle argument. We use a theorem of Worley to show that the set of n s where secπnz is large is sparse enough to ensure that.2 converges. We are grateful to Florian Luca for providing this part of the proof. In summary, we have the following theorem: Theorem. The series in.2 converges absolutely in the following cases: When z = p/q with q odd and s >. 2 When z is algebraic irrational, and s > 2. 3 When z is algebraic irrational, and s = 2. Proof of Theorem, parts and 2. Let z be a rational number with odd denominator in reduced form. It is easy to see that the set of real numbers {secnπz} n N is finite. Let M = max n N secnπz. Then we have secπnz n s M n s. It follows easily from the Weierstrass M-test that.2 converges absolutely for s >. Now we prove the second part of the theorem. By elementary estimates secπnz = cscπnz /2 nz 2 k, 2. n where k n is the integer which minimizes nz 2 k n. Now appeal to the Thue-Siegel-Roth Theorem [Rot55]. In particular, for any algebraic irrational number α, and given ε > 0,
SECANT ZETA FUNCTIONS 3 there exists a constant Cα,ε, such that α p q > Cα,ε q 2+ε. 2.2 If we set α = z, then 2. becomes secπnz n z 2k n+ n+ε. 2n Therefore we have secnπz n s n s ε, and this implies that.2 converges absolutely for s > 2+ε. Since ε is arbitrarily small the result follows. Inorder to prove the thirdpart of Theorem, we requiresome background on continued fractions. Recall that any irrational number z can be represented as an infinite continued fraction z = [a 0 ;a,a 2,...], and the convergents are given by which satisfy [a 0 ;a,...,a l ] = p l, p l = a l p l +p l 2, 2.3 = a l + 2. 2.4 Convergents provide the best possible approximations to algebraic numbers among rational numbers with bounded denominators. In other words, if 0 < q <, then z p q > z p l. 2.5 In addition + > z p l > + +. 2.6 Now we state a weak version of a theorem due to Worley [Wor8, Thm. ]. Theorem 2 Worley. Let z be irrational, k 2, and p/q be a rational approximation to z in reduced form for which z p q < k q 2. Then either p/q is a convergent p l / to z, or where a and b are integers. p q = ap l +bp l a +b, a, b < 2k, Now we can complete the proof of Theorem. The following proof was kindly provided by Florian Luca.
4 MATILDE LALÍN, FRANCIS RODRIGUE, AND MATHEW ROGERS Proof of Theorem, part 3. Let k n be the integer which minimizes nz 2 k n. Let W z denote the set of integers where the quantity is large: { W z = n N : nz } 2 k n logn2. n Then secnπz n 2 secπz + nlogn 2, n W z n=2 which converges. Now assume that n W z. Then z +2k n 2n < logn2 n 2. Consider the convergents of z. Let l be such that 2n <. By Theorem 2 there are at most O log 4 solutions to z p < logn2 2n n 2 2.7 with p Z i.e. consider all values of a, b < 2k = 2log 2. From equations 2.5 and 2.6 we have z p 2n z p l + +. Combining this with equation 2.4 implies nz +2k n n 3 2 n2 + + Hence, if n W z, we find that q 3 l 8 a l+ ++ secnπz n 2 q2 l a l+ +2 ql 3 a l+ql 2 ql 3. q 3 l 8q 2 l a l+ +2. Combining the Thue-Siegel-Roth Theorem equation 2.2 with 2.6, implies that if z is algebraic > q z p l > Cz,ε. l+ q 2+ε l Thus + q +ε l. This allows us to place an upper bound on a l+ : Putting everything together gives the bound and as a result secnπz n 2 n W z a l+ + q ε l secnπz n 2 l= q 4ε l log 4 q ε l, l= q ε l.
SECANT ZETA FUNCTIONS 5 Since + = a l+ + +, we conclude that F l, where F l denotes the lth Fibonacci number. Since the Fibonacci numbers grow exponentially, we have ϕl, where ϕ = + 5. 5 2 Setting ϕ = ϕ ε >, we finally obtain secnπz n 2 n W z l= ϕ l = ϕ <. Thus, it follows that.2 converges absolutely when s = 2. 3. A modular transformation for ψ k z We begin by noting the following trivial properties of ψ s z: ψ s z =ψ s z, 3. ψ s z +2 =ψ s z, 3.2 2 s ψ s 2z =ψ s z+ψ s z +. 3.3 The main goal of this section is to prove that ψ k z also satisfies a modular transformation formula. We start from the partial fractions decompositions of secπx and cscπx, and then perform a convolution trick to obtain an expansion for secπx cscπxz equation 3.5. Differentiating with respect to x then leads to the transformation for ψ k z equation 3.8. This method is originally due to the third author, who used used it to rediscover the Newberger summation rule for Bessel functions [Rog05], [New82]: π J x yj x y sinπx = n= J 2 ny n+x. 3.4 Equation 3.4 follows from applying the convolution trick to partial fractions expansions for J x y and J x y [Lom68]. Lemma. Let χ 4 n denote the Legendre symbol modulo 4. Suppose that x and z are selected appropriately. Then πcscπzxsecπx = zx +8x χ 4 ncscπnz/2 n 2 4x 2 2zx Proof. Recall the classical partial fractions expansions [GR94]: πsecπx =6x 2 πcscπx x =2x k= secπn+/z n 2 z 2 x 2. 3.5 χ 4 n nn 2 4x 2, 3.6 k+ k 2 x 2. 3.7 We assume that z is irrational and algebraic, and x is selected so that the denominators in the sums are never zero.
6 MATILDE LALÍN, FRANCIS RODRIGUE, AND MATHEW ROGERS Both sums converge uniformly. Multiplying the formulas together, expanding via partial fractions, and rearranging the order of summation, we have πsecπx πcscπzx =32zx 3 k+ χ 4 n zx nn 2 4x 2 k 2 z 2 x 2 n k k+ χ 4 n nk 2 z 2 n 2 /4 z 2 n 2 /4 z 2 x 2 k 2 z 2 x 2 =8z 3 x 3 n k =32zx 3 +32zx 3 χ 4 n nn 2 4x 2 k= k= k+ k 2 z 2 x 2 k+ k 2 z 2 n 2 /4 χ 4 n nn 2 4k 2 /z 2 By 3.6 and 3.7 this becomes πsecπx πcscπzx =32x 3 χ 4 n zx n 2 n 2 4x 2 πcscπnz/2 2 nz +2πz 3 x 3 k= k+ k 2 k 2 z 2 x 2 secπk/z. This reduces to 3.5 after several additional applications of 3.6 and 3.7. We can split up the sums because all of the individual terms converge absolutely. For instance, we can prove that k+ secπk/z k 2 z 2 x 2, k= converges absolutely, by showing that the summand is secπk/z k 2, and then applying Theorem for appropriate choices of z. Theorem 3. Let E m denote the Euler numbers, and let B m denote the Bernoulli numbers. Suppose that k 2N. Then for appropriate choices of z: z z +z k ψ k z k ψ k +z z = πik k k [+z 2 m B m E k m m m z m ]. Proof. Recall the Taylor series expansions of cosecant and secant: πxcscπx = 2 m πix m 2B m, secπx = m! πix m E m. m!. 3.8 Expand both sides of 3.5 in a Taylor series with respect to x. Comparing coefficients yields the following identity: πi k k k 2 m B m E k m z m = 2 k χ 4 ncscπnz/2 m n k +z k ψ k +/z.
SECANT ZETA FUNCTIONS 7 Finally, let z +z and z z, and subtract the two results to recover 3.8. The cosecant sums vanish because cscπn+z/2 = cscπn z/2 whenever n is odd. We conclude this subsection with a conjecture on unimodular polynomials. We call a polynomial unimodular if its zeros all lie on the unit circle. We have observed numerically that the polynomials in 3.8 have all of their zeros on the vertical line Rez = 0. Since the linear fractional transformation z = x/+x maps the vertical line to the unit circle, we arrive at the following conjecture: Conjecture 2. We conjecture that the polynomial k k 2 m 2 m 2B m E k m x x m +x k m, 3.9 m is unimodular when k 2N. This new family of polynomials is closely related to the unimodular polynomials introduced in [LR2], [LS3], and [HZ3]. 4. Special values of ψ k z Throughout this section we assume that k is a positive even integer. Let SL 2 Z denote the set of 2 2 integer-valued matrices with determinant equal to one, and let PSL 2 z = SL 2 Z/ I, I. Recall that PSL 2 Z can be identified with the set of linear fractional transformations. The usual group action is a b c d az +b z = cz +d. It is easy to see that multiplying matrices is equivalent to performing compositions of linear fractional transformations. Consider the following matrices in PSL 2 Z: 2 0 A =, B =. 0 2 Equations 3.2 and 3.8 are equivalent to ψ k Az =ψ k z, 4. ψ k Bz =2z + k ψ k z 4.2 k + πik 2 m k B m E k m z + k m[ 2z + m k 2z + k]. m Every matrix C A,B has a factorization of the form C = A j B j 2 A j 3..., so equations 4. and 4.2 together imply that there exists a z-linear relation between ψ k z, ψ k Cz, and π k. Now we outline a strategy to obtain exact evaluations of ψ k z. First select a matrix C, and then find the linear relation between ψ k z, ψ k Cz, and π k. Next choose z so that ψ k z = ψ k Cz. For example, if z = 2j + 2j2j + in 4.2, then Bz = z 4j
8 MATILDE LALÍN, FRANCIS RODRIGUE, AND MATHEW ROGERS whenever j is a non-zero integer. With some work the equation collapses to k ψ k 2j2j + = πik 2 m k B m E k m m m + 2j 2j+ 2j 2j+ k + 2j 2j+ 2j 2j+ m k, 4.3 which holds for j Z \ {0}. Similarly, if we take z = 2j + + 2j2j +, then we arrive at ψ k + k 2j2j + = πik 2 m k B m E k m m m m + 2j+ 2j 2j+ 4.4 2j k k, + 2j+ 2j 2j+ 2j which also holds for j Z \ {0}. The right-hand sides of 4.3 and 4.4 are invariant under the Galois action x x, so both formulas are rational with respect to j. Specializing 4.3 at k = 2 and k = 4 yields ψ 2 2j2j + =3j + π2 6, 75j 2 +46j +6 π 4 ψ 4 2j2j + = 8j +3 80. These results support Conjecture. Further evidence for the conjecture is provided by combining 4.3, 4.4, and 3.3, to obtain identities like ψ 2 8j2j + = π2 6, 4.5 whenever j Z. Ingeneral, it seemstobequitedifficulttoevaluate ψ k jforarbitrarypositiveintegers j. This is due to the fact that 4. and 4.2 restrict the available matrices to a subgroup of PSL 2 Z. If there exists a matrix C A,B which satisfies C j = j, then we can always evaluate ψ k j. We can construct candidate matrices by solving Pell s equation: If we choose C to be given by X 2 jy 2 =. X j Y C = Y X then it follows that C j = X j+jy Y j+x = j, and detc = X 2 jy 2 =. Pell s equation has infinitely many integral solutions when j [HW38], so there are infinitely many choices of C. The main difficulty is to select appropriate values of X and Y so that C factors into products of A s and B s. It is not clear if such a selection is always possible. Notice that A,B Γ 0 2 and so it follows that PSL 2 Z A,B.,
SECANT ZETA FUNCTIONS 9 We conclude this subsection by noting that ψ 2 z = 0 for infinitely many irrational values of z. Setting n = 3j in Proposition below, yields 26j ψ 2 2 = 0, 4.6 3 for any non-zero integer j. It is worth emphasizing that ψ k z is highly discontinuous with respect to z, so these types of results are not surprising. Proposition. Suppose that j and n are integers, and n 0. Then 2j2jn+ ψ 2 = + 3j π 2 n n 6. 4.7 Proof. Setting k = 2 in 4.2 yields Iterating 4.8 gives ψ 2 Bz = 2z + ψ 2z+ z3z2 +4z +2 π 2 2z + 2 6. 4.8 ψ 2 B n z = 2nz + ψ 2z+ nz 3z 2 +4nz +2 2nz + 2 π 2 6. 4.9 The derivation of 4.9 is best accomplished with the aid of a computer algebra system such as Mathematica, because significant telescoping occurs on the right. Now consider the matrix C = A j B n A j 4jn+ 4j2jn + =, 2n 4jn+ and notice that Cz 0 = z 0, where 2j2jn+ z 0 =. n Thus by 4.9 we have ψ 2 z 0 =ψ 2 Cz 0 =ψ 2 B n z 0 +2j = 2nz 0 +2j+ ψ 2z 0 + nz 0 +2j 3z0 +2j 2 +4nz 0 +2j+2 π 2 2nz 0 +2j+ 2 6. We complete the proof by solving for ψ 2 z 0 and simplifying. 5. Speculations and Conclusion Assume that k is a positive even integer. Euler gave the following expression for Bernoulli numbers: B k = 2 2πi kζk. 5. Bernoulli numbers are interesting combinatorial objects, so it is natural to ask if the rational partofψ k j also hasinteresting properties. Thisis anobvious question because ψ k 0 = ζk. For instance, the von Staudt-Clausen theorem gives a complete description of the denominators of Bernoulli numbers: B k = p k + Integer. 5.2 p
0 MATILDE LALÍN, FRANCIS RODRIGUE, AND MATHEW ROGERS Is there an analogue of 5.2 for the rational part of ψ k j? As an example, consider ψ k 6 which can be calculated from equation 4.3. In order to eliminate some trivial factors, we define k k 3+ 6 3 6 β k := 6 πi kψ k 6. By 4.3 we have k β k = 2 m m m k 3+ 6 3 6 B m E k m 3 k m. m 6 The first few values are β 2 = 8/3, β 4 = 508/5, and β 6 = 64896/7, which all have denominators divisible only by primes where p k as hoped for. The first instance where this fails is β 20. The denominator of β 20 equals 5 7, but 7 = 6 does not divide 20. Acknowledgements. The authors wish to thank Shabnam Ahktari for helpful discussions and Florian Luca for providing the proof of part 3 of Theorem. Finally, the authors are grateful to the referee for comments and suggestions that greatly improved the exposition of this article. References [Ber76] B. C. Berndt, Dedekind sums and a paper of G. H. Hardy, J. London Math. Soc. 2 3 976, 29 37. [GR94] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products, Academic Press 994 [HW38] G. H. Hardy and E. M. Wright, An introduction to the theory of numbers, Clarendon Press 938. [HZ3] H. K. Hwang V. Zacharovas, Limit laws of the coefficients of polynomials with only unit roots, Preprint 203, arxiv:30.202. [LR2] M. N. Lalín and M. D. Rogers, Variations of the Ramanujan polynomials and remarks on ζ2j + /π 2j+, to appear in Funct. Approx. Comment. Math.. [LS3] M. N. Lalín and C. J. Smyth, Unimodularity of roots of self-inversive polynomials, Acta Math. Hungar. 38 203, no. -2, 85 0. [Lew9] L. Lewin, Structural properties of polylogarithms, Mathematical Surveys and Monographs, 37. American Mathematical Society, Providence, RI, 99. xviii+42 [Lom68] E. C. Lommel, Studien uber die Bessel schen Functionen, Leipzig 868. [New82] B. S. Newberger, New sum rule for products of Bessel functions with application to plasma physics, J. Math. Phys. 23 982, 278. [Rog05] M. D. Rogers, Partial fractions expansions and identities for products of Bessel functions, J. Math. Phys. 46 2005, 043509. [Rot55] K. F. Roth, Rational approximations to algebraic numbers, Mathematika 2 955, 20; corrigendum, 68. [Wor8] R. T. Worley, Estimating α p/q. J. Austral. Math. Soc. Ser. A 3 98, no. 2, 202 206. Department of Mathematics and Statistics, University of Montreal, Montreal, Canada E-mail address: mlalin@dms.umontreal.ca Department of Mathematics and Statistics, University of Montreal, Montreal, Canada E-mail address: rodriguefrancis@gmail.com Department of Mathematics and Statistics, University of Montreal, Montreal, Canada E-mail address: mathewrogers@gmail.com