arxiv:hep-ph/ v1 18 Nov 1996

Similar documents
Asymptotic Expansions of Feynman Integrals on the Mass Shell in Momenta and Masses

INP MSU 96-34/441 hep-ph/ October 1996 Some techniques for calculating two-loop diagrams 1 Andrei I. Davydychev Institute for Nuclear Physics,

arxiv:hep-ph/ v1 21 Jan 1998

arxiv: v1 [hep-ph] 30 Dec 2015

A magic connection between massive and massless diagrams. Department of Physics, University of Bergen, Allegaten 55, N-5007 Bergen, Norway

Reduction of Feynman integrals to master integrals

PoS(DIS2017)295. Hadronic Higgs boson decay at order α 4 s and α 5 s

Single mass scale diagrams: construction of a basis for the ε-expansion.

Large momentum expansion of two-loop self-energy diagrams with arbitrary masses

arxiv:hep-ph/ v1 30 Oct 2002

Evaluating multiloop Feynman integrals by Mellin-Barnes representation

Evaluating double and triple (?) boxes

Two loop O N f s 2 corrections to the decay width of the Higgs boson to two massive fermions

Tadpole Summation by Dyson-Schwinger Equations

Loop integrals, integration-by-parts and MATAD

RunDec: a Mathematica package for running and decoupling of the strong coupling and quark masses

Forcer: a FORM program for 4-loop massless propagators

arxiv: v2 [hep-ph] 20 Jul 2014

arxiv:hep-ph/ v1 8 Oct 1999

Electroweak accuracy in V-pair production at the LHC

Schematic Project of PhD Thesis: Two-Loop QCD Corrections with the Differential Equations Method

Colour octet potential to three loops

arxiv:hep-ph/ v1 25 Jun 1999

The 4-loop quark mass anomalous dimension and the invariant quark mass.

arxiv:hep-ph/ v3 16 Jan 2008

PUZZLING b QUARK DECAYS: HOW TO ACCOUNT FOR THE CHARM MASS

Two-loop self-energy master integrals on shell

Reduction to Master Integrals. V.A. Smirnov Atrani, September 30 October 05, 2013 p.1

Numerical Evaluation of Multi-loop Integrals

K. Melnikov 1. Institut fur Physik, THEP, Johannes Gutenberg Universitat, Abstract

arxiv:hep-lat/ v1 30 Sep 2005

Calculating four-loop massless propagators with Forcer

A numerical test of differential equations for one- and two-loop sunrise diagrams using configuration space techniques

arxiv:hep-ph/ v2 2 Feb 1998

Five-loop renormalization group functions of O(n)-symmetric φ 4 -theory and ǫ-expansions of critical exponents up to ǫ 5

Triangle diagrams in the Standard Model

1 Running and matching of the QED coupling constant

NNLO antenna subtraction with one hadronic initial state

Computing of Charged Current DIS at three loops

gg! hh in the high energy limit

Ginzburg-Landau Theory of Type II Superconductors

arxiv:hep-ph/ v1 1 Jun 1995

Two-loop Heavy Fermion Corrections to Bhabha Scattering

arxiv: v2 [hep-ph] 4 Jun 2018

Scalar One-Loop Integrals using the Negative-Dimension Approach

Precision Calculations to Top- and Bottom-Yukawa Couplings within the SM and BSM

arxiv: v2 [hep-ph] 2 Apr 2019

Corrections of Order β 3 0α 3 s to the Energy Levels and Wave Function

arxiv:hep-ph/ v3 2 Mar 1999

TVID: Three-loop Vacuum Integrals from Dispersion relations

Numerical Evaluation of Multi-loop Integrals

Numerical Evaluation of Loop Integrals

CP-violating Loop Effects in the Higgs Sector of the MSSM

LIMIT ON MASS DIFFERENCES IN THE WEINBERG MODEL. M. VELTMAN Institute for Theoretical Physics, University of Utrecht, Netherlands

Z+jet production at the LHC: Electroweak radiative corrections

Higgs boson production at the LHC: NNLO partonic cross sections through order ǫ and convolutions with splitting functions to N 3 LO

Two-loop Remainder Functions in N = 4 SYM

arxiv: v1 [hep-ph] 22 Sep 2016

Production of e + e pairs to all orders in Zα for collisions of high-energy muons with heavy nuclei

From Tensor Integral to IBP

Search for Next Generations of Quarks and Leptons at the Tevatron and LHC

Multiloop QCD: 33 years of continuous development

REVIEW REVIEW. Quantum Field Theory II

Quantum Field Theory II

Feynman Integral Calculus

arxiv:hep-ph/ v1 2 Nov 1993

Five-loop massive tadpoles

Review of scalar field theory. Srednicki 5, 9, 10

Precision Measurements in Electron-Positron Annihilation: Theory and Experiment

Systems of differential equations for Feynman Integrals; Schouten identities and canonical bases.

arxiv:hep-ph/ v2 12 Jan 2001

arxiv:hep-ph/ v2 22 Feb 1996

Anomalous dimensions and splitting functions beyond the next-to-next-to-leading order

FIRE4, LiteRed and accompanying tools to solve integration by parts relations

HELAC: a package to compute electroweak helicity amplitudes

REDUCTION OF FEYNMAN GRAPH AMPLITUDES TO A MINIMAL SET OF BASIC INTEGRALS

Physics 217 FINAL EXAM SOLUTIONS Fall u(p,λ) by any method of your choosing.

Two-loop massive fermionic operator matrix elements and intial state QED corrections to e + e γ /Z

Quantum Field Theory II

arxiv: v2 [hep-th] 7 Jul 2016

One-Mass Two-Loop Master Integrals for Mixed

arxiv:hep-ph/ v1 21 Dec 1992

Nonequilibrium dynamics and transport near the chiral phase transition of a quark-meson model

Production of energetic dileptons with small invariant masses. from the quark-gluon plasma. Abstract

Maxwell s equations. electric field charge density. current density

Automating Feynman-diagrammatic calculations

Vacuum Energy and Effective Potentials

Heavy quark masses from Loop Calculations

I. INTRODUCTION In a recent paper [1] we managed to derive closed forms for ladder corrections to selfenergy graphs (rainbows) and vertices, in the co

Numerical multi-loop calculations: tools and applications

XV Mexican Workshop on Particles and Fields

OKHEP The Bjorken Sum Rule in the Analytic Approach to Perturbative QCD

arxiv:hep-th/ v1 2 Jul 1998

arxiv:hep-ph/ v1 25 Sep 2002

arxiv:hep-ph/ v2 2 May 1997

NNLO antenna subtraction with two hadronic initial states

Regularization Physics 230A, Spring 2007, Hitoshi Murayama

arxiv:hep-ph/ v1 29 Aug 2006

Bare Higgs mass and potential at ultraviolet cutoff

MHV Diagrams and Tree Amplitudes of Gluons

Transcription:

TTP96-55 1 MPI/PhT/96-122 hep-ph/9611354 November 1996 arxiv:hep-ph/9611354v1 18 Nov 1996 AUTOMATIC COMPUTATION OF THREE-LOOP TWO-POINT FUNCTIONS IN LARGE MOMENTUM EXPANSION K.G. Chetyrkin a,b, R. Harlander c, J.H. Kühn c and M. Steinhauser a a Max-Planck-Institut für Physik, Werner-Heisenberg-Institut, 80805 Munich, Germany b Institute for Nuclear Research, Russian Academy of Sciences, Moscow 117312, Russia c Institut für Theoretische Teilchenphysik, Universität Karlsruhe, 76128 Karlsruhe, Germany Abstract We discuss the calculation of two-point three-loop functions with an arbitrary number of massive propagators and one large external momentum. The relevant subdiagrams are generated automatically. The resulting massless two-point integrals and massive tadpoles are transformed on-line to FORM-expressions ready to be used by existing FORM packages which calculate them analytically. As an example we compute the quartic mass corrections to the photon polarization function. 1 Introduction. In general it is not possible to compute analytically a threeloop propagator diagram with non-vanishing external momentum and internal particles with different masses. However, often a hierarchy between the mass scales exists. A typical example is the case when the magnitude of the external momentum, q, is much larger than all the masses appearing in the diagram. For this kinematical situation there exists a systematic procedure, the Large 1 The complete paper, including figures, is also available via anonymous ftp at ftp://ttpux2.physik.uni-karlsruhe.de/, or via www at http://www-ttp.physik.uni-karlsruhe.de/cgi-bin/preprints/ Extended version of the contribution to the proceedings of the AIHENP96 workshop in Lausanne, September 1996. Preprint submitted to Elsevier Preprint 5 March 2008

Momentum Expansion (LME), which works on the diagrammatic level and delivers an expansion in m 2 /q 2 for a given dimensionally regularized diagram [1]. Possible physical applications of LME at three-loop level include e.g. the production of t t pairs at e + e colliders and the quark mass effects in the decay rates of Z and Higgs bosons. In this work we discuss a computer-algebra implementation of the LME at three-loop level. Note that the one- and two-loop cases are technically much simpler. The number of various prototype diagrams to be considered is rather small and all necessary general diagrammatic manipulations (except for the final calculation of the very diagrams produced by LME) are quite possible to do by hand as it was demonstrated in [2]. In contrast to the two-loop case an application of LME to even a single three-loop diagram leads to dozens of subdiagrams and corresponding Feynman integrals. This makes the use of computer algebra methods a must. Once the expansion is constructed it is to be followed by usual UV renormalization. As the second step is straightforward we shall only deal with bare Feynman integrals regulated by taking a generic space time dimension D = 4 2ǫ. We start with listing the rules of LME. Given a specific diagram with a large external momentum one gets an asymptotic expansion for the momentum going to infinity by applying the following prescription: (1) Generate all subdiagrams of the initial graph such that (a) they contain both vertices through which the large momentum goes into and out from the initial graph (b) they become one-particle-irreducible when the two vertices are identified. (2) Taylor-expand the integrands of these subdiagrams in all small masses and external momenta generated by removing lines from the initial diagram. (3) In the initial diagram, shrink the subdiagram to a point and insert the result obtained from the expansion in (2). (4) Sum over all terms. We will call the subgraphs from step (1) hard subgraphs or simply subgraphs, while the reduced graph, obtained from step (3), will be called the co-subgraph. For the non-planar three-loop topology all hard subgraphs are shown in Fig. 1. The corresponding co-subgraphs are obtained from the full graph by shrinking the displayed lines to a point. The number of terms generated by LME increases rapidly with the number of loops in the initial diagram. Also, the relation between the expansion momenta in the subgraph and the loop momenta of the co-subgraph becomes non-trivial. Nevertheless, LME provides well-defined rules and therefore it should be possible to automate it. 2

Fig. 1. Possible topologies for the subgraphs of an NO-type diagram. 2 Implementation. We restrict ourselves to three-loop two-point functions containing an arbitrary combination of massive and massless lines. In this section the emphasis will be put on the generation of the terms according to the prescription of Sect. 1. Once they are available, existing program packages written in FORM [3] and based on the method of integration-by-parts [4,5] will be used. The realization of the LME by a computer can be divided into two steps: (a) generation of the relevant subgraphs and the corresponding co-subgraphs including the determination of their topologies, and (b) distribution of momenta in the subgraph and co-subgraph respecting the relations between them. As for step (a), one only considers topologies, disregarding any properties of lines except their relative positions. Especially one neglects their momentum, mass, particle type, etc. For the representation of a diagram we label its vertices by integers and specify lines by their endpoints. In that way, a topology is described by the collection of its lines. To generate the subdiagrams we first note that the full graph is always one of the hard subgraphs. The remaining ones are obtained by going through the following steps: Remove any combination of lines from the initial diagram. Remove the emerging isolated dots and binary vertices. Relabel the remaining vertices by {1,..., #vertices} and build all permutations. Compare with a table listing the basic one- and two-loop topologies. If a subgraph passes the last step, the topology of the corresponding cosubgraph is determined in a similar way. The result of this procedure is a 3

database containing all relevant sub- and co-subgraphs, including information about their topology. Note that no selection criteria bound to line properties have been applied so far. Thus, up to this point the procedure is universal and for each topology this part of the program has to be run only once and for all. Coming to step (b), we now use the database resulting from (a) and specialize to a specific diagram, including masses and momenta. We bind the direction of the momentum carried by a line to the order of the labels representing this line. Furthermore, for each subgraph the small external momenta are guided through it in a way that they touch as few lines as possible. The final result is the full input for the FORM-packages, together with all the necessary administrative files like makefiles etc. 3 Application. We present the result for the photon polarization function in QED up to three loops and mass corrections up to O((m 2 /q 2 ) 2 ). The oneand two-loop results are known since a long time in analytical form [6]. A comparison of these one- and two-loop terms up to O((m 2 /q 2 ) 6 ) calculated by our program package with the expansion of the exact results was successful. At three loops the result reads: Π(q 2 ) = 1 { 20 16π 2 9 4 ( ) ln 3 m + m2 m 2 2 ( ) 2 q 8 + 4 + 8 ln 2 q 2 m 2 + α [ ( ) 5 4ζ(3) ln π 6 m + m2 12 ln 2 q 2 m 2 ( ) m 2 2 ( 2 + + 16ζ(3) 10 ln q 2 3 m 12 q 2 ) ] 2 ln2 m 2 ( ) [ α 2 + 2317 (8 π 432 + ζ(2) ln2 31 ) + 199 ζ(3) + 10ζ(5) 3 144 ( 91 + 24 8 ) 3 ζ(3) ln m 1 ( 2 ln2 3 m + m2 25 + ζ(2) (48 ln2 2 q 2 6 46) + 55 70 95 ζ(3) ζ(5) + ln 3 3 6 m + 5 q 2 ) 2 ln2 m 2 ( ) m 2 2 (923 + q 2 54 + 992 9 ζ(3) 48ζ(4) 20ζ(5) + 8B 4 ( + 745 ) 64 + ζ(2) (96 ln2 92) 36 3 ζ(3) ln m 2 + 157 q 2 6 ln2 m + 28 q 2 ) ]} 2 3 ln3... (1) m 2 where B 4 1.762800 [5]. The imaginary part of Eq. (1) is in agreement with [7] which is an important check for our programs. Terms up to O((m 2 /q 2 ) 6 ) have been calculated by our program package and successfully compared with 4

the results obtained by different methods [8]. To conclude, the manual effort for calculating a three-loop propagator-diagram in the Large Momentum Expansion has been reduced to its minimum. This will find several applications in physics, e.g. the O(α 2 s(m 2 /q 2 ) n )-terms (n = 3, 4,...) of the photon polarization function. The use of these terms as input for the Padé method to get the polarization function in the whole energy range will certainly contribute to a further improvement and stabilization of the results of [8]. From the technical point of view, the program should be extendible also to the hard mass procedure, which will be one of the future projects. References [1] For a recent review see V.A. Smirnov, Mod. Phys. Lett. A 10 (1995) 1485. [2] A.I. Davydychev, V.A. Smirnov and J.B. Tausk, Nucl. Phys. B 410 (1993) 325. [3] J.A.M. Vermaseren, Symbolic Manipulation with FORM (Computer Algebra Netherlands, Amsterdam, 1991). [4] K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. B 192 (1981) 159. [5] D.J. Broadhurst, Z. Phys. C 54 (1992) 54. [6] G. Källén and A. Sabry, K. Dan. Vidensk. Selsk. Mat.-Fys. Medd. 29 (1955) No. 17, see also J. Schwinger, Particles, Sources and Fields, Vol.II, (Addison-Wesley, New York, 1973). [7] K.G. Chetyrkin and J.H. Kühn, Nucl. Phys. B 432 (1994) 337. [8] K.G. Chetyrkin, J.H. Kühn and M. Steinhauser, Phys. Lett. B 371 (1996) 93; K.G. Chetyrkin, J.H. Kühn and M. Steinhauser, MPI/PhT/96-27, TTP-96-13, hep-ph/9606230, to be published in Nucl. Phys. B. 5