Indefinite Integras of Spherica Besse Functions MIT-CTP/487 arxiv:703.0648v [math.ca] 6 Mar 07 Joyon K. Boomfied,, Stephen H. P. Face,, and Zander Moss, Center for Theoretica Physics, Laboratory for Nucear Science, and Department of Physics, Massachusetts Institute of Technoogy, Cambridge, MA 039, USA Dated: March, 07 Highy osciatory integras, such as those invoving Besse functions, are best evauated anayticay as much as possibe, as numerica errors can be difficut to contro. We investigate indefinite integras invoving monomias in mutipying one or two spherica Besse functions of the first kind j with integer order. Cosed-form soutions are presented where possibe, and recursion reations are deveoped that are guaranteed to reduce a integras in this cass to cosed-form soutions. These resuts aow for definite integras over spherica Besse functions to be computed quicky and accuratey. For competeness, we aso present our resuts in terms of ordinary Besse functions, but in genera, the recursion reations do not terminate. I. INTRODUCTION In a mode of eary universe physics, we encountered integras of the form d fj α and d fj k αj β when attempting to compute the covariance between the spherica harmonic modes of a fied. Here, j is the spherica Besse function of the first kind, and we were interested in integers k, 0, and reas α, β. Our functions f were sowy-varying functions of that had been numericay computed at discrete vaues. As a piecewise-poynomia interpoation of f is an appropriate description of such a function, it was natura to consider indefinite integras of the form d n j α and d n j k αj β to numericay compute our integras. Because the spherica Besse functions can osciate widy when α and/or β are arge, it was important that these integras be computed anayticay. We searched the iterature for integration methods for integrands invoving spherica Besse functions. We found a number of resuts regarding definite integras over an infinite range [ 4], but very itte regarding indefinite integras. Resuts for Besse functions with integer order appear to be much more commonpace [5]. NIST s Digita Library of Mathematica Functions [6] was unhepfu, and even the venerabe Gradshteyn and Rhyzik [7] provided itte assistance. Wofram Research [8] had some of the integras that we wanted and Mathematica coud integrate a few more, but essentiay a these resuts were written in terms of generaized hypergeometric p F q functions. As these genera functions are not avaiabe in typica scientific ibraries, we set out to derive our own resuts. This paper presents anaytic resuts for indefinite integras of the form with integer n, integers k, 0 and reas α and β. The resuts have been constructed with an eye towards numerica impementation in a anguage ike C or Python, and hence are written ony in terms of specia functions that are commony avaiabe in scientific ibraries. For some epressions, we find cosed-form resuts, whie for others, we derive recursion reations that decrease the the order of spherica Besse functions appearing in the integras. For joyon@mit.edu face@mit.edu zander@mit.edu A number of our resuts ikey appy more generay, but we do not investigate this.
the cass of integras we consider, a recursion reations terminate with cosed-form soutions. This ony occurs for spherica Besse functions with integer order, as integras invoving j 0 = sinc are anayticay tractabe unike integras invoving J 0. To our knowedge, the recursion reations we derive are new, as are some of the cosed-form soutions that we present. We begin this paper by describing we-known properties of spherica Besse functions. We then construct some preiminary quantities which appear in the genera integras, before treating the singe and product integras in turn. An appendi transates appicabe resuts to ordinary Besse functions. II. OVERVIEW OF SPHERICAL BESSEL FUNCTIONS We begin by presenting we-known resuts regarding spherica Besse functions. Spherica Besse functions are the soutions to the differentia equation d y dy + d d + [ + ]y = 0. 3 They are indeed by the order. Two independent soutions are j and y, the spherica Besse functions of the first and second kind, respectivey. These functions are reated to ordinary Besse functions by π π π j = J +/ and y = Y +/ = + J /. 4 The spherica Besses can be computed from Rayeigh s formuas j = d sin and y = d cos d d. 5 For integer, the spherica Besse functions of the first and second kind are reated to each other by π j + = J / = + y. 6 In this paper, we concern ourseves ony with the spherica Besse functions of the first kind with integer order 0. The first few j are as foows. j 0 = sin j = sin j = 3 = sinc 7a cos sin 3 cos This pattern continues for higher orders, and a j with natura can be written as a sum of sines and cosines mutipied by inverse powers of. We wi restrict ourseves to positive arguments, athough j is defined over the entire compe pane. Negative arguments of j are reated to positive arguments by 7b 7c j = j. 8 A more genera reation appies for compe or non-integer. Spherica Besse functions satisfy two important recursion reations which can be combined to construct further reations. j = j j 9a j = j j 9b
Spherica Besse functions satisfy a cosure reation 0 j kj k d = where δ is the Dirac deta function, as we as an orthogonaity reation j k j d = π k δk k 0 π + δ k for k, N, where δ k is the Kronecker deta. A number of infinite integras over spherica Besse functions are known [3]. π Γ + j d = 0 Γ + j π d = 3 + 0 0 j kj k π k < d = + k > + 4 Here, k < and k > are the smaer and arger of k and k respectivey. The recursion reationships described in the rest of this paper ece when performing integras over many osciations. If one is trying to perform integras over few if any osciations at ow, quadrature methods wi ikey be better than the recursion methods, which can deveop formay canceing divergences at = 0 this is a we-known issue in evauating Besse function recursion reations [5]. From inear regression on 0 00, the first zero of j occurs roughy around 4.75 +.05. 5 Knowedge of the ocation of the first zero of the integrand can hep in determining which computationa method to use. Note that this is an overestimate for very ow ; j 0 has its first zero at = π. It is aso worth knowing the behavior of j near = 0. At = 0, j 0 0 =, whie j = 0 for > 0. The series epansion of j about = 0 is given by π j = + Γ + 3 + O. 6 III. PRELIMINARIES Because spherica Besse functions can be written in terms of sines and cosines divided by powers of their argument, we wi often encounter integras of the form X n = d n sin 7 Y n = d n cos 8 for integer n. It wi be usefu to be abe to compute these efficienty. Integrating by parts, we obtain the foowing two recursion reations. X n = ny n n cos 9 Y n = n sin nx n 0 3
If n > 0, these reations aow us to step down by a singe n, terminating at n = 0 with X 0 = cos and Y 0 = sin. For n < 0, we can invert these recursion reations to aow us to step upwards in n. These recursions terminate at n = with the specia functions X n = n + [n+ sin Y n+ ] Y n = n + [n+ cos + X n+ ] 3 X = Si = Y = Ci = sint dt t 4 cost dt t 5 0 where Si is the sine integra, and Ci is the cosine integra. These specia functions are avaiabe in scientific ibraries such as SciPy and GSL. We now ook cosey at = 0. We have X 0 0 = and Y 0 0 = 0. These can be used to obtain X n 0 = Y n 0 = 0 for a n > 0. Looking towards negative n, we have X 0 = 0, and X n 0 is divergent for n <. Y n 0 is divergent for a n < 0. In the foowing sections, we wi encounter epressions of the form X n α/α n+ and Y n α/α n+. In the regime where α π, we recommend using numerica quadrature. For n 0, one can aso turn to series epansions. X n α α n+ = Y n α α n+ = d n Γn + cosnπ/ sinα = α n+ + α n+ + n α 64 + n + α 4 4 06 + n + Oα6 6 d n Γn + sinnπ/ cosα = α n+ + n+ + n α 3 + n + α4 4 45 + n + Oα6 Note the presence of constant terms that turn on and off with n. In the case of n =, we suggest using epansions for the sine and cosine integras. For ower n, we were unabe to find appropriate epansions, and urge care if numerica resuts are required in this regime. 7 IV. INTEGRALS OF A SINGLE SPHERICAL BESSEL We now turn towards integras of the form I n ; α = d n j α 8 with n Z, N, and α R. The specia case of α = 0 simpifies the integra drasticay, and so we assume α 0. A substitution then yieds I n ; α = α n I n α 9 where we define I n = d n j. 30 4
Using the sma behavior of j, we see that I n is finite at = 0 for + n >. To evauate I n, we begin by ooking at = 0. Because j 0 = sin/, we immediatey have I0 n = d n sin = X n. 3 Of particuar note is I 0 0 = Si. For > 0, we use the recursion reationship to compute I n = j = j j 3 [ ] d n j j = + n I n n j. 33 This aows to decrease by one each step, recursing unti = 0 with I0 n X n. The recursion wi run through, n pairs as i, n i for integer i where 0 i. There are + steps, assuming the recursion doesn t truncate from the + n term vanishing. To truncate, we must have + n odd. The recursion wi then truncate when i = + n /. Due to the range of i, this further requires n +. Mathematica wi integrate I n directy, epressing the resut in terms of a hypergeometric F function. For specia choices of parameters, the hypergeometric function yieds simper functions. These resuts can be directy verified by differentiation and use of the spherica Besse recursion reations 9. I + = + j + 34 I 4+ = 3+ [ + 3j + j +3 ] 35 π I = Γ + j 36 Generay speaking, n = + Ω and n = + Ω for positive even Ω yied resuts simiar to these, but the epressions quicky become more compicated as Ω increases. V. INTEGRALS OF A SPHERICAL BESSEL SQUARED Let us now ook towards the foowing integra. H n ; α = d n j α 37 Again, α = 0 drasticay simpifies the integra, and so we assume α 0. We again perform a substitution to obtain H n ; α = α n H n α 38 where we have defined H n = d n j. 39 From a series epansion of the integrand, H n is finite at = 0 for + n >. 5
We once again start with the = 0 case, writing j 0 = sin/. By further writing sin = cos/, we obtain H0 n = d n j 0 = d n n Y n. 40 If n, this becomes whie for n =, we instead have This is finite in the imit of sma, H n 0 = n n n Y n, 4 H0 = [n Ci]. 4 H 0 = γ + n + 4 + O6 43 where γ is the Euer-Mascheroni constant. We now turn to the case of > 0. Note the foowing reationship, which can be verified by using the recursion reationships 9. [ n n j n j j ] = n + n 3n j + n j n j Performing an indefinite integra over eads to n n j n j j = n + n 3Hn + Hn H n 45 where constants of integration are incuded in the H functions. This resut can be rearranged to sove for H n. H n = H n + n + n 3Hn + n n j n j j 46 This aows us to recurse to ower vaues, which can be repeated unti we get to H0 n, which has been evauated above. Mathematica wi compute H n directy in terms of a hypergeometric F 3 function. For specia vaues of the parameters, the hypergeometric function simpifies to simper functions. Once again, these resuts can be directy verified by differentiation and use of the spherica Besse recursion reations 9. H = H = H = H 4 = H +3 = d j = j j j + j + d j π = 4 + Γ + j + j 4 48 j j j + 49 d j = 3 d 4 j = + 34 + j j + 4 + + 3j + 4 + + + 3j d +3 j = + j + j + 5 4 + 44 47 50 6
Further cosed-form epressions eist for n = Ω, n = + Ω and n = + + Ω with positive even Ω, athough the compeity rapidy increases with Ω. It is usefu to note that these anaytic resuts often aow for an eary termination to the recursion reationships, without having to recurse a the way down to = 0. VI. INTEGRALS OF TWO SPHERICAL BESSELS, SAME ORDER We now turn to integras of two spherica functions that have the same order. K n ; α, β = d n j αj β 5 If α or β = 0, the integra simpifies significanty, whie if α = β, we reduce to a previous case. Hence, we assume α β 0. Note that K n is symmetric in α and β. As for H n, Kn ; α, β is finite at = 0 for + n >. We once again begin with = 0. Writing j 0 = sin/, we obtain K0 n ; α, β = d n sinα sinβ. 53 αβ Noting that the integra can be evauated to yied K0 n ; α, β = [ αβ sinα sinβ = [cosα β cosα + β], 54 α β n Y n α β If α β, then care must be taken in computing these quantities see Section III. For > 0, we again construct a derivative reationship. [ ] n αβ [ nj αj β βj αj β + αj αj β] ] α + β n Y n α + β. 55 = n j αj β αβ n n + 3n j αj β α + β αβ n j αj β 56 This can be derived using the recursion reations 9. Performing an indefinite integra over, we obtain n αβ [ nj αj β βj αj β + αj αj β] = K n ; α, β n n + 3Kn αβ ; α, β α + β αβ Kn ; α, β. 57 Once again, constants of integration are absorbed into the definition of K. Soving for K n ; α, β, we obtain the foowing recursion reation. K n ; α, β = [ α + β K n αβ ; α, β + n n + 3K n ; α, β ] + n n j αj β n [βj αj β + αj αj β] 58 In the imit α = β =, this reduces to the recursion reation 46 for H n. The ony cosed-form soution we coud find for integras of this type is K ; α, β = d j αj β = α β [βj αj β αj αj β]. 59 Once again, this can be checked by differentiation and appication of the recursion reations 9. 7
VII. INTEGRALS OF TWO SPHERICAL BESSELS, DIFFERENT ORDER We now consider our fina integra, where we aow for the order of the Besse functions to be different. L n k; α, β = d n j k αj β 60 If α = 0 or β = 0, the integra once again simpifies drasticay, and so we assume that α and β are nonzero rea numbers. From series epansions of the integrand, L n k ; α, β is finite at = 0 for k + + n >. Our approach in this case is to reduce the arger of the two orders through repeated appication of the recursion reations 9. Assuming that k <, we obtain [ ] L n k; α, β = d n j k α β j β j β = L n β k ; α, β Ln k ; α, β. 6 We can use this repeatedy unti either = k, in which case the K integras appy, or = k +, which we now evauate epicity. Consider the resut n+ j αj β = n n j αj β + β n+ j αj β α n+ j αj β 6 which can be verified by differentiation and appication of the recursion reations 9. Integrating this over yieds n+ j αj β = n L n,; α, β + βk n+ ; α, β αkn+ ; α, β 63 where constants of integration are absorbed into the functions. Rearranging, we then find L n,; α, β = n+ j αj β + αk n+ ; α, β βk n+ n ; α, β. 64 The recursion reations may aso terminate at L n 0; α, β = L n 0; β, α, which we now compute. L n 0; α, β = d n sinα sinβ α β cosβ β = [ ] αβ d n 3 cosα β cosα + β [ ] d n sinα β + sinα + β αβ = [ ] αβ α β n Y n 3α β α + β n Y n 3α + β [ ] αβ α β n X n α β + α + β n X n α + β 65 66 67 If α β, then these quantities wi again need to be computed carefuy see Section III. We were ony abe to find one cosed-form reation for genera L n k ; α, β. α β L k; α, β + [ + kk + ]L 0 k; α, β = d [α β + + kk + ]j k αj β 68 = β j k αj β α j k αj β + k j k αj β 69 Unfortunatey, this is not a particuary usefu formua uness the eact combination depicted here arises. However, specia imits such as α = β or k = do yied usefu resuts, as shown above and beow. 8
If we restrict ourseves to the case of α = β, much more can be derived. Let us begin by defining L n k = d n j k j 70 so that L n k; α, α = α n L n kα. 7 Note that L n k is symmetric under interchange of k and. Again, the basic idea wi be to use a recursion reation to decrease the highest of k or by using Eq. 6. In the case of k =, instead of Eq. 64, we have the simper resut L n = Hn If needed, L n 0 can be computed by L n 0 = n 3 = d = d n 3 d n j j = sin d n sin n 3 cos + n H n n j. 7 cos n sin 73 74 n Y n 3 n X n. 75 L n k can be epressed directy in terms of a rather messy hypergeometric 3 F 4 function. We found that certain specia choices of parameters ead to simpifications in the hypergeometric function that aow the resut to be written quite ceany. As aways, these resuts can be verified by differentiation and appication of recursion reations 9. L 0 k = d j k j = kk + + [j k j j k j + kj k j ] L k = L k k = L k+ k = L k++3 k = d j kj = j kj k + j k+j + k k + j kj + + kk + j k+ j + + j k j + k + + k + + k + k d k π j k j = k++ k + Γ k + Γ + k j k j + j k j k + d k+ j k j = k++3 j k j + j k j k d k++3 j k j = k++4 j k+ j + + j k j k + + Further cosed-form epressions are avaiabe for n = k + Ω, n = k + Ω, n = k + + + Ω and n = k + Ω with positive even Ω, athough once again, the epressions quicky become compicated as Ω increases. 76 77 78 79 80 9
VIII. CONCLUSIONS We have investigated indefinite integras of powers of mutipying one or two spherica Besse functions of the first kind with integer order 0. In the process, we presented a number of cosed-form indefinite integras and constructed recursion reations to aow the numerica computation of a the integras we considered. The recursion techniques that we have presented are highy compementary to numerica quadrature methods. When integrating spherica Besse functions over many osciations, our epressions aow for accurate computations in situations where quadrature methods suffer. On the other hand, our recursion reations have numerica issues at sma, where forma divergences are canceing, which ruins numerica precision. In this regime, there are few if any osciations and quadrature methods ece. In the regime in which we don t need to worry about canceing divergences, we have compared a numerica impementation of our resuts to Mathematica s high-precision numerica integration, with eceent resuts. ACKNOWLEDGMENTS This work is supported in part by the U.S. Department of Energy under grant Contract Number DE- SC00567, and in part by MIT s Undergraduate Research Opportunities Program. Appendi A: Resuts transated into ordinary Besse functions In this appendi, we transate our resuts into epressions invoving ordinary Besse functions. Whie the resuts transate perfecty we, the recursion reations typicay do not terminate. Apart from the recursion reations, a number of these resuts appear in the iterature, e.g., [9].. Singe Besse integras The foowing are transations of Eqs. 33, 34, 35 and 36 respectivey. d n J = + n d + J = + J + d n J n J d +3 J = + [ + J + J +3 ] d J = Γ J A A A3 A4. Besse Squared integras The foowing are transations of Eqs. 46, 47, 48, 49, 50 and 5 respectivey. 0
d n J = d n J + n + n 3 d n J + n n J n J J d J = J 4 J J + + J 4 d J = 4 Γ J + J A7 d J = J J J + A8 d 3 J = + 4 4 + J J + + J 6 + + + J A9 d + J = A5 A6 + J + J + A0 + 3. Two Besses, Same Order The foowing are transations of Eqs. 58 and 59 respectivey. αβ d n J αj β = α + β d n J αj β + n n + 3 d n J αj β + n n J αj β n [βj αj β + αj αj β] A d J αj β = α β [βj αj β αj αj β] A 4. Two Besses, Different Order The foowing are transations of Eqs. 6, 64 and 69 respectivey. n α β d n J k αj β = β d n J k αj β d n J k αj β d n+ J αj β d n J αj β = n+ J αj β + α β d n+ J αj β d J k αj β = k d J kαj β + βj k αj β αj k αj β + k J k αj β A3 A4 A5
The foowing are transations of Eqs. 7, 76, 77, 78, 79 and 80 respectivey. d n J J = + n d n J n J d J kj = k [J k J J k J ] + kj k J d J kj = J kj k + J k+ J + k k + J k J + + kk + [J k+ J + + J k J ] + k + k + + + k + k d k J k J = k+ k + ΓkΓ k J k J + J k J k + d k+ J k J = k+ J k J + J k J k d k++ J k J = k++ J k+ J + + J k J k + + A6 A7 A8 A9 A0 A [] R. Mehrem, J. T. Londergan, and M. H. Macfarane, Anaytic epressions for integras of products of spherica besse functions, Journa of Physics A: Mathematica and Genera 4 99 no. 7, 435. http://stacks.iop.org/0305-4470/4/i=7/a=08. [] L. C. Maimon, On the evauation of the integra over the product of two spherica besse functions, Journa of Mathematica Physics 3 99 no. 3, 64 648. [3] R. Mehrem, The Pane Wave Epansion, Infinite Integras and Identities invoving Spherica Besse Functions, Submitted to: App. Math. Comput. 009, arxiv:0909.0494 [math-ph]. [4] V. I. Fabrikant, Eementary eact evauation of infinite integras of the product of severa spherica besse functions, power and eponentia, Quart. App. Math. 7 03 573 58. [5] W. Rosenheinrich, Tabes of some indefinite integras of besse functions,. http://www.eah-jena.de/rsh/forschung/stoer/besint.pdf. [6] NIST Digita Library of Mathematica Functions. http://dmf.nist.gov/, Reease.0.4 of 06--. http://dmf.nist.gov/. F. W. J. Over, A. B. Ode Daahuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Cark, B. R. Mier and B. V. Saunders, eds. [7] I. S. Gradshteyn, I. M. Ryzhik, A. Jeffrey, and D. Zwiinger, Tabe of integras, series, and products. Esevier, Amsterdam, Boston, Paris, et a., 007. [8] Wofram Research, Spherica besse j integras. Website, 007. http://functions.wofram.com/besse-typefunctions/sphericabessej//0/. [9] Wofram Research, Besse j integras. Website, 007. http://functions.wofram.com/besse-typefunctions/bessej//0/.