# Full text of "Resummation of perturbative QCD by Pade approximants"

## See other formats

TAUP-2451-97 Resummation of perturbative QCD by Pade approximant sQ Einan Gardi School of Physics and Astronomy Raymond and Beverly Sackler Faculty of Exact Sciences Tel-Aviv University, 69978 Tel-Aviv, Israel e-mail: gardi@post.tau.ac.il Abstract In this lecture I present some of the new developments concerning the use of Pade Approximants (PA's) for resumming perturbative series in QCD. It is shown that PA's tend to reduce the renormalization scale and scheme dependence as compared to truncated series. In particular it is proven that in the limit where the (3 function is dominated by the 1-loop contribution, there is an exact symmetry that guarantees invariance of diagonal PA's under changing the renormalization scale. In addition it is shown that in the large (3q approximation diagonal PA's can be interpreted as a systematic method for approximating the flow of momentum in Feynman diagrams. This corresponds to a new multiple scale generalization of the Brodsky-Lepage-Mackenzie (BLM) method to higher orders. I illustrate the method with the Bjorken sum rule and the vacuum polarization function. 1 A lecture given at the Cracow School of Theoretical Physics, Zakopane, May 30 - June 10, 1997 I will talk about resummation of perturbative series in QCD The basic question I deal with is how to use finite order perturbative calculations in QCD to make unambiguous theoretical predictions, with controlled errors. As experiments improve one requires the theoretical predictions to be more accurate. However, in QCD it is very hard to get accurate predictions, ba- sically because the coupling constant is large. This leads to non-negligible non-perturbative effects as well as a badly divergent and renormalization scheme dependent perturbative series. In this talk I will show that PA's which start out as an alternative to a finite order perturbative series having the same formal accuracy, actually have an important advantage over the finite order series^. Through the re- summation of certain all-order effects related to the running of the coupling- constant, PA's become independent of the choice of the renormalization scale and therefore lead to more accurate and more reliable predictions. The ma- terial presented in this lecture appears in greater detail in ref. J|, |J. The outline of the talk is as follows: I will start by introducing the PA's method and the problem of renormalization scale dependence in QCD. I will mention some of the other ideas that were raised to confront the problem of renormalization scale dependence and show how PA's solve it in a most elegant way. Then I will address the question of what higher order effects are summed-up by PA. I will show that there is a direct interpretation of PA's in terms of approximating the momentum distribution of virtual gluons in Feynman diagrams. First, what are PA and how do I use them? I start with an effective charge related to some physical observable, written as a power series in x: S n = x (l + T\x + r 2 x 2 + • • • + r n x n ) (1) where x = a s /7r. A PA is constructed by writing a ratio of two polyno- mials such that when expanded back to a Taylor series, it gives the known coefficients r\ through r n : 1 + CL\X + ... + (InX N 1 + bix + ... + b M x ! P x[ n/m] =x - , ^ , 7T - ZM , N + M = n 2 For a recent review on the nature of perturbative series and resummation techniques see ref. 3 PA have various successful application in physics. Examples of applications to statis- tical physics and quantum field theory are listed in ref. [Q. Applications to QCD appear in ref. fll 1 There is a theorem that for any degree N in the numerator and M in the denominator such that N + M = n there is a unique PA function || . I will mainly deal here with diagonal PA's which are written, in my notation, as x[N/N + 1], having one power of x out of the brackets. I will soon come back to discuss PA's. A renormalized perturbative series in QCD is not expected to give exact predictions for measurable quantities due a few limitations. First, the series is divergent and not even Borel-summable. The resulting ambiguity is related to the existence of non-perturbative effects. Second, at any given order, the partial-sum depends on non-physical parameters, such as the renormalization scale. This also makes the prediction ambiguous. Let us concentrate on the renormalization scale dependence. Usually, when we calculate some observable R that depends on one external mo- mentum Q 2 in perturbation theory, we choose as an expansion parameter the renormalized coupling-constant at the external scale Q 2 . This "natu- ral" choice of fi 2 = Q 2 is, however, quite arbitrary. We can, just as well, use some other expansion parameter y = a s (e t Q 2 )/7c, where t ^ 0. The renormalization group equation -J = f3 x 2 + fax 3 + f3 2 x 4 + ■■■ at determines how the two couplants are related, x = y + (3 t y 2 + (fit 2 + p x t) y 3 + (> 3 t 3 + ^ ftt 2 + ftt) y 4 + ■ ■ ■ (2) and thus how the finite order series can be written in terms of y: S n {t)=y (l + r x y + f 2 y 2 + f 3 y 3 ...+ r n y n ) . (3) The new coefficients fj are different from the original coefficients r^, so as to compensate for the scale shift, such that the total effect is some residual dependence on t which is of the next, uncalculated order. Still, in QCD, since the coupling constant is large, the numerical difference due to the change of scale can be quite large. This limits the predictive power of the theory. Beyond two-loops, there is also the question of scheme dependence which can be parameterized by the higher-order coefficients of the (3 function, 03 and onward. 2 In order to test the significance of this scale and scheme dependence, we studied the polarized Bjorken Sum-Rule. In fig. 1 the Bjorken effective charge at NNLO for Q 2 = 20 GeV 2 is plotted as a function of the renormal- ization group non-physical parameters: the coupling x = a s (fi 2 )/ii and the second coefficient of the /3 function: C 2 = fa/Po- We see that the surface is far from being flat. The same surface is drawn again in fig. 2, but here - as a contour plot. The thick lines are contours of equal effective charge. Large renormaliza- tion scheme dependence corresponds to large higher-order corrections, since these are required to compensate for the scale dependence. This observa- tion makes it clear that we should carefully choose the renormalization scale and scheme that we are using. In fig. 2 one can identify a region of rela- tively low renormalization scale and scheme dependence. Specific scales and schemes are chosen according to different criterions such as the method of Effective Charges 0, the Principal of Minimal Sensitivity M and the BLM scale-setting method |§. For the Bjorken Sum- Rule example (fig. 2), all of the above are located in the central region of low renormalization scale and scheme dependence. Note that in this case, MS , with /i 2 = Q 2 is not a good choice. Let's go back to eq. (0) that describes the scale transformation relating the coupling- const ants x and y defined at two different scales. If we assume that the 1-loop coefficient of the (3 function, j3 Q , is large enough, i.e. Po > PiX { for any i > 1, we can approximate the full relation by one that includes only the leading terms in (3q\ x ~y + (3 ty 2 + Pit 2 y 3 + /? V + • • • This can be written in a closed form: = y l-Poty It is important to realize that in the physical case of QCD with 3 to 5 flavors, this approximation is good. Fig. 3 shows the renormalization scale transformation itself, namely the running coupling constant as a function of the scale. The dashed line is the best we know of the running coupling in 3 QCD (it includes the 4- loop effects), and the solid line is the 1-loop running coupling. The two are quite close and I shall use here the 1-loop formula. I now get to the main point, which is the independence of PA on the renormalization scale. We saw that partial-sums as usually written in per- turbation theory always yield different results in different renormalization scales: S n (0) ± S n (t) where 5„(0) refers to fi 2 = Q 2 as in eq. ([!]), S n (t) refers to fi 2 = e t Q 2 as in eq. (0), and x — yj (1 — Poty)- However, if we construct a diagonal PA from the series in x (eq. ([[])), , s 1 + d\X + ... + a N x N P.WN + 1]{X) - X 1+biX+ ^ +bN+iX N + l and independently, another PA from the series in y (eq. (0)), 5 ( v 1 + aiy + ... + a N y N Py[N/N+i]{y) = y ~ — — ~ — 1 + hy + ... + b N+1 y N+L we will get the same result in both: Px{N/N+l](%) = Py[N/N+l](y)- This is due to the mathematical property of diagonal PA's: they are invariant under homographic transformations of the PA argument (x — > x/(l + Kx), see [|, Q). We know that the all-order result does not depend on the renor- malization scale. The fact that diagonal PA are invariant suggests that they correctly resum certain all-order effects that are related to the running of the coupling-constant. Non-diagonal PA are not exactly invariant. However, on the global level (for large scale shifts t) they always have a reduced scale dependence [|]. Going back to the example we examined above, namely the NNLO Bjorken sum-rule, we show in fig. 4 the a; [0/2] PA. Clearly (compare with the partial- sum of fig. 1) the renormalization scale and scheme dependence is almost completely eliminated! Non-diagonal PA's may be dangerous, since specific renormalization scales and schemes are sometimes particularly deviant, as in the example of the 4 x[l/l] PA for the Bjorken sum-rule shown in fig. 5. Therefore it is best to use a diagonal x[N — 1/N] PA. I now consider the question of what higher-order contributions are summed- up by diagonal PA's. It turns out that we can get some rigorous results [|J if we limit ourselves to the "large /?o" approximation (T^, IT[ |Tj|, where only the leading term in /3q in each perturbative coefficient is taken into account. This approximation corresponds to summing certain higher-order contributions that are related to the exchange of one virtual gluon. I use Neubert's formulation [10|, where resummation is achieved by using the run- ning coupling-constant at the vertices. The resummation integral is then a weighted average of the coupling-constant at all scales: w(k)a s (k 2 )d 4 k = / p(s)x v (e s Q 2 )ds J —oo where w(k) is the Feynman integrand and s = In {j^s)- The superscript V stands for the V-scheme which is the most convenient renormalization scheme for my purposes. While a specific scheme is used here in order to simplify the formulae, it is important to understand that the above resummation integral is scheme-invariant [10|. The function p(s) describes the distribution of momentum of the exchanged gloun, and x v (e s Q 2 ) describes the interaction strength as a function of the momentum. Using a 1-loop formula for x v (e s Q 2 ) I get: Clearly the integral includes contributions from an infinite set of diagrams. The exact distribution function (in the large /3q approximation) has been calculated for a few observables, such as the vacuum-polarization D-function [ |13| , I0| which I shall use here as an example. A representative diagram is the following: 5 In this particular diagram the exchanged gluon is dressed by fermion loops. However, gluonic corrections that are related to the 1-loop running of the coupling are resummed in the above technique just as well. Of course, the resummation integral is not well defined, due to the inte- gration over the infra red pole in the 1-loop formula for the running coupling (Landau pole). This is how infra red renormalons appear in this formulation. I will not deal here with the renormalon ambiguity which cannot be settled completely within perturbative QCD, but rather use the resummation inte- gral to study the PA method. The general methodology is to assume that only a first few coefficients in the perturbative series are known, construct a PA basing on this limited information and then compare the PA with the exact all-order result. We will see below that this comparison can be done also on the level of the momentum distribution function, since PA's can be interpreted as what one obtains by replacing the continuous momentum dis- tribution function with a particular discrete distribution. It was found empirically that the momentum distribution function in the large /3q approximation is a non-negative function in many physical examples [ PH , §]. This justifies a posteriori the probabilistic interpretation implied by the name 'momentum distribution'. If indeed p(s) > for any s, then the resummation integral defines a so-called Hamburger function ||: /OO y /*00 y p(s)-—ds = / — - d<f>(s). (4) -oo 1 + SZ J —oo 1 + SZ where z = /3oX V (Q 2 ) and </>(s) is the indefinite integral of p(s). The pertur- 6 bative coefficients are moments of the distribution function: /oo -oo for i > 0. There is a theorem || that guarantees that for a Hamburger function, a z[N — l/N] PA constructed from the partial-sum: 2N-1 * E fit-*)* 8=0 can be written as: f(z) ~z[N-l/N] =ETT^ i=i 1 "I" </i^ with g« rea/ and r; > for i = 1, 2, ...iV. Through this decomposition of the PA function, together with eq. (f|) one realizes that the PA corresponds to approximating the all-order continuous distribution function by a sum of iV weighted (^-functions: N Pn(s) = E r ^( s ~ i=l or, equivalently, its indefinite integral <fi(s), by a piece- wise constant function composed of N steps: N <I>n(s) = E r ^( s ~ Qi)- i=l Note that Pn(s) is optimal (and unique, of course) since the equation for constructing the PA imply that Pn{s) reproduces the first 2N moments of the distribution function, which we know. Using a diagonal PA of a Hamburger function to identify the optimal scales (ft) and weights (r*) is the basic idea behind the method of Gaussian quadrature for numerical integration . In the Brodsky-Lepage-Mackenzie (BLM) approach || one evaluates the coupling-constant at a scale that corresponds to the average momentum of the exchanged gluon (the BLM scale). This is exactly equivalent to approxi- mating the distribution function by a single 5-function located at its center. 7 This same effect can also be achieved simply by using a a; [0/1] PA for the leading (3q series [§J. In the method described above one uses an x[N — 1/N] PA of the leading (3q series that corresponds to approximating the momen- tum distribution function by a set of iV weighted 5 functions. Therefore it can be viewed as a generalization of the BLM method for higher-orders. In the following, I illustrate the above ideas with the vacuum-polarization D-function. The all-order distribution function p(s) [TIJ is plotted in fig. 6 as a continuous line. We see that there are contributions from both UV scales (positive s) and IR scales (negative s). The different symbols correspond to the locations (qi) and weights (r») of the diagonal PA's poles. For the x[0/l] PA, there is one ^-function at the BLM scale. For the a; [1/2] PA, there are two 5- functions, and so on. In Fig. 7 we see how using x[N — /N] PA corresponds to approximating the integral distribution function 0(s) by a piecewise constant function, composed of N steps. To conclude, we saw that diagonal PA's can be used to resum certain all-order effects that are related to the running of the coupling constant, and thus provide a systematic method for obtaining reliable scale-invariant predictions. There is a rigorous relation between diagonal PA and the mo- mentum distribution of virtual gluons. I stress that this result holds only for a single gluon exchange, i.e. within the large (3q approximation. The way to go beyond this approximation is still unclear. Nevertheless, from the results presented here for the Bjorken Sum-Rule, it is clear that PA's are an important tool for QCD phenomenology. Acknowledgments I thank Marek Karliner, Stan Brodsky, John Ellis and Mark Samuel for a fruitful cooperation and very useful discussions. The research was supported in part by the Israel Science Foundation administered by the Israel Academy of Sciences and Humanities, and by a Grant from the G.I.F., the German- Israeli Foundation for Scientific Research and Development. References [1] J. Fischer, On the Role of Power Expansions in Quantum Field Theory, |hep-ph/9704351 . 8 [2] M.A. Samuel, G. Li and E. Steinfelds, Phys. Rev. D48(1993)869 and Phys. Lett. B323(1994)188; M.A. Samuel and G. Li, Int. J. Th. Phys. 33(1994)1461 and Phys. Lett. B331(1994)114. [3] M.A. Samuel, J. Ellis and M. Karliner, Phys. Rev. Lett. 74(1995)4380; J. Ellis, E. Gardi, M. Karliner and M.A. Samuel, Phys. Lett. B366(1996)268 and Phys. Rev. D54(1996)6986; I. Jack, D.R.T. Jones and M.A. Samuel, Asymptotic Pade Approximants and the SQCD /3 Function, |hep-ph/ 9706249] [4] E. Gardi, Phys. Rev. D56(1997)68. [5] S.J. Brodsky, J. Ellis, E. Gardi, M. Karliner and M.A. Samuel, Pade Approximants, Optimal Renormalization Scales, and Momentum Flow in Feynman Diagrams, hep-ph/9706467 , CERN-TH-97/126. Accepted for publication in Phys. Rev. D. [6] George A. Baker, Jr. and Peter Graves-Morris, Gian-Carlo Rota En- cyclopedia of Mathematics and its Applications (Addison- Wesley, 1981), Vol. 13 and 14. [7] G. Grunberg, Phys. Rev. D29(1984)2315. [8] P.M. Stevenson, Phys. Rev. D23(1981)2916. [9] S.J. Brodsky, G.P. Lepage and P.M. Mackenzie, Phys. Rev. D28(1983)228; [10] M. Neubert, Phys. Rev. D51(1995)5924. [11] M. Beneke and V.M. Braun, Phys. Lett. B348(1995)513. [12] C.N. Lovett-Turner and C.J. Maxwell, Nucl. Phys. B432(1994)147 and Nucl. Phys. B452(1995)188. [13] D.J. Broadhurst, Z. Phys C58(1993)339. [14] M. Beneke, Phys. Lett. B307(1993)154; Nucl. Phys. B405(1993)424. 9 12 16