Regularization-scheme dependence of QCD amplitudes in the massive case

We investigate QCD amplitudes with massive quarks computed in the four-dimensional helicity scheme (FDH) and dimensional reduction at NNLO and describe how they are related to the corresponding amplitudes computed in conventional dimensional regularization. To this end, the scheme dependence of the heavy quark and the velocity-dependent cusp anomalous dimensions is determined using soft-collinear effective theory. The results are checked against explicit computations of massive form factors in FDH at NNLO. Our results complete the description of the scheme dependence of QCD amplitudes at NNLO.


Introduction
The most common procedure to regularize ultraviolet (UV) and infrared (IR) singularities of scattering amplitudes is to apply conventional dimensional regularization (cdr), whereby all relevant quantities are treated as D = 4 − 2ǫ dimensional. In cdr, IR singularities of next-to-next-to leading order (NNLO) scattering amplitudes in massless QCD have a remarkably simple structure [1][2][3][4]. Key ingredients are the cusp anomalous dimension γ cusp and the anomalous dimensions of quarks and gluons, γ q and γ g , respectively.
For practical computations it is sometimes advantageous to apply certain variants of cdr, such as the 't Hooft-Veltman scheme (hv) [5], dimensional reduction (dred) [6] or the four-dimensional helicity scheme (fdh) [7]. This leads to the question how virtual amplitudes computed in these schemes are related to the corresponding amplitudes computed in cdr. In the massless case at NNLO, this question has been answered in ref. [8], where, drawing on earlier results [9][10][11][12][13][14][15][16][17][18], it has been shown that the IR structure of cdr is only modified through changes in the anomalous dimensions. We indicate this regularizationscheme (rs) dependence by the shifts γ cusp → γ RS cusp , γ q → γ RS q and γ g → γ RS g . The explicit expressions of the anomalous dimensions as well as the β functions of the various couplings in the different schemes have been determined at least up to NNLO.

JHEP10(2016)034
In the presence of massive quarks there are additional structures in the IR singularities of QCD amplitudes [19]. Hence, the scheme dependence will also have to be generalized. At NLO the scheme-dependence has been discussed in ref. [20]. The generalization of the scheme dependence at NNLO to QCD amplitudes including massive quarks is the main result of this paper. As we will show, once the scheme-dependent UV renormalization has been carried out, this scheme dependence is contained entirely in two additional anomalous dimensions, the velocity-dependent cusp anomalous dimension γ RS cusp (β) and the anomalous dimension of a heavy quark γ RS Q . In fact, the scheme dependence of γ RS cusp (β) itself is induced solely through the scheme dependence of the cusp anomalous dimension γ RS cusp from the massless case.
With the results presented here it is possible to convert any NNLO QCD amplitude between the four schemes cdr, hv, fdh, and dred. This allows for using whatever scheme is most convenient in the computation of the virtual amplitude and then combine this with the real corrections, typically computed in cdr. In fact, for the generalization to the massive case it is sufficient to consider the difference between the fdh and the hv (or cdr) scheme. If there are no external gluons, fdh is equivalent to dred. Hence, the IR anomalous dimensions are the same, e.g. γ fdh cusp (β) = γ dred cusp (β) and γ fdh Q = γ dred Q . Furthermore, cdr and hv also have the same anomalous dimensions, γ cdr cusp (β) = γ hv cusp (β) and γ cdr Q = γ hv Q . These schemes differ simply in the dimension of the polarization sum of external gluons.
Apart from the four schemes treated in this paper, other possibilities to regularize virtual amplitudes have been considered. The fdh scheme has been adapted to the so-called fdf scheme (four-dimensional formulation) for using unitary-based methods to compute NLO amplitudes [21,22]. There are also proposals to abandon dimensional regularization altogether and perform computations completely in four dimensions in the context of implicit regularization [23][24][25][26], fdr (four-dimensional regularization/renormalization) [27][28][29], and using loop-tree duality to deal with IR singularities at the integrand level [30][31][32].
While this list is by no means exhaustive it shows that despite the impressive technical advances in computing higher-order corrections in cdr there is considerable interest in exploring alternative methods. The results presented here complete the description at NNLO of a first step away from a fully D dimensional treatment of the problem. Apart from allowing to perform computations in fdh and dred, we hope it also helps to understand better the relation between cdr and the different four-dimensional approaches mentioned above. The ultimate goal is, of course, to develop efficient methods to explicitly perform ever more accurate computations.
The paper is organized as follows: in section 2 we briefly review the various schemes, the IR structure of amplitudes and its extension to the massive case. We also discuss the UV renormalization, emphasizing the special features of fdh in the presence of massive quarks. Section 3 is devoted to the computation of γ RS Q and γ RS cusp (β) at NNLO in the fdh scheme. These results are obtained by direct computations using soft-collinear effective theory. In order to obtain an independent test of the scheme dependence of NNLO amplitudes, in section 4 we compare the heavy-quark and heavy-to-light form factors in the fdh and cdr schemes and verify that the results are in agreement with the expected scheme dependence JHEP10(2016)034 obtained from the anomalous dimensions. We also provide a guide on how to actually perform computations in the fdh scheme and show that the modifications compared to cdr are minimal. Finally we present our conclusion in section 5.
2 UV and IR structure of massive QCD

DRED and FDH
As has been shown in a series of papers [8,12,17,18,33], a consistent formulation of the dimensional reduction (dred) and the four-dimensional helicity (fdh) scheme in the framework of massless QCD requires the introduction of three vector spaces. In this work we investigate how this can be extended to the case of massive partons. In doing so we do not consider processes including external vector fields. The names fdh and dred are in the following therefore used synonymously, meaning that whenever a statement about the fdh schemes is made, the same argument also applies in dred. For a detailed discussion and a precise definition of the schemes, of the related vector spaces and their algebraic relations we refer to ref. [14]. Here we only provide the most important characteristics.
In fdh, the underlying quasi 4-dimensional space Q4S with metric g µν is split into a direct sum of the quasi D-dimensional space of cdr with metricĝ µν and a disjoint space Q2ǫS with metricg µν : g µν =ĝ µν +g µν . (2.1) In order to have full control over the contributions originating from Q2ǫS, we define complete contractions of the corresponding metric tensors as As a consequence, arbitrary fdh quantities in general depend on N ǫ . They are in the following denoted by a bar. At the level of the Lagrangian, the structure of the different vector spaces is reflected in a split of the quasi 4-dimensional gluon field into a D-dimensional gluon field and an ǫ-scalar field: A µ =Â µ +Ã µ . The 'particles' associated with these fields are in the following denoted by g andg, respectively. In refs. [15,34,35], it has been shown that because of this split in principle five different couplings need to be distinguished in the bare theory: the gauge coupling α s = g 2 s /(4π), thegqq coupling α e = g 2 e /(4π), and three different quartic g-couplings. However, for the calculations presented in this work it is sufficient to consider only α s and α e .
For later purposes it turns out to be useful to include repeatedly occurring universal factors in the definition of the bare couplings where γ E is the Euler-Mascheroni constant, m is the mass of a heavy fermion, and a i ∈ {a s , a e }. As renormalization prescription for the couplings we use the MS scheme JHEP10(2016)034 throughout this work. The constantsZ α i in fdh are given in e.g. ref. [17]. The perturbative expansion of fdh/dred quantities in terms of the UV renormalized couplings is in the following written as

IR factorization at NNLO in the FDH scheme
In cdr, the IR divergence structure of scattering amplitudes including massive external partons has been investigated up to the two-loop level in ref. [19]. Using a combination of soft-collinear effective theory (SCET) (for an introduction see e.g. ref. [36]) and heavyquark effective theory (HQET) (for an introduction see e.g. ref. [37]) it has been shown that amplitudes with an arbitrary number of massive and massless legs factorize into a hard and a soft function, where the latter depends on both massive and massless Wilson lines. For amplitudes including massive partons, the corresponding IR anomalous dimension has less constraints compared to the massless case and additional color structures arise. Starting from the cdr expression for the IR anomalous dimension, we write the twoparton correlation terms of the respective quantity in fdh as where the capital indices I, J correspond to the massive partons and the angle β IJ is defined as For the definition of the color generators T i , of the kinematic variable s ij , and of the sets {p}, {m} we refer to [19]. In eq. (2.5), the first line corresponds to contributions from massless partons, already discussed in refs. [8,17,18]; the remainder is given by additional terms arising in the massive theory. Suppressing the dependence on the couplings, the complete set of IR anomalous dimensions in fdh/dred is given bȳ whereγg only appears in dred. The quantities in the first line have been computed up to the two-loop level in refs. [8,17,18]; the values ofγ cusp (β IJ ) andγ Q are so far unknown and will be given in section 3. Since there is no difference between the IR anomalous dimensions appearing both in fdh and dred, relation (2.5) also holds in dred.

JHEP10(2016)034
In the massive theory, the IR anomalous dimension also contains three-parton correlation terms which we write in fdh as including the four-velocities of the massive partons v µ I := In refs. [38,39], the functions F 1 and f 2 are given for the case of cdr. Since in fdh these functions do not receive evanescent contributions from the ǫ-scalar up to NNLO, eq. (2.8) is a scheme-independent quantity at this order. Its value in fdh is therefore the same as in cdr.
In analogy to the massless case [8,17], we subtract all IR divergences of QCD loop amplitudes by means of a factorZ which is given by a path-ordered integral overΓ (compare with eqs. (2.8) and (2.12) of ref. [8]). This renormalization factor is given in the effective theory where the heavy quarks have been integrated out. Hence, it is written in terms of α i , the couplings defined in the massless theory. In the massive case, however, we also need to take into account contributions from heavy-quark loops. To reproduce the correct IR behavior of the effective low-energy theory we therefore have to perform a matching of the couplings between the full and the effective theory. For an amplitude describing a process with n external partons then the following relation holds: (2.10) As mentioned above, α i is a coupling in the effective theory, meaning that the heavy quark flavors have been integrated out. It is related to the corresponding coupling of the full theory via the decoupling relation α f i = ζ α i α i . Explicit results for the decoupling constants in the fdh scheme will be given in section 2.4.

Mass renormalization of the ǫ-scalar
In the case of massive fermions there is no symmetry that protects the propagator of the ǫ-scalar from acquiring a mass term ∝ m 2gµν where m is a fermion mass. As a consequence, the ǫ-scalar mass is effectively shifted away from zero, even if the ǫ-scalar is massless at the tree-level. Therefore we have to introduce a mass counterterm δm 2 ǫ in the Lagrangian to restore the initial 'on-shell' condition of a vanishing ǫ-scalar mass [40]. At the one-loop level there is only one diagram that effectively generates a mass term in the ǫ-scalar propagator, see figure 1. To obtain the mass counterterm we need to compute the full one-particle irreducible (1PI) two-point function of the ǫ-scalar, whose JHEP10(2016)034 tensor structure is given by including the dimensionless quantities A and B. The mass counterterm can be extracted by writing the propagator of the ǫ-scalar as (2.12) In order to maintain the ǫ-scalar massless we then require 13) where N H denotes the number of heavy quark flavors and the coupling is defined in eq. (2.3). As a consequence, any time we encounter a massive loop diagram insertion as in figure 1, we add the mass counterterm (2.13) in order to impose the on-shell condition of a massless ǫ-scalar.

Decoupling transformations
The decoupling transformation needed in eq. (2.10) is well known for the gauge coupling.
In order to extend it to α e we apply the procedure described in ref. [41] and build an effective Lagrangian in which the heavy quark flavors have been integrated out. As a consequence, the parameters and fields of the effective theory are in general different from the ones of the full theory. To relate the two theories we introduce decoupling constants in the following way: where g and X stand for parameters and fields of the theory, respectively. In this way we are able to relate the full and the effective bare QCD Lagrangian in terms of the re-scaled parameters and fields

JHEP10(2016)034
The decoupling constants can be obtained from a matching calculation. For ζ 0 A which is related to the gluon field decoupling, for example, we get −ĝ µν whereΠ 0 only contains light degrees of freedom andΠ 0,f receives virtual contributions from the heavy quarks. From eqs. (2.16) we then get Since the l.h.s. does not depend on the kinematics of the process it is possible to consider the special case p = 0. The renormalization of the decoupling constant is done in the usual way by means of the gluon field renormalization constants in the effective and the full theory: ζÂ =ZÂ/Z f A ζ 0 A . The same method also applies to the decoupling of the ǫ-scalar field where, however, according to the discussion in section 2.3 a mass counterterm has to be added in order to maintain the ǫ-scalar massless. In fact, this counterterm is even required to ensure that For the calculations in this work we need the decoupling transformations for α s and α e at the one-loop level which can be obtained from a matching of the gqq andgqq vertices, in analogy to eqs. (2.16) , the (bare) one-loop decoupling constants for g s and g e are entirely given by ζ 0 A and ζ 0 A , respectively. Using (ζ 0 gs ) 2 = ζ 0 αs and ζ αs =Z αs /Z f αs ζ 0 αs and similar for the evanescent coupling we finally obtain for the renormalized decoupling constants of α s and α e . JHEP10(2016)034

Field and mass renormalization of the heavy quarks
To obtain UV-finite Green functions in the fdh scheme we need to perform a renormalization of the heavy quark field and mass, where the corresponding renormalization constants are defined by Extending the standard cdr procedure for obtaining renormalization constants in the onshell (OS) scheme, we write the 1PI self-energy of the heavy quark in fdh as The renormalization constants are then given by To obtain their values we calculated the quantitiesΣ 1 andΣ 2 up to the two-loop level, with sample diagrams shown in figure 2. One point of major importance is that apart from genuine two-loop diagrams we have to include contributions originating from UV (sub)renormalization. This in particular comprises the mass counterterm for the ǫ-scalar given in eq. (2.13), see the r.h.s. of figure 2. In terms of the bare couplings we then get . For later purposes it is convenient to introduce a mass counterterm δm = m − m 0 = m (1 −Z m ) for the heavy quarks. Using eq. (2.23b), a direct calculation ofΣ 1 yields up to the two-loop level. The pure α s terms for N ǫ = 0 correspond to the cdr result.

Field renormalization of the light quarks
In analogy to the previous section we determine the field renormalization of the light quark fields where the corresponding renormalization constant is in the following denoted byZ 2,l . As in the case of heavy quarks,Z 2,l receives contributions from heavy quark loops, see  there is no contribution from the ǫ-scalar mass counterterm. The explicit calculation then yields for the field renormalization of the light quark in the fdh schemē As for the mass counterterm, the pure α s terms are of course not new.

IR anomalous dimensions in the massive case
The aim of this section is to provide all so far unknown IR anomalous dimensions present in the general IR factorization formula (2.5), i.e.γ Q andγ cusp (β). As in the massless case [8], for this we use the SCET framework.

Scheme dependence of the heavy-to-light soft function and γ Q
In ref. [42], it has been shown that the top quark decay factorizes into regions where only soft radiation and (or) radiation collinear to the massless partons are present. More precisely, the factorization consists of a hard function whose renormalization group equation (RGE) depends on the heavy-quark anomalous dimension, a quark jet function, and a soft function. In cdr, the jet and soft functions have been calculated up to the two-loop level in refs. [43] and [44], respectively. In fdh, so far only the jet function is known [8].
The general relation between the corresponding IR anomalous dimensions is given by where γ RS S and γ RS J are the (rs-dependent) anomalous dimensions of the soft and jet function. Eq. (3.1) is a direct consequence of the fact that the RGE of the factorization formula does not depend on the factorization scale. The values ofγ J = γ fdh/dred J andγ q = γ fdh/dred q have been calculated in ref. [8] up to the two-loop level. In order to obtainγ Q = γ fdh/dred Q we therefore have to computeγ S = γ fdh/dred S . Extending the approach of ref. [44], we define the scheme-dependent (bare) soft function as where h v are effective quark fields in HQET (see e.g. ref. [37]), b v are on-shell b-quark states with velocity v, and n is a light-like 4-vector with n · v = 1 and n 2 = 0. The normalization is fixed by b v |h v h v |b v = 1.
For explicit calculations it is useful to express the soft function as a contour integral 3) and to work in Laplace space Since h v and b v are Heisenberg fields, the usual perturbative expansion results in loop diagrams contributing to the heavy quark propagator. As in the massless case, the scheme dependence is related to the UV singularities of such diagrams. At the one-loop level there are no evanescent contributions since the ǫ-scalar does not couple to heavy quark lines, see also ref. [8]. There are exactly three diagrams that induce a scheme dependence of the soft function at the two-loop level. They are shown in figure 4. For the explicit computation we generated the diagrams with QGRAF [45] and applied a tensor reduction of the integrals with Reduze 2 [46], where the master integrals needed in fdh are identical to the ones of cdr given in ref. [44].
In fdh we then get up to the two-loop level Taking the limit N ǫ → 0 in eq. (3.5) we obtain the cdr result which is in agreement with the one given in ref. [44].

JHEP10(2016)034
As for the quark and gluon jet functions [8], all divergences of the soft function can be removed multiplicatively by means of a Z factor To relate Z RS S (Ω, µ) with γ RS S we compare the RGE of the soft function, with the RGE written in terms of γ RS S , where L Ω = ln(Ω 2 /µ 2 ) and the cusp anomalous dimension is known from the massless case [8,17,18]. In fdh, the factorZ S is given by and the coefficients of the β function can be found e.g. in ref. [8]. Imposing minimal subtraction with N ǫ as an independent quantity we can read off the soft anomalous dimension which is scheme independent at the one-loop level. Apart fromγ S it is also possible to extract the already known values of the cusp anomalous dimension as well as the β functions in the fdh scheme, which provides a strong consistency check on the applied procedure. Using the obtained results together with eq. (3.1) we then find for the IR anomalous dimension of the heavy quarks in the fdh scheme. Likeγ S , at NLO it does not depend on N ǫ and is therefore scheme independent, as already found in ref. [20]. However, at NNLO it receives rs-dependent contributions.

JHEP10(2016)034
Coupling of a gluon (left) and an ǫ-scalar (right) to a heavy quark propagator. In the eikonal approximation the latter vanishes.
Eq. (3.12) is the main result of this section. However, for the sake of completeness we give the result of the finite and scheme-independent soft function by setting N ǫ = 2ǫ and taking the subsequent limit ǫ → 0 This result is in agreement with the one given in ref. [44].

Determination ofγ cusp (β)
The velocity-dependent cusp anomalous dimensions can be extracted from the heavy-toheavy soft anomalous dimension Γ hh for the pair production of massive quarks. Using cdr, Γ hh has been calculated in ref. [47] in the framework of the eikonal approximation. This method can also be used to derive the respective quantity in fdh.
In general, the eikonal approximation is suited for describing the emission of soft gluons from partons in a hard scattering process, see the l.h.s. of figure 5. For a vanishing gluon momentum, the Feynman rule for the coupling of a gluon to a massive quark propagator can be reduced tō where in the second line the Dirac equationū(p I )( / p I − m I ) = 0 has been used. Since the Feynman rule (3.15c) does not contain a Dirac matrix anymore, the evaluation of loop contributions is much simpler compared to ordinary QCD.
Extending this to the case of an ǫ-scalar we get Due to the vanishing anticommutator, a direct coupling of ǫ-scalars to massive quark propagators does not exist in the eikonal approximation. Following the approach of ref. [47], the soft anomalous dimension for heavy-quark pair production can be obtained from the UV poles of corresponding eikonal diagrams with oneand two-loop examples shown in figure 6. Since there is no direct coupling of ǫ-scalars to massive quarks, the soft anomalous dimension is scheme independent at the one-loop level. At the two-loop level, however, closed ǫ-scalar loops yield evanescent contributions ∝ α s N ǫ .
In the following, the scalar product of the two outgoing velocity vectors is fixed by v I · v J := − cosh β IJ with β IJ given in eq. (2.6), and the indices of β are suppressed: β IJ =: β. Generalizing eq. (14) of ref. [19] to the case of fdh, the result of the soft anomalous dimension can then be written as (3.17) Using eq. (3.12), it is now possible to extract the velocity-dependent cusp anomalous dimension in fdh which in terms of the renormalized couplings reads Since the terms in the curly brackets do not depend on N ǫ , the scheme dependence of γ cusp (β) is entirely governed by the scheme dependence of the cusp anomalous dimension in the massless case, i.e.γ cusp .

JHEP10(2016)034 4 Guideline for FDH calculations and checks of the results
In order to check the obtained results for the scheme dependence of IR divergences in massive QCD we compute the heavy and the heavy-to-light quark form factor in fdh up to the two-loop level. Apart from a pure check this section is also intended to provide a guideline how practical calculations in the fdh scheme can actually be done. For the two-loop calculations we therefore use the following approach: • At the one-loop level we distinguish the ǫ-scalar from the D-dimensional gluon since the related couplings α s and α e renormalize differently.
• At the two-loop level we use a (quasi) 4-dimensional Lorentz algebra for the evaluation of genuine two-loop diagrams and do not distinguish the ǫ-scalar from the D-dimensional gluon.
• After having applied the UV renormalization we set equal the couplings α s and α e in contributions from one-loop counterterm diagrams.
• Throughout the calculations we identify N ǫ = 2ǫ.
Using this setup it turns out that practical calculations in the fdh scheme are not significantly more complicated than the respective ones in cdr.

Heavy quark form factor
In cdr, the heavy quark form factor has been calculated up to NNLO in ref. [48]. In fdh, the Green function for the interaction of a virtual photon and two massive quarks can be written asV Here and in the following, p 1 and p 2 denote the (outgoing) momenta of the two external quarks with p 2 1 = p 2 2 = m 2 and s = (p 1 + p 2 ) 2 /m 2 . In general, the γ matrices appearing in eq. (4.1b) are scheme-dependent. However, since we are only interested in the structure of F 1 their dimensionality can be chosen arbitrarily. Here and in the following we therefore use D-dimensional γ matrices in the Lorentz decomposition.
The IR anomalous dimensions can be obtained from the heavy-quark form factor,F 1 , which can be extracted from eq. (4.1b) by using an appropriate projection operator. For the proper definition of the projection and other details we refer to ref. [48]. In the fdh scheme, only two diagrams contribute to the form factor at the one-loop level, see figure 7. Using 1-dimensional harmonic polylogarithms [49,50] of the variable JHEP10(2016)034 p 1 p 2 Figure 7. One-loop diagrams contributing to the heavy-quark form factor in fdh.
and notation (2.3) for the couplings, we represent the bare one-loop coefficients of the form factor as 1F To obtain the result at the two-loop level we evaluate the Feynman diagrams (see figure 8) using a quasi 4-dimensional Lorentz algebra. This in particular means that the absolute number of diagrams and master integrals [51,52] we have to evaluate is exactly the same as in cdr. In line with that we do not have to introduce evanescent couplings like α e in the computation of the genuine two-loop diagrams. In the following we give the difference between the UV renormalized form factors in fdh and cdr at the two-loop level. For the renormalization of the couplings, the quark mass, and the fields we use eqs.  results we finally get This difference can be expressed in terms of the IR anomalous dimensions and β functions through eqs. (2.5) and (2.10), in a similar way as shown in ref. [8] for the case of massless partons: where F diff 1 =F ren 10 +F ren 01 − F ren 1 is the difference of the UV renormalized one-loop coefficients, i.e. including a field renormalization of the heavy quarks. The fact that the scheme dependence of the IR divergences related to the heavy-quark form factor can be predicted with the results from sections 2 and 3 constitutes a strong consistency check of the results obtained so far.

Heavy-to-light form factor
The cdr result for the decay process b → u W * → u lν l has been computed at NNLO in refs. [53][54][55][56]. Applying the procedure of the previous section we here extend the calculation to the case of fdh, with sample two-loop diagrams shown in figure 9.
In fdh, the tensor structure of the heavy-to-light form factor can be written as with q = p 1 + p 2 . Again, we are interested in the form factorF 1 which can be extracted by means of a projection operator. Accordingly, the matrixγ µ is treated in D dimensions.
JHEP10(2016)034 Figure 9. Sample two-loop diagrams contributing to the heavy-to-light form factor in fdh.
We compute the bare diagrams up to NNLO and perform the UV renormalization exactly in the same way as described in the previous section, taking into account that here only one leg is massive. Again we have to add a counterterm to subtract the ǫ-scalar mass shift. Using eq. As in the previous section we give the difference between the UV renormalized form factors in fdh and cdr up to the two-loop level:  where L is defined as L = ln µ 2 m 2 . In terms of the IR anomalous dimensions, the β functions, and the factorZ defined in eq. (2.10) this difference is given bȳ  The fact that eq. (4.10) matches with eq. (4.11) constitutes an additional and independent check of our results for the IR anomalous dimensions.

Conclusions
The scheme dependence of massless QCD amplitudes at NNLO had been discussed in ref. [8]. In this paper we complete this study by extending it to amplitudes containing massive quarks. This requires modifications in the UV and IR sector. For the UV part, the presence of heavy quarks modifies the renormalization. In particular, the ǫ-scalar field requires a mass counterterm. Also, the decoupling of α e (the coupling of the ǫ-scalars to the quarks) has to be determined. Furthermore, we have computed the additional contributions required in fdh in the quark mass and the quark wave-function renormalization.
Regarding the IR part, the important result is that the IR structure of massive QCD amplitudes in fdh (and dred) is the same as in cdr (and hv). The only change is in the explicit scheme-dependent expressions of the various anomalous dimensions. In the massive case, there are two additional anomalous dimensions, the velocity-dependent cusp anomalous dimension and the heavy-quark anomalous dimension. We have computed them in the fdh scheme, using a SCET approach.

JHEP10(2016)034
We have checked our results by computing the heavy-quark and heavy-to-light form factor in fdh at NNLO. These results differ from the corresponding expressions in cdr. After UV renormalization, the difference can be reproduced by the scheme dependence of the IR factorization formula. This provides us with a strong consistency check and establishes fdh as a consistent regularization scheme also in the massive case, at least to NNLO.