On the Casimir scaling violation in the cusp anomalous dimension at small angle

We compute the four-loop nf contribution proportional to the quartic Casimir of the QCD cusp anomalous dimension as an expansion for small cusp angle ϕ. This piece is gauge invariant, violates Casimir scaling, and first appears at four loops. It requires the evaluation of genuine non-planar four-loop Feynman integrals. We present results up to Oϕ4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{O}\left({\phi}^4\right) $$\end{document}. One motivation for our calculation is to probe a recent conjecture on the all-order structure of the cusp anomalous dimension. As a byproduct we obtain the four-loop HQET wave function anomalous dimension for this color structure.


Introduction
The cusp anomalous dimension Γ cusp (φ), defined as the anomalous dimension of a Wilson loop with a cusp of angle φ [1], determines the renormalization group evolution of the Isgur-Wise function [2]. In this paper we will mostly be interested in the small φ expansion of Γ cusp (φ). Such an expansion is performed for extracting the CKM matrix element V cb from the B → D * semileptonic decays, see e.g. ref. [3]: the extrapolation of experimental points to φ = 0 is done using the slope and curvature of the Isgur-Wise function, i.e. its φ 2 and φ 4 terms. The cusp anomalous dimension at small angle is also related to real gluon radiation in the case when a heavy quark slightly changes its velocity. This kind of radiation has been considered in refs. [4,5].
The QCD cusp anomalous dimension is currently known to three loops [1,[6][7][8][9] for arbitrary angle φ. Up to this order, the anomalous dimensions for Wilson lines in a given representation R of the color group are related by Casimir scaling: they are given by the quadratic Casimir operator C R times a universal (R-independent) function. Here we explicitly demonstrate that this is not so at four loops. We consider a cusp with a small angle φ, and calculate the first terms of the expansion of its anomalous dimension Γ cusp (φ), namely the φ 2 and φ 4 terms. We consider the specific color structure n f C F,4 with (1.1) The T a F denote the gauge group generators in the fundamental representation, the T a R the ones in the representation R of dimensionality N R = tr R 1, and the round brackets indicate symmetrization, see ref. [10]. This color structure cannot be represented as C R times a universal constant, and thus breaks Casimir scaling.
Non-zero quartic Casimir contributions are known to occur in closely related quantities, such as the static quark anti-quark potential [11,12], which corresponds to the φ → π limit of Γ cusp (φ). Also in N = 4 super Yang-Mills (sYM) theory contributions proportional to the quartic Casimir were found in the Bremsstrahlung function [5], i.e. the φ 2 term, and, very recently, in the light-like limit [13] of Γ cusp .

JHEP10(2017)052
Up to three loops the cusp anomalous dimension has an interesting property [8,9]: when expressed in terms of an effective coupling constant a, which is defined such that the large Minkowskian φ asymptotics, i.e. the light-like limit, of Γ cusp (φ) is given by the firstorder a term only, it becomes a universal function Ω(φ, a) that is independent of the number of fermion or scalar fields in the theory. It has been conjectured in refs. [8,9] that this property holds to all orders of perturbation theory, simply from the intriguing empirical observation at the first three orders. In the present paper we check this conjecture at four loops. The n f C F,4 term we are interested in contains the number of massless fermions n f , and, according to the conjecture, can only arise from some α n s (n > 1) term in a. It cannot be represented as a product of lower-loop color structures, and hence it can only come from the term c n f C F,4 /C R (α s /π) 4 in a/π inserted in the leading term Γ cusp (φ) = C R (a/π)(φ cot φ − 1) + O a 2 . The normalization factor c can be determined from the limit φ → π, where the four-loop Γ cusp is related [8,9] to the three-loop static potential [12,14,15].
We find that the analytic form of our result is different from the conjecture of refs. [8,9]. Interestingly, the numerical values are still surprisingly close to the conjectured ones. While this paper was finalized, the light-like QCD cusp anomalous dimension at four loops has been computed numerically [16]. Its n f C F,4 term is also relatively close, but different from the conjecture. This is in line with our findings here.
As a by-product of our calculation (at φ = 0), we determine the n f C F,4 term in the fourloop anomalous dimension of the HQET heavy-quark field. Currently it is only known to three loops [17,18]. Our result can serve as a non-trivial cross-check of future calculations.
The paper is organized as follows. In section 2, we describe our calculation. In section 3 we compute the heavy quark anomalous dimension and extract from it the QCD on-shell heavy-quark field renormalization constant. In section 4 we present the results for the cusp anomalous dimension to order φ 4 , and compare to the conjecture of refs. [8,9].

Calculation
The QCD cusp anomalous dimension arises from the UV divergences of the Wilson loop where A µ = A a µ T a R is the gluon field, P is the path-ordering operator, the trace is over (color) indices in the representation R of the gauge group. The closed integration contour C has a cusp at a single point and is smooth otherwise. Without loss of generality, we can choose the contour C to consist of two Wilson lines along the directions v µ 1 and v µ 2 with v 2 1 = v 2 2 = 1 that both extend to infinity and end at the cusp point. We denote the angle between them by φ, where φ = 0 corresponds to a Wilson line along v µ 1 = v µ 2 with both ends at infinity, and The open ends of the Wilson lines are considered to be closed at infinity. They can be interpreted as heavy quark lines in HQET with v µ 1 and v µ 2 being the heavy quark velocities. We note that for real v i in Minkowski spacetime, φ = iϕ is purely imaginary and cosh ϕ = v 1 ·v 2 . In this configuration the cusp anomalous dimension was computed through three loops in refs. [1,[6][7][8][9] and we refer to the latter reference for details on the calculational setup. 1 We distinguish two types of HQET Feynman diagrams contributing to the Wilson loop W beyond tree-level: heavy quark self-energy and one-particle-irreducible (cusp) vertex correction diagrams. The sum of the latter depends on the angle φ and is denoted by V (φ). Via a simple Ward identity the self-energy can be related to the vertex correction at φ = 0. We can thus write [6] log

JHEP10(2017)052
where we have introduced the cusp renormalization factor Z. Here and throughout this paper we use dimensional regularization with d = 4 − 2 . The cusp anomalous dimension is then given by We are interested in color structures that generate the quartic Casimir C F,4 defined in eq. (1.1). It first appears in the QCD vertex correction V at four loops and violates Casimir scaling. We will use the equality [10] where the ellipsis in eq. (2.5) stands for terms that can be expressed only in terms of the quadratic Casimirs C R , C F and C A . Equation (2.5) also holds when the order of the adjoint color indices (a, b, c, d) in one of the traces on the right-hand side is interchanged arbitrarily.

JHEP10(2017)052
The quartic Casimir C F,4 occurs in the four-loop contribution proportional to n f . The latter denotes the number of light (massless) fermions in the fundamental representation F . From eq. (2.5) it is clear that the four-loop Feynman diagrams involving C F,4 must have a light fermion loop forming a box that is connected to the Wilson lines via four gluons. There are only six different diagram topologies of that type contributing to V . They are displayed in figure 1. Counting also diagrams with reversed light fermion flow and left-right mirror graphs we arrive at a total of 18 diagrams that contribute to the C F,4 term. 2 Up to the three-loop order an analysis of all color structures and taking into account non-Abelian exponentiation [20,21] makes it possible to rewrite all non-planar integrals in terms of planar integrals only [9]. This is not the case for the non-planar diagrams (b -f ).
Using the HQET building blocks (i = 1, 2) associated with the fermion, heavy quark and gluon lines, respectively, we can write the diagrams of figure 1 in generalized covariant gauge (ξ = 0 corresponds to Feynman gauge) in compact form. The off-shellness δ/2 in the heavy quark propagators serves as an infrared regulator [9] and can be interpreted as an external energy flowing in the opposite direction of the v µ i (indicated by the arrows in figure 1). For the coefficients of the color factor n f C F,4 we have The overall minus sign in the above expressions originates from the closed fermion loop. The sum of all 18 contributions to the n f C F,4 term is gauge invariant and reads 3 We note that there is one more Casimir scaling violating color structure at four loops, namely

JHEP10(2017)052
The gauge invariance can be seen as follows: the n f C F,4 contribution in eq. (2.13) is effectively QED-like, as all diagrams have the same color structure (no relative factors). Now, consider the subdiagrams consisting of the fermion loop and four off-shell gluons attached to it in all possible ways. If we pick out one of the gluon vertices and contract it with the four-momentum of the associated gluon the contributions from the 18 one-loop diagrams add up to zero owing to the Ward-Takahashi identity of QED, see e.g. ref. [22]. Despite being off-shell the gluons are therefore effectively transverse. Thus the (longitudinal) ξ terms in their propagators vanish in the sum of all 18 diagrams. Moreover, renormalization group consistency requires eq. (2.13) to be at most 1/ divergent (for finite φ), as there are no (UV) divergent subdiagrams involved. We will use these properties as strong cross-checks of our calculation.
In this paper we are interested in the expansion of the cusp anomalous dimension for small angle φ. The calculation of the individual terms in this φ expansion is technically simpler than the calculation for arbitrary φ, as there is only one external vector v µ = v µ 1 = v µ 2 and the results are pure numbers. It can therefore be considered as a first step toward the calculation of the full angle-dependent cusp anomalous dimension, but also directly yields some relevant physical information as outlined in section 1.
In practice the Taylor expansion of the expressions in eqs. (2.7)-(2.12) can e.g. be done as follows: we write (in Euclidean spacetime) v 1 = v, v 2 = cos(φ)v + sin(φ)e 2 with v 2 = e 2 2 = 1, v · e 2 = 0 and differentiate the integrands w.r.t. φ. The numerators of the resulting terms include a number of scalar products of the loop momenta k µ i with the unit vector e µ 2 . The denominators are free of such scalar products. Upon integration therefore terms with odd numbers of k i · e 2 , i.e. odd powers of φ, vanish because of the antisymmetry of the integrand. For the terms with even numbers of k i · e 2 we perform a tensor reduction. After evaluation of the Dirac trace we end up with a scalar integrand that only involves the 14 independent scalar products k 2 i , k i · k j and k i · v. The result for each diagram can thus be expressed as a linear combination of integrals with (1) (2) numerator in our problem. Note that for brevity we have suppressed the usual −i0 Feynman prescription in the Q i , which is needed to ensure causality in Minkowski spacetime. The integrals G contributing to the vertex function V have at most 11 propagators. According to its propagator configuration each of them can be assigned to one of the integral topologies displayed in figure 2. The natural way to do this is to map the small angle expanded diagrams D a,d onto topology 1, D b,e onto topology 2, and D c,f onto topology 3. 4 We can now perform an integration-by-parts (IBP) reduction [23] for the integrals of each topology separately. We do this reduction to master integrals (MI) for all relevant integrals up to O φ 4 . To this end, we use the public computer program FIRE5 [24] in combination with LiteRed [25,26]. In this way, we find 32 MI for topology 1, 32 MI for topology 2 and 30 MI for topology 3. Taking into account relations among the MI of different topologies (with less than 11 propagators), the total number of independent MI across the three topologies is 43. Expressing the φ-expanded vertex function V as a linear combination of these 43 MI, we explicitly verify that the ξ-dependence in the coefficients of the MI drops out at O φ 0 and O φ 2 as required by gauge invariance. At order O φ 4 , we restrict ourselves to the Feynman gauge for performance reasons. We note that the extension to higher powers of φ is conceptionally straightforward, and is only limited by computing resources.
As the final step of our calculation we have to solve the MI to sufficiently high order in in order to determine the overall 1/ divergence of V related to the cusp anomalous dimension via eq. (2.4). For the computation of the MI we use the HyperInt package [27]. This code allows to automatically evaluate linearly reducible convergent Feynman integrals in terms of multiple polylogarithms. In our case the latter reduces to transcendental numbers. 5 In order to provide finite Feynman parameter integrals as an input to HyperInt we first switch to a MI basis without divergent subintegrals (except for trivial bubble insertions that can be integrated out.) In practice this is done by inserting a sufficient number of dots on the (off-shell) heavy quark lines of the MI, i.e. by increasing the power of the linear JHEP10(2017)052 propagators by an integer amount. The resulting new basis of MI is then related to the old one by IBP reduction. Next, we Fourier transform to position space and directly integrate out the simple bubble and HQET self-energy subintegrals by hand. This effectively produces integrals with non-integer -dependent propagator powers, but less than four loops. In order to avoid generating new divergent subintegrals in this process, it may be necessary to increase the powers also of some of the involved bubble propagators beforehand. The resulting integrals all have at most one overall UV divergence (∝ 1/ ).
Without loss of generality we can fix the position of the left-most vertex on the Wilson line to the coordinate origin (x 1 = 0) and parametrize the following vertices from left to right along the Wilson line by ρ x 2 , ρ(x 2 + x 3 ), etc., where 0 ≤ x i ≤ 1, i x i = 1. The parameter ρ has the dimension of a length and 0 ≤ ρ ≤ ∞. We can now write down a Feynman parameter representation for the integral over the positions of the non-Wilson-line vertices as a function of the x i and ρ. In particular its dependence on the only dimensionful parameter, ρ, can be easily deduced from dimensional power counting. The position space representation of a HQET Wilson line propagator between the vertices l and m of arbitrary power a reads dq 2π Note that because of the causal −i0 prescription δ can be considered to have an infinitesimally small negative imaginary part. The remaining overall UV divergence, if present, originates from the integration region, where all vertices on the Wilson line are contracted to one point, i.e. where ρ → 0, cf. ref. [9]. In our parametrization, the product of all heavy quark propagators in a diagram together with the ρ-dependence from the Feynman parameter integral thus yields ρ z exp(−iρ δ/2). The power z only depends on and is fixed by the dimensionality of the MI. The possibly divergent ρ integral can therefore be carried out easily. We have thus factored out the UV divergences completely and are left with a convergent integral over Feynman parameters and the x i . Its integrand can now safely be expanded in . The individual terms in this expansion are then computed with HyperInt. Let us illustrate this procedure for the following non-planar nine-propagator MI (number 38 in our list): where We have inserted a dot on the middle heavy quark line in order to render the two-and three-loop subintegrals finite. The first expression in squared brackets in eq. (2.17) arises from the product of the three Wilson line propagators according to eq. (2.16) (for the middle one a = 2). The second term in squared brackets corresponds to the Feynman parameter (α i ) representation of the integration over the two internal space-time vertices (the ones not lying on the Wilson line). Doing the ρ integral yields a factor (δ → δ − i0) , (2.20) which is 1/ (UV) divergent. The α-integrations in eq. (2.17) are projective, see e.g. [28], and we use this freedom to choose α 1 ≡ 1. With this choice the other Feynman parameters are integrated from zero to infinity. After expansion in we evaluate these convergent integrals together with the ones over x 1 and x 2 using HyperInt. The result is Here we have set the IR regulator δ = 1 for convenience. This will not affect the final result for the cusp anomalous dimension. Also the O( ) term in eq. (2.21) will not contribute to Γ cusp (φ) through O φ 4 , as can be verified by inspection of the overall -dependent coefficient of MI 38 in the IBP reduced expressions for the vertex function.
With the methods described above we have computed all MI to the relevant order in . This, in particular, includes the three eleven-propagator MI corresponding to the maximal graphs (with single propagator powers) of the three integral topologies in figure 2. For the latter already the leading term in the expansion turns out to be sufficient. We have checked our analytic results numerically with FIESTA [29]. A list of the results for the MI in electronic form can be found in the supplementary material of the present paper.
For a number of MI it is straightforward to derive four-fold Mellin-Barnes representations by applying the two-fold representation for the heavy-light vertex [30] 6 twice. They can be transformed to a four-fold series, which can be used to obtain the expansions in . We however find this method less convenient than the procedure described above.

The HQET heavy-quark field anomalous dimension
Let us denote the heavy quark self energy (i times the sum of all 1PI heavy quark selfenergy diagrams) by Σ h , and the external heavy quark energy by ω. In our configuration ω = −δ/2. The HQET (off-shell) Ward identity then relates Σ h to the vertex function at zero angle via (3.1)

JHEP10(2017)052
In fact, we have employed this identity already in eq. (2.3). Hence, we can write where Z h is the (MS) heavy quark wave function renormalization factor and V (0) is expressed in terms of the renormalized coupling α s (µ) and gauge parameter ξ(µ). We thus have for the n f C F,4 term of the associated anomalous dimension γ h , which determines the running of the renormalized HQET heavy quark field h through the renormalization group equation As field renormalization is irrelevant to physical observables γ h can depend on the gauge. The gauge dependence has been explicitly shown at three loops, where the complete result is known [17,18]. Unlike other parts, the four-loop n f C F,4 term, however, is gauge invariant for the reason discussed above. For φ = 0 our calculation of the vertex function yields (3.5)