Towards the exact Bremsstrahlung function of ABJM theory

We present the three-loop calculation of the Bremsstrahlung function associated to the 1/2-BPS cusp in ABJM theory, including color subleading corrections. Using the BPS condition we reduce the computation to that of a cusp with vanishing angle. We work within the framework of heavy quark effective theory (HQET) that further simplifies the analytic evaluation of the relevant cusp anomalous dimension in the near-BPS limit. The result passes nontrivial tests, such as exponentiation, and is in agreement with the conjecture made in [1] for the exact expression of the Bremsstrahlung function, based on the relation with fermionic latitude Wilson loops.

H Results for three-loop diagrams 39 1 Introduction Wilson loops (WLs) play an ubiquitous role in gauge theories, both at perturbative and non-perturbative level. Depending on the particular contour, they encode the potential between heavy colored sources and control important aspects of scattering between charged particles. At the same time they provide the basis for the lattice formulation of any gauge theory and represent a general class of observables with deep mathematical meaning. The AdS/CFT correspondence [2][3][4] has triggered further interest for WLs in supersymmetric gauge theories. In fact, they are directly related to fundamental string or brane configurations of the dual theory and are natural observables capable of testing the correspondence itself [5,6]. For example, the vacuum expectation values of some BPS WLs are non-trivial functions of the gauge coupling, potentially providing an interpolation between weak and strong coupling regimes. The exact computation of these quantities is sometimes possible, for instance using localization techniques [7,8], offering non-trivial tests of the AdS/CFT predictions [9,10]. When the theory is conformal, WLs govern the computation of the energy radiated by a moving quark in the low energy limit, the so-called Bremsstrahlung function B(λ) [11]. This function also enters the small angle expansion of the cusp anomalous dimension Γ cusp that, in turn, controls the short distance divergences of a WL in the proximity of a cusp, according to the universal behaviour W L ∼ exp (−Γ cusp log Λ µ ) (where Λ and µ are IR and UV cutoffs, respectively). Precisely, for supersymmetric WLs given as the holonomy of generalized connections that include also couplings to matter, we can consider a cusped WL depending on two parameters, ϕ representing the geometric angle between the two Wilson lines defining the cusp, and the latitude angle θ describing the change in the orientation of the couplings to matter between the two rays [12,13]. At θ 2 = ϕ 2 the cusped WL is BPS, therefore finite, and its anomalous dimension vanishes. For small angles, θ, ϕ 1, the expansion of Γ cusp around the BPS point reads Γ cusp (λ, θ, ϕ) ∼ B(λ)(θ 2 − ϕ 2 ) (1.1) where B(λ) is the Bremsstrahlung function, given as a function of the coupling constant of the theory. In N = 4 SYM theory an exact prescription to extract this non-BPS observable from BPS loops has been given in the seminal paper [11], where results from localization where explicitly used. The very same result can be obtained by starting from an exact set of TBA equations [14][15][16] describing the generalized cusp [17,18], solving the system in the near-BPS limit [19]. Actually, using integrability, one can further obtain a Bremsstrahlung function B L (λ), by considering the insertion of some chiral operator of R-charge L on the tip of the cusp [20]. This last result still calls for a localization explanation (see some progress in this direction in [21]). Moreover, the use of the quantum spectral curve techniques [22,23] has allowed to obtain results away from the BPS point [24]. A proposal for the Bremsstrahlung function in four dimensional N = 2 SYM appeared more recently in [25].
It is quite natural to extend the above investigations to the three-dimensional case, exploring the Bremsstrahlung function in N = 6 superconformal ABJ(M) theories [26,27]. One immediately faces a first difference with the four-dimensional case: In ABJ(M) models not only bosonic but also fermionic matter can be used to build up generalized (super)connections whose holonomy gives rise to supersymmetric loop operators. Precisely, there are two prototypes of supersymmetric Wilson lines, one associated to a generalized gauge connection that includes couplings to bosonic matter only, so preserving 1/6 of the original supersymmetries [28][29][30][31], and another one in which the addition of local couplings to the fermions enhances the operator to be 1/2-BPS [32]. The latter should be dual to the fundamental string on AdS 4 × CP 3 . While the 1/2 BPS operator is cohomologically equivalent (and therefore basically indistinguishable) to a linear combination of 1/6 BPS ones, the fact that they preserve different portions of supersymmetry allows for constructing different non-BPS observables starting from them. Generalized cusps formed with 1/6-BPS rays or 1/2-BPS rays are actually different [33,34] and, consequently, different Bremsstrahlung functions can be defined and potentially evaluated exactly. In particular, in [34] a proposal for the exact Bremsstrahlung function of the 1/6-BPS cusp was put forward, based on the localization result for the 1/6-BPS circular WL, and an extension to the 1/2-BPS case was conjectured.
In [1] an exact formula appeared that expresses the Bremsstrahlung function B 1/2 for 1/2-BPS quark configurations as the derivative of a fermionic latitude WL with respect to the latitude angle. This proposal was suggested by the analogy with the N = 4 SYM case [11] and supported by an explicit two-loop computation that agrees with the result obtained directly from the cusp with 1/2-BPS rays [33]. Given that the fermionic latitude WL is cohomologically equivalent to a bosonic latitude WL [1], B 1/2 can be eventually expressed as in terms of the well-known vev of 1/6-BPS WLs associated to the two U (N ) gauge groups and computed exactly at framing one using localization [8,35,36]. Remarkably, this formula coincides with the proposal of [34], despite being derived from different arguments. Such an expression is also in agreement with a one-loop string computation at strong coupling [37][38][39].
Expanding the matrix model result for W 1/6 and Ŵ 1/6 at weak coupling one can infer that in the planar limit this expression contains odd powers of the 't Hooft coupling λ = N/k only. In principle, there should be a second contribution to (1.2) proportional to the derivative of a multiply n-wound 1/6-BPS circular WL with respect to n. This term should provide the even power expansion in λ. However, using the exact expression for the n-wound WL known from localization, we can verify that both at weak and strong coupling this term vanishes (at least) in the planar limit [1,34]. In particular, up to two loops it vanishes also at finite N .
Based on this observation and using the weak coupling expansion of W 1/6 and Ŵ 1/6 as derived from the matrix model (including the color subleading contributions), eq. (1.2) gives the following prediction for the B 1/2 function, at any finite N 3) The first term has been already checked in [1]. In order to put this conjecture on more solid bases, further checks at higher orders are desirable. In this paper we perform a three-loop computation of B 1/2 directly from its definition in term of the fermionic cusped WL. As a result we obtain exactly the predicted 1/k 3 coefficient in (1.3), so providing a highly non-trivial test of the proposal [1,34]. In particular, at least up to three loops, we have checked that (1.2) is valid not only in the planar limit, but also at finite N . We remark that since we are checking the Bremsstrahlung function through an odd perturbative order, our test is independent of any assumption on the even-λ part that was discarded in (1.2), according to the discussion above.
To simplify the three-loop calculation, we rely on the BPS condition satisfied by the fermionic cusped WL [33] that implies the small angle expansion (1.1). From this equation we can read the Bremsstrahlung function equivalently from the θ or the ϕ expansions of Γ 1/2 (λ, θ, ϕ), setting the other angle to zero. It turns out to be particularly convenient to set ϕ = 0 (flat cusp limit), working with the internal angle θ only. As will be discussed in details, this choice leads to crucial simplifications at several stages of the calculation.
As usual, when computing the cusped WL in the near-BPS limit UV and IR divergences arise. Here we employ dimensional regularization with dimensional reduction (DRED) to regulate the UV divergences and introduce an IR suppression factor to tame the IR divergences coming from the integration region at large distances. The cusp anomalous dimension is then obtained through the usual renormalization group equation for the renormalized cusped WL [40]. The calculation is carried out in momentum space where the Wilson lines are effectively described by non-relativistic, eikonal propagators. This procedure has a physical interpretation in the context of heavy quark effective theory (HQET) [41][42][43].
The plan of the paper is the following. In section 2 we review the construction of 1/2-BPS WLs in ABJM theories and present the proposal for the exact 1/2-BPS Bremsstrahlung function, as derived in [1]. In section 3 we outline the strategy of the three-loop computation, describing in particular the HQET formalism. Section 4 contains a review of the one-and two-loop results that are important in order to check the correct exponentiation. The actual three-loop calculation appears in section 5, while the Bremsstrahlung function is extracted in section 6. Section 7 contains our conclusions. Many technical details, together with the results for the relevant two-and three-loop diagrams are deferred to eight appendices.

The 1/2-BPS cusp in ABJM
We start reviewing the definition of the locally 1/2-BPS cusp in ABJM theory 1 . This is constructed by considering the fermionic WL 2 evaluated along the contour Γ formed by two rays lying in the (1, 2) plane and intersecting at the origin, as illustrated in figure 1. The angle between the rays is π − ϕ, such that for ϕ = 0 they form a continuous straight line. Explicitly, the contour drawn in figure 1 is parametrized as In (2.1) L is a superconnection belonging to the Lie-superalgebra U (N |N ) and given by Generalities on the ABJM theory in our conventions are given in appendices B, C. 2 Here the path-exponential is defined by Edge 2 Edge 1 Figure 1. Eq. (2.2) represents a planar cusp, whose angular extension is given by π − ϕ.
Here Tr denotes the standard matrix trace (and not the super-trace, see [33] for a detailed discussion). The fermionic couplings on each straight-line possess the factorized structure [32,33] η The lowercase index i distinguishes the two different rays forming the cusp. On the first edge, the spinor and R-symmetry factors are respectively given by while on the second one we have (2.6) The θ angle is the counterpart of ϕ in the R-symmetry space. It denotes the angular separation of the two edges in the internal space (i.e. CP 3 ). The two matrices which couple the scalars are determined in terms of the vectors n iM andn N j through the relations M I iJ = δ I J − 2(n I in iJ ). We have respectively As discussed in details in [33] the presence of two angles allows us to find BPS configurations. Explicitly, for θ 2 = ϕ 2 the two straight lines forming the cusp share two Poincarè supercharges and two conformal supercharges. In other words we deal with a globally 1/6−BPS WL. As a consequence, the divergent contributions arising from the cusp geometry are expected to vanish if the BPS condition is satisfied, a fact which was explicitly checked up to second order in perturbation theory [33].

The exact 1/2-BPS Bremsstrahlung function
Away from the BPS point θ 2 = ϕ 2 , the cusped WL defined above suffers in general from UV divergences. At generic angles it has to be renormalized and we expect that the operator possesses an anomalous dimension Γ cusp , according to the universal behaviour where Λ is an IR cutoff and µ stems for the renormalization scale. For the familiar bosonic WL, the divergences associated to a cusp singularity along the contour can be cast in the exponential form (2.8) [40]. This property follows from the presence of general exponentiation theorems for these operators. More precisely, their expectation value can be written as the exponential of the sum of all two-particle irreducible diagrams [42,44]. Extending these results to our case is not straightforward. In fact, the additional couplings between the contour and the fermions appearing in the superconnection may affect the standard analysis presented there. However, despite of these potential technical issues, we expect the trace of the Wilson line in ABJM to respect the usual exponentiation process, so that we can define an anomalous dimension for the cusp according to the standard text-book procedure. Our three-loop result will explicitly confirm the consistency of this picture (see discussion in section 6).
In the limit where the cusp angle is small, the cusp anomalous dimension has a Taylor expansion in even powers of ϕ, whose first coefficient is defined as the Bremsstrahlung function [11]. The 1/2-BPS cusp Γ 1/2 , associated to the 1/2-BPS Wilson lines (2.1) of the ABJM theory, vanishes for ϕ = ±θ and consequently the small angle expansion reads where B 1/2 is a function of the coupling and the number of colors only.
In [1] a conjecture was put forward for the exact expression of B 1/2 . The argument was articulated in a few steps, which we briefly review here.
• First, the Bremsstrahlung function was related to the expectation value of a supersymmetric fermionic circular WL W • F evaluated on latitude contours in the S 2 sphere, that is displaced by an angle θ 0 from the maximal circle, and potentially by an additional internal angle α in the R-symmetry space. Remarkably, the expectation value of such WL seems to depend only on a particular combination of the angles, ν = sin 2α cos θ 0 [1]. The undeformed contour yields the 1/2-BPS circular WL of [32] and corresponds to ν = 1.
Paralleling an analogous derivation for N = 4 SYM [11], the conjecture of [1] states that the Bremsstrahlung function is retrieved by taking derivatives of the WL expectation value with respect to the characteristic parameter ν measuring the deformation Unfortunately, in contrast with the N = 4 SYM case, the latitude WLs have not been yet evaluated via localization, so more work is necessary to derive an exact formula for the Bremsstrahlung function.
• As a second step, the expectation value of the latitude WL W • F was argued to be expressible, via a cohomological equivalence, in terms of bosonic WL operators W • B andŴ • B preserving lower supersymmetry. Using this equivalence, eq. (2.10) can then be rewritten as where Φ B stands for the phase of the WL W • B . Again, in the ν = 1 limit, such WL operators land on the 1/6-BPS circular WL of [29][30][31] and the cohomological equivalence parallels that for 1/2-BPS operators, described in [32].
• The expectation value of W • B operators is not known exactly, either. Nevertheless, it was conjectured that at least the relevant derivative in (2.11) could be traded with the derivative of the expectation value of a multiply wound 1/6-BPS circular WL W 1/6 n with respect to the winding number n, upon proper identification of the latitude and winding parameters n ∂n(ν) ∂ν ν=1 (2.12) Such a trick was also proposed in [34] in the case of the Bremsstrahlung function associated to a cusp constructed with two locally 1/6-BPS rays in ABJM theory. Multiply wound supersymmetric WLs were explored in details in [45], where their expectation value was given using a matrix model average description.
• The expansion of the results for W 1/6 n and Ŵ 1/6 n provides evidence that (2.12) vanishes at all orders, as advocated in [1,34], at least in the planar limit. Under these circumstances, we then arrive at the following simple proposal for the 1/2-BPS Bremsstrahlung function at all orders 3 The Bremsstrahlung function associated to a 1/2-BPS cusped WL is given in terms of the expectation value of the 1/6-BPS circular WL, well-known from localization [8,36]. In [33] it has been checked that at two loops (2.12) vanishes also for finite N . Therefore, up to three loops expression (1.3) provides the expected result for B 1/2 , without requiring the planar limit. This is the prediction that we are going to check with an explicit three-loop calculation of Γ 1/2 .
Somewhat unusually in ABJM theory, such an object seems to possess an odd dependence on the coupling k. Nevertheless, the proposal is in agreement with previously derived results for the cusp anomalous dimension at weak coupling [33], that is up to two-loop order and up to subleading order at strong coupling [37][38][39]. In this paper we put the conjecture on much firmer grounds, by providing an explicit check of it up to three loop order in perturbation theory, including the sub-leading corrections in N , which first appear at this order.

Overview of the computation
Referring to eq. (2.9), the Bremsstrahlung function can be read equivalently from the θ or the ϕ expansions, setting the other angle to 0 We take advantage of this fact for simplifying the evaluation of the Bremsstrahlung function, circumventing the computation of the complete cusp anomalous dimension. In particular, it proves convenient to set ϕ = 0 and work with the internal angle θ only, as shown in the cartoon of figure 2. This strategy brings dramatic simplifications at several stages of the computation, both in the number of diagrams to be taken into account and in their evaluation. φ θ Figure 2. For the 1/2-BPS cusp the small ϕ limit at null internal angle (the red arrow in the cartoon stems for a direction in the internal space) is equivalent to the small θ limit at zero ϕ.

The diagrams
We now discuss how to suitably choose the computational setting in order to pin down the calculation to a relatively small number of Feynman diagrams to be considered at three loops.
First of all, we recall that in Chern-Simons-matter theories remarkable simplifications arise when the WL contour lies in a two-dimensional plane. In fact, this choice triggers the vanishing of several tensor contractions, by virtue of the antisymmetric nature of Chern-Simons gauge propagators (C.2) and cubic vertex (C.5) in their index structure. In particular, whenever an odd number of antisymmetric tensors ε µνρ arises in the algebra of a diagram, their product eventually vanishes. In fact, one can always reduce them to a single Levi-Civita tensor whose indices have to be contracted with three external vectors. The latter all come from the WL contour and hence, lying on a plane, are not linearly independent and always give vanishing expressions when contracted with the ε tensor.
This feature of perturbation theory becomes particularly powerful at odd loop orders. Namely, from ABJM Feynman rules (see appendix C), one can infer that purely bosonic diagrams with gauge interactions vanish identically, by virtue of the aforementioned mechanism. For supersymmetric WLs this singles out contributions containing matter exchanges only. Nonetheless, there are still a number of simplifications.
For example, bosonic 1/6-BPS WL diagrams with scalar exchanges also evaluate to zero, thanks to the vanishing of the trace of their coupling matrices M . Therefore, in ABJM theory bosonic WLs with planar contours, when computed at framing zero, automatically have vanishing expectation values at odd loops [31] As a consequence, only fermionic WLs of the Drukker-Trancanelli type [32] allow for nontrivial contributions at odd orders, as shown in [47]. This is precisely the class of WLs we are computing in this paper.
The arguments provided above imply that in our computation we have to consider only graphs containing at least one matter propagator. This condition substantially narrows down the number of diagrams to be considered.
A further simplification arises from performing the calculation at ϕ = 0, as previously discussed. In fact, taking the contour to lie on a straight line triggers the additional vanishing of diagrams such as the one in figure 3, again by antisymmetry in the contraction of ε tensors with vectors on the line. Finally, aimed at computing the Bremsstrahlung function, one could in principle restrict to diagrams producing factors of cos θ/2 only (see eq. (3.1)). According to the rules in appendix E, these occur only if a fermion line or a scalar bubble stretch between opposite sides of the cusp. Nevertheless, it turns out that we can easily compute the whole cusp at ϕ = 0 including also θ independent terms.
Technically, evaluating the various diagrams entails a few steps. After using the Feynman rules reviewed in appendix C, the integral associated to a generic diagram with fermion exchanges has the following structure dτ 1>2>···ẋ whereẋ i ≡ẋ(τ i ) and we have used the shortening notation The number of internal and contour integrations, spinor structures, derivatives and propagators depend on the particular diagram. Lorentz indices carried by velocities, gamma matrices and derivatives are contracted among themselves or with external ε tensors. If couplings with scalars are also present, Mmatrix structures will appear.
We use identities in appendix E to reduce the η bilinears and the M-matrix structures. We perform the relevant tensor algebra strictly in three dimensions and in an automated manner with a computer program. In this process we drop all the terms containing an odd number of Levi-Civita tensors and reduce all the products of an even number of them to products of metric tensors. We finally arrive at expressions containing scalar products of velocities and derivatives only, where the ultimate step consists in integrating over internal vertices and WL parameters.

The integrals and the HQET formalism
The most striking simplification of working at ϕ = 0 is that the integrals arising from Feynman graphs reduce to propagator-type contributions which evaluate to numbers rather than functions of the angle, considerably reducing the effort needed for their determination.
The integrals could in principle be evaluated directly by solving internal integrations in configuration space and then integrating over the Wilson line parameters. In this case one could make use of the Gegenbauer polynomial x-space technique GPXT [48] to solve the three-dimensional internal integrals along the lines of [49]. However this approach becomes cumbersome and impractical for contributions with more than one internal integration, such as those appearing at three loops.
A more powerful strategy consists in Fourier transforming the integrals to momentum space and perform the contour integrations first. This step casts the various contributions into the form of the non-relativistic Feynman integrals appearing in the effective theory of heavy quarks (HQET) (see [50,51], for a recent application in four dimensions), or rather an Euclidean version thereof, as we perform the whole computation with (+, +, +) signature. This procedure, applied to a single propagator, is sketched pictorially in figure 4. A more detailed application is described in sections 4 and 5.2, when dealing with the example of the one-loop and a particular graph of the three-loop computation. The velocity of the heavy quark is the vector tangent to the WL contour, v =ẋ(τ ).
In general, the resulting integrals suffer from both IR and UV divergences. UV divergences are treated within the framework of dimensional regularization, shifting spacetime integrations to d = 3 − 2 dimensions. As a consequence, the integrals evaluate to Laurent series in the regularization parameter . In order to preserve supersymmetry we work in the dimensional reduction scheme (DRED) [52]. In practice this translates in performing the tensor algebra in the numerators strictly in three dimensions. This prescription can introduce evanescent terms if some tensor integration is attained (in d = 3 − 2 ) and then contractions with three dimensional objects are carried out. Such evanescent terms can be taken into account by suitable regulator dependent factors (see e.g. [53][54][55] for a discussion in the context of WLs in three dimensions). However, in the present case this subtlety is already built-in in our method and does not cause any practical worry. In fact, the index algebra is always handled (in integer dimensions) before any integration, as recalled in the previous section. After this reduction, the only potentially dangerous integrals occur whenever a tensor numerator contains contractions with an ε tensor, but such contributions were argued to vanish and promptly discarded.
In order to tame IR divergences at long distances along the WL contour, following [51] we introduce an exponential damping factor e δ τ (Re(δ) < 0) for each contour integral, enforcing the finiteness of the integrals at large radius. This introduces a residual energy in the heavy quark, shifting the HQET propagators by The parameter δ is customarily set to δ = −1/2, a choice which simplifies the result of the relevant integrals. The final result for the cusp anomalous dimension is scheme independent and is not affected by the particular value of δ, albeit all intermediate steps might be. Up to three loops, HQET propagator integrals were studied in detail in [56,57] for QCD in four dimensions. The analysis performed there reveals that all topologies can be reduced to a set of 1, 2 and 8 master integrals (seven planar and one non-planar) at one, two and three loop order, respectively, by means of integration by parts identities [58,59]. Such a statement is true, independently of the space-time dimensions and applies to our case as well. This reduction step can be made automatic thanks to the Laporta algorithm [60,61] and carried out with one of its available implementations [62][63][64][65][66]. In practice we have used FIRE5 [64], whose C++ version is able to perform the reduction of all the required integrals in a few minutes on a PC.
The relevant master integrals have to be evaluated in d = 3 − 2 dimensions. This task can be carried out straightforwardly relying on previously derived results for such topologies with general integration dimension d and powers of the propagators [56] (some of them are non-integer because arising from the integration of bubble subtopologies). The final evaluation of the master integrals and expansion in series of up to the desired order for the cusp anomalous dimension is spelled out in appendix F.
The procedure described above that reduces the cusp loop integrals to propagatortype integrals for non-relativistic heavy particles, has a very suggestive physical interpretation. The anomalous dimension of the cusped WL captures the renormalization of the current associated to a massive quark, undergoing a transition from a velocity v 1 to a velocity v 2 at an angle ϕ. In fact Wilson lines have been thoroughly employed as a convenient tool for describing the dynamics of such objects. According to this physical picture, the 1/2-BPS Wilson line in ABJM theory is associated to heavy quarks, which are interpreted as N = 3 supermultiplet W bosons, transforming in the fundamental representations of the two gauge groups. These were shown to arise when moving away from the origin of the moduli space [67,68].
The Bremsstrahlung function controls the energy loss by radiation of these massive particles undergoing a deviation by an infinitesimal angle. In particular, the Higgsing procedure mentioned above reveals that this superparticle is able to radiate off pairs of bifundamental fields, both scalars and fermions, quite exotically. The coupling to such matter fields enters the superconnection (2.3) of the 1/2-BPS Wilson line via a parameter θ, which can be interpreted as an angle in the internal SU (4) R R-symmetry space. According to the discussion above, we determine the Bremsstrahlung function by considering the equivalent (but computationally simpler) picture of a heavy quark subject to a small kick in R-symmetry space at fixed vanishing geometrical angle.
In particular, thanks to the ϕ = 0 condition, the integrals involved in the computation of the Bremsstrahlung function are precisely those contributing to the self-energy corrections of a heavy quark, i.e. propagator-type. The presence of the cusp point on the line has the simple effect of increasing the power of a HQET propagator in the diagram. Integration by parts identities can then be employed to reduce the power of the doubled propagator to unity.

The cusp anomalous dimension at zero angle
In the previous subsection we explained how the computation of the expectation value of the cusped WL is performed up to three loops, in the flat cusp limit. In this section we provide the details on how the cusp anomalous dimension and the Bremsstrahlung function can be practically extracted from this result.
As recalled above, the expectation value of the cusped WL suffers from UV and IR divergences. The former determine the renormalization properties of the WL and are the crucial object of our investigation. The latter can be regulated in different ways, either by assigning the heavy quarks a residual energy offsetting them from the mass shell as described above or by computing the loop on a finite interval (−L, L) as done in [33]. Although these divergences are uninteresting for our purposes, different kinds of regularizations can be more or less convenient according to the approach we use to compute the integrals. Therefore, while discussing how to extract the UV divergent behavior of the cusped WL, we will also compare the effectiveness of different IR regularization strategies.
For a cusped WL in a gauge theory UV divergences may have different origins. The first source consists in the divergences of the Lagrangian. In particular, this can cause a running coupling and the presence of a nontrivial β function. The ABJM model is superconformal and such a phenomenon does not occur. Then, divergences associated to the short distance behavior of the Wilson line can arise. In the language of HQET these are the (potentially divergent) radiative corrections to the self-energy of the heavy quarks. Finally, the cusp geometry itself introduces further divergent contributions, which are the ones we are interested in.  The systematics of the renormalization properties of WLs with cusps were discussed in [40,69]. Here we recall some basic and general features.
First of all we can focus our attention only on contributions which are 1PI vertex diagrams in HQET language. For instance if we deal with a configuration of the type represented in figure 5 the subsector governed by the Green function G(t r+1 , ··, t n ) completely decouples from the rest. In fact the corresponding contribution is factorized as If we perform the change of variables t i → t i + t r (with i = r + 1, · · · , n) and we use that G(t r+1 , ··, t n ) is a translationally invariant function we find The factor (B) in (3.5) provides the value of the diagram obtained from the one in figure 5 removing the subsector governed by G(t r+1 , ··, t n ), while the factor (A) is the contribution due to G(t r+1 , ··, t n ) to the vacuum expectation value of a straight-line running from −∞ to 0. In HQET language this factorization is a trivial manifestation of momentum conservation. On the other hand, we have to stress that eq. (3.5) also relies on the judicious choice of the IR regulator which breaks translation invariance along the two rays of the cusp in a controlled way. For instance, if we were to tame the IR divergences by considering a cusp with edges of finite length L, it is easy to see that the above argument would fail and the two contributions would remain in general intertwined.
Because of (3.5) the cusped WL is given by the sum of all the 1PI diagrams, which we shall denote by V The two Wilson line factors in HQET language are identified with the two point functions of the heavy quark and eq. (3.6) is the usual and well-known decomposition of the vacuum expectation value in terms of its 1PI sector and two-point functions.
To single out the cusp anomalous dimension from W (θ, ϕ) , we have first to eliminate the spurious divergences which are due to the fact that the presence of the IR regulator e ±δt mildly breaks gauge invariance 5 . A similar phenomenon also occurs when we consider edges of finite length L [69][70][71][72]: There the gauge invariance is lost because the loop is open. In both cases the gauge-dependent divergent terms can be eliminated by introducing a multiplicative renormalization Z open [69][70][71][72], which, in practice, is equivalent to the subtraction of the straight line or the cusp at θ = ϕ = 0 with V (θ, ϕ) defined in (3.6). In the formalism of [41] Z open corresponds to the wavefunction renormalization of the fermions in the one-dimensional effective description of the Wilson line operator. This factor needs to be eliminated in order to single out the proper renormalization of the relevant operator describing the cusp [40]. Finally we can extract Z cusp by defining the renormalized WL as follows The cusp anomalous dimension is then evaluated through the relation where µ stems for the renormalization scale and we have suppressed the dependence on the angles. The perturbative computation yields V (θ, ϕ) as an expansion in 1 In dimensional regularization every V (i) coefficient is given by a Laurent expansion in . According to the standard textbook dictionary, the cusp anomalous dimension can be read from the residues of the simple poles of Z cusp in the plane This completes the evaluation of the cusp anomalous dimension, which in our case is restricted to the ϕ = 0 limit Γ 1/2 ϕ=0 . Finally the Bremsstrahlung function can be evaluated by taking the double derivative

Review of one-and two-loop results
The cusp anomalous dimension is extracted from the divergent part of the logarithm of the cusped WL expectation value. In the absence of a well-established exponentiation theorem for WLs in ABJM, we do not compute this logarithm directly, rather we consider the whole WL and take the logarithm perturbatively. In particular at three loops, this entails subtracting contributions from the one-and two-loop order corrections. At the order in the expansion to which we compute the three-loop expectation value and its logarithm, namely up to 1/ terms, we need to consider the one-loop correction up to order O( ) and the two-loop contribution up to finite terms in the regulator. In this section we perform such a computation. The relevant 1PI diagrams were already determined in [33], but here we expand them to higher orders in and treat them within the HQET formalism, for consistency with the three-loop results. This produces a mismatch with respect to the results in [33] by scheme dependent terms. However, the leading poles in the regulator are expected to coincide, diagram by diagram. This is indeed the case and provides a handy consistency check.
At one loop, only one non-vanishing diagram arises, which is a fermion exchange across the cusp point, as depicted in figure 6. Throughout this computation, the According to the method outlined above, we Fourier transform this to momentum space. Introducing the exponential IR regulator δ, we obtain In the last line, the squared propagator arises from the two sides of the cusp which degenerate to the same HQET propagator in the ϕ = 0 limit.
We set δ = −1/2 and absorb the imaginary unit in the HQET propagator into the velocity v = iṽ. The resulting imaginary vector withṽ 2 = −1 is suitable for the evaluation of the master integrals in euclidean space, making them manifestly real.
Then we make use of identity (E.3d) to express the η-bilinear in terms of the external velocity v. This gives rise to the θ dependence of the diagram, according to (E.3b), and expresses the numerator of the integrand as k · v. This in turn can be rewritten as an inverse HQET propagator, so that at this stage the diagram can be expressed as a linear combination of scalar HQET one-loop integrals (F.1) where C θ = cos θ 2 . An integration-by-parts identity allows to express everything only in terms of the master integral G 1,1 Finally, using formula (F.7) we expand this integral in up to order O( ) as required by the three-loop computation and obtain the one-loop 1PI HQET vertex correction at zero cusp angle up to an overall factor (e −γ E 16π) k , omitted to keep the expression compact.
At two loops the relevant diagrams are shown in figure 7. In general, each topology of diagrams is associated to a number of independent 1PI configurations, according to the position of the cusp point. Different configurations give rise to different powers of cos θ/2. As an example, in figure 7 we depict the two possible configurations of diagrams with double fermion exchanges, (2) e 1 and (2) e 2 , which give contributions proportional to cos θ/2 and cos 2 θ/2, respectively.
All the diagrams are computed in the HQET framework, following what have been described above for the one-loop case. The integration-by-parts reduction lands on to two master integrals only, which are defined and evaluated in appendix F. The list of results for each individual diagram is given in appendix G. Combining these expressions and expanding them up to order O( 0 ) leads to the following (unrenormalized) expectation value of the 1PI part for the cusped WL at two loops . We stress that no sub-leading corrections in N arise up to two loops.

Diagrams
At three loops it is convenient to classify the diagrams according to the number of fermion and scalar insertions. We provide only a sketch of the possible topologies of diagrams arising in the computation, but each of them can contribute with several different relative positions of the insertion and the cusp points, all giving rise to 1PI configurations, and/or different choices of the gauge vectors, A µ orÂ µ .
The planar diagrams can be organized according to • Diagrams containing scalar insertions, figure 8.
(3) a is the only diagram obtained by contracting three bi-scalar insertions, with no fermion lines. This, along with diagram (3) b , does not actually contribute to the Bremsstrahlung function. In fact they possess a trivial dependence on the internal angle θ, as is immediate to infer using identities (E.3). Nonetheless, we take them into account in order to provide the full result for the cusp at vanishing  Nevertheless, these evaluate to zero after performing the algebra and integrating. Therefore we have not displayed them here.
• Diagrams with no scalar insertions and a fermion line attached to the WL, as depicted in figure 9. Such diagrams arise in a variety of fashions, with up to three internal integration vertices. Several graphs feature lower order corrections to the field propagators.
We recall that the one-loop correction to the scalar two-point function vanishes identically and the one for the fermion line (C.8) is proportional to the difference between the ranks of the two gauge groups, hence it vanishes in the ABJM theory. Diagrams (3) f and (3) i involve the one-loop gluon self-energy correction, whose result is reviewed in appendix C. In our computation we employed an effective propagator in momentum space with an unintegrated scalar bubble, in order to make the corresponding cusp contribution fit into the integral classification (F.1). Diagram (3) j contains the two-loop self-energy correction of the fermion propagator. Its computation in configuration space is spelled out in appendix D.
• Diagrams with four fermion insertions along the line, as drawn in figure 10.  These are all ladder contributions with no internal integrations, making their evaluation easier than the previous contributions. In the cartoon we have only represented all possible 1PI vertex contractions, which are relevant for our computation.
Finally, nonplanar diagrams can be constructed starting from the same field insertions as in graphs (3)

The computation
We tackle the evaluation of three-loop Feynman diagrams by applying the HQET method described in section 3. The integrands related to the diagrams listed in the previous section can be computed using the Feynman rules in appendix C or using the Mathematica R package WiLE [73]. Since the complete calculation is quite long and cumbersome, here we spell out the explicit details of one particular diagram.
As an example, we consider diagram (3) h , selecting in particular the configuration drawn in figure 9 plus the one where the gluon insertion on the contour and the cusp point are interchanged. The two configurations, being symmetric and thus equivalent, yield an overall factor 2. To be definite, we call the sum of these two configurations (3) 1 h . Further possible configurations from the (3) h topology arise, which correspond to other orderings of insertion points on the contour.
Considering configurations (3) 1 h , we may obtain leading and sub-leading contributions, depending whether the gauge lines correspond to A µ orÂ µ .
We start computing the planar contributions. The starting string in configuration space is given by where the constant factor K reads Next, we Fourier transform this expression using (C.1) (which eats up part of K) and arrive to an equivalent HQET-like contribution where we have already introduced the IR regulator δ = −1/2. Next we use (E.3d) to reduce the η-bilinear and extract the corresponding θ dependence. We perform the tensor algebra within the DRED scheme, using in particular a repeated (automatized) application of (A.2) to simplify the product of γ matrices to products of metrics and Levi-Civita tensors. In this process we discard all the contributions proportional to one final ε tensor, as it is doomed to vanish after integration, as explained in section 3.1. The result of this procedure leads to a numerator which can be expressed in terms of scalar products between integrated momenta and external velocities only. It looks like As a final step, we rewrite the above expression in terms of inverse propagators and arrive at a sum of scalar integrals with different powers of the propagators, in a form which FIRE can be fed with. After the FIRE digestion, the diagram evaluates to a simple sum over the master integrals listed in appendix F 1,1,1,1,0,0,0,1 1,1,1,0,0,1,0,1 + 8(d − 1)G 0,1,1,0,1,1,1,0,1 + (32 − 14d)G 1,0,0,0,1,1,1,0,1 (5.4) which can directly be expanded in to the desired order Repeating the exercise for all the diagrams we obtain the results listed in appendix H.
We compute the nonplanar corrections as well, which are produced by two possible sources. At first they can come from planar diagrams, when for instance the choice of the vector boson gives rise to a sub-leading in color contraction of the color matrices. Second, they can be produced by genuinely nonplanar diagrams. We remark that, since we consider color subleading contributions, we also have to take into account the mixed gauge propagator AÂ contribution (C.7c), arising at one-loop from the creation and annihilation of a matter pair.
The first class of contributions arises from diagrams (3) a , (3 figure 12, where we have shown only one particular nonplanar configuration for each topology. Thanks to the linear nature of the HQET propagators as functions of loop momenta, such nonplanar integrals contain linearly dependent propagators. These can be manipulated by partial fractioning (using an automated routine) in such a way to reduce them all to either planar integrals or a single nonplanar topology depicted in figure 15, with arbitrary powers of the propagators. These can be finally reduced by integration-by-parts identities to a single nonplanar master integral, which is defined and computed in appendix F. In particular, after decomposition by partial fractioning, the nonplanar diagrams arising from topologies (3)

The three-loop cusp and the Bremsstrahlung function
We are now ready to provide the full result for the 1/2-BPS cusp anomalous dimension up to three loops, at vanishing geometric angle ϕ = 0. Omitting an overall factor , the sum of all the three-loop contributions to the 1PI vertex function V reads As a first consistency check of the result, we probe exponentiation (2.8). This test was already successfully performed up to two loops in [33] and here we extend the analysis to the next perturbative order. To this end we insert the explicit expressions (4.5, 4.6, 6.1) for V (1) , V (2) , V (3) in the general expansion for the vertex function (3.10). Now, taking the logarithm and expand the resulting expression in powers of 1/k we do find Although the V (2) and V (3) coefficients originally contained higher order poles in 1/ , only single poles appear in (6.2). This proves that quadratic and cubic poles in V (2) (θ) and V (3) (θ), eqs. (4.6, 6.1), are ascribable to powers of lower order divergent contributions to V . Therefore, as expected, the divergent part of the expectation value of the cusped WL exponentiates. In particular, since no nonplanar corrections appear at one-and two-loop orders, their contribution at three loops has only a simple pole in the regulator. We note that the two-loop part in (6.2) is independent of θ. Applying prescription (3.7) to obtain the whole expression for the cusped WL this contribution will then get removed. The final result for the (unrenormalized) cusped WL at ϕ = 0 reads 3) This expression displays uniform transcendentality. Moreover, it satisfies the BPS condition, namely it vanishes for θ = 0, albeit this statement is quite trivial at ϕ = 0, being it an obvious consequence of (3.7).
As a further consistency check, we have also directly verified prescription (3.7) by computing the non-1PI contributions as well.
We renormalize (6.3) by multiplying by a renormalization function Z cusp . Following the steps outlined in section 3.3, from Z cusp we obtain the final result for the anomalous dimension at three loops The Bremsstrahlung function is finally evaluated according to (3.12) and reads This result is in perfect agreement with conjecture (2.13) that was formulated in [1,34]. In fact, using the expectation value for the 1/6-BPS WLs at framing one, as derived from localization and inserting them in the r.h.s. of (2.13) we find exactly expression (6.5). Remarkably, the agreement applies also to the nonplanar part. This is consistent with the fact that up to three loops the arguments of [1] do not rely on restricting to the planar limit.

Conclusions
In this paper we have studied the Bremsstrahlung function for the locally 1/2-BPS cusp in the ABJM model at weak coupling with the aim of proving conjecture (2.13) for its exact determination. We remind that this conjecture is based on linking the Bremsstrahlung function to the expectation value of supersymmetric 1/6-BPS WLs with multiple winding, which are computable exactly via localization [45]. To begin with, it relies on the expectation values of supersymmetric WLs on latitude contours in S 2 (see eq. (2.10)) and hinges on a few steps, some of which lack a rigorous proof. Therefore, in order to substantiate the conjecture, we have performed a precision perturbative check at three loops. Technically, we have taken considerable advantage from the BPS condition for the cusp, in order to make the computation simpler. We have considered the expansion of the cusp anomalous dimension in the small internal angle in the R-symmetry space, at vanishing geometric angle. This setting entails several simplifications at both the level of the diagrams involved, and in their practical evaluation, especially in handling the integrals. These reduce to propagator-type integrals, which we computed by switching to the HQET picture (via Fourier transform to momentum space) and employing integration by parts identities to reduce them to a restricted set of master topologies.
Result (6.5) for the three-loop Bremsstrahlung function agrees with the prediction based on conjecture (2.13). Remarkably, such an agreement extends to the nonplanar part as well. This provides a compelling test in favour of proposal (2.13) and its validity beyond the planar limit. This exact result is interesting per se, but acquires additional appeal in view of a potential integrability based computation of the same quantity. Recent developments on the Quantum Spectral Curve approach [22][23][24]74] in the ABJM model [75,76] seem to point to such a possibility. This would provide a crossed check of the localization based proposal and, by comparison with the integrable description, of the conjecture on the exact expression for the interpolating h function of ABJM [77][78][79][80].
We conclude with comments on future perspectives and extensions of this paper. A first, sensible generalization of the three-loop computation of the cusped 1/2-BPS WL consists in lifting the ranks of the two gauge groups to unequal values, or in other words to compute the same quantity in the ABJ model. From the technical standpoint, this requires modifying the color factors of the diagrams analysed in this article and supplying further graphs which were discarded, because they vanish in the ABJM limit.
Despite the fact that such a result is expected to be straightforward to derive, lower order contributions suggest that the expectation value of the locally 1/2-BPS WL on a cusped contour does not exponentiate [33,81], at least not in a standard fashion. This might hamper the interpretation of the three-loop correction and the extraction of a cusp anomalous dimension from it. Moreover, it is not clear if and how conjecture (2.10) can be extended to the case with different ranks. Still, the ABJ theory is expected to be integrable (it was proven to be so in a particular sector in the limit of [82]), therefore a derivation of its Bremsstrahlung function from integrability is also foreseeable. This, together with a deeper understanding of the ABJ supersymmetric cusp, would grant a firmer handle on the conjecture for the exact interpolating function of the ABJ model [83]. In conclusion, the lifting to the ABJ theory is challenging and requires more thinking, but is certainly an interesting direction to pursue.
In this paper we have found the exact expression of a cusped 1/2-BPS WL at three loops and at finite N , eq. (6.3). A natural extension of this computation is the evaluation of the same 1/2-BPS cusp, but with open geometric angle (ϕ = 0), along the lines of [51]. This task, if performed using the HQET description, demands the evaluation of the 71 master integrals pointed out in [51], in three dimensions. However, this should be simpler than in four dimensions, as the leading divergence of the cusp integrals at loop l is a pole of order l in three dimensions, instead of 2l in four. Consequently, the relevant integrals should be expanded up to terms of transcendentality 2 in the three loop case (compared to transcendental order 5). We have preliminary evidence that some of the master integrals required for the computation do not evaluate to generalized polylogarithms, rather elliptic sectors appear. Elliptic functions are likely to pop up at higher order in the -expansion or they might cancel out when summing all the diagrams, yet their presence can hinder the evaluation of the master integrals. In particular, at a difference with respect to the four dimensional case, we expect that a canonical form [84] does not exist for all of them.

A Spinor and group conventions
We work in euclidean three dimensional space with coordinates x µ = (x 0 , x 1 , x 2 ). The Dirac matrices satisfying the Clifford algebra {γ µ , γ ν } = 2δ µν I are chosen to be The set of matrices (A.1) satisfies satisfies the following set of identities which allows us to simplify the fermionic contributions to the WL. The relevant traces are instead given by Spinor indices are raised and lowered by means of the −tensor: with ε 12 = −ε 12 = 1. In particular, the antisymmetric combination of two spinors can be reduced to scalar contractions: The gamma matrices with two lower indices, (γ µ ) αβ ≡ (γ µ ) γ α ε βγ , are then given by and they obey the useful identity.
Under complex conjugation the gamma matrices transform as follows: [(γ µ ) β α ] * = (γ µ ) β α ≡ βγ (γ µ ) δ γ αδ . As a consequence, the hermitian conjugate of the vector bilinear can be rewritten as follows where we have taken (χ β ) † =χ β and (ψ α ) † =ψ α . The U(N ) generators are defined as T A = (T 0 , T a ), where T 0 = 1 √ N 1 and T a (a = 1, . . . , N 2 − 1) are an orthonormal set of traceless N × N hermitian matrices. The generators are normalized as The structure constant are then defined by In the paper we shall often use the double notation and the fields will carry two indices in the fundamental representation of the gauge groups. An index in the fundamental representation of U(N 1 ) will be generically by the lowercase roman indices i, j, k, . . . , while for an index in the fundamental representation of U(N 2 ) we shall use the hatted lowercase roman indicesî,ĵ,k, . . .

B Basic facts on ABJM action
Here we collect some basic features on the action for general U(N 1 ) k × U(N 2 ) −k ABJ(M) theories. The gauge sector contains two gauge fields (A µ ) i j and (Â µ )ˆiĵ belonging respectively to the adjoint of U(N 1 ) and U(N 2 ). The matter sector instead consists of the complex fields (C I ) iĵ and (C I )ˆi j as well as the fermions (ψ I ) iĵ and (ψ I )ˆi j . The fields (C I ,ψ I ) transform in the (N 1 ,N 2 ) of the gauge group while the couple (C I , ψ I ) belongs to the representation (N 1 , N 2 ). The additional capitol index I = 1, 2, 3, 4 belongs to the R-symmetry group SU (4). In order to quantize the theory at the perturbative level, we introduce the usual Feynman gauge-fixing for both gauge fields and the two corresponding sets of ghosts (c, c) and (c,ĉ). Then the action contains four different contributions Here the invariant SU (4) tensors IJKL and IJKL are defined by setting 1234 = 1234 = 1. The covariant derivatives used in (B.1) are instead given by The action (B.1), in absence of the gauge-fixing terms and of the corresponding ghosts, it is invariant under the following N = 6 SUSY transformations

C Feynman rules
We use the Fourier transform definition From the action (B.1) we read the following Feynman rules 6 • Vector propagators in Landau gauge • Scalar propagator At one-loop, if we work at finite N 1 and N 2 , beside corrections to the tree level vector propagators, a contribution to the mixed AÂ propagator is generated. They read The one-loop fermion propagator reads and is proportional to the difference (N 1 − N 2 ) of the ranks of the gauge groups. Hence it vanishes in the ABJM limit and is not needed in our computation.

D Two-loop self-energy corrections to the fermions
In this appendix we spell out the computation of the two-loop corrections to the fermion two-point functions. Throughout this section we work with different N 1 , N 2 group ranks.
The relevant Feynman diagrams are depicted in figure 13. Considering all the Summing all the contributions we get The reducible two-loop correction of figure 14 gives Each correction to the fermion propagator can be Fourier transformed back to configuration space

E Useful identities on the cusp
We parametrize a point on the line forming the cusp as Simple identities that turn out to be useful along the calculation are The spinor couplings η,η for the fermions and the matrices M for the scalars appearing in the superconnection (2.3) obey the identities (at ϕ = 0) n 1I n I 1 = n 2I n I 2 = 1 (E.3a) n 1I n I 2 = cos Identity (E.3d), and the formulae above it, are widely used in the computation, since they apply whenever a diagram contains a fermion arch. Identities (E.3g) have been used in the diagrams with scalar insertions of figure 8, while identities (E.3j) are relevant for diagram (3) d .

G Results for two-loop diagrams
Here we provide the list of results for the two-loop diagrams of figure 7, considering only the 1PI configurations on the cusp.

H Results for three-loop diagrams
Here we list the results for the HQET-1PI part of the diagrams. The results include the nonplanar corrections and are already given by the sum over different configurations belonging to the same topology, following the classification of section 5.1. A common factor (e −γ E 16π) k 3 is understood.