Next Article in Journal
Raychaudhuri Equation, Geometrical Flows and Geometrical Entropy
Next Article in Special Issue
IDM Benchmarks for the LHC and Future Colliders
Previous Article in Journal
Shock Structure and Relaxation in the Multi-Component Mixture of Euler Fluids
Previous Article in Special Issue
Analysis of the Internal Structure of Hadrons Using Direct Photon Production
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Brief Review of Implicit Regularization and Its Connection with the BPHZ Theorem

by
Dafne Carolina Arias-Perdomo
1,
Adriano Cherchiglia
1,*,
Brigitte Hiller
2 and
Marcos Sampaio
1
1
CCNH Centro de Ciências Naturais e Humanas, Universidade Federal do ABC, Santo André 09210-580, SP, Brazil
2
CFisUC, Department of Physics, University of Coimbra, P-3004-516 Coimbra, Portugal
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(6), 956; https://doi.org/10.3390/sym13060956
Submission received: 29 April 2021 / Revised: 12 May 2021 / Accepted: 20 May 2021 / Published: 27 May 2021
(This article belongs to the Special Issue Higher Order Radiative Corrections in QCD)

Abstract

:
Quantum Field Theory, as the keystone of particle physics, has offered great insights into deciphering the core of Nature. Despite its striking success, by adhering to local interactions, Quantum Field Theory suffers from the appearance of divergent quantities in intermediary steps of the calculation, which encompasses the need for some regularization/renormalization prescription. As an alternative to traditional methods, based on the analytic extension of space–time dimension, frameworks that stay in the physical dimension have emerged; Implicit Regularization is one among them. We briefly review the method, aiming to illustrate how Implicit Regularization complies with the BPHZ theorem, which implies that it respects unitarity and locality to arbitrary loop order. We also pedagogically discuss how the method complies with gauge symmetry using one- and two-loop examples in QED and QCD.

1. Introduction

In Feynman diagram calculations, the building blocks of the perturbative expansion, the transition amplitudes, contain apparent divergences at intermediate stages, and yet all physical quantities we compute in perturbation theory are expected to be finite. A technical problem is immediately posed as one has to perform algebraic operations with divergent quantities. The naive solution is simply to “regulate” the divergences or make them “manifestly finite” so they can make sense at a mathematical level and hope that somehow the physical answer is meaningful when the regularization is lifted. It is important to note that ultraviolet (UV) and infrared (IR) divergences, in the high and low energy domain, respectively, are unavoidable byproducts of the very construction of the quantum field theoretical model as an effective theory. In building Feynman diagram amplitudes using Feynman rules, the product of two distributions is no longer a distribution and is, therefore, ill-defined in the short distance (UV) limit. It does not possess a Fourier transform, for instance [1]. On the other hand, IR divergences cause massless theories and are a consequence of idealizations of the physical situation: taking the region of space–time to be infinite and supposing that massless particles can be detected with infinitely precise energy–momentum resolution. Quantum field theoretical divergences arise in other ways, for instance through the lack of convergence of the perturbation series, which at best is an asymptotic series. Despite all that, the standard model of particle physics is the best-tested physical theory ever. For instance, the theoretical and experimental deviation of the anomalous magnetic moment of the muon a μ has recently been measured to be [2]:
Δ a μ = a μ e x p a μ t h = ( 25.1 ± 5.9 ) × 10 10 .
Most calculations in the standard model including some supersymmetric extensions’ (Analytical continuation in the space–time dimensions clashes with the invariance of an action with respect to supersymmetric transformations. It only holds in general for specific values of the space-time dimension due to the fact that a necessary condition for supersymmetry is equality of Bose and Fermi degrees of freedom) are performed using the variants of dimensional regularization (DREG) although it is well-known that DREG has certain complications with the definition of γ 5 matrices (In d = 4 γ 5 = γ 5 = i γ 0 γ 1 γ 2 γ 3 ) which are only defined in the physical dimension, for instance, in d = 4 we have γ 5 = γ 5 = i γ 0 γ 1 γ 2 γ 3 . and Thus the treatment of the chiral theories is subtle in those methods. In this case, we have recently argued that care must be exercised even when we apply a framework that works essentially in the physical dimension [3,4,5]. The essence of the dimensional regularization is to extend the space–time dimensionality slightly away from the physical dimension and then take the physical limit after the actual calculation [6,7].
UV and IR divergences are conceptually different in the sense that, unlike IR infinities, UV divergences cannot be excused away for an idealization of a physical situation. UV infinities can, however, be renormalized. The renormalization program allows for unambiguously extracting numerical predictions for renormalizable quantum field theories order by order in perturbation theory by redefining the physical constants in the Lagrangian, such as masses, physical fields and coupling constants. Such a program is so successful that renormalizability has become a criterium for a sensible theory after its success in Quantum Electrodynamics. Today, we know that this is a simplistic view. A more sensible way to understand is the effective approach to the problem of UV divergences in which one has to declare explicit restrictions on the domain of energy scales of QFTs and adjust the sensitivity to high-energy phenomena with the tools of renormalization theory [8]. From the practical viewpoint, the problem is even more subtle. At the level of scattering cross-sections, local cancellation of infrared singularities between the so-called real and virtual emission processes in QED and QCD processes mixes UV and IR degrees of freedom in dimensional methods. Whilst finitude is guaranteed by the Kinoshita–Lee–Nauenberg theorem [9,10], which states that suitably defined inclusive quantities will indeed be free of singularities in the massless limit, the physical origin of such cancellations is obscured in DREG as UV and IR divergences can cancel each other out.
DREG is a powerful regularization technique that is convenient not only to make calculations but also to prove theorems to all orders in perturbation theory respecting the relevant Ward–Slavnov–Taylor identities [11,12,13,14]. Any alternative regularization to DREG should be consistent at an arbitrary order in perturbation theory in consonance with the renormalization program of absorbing the infinities into the physical parameters of the theory.
The renormalization program is mathematically established by the Bogoliubov–Parasiuk –Hepp–Zimmermann (BPHZ) theorem [15,16,17,18,19,20]. This scheme was originally developed by Bogoliubov and Parasiuk in terms of a recursive subtraction operation, often called Bogoliubov’s R-operation [15,16,17,18]. This framework makes it possible to subtract overlapping and nested UV divergences in Feynman integrals in a consistent way with perturbation theory . In the BPHZ scheme, the renormalization constants expressed by counterterms at the Lagrangian level are identified with local counterterms at the level of the integrands associated with Feynman graphs. In principle, to render an amplitude UV finite, the BPHZ scheme can be carried out without the need for regularization. In practice, one must adopt a regularization scheme to compute physical quantities. In the presence of a regulator, the BPHZ-scheme provides a consistent way to separate the potentially complicated finite parts of Feynman integrals from the divergent parts. An alternative proof for the finiteness of the renormalized Feynman Integral was given by Zimmermann, by means of the recursive Bogoliubov’s R-operation which leads to a sum over forests of graphs, giving rise to the Zimmermann’s forest formula [15,16,17,18]. The latter is an elegant and comparably simpler proof for the finiteness of Feynman Integrals and works directly in momentum space. A nice review of the BPHZ method can be found in [21].
The purpose of this review is to show that a regularization scheme called Implicit Regularization (IREG) [22], that works entirely in the physical dimension of the model can be implemented to all orders in perturbation theory. The distinguished feature of IREG is that the UV divergences are displayed as loop integrals free of external momentum dependence. While such a program is somewhat trivial at the one-loop level (as one algebraic identity at the integrand level is enough to extract basic divergent loop integrals), at higher-order in perturbation theory, this program is highly non-trivial as it must comply with the BPHZ theorem. A crucial question is whether the divergences at each order in perturbation theory can be expressed in terms of loop integrals according to Bogoliubov’s recursion formula.
In this work, we show that IREG respects unitarity and locality in the BPHZ sense to arbitrary loop order. The renormalization group functions can be obtained without explicit evaluation of the basic divergent integrals by means of a characteristic renormalization constant defined at the one-loop level. We illustrate our framework using a scalar field theory and generalize to abelian and non-abelian theories. We verify that such a program preserves symmetry evoking momentum routing invariance in the loops of Feynman diagrams, which defines a constrained scheme where surface terms are set to vanish.

2. IREG and the BPHZ Algorithm

For simplicity, in the context of this review, we will only consider massless theories and integrals that are free from infrared divergences. We will also be restricted to a space–time with 2 n dimensions, where n is an integer. In general, once a N-loop amplitude of a Feynman graph with L external legs is known, the strategy of IREG is to remove all external momenta from UV divergences, expressing them as a linear combination of basic divergent integrals with one loop momentum only. To fulfill this objective, we are required to perform ( N 1 ) integrations, even though the order in which they must be realized is not immediately clear. In [23], we proposed a systematic procedure to categorize the order of integration which, as a byproduct, displays automatically the counterterms to be subtracted by Bogoliubov’s recursion formula. In this work, we will review its main steps. In order to settle the notation, consider that the integral in k l is the l-th, and to deal with whose general form is given by
I ν 1 ν m = k l A ν 1 ν m ( k l , q i ) i [ ( k l q i ) 2 μ 2 ] ln l 1 k l 2 μ 2 λ 2 ,
where l = 1 N , k l d 2 n k l / ( 2 π ) 2 n (for n integer), q i is an element (or combination of elements) of the set { p 1 , , p L , k l + 1 , , k N } , and μ 2 is an infrared regulator.
We recall that only infrared safe amplitudes are considered, which implies that the limit μ 2 0 is well-defined for the amplitude as a whole, to be taken as the last step of our calculation. The symbol λ stands for an arbitrary non-vanishing parameter that will play the role of the renormalization group scale in the context of IREG. It first appears at one-loop level, surviving to higher-orders due to the use of Equation (10) as we discuss at the end of this section. The function A ν 1 ν m ( k l , q i ) may contain constants and all possible combinations of k l and q i compatible with the Lorentz structure. In the context of gauge theories, it would come from derivative couplings, Dirac algebra, etc. We argue in [24] that the form of A ν 1 ν m ( k l , q i ) , coming from a specific Feynman diagram, is unique after the following steps are taken:
(A)
Internal symmetry group and the usual Dirac algebra must be dealt with first. As extensively discussed in [3], identities only valid in strictly n-dimensional spaces (n integer) such as { γ 5 , γ μ } = 0 must not be used inside divergent amplitudes [3,4,5].
(B)
The requirement of numerator/denominator consistency implies that terms with internal momenta squared in the numerator must be canceled against the denominator. For instance,
k , q k 2 k 2 q 2 ( k q ) 2 | IREG = k , q 1 q 2 ( k q ) 2 | IREG ,
where we consider n = 2, k d 4 k / ( 2 π ) 4 . In the same vein, symmetric integration in divergent amplitudes cannot be enforced. That is,
k k μ 1 k μ 2 m f ( k 2 ) IREG g { μ 1 μ 2 g μ 2 m 1 μ 2 m } ( 2 m ) ! k k 2 m f ( k 2 ) IREG ,
where the curly brackets indicate symmetrisation over Lorentz indices.
After these steps, the resulting multi-loop integrand can be manipulated consistently in the framework of IREG, meaning that (1) each overall-divergent amplitude is separated into a unique finite expression plus a divergent part; (2) power-counting finite expressions are not modified; (3) linearity under the regularization operation R is preserved namely, [ a F + b G ] R = a [ F ] R + b [ G ] R , where F and G are Feynman integrals, a , b are quantities that may only depend on external momenta and/or masses, not the internal loop momenta. Therefore, they can be safely pulled out of the integral. Moreover, the UV content of A n will be cast in terms of well-defined basic divergent integrals, which need not to be evaluated as we will discuss soon.
Given that a normal form for A ν 1 ν m ( k l , q i ) was achieved, we apply the rules of IREG:
(a)
Starting at one loop (which is equivalent to set l = 1 in Equation (2)), we assume an implicit regulator which allows us to remove the external momenta dependence (encoded in p i ) from the UV divergent part of the amplitude by using the identity
1 ( k p i ) 2 μ 2 = j = 0 n i ( k ) 1 ( 1 ) j ( p i 2 2 p i · k ) j ( k 2 μ 2 ) j + 1 + ( 1 ) n i ( k ) ( p i 2 2 p i · k ) n i ( k ) ( k 2 μ 2 ) n i ( k ) ( k p i ) 2 μ 2 ,
in the propagators (for simplicity, we have defined k 1 = k ). As briefly discussed before, μ 0 is a fictitious mass (infrared regulator). It should be emphasized that, since the starting integrals are IR-safe, the infrared regulator will only be needed in intermediate steps of the calculation, canceling in the end result. Therefore, gauge invariance will not be spoiled. After the first step, we can define basic divergent integrals (BDIs) as
I log ( μ 2 ) k 1 ( k 2 μ 2 ) n
I log ( μ 2 ) k 1 ( k 2 μ 2 ) n , I log ν 1 ν 2 r ( μ 2 ) k k ν 1 k ν 2 r ( k 2 μ 2 ) r + n , I quad ( μ 2 ) k 1 ( k 2 μ 2 ) n 1 , I quad ν 1 ν 2 r ( μ 2 ) k k ν 1 k ν 2 r ( k 2 μ 2 ) r + n 1 .
(b)
BDIs with Lorentz indices ν 1 ν 2 r are systematically reduced to linear combinations of BDIs without Lorentz indices (with the same superficial degree of divergence) since we comply with invariance under shifts of the integration momenta and numerator–denominator consistency [3]. Therefore, the total derivatives with respect to the internal momenta must vanish, e.g.,
k k μ k ν ( k 2 μ 2 ) n = 2 n g μ ν 2 n I log ( μ 2 ) I log μ ν ( μ 2 ) = 0 ,
k k μ k ν ( k 2 μ 2 ) n 1 = ( n 1 ) 2 2 ( n 1 ) g μ ν I quad ( μ 2 ) I quad μ ν ( μ 2 ) = 0 .
(c)
After the last step, the divergent part of the amplitude will be given in terms of scalar BDIs only. However, since we still have to take the limit μ 0 , it can be noticed that they are ultraviolet and infrared divergent objects. To isolate these divergences defining a genuine ultraviolet divergent object, we use the identity below
I log ( μ 2 ) = I log ( λ 2 ) + b 2 n ln λ 2 μ 2 , b 2 n i ( 4 π ) n ( 1 ) n Γ ( n ) ,
which introduces λ > 0 as an arbitrary mass scale (renormalization group scale). I quad ( μ 2 ) can be chosen to vanish as μ goes to zero [25]. By adding the divergent part with the finite terms, the limit μ 0 is now well defined since the whole amplitude is power counting infrared convergent from the start. As we will present in our examples, the BDI will be absorbed in the renormalization constants [26] allowing renormalization functions to be obtained using
λ 2 I log ( λ 2 ) λ 2 = b 2 n .
At higher loop-order ( l > 1 in Equation (2)), the procedure is completely analogous. The generalization of the previous formulas are
(a)
After applying in the propagators the identity
1 ( k l q i ) 2 μ 2 = j = 0 n i ( k l ) 1 ( 1 ) j ( q i 2 2 q i · k l ) j ( k l 2 μ 2 ) j + 1 + ( 1 ) n i ( k l ) ( q i 2 2 q i · k l ) n i ( k l ) ( k l 2 μ 2 ) n i ( k l ) ( k l q i ) 2 μ 2 ,
where q i is contained in the set { p 1 , , p L , k l + 1 , , k n } , the UV divergent part of the amplitude is expressed as a linear combination of the objects below (We have already set quadratic divergent BDIs to zero, as previously discussed.)
I log ( l ) ( μ 2 ) k l 1 ( k l 2 μ 2 ) n ln l 1 k l 2 μ 2 λ 2 ,
I log ( l ) ν 1 ν 2 r ( μ 2 ) k l k l ν 1 k l ν 2 r ( k l 2 μ 2 ) r + n ln l 1 k l 2 μ 2 λ 2 .
(b)
As before, higher loop BDIs are reduced to scalar ones by vanishing the total derivatives
k k ν 1 k ν 2 k ν 2 j ( k 2 μ 2 ) n + j 1 ln l 1 ( k 2 μ 2 ) λ 2 = 0 .
For instance,
I log ( l ) μ ν ( μ 2 ) = j = 1 l 1 n j ( l 1 ) ! ( l j ) ! g μ ν 2 I log ( l j + 1 ) ( μ 2 ) .
(c)
We notice, once again, that the BDIs as defined in the last step are UV and IR divergent in the limit μ 0 . To define UV divergent terms only, we apply the identity
I log ( l ) ( μ 2 ) = I log ( l ) ( λ 2 ) b d l ln l μ 2 λ 2 + b d k = 1 n 1 n 1 k j = 1 l 1 ( 1 ) k k j ( l 1 ) ! ( l j ) ! ln l j μ 2 λ 2 .
The μ -dependence will cancel in the amplitude as a whole, since it was IR-safe from the start. As already commented, BDIs can be absorbed in renormalization constants. We take the opportunity to emphasize that a minimal, mass-independent subtraction scheme in IREG amounts to absorb only I log ( l ) ( λ 2 ) . To evaluate renormalization group constants, only derivatives of BDIs with respect to the renormalization scale λ 2 are required [27],
λ 2 I log ( l ) ( λ 2 ) λ 2 = ( l 1 ) I log ( l 1 ) ( λ 2 ) b 2 n α 2 n ( l ) ,
where n 2 , α 4 ( l ) = ( l 1 ) ! , α 6 ( 2 ) = 3 / 2 (a general formula for α 2 n ( l ) can be found in [27]).
Once the rules of IREG are settled, we return to the discussion of our initial problem: how do we identify the order in which the ( n 1 ) integrals must be performed? We will present below a systematic choice that can be done in a way to display the terms to be subtracted by Bogoliubov’s recursion formula implying that the method complies with Lorentz invariance, locality, and unitarity.
The main idea is to adapt identity (12) in such a way that the UV-divergent behaviour of the amplitude as the internal momenta goes to infinity in all possible ways can be clearly identified. For ease of the reader, we consider that q i will only denote external momenta ( p i ) and k is an arbitrary internal momentum. By using the binomial formula, ( p i 2 2 p i · k ) j can be expanded to yield
1 ( k p i ) 2 μ 2 = l = 0 2 ( n i ( k ) 1 ) f l ( k , p i ) + f ¯ ( k , p i ) ,
where we defined,
f l ( k , p i ) j = 0 l / 2 Θ ( B ) l j j ( p i 2 ) j ( 2 p i · k ) l 2 j ( k 2 μ 2 ) l + 1 j , f ¯ ( k , p i ) ( 1 ) n i ( k ) ( p i 2 2 p i · k ) n i ( k ) ( k 2 μ 2 ) n i ( k ) ( k p i ) 2 μ 2 , Θ ( x ) 0 i f x 0 1 i f x > 0 , B n i ( k ) + j l , x m a x { n Z | n x } .
As can be inspected, the terms f l ( k , p i ) behave as k ( l + 2 ) when k by construction while the value of n i ( k ) is chosen to guarantee the UV finitude of f ¯ ( k , p i ) . The above identity is the keystone of our procedure whose application to an arbitrary Feynman amplitude is summarized as follows:
  • Identify which propagators depend on the external momenta, then apply identity (19);
  • Obtain the minimum value of all n j ( k i ) necessary to guarantee the finitude of terms that contain f ¯ ( k i , p j ) as k i in all possible ways;
  • Isolate the UV-divergent terms, allowing a classification in terms of the different ways that the internal momenta approach infinity to be envisaged;
  • Use the rules of IREG, encoded in steps (a)–(c), in the terms identified in step 3 according to their classification;
  • Set aside the divergent terms that contain I log ( l ) ( λ 2 ) and apply the procedure again on the ones that do not.
After step 5, we will obtain two types of terms: or I log ( l ) ( λ 2 ) multiplies an integral or I log ( l ) ( λ 2 ) multiplies only constants and/or polynomials in the external momenta. The first set will amount to the terms to be removed by applying Bogoliubov’s recursion formula, while the latter set will be the typical divergence of the graph, i.e., after subtraction of subdivergences. As emphasized before, the procedure just envisaged will allow IREG to be applied in a systematic way, with the byproduct of identifying the terms to be removed by Bogoliubov’s recursion formula automatically [23].

3. Selected Examples

In this section, we present some selected examples, aiming to pedagogically illustrate how the renormalization procedure can be implemented in IREG up to two-loop. We begin with a scalar theory, moving to more realistic theories afterward (QED and QCD).

3.1. Scalar Theory ϕ 3

We initiate our discussion with a very simple theory, the massless ϕ 3 model defined in 6 dimensions. The choice of dimensions is justified to obtain a more interesting (renormalizable) model, in which only graphs up to three external legs are divergent [28]. The graphs with one external leg have only quadratic divergences which, from the point of view of IREG, could be kept as BDIs. However, they will always cancel out in multiplicatively renormalizable theories [29,30,31], and can be promptly dismissed in massless theories. Therefore, in order to perform the renormalization of the theory up to two-loop order, we need to consider graphs with only two or three external legs corresponding to the renormalization of the propagator and the vertex functions, respectively.
We begin with the one-loop correction for the propagator depicted in Figure 1, whose amplitude reads ( recall we are in six dimensions which implies k d 6 k ( 2 π ) 6 )
Ξ ( 1 ) g 2 2 k 1 k 2 1 ( k p ) 2 = lim μ 2 0 g 2 2 k 1 ( k 2 μ 2 ) 1 [ ( k p ) 2 μ 2 ] , k d 6 k ( 2 π ) 6 .
Notice that, following the rules of IREG, an infrared regulator was introduced in the denominators. The propagator which contains the external momentum can be rewritten in terms of f l and f ¯
Ξ ( 1 ) g 2 = 1 2 k 1 ( k 2 μ 2 ) l = 0 2 ( n ( k ) 1 ) f l ( k , p ) + f ¯ ( k , p ) ,
while n ( k ) is chosen in order to assure the finitude of the term containing f ¯ ( k , p ) . In this specific example, we find by power counting that n ( k ) > 2 , adopting n ( k ) = 3 . The divergent terms can be promptly identified by remembering that f l ( k , p ) behaves like k ( l + 2 )
  • Quadratic divergence
    k f 0 ( k , p ) ( k 2 μ 2 ) = k 1 ( k 2 μ 2 ) 2 ,
  • Linear divergence
    k f 1 ( k , p ) ( k 2 μ 2 ) = k 2 p · k ( k 2 μ 2 ) 3 ,
  • Logarithmic divergence
    k f 2 ( k , p ) ( k 2 μ 2 ) = k 1 ( k 2 μ 2 ) 3 ( 2 p · k ) 2 ( k 2 μ 2 ) p 2 = p 2 3 I log ( μ 2 ) .
For pedagogical reasons, we showed the quadratic and linear divergences, even though they vanish in the limit μ 2 0 . The remaining (UV finite) terms amount to
1 2 k f 3 ( k , p ) + f 4 ( k , p ) + f ¯ ( k , p ) ( k 2 μ 2 ) = 1 2 k 1 ( k 2 μ 2 ) 4 4 p 2 ( p · k ) + p 4 ( p 2 2 p · k ) 3 ( k p ) 2 μ 2 = p 2 b 6 6 ln p 2 μ 2 4 p 2 b 6 9 + O ( μ 2 ) .
After using Equation (10), the limit μ 2 0 is well-defined, and we finally obtain
Ξ ( 1 ) = g 2 p 2 6 I log ( λ 2 ) b 6 ln p 2 λ 2 + 8 b 6 3 .
Similarly, we obtain the amplitude of the one-loop correction for the vertex function shown in Figure 2.
In the amplitude corresponding to V ( 1 ) , h ( p 1 , p 2 ) is a function of p 1 and p 2 vanishing for p 2 = 0 .

3.1.1. Two Loops: Self-Energy Diagrams

After setting the stage with the one-loop graphs, we move to the two-loop contributions. We recall that our main aim is to perform the renormalization of the theory, which implies that only the divergent parts will be kept. Starting with the scalar propagator, the diagrams needed are given in Figure 3.
The amplitude corresponding to P A ( 2 ) is given by
Ξ A ( 2 ) i g 4 = 1 2 k 1 k 2 Δ ( k 1 ) Δ ( k 1 p ) Δ ( k 1 k 2 ) Δ ( k 2 ) Δ ( k 2 p ) , w h e r e Δ ( k i ) 1 k i 2 μ 2 .
As in the one-loop case, we start by rewriting the propagators that depend on the external momenta
k 1 k 2 Δ ( k 1 ) Δ ( k 1 k 2 ) Δ ( k 2 ) l = 0 2 ( n ( k 1 ) 1 ) f l ( k 1 , p ) + f ¯ ( k 1 , p ) m = 0 2 ( n ( k 2 ) 1 ) f m ( k 2 , p ) + f ¯ ( k 2 , p ) .
This time we have two n ( k i ) to be determined, which are chosen to guarantee the finitude of terms containing f ¯ ( k i , p ) as k i in all possible ways. We focus first on n ( k 1 ) . The terms that contain f ¯ ( k 1 , p ) can be compactly written as
k 1 k 2 Δ ( k 1 ) Δ ( k 1 k 2 ) Δ ( k 2 ) f ¯ ( k 1 , p ) Δ ( k 2 p ) .
We want to assure the finitude of the above integral as k 1 . Two cases must be considered:
(1)
Finitude as k 1 and k 2 fixed: n ( k 1 ) > 0 ,
(2)
Finitude as k 1 and k 2 : n ( k 1 ) > 2 ,
which allows us to conclude that n ( k 1 ) should be at least 3. Similarly, we obtain n ( k 2 ) = 3 .
Once the values of n ( k i ) are known, we aim to identify the divergent terms contained in (29) as k 1 and/or k 2 go to infinity. There are three possibilities. We begin analysing the case k 1 and k 2 fixed, which contain divergence terms of the type
k 1 k 2 Δ ( k 1 ) Δ ( k 2 ) Δ ( k 1 k 2 ) f l ( k 1 , p ) m = 0 4 f m ( k 2 , p ) + f ¯ ( k 2 , p ) .
Since f l ( k 1 , p ) goes like k 1 ( l + 2 ) , we find by power counting that the divergent terms are given by
A 1 Ξ k 1 k 2 Δ ( k 1 ) Δ ( k 2 ) Δ ( k 1 k 2 ) f 0 ( k 1 , p ) m = 0 4 f m ( k 2 , p ) + f ¯ ( k 2 , p ) = k 1 k 2 Δ 2 ( k 1 ) Δ ( k 1 k 2 ) Δ ( k 2 ) Δ ( k 2 p ) .
In a similar fashion, the case k 2 and k 1 amounts to
A 2 Ξ k 1 k 2 Δ ( k 1 ) Δ ( k 2 ) Δ ( k 1 k 2 ) f 0 ( k 2 , p ) l = 0 4 f l ( k 1 , p ) + f ¯ ( k 1 , p ) = k 1 k 2 Δ 2 ( k 2 ) Δ ( k 1 k 2 ) Δ ( k 1 ) Δ ( k 1 p ) .
Finally we consider k 1 and k 2 simultaneously. The choice of n ( k i ) = 3 ( i = 1 , 2 ) guarantees that the divergent terms must be of the type
k 1 k 2 Δ ( k 1 ) Δ ( k 2 ) Δ ( k 1 k 2 ) f l ( k 1 , p ) f m ( k 2 , p ) .
Once again, by power counting, we obtain that l and m are constrained by l + m 2 . Cases l = 0 and m = 0 , 1 , 2 are already contained in A 1 Ξ (Equation (32)), while cases m = 0 and l = 0 , 1 , 2 are part of A 2 Ξ (Equation (33)). Thus, we are left only with the case l = m = 1
A 3 Ξ k 1 k 2 Δ ( k 2 ) Δ ( k 1 k 2 ) Δ ( k 1 ) f 1 ( k 1 , p ) f 1 ( k 2 , p ) = k 1 k 2 Δ 3 ( k 1 ) Δ ( k 1 k 2 ) Δ 3 ( k 2 ) ( 2 p · k 1 ) ( 2 p · k 2 ) .
In summary, the divergent terms are the following:
(1)
Case k 1 and k 2 is fixed
A 1 Ξ = k 1 k 2 Δ 2 ( k 1 ) Δ ( k 1 k 2 ) Δ ( k 2 ) Δ ( k 2 p ) ,
(2)
Case k 2 and k 1 is fixed
A 2 Ξ = k 1 k 2 Δ 2 ( k 2 ) Δ ( k 1 k 2 ) Δ ( k 1 ) Δ ( k 1 p ) ,
(3)
Case k 1 and k 2 simultaneously
A 3 Ξ = k 1 k 2 Δ 3 ( k 1 ) Δ ( k 1 k 2 ) Δ 3 ( k 2 ) ( 2 p · k 1 ) ( 2 p · k 2 ) .
Thus, the divergent content of Ξ A ( 2 ) amounts to A 1 Ξ + A 2 Ξ + A 3 Ξ A 4 Ξ . The last term corresponds to the case ( l = m = 0 )
A 4 Ξ k 1 k 2 Δ ( k 2 ) Δ ( k 1 k 2 ) Δ ( k 1 ) f 0 ( k 1 , p ) f 0 ( k 2 , p ) = k 1 k 2 Δ 2 ( k 2 ) Δ ( k 1 k 2 ) Δ 2 ( k 1 ) ,
which must be subtracted since it was counted twice.
The above classification of the divergent terms in different cases (we consider A 4 Ξ as the intersection between the k 1 and k 2 fixed, k 2 and k 1 fixed) is crucial to implement IREG to multi-loop Feynman graphs in a systematic way, since it gives a natural order in which the integrals must be performed. For instance, to evaluate A 1 Ξ further, we must perform the integral in k 1 first. For terms like A 3 Ξ , which are symmetric under k 1 k 2 , one may perform any of the integrals first. We should also emphasize that, as a byproduct, this classification will display the terms to be subtracted by Bogoliubov’s recursion formula automatically as we will show.
Returning to the divergent terms we classified, one may notice that A 1 Ξ and A 2 Ξ have the same structure, and the integral to be dealt with first is
k i Δ 2 ( k i ) Δ ( k i k j ) , i , j = 1 , 2 a n d i j
It is the same amplitude of graph V 1 (Figure 2) by identifying p 1 k j and setting p 2 = 0 . Therefore, we can write
A i Ξ = A ¯ i Ξ + α i Ξ , i , j = 1 , 2 a n d i j A ¯ i Ξ k j Δ ( k j ) Δ ( k j p ) I log ( λ 2 ) , α i Ξ b 6 k j Δ ( k j ) Δ ( k j p ) 2 ln k j 2 μ 2 λ 2 .
We turn to A 3 Ξ . We choose to perform the integral in k 1 first (which is finite), insert the result in the integral in k 2 and use the rules of IREG to obtain
α ¯ 3 Ξ A 3 Ξ = b 6 p 2 I log ( λ 2 ) 3 .
Similarly,
A 4 Ξ = k 2 Δ 2 ( k 2 ) I log ( λ 2 ) b 6 ln k 2 2 μ 2 λ 2 + 2 b 6 = 0
in the limit μ 2 0 .
We will verify that the terms A ¯ i Ξ ( i = 1 , 2 ) are exactly the ones which are going to be subtracted by applying Bogoliubov’s recursion formula. We set them aside for now and evaluate the rest ( α i Ξ ). As usual, after using identity (19) in the propagator that depends on the external momentum, we are able to identify the divergent terms. After taking the limit μ 2 0 , the only one that survives is
α ¯ i Ξ k j Δ ( k j ) f 2 ( k j , p ) b 6 ln k j 2 μ 2 λ 2 + 2 b 6 = b 6 p 2 I log ( 2 ) ( λ 2 ) 3 8 9 I log ( λ 2 ) .
Hence, the divergent content of Ξ A ( 2 ) amounts to
Ξ A ( 2 ) i g 4 1 2 α ¯ 1 Ξ + α ¯ 2 Ξ + α ¯ 3 Ξ + A ¯ 1 Ξ + A ¯ 2 Ξ .
As stated before, the two last terms are exactly the ones that are going to be subtracted after applying Bogoliubov’s recursion formula. Explicitly, the subdivergences of this particular graph are subtracted by the counterterms shown in Figure 4.
The amplitudes of which are, respectively,
i g 4 2 k 2 Δ ( k 2 ) Δ ( k 2 p ) I log ( λ 2 ) = i g 4 2 A ¯ 1 Ξ ,
i g 4 2 k 1 Δ ( k 1 ) Δ ( k 1 p ) I log ( λ 2 ) = i g 4 2 A ¯ 2 Ξ .
Notice that we are adopting a minimal subtraction scheme which, in the context of IREG, corresponds to the subtraction of basic divergent integrals [32]. Therefore, after the subtraction of subdivergences, we obtain
Ξ ¯ A ( 2 ) i g 4 b 6 p 2 6 [ 2 I log ( 2 ) ( λ 2 ) 13 3 I log ( λ 2 ) + finite ] .
We consider next the two-loop nested graph ( P B ( 2 ) ) with an amplitude of
Ξ B ( 2 ) i g 4 = 1 2 k 1 k 2 Δ 2 ( k 1 ) Δ ( k 1 p ) Δ ( k 2 ) Δ ( k 1 k 2 ) .
Once again, identity (19) is applied in the propagators that depend on the external momentum, allowing us to choose n ( k 1 ) = 3 to assure that terms containing f ¯ ( k 1 , p ) are finite as k 1 in all possible ways. We proceed to classify the divergent terms. It is easy to see that the case k 1 and k 2 fixed does not have any divergent term, while the case k 2 and k 1 fixed does given below
k 1 k 2 Δ 2 ( k 1 ) Δ ( k 1 k 2 ) Δ ( k 2 ) l = 0 4 f l ( k 1 , p ) + f ¯ ( k 1 , p ) = k 1 k 2 Δ 2 ( k 1 ) Δ ( k 1 k 2 ) Δ ( k 2 ) Δ ( k 1 p ) .
For definiteness, we denote the above integral B 1 Ξ . It is the only term that we have to deal with (the divergent terms from the case k 1 and k 2 simultaneously are contained in the above integral). It should be noted that, although it is the original amplitude, we have now a natural order to implement IREG. After a straightforward use of the rules in the integral in k 2 , we have
B 1 Ξ = B ¯ 1 Ξ + β 1 Ξ , B ¯ 1 Ξ k 1 Δ ( k 1 ) Δ ( k 1 p ) I log 3 ( λ 2 ) , β 1 Ξ b 6 3 k 1 Δ ( k 1 ) Δ ( k 1 p ) ln k 1 2 μ 2 λ 2 8 3
By applying the procedure again in β 1 Ξ , we can obtain the following divergent terms
β ¯ 1 Ξ b 6 3 k 1 Δ ( k 1 ) l = 0 2 f l ( k 2 , p ) ln k 1 2 μ 2 λ 2 8 3 ] = b 6 p 2 9 I log ( 2 ) ( λ 2 ) 10 3 I log ( λ 2 ) .
Thus, the divergent content of Ξ B ( 2 ) is given by
Ξ B ( 2 ) i g 4 1 2 β ¯ 1 Ξ + B ¯ 1 Ξ .
As before, the last term is the one to be removed by an application of Bogoliubov’s recursion formula since the counterterm we consider in this case is shown in Figure 5.
After removing the subdivergence, we finally obtain
Ξ ¯ B ( 2 ) i g 4 b 6 p 2 18 I log ( 2 ) ( λ 2 ) + 10 3 I log ( λ 2 ) + f i n i t e .
Finally, we are able to write down the divergent part of the two-point function at two-loop order
Ξ ¯ div ( 2 ) Ξ ¯ A ( 2 ) + Ξ ¯ B ( 2 ) div = i g 4 p 2 6 5 b 6 3 I log ( 2 ) ( λ 2 ) 29 b 6 9 I log ( λ 2 ) .

3.1.2. Two-Loop Vertex Renormalization

We consider now the renormalization of the vertex. The Feynman diagrams to be evaluated are depicted in Figure 6.
The amplitude corresponding to V A ( 2 ) is given by
Λ A ( 2 ) i g 5 k 1 k 2 Δ ( k 1 ) Δ ( k 2 k 1 ) i = 1 2 Δ ( k i p 1 ) Δ ( k i p 2 ) .
Our first task is to obtain n j ( k i ) . We apply the same procedure as before, choosing n 1 ( k 1 ) = n 2 ( k 1 ) = n 1 ( k 2 ) = n 2 ( k 2 ) = 1 which amounts to
k 1 k 2 Δ ( k 1 ) Δ ( k 2 k 1 ) f 0 ( k 1 , p 1 ) + f ¯ ( k 1 , p 1 ) f 0 ( k 1 , p 2 ) + f ¯ ( k 1 , p 2 ) × f 0 ( k 2 , p 1 ) + f ¯ ( k 2 , p 1 ) f 0 ( k 2 , p 2 ) + f ¯ ( k 2 , p 2 ) .
The divergent terms come only from the case k 2 and k 1 fixed, yielding
A 1 Λ k 1 k 2 Δ ( k 1 ) Δ ( k 2 k 1 ) f 0 ( k 2 , p 1 ) f 0 ( k 2 , p 2 ) f 0 ( k 1 , p 1 ) + f ¯ ( k 1 , p 1 ) f 0 ( k 1 , p 2 ) + f ¯ ( k 1 , p 2 ) = k 1 k 2 Δ ( k 1 ) Δ ( k 2 k 1 ) Δ 2 ( k 2 ) Δ ( k 1 p 1 ) Δ ( k 1 p 2 ) = A ¯ 1 Λ + α 1 Λ , A ¯ 1 Λ k 1 Δ ( k 1 ) Δ ( k 1 p 1 ) Δ ( k 1 p 2 ) I log ( λ 2 ) , α 1 Λ b 6 k 1 Δ ( k 1 ) Δ ( k 1 p 1 ) Δ ( k 1 p 2 ) 2 ln k 1 2 μ 2 λ 2 .
By applying the procedure in α 1 Λ , we obtain the divergent term below
α ¯ 1 Λ b 6 k 1 Δ ( k 1 ) f 0 ( k 1 , p 1 ) f 0 ( k 1 , p 2 ) 2 ln k 1 2 μ 2 λ 2 = 2 b 6 I log ( λ 2 ) b 6 I log ( 2 ) ( λ 2 ) .
Hence, the divergent content of Λ A ( 2 ) is given by
Λ A ( 2 ) | div i g 5 α ¯ 1 Λ + A ¯ 1 Λ
where the last term is removed by adding the counterterm shown in Figure 7.
Therefore, after the subtraction of the subdivergence, we have
Λ ¯ A ( 2 ) i g 5 b 6 I log ( 2 ) ( λ 2 ) + 2 I log ( λ 2 ) + finite .
We move to graph V B ( 2 ) , with an amplitude of
Λ B ( 2 ) i g 5 1 2 k 1 k 2 Δ 2 ( k 1 ) Δ ( k 1 p 1 ) Δ ( k 1 p 2 ) Δ ( k 2 ) Δ ( k 2 k 1 ) .
After choosing n 1 ( k 1 ) = n 2 ( k 1 ) = 1 , we notice that the divergent terms are all contained in the case k 2 and k 1 fixed, yielding
B 1 Λ k 1 k 2 Δ 2 ( k 1 ) Δ ( k 2 ) Δ ( k 2 k 1 ) f 0 ( k 1 , p 1 ) + f ¯ ( k 1 , p 1 ) f 0 ( k 1 , p 2 ) + f ¯ ( k 1 , p 2 ) = k 1 k 2 Δ 2 ( k 1 ) Δ ( k 2 ) Δ ( k 2 k 1 ) Δ ( k 1 p 1 ) Δ ( k 1 p 2 ) = B ¯ 1 Λ + β 1 Λ , B ¯ 1 Λ k 1 Δ ( k 1 ) Δ ( k 1 p 1 ) Δ ( k 1 p 2 ) I log 3 ( λ 2 ) , β 1 Λ k 1 Δ ( k 1 ) Δ ( k 1 p 1 ) Δ ( k 1 p 2 ) b 6 3 ln k 1 2 μ 2 λ 2 8 b 6 9 .
Repeating the procedure in β 1 Λ yields
β ¯ 1 Λ k 1 Δ ( k 1 ) i = 1 2 f 0 ( k 1 , p i ) b 6 3 ln k 1 2 μ 2 λ 2 8 b 6 9 = b 6 3 I log ( 2 ) ( λ 2 ) + 8 3 I log ( λ 2 )
Thus, the divergent content of Λ B ( 2 ) is given by
Λ B ( 2 ) | div i g 5 2 ( β ¯ 1 Λ + B ¯ 1 Λ ) .
The last term is going to be removed after applying Bogoliubov’s recursion formula since the counterterm for this graph is given in Figure 8.
After removing the subdivergence, we obtain
Λ ¯ B ( 2 ) i g 5 b 6 6 I log ( 2 ) ( λ 2 ) 4 b 6 9 I log ( λ 2 ) + finite .
Finally, we evaluate graph V C ( 2 ) . Denoting its amplitude as Λ C ( 2 ) we have
Λ C ( 2 ) i g 5 1 2 k 1 k 2 Δ ( k 1 ) Δ ( k 2 ) Δ ( k 1 p 1 ) Δ ( k 2 p 2 ) Δ ( k 1 + k 2 p 1 ) Δ ( k 1 + k 2 p 2 ) .
As usual, we choose n 1 ( k 1 ) = n 2 ( k 2 ) = n 1 ( k 1 + k 2 ) = n 2 ( k 1 + k 2 ) = 1 to find that the only divergent term comes from the case k 1 and k 2 , simultaneously yielding
C 1 Λ k 1 k 2 f 0 ( k 1 , p 1 ) f 0 ( k 2 , p 2 ) f 0 ( k 1 + k 2 , p 1 ) f 0 ( k 1 + k 2 , p 2 ) = k 1 k 2 Δ 2 ( k 1 ) Δ 2 ( k 2 ) Δ 2 ( k 1 + k 2 ) .
After performing the integral over k 2 , the rules of IREG can be applied to yield
Λ ¯ C ( 2 ) Λ C ( 2 ) = i g 5 b 6 I log ( λ 2 ) + f i n i t e .
Grouping all the results, we obtain that the divergent part of the three-point function at two-loop order is given by
Λ ¯ div ( 2 ) ( Λ ¯ A ( 2 ) + Λ ¯ B ( 2 ) + Λ ¯ C ( 2 ) ) div = i g 5 5 b 6 2 I log ( 2 ) ( λ 2 ) 17 b 6 3 I log ( λ 2 ) .

3.1.3. Two-Loop Renormalization Group Functions

In summary, we can obtain the one- and two-loop counterterms to the propagator and vertex function in a minimal subtraction scheme as
Ξ ct = i g 2 6 I log ( λ 2 ) g 4 6 5 b 6 3 I log ( 2 ) ( λ 2 ) 29 b 6 9 I log ( λ 2 ) ;
Λ ct = i g 2 I log ( λ 2 ) g 4 5 b 6 2 I log ( 2 ) ( λ 2 ) 17 b 6 3 I log ( λ 2 ) .
Given the counterterms, it is an easy task to obtain the renormalization group functions. For completeness, we present the usual definitions
ϕ o Z ϕ 1 2 ϕ , g o Z g g , Ξ ct Z ϕ 1 , Λ ct Z g Z ϕ 3 2 1 , γ λ ln Z ϕ λ , β g λ ln Z g λ ,
where λ plays the role of the renormalization group scale in the context of IREG. We finally obtain
γ = g 2 6 ( 4 π ) 3 + 13 g 4 216 ( 4 π ) 6 + O ( g 6 ) ,
β = 3 g 3 4 ( 4 π ) 3 125 g 5 144 ( 4 π ) 6 + O ( g 6 ) .
This result agrees with the one in the literature [33]. A treatment for massive theories can be performed in a similar fashion; for further details, we refer the reader to [23]. Before moving to gauge theories, we emphasize that the algorithm applied in this section can be generalized to arbitrary loop order. We provide in [23] examples at 4- and n-loop order in the context of the ϕ 3 model.

3.2. Gauge Theories

In this section, we will apply our procedure to massless gauge theories up to two-loop order. Since the complexity and number of diagrams is far greater than the example we provided for the ϕ 3 model, we will not present all the details. Our main aim is to discuss the importance of defining a normal form (as stated in Section 2) and collect known results. We will also only be interested in the UV behavior, meaning that all results consider only off-shell external momenta, avoiding the appearance of IR divergences.
We begin with QED at 1-loop level. The divergent diagrams represent the radiative correction to the photon and electron propagators, as well as the vertex diagram. In the context of IREG, they have been computed in [34,35]. We will comment on the photon and electron propagator in more detail, since they illustrate some interesting features, while for the vertex, we will just quote the result.
Starting with the photon propagator in the massless limit, after a standard application of Feynman rules, we obtain the amplitude
i Π μ ν ( p ) = ( ) ( i e ) 2 k T r γ μ i ( k ) γ ν i ( k p ) .
In the context of IREG, our first step is to perform Dirac algebra in order to obtain a normal form compatible with gauge symmetry [3,24]. For pedagogical purposes, we will skip this step and discuss the consequences. Removing the Dirac trace outside the integral, the amplitude can be written as
i Π μ ν ( p ) = ( e 2 ) T r γ μ γ α γ ν γ β ( I α β I α p β ) w h e r e I α 1 α n = k k α 1 k α n k 2 ( k p ) 2 .
It is an easy task to apply the IREG rules to these integrals to obtain
i Π μ ν ( e 2 ) = 4 3 I log ( λ 2 ) b ln p 2 λ 2 + 5 3 b ( g μ ν p 2 p μ p ν ) + 2 b 3 g μ ν
As can be clearly seen, the result is NOT gauge invariant. The offending term comes from dismissing the first step on defining a normal form. Actually, by performing the trace before applying the regularization rules, one would obtain a term such as
k k 2 k 2 ( k p ) 2 = k 1 ( k p ) 2 = lim μ 2 0 k 1 ( k p ) 2 μ 2 = lim μ 2 0 I quad ( μ 2 ) = 0
while by performing the traces afterwards, one obtains
g α β k k α k β k 2 ( k p ) 2 = g α β p α p β 3 g α β p 2 12 I log ( λ 2 ) b ln p 2 λ 2 + 13 b 6 g α β b p 2 24 = b p 2 6
where we have already discarded I quad ( μ 2 ) . Notice that the difference in both results has its origin in the fact that symmetric integration is not allowed in divergent integrals [3] in dimension-specific methods such as IREG. More importantly, the difference among Equations (79) and (80) is exactly the offending term in Equation (78), since the Dirac trace will contribute with a global factor of (−4).
Therefore, after the correct application of the rules of Section 2, one obtains
i Π μ ν = 4 e 2 3 I log ( λ 2 ) b ln p 2 λ 2 + 5 3 b ( g μ ν p 2 p μ p ν )
which is transverse as expected.
We move to the electron self-energy whose amplitude (in Feynman gauge) is given by
i Σ ( p ) = ( i e ) 2 k i k 2 γ ν i ( p k ) γ ν .
The application of IREG to it is straightforward, remembering that one must first perform the Dirac algebra
i Σ ( p ) = 2 e 2 k p k k 2 ( k p ) 2 .
The interesting point to be noted here will be the difference in the finite part in contrast to the dimensional regularization result. The reason boils down to the difference in the definition of the Dirac algebra in the different methods. For instance, in IREG,
γ β p γ β = 2 p ,
while in DReg,
γ β p γ β = ( 2 d ) γ μ = 2 ( 1 ϵ ) p .
In the context of DReg, the evanescent part [36,37] will hit poles in ϵ , generating more finite terms in contrast to the result in IREG. Nevertheless, our result is identical to the one of Dimensional Reduction after reducing the ϵ -scalar coupling to the standard gauge coupling [36]. Thus, we obtain
i Σ ( p ) = ( e ) 2 p I log ( λ 2 ) b ln p 2 λ 2 + 2 b .
Finally, the vertex function can be evaluated in a similar manner, we just quote its divergent part
i Λ μ = e 3 γ μ I log ( λ 2 ) .

3.2.1. One-Loop renormalization for QED and QCD

Proceeding to the renormalization program, the definitions in QED are
ψ 0 = Z 2 ψ , A μ 0 = Z 3 A μ , e 0 = Z e e .
With this, it is straightforward to obtain the renormalization functions
Λ c t = Z 1 1 , Σ c t = Z 2 1 , Π c t = Z 3 1 , Z 1 = Z 2 , Z 3 = 1 + 4 3 i e 2 I log ( λ 2 )
and check that the Ward identities are satisfied.
For QCD, the analysis is analogous. The usual definitions are
A 0 μ a = Z 3 1 / 2 A μ a , c 0 a = Z ˜ 3 1 / 2 c a , ψ 0 = Z 2 1 / 2 ψ , g 0 = Z g g , Z 1 Z g Z 3 1 / 2 , Z 4 Z g 2 Z 3 2 , Z ˜ 1 Z g Z ˜ 3 Z 3 1 / 2 , Z 1 F Z g Z 2 Z 3 1 / 2 .
Given the large number of diagrams, we will just present the renormalization functions at 1-loop order, and refer the reader to [29]
Z 3 = 1 i g 2 5 3 C 2 ( G ) 4 3 n f C ( r ) I log ( λ 2 ) , Z 2 = 1 + i g 2 C 2 ( r ) I log ( λ 2 ) , Z ˜ 3 = 1 i g 2 2 C 2 ( G ) I log ( λ 2 ) , Z 1 F = 1 + i g 2 ( C 2 ( G ) C 2 ( r ) ) I log ( λ 2 ) , Z ˜ 1 = 1 + i g 2 2 C 2 ( G ) I log ( λ 2 ) , Z 1 = 1 + i g 2 2 3 C 2 ( G ) + 4 3 C ( r ) n f I log ( λ 2 ) , Z 4 = 1 + i g 2 1 3 C 2 ( G ) + 4 3 C ( r ) n f I log ( λ 2 ) .
Once again, the Slavnov–Taylor identities can be easily checked
Z 1 Z 3 = Z ˜ 1 Z ˜ 3 = Z 1 F Z 2 = Z 4 Z 1 = 1 + i g 2 C 2 ( G ) I log ( λ 2 )

3.2.2. Two-Loop β Functions for QED and QCD

In the context of this review, we will focus on the calculation of the two-loop coefficient of the β -function gauge coupling both in QED and QCD. Since these coefficients are universal if a renormalization substraction scheme independent of the mass is chosen (which is the case of IREG), we must recover the well-known results [38,39]. Given this objective, we will restrict ourselves to two-point functions only, which amounts to the following topologies depicted in Figure 9.
The reason of this restriction is as follows: in the context of QED, the knowledge of the photon self-energy is enough to determine the gauge coupling β -function. In QCD, if one applies the background field method [40], this statement is also true since Z g = Z A 1 / 2 . Notice that tadpoles diagrams were already removed, since they will amount to scaleless quadratic divergences (encoded as I quad ( l ) ( μ 2 ) ), which vanish in massless theories. Given the general topologies shown in Figure 9, we need to to fill them with the particle content of the theory under study. For instance, in QED, since four-point interactions do not occur, one can realize only topologies T2 and T3. This task was automatically performed by [41]. After building the amplitudes in all theories, we designed in-house routines to perform Dirac and Lorentz algebra in 4-dimensions in the amplitude as a whole, making use of FormCalc [42].
This is crucial to implement a normal form that respects gauge invariance at two-loop level, as we briefly discuss. In order to illustrate this point, consider the diagram of Figure 10 , whose amplitude (in the Feynman gauge) is schematically given by
A k Π μ ν α β ( k , p ) l F α β ( q , k , p ) q 2 ( q k ) 2 ,
where we denote q as the internal momentum of the sub diagram (gluon loop), k is the internal momentum of the complete diagram, and p is the external momentum. As we presented in Section 2, one must first perform Dirac algebra, Lorentz contractions, etc., before applying the rules of IREG. Therefore, if F α β contains a term such as q α q β while Π μ ν α β ( k , p ) has a term such as g α β k μ k ν , one will obtain
B k k μ k ν k 4 ( k p ) 2 q q 2 q 2 ( q k ) 2 = k k μ k ν k 4 ( k p ) 2 q 1 ( q k ) 2 = 0
Notice the similarity between the equation above and Equation (79). On the other hand, if one opts to perform the contraction with g α β afterwards (see Equation (80)),
g α β q q α q β q 2 ( q k ) 2 = b k 2 6 ,
a non-null result to B will be obtained. As in the 1-loop case of QED, the second choice will result in the breaking of gauge invariance.
We proceed to a collection of results found in [24] in which we computed the two-loop correction to Z A , the renormalization function of the external gauge boson (the photon for QED, the background gluon field for QCD). Defining
Z A = 1 + g 2 ( 4 π ) 2 Z A ( 1 ) + g 4 ( 4 π ) 4 Z A ( 2 ) ,
one obtains for QED
Z A ( 1 ) = 4 3 b I log ( λ 2 ) ,
Z A ( 2 ) = 4 b I log ( λ 2 ) ,
and for QCD
Z A ( 1 ) = 11 b 2 3 b n f I log ( λ 2 ) ,
Z A ( 2 ) = 54 b 2 I log 2 ( λ 2 ) 2 b I log ( 2 ) ( λ 2 ) + 210 b 38 3 b n f I log ( λ 2 ) .
We have written the result for the SU(3) group, while keeping the number of fermions as a free parameter n f . As standard, by defining the β -function by
β = g β 0 g 4 π 2 + β 1 g 4 π 4 ;
one obtains
β 0 = 1 2 λ λ Z A ( 1 ) , β 1 = 1 2 λ λ Z A ( 2 ) .
Finally, using the results shown in Equations (97)–(99), we obtain the well-known one- and two-loop contributions for the gauge β coupling in QED and QCD [38,39]
QED : β 0 = 4 3 ; β 1 = 4 ;
QCD : β 0 = 11 2 3 n f ; β 1 = 102 38 3 n f .
Notice that the renormalization group β function has been obtained in a gauge-invariant fashion without explicitly evaluating the UV divergencies. This is in contrast with other schemes that operate in the physical dimension such as differential renormalization [1] in which divergent expressions are replaced by finite. Moreover, in IREG, the basic divergent objects depend explicitly on a renormalization scale λ whereas the back-bones of ultraviolet divergencies in dimensional methods are terms proportional to 1 / ϵ n , n = 1 , 2 , , as ϵ 0 .

4. Concluding Remarks

The purpose of this review was to present in a pedagogical way how the IREG method is implemented to comply with the powerful framework of BPHZ, which is based on the fundamental principles of quantum field theory, unitarity, causality and locality. An algorithm has been shown that delivers the integrals involved in multi-loop amplitudes being decomposed in structures that are identified as the counterterms and divergencies of the order according to the BPHZ scheme. Various examples, ranging from the cubic scalar theory in six space time dimensions, to QED and QCD in the background field method have been worked out, highlighting the procedure. A further benefit of the method is that it automatically delivers all the necessary ingredients to obtain renormalization group functions, of which we have presented the beta functions to two-loop order of the above-mentioned theories, with known universal coefficients.

Author Contributions

Conceptualization, A.C., B.H., M.S; methodology, A.C., B.H., M.S; software, D.C.A.-P., A.C.; validation, D.C.A.-P.; formal analysis, A.C., B.H., M.S; investigation, D.C.A.-P., A.C., B.H., M.S; writing—original draft preparation, D.C.A.-P., A.C., B.H., M.S; writing—review and editing, D.C.A.-P., A.C., B.H., M.S; funding acquisition, A.C., B.H., M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This publication is based upon work from COST Action CA16201 PARTICLEFACE, supported by COST (European Cooperation in Science and Technology, www.cost.eu). M.S. acknowledges a research grant from CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico—303482/2017-6). We acknowledge support from Fundação para a Ciência e Tecnologia (FCT) through the projects UID/FIS/04564/2020 and CERN/FIS-COM/0035/2019, and the networking support by the COST Action CA16201. This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES)—Finance Code 001.

Data Availability Statement

Not Applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Freedman, D.Z.; Johnson, K.; Latorre, J.I. Differential regularisation and renormalisation: A New method of calculation in quantum field theory. Nucl. Phys. B 1992, 371, 353–414. [Google Scholar]
  2. Abi, B.; Albahri, T.; Al-Kilani, S.; Allspach, D.; Alonzi, L.P.; Anastasi, A.; Anisenkov, A.; Azfar, F.; Badgley, K.; Baeßler, S.; et al. Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm. Phys. Rev. Lett. 2021, 126, 141801. [Google Scholar] [CrossRef]
  3. Bruque, A.M.; Cherchiglia, A.L.; Pérez-Victoria, M. Dimensional regularisation vs methods in fixed dimension with and without γ5. J. High Energy Phys. 2018, 1808, 109. [Google Scholar]
  4. Viglioni, A.C.D.; Cherchiglia, A.L.; Vieira, A.R.; Hiller, B.; Sampaio, M. γ5 algebra ambiguities in Feynman amplitudes: Momentum routing invariance and anomalies in D = 4 and D = 2. Phys. Rev. D 2016, 94, 065023. [Google Scholar]
  5. Porto, J.S.; Vieira, A.R.; Cherchiglia, A.L.; Sampaio, M.; Hiller, B. On the Bose symmetry and the left- and right-chiral anomalies. Eur. Phys. J. C 2018, 78, 160. [Google Scholar]
  6. ’t Hooft, G.; Veltman, M.J.G. Regularization and Renormalization of Gauge Fields. Nucl. Phys. B 1972, 44, 189–213. [Google Scholar] [CrossRef] [Green Version]
  7. Bollini, C.G.; Giambiagi, J.J. Lowest order divergent graphs in nu-dimensional space. Phys. Lett. B 1972, 40, 566–568. [Google Scholar] [CrossRef]
  8. Rivat, S.; Grinbaum, A. Philosophical foundations of effective field theories. Eur. Phys. J. A 2020, 56, 90. [Google Scholar] [CrossRef]
  9. Kinoshita, T. Mass singularities of Feynman amplitudes. J. Math. Phys. 1962, 3, 650. [Google Scholar]
  10. Lee, T.D.; Nauenberg, M. Degenerate Systems and Mass Singularities. Phys. Rev. 1964, 133, B1549. [Google Scholar]
  11. Breitenlohner, P.; Maison, D. Dimensional Renormalization and the Action Principle. Commun. Math. Phys. 1977, 52, 11. [Google Scholar]
  12. Breitenlohner, P.; Maison, D. Dimensionally Renormalized Green’s Functions for Theories with Massless Particles—1. Commun. Math. Phys. 1977, 52, 39. [Google Scholar]
  13. Breitenlohner, P.; Maison, D. Dimensionally Renormalized Green’s Functions for Theories with Massless Particles—2. Commun. Math. Phys. 1977, 52, 55. [Google Scholar]
  14. Bonneau, G. Zimmermann Identities And Renormalization Group Equation in Dimensional Renormalization. Nucl. Phys. B 1980, 167, 261. [Google Scholar]
  15. Bogoliubov, N.; Parasiuk, O. On the Multiplication of the causal function in the quantum theory of fields. Acta Math. 1957, 97, 227. [Google Scholar]
  16. Hepp, K. Proof of the Bogolyubov-Parasiuk theorem on renormalisation. Commun. Math. Phys. 1966, 2, 301. [Google Scholar]
  17. Zimmermann, W. Local field equation for A4 coupling in renormalized perturbation theory. Commun. Math. Phys. 1967, 6, 161. [Google Scholar]
  18. Zimmermann, W. Convergence of Bogoliubov’s method of renormalisation in momentum space. Comm. Math. Phys. 1969, 15, 208–234. [Google Scholar]
  19. Piguet, O.; Sorella, S.P. Algebraic Renormalisation: Perturbative Renormalisation, Symmetries and Anomalies; Lecture Notes in Physics Monographs; Springer: Berlin, Germany, 1995; Volume 28, p. 1. [Google Scholar]
  20. Epstein, H.; Glaser, V. The role of locality in perturbation theory. Ann. Inst. Henri Poincaré Sect. A 1973, 19, 211. [Google Scholar]
  21. Herzog, F. Zimmermann’s forest formula, infrared divergences and the QCD beta function. Nucl. Phys. B 2018, 926, 370–380. [Google Scholar] [CrossRef]
  22. Battistel, O.A.; Mota, A.L.; Nemes, M.C. Consistency conditions for 4-D regularizations. Mod. Phys. Lett. A 1998, 13, 1597. [Google Scholar]
  23. Cherchiglia, A.; Sampaio, M.; Nemes, M.C. Systematic Implementation of Implicit regularisation for Multi-Loop Feynman Diagrams. Int. J. Mod. Phys. A 2011, 26, 1. [Google Scholar]
  24. Cherchiglia, A.; Arias-Perdomo, D.C.; Vieira, A.R.; Sampaio, M.; Hiller, B. Two-loop renormalisation of gauge theories in 4D Implicit Regularisation: Transition rules to dimensional methods. arXiv 2020, arXiv:2006.10951. [Google Scholar]
  25. Dias, E.W. Generalização do procedimento de regularização implícita para ordens superiores em teorias de calibre abelianas. Ph.D. Thesis, Federal University of Minas Gerais, Belo Horizonte, Brazil, 2008. [Google Scholar]
  26. Brito, L.C.; Sampaio, M.; Fargnoli, H.; Nemes, M.C. Systematisation of Basic Divergent Integrals in Perturbation Theory and renormalisation Group Functions. Phys. Lett. B 2009, 673, 220. [Google Scholar]
  27. Ferreira, L.; Cherchiglia, A.; Nemes, M.C.; Hiller, B.; Sampaio, M. Momentum routing invariance in Feynman diagrams and quantum symmetry breakings. Phys. Rev. D 2012, 86, 025016. [Google Scholar]
  28. Muta, T. Foundations of QCD; World Scientific: Singapore, 1987. [Google Scholar]
  29. Sampaio, M.; Baêta Scarpelli, A.P.; Ottoni, J.E.; Nemes, M.C. Implicit regularisation and renormalisation of QCD. Int. J. Theor. Phys. 2006, 45, 436. [Google Scholar]
  30. Cherchiglia, A.; Cabral, L.A.; Nemes, M.C.; Sampaio, M. (Un)determined finite regularisation dependent quantum corrections: the Higgs boson decay into two photons and the two photon scattering examples. Phys. Rev. D 2013, 87, 065011. [Google Scholar]
  31. Cherchiglia, A.L.; Vieira, A.R.; Hiller, B.; Scarpelli, A.P.B.; Sampaio, M. , Guises and disguises of quadratic divergences. Ann. Phys. 2014, 351, 751. [Google Scholar]
  32. Sampaio, M.; Scarpelli, A.P.B.; Hiller, B.; Brizola, A.; Nemes, M.C.; Gobira, S. Comparing implicit, differential, dimensional and BPHZ renormalization. Phys. Rev. D 2002, 65, 125023. [Google Scholar]
  33. Macfarlane, A.J.; Woo, G. ϕ3 Theory in Six Dimensions and the Renormalization Group. Nucl. Phys. B 1974, 77, 91–108. [Google Scholar]
  34. Battistel, O.A. Uma estratégia para manipulações e cálculos envolvendo divergências em TQC. Ph.D. Thesis, Federal University of Minas Gerais, Belo Horizonte, Brazil, 2000. [Google Scholar]
  35. Dias, E.W.; Scarpelli, A.P.B.; Sampaio, M.; Nemes, M.C. Implicit regularisation beyond one loop order: Gauge field theories. Eur. Phys. J. C 2008, 55, 667. [Google Scholar]
  36. Gnendiger, C.; Signer, A.; Stöckinger, D.; Broggio, A.; Cherchiglia, A.L.; Driencourt-Mangin, F.; Fazio, A.R.; Hiller, B.; Mastrolia, P.; Peraro, T.; et al. To d, or not to d: Recent developments and comparisons of regularisation schemes. Eur. Phys. J. C 2017, 77, 471. [Google Scholar]
  37. Bobadilla, W.J.T.; Sborlini, G.F.R.; Banerjee, P.; Catani, S.; Cherchiglia, A.L.; Cieri, L.; Dhani, P.K.; Driencourt-Mangin, F.; Engel, T.; Ferrera, G.; et al. May the four be with you: Novel IR-subtraction methods to tackle NNLO calculations. Eur. Phys. J. C 2021, 81, 250. [Google Scholar] [CrossRef]
  38. Vladimirov, A.A. Method for Computing Renormalization Group Functions in Dimensional Renormalization Scheme. Theor. Math. Phys. 1980, 43, 417. [Google Scholar] [CrossRef]
  39. Larin, S.A.; Vermaseren, J.A.M. The Three loop QCD Beta function and anomalous dimensions. Phys. Lett. B 1993, 303, 334–336. [Google Scholar] [CrossRef] [Green Version]
  40. Abbott, L.F. The Background Field Method Beyond One Loop. Nucl. Phys. B 1981, 185. [Google Scholar]
  41. Hahn, T. Generating Feynman diagrams and amplitudes with FeynArts 3. Comput. Phys. Commun. 2001, 140, 418–431. [Google Scholar]
  42. Hahn, T.; Perez-Victoria, M. Automatized one loop calculations in four-dimensions and D-dimensions. Comput. Phys. Commun. 1999, 118, 153–165. [Google Scholar]
Figure 1. Graph P ( 1 ) .
Figure 1. Graph P ( 1 ) .
Symmetry 13 00956 g001
Figure 2. Graph V ( 1 ) .
Figure 2. Graph V ( 1 ) .
Symmetry 13 00956 g002
Figure 3. Graphs P A ( 2 ) and P B ( 2 ) , respectively.
Figure 3. Graphs P A ( 2 ) and P B ( 2 ) , respectively.
Symmetry 13 00956 g003
Figure 4. Counterterms for P A ( 2 ) .
Figure 4. Counterterms for P A ( 2 ) .
Symmetry 13 00956 g004
Figure 5. Counterterm for P B ( 2 ) .
Figure 5. Counterterm for P B ( 2 ) .
Symmetry 13 00956 g005
Figure 6. Graphs V A ( 2 ) , V B ( 2 ) and V C ( 2 ) respectively.
Figure 6. Graphs V A ( 2 ) , V B ( 2 ) and V C ( 2 ) respectively.
Symmetry 13 00956 g006
Figure 7. Counterterm for V A ( 2 ) .
Figure 7. Counterterm for V A ( 2 ) .
Symmetry 13 00956 g007
Figure 8. Counterterm for V B ( 2 ) .
Figure 8. Counterterm for V B ( 2 ) .
Symmetry 13 00956 g008
Figure 9. Two-loop topologies for two-point functions.
Figure 9. Two-loop topologies for two-point functions.
Symmetry 13 00956 g009
Figure 10. Two-loop diagram which contains a gluon loop as sub-diagram.
Figure 10. Two-loop diagram which contains a gluon loop as sub-diagram.
Symmetry 13 00956 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Arias-Perdomo, D.C.; Cherchiglia, A.; Hiller, B.; Sampaio, M. A Brief Review of Implicit Regularization and Its Connection with the BPHZ Theorem. Symmetry 2021, 13, 956. https://doi.org/10.3390/sym13060956

AMA Style

Arias-Perdomo DC, Cherchiglia A, Hiller B, Sampaio M. A Brief Review of Implicit Regularization and Its Connection with the BPHZ Theorem. Symmetry. 2021; 13(6):956. https://doi.org/10.3390/sym13060956

Chicago/Turabian Style

Arias-Perdomo, Dafne Carolina, Adriano Cherchiglia, Brigitte Hiller, and Marcos Sampaio. 2021. "A Brief Review of Implicit Regularization and Its Connection with the BPHZ Theorem" Symmetry 13, no. 6: 956. https://doi.org/10.3390/sym13060956

APA Style

Arias-Perdomo, D. C., Cherchiglia, A., Hiller, B., & Sampaio, M. (2021). A Brief Review of Implicit Regularization and Its Connection with the BPHZ Theorem. Symmetry, 13(6), 956. https://doi.org/10.3390/sym13060956

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop