Computig the Iverse Melli Trasform of Holoomic Sequeces usig Kovacic s Algorithm arxiv:8.9v [cs.sc] 2 Ja 28 Research Istitute for Symbolic Computatio RISC) Johaes Kepler Uiversity Liz, Alteberger Straße 9, A- Liz, Austria E-mail: jablige@risc.jku.at We describe how the extesio of a solver for liear differetial equatios by Kovacic s algorithm helps to improve a method to compute the iverse Melli trasform of holoomic sequeces. The method is implemeted i the computer algebra package HarmoicSums. th Iteratioal Symposium o Radiative Correctios 2-29 September, 27 St. Gilge, Austria Speaker. This work has bee supported i part by the Austria Sciece Fud FWF) grat SFB F5 F59-N5). c Copyright owed by the authors) uder the terms of the Creative Commos Attributio-NoCommercial-NoDerivatives. Iteratioal Licese CC BY-NC-ND.). https://pos.sissa.it/
. Itroductio There have bee several methods proposed to compute the iverse Melli trasform of special sequeces, for istace i [] a algorithm usig rewrite rules) to compute the iverse Melli trasform of harmoic sums was stated. This algorithm was exteded i [] to geeralized harmoic sums such as S-sums ad cyclotomic sums. A differet approach to compute iverse Melli trasforms of biomial sums was described i []. I [2] a method to compute the iverse Melli trasform of geeral holoomic sequeces was described. That method uses holoomic closure properties ad was implemeted i the computer algebra packageharmoicsums [,, 5,, 7]. I the frame of the method a liear differetial equatio has to be solved. So far the differetial equatios solver of HarmoicSums was oly able to fid d Alembertia solutios [8]. Recetly the solver was geeralized ad therefore more geeral iverse Melli trasforms ca be computed. I the followig we repeat importat defiitios ad properties compare [2,, ]). Let K be a field of characteristic. A fuctio f = fx) is called holoomic or D-fiite) if there exist polyomials p d x), p d x),..., p x) K[x] ot all p i beig ) such that the followig holoomic differetial equatio holds: p d x) f d) x)+ + p x) f x)+ p x) fx) =..) We emphasize that the class of holoomic fuctios is rather large due to its closure properties. Namely, if we are give two such differetial equatios that cotai holoomic fuctios fx) ad gx) as solutios, oe ca compute holoomic differetial equatios that cotai fx) + gx), fx)gx) or x fy)dy as solutios. I other words ay compositio of these operatios over kow holoomic fuctios fx) ad gx) is agai a holoomic fuctio hx). I particular, if for the ier buildig blocks fx) ad gx) the holoomic differetial equatios are give, also the holoomic differetial equatio of hx) ca be computed. Of special importace is the coectio to recurrece relatios. A sequece f ) with f K is called holoomic or P-fiite) if there exist polyomials p d ), p d ),..., p ) K[] ot all p i beig ) such that a holoomic recurrece p d ) f +d + + p ) f + + p ) f =.2) holds for all N from a certai poit o). I the followig we utilize the fact that holoomic fuctios are precisely the geeratig fuctios of holoomic sequeces: if fx) is holoomic, the the coefficiets f of the formal power series expasio fx)= = form a holoomic sequece. Coversely, for a give holoomic sequece f ), the fuctio defied by the above sum i.e., its geeratig fuctio) is holoomic this is true i the sese of formal power series, eve if the sum has a zero radius of covergece). Note that give a holoomic differetial equatio for a holoomic fuctio fx) it is straightforward to costruct a holoomic recurrece for the coefficiets of its power series expasio. For a recet overview of this holoomic machiery ad further literature we refer to []. A additioal property of f x
holoomic fuctios was give for example i [2] ad []: if the Melli trasform of a holoomic fuctio M[ fx)]) := x fx)dx=f)..) of a holoomic fuctio is defied i.e., the itegral x fx)dx exists, the it is a holoomic sequece. Coversely, if the Melli trasform M[ fx)]) of a fuctio fx) is holoomic, the also the fuctio fx) is holoomic. I this article we will report a extesio of HarmoicSums based o Kovacic s algorithm [2] that supports the user to calculate the iverse Melli trasform i terms of iterated itegrals that exceed the class of d Alembertia solutios. The paper is orgaized as follows. I Sectio 2 we revisit the method to compute the iverse Melli trasform of holoomic fuctios from [2], while i Sectio we explai the geeralizatio of the method ad i Sectio we give some examples. 2. The Iverse Melli Trasform of Holoomic Sequeces I the followig, we deal with the followig problem: Give a holoomic sequece F). Fid, wheever possible, a holoomic fuctio fx) such that for all N from a certai poit o) we have M[ fx)]) = F). I [2] a procedure was described to compute a differetial equatio for fx) give a holoomic recurrece for M[ fx)]). Give this procedure the followig method to compute the iverse Melli trasform of holoomic sequeces was proposed i [2]:. Compute a holoomic recurrece for M[ fx)]). 2. Use the method metioed above to compute a holoomic differetial equatio for fx).. Compute a liear idepedet set of solutios of the holoomic differetial equatio usig HarmoicSums).. Compute iitial values for M[ fx)]). 5. Combie the iitial values ad the solutios to get a closed form represetatio for fx). I our applicatios we usually apply this method o expressios i terms of ested sums, however as log as there is a method to compute the holoomic recurrece for a give expressio i.e., item ca be performed) this proposed method ca be used. Aother possible iput would be a holoomic recurrece together with sufficiet iitial values. Note that util recetly HarmoicSums could fid all solutios of holoomic differetial equatios that ca be expressed i terms of iterated itegrals over hyperexpoetial alphabets [, 9,, ]; these solutios are called d Alembertia solutios [8]. Hece as log as such solutios were sufficiet to solve the differetial equatio i item we succeeded to compute fx). I case d Alembertia solutios do ot suffice to solve the differetial equatio i item, we have to exted the solver for differetial equatios. 2
. Beyod d Alembertia solutios of liear differetial equatios Util recetly oly d Alembertia solutios of liear differetial equatios could be foud i usig HarmoicSums compare [2]), but i order to treat more geeral problems the differetial equatio solver had to be exteded. I [2] a algorithm to solve secod order liear homogeeous differetial equatios is described. We will refer to this algorithm as Kovacic s algorithm. Cosider the holoomic differetial equatio p i x) C[x]) Kovacic s algorithm decides whether.) p 2 x) f x)+ p x) f x)+ p x) fx)=..) has a solutio of the form e ω where ω Cx); has a solutio of the form e ω where ω is algebraic over Cx) of degree 2 ad the previous case does ot hold; all solutios are algebraic overcx) ad the previous cases do ot hold; has o such solutios; ad fids the solutios if they exist. Note that the solutios Kovacic s algorithm ca fid are called Liouvillia solutios []. I case Kovacic s algorithm fids a solutio, it is straightforward to compute a secod solutio which will agai be Liouvillia. This algorithm was implemeted iharmoicsums. Example. Cosider the followig differetial equatio: x 89x+7x 2) + 9x 2 2 7x+59x 2) D x + 8x x+27x 2) ) D 2 x fx)=, with the implemetatio of Kovacic s algorithm i HarmoicSums we fid the followig two solutios: x x 2+2 x x f x) =, x) x 5 +27x) + x x 2 2 x x f 2 x) =. x) x 5 +27x). Composig solutios Suppose we are give the liear differetial equatio q i, p i C[x];d > ) ) q d x)d d x + + q x) fx) =,.2) which factorizes liearly ito d first-order factors. The this yields d liearly idepedet solutios of the form f2 x) f x), f x) f x) dx, f f2 x) f x) x) f x) f 2 x) dxdx,..., f f2 f d x) x) f x) f d x) dx dx,
where the f i are hyperexpoetial fuctios i.e., D x f i x) f i x) Kx) ). These solutios are also called d Alembertia solutios of.2), compare [, 8]. Now suppose that a give differetial equatio does ot factorize liearly, but cotais i betwee secod-order factors, which ca be solved, e.g., by Kovacic s algorithm. Let the followig differetial equatio correspod to a secod order factor: p2 x)d 2 x+p x)d x + p x) ) fx) =,.) the we ca compose the solutios of the first order ad secod order factors as follows. Let sx) be solutio of.2) ad let g x) ad g 2 x) be solutios of.). The g x) sx), sx) sx) dx ad sx) g2 x) sx) dx are solutios of p2 x)d 2 x + p x)d x + p x) ) q d x)d d x x)) + + q fx)=. I additio, if we defie wx) := p 2 x)g x)g 2x) g x)g 2 x)) the g x), g 2 x) ad g x) sx)wx)g 2 x)dx g 2 x) sx)wx)g x)dx are solutios of q d x)d d x + + q x)) p2 x)d 2 x+p x)d x + p x) ) fx)=.. Examples Example 2. We wat to compute the iverse Melli trasform of ) ) f :=. 27 We fid that which leads to the differetial equatio 2+)+2) f + 9+)2+) f + =, 27x ) fx)+9x7x ) f x)+8x 2 x ) f x)=, for which we fid with the help of Kovacic s algorithm the geeral solutio 2+2 x x 2 2 x x sx)= c + c 2 x x 2, x x 2 for some costats c ad c 2. I order to determie these costats we compute x sx)dx = 9 c + 8π ) + 8log2) + 9 c 2 + 8π ) 8log2),
x 2 sx)dx = 8 c 7+ 2π ) + 2log2) + 8 c 2 7+ 2π ) 2log2). Sice f = /9 ad f 2 = 8/2 we ca deduce that c = c 2 = f = [ 2 2 x x+ π M xx 2/ i= π ad hece ] 2+2 x x Example. Durig the Computatio of the iverse Melli trasform of ) i i i we have to solve the followig differetial equatio: = 2727x ) fx)+ x 2 2x+ ) f x)+9x 5x 2 298x+ ) f x) +8x )x 2 27x ) f ) x). We are able to fid the geeral solutio of that differetial equatio usigharmoicsums: x c 27x + c 2 Give this geeral solutio we fid: ) ) i 27 + = i i i= τ) 2/ + τ τ dτ π π 27x ) x 27 ) x + c x ) x 27 27x dx 27 x x ). τ+ τ) 2/ dτ τ 27x ) 2/ dτ τ + τ τ 27x dx ) 2/ dτ + τ τ τ ) 27x Fially, we list several examples that could be computed usigharmoicsums: Example. ) = 2 2 2π ) = 2 2 ) = 27) x + x+ x ) x+ dx, x xx x + x + x ) x + x xx dx, x + x + x ) ) 2/ x + dx, x xx dx.. 5
i= i= i= ) i = i ) i 2i i ) i i i 2 = 27 ) + π 2 = + x x π x 27 ) + π x ) ) ) 27 2 2 x x+ 2+2 x x ) x 27 27x x ) x 27 27x logx)dx log x ) dx, x27x ) + y + ) y x y + dydx y yy / x dx y y 27 z+ + z dzdydx, zz 2/ ) x ) ) 27 27x dx. Ackowledgemets I wat to thak C. Scheider, C. Raab ad J. Blümlei for useful discussios. Refereces [] J. Abliger, The package HarmoicSums: Computer Algebra ad Aalytic aspects of Nested Sums, PoS LL 2 2) 9 [arxiv:7.8 [cs.sc]]. [2] J. Abliger, Iverse Melli Trasform of Holoomic Sequeces, PoS LL 2 2) 7 [arxiv:.285 [cs.sc]]. [] J. Abliger, Computer Algebra Algorithms for Special Fuctios i Particle Physics, PhD Thesis, J. Kepler Uiversity Liz, 22. arxiv:5.87 [math-ph]. [] J. Abliger, J. Blümlei, C. G. Raab ad C. Scheider, Iterated Biomial Sums ad their Associated Iterated Itegrals, J. Math. Phys. 55 2) 2 doi:./.98 [arxiv:7.822 [hep-th]]. [5] J. Abliger, J. Blümlei ad C. Scheider, Harmoic Sums ad Polylogarithms Geerated by Cyclotomic Polyomials, J. Math. Phys. 52 2) 2 doi:./.2972 [arxiv:5. [math-ph]]. [] J. Abliger, J. Blümlei ad C. Scheider, Aalytic ad Algorithmic Aspects of Geeralized Harmoic Sums ad Polylogarithms, J. Math. Phys. 5 2) 82 doi:./.87 [arxiv:2.78 [math-ph]]. [7] J. Abliger, A Computer Algebra Toolbox for Harmoic Sums Related to Particle Physics, arxiv:.7 [math-ph]. [8] S.A. Abramov ad M. Petkovšek, D Alembertia solutios of liear differetial ad differece equatios, i proceedigs of ISSAC 9, ACM Press, 99. [9] M. Brostei, Liear Ordiary Differetial Equatios: breakig through the order 2 barrier, i proceedigs of ISSAC 92, ACM Press, 992. [] P.A. Hedriks ad M.F. Siger, Solvig differece equatios i fiite terms, J. Symbolic Comput. 27, 999. [] M. Kauers ad P. Paule, The Cocrete Tetrahedro, Text ad Moographs i Symbolic Computatio, Spriger, Wie, 2.
[2] J.J. Kovacic, A algorithm for solvig secod order liear homogeeous differetial equatios, J. Symbolic Comput. 2, 98. [] M. Petkovšek, Hypergeometric solutios of liear recurreces with polyomial coefficiets, J. Symbolic Comput., 992. [] E. Remiddi ad J.A.M. Vermasere, Harmoic Polylogarithms, It. J. Mod. Phys. A5, 2. [arxiv:hep-ph/99527v] 7