arxiv: v1 [math.nt] 12 Sep 2018

Transcription

1 arxiv: v [mathnt] Sep 08 CONNECTION COEFFICIENTS FOR HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: A RANDOM WALK APPROACH LIN JIU AND CHRISTOPHE VIGNAT Abstract We consider the use of random walks as an approach to obtain connection coefficients for higher-order Bernoulli and Euler polynomials In particular, we consider the cases of a -dimensional linear reflected Brownian motion and of a 3-dimensional Bessel process Considering the successive hitting times of two, three, and four fixed levels by these random walks yields non-trivial identities that involve higher-order Bernoulli and Euler polynomials Introduction Different types of random walks have been studied in the literature, together with their connections to different fields of mathematics and physics; for a modern introduction, see for example [6] In this paper, we focus on two specific random walks: the -dimensional linear reflected Brownian motion and the 3-dimensional Bessel process It would seem that there is no relation between these processes and special functions such as Bernoulli and Euler polynomials that appear mostly in number theory and combinatorics However, we will show how the study of the hitting times of these two processes allows us to derive non-elementary identities for higher-order Bernoulli and Euler polynomials, denoted by B n p x and E n p x, respectively These polynomials are defined through their generating functions as p t p e t e xt B n p tn x and n! e t e xt + n0 n0 E n p tn x n! The special case p yields the usual Bernoulli and Euler polynomials: B n x B n x and E n x E n x; in addition, Bernoulli and Euler numbers are the special evaluations B n B n 0 and E n n E n / One of the key tools of this article is the interpretation of the polynomialsb n p x and E n p x as probabilistic moments of certain random variables that are related to the hitting times of some random walks, which will be introduced in Section This study arises from early work on the higher-order Euler polynomials In a previous article [5, eq 38, p 78], we obtained the following expansion for the usual Euler polynomials as a linear combination of higher-order Euler polynomials: for any positive integer N, E n x N n p N l N l E n l +Nx, ln

2 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK where the positive coefficients p N l are defined by the generating function T N /t p N l t l, via the N-th Chebychev polynomial T N t These coefficients p N l also appear as transition probabilities in the context of a random walk over a finite number of sites [5, Note 48, p 787] This characterization provides an alternate interpretation of in terms of a Markov process For the two random walks considered in this paper, we shall focus on their consecutive hitting times: namely, for a process starting from the origin at time t 0, we consider the first epochs - called hitting times - at which the process reaches successive levels 0 a 0 < a < a < < a N In particular, we study some special cases for N, 3, and 4 From the moment generating functions of these hitting times, we derive identities involving B n p x and E n p x which will be presented in Section 3 and Section 4 l0 Probabilistic Preliminaries and Principle The symbols B, E and U We shall fully make use of the classical umbral symbols B, E, and U, defined as follows - for an introduction to the classical umbral calculus, see for example [4, 9] The Bernoulli symbol B satisfies the evaluation rule B+x n B n x Equivalently, this symbol B can be interpreted as a random variable [3, Thm 3, p 384]: if L B is a random variable distributed according to the square hyperbolic secant density p B t πsech πt/, then, for any suitable function f, f x+b E [ f x+il B ] R f x+it p B tdt In particular, choosing f x x n produces the Bernoulli polynomial, so that B coincides with the random variable il B / More generally, the p-th order Bernoulli polynomial can be expressed as B p n x x+b + +B p n, where B i p i is a set of p independent umbral symbols or random variables that satisfy the two following evaluation rules: : if B and B are independent symbols then : if B B, B n Bm B nb m ; B n Bm Bn+m B n+m For simplicity, we shall denote the sum of p independent Bernoulli symbols as B p B + +B p From the generating function, we deduce that e Bt t e t, etb+ t sinht, and p +p t etb sinht p

3 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK 3 Similarly, we let L E be a random variable distributed according to the hyperbolic secant density p E t sechπt, and define the Euler symbol E by [ f x+e E f x+il E ] f x+it p E tdt Then, denoting E p E + +E p, the sum of p independent Euler symbols, the Euler polynomial of order p is expressed as Also, we have E p n x [x+e p] n e te e t +, ete+ secht, and From the generating function R e te p +p sech p t e Bt t e t t e t e t + etb+e, we deduce that the symbols B and E satisfy the rule B B +E, in the sense that for any suitable function, f x+b f x+b+e The last useful symbol is the uniform symbol U with the evaluation rule U n n+ It can be interpreted as a random variable U uniformly distributed over the interval [0,], and for any suitable function f, f x+u E[f x+u] 0 f x+tdt Thus, we have in terms of generating functions e tu t n n+! et, e tu sinht p, and e p sinht tup t t t n0 where U p U + +U p denotes the sum of p independent uniform symbols An important link between the Bernoulli symbol B and the uniform symbol U is deduced from the identity e tu+b e tu e tb t e t et ; t this shows that, for any suitable function f, 3 f x+b +U f x, so that the actions of these two symbols cancel each other In what follows, we will use independent copies of Bernoulli, Euler and uniform symbols that, for example in the uniform case, they will be denoted by U, U, and U p,u p, in order to be distinguished

4 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK 4 Level sites with loop and loops The general setting of this paper is thus a random walk starting from the origin in R or R 3 and hitting some defined levels -either some points on the line or some sphere in the 3-dimensional space, called sites Looping back and forth, each random walk can reach each site multiple times, and we are interested in the first time the process reaches each of these sites Before we derive identities, we shall first consider the contribution of loops to the hitting times of these sites Consider one possible loop between sites a and b, with a < b Let φ a b be the moment generating function of the hitting time of site b staring from site a, without hitting any other site if any; and also let φ b a be the counterpart from b to a of φ a b Let us further denote by I a,b φ a b φ b a ; the random walk can loop an arbitrary number of times k 0 between sites a and b without visiting any other sites so that the overall contribution of these visits to the generating function is 4 I k a,b I a,b Next, consider one possible loop between a and b, and another between c and d, with a < b c < d Using similar notations as above, we define I a,b φ a b φ b a and I c,d φ c d φ d c The possible contributions are as follows: k loops between sites a and b followed by l loops between sites c and d, with k,l 0,,, contributing kl 0 I k a,bi l c,d I a,b I c,d ; k loops between sites a and b followed by l loops between sites c and d, then followed by k loops between sites a and b and finally l loops between sites c and d, with k, l nonnegative and k, l positive, contributing k l 0,k,l I k a,b Il c,d Ik a,b Il cd I a,b I c,d I a,b I c,d ; the general term will consist of k loops between sites a and b followed by l loops between sites c and d and so on followed by k n loops between sites a and b and l n loops between sites c and d, with k, l n nonnegative and the other indices being positive, overall contributing to the generating function I n a,b In c,d I a,b n I c,d n Therefore, the overall contribution of all possible loops over sites a and b followed by loops over sites c and d is I n a,b 5 In c,d I a,b n I c,d n I a,b +I c,d I a,b +I c,d k n

5 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK 5 3 One-Dimensional Reflected Brownian Motion 3 Introduction Consider the -dimensional reflected Brownian motion on R + For simplicity, we let φ r s z be the generating function of the hitting time of site s starting from site r; φ r s tz be the generating function of the hitting time of site s starting from site r without reaching site t In the case of three consecutive sites a, b and c with 0 < a < b < c, the generating functions of the corresponding hitting times can be found in [, p 98 and p 355]: with w z, φ a b z coshaw coshbw sinhc bw φ b a cz sinhc aw sinhb aw φ b c az sinhc aw Let us now label N + sites on the real positive line as 0 a 0 < a < a < < a N 3 The case of 3 sites We consider here 3 sites 0 a 0 < a < a ; the study of hitting times yields the following identity Theorem Let 0 < a < a, then we have, for any positive integer n, the expansion of Euler polynomials as a convex combination of higher-order Bernoulli polynomials as x x E n E n 34 n+ where the coefficients + 3 a a a a a n a n a n a + p k B k+ n p k a a k a a x + a + k, 4a 4a are the probability weights of a geometric distribution with parameter a /a Proof From formulas above, we have φ 0 a z secha w, φ 0 a z secha w, φ a a 0z sinha w sinha w, and φ a 0 z sinha a w a sinha w

6 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK 6 The process includes one possible loop between sites 0 and a, so that by 4, we have 35 φ 0 a z φ 0 a zφ a a 0z φ0 a zφ a 0 z k a secha w sinha w sinha w sinha w cosha wsinha w secha w sinha k a w sinha w sinha aw cosha wsinha w sinha w cosha wsinha w sinha a w, which coincides with φ 0 a z secha w, since sinha a w sinha wcosha w cosha wsinha w Meanwhile, we transform 35 into a form involving Bernoulli, uniform and Euler symbols as follows: e awe+ a a e awu e awb+ e awek+ +k+ e awu +awb+ k a a e k a awuk e a wb k +k a p k e awek+ +k++a a wu k k+a wb k +k Multiplying both sides by e xw and comparing coefficients of w n yield x+a E +a n [ p k x+a U +a B++a E k+ +k + +a a U k k +a ] n B k +k n, p k [x+a a B+a E k+ +a B k +a +a k] where, in the last step, we have used the fact that B and U cancel each other according to 3 Now, the substitution x x+a a U yields x+a a U +a E +a n n p k x+a E k+ +a B k+ +a +a k

7 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK 7 Here, the left-hand side can be computed as [ n a n E + x + a + a ] n x U n a n a E n + a + x n a n E n + 0 a + n a n x [E n a + 3 a a a U a a E n a a a n+ while the right-hand side is, by, p k x+a E k+ +a B k+ +a +a k n p k x+4a B k+ +a +a k x p k 4 n a n B n k+ + a + k 4a 4a Further simplification completes the proof n t x a + Remark For uniformly spaced levels a a / a/, identity 34 collapses to the trivial identity The case of 4 sites We consider now the case of four sites 0 a 0 < a < a < a 3 with two possible loops, one between a 0 and a and one between a and a For the sake of simplicity, we consider only the case where all sites uniformly distributed, with a i i for i 0,,, 3, and obtain the following identity Theorem For any positive integer n, the Euler polynomial of degree n can be expressed as a linear combination of high-order Euler polynomials of the same degree as 3 k n 36 E n x 3x+k 4 k+ek+3 n Proof By 5 and the basic identity sinhw sinhwcoshw, we have φ 0 3 z φ 0 zφ 0zφ 3 z k φ 0 zφ 0 z+φ 0zφ 3z sinhw sechw sinhw k 3 sech k+3 w, 4 4 sechw sinhw sinhw + sinhw sinhw which is equal to sech3w φ 0 3 z Therefore, we deduce e 3wE+ 3 k k+3 +k+3, 4 k+ewe dt ] k sinhw sinhw ;

8 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK 8 namely, e 6wE 3 k k+3 +wk 4 k+ewe Multiplying both sides by e xw and comparing coefficients of w n, we obtain x+6e n where the left-hand side is x+6e n 6 n E + x 6 3 k 4 k+ x+k +E k+3 n, n 6 n E n x 6 ; and the right-hand side is 3 k 4 k+ x+k+e k+3 n 3 k x 4 k+n +k +Ek+3 n 3 k x 4 k+n E n k+3 +k Simplification completes the proof In the general case of 4 sites that are not uniformly spaced, we obtain a much more complicated result as follows: Theorem 3 For any 0 < a < a < a 3 and arbitrary positive integer n, we have x E n + a 3 with the coefficients and k l0 q k,l a a 3 n x E n l +U k l + a a B k l + r k,l a a q k,l a + a a a B + a 3 a, k a a l+ a k l+ a 3 a k l l a k+ a 3 a k l+, r k,l a 3 +k la +3l k +a a B + a a a U l

9 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK 9 Proof Apply 5 to obtain φ 0 a3 z φ 0 a zφ a a 0zφ a a 3 a z k φ 0 a zφ a 0 z+φ a a a 0zφ a a z a3 secha w sinha w sinha w sinha a w sinha 3 a w secha w sinha a w + sinha w sinha w sinha w sinha 3 a w sinha 3 a w k k sinh k+ sech l+ a wsinh l+ a a w a w l l0 sinhk l+ a wsinh k l a 3 a w sinh k l+ a 3 a w Now, we multiply both sides by e xw and compare coefficients of w n, with the notation and to obtain x+a 3 E +a 3 n q k,l k a a l+ a k l+ a 3 a k l l a k+ a 3 a k l+ r k,l a 3 +k la +3l k +a, k l0 q k,l [ x+a a B+a 3 a B +a E l +a a U l +a U k l +a a B k l +r k,l ] n 4 Bessel Process in R 3 4 Introduction We consider now a Bessel process R 3 t in R 3 Using similar notations φ a b z and φ a b cz as in the previous section and considering sites a < b < c that are now concentric spheres of radii a, b, and c, we will need the following formulas from [, pp ]: Remark One can easily check that φ a b z bsinhaw asinhbw, csinhb aw φ b c az bsinhc aw, asinhc bw φ b a cz csinhc aw φ a 0 a z 0 sinha a w a sinha 0w 0, k

10 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK 0 so that the first possible loop is now between a and a, different from the - dimensional case, where the first possible loop was between the origin 0 and the first site a 4 The case N 3 Theorem 4 In the case of three concentric spheres of arbitrary radii 0 < a < a < a 3, we deduce the following expression of a Bernoulli polynomial of degree n in terms of higher-order Bernoulli polynomials B n x+a3 a 3 with and a a 3 n p k B k+ n + a 3 a B k+ a ] n [ x a +β k + a a a U + a 3 a a U k + a a U k p k αα k, β k a 3 +ka ka a 3 α a 3 a a a 3 a a The special case of uniformly spaced radii is as follows Corollary 5 In the case of three spheres of radii a,a, and a 3 3, we obtain the identity 44 3 n+ [ x B n+ n x B n ] 3 4 k E n k+ 4 x+3+k The specialization to the case x 0 and n odd of this identity is as follows Corollary 6 The even Bernoulli number B m can be expressed as a convex combination of high-order Euler polynomials as follows m k B m m 3 m E k+ m k Proof Take x 0 and n m in 44 By Entries 447 and 443 of [8], the left-hand side becomes 3 m [ ] [ 5 B m B m 3m m 3 m B m + m ] B m m 6 m 3m m 3 m B m + m 3 m 3 m B m ; 4m while the right-hand side is k E k+ m k + 3 And we obtain, after simplification, the desired result

11 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK Proof of Theorem 4 Remarking that lim x 0 [sinhxt/x] t, from 4, 4 and 43, we have so that φ a0 a 3 z φ a0 a z aw sinha, w φ a0 a 3 z a3w sinha, 3w a sinha φ a a a0 z w a sinha, w a sinha φ a a a3 z 3 a w a sinha, 3 a w a3 sinha φ a a 3 a z a w a sinha, 3 a w φ a0 a zφ a a z φ a0 a a zφ a3 a a z k φa a a0 3 a z a w a sinha w a3 sinha a w sinha w a sinha w sinha3 a w sinha 3 a w sinha w sinha w a 3 wsinha a w a sinha 3 a w k sinh k a 3 a wsinh k a w sinh k+ a 3 a wsinh k+ a w Simple algebra shows that this expression coincides with a 3 w φ a0 a 3 sinha 3 w Then, we obtain where e a3wb+ e wa au { r k e w[a3 auk k+a U k k] e w[a3 abk+ +k++a B k+ +k+] r k a 3a a a 3 a k a k a 3 a k+ a k+ Similarly, we multiply both sides by e wx and look at the coefficients of w n For simplicity, we let s k a 3 +ka ka, so that a 3 n B n x+a3 a 3 The special case a i i, with r k [ x+a 3 B+a 3 +s k +a a U +a 3 a U k +a U k +a 3 a B k+ +a B k+ ] n β k s k 3+k and ρ k r k aii aii 3 4 } k, 4

12 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK produces 6 n B n x+3 6 k 0 k 0 k 0 k 0 k 0 ρ k x+3+k +U +U k +U k +4B k+ +4B k+ n ρ k x+3+k +U k+ +4B k+ n ρ k x+3+k +U k+ +B k+ +E k+ n ρ k x+3+k +B+E k+ n x+3+k ρ k n E n k+ +B, or x 3 n B n 6 + x+3+k ρ k E n k+ +B k 0 Replacing x by x+u completes the proof 43 The case N 4 We have now two possible loops among the sites with radii a, a 3 and a 4 Following the notation in Subsection, we have and similarly, I a,a φ a a a0 zφ a a a3 z sinha wsinha 3 a w sinha wsinha 3 a w I a,a 3 φ a a 3 a zφ a 3 a a4 z sinha a wsinha 4 a 3 w sinha 3 a wsinha 4 a w In order to obtain a simple expression for I a,a + I a,a 3, we further assume that the levels are uniformly distributed, ie, a i i, for i 0,,, 3, 4, so that This produces the following identity I a,a +I a,a 3 sinh w sinh w sech w Theorem 7 For any positive integer n, the Bernoulli polynomial of degree n can be expanded as the convex combination of higher-order Euler polynomials of the same degree as x+4 45 B n 6 3 n Proof From kek+ n x+k +3 φ 0 4 z φ 0 zφ 0zφ 3 zφ 3 4 z I a,a +I a,a 3 k,

13 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK 3 namely, 4w sinh4w we deduce that w sinhw sinhw sinhw w sinhw sech k+ w k, e 4wB+ e wb + 3sinhw sinhw 4sinhw 3sinhw k+ +k+ kewe After multiplying by e xw, coefficients of w n on both sides yield [ x+b x+8b+4 n ++E k+ +k + ] n k Apply the substitution x x+u, to get the left-hand side as n x+6b+4 n x+4 x+4 6 n +B 6 n B n, 6 6 and the right-hand side as [ x+e k+ +k +3 ] n k n Further simplification completes the proof n [ ] E k+ + x+k+3 n E k+ n k x+k+3 k sech k w k Remark Note the analogy between 45 and 36 whereas they are obtained from two different setups 5 Conclusion We have shown how the setup of random processes allows to obtain non trivial identities among higher-order Bernoulli and Euler polynomials The underlying principle of this approach is that these special functions appear naturally in the generating functions of the hitting times of these random processes Several remarks are in order at this point: - the identities obtained from this approach are not of the usual, convolutional type - see for example [] for such identities They are rather connection-type identities between the usual Bernoulli and Euler polynomials and their higher-order counterparts; - these inherently involve a mixture of higher-order Bernoulli and Euler polynomials; - the interest of this approach is that each term in the obtained decomposition can be physically related to a relevant object, namely one loop in a trajectory of a random process - this work should be considered as a first approach only to a more general project in which the richness of the possible setups for random walks is expected to generate a number of non-trivial identities about more general special functions

14 HIGHER-ORDER BERNOULLI AND EULER POLYNOMIALS: RANDOM WALK 4 Acknowledgement The authors would like to thank Prof P Salminen and Prof B Q Ta for their valuable suggestions on this work They also thank Prof K Dilcher for the organization of the Eighteenth International Conference on Fibonacci Numbers and their Applications References [] A N Borodin and P Salminen, Handbook of Brownian Motion - Facts and Formulae, Birkhauser, 05 [] K Dilcher, Sums of products of Bernoulli numbers, J Number Theory, , 3 4 [3] A Dixit, V H Moll, and C Vignat, The Zagier modification of Bernoulli numbers and a polynomial extension Part I, Ramanujan J 33 04, [4] I Gessel, Applications of the classical umbral calculus, Algebra Universalis, , [5] L Jiu, V H Moll, and C Vignat, Identities for generalized Euler polynoimials, Integral Transforms Spec Funct 5 04, [6] G Lawler, V Limic, Random Walk: a Modern Introduction, Cambridge University Press, 00 [7] L B Klebanov, A V Kakosyan, S T Rachev, and G Temnov, On a class of distributions stable under random summation, J Appl Probab 49 0, [8] F W J Olver, D W Lozier, R F Boisvert and C W Clark, NIST Handbook of Mathematical Functions, Cambridge University Press, 00 [9] S Roman, The Umbral Calculus, Academic Press, 984 Department of Mathematics and Statistics, Dalhousie University, Halifax, NS, B3H 4R, Canada address: LinJiu@dalca Department of Mathematics, Tulane University, New Orleans, LA, 708, USA address: cvignat@tulaneedu