arxiv: v1 [math.ca] 1 Jan 2019

Transcription

1 An expnsion for the sum of product of n exponentil nd Bessel function R. B. Pris rxiv:90.004v [mth.ca] Jn 09 Division of Computing nd Mthemtics, Aerty University, Dundee DD HG, UK Astrct We exmine convergent representtions for the sum of decying exponentil nd Bessel function in the form e n Jν(n, nν ( where J ν(x is the Bessel function of the first kind of order ν > nd, re positive prmeters. By mens of doule Mellin-Brnes integrl representtion we otin convergent symptotic expnsion tht enles the evlution of this sum in the limit 0 with < π fixed. A similr result is derived for the sum when the Bessel function is replced y the modified Bessel function K ν(x. The lternting versions of these sums re lso mentioned. Mthemtics Suject Clssifiction: 33C05, 33C0, 33C0, 4A30, 4A60 Keywords: Bessel functions, Mellin-Brnes integrl, symptotic expnsion. Introduction We consider the sums e n S J,ν (, = ( J e n ν(n, nd S K,ν (, = nν ( K ν(n, (. nν where J ν (x is the Bessel function of the first kind nd K ν (x the modified Bessel function. The prmeters, re ssumed to e positive nd the order ν > for S J,ν(, nd ν 0 for S K,ν (,. Convergent expnsions re derived for these sums tht involve the polylogrithm function Li s (x of negtive integer order s when > nd the Riemnn zet function when ; the ltter cse requires, (0,π to secure convergence. The evlution of the lternting versions of (. follows in strightforwrd mnner in terms of these sums with prmeters, nd,. Our min interest in the sums in (. is the limit s 0 when convergence of the first sum ecomes slow. When = 0, the sum S J,ν (0, hs een considered y Tri cković et l. in [4], where pproches using Poisson s summtion formul nd Bessel s integrl were employed to derive convergent expnsions. The sums S J,ν (0, nd S K,ν (0, hve een discussed in [] using Mellin trnsform pproch. With the presence of the exponentil fctor in the ove sums it is found tht doule Mellin-Brnes integrl is required to derive n expnsion s 0. The tretment of such integrls hs een discussed in [3, Chpter 7] in the context of multi-dimensionl Lplce-type integrls. E-mil ddress: r.pris@erty.c.uk

2 R.B. Pris. A series representtion for S J,ν (, vlid when We consider the sum S J,ν (, = e n ( nν J ν(n ( > 0, > 0, (. where J ν (x denotes the usul Bessel function of rel order ν >. In the cse = 0, the sum hs the vlue [, Theorem, (.9] S J,ν (0, = Γ(ν + Γ(+ν (ν >. (. Expnding the Bessel function s series, we hve S J,ν (, = e n ( k ( nk k!γ(+ν +k = ( k ( k k!γ(+ν +k e n n k upon interchnging the order of summtion. If we introduce the polylogrithm function Li s (x defined y Li s (z = n s z n [, (5..0], we then otin S J,ν (, = ( k ( k k!γ(+ν +k Li k(e. (.3 The polylogrithm function of negtive integer order n, n = 0,,,... is { Li n (e = n! n + {( +πki k { = n! n + ( k n+ +( n+ } πki }, +4π k n+ cos[(n+φ k ] with φ k := rctn(πk/. For lrge n with fixed we hve Li n (e n! n (n. The lte terms in the sum (.3 consequently possess the ehviour ( k ( k Γ(k + = ( k Γ(k + k!γ(+ν +k Γ(+ν +k ( k ( k = ( O( k k ν / s k, so tht the sum (.3 converges when >, nd when = provided ν >. Clcultions with Mthemtic confirm the vlidity of this sttement, where it should e noted tht Li s (z is uilt-in lirry function. 3. A series representtion of S J,ν (, for 0 To del with the cse <, nd in prticulr to otin n expnsion for S J,ν (, s 0, we employ the Chen-Mellin integrl [3, p. 89] e x = πi c+ i c i Γ(sx s ds ( rg x < π (3.

3 A modified Bessel function sum 3 nd the Mellin-Brnes integrl [, (0.9.] ( x ν J ν (x = c+ i Γ(s πi c i Γ(+ν s ( x s ds (x > 0, ν > 0, (3. where c > 0 in oth cses so tht the integrtion pth lies to the right of the poles of the integrnd situted t s = 0,,,.... Comintion of these two results leds to the doule Mellin-Brnes integrl representtion for S J,ν (, given y, provided ν > 0, S J,ν (, = ( πi i i Γ(sΓ(t Γ(+ν +s t ( s n n s t dsdt = ( πi i i Γ(sΓ(t Γ(+ν +s t ( s ζ(s+tdsdt, (3.3 where ζ(s is the Riemnn zet function. The integrtion contours re indented to pss to the right of the pole of the zet function nd those of Γ(s nd Γ(t. Let us first consider the poles t s = n, t = m (m,n = 0,,,... tht result from the gmm functions in the numertor of (3.3. Displcement of the integrtion pths to the left produces the forml series T = = ( m+n m!n! ζ(0 Γ(+ν m ( n Γ(+ν +n ζ( n m ( n m+ ( n (m+!n! Γ(+ν +n ζ( m n, where even vlues of m (prt from m = 0 when n = 0 do not contriute to the doule sum on ccount of the trivil zeros of ζ(s t s =, 4,.... Using the functionl reltion for ζ(s [, (5.4.] ζ(s = s π s ζ( sγ( ssin πs, the vlue ζ(0 = nd the dupliction formul Γ(z = z π / Γ(zΓ(z +, we find tht T = Γ(+ν + π where = Γ(+ν + π ( n cosπ(m+n (m+!n! ( m ( n Γ(m+n+ ζ(m+n+ Γ(n++ν π π ( m+ ( n ( m A m,n, (3.4 π π A m,n := Γ(m+n+Γ(m+n+ 3 m!n!γ(m+ 3 ζ(m+n+. (3.5 Γ(n++ν Appliction of Stirling s formul for the gmm function, together with the fct tht ζ(m + n + = O( for lrge m, n, shows tht A m,n hs only lgeric growth in m nd n s m,n. Consequently the doule sum in (3.4 converges when < π nd < π. We now del with the pole of the zet function when s + t = with unit residue. Setting t = s, we find the contriution i ( s Γ(sΓ( s ds πi i Γ(+ν +s If the integtion pth in (3. is ent ck into loop tht encloses the poles of Γ(s with endpoints t infinity in R(s < 0, then the integrl holds without restriction on ν; see [3, p. 5].

4 4 R.B. Pris where the integrtion pth is indented to seprte the poles of Γ(s nd Γ( s. Evlution of the residues of Γ( s t s = + k, k = 0,,,... leds to the contriution T = ( k Γ( k + ( k k!γ( +ν k. (3.6 The terms in this sum ehve like ( k Γ( k + ν ( k ( ( k πγ( k + cosπ(ν k = ( k O k ν / (k, so tht (3.6 converges for <, nd for = provided ν >. By seprting terms with even nd odd k, we cn express T in terms of well-known specil functions to yield T = Γ(ν + (+ ν Γ(ν F (, ν; 3 ;. (3.7 Alterntively, the hypergeometric function cn e written in different form using Euler s trnsformtion [, (5.8.] to yield T = Γ(ν + (+ To summrise, we hve the following expnsion. Theorem. Let ν > S J,ν (, = Γ(ν + ν ( + Γ(ν F (, +ν; 3 ; +. nd. Then the following convergent expnsion holds (+ + π ν Γ(+ν Γ(ν F (, ν; ( ( m A m,n π provided < π nd < π. The coefficients A m,n re defined in (3.5,. 3 ; m ( n, (3.8 π When = 0 it is seen tht (3.8 correctly reduces to (.. The ove expnsion cn e rrnged in form suitle for symptotic clcultions in the limit 0. We note tht the series representtion of the hypergeometric function in (3.7 is lredy in the form of convergent symptotic expnsion s 0. We find fter routine mnipultion tht T = Γ(+ν + ( k+ ( k C k, π π where Consequently we otin C k = n=0 ( n A k,n ( < π. (3.9 π Theorem. Let ν > nd < π. Then the following (convergent symptotic expnsion holds for 0 S J,ν (, = Γ(ν + (+ ν Γ(+ν Γ(ν F (, ν; 3 ;

5 A modified Bessel function sum 5 + π ( k+ ( k C k, (3.0 π where the coefficients C k re defined in (3.9. Finlly, the lternting version of (. is given y ( n e n Ŝ J,ν (, = ( J ν (n. nν It is esily seen tht Ŝ J,ν (, = S J,ν (, S J,ν (,, from which the expnsion s 0 cn e otined from Theorem. 4. The expnsion involving the modified Bessel function The sme procedure cn e rought to er on the sum involving the modified Bessel function K ν (z given y e n S K,ν (, = ( K ν(n ( > 0, > 0, (4. nν where we consider only ν 0, since K ν (x = K ν (x. We remove from considertion the cse when ν hs hlf-integer vlues when the ove sum reduces to sum of exponentils; see the ppendix. From the stndrd definition K ν (x = π(i ν (x I ν (x/(sinπν vlid for non-integer ν nd the series expnsions of I ±ν (x, we otin S K,ν (, = π sinπν { ( k ν k!γ( ν+k Li k+ν(e when >. Use of the integrl representtion [, (0.3.3] ( x ν K ν (x = 4πi c+ i c i Γ(sΓ(s ν( x s ds ( } k k!γ(+ν+k Li k(e ( rg x < π, where c > mx{0,ν}, together with (3., shows tht S K,ν (, cn e expressed s the doule Mellin-Brnes integrl S K,ν (, = ( i Γ(sΓ(s νγ(t t ( πi s ζ(s+tdsdt, (4.3 i where, s in (3.3, the integrtion contours re indented to pss to the right of the pole of ζ(s+t nd those of the gmm functions. Evlution of the residues t s = n, t = m nd s = n+ν, t = m (m,n = 0,,,... yields, provided ν 0,,, 3,,..., T = = π sinπν { π sinπν ( m m ( n m!n! Γ(+ν π { ( ν ζ( n m+ν Γ(n+ ν ζ( n m } Γ(n++ν { ( m+ ( n ( m+n A m,n π π ( m ( n ν }} +B m,n π π There is n error in the sector of vlidity in [, (0.3.3] nd lso in [3, p. 4]. (4. (4.4

6 6 R.B. Pris for < π nd < π. The coefficients A m,n re defined in (3.5 nd B m,n = Γ( m+n+ νγ( m+n+ ν Γ( m+ Γ( m+n!γ(n+ ν ζ(m+n+ ν cscπ( m ν. (4.5 To del with the pole of the zet function when s+t =, we otin from (4.3 i ( s Γ(sΓ(s νγ( s ds, 4πi i where the integrtion pth is indented to seprte the poles of Γ(s nd Γ(s ν on the left from those of Γ( s on the right (ν,, 3,,.... Evlution of the residues on the right yields the contriution T = ( k ( k Γ( k! k + Γ( k + ν (4.6 ν π = Γ( ( ν Γ( ν < 3 F (, ν; ;. The lte terms in this expnsion possess the ehviour ( k ( k ν / (/ k, so tht the series in (4.6 converges for (ν 0. Then we hve the following expnsion. Theorem 3. Let ν 0,,, 3,,... nd. Then the following convergent expnsion holds S K,ν (, = Γ( ν sinπν ( ν πcscπν + 4Γ(+ν Γ( ν F (, ν; 3 ; { ( m+ ( n ( m ( n ν } ( m+n A m,n +B m,n π π π π provided < π nd < π. The coefficients A m,n nd B m,n re defined in (3.5 nd (4.5. The doule sum ppering in (4.7 cn e written in the lterntive (symptotic form sinπν where the coefficients D k re defined y with C k defined in (3.9 nd (4.7 ( k ( k D k, (4.8 π D k = C k, D k+ = C k+ +( k C k (k = 0,,,... C k = ( n ν ( n B k,n ( < π. π n=0 In the cse of integer vlues of ν we need to tke the limiting vlue of the doule sum T in (4.4. For exmple, if ν = 0 we hve upon setting ν = ǫ π T = lim ǫ 0 sinπǫ ( m m ( n m!n! { ζ( n m+ǫ( ǫ ζ( n m } Γ(n+ ǫ Γ(n++ǫ

7 = A modified Bessel function sum 7 ( m m ( } {ζ n m!(n! ( n m+ζ( n m{ψ(n+ log }. Upon noting tht the hypergeometric function in (4.7 reduces to n rcsin when ν = 0, we then find e n K 0 (n = ( m m ( } {ζ n m!(n! ( n m+ζ( n m{ψ(n+ log } + π rcsin(/ (/ ( < provided < π nd < π. Finlly, the lternting version of (4. is given y It is esily seen tht Ŝ K,ν (, = ( =. ( n e n ( K ν (n. nν Ŝ K,ν (, = S K,ν (, S K,ν (,, from which the expnsion s 0 cn e otined from Theorem 3. (4.9 Appendix: The cse of S K,ν (, when ν hs hlf-integer vlues When ν = m+ we hve x m Km+ (x = e x x m+ m c k (m x k (m = 0,,,..., where the c k (m re known constnts with c 0 (m =. Then S K,ν (, reduces to S K,m+ (, = m m+ m c k (m k e (+n n m+k+. The sum n n w e αn for integer w hs een discussed in [3, (4..0], where it is shown tht n w e αn = ( αw Γ(w {γ log α+ψ(w}+ ( k ζ(w kα k k! (α < π, where γ is the Euler-Mscheroni constnt, ψ(x is the logrithmic derivtive of the gmm function nd the prime on the summtion sign denotes the omission of the term corresponding to k = w. From this lst result it is then possile to construct the expnsion of S K,ν (, when ν = m+. The simplest cse when ν = yields when + < π. { π S K, (, = k= ( k } ζ( k(+ k log(+ k! (A.

8 8 R.B. Pris References [] F.W.J. Olver, D.W. Lozier, R.F. Boisvert nd C.W. Clrk (eds., NIST Hndook of Mthemticl Functions, Cmridge University Press, Cmridge, 00. [] R.B. Pris, The evlution of single Bessel function sums, Mth. Aetern 8 ( [3] R.B. Pris nd D. Kminski, Asymptotics nd Mellin-Brnes Integrls, Cmridge University Press, Cmridge, 00. [4] S.B. Tri cković, M.V. Vidnović nd M.S. Stnković, On the summtion of series in terms of Bessel functions, J. Anl. Appl. 5 (