Abstract

This review is devoted to collecting some results on the high spin expansion of (minimal) anomalous dimension. Thanks to the recent rationale on integrability, planar super Yang-Mills theory (or its string counterpart) represents a very practicable field. Here the attention will be restricted to its sector of twist operators, although the analysis tools are quite general (in integrable theories). Some structures and ideas turn out to be general also for other sectors or gauge theories.

1. Framework and Beyond

We will move our investigation within the maximally supersymmetric gauge theory in planar limit, that is, for number of colours and coupling , so that the `t Hooft coupling may stay finite. Among the different sectors (perturbatively closed under renormalisation), we also pick up the twist sector, spanned by local composite operators of trace form where is the (light-cone) covariant derivative acting in all the possible ways on the complex bosonic fields . Trace ensures, of course, gauge invariance. The Lorentz spin of these operators is and is the -charge which also coincides with the twist (classical dimension minus the spin). Besides, this sector may be described—thanks to the AdS/CFT correspondence [13]—by spinning folded closed strings on spacetime with and angular momenta and , respectively [4, 5].

As being in a conformal model, suitable superpositions of operators form dilatation operator eigenvectors with definite dimensions (eigenvalues), which are made up of a classical part plus an anomalous one. For instance, in the sector (1.2) this spectral problem shows up dimensions where is the anomalous part. According to the AdS/CFT strong/weak coupling duality, the set of anomalous dimensions of composite operators in SYM coincides with the energy spectrum of the string theory ([15] and references therein), although the perturbative regimes are interchanged. The highly nontrivial problem of evaluating the anomalous part in SYM was greatly simplified by the discovery of integrability in the purely bosonic sector at one loop [6]. Later on, this fact has been extended to all the gauge theory sectors and at all loops in a way which shows up integrability in a weaker sense but still furnishes the investigators many powerful tools [710]. More in detail, any operator (e.g., of the form (1.2)) has been thought of as a state of a “spin chain,” whose Hamiltonian is, of course, the dilatation operator itself, although the latter does not have an explicit expression of the spin chain form, but for the first few loops. Nevertheless, the large size (asymptotic) spectrum has turned out to be exactly described by certain Asymptotic Bethe Ansatz-like equations (the so-called Beisert-Staudacher equations, cf. [711] and references therein). In other words, the anomalous dimensions coincide with the energies given by the Bethe Ansatz solutions (or roots): this is, of course, a great simplification of the initial spectral problem.

Unfortunately, this works only for infinitely long operators: anomalous dimensions of operators with finite length depend not only on Asymptotic Bethe Ansatz (ABA) data but also on finite size “wrapping” corrections: in the perturbative expansion wrapping effects are observed [12, 13] starting from the order . Recent progress [1417] has shown that a set of Thermodynamic Bethe Ansatz (TBA) equations provides a basis for exact (any length at any coupling) predictions for anomalous dimensions of planar SYM.

However, in the sector of SYM relevance of wrapping effects seems to be reduced, even for short operators, if one goes to the high spin limit. For instance, findings of [18, 19] have showed that, at least at twist two and up to five loops, wrapping corrections start contributing at order . This fact has pushed the idea of applying ABA techniques to the study of the high spin limit of twist operators, which—on the other hand—had already received much attention in the past literature. Indeed, the high spin behaviour of the anomalous dimension shows a Sudakov behaviour determined by the so-called universal (since it does not depend on or the flavour) scaling function, [11, 2023]. Actually, this behaviour is more general than in planar  SYM and the adjective scaling is due to the linear value of the cusp anomalous dimension with coefficient while the cusp angle tends to infinity [24]. (Polyakov noticed as first that for cusped Wilson loop vacuum expectation value the charge renormalisation is not enough as in the noncusped case, because of an extra logarithmic divergence due to the high bremsstrahlung at the cusp. He was led to consider cusps because of their importance in the loop dynamics (in Euclidian space-time) [25].) This large angle behaviour is due to the dominance by the lowest twist () in the renormalisation of the vacuum expectation value of a cusped Wilson loop with a very large angle. In the end, for an infinite angle cusp (i.e., with one light-cone segment) equals twice the cusp anomalous dimension of a light-cone Wilson loop [26]. Additionally, the one-loop problem (and thus ) stays exactly the same for twist operators in QCD as long as the partonic helicities are aligned [27, 28]; this fact also has justified partially the great interest on the twist operators (1.2).

In general, the high spin limit of anomalous dimensions of twist operators goes on as a series of logarithmic (inverse) powers that is, looks, at this order, like an expansion in the large “size” . (with , we indicate terms going to zero faster than any inverse powers of ) Recently [11], the leading term (i.e., ) was obtained—in the hypothesis of being wrapping free—from the solution of a linear integral equation directly derived from the ABA via the root density approach. Moreover, was carefully studied and tested both in the weak [11, 21] and strong coupling limit [2935].

The subleading (constant) contribution received also much attention. In the ABA framework, it was shown [36] to come from the solution of a nonlinear integral equation (NLIE). Then, it was obtained [37] starting from a linear integral equation (LIE). Explicit weak and strong coupling expansions performed using the LIE of [37] are present in [38, 39] and agree with some string theory computations [4042]. Importantly, after results of [18, 19] (cf. above) it may be inferred that both and are exactly given by this approach based on the ABA (without wrapping corrections). Besides, both terms fix the coefficients via the reciprocity relation which may still be, consistently with [18, 19], wrapping free.

The latter reasoning would support that even sublogarithmic terms from ABA be exact, but it is even more unclear if and at what extent this might be trusted. In [43] we studied and computed in a systematic way the ABA contribution to them, by using a set of integral equations. Possible wrapping corrections to our results are still to be determined. In this respect, in addition to TBA findings, also results on the string side of the correspondence in the spirit of [4042] could be of fundamental importance. Nevertheless, the definitive proof on the absence of wrapping corrections should only come from exact (maybe TBA) calculations.

Wrapping effects should be negligible if one takes the infinite twist limit and one restricts to the study of the scaling functions , appearing in the expansion (This should be true at least for small values of . But, again, this reasoning is not rigorous.) For this reason, a reasonable amount of activity was also devoted to the study of limit (1.7), using ABA techniques. Results concerning at weak coupling are present in [36]. Strong coupling behaviour of was studied in [4447] by relying on linear integral equations. A detailed study of , , can be found in [43, 48].

In this review we want to give a summary of our activity in the framework of high spin limit of twist operators. We will first give a general description of the method we have used and whose main tools are integral equations derived from the ABA equations. Then, we will briefly report the most important results we obtained.

This review is organised as follows. In Section 2 we review the ABA equations for twist operators. We illustrate the properties of the minimal anomalous dimension operators and explain a technique which allows exact computations of ABA contributions to the anomalous dimension by using the so-called Nonlinear Integral Equation (NLIE). In Section 3 we specialise to the high spin limit and show that, if one neglects terms, asymptotic anomalous dimension can be computed by relying on integral equations. In Section 4 we study the high spin limit at fixed twist. In Section 5 results in the scaling limit (1.7) are discussed.

2. All-Loop ABA and the (N)LIE

Let us recall the Asymptotic Bethe Ansatz equations [711] for the sector: where being the 't Hooft coupling. The so-called dressing factor [11, 4951] is given by the functions being and being the density of the th charge: It is now clear that configurations of Bethe roots, that is, solutions of (2.1), and the corresponding eigenvalues of the energy are related, respectively, to composite operators and their anomalous dimensions in the sector of SYM.

In the sector states of twist are described by an even number of real Bethe roots which satisfy (2.1). Bethe roots localize in an interval of the real line. In addition to Bethe roots, also real “holes” [11, 20, 36, 37, 52] are present (a better explanation of the nature of holes will be given in the following). For any state, two holes reside outside the interval and the remaining holes lie inside this interval. We will indicate with these “internal” holes.

In this paper we will focus on the minimal anomalous dimension state. For such a state the positions of both roots and holes are symmetric with respect to the origin. For what concerns the internal holes, they all concentrate near the origin, with no roots lying in between.

In general, an efficient way to treat states described by solutions of a nonlinear set of Bethe Ansatz equations consists in writing a nonlinear integral equation, which is completely equivalent to them. The nonlinear integral equation is satisfied by the counting function , which in the case (2.1) reads as where It follows from its definition that the counting function , as a function of the real variable , is a monotonously decreasing function and that Therefore, there are real points such that . It is a simple consequence of the definition of that of them coincide with the Bethe roots. The remaining points are called “holes” and their role will be of fundamental importance in what follows. As we anticipated before, for the minimal anomalous dimension state the internal holes concentrate near the origin; that is, their positions are determined by the relations Excited states are obtained by making different choices for the “quantum numbers” . It follows from the structure of (2.9) that, for any state, depend in a nonlinear way on the counting function .

The nonlinear integral equation for is written by using a modification of the standard ideas underlying the procedure concerning the excited state NLIE [5357], this modification being dictated by the physical situation with two “important” external holes. Suppose that in the interval of the real line Bethe roots and holes are present. Then, by using Cauchy theorem we can express a sum over the Bethe roots of an observable as The right-hand side of (2.10) does not depend on as far as no poles of the integrands lie in the region , . In addition, must be kept sufficiently small, in such a way that where belongs to the integration contour of the last three integral terms in (2.10). This condition is assured by the monotonicity of the counting function, that is, in our case , . We apply (2.10) to the sum over the Bethe roots contained in (2.6): What we have obtained is a nonlinear integral equation—for the counting function —which describes—in a way which is alternative to the Bethe Ansatz equations—the minimal anomalous dimension state. Since Bethe roots are localised in an interval of the real axis (i.e., ), NLIE (2.12) is different from the equation introduced in [5357], which contains integrations over the entire real axis. This difference reveals crucial in the specific problem of twist operators in the sector. Indeed, important simplifications to the structure of (2.12) arise in the high spin limit: if we decide to neglect terms going to zero faster than any inverse power of , then it is possible to get rid of the last three nonlinear terms in (2.12). We are now going to show this important property.

Coming back to (2.10), we first use (2.11) in order to replace all the with We obtain In order to evaluate the nonlinear terms in this expression, we first assume that it is possible to exchange the series with the integrals. Then, we use the following formulae: We need to remark that results (2.15) are correct if the above infinite sums make sense either as convergent or as asymptotic series. In our case, we can use (2.15), since the series we will get are asymptotic.

Using (2.15) in the evaluation of the above nonlinear terms, the dependence on cancels out, as it should be, and we are left with Now, we are allowed to choose in such a way that : therefore, since is bounded, all the terms in (2.16) proportional to the sine function are zero. Thus, we are left only with terms containing the cosine function, that is, after summing over , where is the Bernoulli polynomial. Relation (2.17) allows to write the final formula: The following expressions give an idea of the form of the first nonlinear terms appearing in (2.18). Next step is to apply (2.17) to the nonlinear integral terms (the last three ones) contained in (2.12). We have to replace with . From the form of (2.17) we realise that such nonlinear terms are proportional to derivatives When , and for such derivatives On the other hand, for the derivatives of the counting function we can give the estimates Putting all together, we conclude that when the spin the nonlinear integral terms contained in (2.12) are . Therefore, if we decide to neglect terms and to focus only on corrections , we are entitled not to consider all the nonlinear integral terms in (2.12). In this case we are left with and this equation has to be solved together with condition (2.9), which fixes the holes positions in terms of the counting function . It is important to remark that within our approximations nonlinearity with respect to enters only through (2.9).

Equations (2.23) and (2.9) are our starting point for the study of the high spin limit. They will be worked out in next section.

3. Integral Equations for the Logarithmic Terms

Let us consider integral equation (2.23) satisfied by the counting function in the high spin limit. Passing to derivatives, we define We have We now decompose such equation in its one loop (with index ) and higher than one loop (with index ) contributions. We set , , , where and where, after neglecting quantities which are , and satisfy In (3.4) and (3.5) the notation stands for the one loop component of the position of the th internal hole. In order to solve the one loop (3.4), we consider (see 3.52 of [37]) a function , whose Fourier transform reads Function (3.6) satisfies the following important property [37]: which allows to extend to the whole real axis the integrations involving : a look at many loops (3.5) shows that this is just the property we need in order to try to solve it.

It follows from results in [21] and from numerical checks that we can extend to the entire real axis the integrations involving : what we are missing are terms. Putting all these pieces together, the equation satisfied by is It is now convenient to pass to Fourier transforms. We have Therefore, the Fourier transform of (3.8) reads Inserting in (3.10) the expression (3.6) for , introducing the “magic” kernel , defined in [11] as and restricting to , we finally get the integral equation: Since we are interested in the computation of the anomalous dimension, we go back to the general formula (2.18) for the evaluation of observables: If we specialise to the energy (anomalous dimension), we have Again, in the large spin limit, one has (with ) together with the estimate (2.22) for the counting function. It follows that the nonlinear terms contained in the expression (3.13) for the energy are . Therefore, if we neglect terms, we are allowed to work with only the linear expression: which we find convenient to write in terms of the one loop and higher than one loop densities: Extending the domains of integration to the real axis does not give problems, since we are neglecting terms: we get Passing to Fourier transforms we have and, consequently, Comparing (3.20) with (3.10), we see that which extends the Kotikov-Lipatov relation [58] to all the sublogarithmic , , contributions and allows to compute the high spin anomalous dimension from the higher than one loop density.

For computational reasons, it is more convenient to use the function The function (3.22) satisfies the integral equation (for ): where . We remark that such equation depends on the position of the holes (but not on both , ) and that the quantity is also an unknown of the problem and has to be determined by solving (2.9): which introduces nonlinear effects in the equation for the density . Equation (3.23) is exact if we neglect—in the high spin limit—terms which are and is the main integral equation of this paper. (Actually, for twist two there is a simple way to get also the and terms from the first of (3.23): it is sufficient to replace in the argument of the function the quantity , which represent the leading contribution to the position of the external hole at large , with the more accurate estimate [59, 60] . We then get the same result as of [38].) It will be our starting point in order to investigate the high spin limit of (minimal) anomalous dimension of twist operators in the sector.

We will study two cases. First, we will consider the limit: Then, we will focus on the scaling limit (1.7) [20]: The case (3.25) will be reported in Section 4, the case (1.7) in Section 5.

4. High Spin at Fixed Twist

In this section we report our results on the fixed twist case (3.25). We naturally make a move from (3.23), in which the various contributions to the known (forcing) term are separated according to their power of . As a consequence of the structure of the forcing term, the high spin expansion of the function goes on as a series of logarithmic (inverse) powers (1.6): Consequently, the anomalous dimension at high spin follows the same fate: Coming back to , the term proportional to satisfies the BES linear integral equation: whose solution determines the universal scaling function through . The BES equation was introduced in [11, 21] and thoroughly studied in the weak [11, 21] and the strong coupling [2935] limit. The interested reader can refer to these papers for a detailed study of (4.3).

The four subsequent terms—independent of —appear in the linear integral equation for the density which determines the virtual scaling function . This equation was written in [24, equation 4.11], and then it was reobtained in [38] and used there and in our contemporaneous paper [39] (where it appears as equation 3.3). In notations used in this paper it reads () the virtual scaling function being As in the case of , an explicit expression for , interpolating from weak to strong coupling, has not been found yet. What we did is expanding (4.4) in a systematic way for small , thus getting the weak coupling convergent series for (formula of [39]): In addition, in [39] we performed the strong coupling analysis, by means of analytical and numerical computations. For the first leading terms (formula () of [39]) we provided the following outcome: where , , is the Catalan's constant and (see also [38]) Importantly, in this asymptotic expansion we could show that the only trace of the twist comes up in the piece and thus cancels out completely, at order , in the asymptotic (large ) expansion of . It follows that the constant term (i.e., ) in at order is for any twist: this allowed successful comparisons with string theory results [4042], which does not distinguish between null and small values of .

An alternative analysis of the strong coupling limit can be found in [38]. This computation was done adapting the method introduced in [34] for .

For what concerns , that is, the , , contributions to the anomalous dimensions, they are “driven” by the holes depending parts of (3.23). As a consequence of (3.21), one has , where satisfies the integral equation (): with the symbol standing for the component proportional to . The analysis of this case was done in detail in the paper [43]. We report here the main results.

In order to fulfil condition (2.9), the position of the internal holes has to expand in inverse powers of : Introducing the (even) derivatives in zero of the (even) function , developing them in powers of , and imposing the condition (2.9) for the holes, we eventually get the following recursive equation for the unknowns : where the contained in the second term of the right-hand side are constrained by the conditions , and the ones in the third term by , .

The next step is the Neumann expansion for the functions (in the domain ): Straightforward but lengthy calculations, originating from (3.23), lead to the conclusion that the Neumann modes satisfy the system (We use the notation: ) In (4.14) appears as a coefficient in the high spin expansion of the internal holes-depending function appearing in (4.9). As a consequence of (4.12), depends on the various coefficients of (4.10) as In order to study the system (4.14), it is useful to introduce the “reduced coefficients" , defined in and of [47] as solutions of the reduced systems: with the explicit (i.e., not depending on the densities ) forcing terms: Indeed, solutions to the systems (4.14) are linear combinations of the various with coefficients depending on . In particular, for what concerns , its expression in terms of and reads The various are written in terms of the densities with the help of (4.12), leaving eventually as depending on and , with . This last property is very important, since it makes possible to build up a recursive calculation scheme for the , opening the way to push the computation up to the desired order in .

As an example in paper [43] we gave the following exact results for the first : Due to the aforementioned recursive properties, such expressions can be explicitly computed in the weak and in the strong coupling limit. Weak coupling expansions are presented in Appendix A of [43]. The strong coupling leading term is given in Section of [43].

5. High Spin and Large Twist

We now want to study the anomalous dimension in the limit (1.7): Since the number of internal holes becomes infinite, we need to treat the sum over them contained in (3.23) by means of the technique discussed in Section 2. Before doing that, we introduce two points which separate the internal holes and the Bethe roots: since the number of internal holes equals , the “separator” has to satisfy the following relation: Now we can specialise formula (2.18) to the sum over the internal holes . Remember that in no Bethe roots are present: therefore the right-hand side of (2.18) is zero. We get where denotes the th derivative of the density in . This expression has to be inserted in (3.23) and worked out together with condition (5.2): we get in general a nonlinear integral equation for the quantity , related to the Fourier transform of the density of Bethe roots and internal holes through (3.22).

This joined analysis of (3.23) and (5.2) gets simplified in the limit (1.7). Indeed, it follows from the structure of (3.23) that, in the case of limit (1.7), the function expands as [36, 43, 48] and, correspondingly, the anomalous dimension behaves as (1.8): The functions and are called generalised scaling functions. In particular coincides with the universal scaling function .

Similarly, the separator between internal holes and Bethe roots, , enjoys the following scaling in the limit (1.7): The constants are connected to the various components in which the density and its derivatives in zero expand in the limit (1.7): by means of (5.2): For instance, using (5.7) and (5.8) in (5.2), we get, for the first , Therefore, once we take the limit (1.7) and use condition (5.2), we easily get that (3.23) splits in a set of equations, one for every function . Remember that coincides with the BES density (function of last section): therefore, it satisfies the BES linear equation (4.3), that is, Another particular case is the function : this satisfies the linear equation and determines the generalised scaling function . Equation (5.12) was studied in paper [44].

The last case that has to be treated separately is . However, in this case, the density is obtained from ( contribution at fixed twist, see last section) and , by means of the equality In the three cases we have discussed up to now, the dynamics of the internal holes is not relevant: for the required approximations the internal holes can be supposed all lying at the origin; that is, the sum (5.3) can be considered as vanishing.

The particular form of the expansion (5.7) with coefficients given by (5.10) is however of fundamental importance for all the remaining cases. They can be studied together, by means of the integral equation in which the sum over the internal holes is evaluated through (5.4) and the “separator” is expanded as in (5.7), with coefficients depending on the quantities and determined by the use of (5.2).

From the methodological point of view, we need to distinguish two cases. The first one covers the values and . In this case only the first term in the right-hand side of (5.4), which is linear in the density, is relevant. Therefore in the final equation nonlinearity comes only from the nonlinear dependence of on the density and their derivatives in zero: For instance, using (5.10) we get for the first cases In the remaining cases, that is, , the evaluation of the sum over the internal holes (5.3) involves also terms explicitly nonlinear in the density. Thanks to formulae (2.18) and (5.4), however, everything is under control and, for instance, for the first values of we get In both cases, next steps are the usual ones and details can be found in [43, 47, 48]. First, we perform a Neumann expansion for the even functions , in the domain : This implies that the generalised scaling functions are expressed as The Neumann modes satisfy linear systems and are linear combinations of the “reduced” coefficients which are solutions of the systems (4.18). The coefficients driving such linear combinations depend on the densities and their derivatives in zero, , with . This property allows to find, step by step, exact expressions for the generalised scaling functions in terms of and . For the first of them we get the following results: Explicit expressions in the weak and strong coupling limit for the various can be given without much ado, because of the iterative structure of the relations which determine them. They can be found in [43, 47, 48].

In particular, strong coupling limit of the anomalous dimension is of interest since it can be checked against string theory calculations. For what concerns the function performing such a check, is not a difficult task, after results by Alday and Maldacena [23]. Introducing (at large ) the quantity in [23], it was proved that in the limit (1.7), when , , with fixed, the quantity has to coincide with the energy density of the ground state of the nonlinear sigma model with mass gap . When we are in the nonperturbative regime of the nonlinear sigma model. In this case the energy density can be computed by using Bethe Ansatz-related techniques. This computation has been systematically performed in [61]. In order to have agreement between our calculations for and computations of [61] for the nonlinear sigma model, we must have that the quantities computed in that paper have to be related to by the relation . Our results [39] for the strong coupling limit of and are where we used the compact notations and agree with corresponding formulae contained in [61]. (In [39] results for the strong coupling limit of and checks with nonlinear sigma model results were performed up to )

6. Conclusions

In this paper we have discussed our activity on the study of minimal anomalous dimension of twist operators in the sector of SYM at high spin. We preferred to give a general description of the methods we used: the reader can refer to the original papers for details and for a complete list of our results.

The main tool we used is integral equation (3.23) for the density of Bethe roots. This equation coincides with the exact nonlinear integral equation equivalent to the Asymptotic Bethe Ansatz equations of the sector if we neglect terms going to zero faster than any inverse power of the logarithm of the spin. Equation (3.23) is linear, apart from the nonlinear dependence of the internal holes positions on the counting function . These nonlinear effects are taken into account by inverting relation (2.9): in the high spin limit such inversion is feasible because of the recursive properties of (4.12) determining the various . This allows to give explicit expressions (4.20) for the coefficients of the high spin expansion of the anomalous dimension.

When the twist goes to infinity, the internal holes are described by their density and contribute to integral equation (3.23) introducing explicitly nonlinear terms. These nonlinear terms are evaluated (see Section 5) by using properties and techniques of the nonlinear integral equation, which are discussed in Section 2. Their contribution to main equation (3.23) is obtained by applying formula (2.18) and is reported in (5.16) and (5.17). In this case also, recursive properties of the equations determining the high spin limit of the density are crucial in order to obtain explicit expressions for the anomalous dimension.

Acknowledgments

The authors thank D. Bombardelli, G. Infusino, and P. Grinza for collaboration on the various papers reported here. They acknowledge the INFN Grant Iniziative specifiche FI11 and PI14, the Italian University PRIN 2007JHLPEZ “Fisica Statistica dei Sistemi Fortemente Correlati all'Equilibrio e Fuori Equilibrio: Risultati Esatti e Metodi di Teoria dei Campi” for travel financial support.