arxiv: v1 [math.ca] 25 Mar 2015

Transcription

1 The asymptotic expansion of a generalisation of the Euler-Jacobi series arxiv: v [math.ca] 25 Mar 205 R. B. Paris Division of Computing and Mathematics, University of Abertay Dundee, Dundee DD HG, UK Abstract We consider the asymptotic expansion of the sum S p (a;w) = n= e anp n w as a 0 in arg a < 2π for arbitrary finite p > and w > 0. Our attention is concentratedmainlyonthecasewhen pandw areboth evenintegers, wherethe expansion consists of a finite algebraic expansion together with a sequence of increasingly subdominant exponential expansions. This exponentially small component produces a transformationfor S p (a;w) analogousto the well-knownpoisson-jacobi transformation for the sum with p = 2 and w = 0. Numerical results are given to illustrate the accuracy of the expansion obtained. Mathematics Subject Classification: 30E5, 33B0, 34E05, 4A30 Keywords: Euler-Jacobi series, Poisson-Jacobi transformation, asymptotic expansion, inverse factorial expansion. Introduction We consider the asymptotic expansion of the sum S p (a;w) = n= e anp n w (.) as the parameter a 0 in arg a < 2π, where p > 0 and, for convenience, w will be supposed throughout to be real and positive. When w = 0, this sum is known as the Euler-Jacobi series and when a = 0 then S p (0;w) reduces to the Riemann zeta function ζ(w) (provided R(w) > ). Consequently, the series in (.) can also be viewed as a smoothed Dirichlet series for ζ(w). The asymptotics of S p (a;w) as a 0+ for p a rational fraction and w < 0 was considered by Ramanujan and is discussed in [, Chapter 5]. It was shown that the expansion in this case consisted of an asymptotic sum involving the Riemann zeta function. A hypergeometric function approach for the Euler-Jacobi series when p is

2 2 R. B. Paris a rational fraction and w = 0 has been discussed at length in the monograph [4], where great emphasis was placed on obtaining exponentially small expansions. The case w = 0 and arbitrary p > 0, thus generalising the work in [4], has been investigated in [6, 8.] also for a 0+ using both a Mellin-Barnes integral approach and also a saddle point analysis of a Laplace-type integral representation. A similar Mellin-Barnes integral approach for this latter case has also been independently considered in [2]. The saddle point approach is instructive for understanding the appearance of exponentially small terms in the expansion as the parameter p increases. It was established in [6, 8.] that an additional exponentially small contribution appears when p = 2 + 4k, k = 0,,2,... Indeed, the appearance of the first exponentially small expansion as p passes through the classical value p = 2 was demonstrated to be associated with a Stokes phenomenon [6, 8..7]. In the case p = 2, w = 0 the behaviour of the sum S 2 (a;0) as a 0 can be obtained from the classical Poisson-Jacobi transformation given by S 2 (a;0) = e an2 = 2 n= π a 2 + π a e π2 n 2 /a n= (.2) valid for all values of a in arg a < 2π. This well-known transformation relates a sum of Gaussian exponentials involving the parameter a to a similar sum with parameter π 2 /a. In the small-a limit, the convergence of the sum on the left-hand side becomes slow, whereas the sum on the right-hand side converges rapidly. Various proofs of (.2) exist in the literature; see, for example, [6, p. 20], [8, p. 60] and [9, p. 24]. ThedominantasymptoticexpansionofS p (a;w) asa 0in arg a < 2 π forgeneral p > 0 and w > 0 is relatively straightforward and is found to consist, in general, of a single term proportional to a (w )/p, together with a series in ascending powers of a with coefficients involving the Riemann zeta function (the algebraic expansion). When 0 < p, the expansion is convergent and the result is exact; when p > the expansion is asymptotic as a 0. The most interesting case arises when p and w are both even integers. The above-mentioned algebraic expansion then terminates after a finite number of terms, and it becomes essential for accurate estimation to also include a subdominant sequence of exponentially small expansions. This exponentially small component produces a transformation for S p (a;w) analogous to the Poisson-Jacobi transformation in (.2), but valid as a 0 in arg a < 2π. This similarly involves a finite sequence of series similar to (.2) with a in the exponential replaced by an inverse power of a, but with each term decorated by an asymptotic series in ascending powers of a /(p ). The approach we employ in this paper is based on a Mellin-Barnes integral representation for S p (a;w) and is similar to that described in [6, 8..4]. An algorithm for the determination of the coefficients in the exponentially small asymptotic series is described. The case p = 2 when w is an even integer has been recently discussed in [7], where the coefficients in the decorating asymptotic series can be given in closed form. An application of the series when p = 2, with w = 2 and w = 4, has arisen in the geological problem of thermochronometry in spherical geometry [0]. 2. An expansion for S p (a;w) as a 0 when w, p 2,4,... We examine the expansion of the series S p (a;w) defined in (.) as a 0 in the sector arg a < 2π, where p > 0 and for simplicity in presentation we shall assume

3 Poisson-Jacobi-type transformation 3 throughout real values of w > 0. The case p =, w = 0 may be excluded from our consideration since the series in this case is summable as a geometric progression. Our starting point is the well-known Cahen-Mellin integral (see, for example, [6, 3.3.]) z α e z = c+ i Γ(α s)z s ds (z 0, arg z < 2πi 2π), (2.) c i where c < R(α) so that the integration path passes to the left of all the poles of Γ(α s) situated at s = k+α (k = 0,,2,...). Then, it follows that S p (a;w) = e anp n n= w = n= c+ i = 2πi c i n w 2πi c+ i c i Γ( s)(an p ) s ds Γ( s)ζ(w ps)a s ds ( arg a < 2π), (2.2) upon reversal of the order of summation and integration, which is justified when c > max{0,(w )/p}, and evaluation of the inner sum in terms of the Riemann zeta function. The integrand in (2.2) possesses simple poles at s = k (k = 0,,2,...) and s = s 0 (w )/p, except if s 0 = M (that is, w = pm + ), where M is a non-negative integer, when the pole at s = s 0 is double. The residue at the double pole is obtained by making use of the fact that ζ(s) /(s ) + γ in the neighbourhood of s =, where γ is Euler s constant, to find ( a) M {γ M! p log a+ } p ψ(m +) (M = 0,,2,...), where ψ(x) is the logarithmic derivative of the gamma function. The case when w and p are even positive integers requires a separate treatment which is discussed in Section The case 0 < p < We first consider the case 0 < p <. The integration path in (2.4) can be made to coincide with the imaginary s-axis together with a suitable indentation to lie to the left of the poles at s = 0 and s = s 0 (when 0 < w < ). Then use of the functional relation for ζ(s) given by [9, p. 269] shows that the integrand can be written as ζ(s) = 2 s π s ζ( s)γ( s)sin 2πs (2.3) (2π) w ζ( w+ps) Γ( w +ps) Γ(+s) sin 2 π(ps w) sinπs a s (2π) ps. With s = Re iθ, where R and is chosen so that the arc passes between the poles on the positive real axis, the logarithm of the dominant real part of the integrand is controlled by (p )Rcosθ log R+O(R). When θ < 2 π, this last expression tends to when 0 < p <. Consequently, the integration path can be bent back to enclose the poles of the integrand to yield the

4 4 R. B. Paris convergent result ( ) k S p (a;w) = J p (a;w)+ ζ(w kp)a k k! (0 < p < ), (2.4) where J p (a;w) = ( ) w p Γ a (w )/p (w pm +) p ( a) M {γ M! p log a+ ψ(m +)} (w = pm +) p and the prime on the sum over k denotes the omission of the term corresponding to k = M when w = pm +. When p =, the sum S (a;w) is given by (2.4) but now the sum over k on the right-hand side converges when a < 2π; see [6, 4.2.2] for details. If we let p =, w = 0 then use of the facts that ζ( 2k) = B 2k /(2k) and ζ( 2k) = 0, where B 2k are even-order Bernoulli numbers, shows that S (a;0) = a 2 + a k= B 2k (2k)! a2k ( a < 2π), which correctly reduces to the trivial summation /(e a ) by application of [5, Eq. (24.2.)]. 2.2 The case p > When p >, the integration path in (2.2) cannot be bent back over the poles and we proceed in a similar manner to that described for the case w = 0 in [6, 8..4]. Consider the integral taken roundthe rectangular contour with vertices at c±it, c ±it, where c > 0. The contribution from the upper and lower sides s = σ ± it, c σ c, vanishes as T provided arg a < 2π, since from the behaviour Γ(σ ±it) = O(t σ 2 e 2 πt ), ζ(σ ±it) = O(t µ(σ) log A t) (t ), where for σ and t real µ(σ) = 0 (σ > ), 2 2 σ (0 σ ), 2 σ (σ < 0), A = (0 σ ), A = 0 otherwise, the modulus of the integrand is controlled by O(T σ+µ(σ) 2 log Te T ), with = 2π arg a. Displacement of the integration path to the right over a finite set of poles then yields (provided w and p are not even integers) N ( ) k S p (a;w) = J p (a;w)+ ζ(w kp)a k +R N, (2.5) k! where N is a positive integer such that N > s and the prime on the sum over k again denotes the omission of the term corresponding to k = M when w = pm +.

5 Poisson-Jacobi-type transformation 5 The remainder R N is given by R N = c+ i 2πi = (2π)w 2πi c i c+ i c i Γ( s)ζ(w ps)a s ds ζ( w+ps) Γ( w +ps) Γ(+s) sin 2 π(ps w) sinπs a s (2π) psds where c = N 2 and the second expression follows from (2.3). Upon use of the result ζ(x+iy) < ζ(x) when x >, we obtain the bound where φ = arg a and ( ) a N R N < (2π) w ζ(pn p(s 0 2 )) 2 (2π) p e φt F(t)dt, (2.6) F(t) = Γ( w+ps) Γ(+s) cosh 2 πpt cosh πt (s = N 2 +it). On the integration path, F(t) is regular and satisfies F(t) = O(e 2 π t ) as t ±. Hence the integral in (2.6) is convergent and independent of a provided φ < 2π. It then follows that R N = O(a N 2 ) (a 0 in arg a < 2 π). The expansion (2.5) is the dominant algebraic expansion associated with S p (a;w) valid as a 0 in arg a < 2π, provided w (> 0) and p are not even integers when the sum in (2.5) is finite. The same analysis can be applied to the case with non-positive w to yield the Berndt-Ramanujan result [, Theorem 3., p. 306] S p (a; w) = ( ) +w N p Γ a (+w)/p ( ) k + ζ( w kp)a k +O(a N 2) (w 0) p k! as a 0 in arg a < 2π. The reflection formula (2.3) can be employed to convert the argument of the zeta function to a positive form. 3. The expansion of S p (a;w) when w and p are even integers Throughout this section we let w and p be even positive integers, with w = 2m where m =,2,... In this case s 0 = (2m )/p, which cannot equal an integer and so no double pole can arise. More importantly, there is now only a finite set of poles of the integrand in (2.2) at s = s 0 and s = 0,,2,...,K, where K = w/p, since the poles of Γ( s) at s = K+k (k =,2,...) are cancelled by the trivial zeros of the zeta function ζ(s) at s = 2, 4,... This has the consequence that the integrand is holomorphic in R(s) > max{s 0,K}, so that further displacement of the contour can produce no additional algebraic terms in the expansion of S p (a;w). Thus, we find from (2.4) upon displacement of the integration path to the right over the poles of the integrand S p (a;w) = ( ) w p Γ a (w )/p + p K ( ) k ζ(w pk)a k +( ) m (2π) w I L, (3.) k!

6 6 R. B. Paris where I L = ( )m (2π) w 2πi L Γ( s)ζ(w ps)a s ds (3.2) and L denotes a path parallel to the imaginary axis with R(s) > (w/p) + δ, with δ denoting an arbitrary positive quantity. This is easily seen to satisfy the requirement R(s) > max{s 0,K 0 } necessary for the validity of (3.). We now employ the functional relation for ζ(s) in (2.3) to convert the argument of the zeta function in (3.2) into one with real part greater than unity. The integral in (3.2) can then be written in the form I L = Γ( w +ps) sin 2 ζ( w +ps) πps 2πi L Γ(+s) sinπs χ s ds, where χ = (2π) p a. In [6, 8..4], the zeta function appearing in the above integrand was written as an infinite series. Here we follow a suggestion made in [2] and retain this function in the integrand; see also [4, 8]. Making use of the expansion (see, for example, [6, p. 368]) sin 2 πps sinπs { N = 2 cosπ( 0 (p/2 even) 2 p 2r )s+ (p/2 odd) r=0, N = [ 4 p], where square brackets denote the nearest integer part, we obtain I L = N r=0 {J + r +J r }+ { 0 (p/2 even) J (p/2 odd). (3.3) Here we have defined the integrals J ± r and J by and J r ± = Γ( w +ps) ζ( w +ps) (χe πi(p/2 2r ) ) s ds (3.4) 2πi L Γ(+s) J = ζ( w+ps) Γ( w+ps) χ s ds. (3.5) 2πi L Γ(+s) 3. Asymptotic evaluation of J ± r and J The integrals J ± r and J have no poles in the half-plane R(s) > (w/p) +δ, so that we can displace the path L as far to the right as we please. On such a displaced path s is everywhere large. Let M denote an arbitrary positive integer. The ratio of gamma functions appearing in (3.4) and (3.5) may then be expanded by making use of the result (for p > ) given in [6, p. 53] Γ( w +ps) Γ(+s) = A { M 2π (hκκ ) s } ( ) j c j Γ(κs+ϑ j)+ρ M (s)γ(κs+ϑ M), (3.6) where c 0 =, ρ M (s) = O() as s in arg s < π and κ = p, h = p p, ϑ = 2 w, A = (2π) 2 κ 2 ϑ p ϑ. (3.7) The nearest integer part corresponds to [x] = N when x is in the interval (N 2,N + 2 ]. Note that when p = 2, we have N = 0 and the above expansion contains no information.

7 Poisson-Jacobi-type transformation 7 The coefficients c j c j (w,p) (0 j 4) are listed in [6, pp ] where an algorithm for their determination is described; see Section 4 for details. Substitution of the expansion (3.6) into the integrals J ± r in (3.4) then produces J ± r = A 2π M ( ) j c j 2πi = A M (Xe πiψr ) ϑ j( )j c j 2πκ 2πi L ζ( w+ps)γ(κs+ϑ j)(xe πiψr ) κs ds+r ± M,r L Γ(u)ζ(qu+λ j )(Xe πiψr ) u du+r ± M,r. Here we have made the change of variable u κs+ϑ j, with L denoting the modified integration path, and have defined X := κ(hχ) /κ, ψ r := 2 p 2r κ, q := p κ, λ j := + κ (w +p(j 2 )) (3.8) together with the remainders R ± M,r = A ρ M (s)ζ( w+ps)γ(κs+ϑ M)(Xe πiψr ) κs ds. (3.9) 4πi L We note that λ j > 0 for j 0 when w > 0 and p 2. The above integrals appearing in J r ± may now be evaluated by means of (2.2), when we replace s by s and allow the integration path (c i,c + i) to coincide with the path L, to yield the sum S q (Xe πiψr ;λ j ) as defined in (.). This evaluation is valid provided that the variable Xe πiψr satisfies the convergence condition arg(xe πiψr ) < 2π; that is arg a κ πψ r < 2 π (0 r N ). It is routine to verify that these conditions are met when arg a < 2 π. Thus we find J r ± = A M ( ) j c j (Xe πiψr ) ϑ j S q (Xe πiψr ;λ j )+R ± M,r 2πκ. (3.0) Bounds for the remainders of the type R ± M,r have been considered in [6, p. 7, Lemma 2.7]; see also [3, 0.]. The integration path in (3.9) is such that R( w +ps) >, so that we may employ the bound ζ(x+iy) ζ(x) for x >. A slight modification of Lemma 2.7 in [6, p. 7] then shows that R ± M,r = O (X ϑ M e Xe πiψr) (3.) as a 0 in the sector arg a < 2 π. An analogous procedure applied to J in (3.5) shows that as a 0 in arg a < 2 π. J = A M ( ) j c j X ϑ j S q (X;λ j )+O(X ϑ M e X ) (3.2) 2πκ 3.2 The expansion of S p (a;w) The expansion of I L as a 0 in arg a < 2 π then follows from (3.), (3.3), (3.0) and (3.2). We obtain the following theorem.

8 8 R. B. Paris Theorem. Let m and M be positive integers. Then, when w = 2m and p is also an even positive integer, with K = w/p and N = [ 4p] (with square brackets denoting the nearest integer part), we have the expansion valid as a 0 in arg a < 2 π with S p (a;w) = ( ) w p Γ a (w )/p + p I L = N r=0 K ( ) k ζ(w pk)a k +( ) m (2π) w I L, (3.3) k! E r (a;w,p)+δ pp Ê N (a;w,p), (3.4) where δ pp is the Kronecker symbol with p = 4N +2. The sums E r (a;w,p) are given by 2 E r (a;w,p) = A 2πκ M ( ) j c j (Xe πiψr ) ϑ j S q (Xe πiψr ;λ j )+R M,r (3.5) ± for 0 r N, where X = κ(h(2π) p /a) /κ, ψ r = ( 2p 2r )/κ (0 r N ), λ j = +(w +p(j 2 ))/κ, q = p/κ and the parameters κ, h, ϑ and A are defined in (3.7). The leading coefficient c 0 = and c j c j (w,p) (j ) are discussed in Section 4. The sum ÊN(a;w,p) is also given by (3.3) when we put ψ N 0 and omit the summation ±. The remainders R M,r satisfy the bound R M,r = O(max{X ϑ M e Xe±πiψr }) (0 r N). It is seen from (3.3), (3.4) and (3.5) that the sum S p (a;w) has been expressed in termsofthesumss q (Xe πiψr ;λ j ), whichinvolve thereciprocal poweroftheasymptotic variable a scaling like a /κ. Thus, as a 0 the argument X. It is obvious from the definition in (.) (when q > 0) that S q (z;λ j ) e z (z in arg z < 2 π), so that the E r (a;w,p) represent a series of exponentially small expansions of increasing subdominance in the small-a limit. In addition, the number of exponentially small expansions increases by one each time p increases by 4. By means of a saddle-point analysis in the case w = 0, this was demonstrated to correspond to a Stokes phenomenon when p was allowed to vary continuously through the values p = 2,6,0,...; see [6, 8..2, 8..7]. Finally, we remark that the exponents p and q are conjugate exponents [2], since p + q =. When a is a real parameter, the expansion in Theorem can be expressed in a different form by using (.) to represent the S q (Xe πiψr ;λ j ) as infinite sums. Then from (3.5) we obtain the following theorem: Theorem 2. Let w and p be even positive integers, N = [ 4p] and M be a positive integer. Then, the exponentially small expansions in (3.5) valid as a 0+ can be written in the form E r (a;w,p) = A n w Xn ϑ e Xncosπψr Υ n,r (0 r N ), (3.6) πκ n= 2 The symbol signifies that the series with ± signs are to be added. ±

9 Poisson-Jacobi-type transformation 9 Ê N (a;w,p) = A 2πκ The Υ n,r (0 r N) have the asymptotic expansions Υ n,r = M n w Xn ϑ e Xn Υ n,n. (3.7) n= ( ) j c j X j n cos[x nsinπψ r +π(j ϑ)ψ r ]+O(X M n ), (3.8) where X n = Xn p/κ = κ(h(2πn) p /a) /κ, ψ r = ( 2 p 2r )/κ, ψ N 0 and the other quantities are as defined in Theorem. The result in (3.3) and (3.4) is the analogue of the Poisson-Jacobi transformation in (.2) corresponding to w = 0, p = 2. In this latter case, N = 0 and, from (3.4), I L = E 0 (a;0,2). The ratio of gamma functions in (3.6) is replaced by the single gamma function Γ(s+ 2 ) by the duplication formula for the gamma function, with the result that c 0 =, c j = 0 (j ) and consequently Υ n,n = for all n. Then, since K = 0 and ζ(0) = 2, (3.3), (3.4) and (3.7) reduce to (.2). The resulting expansion is valid for all values of the parameter a (not just a 0) satisfying arg a < 2π. When w = 0, p = 2m, the expansions (3.6) and (3.7) reduce to those given in [7]; see Section 5. We observe that the n-dependence in the sums E r (a;w,p) from the factor n w Xn ϑ is given by n w n pϑ/κ = n (2w+p 2)/(2κ). Since p 2 and w > 0, this is seen to correspond to a negative power of n. 4. The coefficients c j We describe an algorithm for the computation of the coefficients c j c j (w,p) that appear in the exponentially small expansions E r (a;w,p) in (3.5) and (3.8). The expressionfor theratio oftwo gamma functionsin(3.6), withα w forconvenience, takes the form Γ(α+ps) Γ(+s)Γ(κs+ϑ) = A { M 2π (hκκ ) s c j ρ M (s) + ( κs ϑ) j ( κs ϑ) M where the parameters κ, h, ϑ and A are defined in (3.7) and (α) j = Γ(α + j)/γ(α) is the Pochhammer symbol. If we introduce the scaled gamma function Γ (z) = Γ(z)/( 2πz z 2e z ), then we have Γ(βs+γ) = Γ (βs+γ)(2π) 2 e βs (βs) βs+γ 2 e(βs;γ), }, where e(βs;γ) := exp [(βs+γ 2 )log(+ γ ] βs ) γ. The above ratio of gamma functions may therefore be rewritten as R(s)G(s) = M c j ( κs ϑ) j + ρ M (s) ( κs ϑ) M (4.)

10 0 R. B. Paris as s in arg s < π, where R(s) = e(ps; α) e(s;)e(κs;ϑ), G(s) = Γ (α+ps) Γ (+s)γ (κs+ϑ). We now let ξ := (κs) and expand R(s) and G(s) for ξ 0 making use of the well-known expansion [6, p. 7] Γ (z) ( ) k γ k z k ( z ; arg z < π), where γ k are the Stirling coefficients, with γ 0 =, γ = 2, γ 2 = 288, γ 3 = , γ 4 = ,.... After some straightforward algebra we find that R(s) = + ξ 2 { α(α )κ p } ϑ(ϑ ) +O(ξ 2 ), G(s) = + ξ ( p ) +O(ξ 2 ), 2 p so that upon equating coefficients of ξ in (4.) we can obtain c. The higher coefficients can be obtained by matching coefficients recursively with the aid of Mathematica to find [6, p. 47] c 2 = c 0 =, c = 24p (2 5p+2p2 2w +2pw +2w 2 ), 52p 2(4+28p 87p2 +28p 3 +4p 4 +48w 26pw +24p 2 w +44p 3 w 96w 2 20pw p 2 w 2 96w pw 3 +44w 4 ),.... (4.2) The rapidly increasing complexity of the coefficients with j 3 prevents their presentation. However, this procedure is found to work well in specific cases when the various parameters have numerical values, where up to a maximum of 00 coefficients have been so calculated. In Table we present the values 3 of the coefficients c j for j 8 in the specific examples considered in Section 5. When p = 2, useof the duplication formulashows that the ratio of gamma functions in (3.6) becomes Γ( 2 2 w+s)γ( 2 w+s) Γ(+s) = M ( ) j c j Γ(s+ 2 w j)+ρ M(s)Γ(s+ 2 w M). The coefficients c j in this case can be expressed in closed form as [6, p. 53] c j = ( 2 w) j( w) j j! = 2 2j (w) 2j j! (p = 2). (4.3) Finally, we mention that when w = 0 (corresponding to the Euler-Jacobi series) the coefficients c j are listed for p 2 and j 8 in [6, p. 374]. 3 In the tables we write the values as x(y) instead of x 0 y.

11 Poisson-Jacobi-type transformation Table : The coefficients c j ( j 8) for different p and w. j p = 4, w = 2 p = 4, w = 4 p = 6, w = 2 p = 6, w = (0) (0) (0) (0) (0) () (0) () () ().38009() () () (2) () (2) (2) (3) (2) (3) (3) (4) (3) (4) (4) (5) (4) (5) (5) (6).39672(5) (6) 5. Numerical results and concluding remarks We present some examples of the expansion of S p (a;w) given in Theorem when p and w = 2m are even integers. For convenience in presentation, we extract the factor e z from the sum S q (z;λ j ) by writing S q (z;λ j ) = e z e z(nq ) Ŝ 2 (z;λ j ), Ŝ q (z;λ j ) := n n= λ. j It follows that, when q > 0, Ŝq(z;λ j ) = O() as z in arg z < 2 π. Example. In the case p = 2, we have κ =, q = 2, N = 0, K = m, p = 2, ψ 0 0 and X = π 2 /a. The quantity δ pp = so that the exponentially small component of S 2 (a;2m) consists of the single term Ê0(a;2m,2). From (3.3), (3.4), (3.5) and (4.3) we therefore find = ( ) m ( a π S 2 (a;2m) ( ) 2m 2 Γ m ( ) k a m 2 ζ(2m 2k)a k 2 k! ) 2m { M 2 e π 2 /a ( ) j (2m) 2j j! ( ) a jŝ2 4π 2 (π 2 /a;2m+2j)+o(a )} M (5.) as a 0 in the sector arg a < 2π. The expansion in this case has been given in an equivalent form in [7]. Example 2. When p = 4, we have κ = 3, q = 4 3, N =, K = 2 m, p = 6 and X = 3( 2 π)4/3 a /3, λ j = 3 (2m+4j +). The quantity δ pp = 0 so that there is the single exponentially small term E 0 (a;2m,4) with ψ 0 = 3. Then we find the expansion as a 0 in arg a < 2π given by S 4 (a;2m) ( ) 2m 4 Γ a (2m )/4 4 K ( ) k ζ(2m 4k)a k k!

12 2 R. B. Paris ( ) 2a (4m )/6 { M = ( ) m e Xe 3 πi 3 πiϑ π ± ( ) j c j (Xe 3 πi ) j Ŝ4 3 } (Xe 3 πi ;λ j )+O(a M/3 ) where the coefficients c j c j (2m,4) can be obtained from (4.2) and [6, p. 47] as c 0 =, c = 48 (7+36m+24m2 ), c 2 = 4608 ( m+704m m m 4 ), c 3 = ( m m m m m m 6 ),.... These coefficients are listed in Table for j 8 when m = and m = 2. Example 3. When p = 6, we have κ = 5, q = 6 5, N =, K = 3m and X = 5( 3 π)6/5 a /5, λ j = 5 (2m+6j +2). (5.2) In this case p = 6, so that δ pp = and there are now two exponentially small expansions E 0 (a;2m,6), with ψ 0 = 2 5, and Ê(a;2m,6). Then, as a 0 in arg a < π, we have the expansion 2 S 6 (a;2m) ( ) 2m 6 Γ a (2m )/6 6 ( ) 3a (4m )/0 { M = ( ) m e Xe 2 5 πi 2 5 πiϑ π ± ( ) 3a (4m )/0 { M +( ) m e X ( ) j c j π X j The first few coefficients c j c j (2m,6) are K ( ) k ζ(2m 6k)a k k! ( ) j c j (Xe 2 5 πi ) j Ŝ6 5 c 0 =, c = 36 (+30m+2m2 ), } (Xe 2 5 πi ;λ j )+O(a M/5 ) } Ŝ6(X;λ j )+O(a M/5 ). (5.3) 5 c 2 = 2592 ( m+46m2 +392m 3 +44m 4 ), c 3 = ( m m m m m m 6 ). These coefficients are listed in Table for j 8 when m = and m = 2. We show the results of numerical calculations to demonstrate the achievable accuracy of the expansion in Theorem. We define the difference between S p (a;w) and the finite algebraic expansion by S p,w S p,w (a) := S p (a;w) ( ) w p Γ a (w )/p p K ( ) k ζ(w pk)a k. (5.4) k!

13 Poisson-Jacobi-type transformation 3 Table 2: Values of the absolute error in the computation of S p(a;w) defined in (5.4) using the expansions (5.) and (5.2). The value of the index j 0 corresponds to optimal truncation of the subdominant expansion E 0(a;w,p). p = 2, w = 2 p = 2, w = 4 a S p,w S p,w E 0 j 0 S p,w S p,w E 0 j ( 06) 6.637( 09) ( 07) 3.642( 08) ( 07) 8.089( 2) 9.296( 09) 4.659( ) ( 0).260( 7) ( 2) 7.635( 7) ( 24).542( 43) ( 26) 9.830( 43) ( 46).486( 86) ( 49) 9.63( 86) 95 p = 4, w = 2 p = 4, w = 4 a S p,w S p,w E 0 j 0 S p,w S p,w E 0 j ( 03).329( 06) ( 04) 8.742( 06) ( 04) 2.749( 08) 4.805( 05) 4.879( 07) ( 05) 2.56( 0) ( 06).420( 09) ( 09) 4.62( 7) ( 09) 3.04( 6) ( 4).033( 36) ( 6) 6.799( 36) 5 In Table 2 we present the absolute error in the computation of S p (a;w) for different values of the parameter a in the two cases p = 2 and p = 4, with w = 2 and w = 4 using the expansions given in (5.) and (5.2). The first column in each entry displays the absolute valueof S p,w ; that is, theaccuracy achievable with justthealgebraic expansion and no subdominant exponential terms. The second column shows the absolute error when the single optimally truncated exponential expansion E 0 (a;w,p) (denoted by E 0 in the table) is included. The optimal truncation index j 0, corresponding to truncation of the exponential expansion E 0 (a;w,p) at, or near, the least term in magnitude, is indicated in the final column. The situation when there is only a single subdominant exponentially small expansion present is straightforward: this sum is truncated at some suitable point thereby introducing a truncation error. If truncation is optimal, then the resulting error is exponentially more recessive than the parent exponential expansion. However, in the case of two, or more, exponential expansions of different degrees of subdominance (corresponding to p 6) the situation is not so obvious. It is not clear, without further investigation, how the error from the truncated leading exponential series compares with the contribution from the next series. We illustrate this by considering the case p = 6 and w = 2 given in (5.3). In Table 3 we present the absolute error in the computation of S 6 (a;2) as a function of the parameter a. We show, in order, the value of S 6,2 and the absolute error in S 6,2 E 0 (a;2,6) when the leading subdominant exponential expansion E 0 (a;2,6) is optimally truncated at index j 0. The fourth column gives the absolute error when the first few terms of the second exponential expansion Ê(a;2,6) are included (for brevity in the table these exponential expansions are labelled E 0 and E, and their

14 4 R. B. Paris Table 3: Values of the absolute error in the computation of S 6(a;2) defined by (5.4) using the expansion (5.3). The value of the index j 0 corresponds to optimal truncation of the expansion E 0(a;2,6). a S 6,2 S 6,2 E 0 j 0 S 6,2 E 0, Min E 0 E (j = 0) ( 02) 3.780( 05) ( 05) 5.095( 05) ( 03) 3.037( 05) 8.200( 05).729( 05).9( 05) ( 04).93( 07) ( 08).228( 07).904( 07) ( 03).099( 08) ( 09) 9.090( 09) 2.48( 08) ( 04) 3.452( ) ( 2) 3.757( 2) 4.053( ) ( 07) 8.80( 7) ( 9) 3.024( 9) 9.20( 7) ( 0).522( 25) 5.963( 30).964( 30).564( 25) sum is denoted by E 0, ). The final two columns show the values of the least term (including prefactors) in E 0 (a;2,6) at optimal truncation and the values of the leading term (j = 0) of the sub-subdominant expansion Ê(a;2,6). A cursoryinspection of Table 3 shows that for a 0. theleading term of Ê(a;2,6) islessthantheminimumtermofe 0 (a;2,6) andconsequentlythatinclusionofê(a;2,6) cannot improve the accuracy. For a< 0.0, the reverse is true: the leading terms of Ê (a;2,6) aregreaterthantheminimumtermofe 0 (a;2,6) andtheirinclusiontherefore increases the overall accuracy. However, it is clear that in both cases the final accuracy achievable is limited by the optimal truncation of the leading subdominant expansion E 0 (a;2,6). Further improvement in the accuracy would require a hyperasymptotic treatment in order to deal with the divergent tails of E 0 (a;2,6) and Ê(a;2,6). A possible hyperasymptotic scheme for the Euler-Jacobi series with p = 3 and w = 0 has been discussed in [4, 8]. Finally we remark that the asymptotics of the alternating version of (.) can be deduced from the result in Theorem by making use of the identity ( ) ne anp n n= w = 2 w S p (2 p a;w) S p (a;w). References [] B. C. Berndt, Ramanujan s Notebooks, Part II, Springer-Verlag, New York, 985. [2] J. Boersma, Private communication (2000). [3] B. L. J. Braaksma, Asymptotic expansions and analytic continuations for a class of Barnes integrals, Compas. Math. 5 (963) [4] V. Kowalenko, N. E. Frankel, M. L. Glasser and T. Taucher, Generalised Euler-Jacobi Inversion Formula and Asymptotics Beyond All Orders, London Math. Soc. Lecture Notes Series 24, Cambridge University Press, Cambridge, 995. [5] F. W. J. Olver, D. W. Lozier, R. F. Boisvert and C. W. Clark (eds.), NIST Handbook of Mathematical Functions, Cambridge University Press, Cambridge, 200.

15 Poisson-Jacobi-type transformation 5 [6] R. B. Paris and D. Kaminski, Asymptotics and Mellin-Barnes Integrals, Cambridge University Press, Cambridge, 200. [7] R. B. Paris, A Poisson-Jacobi-type transformation for the sum n= n 2m exp( an 2 ) for positive integer m. 205 arxiv:50: [8] E. C. Titchmarsh, Introduction to the Theory of Fourier Integrals, Oxford University Press, Oxford, 975. [9] E. T. Whittaker and G. N. Watson, Modern Analysis, Cambridge University Press, Cambridge, 952. [0] R. A. Wolf, K. A. Farley and D. M. Kass, Modeling of the temperature sensitivity of the apatite (U-Th)/He thermochronometer, Chem. Geology 48 (998) 05 4.