QCD corrections to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H \rightarrow \textit{gg}$$\end{document}H→gg in FDR

I apply FDR—a recently introduced four-dimensional approach to quantum field theories (QFTs)—to the computation of the NLO QCD corrections to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H \rightarrow gg$$\end{document}H→gg in the large top mass limit. The calculation involves all key ingredients of QCD—namely ultraviolet, infrared, and collinear divergences, besides \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _S$$\end{document}αS renormalization—and paves the way for successful use of FDR in massless one-loop QFT computations. I show in detail how the correct result emerges in FDR, and discuss the translation rules to dimensional regularization.


Introduction
Many of the difficulties of higher-order calculations in QFT can be traced back to the treatment, in the framework of dimensional regularization (DR) [1], of the infinities arising in the intermediate steps of the computation. Ultraviolet (UV), infrared (IR), and collinear (CL) divergences are first dimensionally regulated, and then renormalized awayin the UV case-or canceled by combining virtual and real contributions, or reabsorbed in the collinear behavior of the initial state parton densities. In order to attack the problem numerically, it is often necessary to subtract and add back approximations of the IR/CL singular structures. At one loop, several well-tested subtraction procedures have been introduced in the last two decades [2][3][4][5][6][7][8]. At two loops and beyond, the situation is more involved, but progress is under way [9][10][11][12][13][14].
The first obvious ingredient, which may lead to a significant simplification in the above picture, is a computational procedure in which all parts of the calculation can be directly treated in four dimensions. As for the virtual contribution, the FDR approach has been recently introduced in reference [15], which allows a subtraction of the UV divergences a e-mail: pittau@ugr.es at the level of the integrand, leaving a four-dimensional integration over the loop momenta. In the same work, the use of FDR as an IR regulator in QED with massive fermions is also suggested.
In this paper, I present the first application of the FDR ideas in the context of fully massless QCD, where the issues related to gauge invariance are much more subtle than in the QED case. I concentrate, in particular, on the calculation of the O(α S ) gluonic corrections to the H → gg decay in the m top → ∞ limit, and re-derive the well-known fully inclusive result [16,17] where is the lowest order contribution, with N F = 0 in α S (M 2 H ), since only gluons are considered.
Despite its simplicity, all key ingredients of massless QCD are present in this process, such as the simultaneous occurrence of IR/CL divergences and UV renormalization. The fact that the correct expression is reproduced shows that FDR is a valid and consistent approach in massless QFTs, and it gives confidence in its potential to simplify multi-leg/loop computations.
The outline of the paper is as follows. Section 2 provides the set-up of the calculation. In Sect. 3, I review the FDR treatment of the UV divergences and discuss its interplay with the IR and CL infinities. Section 4 presents the FDR computation of the virtual part, while Sect. 5 deals with the real contribution and its merging with the oneloop piece. The connection between FDR and DR is discussed in Sect. 6 and the final conclusions are drawn in Sect. 7.
The gray blobs in V 6 and V 7 represent gluon wave-function corrections and the dashed line stands for the Higgs field. R 1 ( p i , p j , p k ) corresponds to three diagrams with permuted gluons

The model for H → gg
The effective interaction of one Higgs field H with two, three and four gluons-mediated by an infinitely heavy top loopis described by the Lagrangian [18,19] where and v is the vacuum expectation value, v 2 = (G F √ 2) −1 . The corresponding Feynman rules are given in [20], and the diagrams for the decay rate Γ (H → gg) are drawn in Fig. 1.
There are five graphs contributing to the virtual part Γ Vwithout counting gluon wave-function corrections-and four diagrams for the real radiation Γ R . In the following, I separately compute, in FDR, the two pieces, showing how IR/CL divergences drop in the sum

FDR versus infinities
The FDR subtraction of UV infinities is better illustrated with an explicit example. Consider the one-loop quadratically divergent rank-two tensor with Its UV convergence can be improved by first deforming the propagators by a vanishing amount μ 21 and then by repeatedly using the identity wherē Note that the prescription in Eq. (8) avoids possible infrared divergences in the r.h.s. of Eq. (9). The integrand in Eq. (6) can then be rewritten as where the terms in square brackets are UV divergent but depend only on μ 2 . The FDR definition of the integral in Eq. (6) is obtained by integrating the expansion in Eq. (11), after dropping the divergent pieces, and taking the physical limit μ → 0: The r.h.s. of Eq. (12) corresponds to a well-defined fourdimensional integral, in which all UV divergences are explicitly subtracted. Furthermore, IR and CL divergences get also regulated by the propagator deformation. The gaugeinvariance properties of this definition are discussed in detail in [15,21]. In the rest of this section, I mostly concentrate on CL and IR infinities, and, in particular, on the matching between virtual and real contributions. A convenient starting point to study the CL singularities is the fully massless limit of Eq. (12): which vanishes, after tensor decomposition, since p 2 1 = 0. Analogously, one proves that Those results coincide with DR-in which scale-less integrals are zero [22]-and are due to a cancelation between two ln(μ 2 ) of CL and UV origin, respectively. For example, where the first logarithm develops a CL singularity in the limit p 2 → 0. Thus, the virtual CL infinities, generated by 1 → 2 splittings of massless particles, are naturally regulated by the μ 2 -deformed propagators inside the loop, while the external momenta remain massless, as illustrated in Fig. 2 a and c. The real counterpart of this procedure is exemplified in Fig. 2b and d, and corresponds to a phase space in which all would be massless external particles are given a common mass μ and the internal ones stay massless. In other words, one has to replace 2 (16) in any possible singular denominator of the real matrix element squared, integrate over the aforementioned massive phase space, and take the limit μ → 0. As for the IR divergences, the reasoning follows the same lines. For example, the only IR/CL divergent scalar one-loop three-point function generated by the cut in Fig. 2e is 3 with By denoting one computes which is indeed fully matched by the inclusive real contribution in Fig. 2f, as will be shown in Sect. 6. In the following, I use the described approach to UV/CL/IR infinities to compute Γ V (H → gg) and Γ R (H → ggg).

The virtual part Γ V (H → gg)
The calculation is greatly simplified by Eqs. (13) and (14). In fact, only the diagrams V 1 and V 2 in Fig. 1 contribute-as in DR-and the gluon wave function corrections vanish. One computes This simple expression is obtained after a standard Passarino-Veltman [23] decomposition, the only subtlety being the FDR treatment of μ 2 [15,21]: for consistency with Eq. (8), 2 See Fig. 2b. 3 The FDR integration corresponds to a normal integration, in this case, because C(s) is UV finite. a q 2 appearing in the numerator of a diagram should also be deformed, and integrals involving μ 2 , such as require the same integrand expansion as if μ 2 = q 2 . For example, from Eq. (12), The final result follows by inserting Eq. (20) into (21),

The real radiation Γ R (H → ggg) and the fully inclusive result
The unpolarized matrix element squared, derived from the real emission diagrams in Fig. 1 where s i j = ( p i + p j ) 2 . This expression is obtained from the massless result with the replacement 2( p i · p j ) → s i j , in accordance with Eq. (16). As described in Sect. 3, in order to match the virtual IR/CL singularities, |M| 2 should be integrated over a massive three-gluon phase space with p 2 i = μ 2 , which can be parametrized as where √ s is the Higgs mass. It is convenient to introduce the dimensionless variables with μ 0 given in Eq. (19), in terms of which, using the condition, all IR/CL divergent bremsstrahlung integrals can be reduced to the following ones: where the integration region reads, in the fully inclusive case, with Thus Finally, one computes, up to terms which vanish in the limit μ 0 → 0, and so that Adding this to Eq. (25), and accounting for the finite renormalization term in Eq. (4), one obtains All CL/IR ln(μ 2 ) and ln 2 (μ 2 ) cancel in Eq. (37), so that the remaining μ is directly interpreted as the renormalization scale. This is a typical procedure in FDR: since the UV infinities are subtracted from the very beginning, the unphysical left-over μ dependence is eliminated, on the perturbative level one is working at, by a finite renormalization, which fixes the bare parameters in terms of the observables [24]. This is obtained, in the case at hand, by simply replacing 4 in Eq. (37). Then the logarithm is reabsorbed in the gluonic running of the strong coupling constant, and Eq. (1) follows.

FDR versus DR
In this section, I discuss the transition rules between FDR and DR. This is particularly important in QCD, where NLO calculations have to be matched with the running of α S and parton densities, conventionally derived in DR. I consider UV, CL, and IR divergences in turn, showing the equivalence of FDR with the dimensional reduction [25] version of DR, widely used in supersymmetric theories. I start by establishing the connection between the 1/ DR regulator and the ln(μ 2 ) appearing in FDR. As for the UV infinities, it is sufficient to compare the FDR and DR variants of any divergent integral. For instance, the DR counterpart of Eq. (15) (with p 2 = 0) reads where Thus, DR and FDR UV regulators are linked through the simple MS replacement → ln(μ 2 ).
CL virtual singularities follow the same pattern, as can be inferred from the exact UV/CL cancelation in Eqs. (13) and (14). As a consistency check, the DR version of J p reads which indeed coincides with Eq. (35) if = ln(μ 2 ). Finally, the ln 2 (μ 2 ) terms-generated by overlapping IR/CL singularities-drop, together with the full constant part, when adding virtual and fully inclusive real contributions, which can be traced back to the following relation: between Eqs. (20) and (34). An easy calculation shows that the same happens in DR. In fact where p 2 1 = p 2 2 = 0 and s = −2( p 1 · p 2 ), and In summary, Eq. (41) is the only relation needed between the two regulators. However, an important difference between DR and FDR follows from self-contractions of metric tensors coming from the Feynman rules. In DR g αβ g αβ = n, while g αβ g αβ = 4 in FDR. This, together with Eq. (41), and the FDR treatment of μ 2 discussed in Sect. 4, makes explicit the equivalence between FDR and dimensional reduction in the MS scheme. Having established this, all the well-known transition rules between dimensional reduction and DR [26,27] can be directly applied to FDR. In the case of Eq. (37), it turns out that the expression is the same in both dimensional reduction (or FDR) and DR. Therefore, the correct strong coupling constant to be used is the customary α S (μ 2 ) in the MS scheme, proving that the FDR result coincides with Eq. (1).

Conclusions
I have presented an FDR calculation of the gluonic QCD corrections to H → gg in the large top effective theory, demonstrating that ultraviolet, collinear, and infrared divergences can be simultaneously and successfully regulated in four dimensions. I have proved the equivalence, at the oneloop level, of dimensional reduction and FDR, making the latter approach attractive also in supersymmetric calculations, where the fermionic and bosonic sectors must share the same number of degrees of freedom.
The advantage of directly working in the four-dimensional Minkowsky space is expected to lead to considerable simplifications in higher-order QFT computations, especially in connection with numerical techniques. This issue, together with the extension of FDR to more loops, is currently under study.