## Fiesta.dvi

**Applications of FIESTA**
*Institut für Theoretische Teilchenphysik, Karlsruhe Institute of Technology (KIT), D-76128Karlsruhe, GermanyE-mail: *
**A.V. Smirnov**

*Scientific Research Computing Center, Moscow State University, 119992 Moscow, Russia*

E-mail:
Sector decomposition in its practical aspect is a constructive method used to evaluate Feynman
integrals numerically. We present a new program performing the sector decomposition and in-
tegrating the expression afterwards. The program can be also used in order to expand Feynman
integrals automatically in limits of momenta and masses with the use of sector decompositions
and Mellin–Barnes representations. The program is parallelizable on modern multicore com-
puters and even on multiple computers. Also we demonstrate some new numerical results for
four-loop massless propagator master integrals.

*13th International Workshop on Advanced Computing and Analysis Techniques in Physics ResearchFebruary 22-27, 2010Jaipur, India*
Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence.

*Applications of FIESTA*
Originally sector decomposition was used as a tool for analyzing the convergence and prov-
ing theorems on renormalization and asymptotic expansions of Feynman integrals [, .

After [], the sector decomposition approach has become an efficient tool for numerical evaluating
Feynman integrals (see Ref. for a recent review). At present, there are two public codes per-
forming the sector decomposition [and . The latter one was developed by the present authors;
it is named FIESTA which stands for "Feynman Integral Evaluation by a Sector decomposiTion
Approach". Recently FIESTA was greatly improved in various aspects [].

During the last year FIESTA was widely used, some of application are listed in [. In ]
we used FIESTA in order to confirm numerically the recent analytic results for

*master integrals*
(MI's) for four-loop massless propagators which recently were analytically evaluated in . Here
we provide some more numerical results for extra orders in epsilon expansions for these MI's.

**2. Theoretical background and software structure**
FIESTA calculates Feynman integrals with the sector decomposition approach. After per-
forming Dirac and Lorentz algebra one is left with a scalar dimensionally regularized Feynman
integral []

*F*(

*a*1, . . ,

*an*) R · · · R d

*dk*1.d

*dkl*
, where

*d *= 4 − 2ε is the space-time dimension,

*a*
*E *1 .

*Ean*
indices,

*l *is the number of loops and 1/

*En *are propagators. We work in Minkowski space wherethe standard propagators are the form 1/(

*m*2 −

*p*2 −

*i*0). Other propagators are permitted, see .

Substituting 1
dαα

*ai*−1

*e*−

*iEi *, after usual tricks , performing the decomposition of
the integration region into the so-called

*primary sectors *and making a variable replacement,
one results in a linear combination of integrals R 1

*i *. .

*dxn*′

*j*=1

*x j*
If the functions

*UA*−(

*l*+1)

*d*/2 had no singularities in ε, one would be able to perform the expansion
in ε and perform the numerical integration afterwards. However, in general one has to resolvethe singularities first. Thus, one starts a process the sector decomposition aiming to end with a
sum of similar expressions, but with new functions

*U *and

*F *which have no singularities (all the
singularities are now due to the part ∏

*nj*=1

*x*′

*a*′

*j*−1). The way sector decomposition is performed is
called a

*sector decomposition strategy *([) and is an essential part of the algorithm (let us
also mention a geometrical approach to sector decomposition [] which is rather complicated in
implementation as a strategy on a computer but promises to be the optimal one).

After the sector decomposition one resolves the singularities by evaluating the first terms of
the Taylor series: in those terms one integration is taken analytically. Afterwards the ε-expansioncan be performed and finally one can do the numerical integration.

FIESTA is written in Mathematica [and C. The user is not supposed to use the C
part directly as it is launched from Mathematica via the Mathlink protocol in order to per-
form a numerical integration. To run FIESTA, the user has to load the FIESTA package into
Mathematica 6 or 7. In order to evaluate a Feynman integral one has to use the command
loop_momenta is a list of all loop momenta, propagators is a list of all propagators, subst
is a list of substitutions for external momenta, masses and other values. For example,

*Applications of FIESTA*
SDEvaluate[UF[{k},{-k2,-(k+p1)2,-(k+p1+p2)2,-(k+p1+p2+p4)2},{p21 →0,p22 →0,p24 →0, p1 p2 →-s/2,p2 p4 →-t/2,p1 p4 →-(s+t)/2,s→-3,t→-1}], {1,1,1,1},0]
evaluats the massless on-shell box diagram with Mandelstam variables equal to −3 and −1.

**3. Numerical results for four-loop massless propagators**
**Figure 1: ***M*61–

*M*43: the thirteen complicated four-loop master integrals according to [. The two MI's

*M*52 and

*M*43 can be identically expressed through the three-loop nonplanar MI

*N*0.

In [a full set of four-loop massless propagator-like MI's was identified. There are 28 inde-
pendent MI's. Analytical results for these integrals were obtained in . The most complicated
MI's are demonstrated on Fig. ε

*m *after

*Mi j *stands for the maximal term in ε-expansion of

*Mi j*which one needs to know for evaluation of the contribution of the integral to the final result for a
four-loop integral after reduction is done, see []. Two of the complicated integrals (

*M*43 and

*M*52)are related by a simple factor with the three-loop MI

*N*0 [so it is enough to evaluate remainingeleven complicated MI's

*M*61–

*M*36 as well as first three terms of the ε-expansion of

*N*0.

We calculated them (for

*q*2 = −1) using FIESTA with the Cuba[Vegas integrator and
1 500 000 sampling points for integration. Our results alongside with the corresponding analytical
expressions (transformed to the numerical form) from [look like follows1:
: 0.08333 ± 0 (0.08333); ε−3: 0.916667 ± 0.000018 (0.91666); ε−2: 5.64251 ± 0.00022
(5.6425109); ε−1: 27.6413 ± 0.00077 (27.6412581); ε0: 98.638 ± 0.0034 (98.637928); ε1:342.736 ± 0.012 (342.7349920); ε2: 857.88 ± 0.048 (857.8735165); ε3: 2659.84 ± 0.19(2659.825402); ε4: 4344.28 ± 0.75 (unknown); ε5: 17483.1 ± 5.7 (unknown).

: 0.601028 ± 0.000012 (0.601028); ε−1: 7.4231 ± 0.00024 (7.423055); ε0: 44.9127 ±
0.00073 (44.91255); ε1: 217.023 ± 0.0037 (217.0209); ε2: 780.436 ± 0.013 (780.432); ε3:2678.13 ± 0.053 (unknown); ε4: 7195.9 ± 0.3 (unknown).

1Please, note that the overall normalization used by FIESTA is

*different *from the one employed by the authors of

*Applications of FIESTA*
: 5.184645 ± 0.000042 (5.1846388); ε0: 38.8948 ± 0.00039 (38.8946741); ε1: 240.069
± 0.0019 (240.0684359); ε2: 948.623 ± 0.0091 (unknown); ε3: 3679.77± 0.06 (unknown).

: 20.73860 ± 0.00023(20.7385551); ε0: 102.033 ± 0.003 (102.0326759); ε1: 761.60 ±
0.011 (761.5969858); ε2: 2326.18 ± 0.062 (unknown); ε3: 12273.6 ± 0.4 (unknown).

: 20.73860 ± 0.00024 (20.7385551); ε0: 145.381 ± 0.0029 (145.3808999); ε1: 985.91
± 0.014 (985.9082306); ε2: 3930.65 ± 0.076 (unknown); ε3: 17486.6 ± 0.6 (unknown).

ε0: 55.58537 ± 0.00031 (55.5852539); ε1: 175.325 ± 0.004 (unknown); ε2: 1496.52 ± 0.02
: 52.0181 ± 0.0003 (52.0178687); ε1: 175.50 ± 0.0036 (175.496447); ε2: 1475.272 ±
0.0098 (unknown); ε3: 2623.5 ± 0.1 (unknown).

: -5.184651 ± 0.000048 (-5.184638); ε0: -32.0962 ± 0.00057 (-32.09614); ε1: -91.158
± 0.0052 (-91.1614); ε2: 119.06 ± 0.043 (unknown); ε3: 2768.6 ± 0.45 (unknown).

: 20.73857 ± 0.00026 (20.7385551); ε1: 190.60 ± 0.0023 (190.600238); ε2: 1049.20 ±
0.014 (1049.194196); ε3: 4423.84 ± 0.072 (unknown); ε4: 16028.8 ± 0.5 (unknown).

: -10.36931 ± 0.00006 (-10.3692776); ε0: -70.990 ± 0.0011 (-70.99081719); ε1: -
21.650 ± 0.013 (-21.663005); ε2: 2832.69 ± 0.096 (unknown).

: -10.36933 ± 0.00006 (-10.36927); ε0: -58.6187 ± 0.0013(-58.6210); ε1: 244.681 ±
0.015 (unknown).

: -5.18467 ± 0.000042 (-5.184638); ε0: 14.3989 ± 0.00081 (14.39739); ε1: 739.979 ±
0.0099 (unknown).

Here for each MI we provide our numerical result for coefficients of ε-expansion in comparison (inparentheses) with the known from [analiycal results (if any). As we can see, our calculations
reproduce the result of ] with 3-4 correct digits. The extra terms in the ε-expansion of each MIwhich are currently unavailable analytically but are necessary for future five-loop calculations.

Usually, analytical evaluation of multiloop MI is a kind of art. It requires a lot of efforts (and
CPU time). In many situations, independent checkup is hardly any possible in reasonable time.

That is why the simple in use tools for numerical evaluation like FIESTA are important.

*Acknowledgments. *This work was supported in part by DFG through SBF/TR 9 and the Rus-
sian Foundation for Basic Research through grant 08-02-01451.

[1] K. Hepp, Commun. Math. Phys.

**2 **(1966) 301.

*Applications of FIESTA*
[2] E.R. Speer, J. Math. Phys.,

**9 **(1968) 1404; M.C. Bergère and J.B. Zuber, Commun. Math. Phys.

**35**
(1974) 113; M.C. Bergère and Y.M. Lam, J. Math. Phys.

**17 **(1976) 1546; O.I. Zavialov,

*Renormalized*

quantum field theory, Kluwer Academic Publishers, Dodrecht (1990); V.A. Smirnov, Commun. Math.

Phys.

**134 **(1990) 109.

[3] P. Breitenlohner and D. Maison, Commun. Math. Phys.

**52 **(1977) 11; 39,55;

[4] M.C. Bergère, C. de Calan and A.P.C. Malbouisson, Commun. Math. Phys.

**62 **(1978) 137;

K. Pohlmeyer, J. Math. Phys.

**23 **(1982) 2511.

[5] V.A. Smirnov,

*Applied asymptotic expansions in momenta and masses*, STMP

**177**, Springer, Berlin,

Heidelberg (2002).

[6] T. Binoth and G. Heinrich, Nucl. Phys. B,

**585 **(2000) 741; Nucl. Phys. B,

**680 **(2004) 375; Nucl.

Phys. B,

**693 **(2004) 134.

[7] G. Heinrich, Int. J. of Modern Phys. A,

**23 **(2008) 10. [arXiv:0803.4177].

[8] C. Bogner and S. Weinzierl, Comput. Phys. Commun.

**178 **(2008) 596 [arXiv:0709.4092 [hep-ph]];

Nucl. Phys. Proc. Suppl.

**183 **(2008) 256 [arXiv:0806.4307 [hep-ph]].

[9] A.V. Smirnov and M.N. Tentyukov, Comput. Phys. Commun.

**180 **(2009) 735 [arXiv:0807.4129

[10] A.V. Smirnov, V.A. Smirnov and M.N. Tentyukov,[arXiv:0912.0158 [hep-ph]].

[11] A. V. Smirnov, V. A. Smirnov and M. Steinhauser, Phys. Lett. B

**668**, 293 (2008) [arXiv:0809.1927

[hep-ph]]; R. Bonciani and A. Ferroglia, JHEP

**0811**, 065 (2008) [arXiv:0809.4687 [hep-ph]];

Y. Kiyo, D. Seidel and M. Steinhauser, JHEP

**0901**, 038 (2009) [arXiv:0810.1597 [hep-ph]];G. Bell,

Nucl. Phys. B

**812**, 264 (2009) [arXiv:0810.5695 [hep-ph]];V. N. Velizhanin, arXiv:0811.0607

[hep-th];T. Ueda and J. Fujimoto, arXiv:0902.2656 [hep-ph];D. Seidel, arXiv:0902.3267

[hep-ph];;G. Heinrich, T. Huber, D. A. Kosower and V. A. Smirnov, Phys. Lett. B

**678**, 359 (2009)

[arXiv:0902.3512 [hep-ph]];P. A. Baikov, K. G. Chetyrkin, A. V. Smirnov, V. A. Smirnov and

M. Steinhauser, Phys. Rev. Lett.

**102**, 212002 (2009) [arXiv:0902.3519 [hep-ph]];J. Gluza, K. Kajda,

T. Riemann and V. Yundin, PoS A

**CAT08**, 124 (2008) [arXiv:0902.4830 [hep-ph]];S. Bekavac,

A. G. Grozin, D. Seidel and V. A. Smirnov, Nucl. Phys. B

**819**, 183 (2009) [arXiv:0903.4760

[hep-ph]];R. Bonciani, A. Ferroglia, T. Gehrmann and C. Studerus, JHEP

**0908**, 067 (2009)

[arXiv:0906.3671 [hep-ph]];M. Czakon, A. Mitov and G. Sterman, Phys. Rev. D

**80**, 074017 (2009)

[arXiv:0907.1790 [hep-ph]];A. Ferroglia, M. Neubert, B. D. Pecjak and L. L. Yang, arXiv:0907.4791

[hep-ph]; JHEP

**0911**, 062 (2009) [arXiv:0908.3676 [hep-ph]];S. Bekavac, A. G. Grozin, P. Marquard,

J. H. Piclum, D. Seidel and M. Steinhauser, arXiv:0911.3356 [hep-ph];M. Dowling, J. Mondejar,

J. H. Piclum and A. Czarnecki, arXiv:0911.4078 [hep-ph];A. V. Smirnov, V. A. Smirnov and

M. Steinhauser, [arXiv:0911.4742 [hep-ph]].

[12] A.V. Smirnov, M. Tentyukov, [arXiv:1004.1149 [hep-ph]].

[13] P.A. Baikov and K.G. Chetyrkin, [arXiv:1004.1153 [hep-ph]].

[14] G. 't Hooft and M. Veltman, Nucl. Phys. B 44 (1972) 189;C.G. Bollini and J.J. Giambiagi, Nuovo
Cim. 12 B (1972) 20.

[15] T. Kaneko and T. Ueda, [arXiv:0908.2897 [hep-ph]].

[16] T. Hahn, Comput. Phys. Commun. 168 (2005) 78, arXiv: hep-ph/0404043
[17] Wolfram Research, Inc., Mathematica, Version 7.0, Champaign, IL (2008).

[18] P. A. Baikov, Phys. Lett. B

**634**, 325 (2006) [arXiv:hep-ph/0507053].

Source: http://pos.sissa.it/archive/conferences/093/081/ACAT2010_081.pdf

Substance Withdrawal ManagementGuidelines for medical and nursing practitioners in primary health, specialist addiction, custodial and general hospital settings Substance Withdrawal ManagementGuidelines for medical and nursing practitioners in primary health, specialist addiction, custodial and general hospital settings The main authors of this guideline were Ashley Koning, Michelle Fowler, Moira Gilmour and Vanessa Caldwell. Contributions by the members of the Matua Rak_i team and the Withdrawal Management Reference Group were greatly appreciated.

tce CHEMENGERS WHO CHANGED THE WORLD EFINING, reactor design, bulk plastics – the kind of world- " no one initially realised that progesterone, changing developments featured especially when combined with oestrogen, thus far in this series include many of the stereotypical activities one would imagine