D i k B(n;D)= X LX ;d= l= A (il) d n d d+ (A ) (kl)( + )B(n;D); where c B(:::;n c;:::) = B(:::;n c ;:::), in prticulr c n = n c. Using the reltions [n

Similar documents
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

5.7 Improper Integrals

Definite integral. Mathematics FRDIS MENDELU

Lecture 14: Quadrature

Definite integral. Mathematics FRDIS MENDELU. Simona Fišnarová (Mendel University) Definite integral MENDELU 1 / 30

Lecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

MA 124 January 18, Derivatives are. Integrals are.

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

Introduction to Determinants. Remarks. Remarks. The determinant applies in the case of square matrices

14.4. Lengths of curves and surfaces of revolution. Introduction. Prerequisites. Learning Outcomes

Math 8 Winter 2015 Applications of Integration

Review of Calculus, cont d

Abstract inner product spaces

4.4 Areas, Integrals and Antiderivatives

8 Laplace s Method and Local Limit Theorems

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as

MAC-solutions of the nonexistent solutions of mathematical physics

Problem Set 3 Solutions

Improper Integrals, and Differential Equations

The Regulated and Riemann Integrals

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

Construction of Gauss Quadrature Rules

Overview of Calculus I

Math 1B, lecture 4: Error bounds for numerical methods

Summary: Method of Separation of Variables

Math& 152 Section Integration by Parts

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

We know that if f is a continuous nonnegative function on the interval [a, b], then b

Chapter 1. Basic Concepts

The graphs of Rational Functions

MAA 4212 Improper Integrals

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

Week 10: Line Integrals

Section 4.8. D v(t j 1 ) t. (4.8.1) j=1

Lecture 1. Functional series. Pointwise and uniform convergence.

Consequently, the temperature must be the same at each point in the cross section at x. Let:

p(t) dt + i 1 re it ireit dt =

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Riemann Sums and Riemann Integrals

Math 113 Exam 2 Practice

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Lecture 19: Continuous Least Squares Approximation

The Fundamental Theorem of Calculus. The Total Change Theorem and the Area Under a Curve.

Review of basic calculus

Orthogonal Polynomials and Least-Squares Approximations to Functions

Best Approximation in the 2-norm

Chapter 3 MATRIX. In this chapter: 3.1 MATRIX NOTATION AND TERMINOLOGY

Riemann Sums and Riemann Integrals

Math 5440 Problem Set 3 Solutions

Lesson 1: Quadratic Equations

MATH , Calculus 2, Fall 2018

. Double-angle formulas. Your answer should involve trig functions of θ, and not of 2θ. sin 2 (θ) =

Best Approximation. Chapter The General Case

4181H Problem Set 11 Selected Solutions. Chapter 19. n(log x) n 1 1 x x dx,

Numerical Integration

Math 5440 Problem Set 3 Solutions

Math 554 Integration

Stuff You Need to Know From Calculus

Name Solutions to Test 3 November 8, 2017

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER /2019

4. Calculus of Variations

5.2 Volumes: Disks and Washers

Recitation 3: More Applications of the Derivative

potentials A z, F z TE z Modes We use the e j z z =0 we can simply say that the x dependence of E y (1)

Chapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

Field-Induced Axion Luminosity of Photon Gas via a-interaction N.V. Mikheev, A.Ya. Parkhomenko and L.A. Vassilevskaya Yaroslavl State (Demidov) Univer

Math 31S. Rumbos Fall Solutions to Assignment #16

New Expansion and Infinite Series

221B Lecture Notes WKB Method

Numerical Analysis: Trapezoidal and Simpson s Rule

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

Math 3B: Lecture 9. Noah White. October 18, 2017

Heat flux and total heat

Section 7.1 Integration by Substitution

A5682: Introduction to Cosmology Course Notes. 4. Cosmic Dynamics: The Friedmann Equation. = GM s

4 The dynamical FRW universe

Jackson 2.7 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:

ODE: Existence and Uniqueness of a Solution

a a a a a a a a a a a a a a a a a a a a a a a a In this section, we introduce a general formula for computing determinants.

Numerical Analysis. 10th ed. R L Burden, J D Faires, and A M Burden

Introduction to the Calculus of Variations

AMATH 731: Applied Functional Analysis Fall Additional notes on Fréchet derivatives

1.9 C 2 inner variations

Module 6: LINEAR TRANSFORMATIONS

Calculus II: Integrations and Series

Lecture 20: Numerical Integration III

Math 116 Calculus II

ENGI 3424 Engineering Mathematics Five Tutorial Examples of Partial Fractions

Practice final exam solutions

How can we approximate the area of a region in the plane? What is an interpretation of the area under the graph of a velocity function?

Theoretical foundations of Gaussian quadrature

Reversals of Signal-Posterior Monotonicity for Any Bounded Prior

Transcription:

Explicit solutions of the multi{loop integrl recurrence reltions nd its ppliction? P. A. BAKOV ; nstitute of ucler Physics, Moscow Stte University, Moscow 9899, Russi The pproch to the constructing explicit solutions of the recurrence reltions for multi{loop integrls re suggested. The resulting formuls demonstrte high eciency, t lest for 3{loop vcuum integrls cse. They lso produce new type of recurrence reltions over the spce{time dimension. Vcuum cse Recently [,] new pproch to implement recurrence reltions [3] for the Feynmn integrls ws proposed. n this work we extend the generl formuls for the solutions of the recurrence reltions to the multi{loop cse. Let us consider rst vcuum L-loop integrls with = L(L + )= denomintors (so tht one cn express through them ny sclr product of loop moment) of rbitrry degrees: B(n;D)=m n i D = X ij LD d D p :::d D p L D n :::D n A (ij) p i p j m ; p k p l = X = ; () (A ) (kl)(d + m ): () The recurrence reltions tht result from integrtion by prts, by letting (@=@p i ) p k ct on the integrnd [3], re:? Tlk presented t the AHEP'96 (EPFL-UL Lusnne, Sept. -6 996) Supported in prt by the RBRF (grnt 96{0{00654), TAS (grnt 93-80- ext); e-mil: bikov@theory.npi.msu.su Preprint submitted to Elsevier Preprint 9 ovember 996

D i k B(n;D)= X LX ;d= l= A (il) d n d d+ (A ) (kl)( + )B(n;D); where c B(:::;n c;:::) = B(:::;n c ;:::), in prticulr c n = n c. Using the reltions [n d d+ ; ]= d ; X = A (il) (A ) (kj) = (i (k l) j) ; they cn be represented s D L i k B(n;D)= X LX ;d= l= (A ) (kl)( + )A (il) d n d d+ B(n;D): (3) The common wy of using these reltions is step{by{step reexpression of the integrl () with some vlues of n i through set of integrls with shifted vlues of n i, with the nl gol to reduce this set to liner combintion of severl "mster" integrls k (D) with some "coecient functions" F k (n;d):. B(n;D)= X k F k (n;d) k (D) evertheless, to nd proper combintions of these reltions nd proper sequence of its use is the mtter of rt even for the tree{loop integrls with one mss [4]. Then, even in cses when such procedures were constructed, they led to very time nd memory consuming clcultion becuse of lrge reproduction rte t every recursion step. nsted, let us construct the F k (n;d) directly s solutions of the given recurrence reltions. ote, tht if we nd ny set of the solutions, we could construct F k (n;d) s their liner combintions. Let us try the solution of (3) in the following form: f k (n) = ({) dx dx x n xn g(x ) (4) where integrl symbols denote subsequent complex integrtions with contours which will be described lter. Acting by some opertor O i ( ;n + ) (ll decresing opertors should be plced to the left) on (4) nd performing the integrtion by prts one gets (s.t. re surfce terms):

O i ( ;n + )f k (n)= ({) dx dx x n xn O i (x ;@ )g(x ) + (s. t.): So, if we choose the g(x ) s the solution of O i (x ;@ )g(x ) = 0 nd cncel the surfce terms by proper choosing of integrtion contours (for exmple, closed or ended in the zero points) we nd tht (4) is solution of reltions O i ( ;n + )f k (n) = 0, nd dierent choices of contours correspond to dierent solutions. The dierentil equtions for (3) hve the solution g(x )=P(x + ) (D L )=, where P (x ) = det( X = (A ) (kl) x ) is the polynomil in x of degree L, so we get the desirble solutions of (3): f k (n;d)= ({) dx dx x n xn det((a ) (kl)(x + )) D L : (5) Finlly, let us derive from (3) the recurrence reltions with D-shifts. ote tht if f k (n i ;D) is solution of (3), then by direct substitution to (3) one cn check tht P ( + )f k (n i ;D ) lso is solution. Hence, if f k (n i ;D)is complete set of solutions, then f k (n i ;D)= X n S k n(d)p( + i )f n (n i ;D ); where the coecients of mixing mtrix S is numbers, tht is do not ct on n i. For the solutions (5) the S is the unit mtrix (the incresing of D by leds to ppering of fctor P (x ) in the integrnd of (5)), but the desire to come to some specic set of mster integrls my led to nontrivil mixing. These reltions look dierent from recently proposed in [5], lthough further investigtions cn give some connections with them. To check the eciency of this pproch we evluted (using REDUCE) the rst 5 moments in the smll q expnsion of the 3-loop QED photon vcuum polriztion. The 3-loop contribution to the moments re expressed through bout 0 5 three{loop sclr vcuum integrls with four mssive nd two mssless lines. The integrl (5) in this cse cn be solved to nite sums of the Pochhmmer's symbols (see []). Moreover, it is not necessry to evlute these integrls seprtely. nsted, we evluted few integrls of (5) type, 3

but with P D= producted by long polynomil in x i (the results see in [,6], they re in greement with QCD clcultions [7] mde by FORM). The comprison with the recursive pproch shows resonble progress: the common wy used in [6] demnds severl CPU hours on DEC-Alph to clculte full D dependence of the rst moment, nd further clcultions becme possible only fter trunction in (D= ). n the present pproch the full D clcultion for the rst moment demnds few minutes on PC. on{vcuum cse Suppose tht integrls () depend on R externl moment p i (L<iL+R). The number of the denomentors re now = L(L +)=+LR, nd the number of dditionl ("externl") invrints re = R(R +)=. Let us expnd the integrls in forml seri over "denomintor{like" objects D of () type with = +; ::; +, depending on externl moment only: B(n l;(l=;:::; );p k;(k=l+;:::;l+r) )= d D p :::d D p L D n :::Dn = = X n i (i> ) m n i+ +LD b(n i ; (i=;:::; + ) ) Y+ i= + D n i i : (6) We dene such generl expnsion in order to write the recurrence reltions in compct form, in prctice the coecients A (ij) nd my bevery simple. The expnsion with negtive n i corresponds to the lrge moment expnsion, with positive ones to the expnsion ner points m. The n i cn lso be noninteger, but with unit shifts. Acting by (@=@p i ) p k, (i = ;:::;L;k = ;:::;L + R) on the integrnd we get recurrence reltions. The dditionl reltions we get cting by p k (@=@p i ), (i; k = L +;:::;L+R) on both sides of (6). These new reltions look like the old ones with only exception tht they hve no terms proportionl to spce{time dimention D. The complete set of recurrence reltions is now ((D L R ) i k (D R ) ^ i k) b(n;d)= X+ = ;d= X L+R l= (A ) (kl)( + )A (il) d n d d+ b(n;d); 4

where ^ i k=( i k if i; k > L, else 0). The corresponding dierentil equtions hve the solution g(x )=g 0 (x + ), where g 0 (x) = det (A ) (kl) x D L R det 0 (A ) (kl) x D R ; (7) nd det 0 denotes the minor with k; l > L. So, one cn use the representtion (4), but the problem of resolving it to explicit formuls demnds futher investigtions. Finlly note tht one cn formlly obtin the formuls (5, 7) by "chnge of integrtion vribles" from loop moment to "denomentor{like objects" D. The weight function for this chnge is d D p d D p L Y i (D i =m x i ) / det((a ) (kl)(x + )) D L : References [] P. A. Bikov, Phys. Lett. B385 (996) 404, hep-ph/960367. [] P. A. Bikov, hep-ph/960454. [3] K. G. Chetyrkin, F. V. Tkchov, ucl. Phys. B9 (98) 59; F. V. Tkchov, Phys. Lett. B00 (98) 65. [4] D. J. Brodhurst,. Phys. C54 (99) 599; L. V. Avdeev, Comp. Phys. Commun. 98 (996) 5. [5] O. V. Trsov, Phys. Rev.D54 (996) 6479. [6] P. A. Bikov, D. J. Brodhurst Preprint OUT-40-54, hep-ph/9504398, P- 95-3/377, Proc. of 4th nt. Workshop on Softwre Engineering nd Articil ntelligence for High Energy nd ucler Physics, Pis, tly, 995, eds. B. Denby, D. Perret{Gllix, 67. [7] K. G. Chetyrkin, J. H. Kuhn, M. Steinhuser, TTP-96-3 (996), hep-ph/960630. 5