arxiv: v1 [hep-th] 16 Feb 2018

Transcription

1 Fast Summation of Diverent Series and Resurent Transseries in Quantum Field Theories from Meijer-G Approximants arxiv: v1 [hep-th] 16 Feb 2018 Héctor Mera, 1, 2, 3, Thomas G. Pedersen, 2, 3 and Branislav K. Nikolić 1 1 Department of Physics and Astronomy, University of Delaware, Newark, DE , USA 2 Department of Physics and Nanotechnoloy, Aalbor University, DK-9220 Aalbor Øst, Denmark 3 Center for Nanostructured Graphene CNG, DK-9220 Aalbor Øst, Denmark Dated: February 19, 2018 We demonstrate that a Meijer-G-function-based resummation approach can be successfully applied to approximate the Borel sum of diverent series, and thus to approximate the Borel-Écalle summation of resurent transseries in quantum field theory QFT. The proposed method is shown to vastly outperform the conventional Borel-Padé and Borel-Padé-Écalle summation methods. The resultin Meijer-G approximants are easily parameterized by means of a hypereometric ansatz and can be thouht of as a eneralization to arbitrary order of the Borel-Hypereometric method [Mera et al. Phys. Rev. Lett. 115, ]. Here we illustrate the ability of this technique in various examples from QFT, traditionally employed as benchmark models for resummation, such as: 0-dimensional φ 4 theory, φ 4 with deenerate minima, self-interactin QFT in 0-dimensions, and the computation of one- and two-instanton contributions in the quantum-mechanical double-well problem. I. INTRODUCTION Perturbative expansions in quantum mechanics, quantum field theory QFT and strin field theory often have zero radius of converence, i.e., they are asymptotic[1, 2]. Optimal truncation of such series to a small number of terms can provide experimentally relevant results but only at sufficiently small couplin constant. This is exemplified by the asymptotic series for Stark effect of atoms and molecules in an external electric field where first and second order terms match measurements well, but only for very weak electric fields [3], or by hih precision calculations of multiloop Feynman diarams in quantum electrodynamics [4] where the fine-structure constant is small. On the other hand, extractin physically relevant information from asymptotic series at larer couplin constants calls, almost invariably, for resummation techniques as exemplified by the Stark effect [5, 6] and field assisted excitonic ionization in novel materials [7], anharmonic oscillators in quantum mechanics [8, 9], φ 4 theory in QFT [10], quantum chromodynamics [11], strin perturbation theory [12, 13], and diarammatic Monte Carlo in condensed matter physics [14]. Given the ubiquity of diverent series in physics, research on summation techniques remains an active research area [15, 16]. Conventional resummation is, however, not sufficient in the presence of the so-called Stokes phenomenon where different asymptotic expansions hold in different reions of the plane made up of complexified expansion parameter values [17 23]. Thus the Stokes phenomenon requires enerally distinct resummations in each of these reions. This complexity is captured by resurent transseries [21, 24 27], which include both analytic polynomial terms and nonanalytic exponential and hypereometric2f1@mail.com loarithmic terms. In principle, resurent transseries offer a non-perturbative framework to reconstruct the oriinal function, which has led to recent viorous efforts to examine their promise in physically relevant examples, where different sectors are enerated by non-perturbative semiclassical effects like instantons[22, 23, 27]. Nonetheless, the resurent transseries also need to be resummed in order to obtain a sensible result. However, the conventional Borel-Padé-Écalle resummation used for this purpose typically requires a lare number of terms [23, 28] from each sector in order to obtain reasonably accurate results beyond the weak couplin reime. This makes it useless for problems in QFT [4, 10] or many-body perturbation theory in condensed matter physics [14] where 10 orders are available at best. Here we introduce a new alorithm which replaces the standard Padé approximants in the Borel plane by more eneral and flexible hypereometric functions of, in principle, arbitrarily hih order, thereby achievin reat converence acceleration towards the exact sum of a resurent transseries. Hence, in this work, we replace the conventional Borel-Padé-Écalle summation by Borel-Hypereometric-Écalle summation whose approximants admit a representation in terms of Meijer- G functions which are easily parameterized. Prior to oin into technical details, we hihliht the power of our alorithm by notin that our approximants convere to the exact non-perturbative, ambiuity-free, partition function for 0-dimensional self-interactin QFT with just five orders, whereas Borel-Padé-Écalle resummation in Ref. 28 needed tens of terms to find a ood approximation at intermediate couplin strenth. In practice only a handful of expansion coefficients are typically available: for instance the ɛ-expansion for the On- symmetric φ 4 theory is only known to six-loop order [10], while the five-loop QCD beta function and anomalous dimension have been calculated only very recently [29]. It

2 2 is then clear that, for a summation technique to be practical, it needs to be able to return accurate estimates of the sum of a diverent series with only a few coefficients. The hih-accuracy at low orders of the Meijer-G approximants introduced in this work makes them suitable for practical applications. The paper is oranized as follows. In Sec. II we introduce the alorithm to transform the expansion coefficients of a diverent expansion into a table of Meijer- G approximants. In Section III we apply our alorithm to sum four well-known examples in QFT: 0- dimensional φ 4 theory, φ 4 theory with deenerate minima, self-interactin 0-dimensional QFT and finally the computation of the one- and two-instanton contributions in the double-well problem of quantum mechanics. While in the first of these examples we deal with a Borelsummable diverent expansion, the other three cases offer an opportunity to demonstrate the efficiency of our approach for the summation of transseries/multi-instanton expansions, i.e., those cases where the perturbation expansion is not Borel-summable. The Meijer-G summation method is shown to work well in all of these cases, providin a fast way to evaluatin the Borel sum of a diverent series and massively outperformin the Borel- Padé and Borel-Padé-Écalle approach. In the first three examples the Meijer-G approximants convere to the exact result at five-loop order, while in the latter case the converence is slower althouh fast when compared to Borel-Padé-Écalle approaches. Finally in Sec. IV we discuss the advantaes and disadvantaes of our approach as an alternative to traditional Borel-Padé techniques. We end with the conclusions in Sec. V. II. MEIJER-G APPROXIMANTS An efficient resummation technique should be capable of takin a handful of coefficients and return an accurate estimate of the sum of a diverent perturbation series. In this Section such a technique is introduced, providin an alorithm that efficiently transforms the coefficients of a diverent perturbation expansion into a table of Meijer-G functions which serve as approximants. For completeness we first briefly review the traditional Borel- Padé resummation, emphasisin the inherent difficulties faced by such method when it comes to summation on the cut. We will then briefly review alternatives to Padé and Borel-Padé that make use of analytic continuation functions with a built-in branch cut, in particular the hypereometric resummation method we introduced in Refs After reviewin these other approaches we finally introduce the alorithm to calculate Meijer-G approximants and hihliht various properties that make this method extremely well-suited to yield inexpensive and yet accurate low order approximations to the sum of a diverent series. A. Borel-Padé Resummation When it comes to sum diverent series, Borel-Padé has become the dominant approach [1, 15, 30]. There are various reasons for the popularity of this approach but it can be arued that the most important of these is its alorithmic simplicity. The Borel-Padé approach is in essence a simple recipe to transform the coefficients of a diverent expansion into a table of approximants, which approximate the Borel sum of a diverent series and are typically evaluated by numerical contour interation. Another advantae of this approach is that it relies on widely studied approaches: conditions for Borelsummability are by now well understood [30 33] and the properties of Padé approximants used at a crucial step in the alorithm are also very well known as they have been studied in depth for decades [34]. Given a diverent expansion Z n=0 z n n, where z n are the expansion coefficients and is the expansion parameter the couplin, the alorithm to calculate the Borel sum of a diverent series can be summarized as follows: Step 1: Borel transform. Calculate the Boreltransformed coefficients: b n = z n /n!. Step 2: Summation in the Borel plane. Sum the series Bτ n=0 b nτ n ; this series is called the Boreltransformed series. The complex-τ plane is known as the Borel plane. Step 3: Laplace transform. The Borel sum of the series, Z B is iven by the Laplace transform Z B = 0 e τ Bτdτ. 1 The rationale behind the Borel summation method is simple: the coefficients of a series with zero radius of converence typically row factorially at lare orders. In the first step such factorial rowth is removed and the Borel-transformed series is more tractable since it has a finite non-zero radius of converence. By summin the Borel-transformed series one finds the function Bτ and then the Laplace transform can be calculated numerically to find the Borel sum Z B. The Borel-Padé summation method is a practical alorithm to find, in principle, increasinly accurate approximations to the Borel sum. In this approach the second step above is specialized to Padé summation: Step 2: Padé summation in the Borel plane. Use Padé approximants to approximately sum the series Bτ n=0 b nτ n, from the knowlede of the partial sums to order N, i.e., N n=0 b nτ n. In Padé summation one approximates Bτ by a rational function of τ B L/M τ = L n=0 p nτ n 1 + M n=1 q nτ n. 2 where L + M = N. Here the coefficients p n and q n are found by equatin order by order the Taylor series of B L/M τ around τ = 0 with the asymptotic expansion

3 3 of Bτ, up to the desired order. These approximants are then used in step 3 of the alorithm to evaluate the Laplace transform and therefore to find the L/M-Borel- Padé approximant to the Borel sum of Z, Z B,L/M, as Z B,L/M = 0 e τ B L/M τdτ. 3 Clearly the Borel-Padé method returns a table of approximants: for instance by knowin the partial sums to second order N = 2 one can calculate B 1/0, B 0/1, B 2/0, B 1/1, and B 0/2. There are, however, various types of problems for which Borel-Padé approximants are not well suited. Padé approximants are rational functions whose built-in sinularities are poles. It is well-known that in many situations the converence of the Borel-transformed series is limited by the branch point sinularity closest to the oriin, the terminal point of a dense line of poles known as a branch cut. In such scenarios Borel-Padé approximants can convere very slowly since many poles may be needed to properly mimic a branch cut. Thus, to overcome this difficulty with Borel-Padé one has to replace Padé approximants, adoptin instead as approximants functions that are able to mimic branch cuts in the Borel plane. The advantaes of that stratey are hihlihted in Fi. 1 where we show a domain colorin plot of a function Bτ in the Borel plane correspondin to the partition function of 0-dimensional φ 4 for = 1. We show the 20/20 Padé approximant, B 20/20 τ top panel toether with the exact Bτ bottom panel. We also show a hypereometric 2 F 1 approximant middle panel in the Borel plane which corresponds to the third order Meijer-G approximant, introduced in Sec. II C. The Padé approximant in the Borel plane attempts to reproduce the branch cut by brute force placin poles next to each other alon the neative τ axis. The Meijer-G approximant has a builtin branch cut in the Borel plane and it can thus very accurately mimic the branch cut usin only three orders of perturbation theory. One has to carefully look at the the details in Fi. 1 to see the very minor differences between the exact Bτ and its third-order hypereometric approximant. B. Hypereometric Resummation Recently we have introduced hypereometric resummation a technique that enables summation on the cut usin only a small number of expansion coefficients [5 7, 35, 36]. Various flavours of this technique were applied to a variety of problems with ood results: in particular it was shown how one could use low order data to derive accurate approximations to the decay rate in Stark-type problems [5 7]. Typically the idea is to use hypereometric functions to analytically continue diverent series. Oriinally the technique was devised by noticin the shortcomins of Padé approximants when FIG. 1. Domain colorin plot of approximants in the Borel plane calculated for the partition function of φ 4 theory in zero dimensions for = 1. Branch cuts are represented by a discontinuous chane of color. Around poles one sees radual colour chanes, intersectin black lines and concentric white lines. Top panel: 20/20 Padé approximant. Middle panel: hypereometric 2F 1 approximant. Bottom panel: Exact. The Padé approximant mimics the exact branch cut by placin poles next to each other and hih orders are needed to accurately model the cut. In contrast, the third order 2F 1 approximant with a built-in branch cut very nearly reproduces the exact branch cut. applied to problems in nonequilibrium nanoelectronics modellin within the nonequilibrium Green s function theory [37 40]. As pointed out in Sec. II A one seeks to substitute Padé approximants by more eneral functions that are equipped with branch cuts and can work well in those cases where Padé approximants do not work or convere very slowly. Hypereometric functions looked particularly promisin as they are endowed with a branch cut and eneralize eometric series which can be summed exactly by Padé approximants. For instance iven a diverent function Z n=0 z n n, with normalized zeroth order coefficient z 0 = 1, one can attempt to find a hypereometric function[41, 42], 2F 1 h 1, h 2, h 3, h 4, such that Z = 2 F 1 h 1, h 2, h 3, h 4 + O 5. 4 Similar ideas were already considered in Refs. 30, 43 45, where the authors used it as part of an alorithm to find converent stron couplin expansions from diverent weak couplin ones. Similarly, in Refs. 46 and 47, the authors used products of hypereometric functions as approximants. Foreshadowin all of these works are

4 4 contributions by Stilliner and coworkers [48, 49], where early Borel-Hypereometric approximants typically 1 F 0 hypereometric functions in the Borel plane are considered, as well as early versions of self-similar factor approximants and exponential-borel approximants. Of the works mentioned, the work of Stilliner is the closest in philosophy to our own work, as he also considers the ratio test of series converence as a startin point as we do in Sec. III C and we did in earlier work [35]. The hypereometric approximants we introduced in Refs. 5 7, and 35 have a number of clear limitations, which we overcome in this work and which we enumerate below: 1. Hypereometric resummation is uncontrolled: Hypereometric 2 F 1 approximants of third, fourth and fifth order can be constructed in various ways, but in Refs. 5 7, and 35 we did not ive a recipe for constructin a table of hypereometric approximants. How does one parameterize and build, say, a 21st order hypereometric approximant? In order to have control over any approximations one develops, it is fundamental to be able to increase the order of the approximation and to study the converence properties of the approximants. In this work we provide a set of approximants that can, in principle, be computed at any order. 2. Difficult parameterization at lare orders: A natural attempt to provide a eneralization to arbitrary orders would be to state that eneral hypereometric functions q F p [41, 42] constitute the approximant space. These functions contain p + q + 1 parameters, p F q h 1,..., h p ; h p+1,..., h p+q, h p+q+1 τ, that need to be calculated by equatin each order of the asymptotic expansion that one seeks to sum with the correspondin order of the expansion of the hypereometric approximant. However one faces a deeneracy problem as all the hypereometric functions obtained by permutin elements of each of the parameter sets h 1,..., h p and h p+1,..., h p+q are one and the same. Hence when determinin the parameters h i the computational time rows factorially with order there is a factorially lare number of solutions, all of which correspond to the same hypereometric function. So calculations at the six-loop order are already very expensive and nearly impossible at the 7th loop order. The same problem is found with other approximants; for instance, self-similar factor approximants [46, 47] have rarely been computed beyond the sixth loop order. In this work we put forward an alorithm that enables fast parameterization of hypereometric approximants in the Borel plane, at arbitrarily hih orders. 3. Inaccuracies for expansions with zero radius of converence: It was noted that for very small couplins the hypereometric approximants as described in Ref. 5 ave inaccurate results. While these inaccuracies were exponentially small, they were conceptually important. The reason for them was that the radius of converence of hypereometric p+1 F p functions is not zero, while we were applyin the hypereometric approximants to problems with zero radius of converence. For instance an approximant iven by Eq. 4 above has radius of converence c = 1/h 4. When applied to series with zero radius of converence the value of h 4 was typically found to be very lare but finite. Therefore the hypereometric approximants had a Taylor series with a tiny, but non-zero, radius of converence. Hence, in Refs. 6 and 7 we used a different parameterization of the approximants, which positioned the tip of the hypereometric branch cut exactly at the oriin and therefore alleviated this problem. However there we did not come up with a clear approach to compute similar approximants of hiher order. The approximants derived in this work can have zero radius of converence, therefore alleviatin this difficulty found in our previous approach. Clearly there is a wide variety of problems for which Borel-Padé resummation can be improved. However attempts to improve it can easily fall into severall pitfalls. In the case of hypereometric resummation these were difficulties in both extendin the approach to arbitrarily hih orders and dealin with series with zero radius of converence. These difficulties are larely surpassed by the Meijer-G approximants we introduce next. C. Meijer-G Resummation We now present an alorithm to transform the low order coefficients of a diverent perturbation expansion, Z n=0 z n n with normalized coefficients, z 0 = 1, into a table of approximants to its Borel sum. The alorithm consists of four easy steps. For ease of presentation, in this work we will compute mostly odd-order approximants, ivin a short description of the alorithm to compute even-order approximants below. Step 1: Borel transform. Imaine that we know only N coefficients, z 0, z 1,..., z N. In this step we compute the Borel-transformed coefficients b n = z n /n! toether N ratios of consecutive Borel-transformed coefficients, rn = b n+1 /b n. Here we will assume that N is an odd number. Step 2: Hypereometric ansatz. We make the ansatz that rn is a rational function of n. Thus we define a rational function of n, r N n as r N n = l m=0 p mn m 1 + l m=1 q mn m, 5 where l = N 1/2 and the N unknown parameters p m and q m are uniquely determined by the N input ratios

5 by means of N equations, rn = b n+1 b n = r N n, 0 n N 1. 6 It should be noted that this is a system of N linear equations with N unkowns q m and p m which can be easily solved by a computer. The hypereometric ansatz is the crucial step that allows an extremely fast parameterization of lare order Meijer-G approximants. Step 3: Hypereometric approximants in the Borel plane. In this step we undertake the parameterization of hypereometric approximants in the Borel plane. To do this for N > 1, we use the calculated p m and q m to set two equations 1 + l p m x m = 0, m=0 l q m y m = 0, 7 m=1 which yield two solution vectors x 1,..., x l and y 1,..., y l. We refer to these vectors as hypereometric vectors. It follows from the definition of hypereometric functions that the hypereometric vectors determined in this way uniquely determine the hypereometric function B N τ l+1 F l x, y, p l q l τ, 8 where x = 1, x 1,..., x l, y = y 1,..., y l and l+1f l is a eneralized hypereometric function [41]. The function B N τ is the hypereometric approximant in the Borel plane; it provides an N-th order approximation to the sum of the Borel-transformed series. Thanks to the hypereometric ansatz and the hypereometric vector equations, the function B N τ can be easily parameterized for arbitrary lare N. Step 4: Meijer-G approximants. In this last step we need to reinstate the n! removed from the expansion coefficients by means of the Borel transform. This is achieved by means of the Laplace transform Z B,N 0 e τ B N τdτ, 9 which ives the desired approximantion to the Borel sum of the asymptotic expansion of Z in the complexified -plane. This expression admits the representation Z B,N = Πl i=1 Γ y i Π l i=1 Γ x i G l+2,1 l+1,l+2 1, y 1,..., y l 1,1, x 1,..., x l q l, p l 10 where Γx is Euler s Gamma function and Gp,q m,n a1,...,a p b 1,...,b q z is Meijer s G-function MGF [41, 50, 51]. This alorithm then transforms N available input coefficients z n into a table of Meijer-G functions, which approximate the Borel sum of Z. We once aain emphasize how easily one can parameterize these extremely complex functions: all what was needed was the hypereometric ansatz and the resultin hypereometric vectors. Once these are determined the Meijer-G approximants can be parameterized for arbitrarily lare orders. There are various remarks that we would like to add before movin onto the practical application of this method. The even order approximants can be computed in exactly the same way, but one starts from a oncesubtracted series, for instance from the once-subtracted Borel-transformed series n=0 b nτ n 1/b 1 τ and the B N are calculated as above but with l = N/2. The Laplace transform then ives the even order approximants. It is also easy to compute approximations to the Mitta-Leffler ML sum [52], which is a eneralization of the Borel sum where one can remove superfactorial asymptotic coefficient rowth as opposed to factorial by usin a ML transform of the coefficients as opposed to a Borel transform; in step 1 the ML-transformed coefficients are b n = a n /Γαn + β, where α and β are real numbers, and the odd approximants in step 4 are iven by f N λ = e τ τ β 1+α B 0 N λτ α dτ. When α = β = 1 ML summation is equivalent to Borel summation. These approximants also admit a MGF representation not shown. Furthermore, further eneralizations of Borel and ML summation methods, where essentially arbitrary asymptotic coefficient rowth is removed can be easily arrived at. It follows from the differential equation satisfied by MGFs that the Meijer-G approximants iven by Eq. 10 provide a reularizin analytic continuation of the diverent hypereometric functions [50]. Such hypereometric functions have zero radius of converence; the fact that Meijer-G approximants are able to sum them clearly illustrates the potential of these functions for the summation of diverent series. l+2f l y 1,..., y l ; 1, 1, x 1,..., x l ; p lλ q l III. APPLICATIONS IN QUANTUM FIELD THEORY In this section Meijer-G approximants are used to sum partition functions in QFT usin only a few orders of perturbation theory. In particular, in Sec. III A we will consider the summation of the partition function in φ 4 theory, which is Borel summable. We will also consider the summation of non-borel summable series by considerin φ 4 theory with deenerate minima [53] in Sec. III B, as well as a self-interactin QFT model in zero dimensions [28] in Sec. III C. In such cases the standard perturbation series needs to be upraded to a transseries. These examples are often utilized for benchmarkin new summation techniques, since they contain hihly diverent series which are very difficult to sum. It will be shown that, remarkably, these cases constitute a bestcase scenario for the application of the Meijer-G summation technique developed above. This is so because the 5

6 6 Meijer-G approximants convere in all of these cases at order N = 5. This means that by applyin the Meijer-G summation procedure we arrived at a closedform analytic expression; in other words, these partition functions belon to the space of functions reproducible by our summation technique and thus they are easily summable by means of Meijer-G approximants. The reason for this is that in these cases the hypereometric ansatz turns out to be exact, and the Borel-transformed series sums exactly to a hypereometric function of the form 3 F 2 1, h 1, h 2 ; h 3, h 4, h 5 τ [or, equivalently, to a hypereometric function of the form 2 F 1 h 1, h 2, h 3, h 4 τ]. Therefore, in Sec. III D we complement our study by considerin the summation of the transseries expansion for the double-well potential in quantum mechanics, includin the one- and two-instanton contributions [27], which constitutes a more challenin example since the approximants exhibit hihly non-trivial converence properties. A. Partition function in φ 4 Theory The partition function in zero-dimensional φ 4 theory is iven by Z = 1 e φ2 /2 φ 4 /4! dφ, 11 2π which for Re[] > 0 can be written as 3 Z = 2π e 3 4 K , 12 where K ν x is a modified Bessel function of the second kind. This partition function is commonly used to benchmark resummation techniques see Refs. 15 and 54 for two recent examples. The first few terms of the asymptotic expansion about = 0 read Z The expansion coefficients row factorially at lare orders and thus this expansion has zero radius of converence. Here the calculation of Z by direct resummation of the asymptotic expansion serves two main purposes. On the one hand it provides a simple test system for benchmarkin aainst Borel-Padé. On the other hand the parameterization of the approximants can be performed analytically at low orders by followin our alorithm and thus is valuable from a didactical perspective. Let us calculate the first-order and third-order Meijer G approximant for this problem analytically, by runnin the above-iven alorithm explicitely. We start with the first order calculation and proceed step by step. 1. Borel transform: we calculate the Borel transformed coefficients b n = z n /n!. In a first-order calculation we have just two coefficients, b 0 = 1 and b 1 = 1/8. 2. Hypereometric ansatz: we compute the ratios between consecutive coefficients as a function of n, rn and we approximate the ratio by a rational function of n. In this case we have only one ratio r0 = b 1 /b 0 = 1/8 and the rational function that approximates rn is just a constant: rn r0, n. 3. Hypereometric approximant in the Borel plane: in this case the hypereometric vectors are empty, x 1,..., x l = {} and y 1,..., y l = {} and therefore x = 1 and y = {}. Since l = N 1/2 = 0 for N = 1, the first order hypereometric approximant in the Borel plane is iven by B H,N=1 τ = 1 F 0 1, r0τ. 14 This 1 F 0 hypereometric function is just the 0/1 Padé approximant to the Borel-transformed series, i.e., B H,1 τ = r0τ. 15 It should be clear that the hypereometric vectors are empty, and that the above-iven hypereometric function in the Borel plane reduces to the eometric case. Therefore our first-order hypereometric approximant coincides with the 0/1 Padé approximant. 4. Meijer-G approximant: once the hypereometric approximant in the Borel plane is found, we can use Eq. 10 to immediately write down the correspondin Meijer-G approximant in the complexified- plane, which reads Z B,1 = 1 8 G 2,1 1,2 0 0, This Meijer-G approximant is just the 0/1 Borel- Padé approximant. This is an interestin aspect of hypereometric and Borel-hypereometric resummation: to first order hypereometric approximants are just 0/1 Padé approximants; translatin this observation to the Borel plane shows that the correspondin Meijer-G Borelhypereometric approximant is just the 0/1 Borel- Padé approximant. Indeed, e τ Z B,1 = dτ = r r0τ U 1, 1, r0, 17 is both the 0/1 Borel-Padé approximant and the first order Meijer-G approximant. The conclusion from this first example calculation is that our first-order hypereometric-borel Meijer-G approximant is the first-order Borel-Padé approximant. This illustrates a eneral property of Meijer-G approximations

7 7 to the Borel sum: their first order is just identical to a first-order Borel-Padé approximant. The interested reader can now look to the closely related approach put forward in Ref. 15 both approaches take very different routes beyond first order. Next we run aain our alorithm to obtain the third order approximant. While the first order Meijer-G approximant is identical to the first order Borel-Padé approximant we will see shortly that the third order Meijer-G approximant is substantially more accurate than Borel- Pade approximants of the same and much hiher orders. 1. Borel transform. In this case we have four Boreltransformed coefficients: b 0 = 1, b 1 = 1/8, b 2 = 35/768 and b 3 = 385/ Hypereometric Ansatz. Here we have three ratios: r0 = 1/8, r1 = 35/96, and r3 = 11/24. We approximate rn as rn = r 3 n where r 3 n = p 0 + p 1 n 1 + q 1 n, 18 and use the known ratios, r0, r1 and r2 to find p 0, p 1 and q 1 by requirin rn = r 3 n, n = 0, 1, 2, 19 which leads to three equations r0 = 1 8 = p 0, 20 r1 = = p 0 + p q 1, 21 r2 = = p 0 + 2p q 1, 22 with three unknowns, p 0, p 1 and q 1. These equations yield a solution p 0 = 1 8, p 1 = , q 1 = Therefore our third order rational approximation to the ratios rn reads rn r 3 n = 1/8 113n/ n/9 3. Hypereometric Approximant in the Borel plane. To build the hypereometric approximants in the Borel plane we need to find the hypereometric vectors, that is to find the values of x and y that solve these two equations p 0 + p 1 x = 0, q 1 y = 0, 26 which yield x 1 = p 0 /p 1 and y 1 = 1/q 1 and thus find the vectors x = 1, p 0 /p 1 and y = 1/q Z Exact BP5/5 BP10/10 BP20/20 Meijer-G N=3 Meijer-G N= FIG. 2. Z for 0-dimensional φ 4 theory calculated usin MGFs and Borel-Padé approximants. For lare values of the third-order Meijer-G approximant filled inverted trianles and the fifth order Mejer-G approximant empty circles are compared with hiher-order 5/5 empty trianles, 10/10 filled diamonds and 20/20 empty squares Borel-Padé approximants and with the exact Z solid line with filled circles. The third order Meijer-G approximant is more accurate that 5/5 and 10/10 Borel-Padé approximants, and slihtly less accurate than the 20/20 Borel-Padé approximant. The fifthorder Meijer-G approximant is exact and thus more accurate than any Borel-Padé approximant. Hence the third order hypereometric approximant in the Borel plane is B H,3 τ = 2 F 1 1, p 0 p 1, 1 q 1, p 1 q 1 τ. 27 Substitutin values of p i and q i we et B H,3 τ = 2 F 1 1, , 9 7, 113 τ Meijer-G approximant. We read-off the Meijer-G approximant directly from the hypereometric vectors to et a third order approximant that can be compactly written as Z B,3 = Γ9/7 Γ27/113 G 3,1 2,3 1,9/7 1,1,27/ One can easily calculate hiher order Meijer-G approximants. It turns out that the fifth order Meijer-G approximant is convered and equal to the exact Borel sum, i.e., Z B,5 = Z, 30 and all hiher order Meijer-G approximants are also equal to Z, i.e., Z B,5 = Z B,7 = = Z. 31

8 8 Method = 1 = 10 = 100 BP2/ i i i BP5/ i i i BP10/ i i i BP20/ i i i Meijer-G N = i i i Meijer-G N = i i i Exact i i i TABLE I. Z for 0-dimensional φ 4 theory evaluated on the cut, for = 1, = 10 and = 100 usin Borel-Padé approximants of orders 4 BP2/2, 10 BP5/5, 20 BP10/10 and 40 BP20/20, compared with third and fifth order Meijer-G approximants and with the exact value. All resuts are iven to six sinificant diits. The accuracy of the third order Meijer-G approximant is comparable to that of BP approximants of hiher order at weak couplins = 1 and reater at lare couplins = 10 and = 100. The fifth order Meijer-G approximant is exact. What is happenin is that the rational approximations used in the hypereometric ansatz have convered at fifth order, i.e., the ratio between consecutive Boreltransformed coefficients is indeed a rational function of n, and, in fact, of the form specifically rn = p 0 + p 1 n + p 2 n q 1 n + q 2 n 2, 32 rn = 1/8 2n/3 2n2 / n + n 2, 33 which reproduces the ratios between Borel-transformed coefficients up to arbitrarily hih orders. Approximatin these ratios by rational functions of hiher order, such as r 7 n, one finds the same rational function once aain. Hence the approximants of order fiver or hiher are convered. We now compare the performance of Meijer-G approximants with that of Borel-Padé approximants. In Fi. 2 we compare the 5/5, 10/10 and 20/20 Borel-Padé approximants of orders 10, 20 and 40 respectively with the third and fifth order Meijer-G approximant for lare values of. It is clear that the third order Meijer-G approximant is more accurate than the 10/10 Borel-Padé approximant, but less accurate than the 40-order 20/20 Borel-Padé approximant. The fifth order Meijer-G approximant is exact and it is therefore more accurate than any Borel-Padé approximant. It is instructive to compare Z with Z 3 for < 0. For instance while Z 10 + iɛ = i , 34 Z iɛ = i , 35 which demonstrates the reat accuracy of Meijer-G approximants for rather lare neative couplins; in particular the imainary part of Z is reproduced within a percent. In Table I we compare aain the third and fifth order order Meijer-G approximants aainst the Borel- Padé approximants this time for neative couplins. We see that for = 1 all approximants are very accurate. The third order Meijer-G approximant is more accurate than the 2/2 Borel-Padé approximant, but less accurate than the other Borel-Padé approximants shown. For = 10 the Meijer-G approximant is already more accurate than 2/2, 5/5 and 10/10 Borel-Padé approximants. Finally for = 100 the Meijer-G approximant is more accurate than all the Borel-Padé approximants reported. The fifth order Meijer-G approximant reproduces the exact result. These findins demonstrate how to easily build accurate Meijer-G approximations to the sum of a Borelsummable diverent series, and are confirmed by calculations for the quartic anharmonic oscillator and the Heisenber-Euler action for QED, which will be shown elsewhere. These approximants have reat potential for the summation of Borel-summable series and should be accurate in cases where the converence of the Boreltransformed series is limited by a branch point sinularity. Here we have considered the case of 0-dimensional φ 4 theory, which is a challenin example for standard approximants which typically require a lare number of coefficients in order to yield accurate estimates for Z; in contrast Meijer-G approximants convere rapidly to the exact answer. In many cases, however, perturbation expansions are not Borel-summable. In such cases one needs to sum a resurent transseries. Below it is shown that Meijer-G approximants can also be used to provide economical and accurate approximations to the sum of such expansions. B. Deenerate Vacua Marucho [53] considered a partition function of the form Z = 1 e φ2 1 φ 2 /2 dφ, > π

9 9 The term φ 2 /2 in the exponent may be rearded as a freefield = 0 and a traditional perturbative approach, such as diarammatics, results in an asymptotic expansion in powers of. The first few terms of this expansion are Z , 37 and its eneral term can be found in the paper in Ref. 53. Hence we are once aain in a situation where we happen to know all of the coefficients. We stress that this only happens in toy models and expansion coefficients are rarely known at lare orders. With these coefficients we can easily run our alorithm and find the odd order Meijer-G approximants, which read Z B,3 = Z B,5 = G 3,1 2,3 G 4,1 3,4 Z B,7 = G 5,1 4,3. 0, ,0,0 0,0,0 0.25, 0.75,0,0 0,0,0,0 0.25, 0.75,0,0, , , These Meijer-G approximants to the Borel sum are fully convered at order N = 5 and larer; by increasin the order of the approximants we et the same function aain and aain: Z B,5 = Z B,7 =. This means that we have found the exact Borel sum by means of our approach; the Borel sum of Z belons to the space of reproducible functions associated with Borel-Hypereometric resummation. The convered approximants aree with the exact Borel sum [53] iven by. Z B i0.1e / K 1/ x Our approximants actually find the exact Borel sum for this problem with just a few orders. In contrast Borel- Padé would have required very hih order coefficients to produce accurate results up to modest values of the couplin. Unfortunately for the present problem havin very accurate approximations to the Borel sum turns out to be insufficient to accurately approximate Z. The problem is of course that Z is not Borel-summable. This is clearly seen in Fiure 3 which shows the real and imainary parts of Z B,3 and Z B,5 toether with the exact Z. The Borel sum fails to accurately represent Z; in Fi. 3a we see that the real part of the Borel sum is off by a constant factor of 2; in Fi. 3b we see that the Borel sum posseses an imainary part, while Z is manifestly real. Note that the N = 3 approximant is nearly convered while the N = 5 approximant is fully convered to the exact value of the Borel sum, iven in Ref. 53. When performin a Borel summation it is not uncommon to obtain imainary parts. While in some cases it is rather straihtforward to ive them a physical meanin, in the present case it is not so clear and furthermore such physical interpretation is not our objective which is to accurately approximate Z. This imainary part is called non-perturbative ambiuity. It is ambiuous because it chanes sin dependin on whether one evaluates Z B + iɛ or Z B iɛ, i.e. Im[Z B + iɛ] = Im[Z B iɛ]. Non-Borel summability for this problem is explained in Ref. 53. For completeness we ive our version here see Fi. 4. For very small but non-zero couplins the interand in Z contains contributions from two identical aussians, each of which sits on a different saddle, one at x = 0 and the other at x = 1/. In the literature the former is referred to as the perturbative saddle while the latter is known as the nonperturbative saddle. If the area under each aussian is A, then Z 0 + = 2A. When doin perturbation theory around = 0 one is effectively takin into account contributions only from the perturbative saddle and thus the perturbative estimate will be Z 0 + = A, which is wron by a factor of two. The Meijer-G approximants succeed in summin up the perturbation expansion around = 0 from only a few terms, but that is not enouh. The Borel sum misses the contribution from the second aussian entirely and thus underestimates the real part of Z by a factor of two. It is then clear that Z has an essential sinularity at = 0. Perturbation theory is not sufficient to accurately evaluate Z even in the limit 0. In order to obtain an accurate estimate of Z we need to uprade the perturbation expansion to a resurent transseries which, for the present problem, is of the form Z = ±i e F 0 1/4, 3/4, F 0 1/4, 3/4, = ±i e , where the 2 F 0 factors are diverent hypereometric series and the upper sin is for Im > 0, while the lower sin is for Im < 0. This transseries contains two diverent series that need to be summed. The so-called Borel-Padé-Écalle method uses the Borel-Padé technique to sum each of these diverent series [17 20, 22, 28]. Here we use instead Meijer-G approximants to sum each of the diverent series that appear in the transseries expansion.

10 10 are[z] Z Re[Z B,3] Re[Z B,5] b Im[Z] Im[Z B,3] Im[Z B,5] FIG. 3. Partition function, Z, for φ 4 theory with deenerate minima. a The real part of the exact Z solid is compared with the real part of the third- Z B,3, filled circles and fifth-order Z B,5, empty trianles Meijer-G approximants to the Borel sum of Z. The Meijer-G approximants are well convered for the rane of values of shown and are not clearly distinuishable on the scale of the plot. The fifth order approximant is convered to the exact Borel sum iven by Marucho [53]. The real part of the Borel sum underestimates Z by a factor of two, Re[Z B,5] = Z/2. b Imainary parts of the Meijer-G approximants are not zero, while Z is manifestly real, Im[Z B,5] 0 while Im[Z] = 0. Perturbation theory for φ 4 with deenerate minima is not Borel-summable. e - = A A a Z Z Re[Z B,1,3]+Re[Z B,2,3] Re[Z B,1,5]+Re[Z B,2,5] FIG. 4. Summability problem in φ 4 theory with deenerate minima. Exponentiatin a symmetric double-well with = yields two equal-area aussians, each centered on its respective well. The left well is called perturbative saddle and the riht well is called non-perturbative saddle. For any > 0, no matter how small, Z is the total area under the aussians. Borel summation of perturbation theory yields exactly Z/2 plus an imainary part called the nonperturbative ambiuity, which is exponentially supressed as 0. The Borel sum of the standard diverent perturbation theory of Z accounts only for contributions comin from the perturbative saddle. Such an approach can then be described as the Borel- Hypereometric-Écalle summation method. Meijer-G approximants are extremely-well suited to sum this expansion exactly since as discussed in in Section II they are reularizin analytic continuations of the diverent hypereometric series n+2 F n. Therefore the exact Borel sum of each of the two diverent hypereometric functions that enter Eq.39 can be obtained directly from Meijer-G approximants of order N 5. Definin Z 1 ±i e 1 2F 0 1/4, 3/4, 32, 40 Z F 0 1/4, 3/4, 32, 41 we can find their respective N th order Meijer-G approximations to their Borel sums, denoted Z B,1,N and Z B,2,N b Im[Z] Im[Z B,1,3]+Im[Z B,2,3] Im[Z B,1,5]+Im[Z B,2,5] FIG. 5. Transeries summation by Meijer-G approximants for φ 4 theory with deenerate minima. a Exact real part of Z versus third- and fifth-order summed-up transseries. b Imainary part of Z as computed by third order and fifth order Meijer-G approximants; the third order Meijer-G approximants nearly cancel the non-perturbtaive ambiuity compare with Fi. 3b; the fifth order approximants are exact and thus have zero imainary part.

11 11 The third order approximants are Z B,1,3 = ±e 1 i G 3,1 2,3 Z B,2,3 = G 3,1 2,3 0, ,0, , , 43 0, ,0,0 and the third order approximation to the full transseries is Z = Z B,1,3 + Z B,2,3 + O while the the fifth order approximant is convered to the exact value. The fifth order approximantion to the transseries actually reconstructs Z exactly: Z = Z B,1,5 + Z B,2,5. 45 Hence we use Meijer-G functions to sum up the 2F 0 diverent hypereometric series that appear in the transseries iven by Eq.39. In Fiure 5 we compare the exact Z with the results of third- and fifth-order Borel- Hypereometric-Écalle summation. In Fi. 5a we see that third order and fifth order Meijer-G approximants are excellent approximations to the exact ReZ; the factor of two discrepancy has been removed by proper transseries summation. In Fi. 5b we see that the nonperturbative ambiuity has also been cancelled; while for third-order Meijer-G approximants this cancellation is not complete the size of the ambiuity has been reduced by two orders of manitude, as can be seen by comparin Fi. 5b with Fi. 3b. The cancellation is exact for fifth-order Meijer-G approximants that sum the transseries exactly. Table II shows both the nonperturbative ambiuity and the missin factor 2 in the Meijer-G approximants to the Borel sum of Z Z B,3 and Z B,5 for selected values of ranin from intermediate to very stron couplins. When the Borel- Hypereometric-Écalle summation is adopted the ambiuity is removed and the factor 2 is restored. The third order approximant, while not exact, is remarkably accurate. Usin the interal considered by Marucho [53] we have illustrated the application of Meijer-G approximants to transseries summation. Such approach can be thouht of as Borel-Hypereometric-Écalle summation, as opposed to Borel-Padé-Écalle and consists in usin Meijer- G approximants or hypereometric approximants on the Borel plane; or rational approximations for the ratio between consecutive coefficients to sum the diverent series that contribute to the resurent transseries. In the present case we found that the transseries can be summed exactly with fifth order Meijer-G approximants. While not exact, the third order approximation provides excellent approximations to Z and removes most of the non-perturbative ambiuity. C. Self-interactin QFT Ref. 28 considers the followin partition function of a 0+0-dimensional self-interactin QFT Z = 1 2π π/2 π/2 e sinφ2 2 dφ = π e I0, 4 46 where I 0 x is a modified Bessel function of the first kind. Z has the asymptotic expansion Z Usin the techniques put forward in Section II we once aain obtain a converent set of Meijer-G approximantions to the Borel sum of this series. The convered Meijer-G approximant to the Borel sum is of fifth order and reads Z B,5 = 1 4π G 4,1 0,0,0 3,4 1/2, 1/2,0, = i 1 2π e 4 K As in the previous example, increasin the order of the approximant yields the same Meijer-G function. The convered approximants can be compared with the very accurate, but not exact, third-order Meijer-G approximant Z B,3 = 26 Γ ,6/13 2, /53,0, Γ G 3,1 Evaluatin Z and Z B,5 numerically for selected values of we see that but while Re[Z B,5 ] = Re[Z], 50 Im[Z B,5 ] 0, 51 Im[Z] = Thus we find the non-perturbative ambiuity once aain. For example, evaluatin the third- and fifth-order approximants at = 1 we obtain while Z B,3 = 1 ± iɛ = ± i , Z B,5 = 1 ± iɛ = ± i , Z = 1 =

12 12 Approximant = 1 = 10 = 100 Z B, i i i Z B, i i i Z B,1,3 + Z B,2, i i i Z B,1,5 + Z B,2, Exact TABLE II. Z for 0-dimensional φ 4 theory with deenerate minima evaluated for = 1, = 10 and = 100 usin the third and fifth order Meijer-G approximants to the Borel sum, Z B,3 and Z B,5, as well as the third and fifth order Meijer-G approximations to the resurent transseries, Z B,1,3 + Z B,2,3 and Z B,1,5 + Z B,2,5. We also include the exact result obtained by evaluatin Z numerically in Eq. 36. The approximations to the Borel sum fail to reproduce both the real part and the imainary part of the exact result. In contrast the Meijer-G resummation of the resurent transseries for this problem ives a better account of the exact Z: the third-order Meijer-G approximant is in excellent areement with the exact result, nearly cancellin the non-perturbative ambiuity; the fifth order Meijer-G approximant is exact. All quantities are iven to 12 sinificant diits. Repeatin at the very lare value of = 100 we find Z B,3 = 100 ± iɛ = ± i , Z B,5 = 100 ± iɛ = ± i , while Z = 100 = This comparison shows that the third order Meijer-G approximants are excellent approximations to the exact Borel sum, and that fifth order Meijer-G approximants reproduce the exact Borel sum. Finally it is clear that the series is not Borel summable: usin the Borel sum to estimate the value of Z results in a fictitious imainary part with no clear physical interpretation. The Borel sum fails to reproduce the analytic structure of Z on the complexified -plane. Cherman et al. [28] provide a very clear explanation of the appearence of the non-perturbative ambiuity in this problem, which echoes the discussion iven above for the case of φ 4 theory with deenerate vacua. In the present case, we have also an action sinφ2 2 with a perturbative saddle at φ = 0, with action S 1 = 0, as well as a non-perturbative saddle at φ = π/2 where the action evaluates to 1/4 S 2 = 1/4. In order to properly evaluate Z we thus need to sum a transseries expansion of the form Z, σ 0, σ 1 σ 0 e S0 n=0 c n,0 n + σ 1 e S1 n=0 c n,1 n, 55 where c n,0 c n,1 is the n-th order expansion coefficient around the perturbative non-perturbative saddle, and σ 0 and σ 1 are the correspondin transseries parameters. The transseries parameters have been found to be iven by σ 0 = 1 and σ 1 = i, for Im > 0 and σ 1 = i, for Im < 0. The diverent series in the transseries expansion can both be summed exactly by Borel-Hypereometric summation; the Borel-transformed series are both hypereometric series and n=0 n=0 c n,0 1 n! τ n 2 F 1 2, 1, 1; 2τ, 56 2 c n,1 1 n! τ n 2 F 1 2, 1, 1; 2τ, 57 2 and thus belon to the space of functions that can be reconstructed by our technique. In particular, third order Meijer-G approximants cannot be exact, as the hypereometric approximants in the Borel plane are of the form 2F 1 1, h 1, h 2, h 3 τ rather than 2 F 1 1/2, h 1, h 2, h 3 τ. But fifth order Meijer-G approximants can be exact since 1 2F 1 2, 1 1 2, 1; 2τ = 3 F 2 2, 1 2, 1, 1, 1; 2τ, 58 and therefore the Borel transformed series can be summed up exactly with fifth order Meijer-G approximants. In Fiures 6a-b we compare Meijer-G or Borel-Hypereometric approximants with the exact Z. In Fi 6a we show the real part of Z: the third order approximant already ives an excellent approximation to Re[Z] while the fifth order approximant reproduces its exact value. In Fi 6b we see that both approximants have a non-perturbative ambiuity, which is aain well approximated by the third order approximant. The non-perturbative ambiuity can be removed by the transseries summation by Meijer-G approximants, as shown in Fiure 6c-d. The real part of Z, shown in Fi 6c, is unaffected but the non-perturbative ambiuity has been effectively removed. Comparin Fi. 6b with Fi. 6d we see that thirdorder transseries summation dramatically reduced the size of the non-perturbative ambiuity. For testin and pedaoical purposes we have chosen as a first application of our method the problem of approximatin partition functions in 0-dimensional QFT. For