Acessibilidade / Reportar erro

A local non-Abelian gauge invariant action stemming from the nonlocal operator FmuN(D²)-1FmuN

Abstract

We report on the nonlocal gauge invariant operator of dimension two, FµN (D²)-1 FµN. We are able to localize this operator by introducing a suitable set of (anti)commuting antisymmetric tensor fields. Starting from this, we succeed in constructing a local gauge invariant action containing a mass parameter, and we prove the renormalizability to all orders of perturbation theory of this action in the linear covariant gauges using the algebraic renormalization technique. We point out the existence of a nilpotent BRST symmetry. Despite the additional (anti)commuting tensor fields and coupling constants, we prove that our model in the limit of vanishing mass is equivalent with ordinary massless Yang-Mills theories by making use of an extra symmetry in the massless case. We also present explicit renormalization group functions at two loop order in the ${\overline{\mbox{MS}}}$ scheme.

Yang-Mills gauge theory; BRST symmetry; Renormalization; Mass


A local non-Abelian gauge invariant action stemming from the nonlocal operator FmN(D2)-1FmN

D. DudalI,* * Postdoctoral fellow of the Special Research Fund of Ghent University. ; M. A. L. CapriII; J. A. GraceyIII; V. E. R. LemesII; R. F. SobreiroII; S. P. SorellaII,† † Work supported by FAPERJ, Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro, under the program Cientista do Nosso Estado, E- 26/151.947/2004. ; H. VerscheldeI

IGhent University Department of Mathematical Physics and Astronomy, Krijgslaan 281-S9, B-9000 Gent, Belgium

IIUERJ - Universidade do Estado do Rio de Janeiro, Rua São Francisco Xavier 524, 20550-013 Maracanã, Rio de Janeiro, Brazil

IIITheoretical Physics Division, Department of Mathematical Sciences, University of Liverpool, P.O. Box 147, Liverpool, L69 3BX, United Kingdom

ABSTRACT

We report on the nonlocal gauge invariant operator of dimension two, FµN (D2)-1FµN. We are able to localize this operator by introducing a suitable set of (anti)commuting antisymmetric tensor fields. Starting from this, we succeed in constructing a local gauge invariant action containing a mass parameter, and we prove the renormalizability to all orders of perturbation theory of this action in the linear covariant gauges using the algebraic renormalization technique. We point out the existence of a nilpotent BRST symmetry. Despite the additional (anti)commuting tensor fields and coupling constants, we prove that our model in the limit of vanishing mass is equivalent with ordinary massless Yang-Mills theories by making use of an extra symmetry in the massless case. We also present explicit renormalization group functions at two loop order in the ${\overline{\mbox{MS}}}$ scheme.

Keywords: Yang-Mills gauge theory; BRST symmetry; Renormalization; Mass

I. INTRODUCTION

We shall consider pure Euclidean SU(N) Yang-Mills theories with action

where , a = 1,...,N2 - 1 is the gauge boson field, with associated field strength

The theory (1) is invariant with respect to the local gauge transformations

with

denoting the adjoint covariant derivative.

As it is well known, the theory (1) is asymptotically free [1, 2], i.e. the coupling becomes smaller at higher energies and vice versa. At very high energies, the interaction is weak and the gluons can be considered as almost free particles. However, in spite of the progress in the last decades, we still lack a satisfactory understanding of the behaviour of Yang-Mills theories in the low energy regime. Here the coupling constant of the theory is large and nonperturbative effects have to be taken into account.

The introduction of condensates, i.e. the (integrated) vacuum expectation value of certain operators, allows one to parametrize certain nonperturbative effects arising from the infrared sector of e.g. the theory described by (1). Via the Operator Product Expansion (OPE) (viz. short distance expansion), which is applicable to local operators, one can relate these condensates to power corrections which give nonperturbative information in addition to the perturbatively calculable results. If one wants to consider the possible effects of condensates on physical quantities in a gauge theory, quite clearly only gauge invariant operators should be considered. The most famous example is the dimension 4 gluon condensate áasñ, giving rise to power corrections. Via the SVZ (Shifman-Vainshtein-Zakharov) sum rules [3], one can extract phenomenological estimates for áasñ.

In recent years, a great deal of interest arose in dimension 2 condensates in gauge theories. Most attention was paid to the gluon condensate áñ in the Landau gauge, due to the work of [4, 5], as the quantity

which is gauge invariant due to the minimization along the gauge orbits, could be physically relevant. In fact, as shown in [4, 5] in the case of compact QED, the quantity áñ seems to be useful in order to detect the presence of nontrivial field configurations like monopoles. One can show that can be written as an infinite series of nonlocal terms, see [6, 7] and references therein, namely

Since the operator is nonlocal, it falls beyond the applicability of the OPE annex sum rules, which refer to local operators.

However, in the Landau gauge, ¶m Aµ = 0, all nonlocal terms of expression (6) drop out, so that reduces to the local operator , hence the interest in the Landau gauge and its dimension two gluon condensate áñ. A complication is that the explicit determination of the absolute minimum of along its gauge orbit, and moreover of its vacuum expectation value, is a very delicate issue intimately related to the problem of Gribov copies [8-13].

Nevertheless, some nontrivial results were proven concerning the operator . In particular, we mention its multiplicative renormalizability to all orders of perturbation theory, in addition to an interesting and numerically verified relation concerning its anomalous dimension [14, 15]. An effective potential approach consistent with renormalizability and renormalization group requirements for local composite operators (LCO) has also been worked out for this operator, giving further evidence of a nonvanishing condensate áñ ¹ 0, which lowers the nonperturbative vacuum energy [16]. The LCO method yields an effective gluon mass squared ~ áñ of a few hundred MeV [16-19].

In [20], it was already argued that gauge (in)variant condensates could also influence gauge variant quantities such as the gluon propagator. An OPE argument based on lattice simulations in the Landau gauge has indeed provided evidence that a condensate áñ could account for quadratic power corrections of the form ~ , reported in the running of the coupling constant as well as in the gluon propagator, see e.g. [21-26]. This OPE approach allows one to obtain an estimate for the soft part áñIR originating from the infrared sector. The OPE can also be employed to relate this condensate to an effective gluon mass [21].

The presence of mass parameters in the gluon propagator have also been advocated from the lattice perspective several times [27-32], whilst effective gluon masses also found phenomenological use [33, 34].

A somewhat weak point about the operator is that it is unclear how to deal with it in gauges other than the Landau gauge. Till now, it seems hopeless to prove its renormalizability out of the Landau gauge. In fact, at the classical level, adding (5) to the Yang-Mills action is equivalent to add the so-called Stueckelberg action, which is known to be not renormalizable [35, 36]. We refer to [7] for details and references. As already mentioned, also the OPE becomes useless outside the Landau gauge for this particular operator.

In other gauges, there can be found other renormalizable local operators, which condense and give rise to a dynamical gluon mass. Next to the Landau gauge [14, 16-19, 37, 38] the maximal Abelian [39, 40], linear covariant [41-43] and Curci-Ferrari gauges [44, 45] have been investigated in the past.

The relevant operators in these other gauges are however gauge variant, and as a consequence also the effective gluon mass. From this perspective, it is worthwhile to find out whether a gauge invariant framework might be found for a dynamical mass, and related to it for power corrections.

In order to have a starting point, we need a dimension 2 operator that is gauge invariant. This necessarily implies a nonlocal operator, since gauge invariant local operators of dimension 2 do not exist. We would also need a consistent calculational framework, which requires an action only containing local terms. Therefore, we should find an operator that can be localized by means of a finite set of auxiliary fields, in such a way that the local gauge invariance is respected. As looks a bit hopeless from this viewpoint as it is a infinite series of nonlocal terms (6), we moved our attention instead to the nonlocal gauge invariant operator

This operator caught already some attention in 3 dimensional gauge theories in relation to a dynamical mass generation [46].

In the following sections, we shall show that the operator (7) can be localized, giving rise to a local, classically gauge invariant action. Afterwards, we discuss how to investigate the renormalizability of the action once quantized. Eventually, we need to introduce a slightly more general classical action in order to obtain a quantum action that is renormalizable to all orders of perturbation theory. In the case of vanishing mass, the equivalence of our action with usual Yang-Mills theories can be shown. We shall point out the existence of a naturally extended version of the usual BRST symmetry. Before turning to conclusions, we explicitly give various renormalization group functions, verifying the renormalizability at the practical level.

II. CONSTRUCTION OF THE ACTION AND ITS RENORMALIZABILITY AT THE QUANTUM LEVEL

A. The action at the classical level

We can add the operator (7) to the Yang-Mills action as a mass term via

with

As we have discussed in [7], the action (8) can be localized by introducing a pair of complex bosonic antisymmetric tensor fields, (), and a pair of complex anticommuting antisymmetric tensor fields, (), belonging to the adjoint representation, according to which

Therefore, we obtain a classical local action which reads

where

which is left invariant by the gauge transformations

B. The action at the quantum level

In order to discuss the renormalizability of (11), we relied on a method introduced by Zwanziger in [47, 48]. Instead of using (11) with m coupled to the composite operators and , we introduce 2 suitable external sources Vrsmn and rsmn and replace Sm by

At the end, the sources Vsrmn(x), srmn(x) are required to attain their physical value, namely

so that (14) reduces to Sm in the physical limit.

From now on, we assume the linear covariant gauge fixing, implemented through

In [7], we wrote down a list of symmetries enjoyed by the action

i.e. in absence of the sources. Let us only mention here the BRST symmetry, generated by the nilpotent transformation s given by

It turns out that one can introduce all the necessary external sources in a way consistent with the starting symmetries. This allows to write down several Ward identities by which the most general counterterm is restricted using the algebraic renormalization formalism [7, 49]. After a very cumbersome analysis, it turns out that the action (11) must be modified to

with

in order to have renormalizability to all orders of perturbation theory. We notice that we had to introduce a new invariant quartic tensor coupling labcd, subject to the generalized Jacobi identity

and the symmetry constraints

as well as two new mass couplings l1 and l3. Without the new couplings, i.e. when l1º 0, l3º 0, labcdº 0, the previous action would not be renormalizable. We refer to [7, 50] for all the details. We also notice that the novel fields , and are no longer appearing at most quadratically. As it should be expected, the classical action Scl is still gauge invariant w.r.t. (13).

III. FURTHER PROPERTIES OF THE ACTION

A. Existence of a BRST symmetry with a nilpotent charge

The BRST transformation (18) no longer generates a symmetry of the action Sphys. However, we are able to define a natural generalization of the usual BRST symmetry that does constitute an invariance of the gauge fixed action (19). Indeed, after inspection, one shall find that

with

Hence, the action Sphys is invariant with respect to a nilpotent BRST transformation . We obtained thus a gauge field theory, described by the action Sphys, (19), containing a mass term, and which has the property of being renormalizable, while nevertheless a nilpotent BRST transformation expressing the gauge invariance after gauge fixing exists simultaneously. It is clear that stands for the usual BRST transformation, well known from literature, on the original Yang-Mills fields, whereas the gauge fixing part Sgf given in (20) can be written as a -variation, ensuring that the gauge invariant physical operators shall not depend on the choice of the gauge parameter [49].

B. Existence of a ''supersymmetry'' when m º 0

We define a nilpotent (anti-commuting) transformation ds as

Then one easily verifies that (25) generates a ''supersymmetry'' of the action since

with

Taking another look at the transformations s and , respectively given by (18) and (24), one recognizes that

Since ds is a nilpotent operator, it possesses its own cohomology, which is easily identified with polynomials in the original Yang-Mills fields {, , , ª}. The auxiliary tensor fields, {, }, do not belong to the cohomology of ds, because they form pairs of doublets [49]. This fact can be brought to use in the following subsection.

C. Equivalence with Yang-Mills gauge theory when m º 0

If the mass m º 0, we would expect that the action (19) would be equivalent with the usual Yang-Mills gauge theory, since in the nonlocal formulation (9), we would have introduced ''nothing''. In the local renormalizable formulation (19), this would also be trivially true when labcdº 0 as then we would only have added a -although quite complicated- unity to the Yang-Mills action. Unfortunately, since renormalization forbids setting labcd = 0, we must find another argument to relate to the usual Yang-Mills gauge theory. As proven in [50], the ''supersymmetry'' ds of (25) can be used to show that

where

is a generic Yang-Mills functional. The expectation value of any Yang-Mills Green function, constructed from the fields {,,a,} and calculated with the original (gauge fixed) Yang-Mills action SYM + Sgf, is thus identical to the one calculated with the massless action , where it is of course assumed that the gauge freedom of both actions has been fixed by an identical gauge fixing.

The foregoing result also reflects on the renormalization group functions. As usual, we employ a massless renormalization scheme known as the scheme. As a consequence, we can set m = 0 in order to extract the ultraviolet behaviour. Using (29), we conclude that all the renormalization group functions of the original Yang-Mills quantities are not affected by the presence of the extra fields or couplings. This fact shall be explicitly verified in the next section.

IV. EXPLICIT RENORMALIZATION AT TWO LOOP ORDER

Having proven the renormalizability of the action (19), we shall now compute explicitly the two loop anomalous dimension of the fields and the one loop b-function of the tensor coupling labcd. The corresponding details can be found in [7] for one loop results, while two loop results are discussed in [50].

We have regarded the mass operator as an insertion and split the Lagrangian into a free piece involving massless fields with the remainder being transported to the interaction Lagrangian. To renormalize the operator, we insert it into a massless Green function, after the fields and couplings have been renormalized in the massless Lagrangian. An attractive feature of the massless field approach is that we can use the MINCER algorithm to perform the actual computations. This algorithm, [51], written in the symbolic manipulation language FORM, [52, 53], is devised to extract the divergences from massless 2-point functions. The propagators of the massless fields in an arbitrary linear covariant gauge are

where p is the momentum. The necessary Feynman diagrams were generated automatically with QGRAF [54].

We first checked that the same two loop anomalous dimensions emerge for the gluon, Faddeev-Popov ghost and quarks in an arbitrary linear covariant gauge as when the extra localizing fields are absent. It was also explicitly verified that the correct coupling constant renormalization constant is found. These results are in agreement with the general argument of the previous subsection.

We have implemented the properties (21) and (22) of the labcd coupling in a FORM module, while it was assumed that

which follows from the fact that there is only one rank 2 invariant tensor in a classical Lie group.

At two loops in the scheme, we find that

where NA is the dimension of the adjoint representation of the colour group, a = and we have also absorbed a factor of into labcd here and in later anomalous dimensions. These anomalous dimensions are consistent with the general observation that these fields must have the same renormalization constants, in agreement with the output of the Ward identities [7]. A check on (33) is that after the renormalization of the 3-point gluon vertex, the correct gauge parameter independent coupling constant renormalization constant emerges.

We also determined the one loop b-function for the labcd couplings. As this is present in a quartic interaction it means that to deduce its renormalization constant, we need to consider a 4-point function. However, in such a situation the MINCER algorithm is not applicable since two external momenta have to be nullified and this will lead to spurious infrared infinities which could potentially corrupt the renormalization constant. Therefore, for this renormalization only, we have resorted to using a temporary mass regularization introduced into the computation using the algorithm of [55] and implemented in Form. Consequently, we find the gauge parameter independent anomalous dimension

from both the and vertices where is the totally symmetric rank four tensor defined by

with denoting the group generator in the adjoint representation, [56]. Producing the same expression for both these 4-point functions, aside from the gauge independence, is a strong check on their correctness as well as the correct implementation of the group theory. Moreover, as it should be, babcd enjoys the same symmetry properties as the tensor labcd, summarized in (22).

We notice that labcd = 0 is not a fixed point due to the extra labcd-independent terms. If we had not included the labcd-interaction term in the original action, then such a term would inevitably be generated at one loop through quantum corrections, meaning that in this case there would have been a breakdown of renormalizability.

Finally, we turn to the two loop renormalization of the mass m. The corresponding operator can be read off from (19) and is given by

We insert into a - 2-point function and deduce the appropriate renormalization constant, leading to the anomalous dimension

as the two loop anomalous dimension. The gauge parameter independence is again a good check, as the operator is gauge invariant.

V. CONCLUSIONS

We added a nonlocal mass term (9) to the Yang-Mills action (1), and starting from this, we succeeded in constructing a renormalizable massive gauge model, which is gauge invariant at the classical level and when quantized it enjoys a nilpotent BRST symmetry. This BRST symmetry ensures that the expectation value of gauge invariant operator is gauge parameter independent. We have also proven the equivalence of the massless version of our model with Yang-Mills gauge theories making use of a ''supersymmetry'' existing between the extra fields in that case. We presented explicit two loop renormalization functions, thereby verifying that the anomalous dimensions of the original Yang-Mills quantities remain unchanged.

Many things could be investigated in the future concerning the gauge model described by (19).

At the perturbative level, it could be investigated which (asymptotic) states belong to a physical subspace of the model, and in addition one should find out whether this physical subspace can be endowed with a positive norm, which would imply unitarity. The nilpotent BRST symmetry (24) might be useful for this.

The model (19) is also asymptotically free, implying that at low energies nonperturbative effects, such as confinement, could set in. Proving and understanding the possible confinement mechanism in this model is probably as difficult as for usual Yang-Mills gauge theories.

It would also be interesting to find out whether a dynamically generated term m(-B)F might emerge, which in turn could influence the gluon Green functions. This might also be relevant in the context of gauge invariant power corrections, an issue that recently has also attracted attention from the gauge/string duality side, the so-called AdS/QCD [57, 58].

Acknowledgments

D. Dudal and H. Verschelde are grateful to the organizers of the stimulating IRQCD06 conference. The Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq-Brazil), the Faperj, Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro, the SR2-UERJ and the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) are gratefully acknowledged for financial support.

[9] Semenov-Tyan-Shanskii and V.A. Franke, Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V.A. Steklov AN SSSR, Vol. 120, 159 (1982). English translation: New York: Plenum Press 1986.

[25] S. Furui and H. Nakajima, hep-lat/0503029.

[52] J. A. M. Vermaseren, math-ph/0010025.

[58] C. Csaki and M. Reece, hep-ph/0608266.

Received on 12 September, 2006

  • [1] D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973)
  • [2] H. D. Politzer, Phys. Rev. Lett. 30, 1346 (1973).
  • [3] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979).
  • [4] F. V. Gubarev, L. Stodolsky, and V. I. Zakharov, Phys. Rev. Lett. 86, 2220 (2001).
  • [5] F. V. Gubarev, and V. I. Zakharov, Phys. Lett. B 501, 28 (2001).
  • [6] M. Lavelle and D. McMullan, Phys. Rept. 279, 1 (1997).
  • [7] M. A. L. Capri, D. Dudal, J. A. Gracey, V. E. R. Lemes, R. F. Sobreiro, S. P. Sorella, and H. Verschelde, Phys. Rev. D 72, 105016 (2005).
  • [8] V. N. Gribov, Nucl. Phys. B 139, 1 (1978).
  • [10] D. Zwanziger, Nucl. Phys. B 345, 461 (1990).
  • [11] G. Dell'Antonio and D. Zwanziger, Nucl. Phys. B 326, 333 (1989).
  • [12] G. Dell'Antonio and D. Zwanziger, Commun. Math. Phys. 138, 291 (1991).
  • [13] P. van Baal, Nucl. Phys. B 369, 259 (1992).
  • [14] D. Dudal, H. Verschelde, and S. P. Sorella, Phys. Lett. B 555, 126 (2003).
  • [15] J. A. Gracey, Phys. Lett. B 552, 101 (2003).
  • [16] H. Verschelde, K. Knecht, K. Van Acoleyen, and M. Vanderkelen, Phys. Lett. B 516, 307 (2001).
  • [17] R. E. Browne and J. A. Gracey, Phys. Lett. B 597, 368 (2004).
  • [18] J. A. Gracey, Eur. Phys. J. C 39, 61 (2005).
  • [19] R. E. Browne and J. A. Gracey, JHEP 0311, 029 (2003).
  • [20] M. J. Lavelle and M. Schaden, Phys. Lett. B 208, 297 (1988).
  • [21] K. I. Kondo, Phys. Lett. B 514, 335 (2001).
  • [22] P. Boucaud, A. Le Yaouanc, J. P. Leroy, J. Micheli, O. Pene, and J. Rodriguez-Quintero, Phys. Rev. D 63, 114003 (2001).
  • [23] P. Boucaud, J. P. Leroy, A. Le Yaouanc, J. Micheli, O. Pene, F. De Soto, A. Donini, H. Moutard, and J. Rodriguez-Quintero, Phys. Rev. D 66, 034504 (2002).
  • [24] E. Ruiz Arriola, P. O. Bowman, and W. Broniowski, Phys. Rev. D 70, 097505 (2004).
  • [26] S. Furui and H. Nakajima, Phys. Rev. D 73, 074503 (2006).
  • [27] P. Marenzoni, G. Martinelli, and N. Stella, Nucl. Phys. B 455, 339 (1995).
  • [28] D. B. Leinweber, J. I. Skullerud, A. G. Williams, and C. Parrinello [UKQCD Collaboration], Phys. Rev. D 60, 094507 (1999) [Erratum-ibid. D 61 (2000) 079901].
  • [29] K. Langfeld, H. Reinhardt, and J. Gattnar, Nucl. Phys. B 621, 131 (2002).
  • [30] K. Amemiya and H. Suganuma, Phys. Rev. D 60, 114509 (1999).
  • [31] V. G. Bornyakov, M. N. Chernodub, F. V. Gubarev, S. M. Morozov, and M. I. Polikarpov, Phys. Lett. B 559, 214 (2003).
  • [32] A. Cucchieri, T. Mendes, and A. R. Taurines, Phys. Rev. D 71, 051902 (2005).
  • [33] G. Parisi and R. Petronzio, Phys. Lett. B 94, 51 (1980).
  • [34] J. H. Field, Phys. Rev. D 66, 013013 (2002).
  • [35] H. Ruegg and M. Ruiz-Altaba, Int. J. Mod. Phys. A 19, 3265 (2004).
  • [36] H. van Dam and M. J. G. Veltman, Nucl. Phys. B 22, 397 (1970).
  • [37] R. F. Sobreiro, S. P. Sorella, D. Dudal, and H. Verschelde, Phys. Lett. B 590, 265 (2004).
  • [38] D. Dudal, R. F. Sobreiro, S. P. Sorella, and H. Verschelde, Phys. Rev. D 72, 01401 (2005).
  • [39] D. Dudal, J. A. Gracey, V. E. R. Lemes, M. S. Sarandy, R. F. Sobreiro, S. P. Sorella, and H. Verschelde, Phys. Rev. D 70, 114038 (2004).
  • [40] M. A. L. Capri, V. E. R. Lemes, R. F. Sobreiro, S. P. Sorella, and R. Thibes, Phys. Rev. D 72, 085021 (2005).
  • [41] D. Dudal, H. Verschelde, V. E. R. Lemes, M. S. Sarandy, R. F. Sobreiro, S. P. Sorella, and J. A. Gracey, Phys. Lett. B 574, 325 (2003).
  • [42] D. Dudal, H. Verschelde, J. A. Gracey, V. E. R. Lemes, M. S. Sarandy, R. F. Sobreiro, and S. P. Sorella, JHEP 0401, 044 (2004).
  • [43] R. F. Sobreiro and S. P. Sorella, JHEP 0506, 054 (2005).
  • [44] D. Dudal, H. Verschelde, V. E. R. Lemes, M. S. Sarandy, S. P. Sorella, and M. Picariello, Annals Phys. 308, 62 (2003).
  • [45] D. Dudal, H. Verschelde, V. E. R. Lemes, M. S. Sarandy, R. F. Sobreiro, S. P. Sorella, M. Picariello, and J. A. Gracey, Phys. Lett. B 569, 57 (2003).
  • [46] R. Jackiw and S. Y. Pi, Phys. Lett. B 403, 297 (1997).
  • [47] D. Zwanziger, Nucl. Phys. B 323, 513 (1989).
  • [48] D. Zwanziger, Nucl. Phys. B 399, 477 (1993).
  • [49] O. Piguet and S. P. Sorella, Lect. Notes Phys. M28, 1 (1995).
  • [50] M. A. L. Capri, D. Dudal, J. A. Gracey, V. E. R. Lemes, R. F. Sobreiro, S. P. Sorella, and H. Verschelde, Phys. Rev. D 74, 045008 (2006).
  • [51] S. G. Gorishnii, S. A. Larin, L. R. Surguladze, and F. V. Tkachov, Comput. Phys. Commun. 55, 381 (1989).
  • [53] S. A. Larin, F. V. Tkachov, and J. A. M. Vermaseren, The FORM version of MINCER, NIKHEF-H-91-18.
  • [54] P. Nogueira, J. Comput. Phys. 105, 279 (1993).
  • [55] K. G. Chetyrkin, M. Misiak, and M. Munz, Nucl. Phys. B 518, 473 (1998).
  • [56] T. van Ritbergen, A. N. Schellekens, and J. A. M. Vermaseren, Int. J. Mod. Phys. A 14, 41 (1999).
  • [57] O. Andreev, Phys. Rev. D 73, 107901 (2006).
  • *
    Postdoctoral fellow of the Special Research Fund of Ghent University.
  • †
    Work supported by FAPERJ, Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro, under the program Cientista do Nosso Estado, E- 26/151.947/2004.
  • Publication Dates

    • Publication in this collection
      11 May 2007
    • Date of issue
      Mar 2007

    History

    • Received
      12 Sept 2006
    Sociedade Brasileira de Física Caixa Postal 66328, 05315-970 São Paulo SP - Brazil, Tel.: +55 11 3091-6922, Fax: (55 11) 3816-2063 - São Paulo - SP - Brazil
    E-mail: sbfisica@sbfisica.org.br