arxiv: v1 [gr-qc] 13 Aug 2007

Transcription

1 Anomalously small wave tails in higher dimensions Piotr Bizoń, 1 Tadeusz Chmaj,,3 Andrzej Rostworowski 1 1 M. Smoluchowski Institute of Physics, Jagiellonian University, Kraków, Pol H. Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Kraków, Pol 3 Cracow University of Technology, Kraków, Pol Dated: November 4, 018 We consider the late-time tails of spherical waves propagating on even-dimensional Minkowski spacetime under the influence of a long range radial potential. We show that in six higher even dimensions there exist exceptional potentials for which the tail has an anomalously small amplitude fast decay. Along the way we clarify amend some confounding arguments statements in the literature of the subject. arxiv: v1 [gr-qc] 13 Aug 007 I. INTRODUCTION It is well known that sharp propagation of free waves along light cones in even-dimensional flat spacetimes, known as Huygens property, is blurred by the presence of a potential. Physically, the spreading of waves inside the light cone is caused by the backscattering off the potential. If the potential falls off exponentially or faster at spatial infinity, then the backscattered waves decay exponentially in time, while the long range potentials with an algebraic fall-off give rise to tails which decay polynomially in 1/t. The precise description of these tails is an important issue in scattering theory. There are two main approaches to this problem in the literature. On the one h, there are mathematical results in the form of various decay estimates. These results are rigorous, however they rarely give optimal decay rates inside the light cone provide very poor information about the amplitudes of tails. The notable exception is the work of Strauss Tsutaya [1] recently strengthened by Szpak [] where the optimal pointwise decay estimate for the tail was proved in four dimensions. Unfortunately, to the best of our knowledge, there is no analogous result in higher dimensions. On the other h, there are non-rigorous results in the physics literature based on perturbation theory. The most complete work in this category was done by Ching et al. [3] who derived first-order approximations of the tails for radial potentials. Although these results were originally formulated for partial waves in four dimensions, they can be easily translated to spherical waves in higher dimensions. Ching et al. noticed that there are exceptional potentials for which the first-order tail vanishes, however they did not pursue their analysis to the second order, apart from giving some dimensional arguments. The main purpose of this paper is to analyze the tails for such exceptional potentials in more detail. One of the physical motivations behind our work stems from the fact that this kind of potentials arise in the study of linearized perturbations of higher even-dimensional Schwarzschild black holes. The behavior of tails on the Schwarzschild background in well known in four dimensions see [4], [5], [3], [6], [7], [8], but not in higher even dimensions despite statements to the contrary in the literature [9]. Although our analysis is restricted to the flat background, it sheds some light on the problem of tails on the black hole background because the properties of tails are to some extent independent of what happens in the central region. The rest of the paper is organized as follows. In section we construct the iterative scheme for the perturbation expansion of a spherically symmetric solution of the linear wave equation with a potential. This scheme is applied in section 3 to derive the first second-order approximations of the tails for radial potentials which fall off as pure inverse-power at infinity. In section 4 we discuss the modifications caused by subleading terms in the potential. Section 5 contains numerical evidence confirming the analytic formulae from sections 3 4. Finally, in section 6 we give a heuristic argument to predict the behavior of tails outside Schwarzschild black holes in higher even dimensions. Technical details of most calculations are given in the appendix. Throughout the paper we use the succinct notation summation technics from the excellent book by Graham et al. [14]. In particular, we shall frequently use the following abbreviations x 0 := 1, x k := x x 1 x k 1, k > 0, 1 x 0 := 1, x k := x x+1 x+k 1, k > 0. II. ITERATIVE SCHEME We consider the wave equation with a potential in even-dimensional Minkowski spacetime R d+1 tφ φ+λvφ = 0. 3

2 The prefactor λ is introduced for convenience - throughout the paper we assume that λ is small which allows us to use it as the perturbation parameter. The precise assumptions about the fall-off of the potential will be formulated below. We restrict attention to spherical symmetry, i.e., we assume that φ = φt,r V = Vr. Then, equation 3 becomes Lφ+λVrφ = 0, L := t r d 1 r. 4 r We are interested in the late-time behavior of φt, r for smooth compactly supported or exponentially localized initial data. φ0,r = fr, t φ0,r = gr. 5 To determine the asymptotic behavior of solutions we define the perturbative expansion Born series φ = n=0λ n φ n, 6 where φ 0 satisfies initial data 5 all φ n with n > 0 have zero initial data. Substituting this expansion into equation 4 we get the iterative scheme Lφ n = Vφ n 1, φ 1 = 0, 7 which can be solved recursively. The zeroth-order solution is given by the general regular solution of the free radial wave equation which is a superposition of outgoing ingoing waves [10] where φ ret 0 t,r = 1 r l+1 l l k! k!l k! k=0 φ 0 t,r = φ ret 0 t,r+φadv 0 t,r, 8 a k u, φadv v u l k 0 t,r = 1 r l+1 l k+1 l k! a k v 1, 9 k!l k! v u l k u = t r, v = t+r are the retarded advanced times, respectively. Here in the following, instead of d, we use the index l defined by d = l +3 remember that we consider only odd space dimensions d. Note that for compactly supported initial data the generating function ax can be chosen to have compact support as well this condition determines ax uniquely. To solve equation 7 for the higher-order perturbations we use the Duhamel representation for the solution of the inhomogeneous equation Lφ = Nt,r with zero initial data φt,r = 1 t r l+1 0 dτ t+r τ t r τ k=0 ρ l+1 P l µnτ,ρdρ, 10 where P l µ are Legendre polynomials of degree l µ = r + ρ t τ /rρ note that 1 µ 1 within the integration range. This formula can be readily obtained by integrating out the angular variables in the stard formula φ = G ret N where G ret t,x = π l+1 1 Θtδ l t x is the retarded Green s function of the wave operator in d+1 dimensions see, for example, [11]. It is convenient to express 10 in terms of null coordinates η = τ ρ ξ = τ +ρ φt,r = 1 l+3 r l+1 t+r t r t r dξ ξ η l+1 P l µnη,ξdη, 11 ξ where now µ = r +ξ t/rξ η. Using this representation we can rewrite the iterative scheme 7 in the integral form 1 φ n t,r = l+3 r l+1 t+r t r t r dξ ξ η l+1 P l µvρη,ξφ n 1 η,ξdη. 1 ξ This master equation will be applied below to evaluate the first two iterates for a special class of potentials. It is natural to expect that for sufficiently small λ these iterates provide good approximations of the true solution.

3 3 III. PURE INVERSE-POWER POTENTIALS AT INFINITY In this section we consider the simple case below referred to as type I when the potential is exactly Vr = r α for r greater than some r 0 > 0. We assume that α >. The modifications caused by subleading corrections to the pure inverse-power decay of the potential will be discussed in section 4. A. Generic case We wish to evaluate the first iterate φ 1 t,r near timelike infinity, i.e, for r = const t. Thanks to the fact that φ 0 η,ξ has compact support we may interchange the order of integration in 1 drop the advanced part of φ 0 η,ξ to obtain φ 1 t,r = α l+3 r l+1 dη t+r t r ξ η l+1 α P l µφ ret 0 η,ξdξ, 13 where we have substituted V = α ξ η α. Plugging 9 into 13, after a long calculation see appendix A for the technical details, we get lα α 3 l φ 1 t,r = α+3l 1 α dη aη [ r ] α 1+l 1 n n l+n! α n α 1 l 1 n!l+n+1! 0 n α / n m n+l 1 m n l 1 l α +1 m m r m n α m. 14 Asymptotic expansion of 14 near timelike infinity yields the following first-order approximation of the tail φt,r λφ 1 t,r = λ Cl,α [ t α+l A+α+l B ] 1 +O t t, 15 where Cl,α = α+l 1 l+1!! A = aηdη, B = lα α 3 l, 16 aηηdη. 17 In general A 0 the tail decays as t α l, however there are nongeneric initial data for which A = 0 then the tail decays as t α l 1 ; in particular this happens for time symmetric initial data for which ax is an odd function. Remark 1. It is easy to check that if the function φt,r satisfies equation 4, then the function ψ = r l+1 φ satisfies the radial wave equation for the lth multipole t r +ll+1/r ψ +λvrψ = The late-time tails for this equation were studied by Ching et al. [3] who derived the formula equivalent to 15 via the Fourier transform methods.

4 4 B. Exceptional case It follows from 15 that if α is an odd integer satisfying 3 α l +1, then φ 1 t,r vanishes identically due to factor α 3 l in 16 there is no polynomial tail whatsoever in the first order. Thus, in order to compute the tail in this exceptional case we need to go the second order of the perturbation expansion. Using 1 proceeding as above we get the second iterate φ t,r = α l+3 r l+1 dη where φ ret 1 is the outgoing solution of the inhomogeneous equation t+r t r ξ η l+1 α P l µφ ret 1 η,ξdξ, 19 Lφ 1 = Vφ 0. 0 In general φ 1 is a sum of the solution of the homogeneous equation the particular solution of the inhomogeneous equation. The homogeneous part has the form 9 with a different generating function than a, but still compactly supported, thus for the same reason as above it gives no contribution to the tail. The particular solution of the inhomogeneous equation 0 reads l α/+1/ φ NH 1 l = α 1r α+l l α/+1/ q q α/ q α q q=0 φ H l 1 q r q, 1 where φ H l 1 q denotes the solution of the homogeneous equation with d = l 1 q+3 the same generating function a as in φ 0 see 9. The formula 1 can be easily derived by the method of undetermined coefficients we emphasize that this formula is valid only for odd α satisfying 3 α l + 1. Substituting 1 into 19, after a long calculation see appendix A for the technical details, we obtain the following asymptotic behavior near timelike infinity φt,r λ φ t,r = λ Dl,α [A+α+l 1 Bt ] 1t t +O α+l 1, where the coefficients A B are defined in 17 Dl,α = α+l l+1!! α 3 α 5 l 1 α +l α 1 l 1 l+α/ 1/, α/, α, 1 F α, α, α l 1/ Here F sts for the generalized hypergeometric function F 1. 3 a1,..., a m z = a k 1,..., ak m z k b 1,..., b n k 0 b k 1,..., bk k!. 4 n TABLE I: The first few coefficients Dl,α α l /5 40/3 3 96/ /5 4 64/7 1790/ / /7 We remark that the behavior Oλ t l+α 1 of the tail was conjectured before by Ching et al. [3] on the basis of dimensional analysis.

5 5 IV. GENERAL POLYNOMIALLY DECAYING POTENTIALS In this section we analyze how the presence of subleading corrections to the pure inverse-power asymptotic behavior of the potential affects the results obtained in section 3. We restrict ourselves to the most interesting common case below referred to as type II when near infinity Vr = 1 r α 1+ βr γ +o If Cα,l 0, then the dominant behavior of the tail is of course the same as in 15: 1 r α+γ, γ > 0. 5 φt,r λacl,αt α+l. 6 However, in the exceptional case, when Cα,l = 0, the situation is more delicate. As we showed above, in this case there is the second-order contribution to the tail given by φ t,r ADl,αt α+l 1. 7 In contrast to the type I case where the first-order tail vanishes identically, in the type II case the subleading term in the potential produces the first-order contribution which is given by 15 with α replaced by α+γ: φ 1 t,r βacl,α+γt α+γ+l, 8 assuming that α+γ is not an odd integer d otherwise one has to repeat the analysis for the next subleading term in the potential. Now, comparing the decay rates in 7 8 we conclude that the leading asymptotics of the tail is given by the first-order term λφ 1 t,r if γ α we call it subtype IIa, but otherwise, i.e. for γ > α subtype IIb, the second-order term λ φ t,r is dominant for t. Remark. In the context of equation 18 a formula analogous to 8 was obtained by Hod who studied tails in the presence of subleading terms in the potential see subgroup IIIb in [1]. However, Hod s analysis, restricted to the first-order approximation, was inconclusive because, as we just have shown, without the second-order formula 7 one is not in position to make assertions about the dominant behavior of the tail. V. NUMERICS In order to verify the above analytic predictions we solved numerically the initial value problem 4-5 for various potentials initial data. Our numerical algorithm is based on the method of lines with finite differencing in space explicit fourth-order accurate Runge-Kutta time stepping. As was pointed out in [3], a reliable numerical computation of tails requires high-order finite-difference schemes, since otherwise the ghost potentials generated by discretization errors produce artificial tails which might mask the genuine behavior. The minimal order of spatial finite-difference operators depends on the fall-off of the potential for the cases presented below the fourth-order accuracy was sufficient, but for the faster decaying potentials a higher-order accuracy is needed. To eliminate highfrequency numerical instabilities we added a small amount of Kreiss-Oliger artificial dissipation All computations were performed using quadruple precision which was essential in suppressing round-off errors at late times. The numerical results presented here were produced for initial data of the form φ0,r = exp r, t φ0,r = exp r. 9

6 6 As follows from 9 the generating function for these data is We considered the following potentials ax = l+ 1 xexp x, hence A = axdx = π/ l+. 30 tanh α+ r r Vr = α type I 31a tanh α+ r r α 1+ tanhγ r r γ type II, 31b for various values of α γ. The regularizing factor tanhr introduces exponentially decaying corrections to the pure inverse-power behavior at infinity but such corrections do not affect the polynomial tails. The numerical verification of the formulae 15, 8, is shown in tables II III. The observed decay rates agree perfectly with analytic predictions, while small errors in the amplitudes are due to neglected higher-order terms in the perturbation expansion. TABLE II: The generic case: numerical verification of the analytic formula 15 for the potential 31a λ = 0.1 initial data 9. Comparing the second column of this table corresponding to α = 3.01 with the last column of table III one can see the discontinuity of the decay rate at α = 3 for d = 5 7. α =.5 α = 3.01 α = 4 Theory Numerics Theory Numerics Theory Numerics d = 3 d = 5 d = 7 Exponent Amplitude Exponent Amplitude Exponent Amplitude TABLE III: The exceptional case: comparison of analytic numerical parameters of the tails for the potential 31b the first two columns 31a the third column with α = 3, λ = 0.1, initial data 9. The analytic results are given by the formula 8 for the subtype IIa potential, by the formula 7 for the type I IIb potentials. Note that although the dominant tails for the type I the subtype IIb potentials are theoretically the same, in the case IIb there is an additional first order error due to the subdominant term Oλt l+α+γ which accounts for a slight difference in numerical accuracy between these two cases. γ = 1/ subtype IIa γ = 1.75 subtype IIb type I Theory Numerics Theory Numerics Theory Numerics d = 5 d = 7 d = 9 Exponent Amplitude Exponent Amplitude Exponent Amplitude

7 7 VI. SCHWARZSCHILD BACKGROUND Consider the evolution of the massless scalar field outside the d + 1 dimensional Schwarzschild black hole ds = 1 1 r d dt r d dr +r dω d 1, 3 where dω d 1 is the round metric on the unit sphere Sd 1 d 5 is odd. Here we use units in which the horizon radius is at r = 1. Introducing the tortoise coordinate x, defined by dr/dx = 1 1/r d, decomposing the scalar field into multipoles, one obtains the following reduced wave equation for the jth multipole [13] tψ xψ +Uxψ = 0, U = 1 1 j +d 3j +d 1 r d 4r + d 1 4r d. 33 Note that 33 is the 1+1 dimensional wave equation on the whole axis < x <. For large positive x we have r = x+ 1 1 d 3x d 3 d 1 1 +O d 5d 3 xd 5 x 3d 7, 34 which implies that j +d 3j +d 1 Ux = 4x +Vx, Vx = a x d + b 1 +O xd with d 1jj +d a = d 3 x 3d 4 as x, 35 b = d 3d 3d d 1 4jj +d 1+dd 3 4d 5d For large negative x near the horizon the potential is exponentially small, so one expects that the backscattering off the left edge of the potential can be neglected. If so, the decay rate but not the amplitude! should follow from the analysis of section 4. Comparing equation 33 for large positive x to equation 18 with the potential 5 using 35 we find that l = j +d 3/ the potential V is of the subtype IIa with α = d γ = d. Thus, applying 8 we get the first-order tail ψt,x t j+3d Remark 3. Late-time tails outside higher dimensional Schwarzschild black holes were studied in [9], however in the even-dimensional case the reasoning presented there is not correct, even though the result agrees with 37. The reason is that the analysis of [9] is based on the application of Ching et al. conjecture about the decay of the second-order tail t l+α which for l = j+d 3/ α = d gives t j+3d 5. Unfortunately, this conjecture does not apply to the problem at h. For j = 0 this is evident because the leading term in V proportional to x d vanishes since by 36 a = 0, while the subleading term proportional to x d is of generic type. For j > 0 this follows from the fact that the potential is of the subtype IIa. Thus, for all j 0 the dominant first-order contribution to the tail comes from the subleading term in the potential. The agreement of the decay rate obtained in [9] with 37 is accidental due to the fact that the subdominant term in 35 not considered in [9] is on a borderline between subtypes IIa IIb. Admittedly, the hwaving argument leading to 37 is far from satisfactory. Unfortunately, we have not been able to carry over the analysis from sections -4 in the case of equation 33. There are two difficulties in this respect. First, in contrast to the spherical case, Huygens principle is not valid for the free wave equation in 1+1 dimensions. Second, there is no natural small parameter in the problem. In the impressive tour de force work [7] Barack showed how to overcome these difficulties for a restricted class of initial data in four dimensions. It would be interesting to generalize Barack s approach to higher even-dimensional Schwarzschild spacetimes.

8 8 Acknowledgments: PB thanks Nikodem Szpak for helpful discussions Leor Barack for clarifying some details of the paper [7]. AR thanks Prof. Bernd Brügmann for hospitality in his group at FSU Jena, where a part of this work was done. This research was supported in part by the MNII grant 1PO3B019 grant 189/6. PR UE/007/7. APPENDIX A Throughout the appendix we use the notation of [14] in which the square bracket around a logical expression returns a value 1 if the expression is true a value 0 if the expression is false: 1 if condition = true [condition] = 0 if condition = false In order to derive the asymptotic behavior of the iterates near timelike infinity fixed r t we need to evaluate the following expression where Ft,r; β, L = β 4r l+1 L k=0 c L,k dη t+r t r dξ P l µ ξ η β+l kak η, A1 c L,k = L k! k!l k! A µ = ξ t+r rξ η. A3 From 9 13 we have φ 1 t,r = Ft,r; α, l, A4 from 19 1 we have φ t,r = l α/+1/ 1 α 1r α+l l α/+1/ q q α/ q α q Ft,r; α 1+q, l 1 q. A5 q=0 Since aη has compact support, it is advantageous to begin with integration by parts dη P l µ ξ η β+l kak η = dη 1 k dk dη k P l µ ξ η β+l k aη. For µ as defined in A3 for any function gµ the following identity holds d k gµ k dη k ξ η β = k r β +k 1 k j t ξ j g j µ, A6 j r ξ η β+k+j hence Ft,r; β, L = β 4r l+1 dη aη j=0 0 j k L 1 k k c L,k β +L 1 k j 1 t+r r t ξ j j r j dξ ξ η β+l+j P j l µ. A7 t r

9 9 The sum over k can be evaluated explicitly L 1 k k L k! j k!l k! β +L 1k j = 1 L L j Let us define k=j I := 1 r j t+r t r dξ r t ξ j ξ η β+l+j P j l µ. β L j. A8 A9 Changing the integration variable from ξ to µ integrating by parts, we get Using the identity [15] I = 1 j rj+1 β +L j [ r ] β 1+L µ k = +1 1 l=k,k,k 4,... dµp l µ dj dµ j k l/ k l ] exping [1 µ dj β +L j dµ j j 1 r µ in Taylor series we get I = 1 j rj+1 β +L j [ r ] β 1+L j+n/ m=0 j m = rl+1 β +L l [ r ] β 1+L j+l+n/ m=0 j m β +L n=0 β +L j j +n m = rl+1 β +L l [ r ] β 1+L j+l+n/ m=0 j m [1 µ j 1 r ] β +L j µ. A10 l+1k! P l µ, A11!k +l+1!! +1 j +n j dµp l µµ n β +L l/ n=0 β +L j j +l+n m β +L l/ n=0 β +L j j +l+n m 1 1 j+n+m r 1 j+n m +1 j +l+n j dµp l µµ l+n 1 l+m r Collecting the results of A8, A10, A1 plugging them into A7 we get Ft,r; β, L = β+l+1 4 where Sβ,L = L β j +l+n L j j j=0 β +L l dη aη [ r ] β 1+L j+l+n/ m=0 j+n m j +l+n j l+1l+n!l+n! n!l+n+1! j+n m r 1 l+m. A1 β +L l/ n=0 1 m j β +L j m j +l+n m l +n!l+n! n!l+n+1! 1L+l L!Sβ,L, A13 j+n m r. A14

10 10 1. First-order approximation To evaluate the first iterate φ 1 t,r we apply the formula A13 with β = α L = l. Then l Sα,l = α l+n+j j+n 1 j+n m j α +l j m r, A15 l j j j +n m l j +m j=0 m=j l/ where we shifted the summation index m j +n m. Next, we interchange the order of summation according to [0 j][j l][m n j][j m+l] [ l m < 0][0 j l+m] + [0 m < n][0 j l] + [n m l+n][m n j l], convert the sum over j into a generalized hypergeometric function [14]. Defining t j = 1 j+n m α l+n+j j α +l j, l j j j +n m l j +m we see that t 0 0 iff n = m, thus the sums for [ l m < 0] [0 m < n] do not contribute to A15 we are left with Sα,l = n+l m l+n m r j=0 1 j α l+n m j j +m n α +l+n m j j l+n+m j l+n+m+j j +m n, A16 where we shifted the summation index j j +m n. Defining t j = 1 j α l+n+m+j j +m n α +l+n m j l+n m j j +m n j l+n+m j, we see that hence t 0 = α l+n m l+n m! l+n+m! m n!l+n! α +l+n ml+n+m l+n+m! t j+1 j l+n mj l+n+mj +l +n+m+1 = t j j +α 1 l+n mj + α l+n mj +1, Sα,l = n+l m r α l+n m α +l+n ml+n+m l+n m! m n!l +n! l+n m, l+n+m, l+n+m+1 F α 1 l+n m, α l+n m 1 = n+l m r 1+l+n m π α l+n m α +l+n ml+n+m l+n m! m n!l +n! α l +n mα 1 l+n m α 3 +m α l+n α +m α 1 l+n, A17

11 11 where in the last equation we used the identity a+1, a, b+c 1/ F b, c 1 = b+c π Substituting into A17 we get Sα,l = n+l b a c a bc 1+b+a α l+n m α 1 l+n m = α 1, 1+c+a. α +l+n m l+n+m α l+n m = 1 l+n+m α++m m r 1 l+n+m 1+l+n m π l+n m!m n!l +n! The last equation can be still simplified due to the identity We have α α 3 α 1 α++m α +m α 3 +m α l+n α 1 l+n. α 1 α+ α l α 1 l = 1l π α++m = α+ m α+, α 3 +m = α 3 m α 3 α l n α 1 l n α +m = α l α l 1 n, = α 1 l α 1 l 1 n, = α m α, lα α 3 l. A18, so finally Sα,l = n+l α+ m α 3 m = m α m, +1 m r 1 n+m l+n l+n m!m n!l +n! lα α 3 α +1 m α l 1 n α 1 l 1 n α m. A19 l Plugging A19 into A13 with β = α L = l we get the expression 14.

12 1. Second-order approximation The calculation in the second order β = α 1+q L = l 1 q is only a slight modification of what we have already done in the first order. Following the same steps which led us from A15 to A18 we get n+l m r 1 l+n+m 1+L+n m π Sβ,L = L+n m!m n!l +n! β + l L +m β 3 + l L +m The last equation can be simplified due to the identity β + l L β 3 + l L β 3 l L = 1l π which for L = l reduces to A18. We have hence Sβ,L = n+l β 1 β ++l L+m β l+l +n β 1 β ++l L β L β +l L β 1. A0 l+l +n β 1 L β l L, A1 β ++l L+m = β ++l L m β ++l L, β 3 + l L +m = β 3 + l L m β 3 + l L, β β n = n, β 1 β 1 β 1 n = n, β 1 1 β + l L β +m = + l L m β + l L, β ++l L m m = m β β 3 + l L + l L m +1, m r 1 n+m L+n L+n m!m n!l +n! L L β 3 l L β +l L β β l L + l L +1 Plugging A into A13 we get m β β + l L α 5 l 1 q α +l l 1 q α 3 1+q Ft,r; α 1+q, l 1 q = 1 q α+3l q 4 α 4 dη aη [ r ] α 3+l α 1 n n l+n! α 3 n n+l 1 q n!l+n+1! α+1 ln l+q 1 m l 1 q m n n=0 A n n 1 β 1 1 m. α+m α+q m A3 m r.

13 Substituting this into A5 exping in 1/t we have 1 φ t,r = α 1 α+l 4l+1!! l α 1/ 1 q l p q q α/ q q=0 [ 1 t α+l A+α+l 1 B t α q α 5 1 +O t ] l 1 q α +l l 1 q α 3 1+q 13, A4 with A B defined in 17. Converting the sum over q into the generalized hypergeometric function we get. [1] W. Strauss K. Tsutaya, Discrete Cont. Dyn. Sys. 3, [] N. Szpak, arxiv: [math-ph] [3] E. S. C. Ching et al., Phys. Rev. D5, [4] R. Price, Phys. Rev. D5, [5] E. W. Leaver, Phys. Rev. D34, [6] C. Gundlach, R. Price J. Pullin, Phys. Rev. D49, [7] L. Barack, Phys. Rev. D59, [8] M. Dafermos I. Rodnianski, Invent. Math. 16, [9] V. Cardoso et al., Phys. Rev. D68, [10] J. G. Kingston, Quart. Appl. Math. 46, [11] H. Lindblad C. D. Sogge, Amer. J. Math. 118, [1] S. Hod, Class. Quantum Grav. 18, [13] A. Ishibashi H. Kodama, Prog. Theor. Phys. 110, [14] R. L. Graham, D. E. Knuth O. Patashnik, Concrete Mathematics Reading, Massachusetts: Addison-Wesley, [15]