Table of Contents Author Guidelines Submit a Manuscript
Advances in High Energy Physics
Volume 2010, Article ID 703064, 31 pages
Review Article

Scattering Amplitudes, the AdS/CFT Correspondence, Minimal Surfaces, and Integrability

School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540, USA

Received 22 March 2010; Accepted 8 April 2010

Academic Editor: Carlos Nunez

Copyright © 2010 Luis F. Alday. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


We focus on the computation of scattering amplitudes of planar maximally supersymmetric Yang-Mill in four dimensions at strong coupling by means of the AdS/CFT correspondence and explain how the problem boils down to the computation of minimal surfaces in AdS in the first part of this paper. In the second part of this review we explain how integrability allows to give a solution to the problem in terms of a set of integral equations. The intention of the review is to give a pedagogical, rather than very detailed, exposition.

1. Introduction

The aim of this review is to study gluon scattering amplitudes of four-dimensional planar maximally supersymmetric Yang-Mills (MSYMs). We hope that the study of such amplitudes would teach us something about scattering amplitudes of QCD, but at the same time they are much more tractable. The reason for such tractability is twofold. On one hand, perturbative computations are much simpler than in QCD, due to the high degree of symmetry. In fact enormous progress has been made in the last few years. On the other hand, the strong coupling regime of the theory can be studied by means of the AdS/CFT duality, by studying a weakly coupled string sigma-model.

The aim of these lectures is to explain how to use the AdS/CFT duality in order to compute gluon scattering amplitudes of planar MSYM at strong coupling.

In the next section we briefly describe some weak coupling perturbative results for scattering amplitudes in planar MSYM. Scattering amplitudes of massless gauge theories in four dimensions suffer from IR divergences; hence, a regulator needs to be introduced in order to define them properly. A convenient choice is dimensional regularization, in which the theory is studied in dimensions. As we will see, perturbative computations suggest an ansatz, the so-called BDS ansatz, according to which the amplitude exponentiates. By now, we know that the BDS ansatz is not correct at higher loops and for a large number of particles. Still, it was a useful guidance for recent developments, and perturbative results are often expressed in terms of deviations from this ansatz.

In Section 3 we proceed and explain how the AdS/CFT duality can be used to compute scattering amplitudes at strong coupling. As in the weak coupling computation, amplitudes are IR divergent as a regulator has to be introduced. In order to set up the computation we introduce a -brane as regulator. Actual computations, however, are done by using the analogous of dimensional regularization, since we are interested in comparing our results to the expectations from the perturbative regime. The computation of scattering amplitudes at strong coupling reduces to a minimal area problem in AdS. We show in detail how the prescription works for the scattering of four gluons.

In Section 4 we focus on the mathematical problem of computing the area minimal surfaces in AdS relevant to the study of scattering amplitudes. This lead us to study classical strings on AdS. By means of a reduction, sometimes called Polhmeyer reduction, one can focus on the physical degrees of freedom of the problem. The problem possesses integrability, which in practical terms means that we can introduce a spectral parameter and study how physical quantities depend on such parameter. In particular, this allows to give a solution of the problem in terms of a set of integral equations that have the form of thermodynamic Bethe ansatz (TBA) equations. Furthermore, we argue that the area of the minimal surfaces, directly related to the scattering amplitude at strong coupling, is nothing but the free energy of the above-mentioned TBA system.

In Section 5 we end up with some conclusions and a list of open problems. Finally, we include an appendix with a Mathematica code solving numerically the TBA integral equations for a particular example.

The present paper is not intended to be a complete/fully precise exposition of all the recent results regarding scattering amplitudes in planar maximally supersymmetric Yang-Mills. Rather, we attempt to explain in a pedagogical manner the computation at strong coupling and derive the final solution in the most straightforward possible manner.

2. Perturbative Scattering Amplitudes of Planar MSYM

In this section we briefly discuss the basics regarding the perturbative computation of planar gluon scattering amplitudes of MSYM. For more complete reviews see [1, 2] and references therein.

2.1. Basics

Being massless particles of spin one in four dimensions, gluon states are characterized by their helicity, , four momenta, , with , and color indices in the adjoint representation of the gauge group . Generic scattering amplitudes depend in a complicated manner on these.

In the planar limit it is useful to write the amplitudes in the so-called color decomposed form where we denoted by the -loop, -point amplitude. The sum runs over noncyclic permutations, denotes the number of colors and the gauge theory coupling constant. Furthermore, we denoted by the gauge group generators in the adjoint representation. In the nonplanar case, there are also contributions with multi-traces, suppressed in the large limit. This decomposition separates the color structure from the kinematics. The term proportional to a given color structure is the so-called leading color ordered amplitude , which will only depend on the momenta and the helicities of the particles undergoing the scattering.

As already mentioned, scattering amplitudes of massless particles in four dimensions suffer from infrared divergences. In order to see this more clearly, consider the scattering of four gluons at one loop. The amplitude contains a factor of the form where is the momentum running along the loop and are the momenta of the external particles. We recognize two kinds of divergences. First, from the region , these are called soft divergences, since they are due to the interchange of gluons with very low (soft) momenta among external gluons. The second class of divergences comes from the region and are called collinear divergences, since in this case the momentum of the gluon interchanged is parallel to the momentum of one of the external gluons. Notice that in addition we can have with , so soft and collinear effects can combine, giving an enhanced divergence.

Such amplitude can be regularized in several ways. A popular choice is to use a version of dimensional regularization and consider the theory in dimensions. (Another possibility is to go to the Coulomb phase of the theory, staying in four dimensions [34]. This is equivalent to giving a small mass to the propagators in the outer loop when computing scattering amplitudes.) More precisely, we still consider four-dimensional external states but consider the momentum running inside the loop in a dimension slightly higher than four ( is negative) (In practice, we do the computation in four dimensions and once we arrive to the final result (2.2), we replace .)Once a regulator is introduced, the amplitudes are finite, but the price we have to pay is an explicit dependence on such regulator. (and the consequent breaking of some of the symmetries is present in the original problem.) In the case of dimensional regularization, IR divergences manifest themselves as poles in . From the explicit one loop expression (2.2) and the discussion above, we see that soft and collinear divergences combine to give a divergent factor 1. In the general case of loops, we expect the same phenomenon to occur for each subloop, hence we expect a divergence of the form As already mentioned, color-ordered amplitudes depend also on the helicities of the external gluons. As a consequence of supersymmetry, scattering amplitudes of MSYM satisfy supersymmetric Ward identities. In particular, these identities imply that certain amplitudes vanish at all loops, namely, an amplitude vanish if all helicities, or all but one, are plus The first nontrivial amplitude is the so-called maximally helicity violating (MHV) amplitude, and it has the form , namely, two external states have minus helicity and the rest plus. A great simplicity of studying MHV amplitudes is given by the fact that they contain a single Lorentz structure, which is already captured by the tree level amplitude, to all orders in perturbation theory. It is then convenient to factor out the tree level amplitude and study the so-called reduced amplitude which depends only on the kinematical invariants and the regulator. So, in summary, we consider the planar limit and furthermore we focus on MHV amplitudes only. The reduced amplitude will be simply a scalar function of the kinematical invariants (and of course ).

2.2. The BDS Ansatz

In the last years, a bast technology has been developed for the computation of perturbative scattering amplitudes; see, for instance, [1]. Such technology made it possible for the computation of scattering amplitudes for four (and to some extent for five) gluons to rather high order in the loop expansion. Bern, Dixon, and Smirnov (BDS) noticed in [3], after previous hints in [4], that all known MHV amplitudes satisfied a recursion relation, allowing to write amplitudes at a given order in terms of lower-order amplitudes. These recursions relations could be encoded by the following ansatz, which was then proposed by BDS for MHV amplitudes of any number of gluons at any loop order: where is proportional to the 't Hooft coupling constant and keeps track of the perturbation order. We have also introduced an IR scale , in order to account for the dimensionfullness of the coupling constant in dimensions different from four.

The above relation, from now on called BDS ansatz, expresses the (log of the) amplitude in terms of the one loop amplitude . Note that this relation is highly nontrivial, since the constants and do not depend on the number of gluons or the kinematics. Notice that in addition to the constant piece, there can be noniterating terms of order . Due to the presence of these noniterating terms, it is more proper to say that the BDS ansatz is an ansatz for the log of the amplitude rather than for the amplitude itself, since these terms of order could hit the divergent term in order to give finite contributions when expanding the exponential.

The structure of the IR divergent terms in (2.7) was determined in [5], where it was shown that indeed IR divergences factorize and exponentiate as in (2.7). The highly nontrivial content of the ansatz is that the finite pieces also exponentiate in a similar manner and are given, at all loops, by the same functions that characterize the IR divergent terms (times the finite pieces at one loop).

We will be particularly interested in the scattering of four gluons; in that case the BDS ansatz reduces to where and are the usual Mandelstan variables for the scattering of four particles. The amplitude has two pieces, a divergent one and a finite one, and is characterized by two functions. The so-called cusp anomalous dimension controls both the leading divergent piece and the finite piece while the so-called collinear anomalous dimension, , controls the subleading divergent pole. Note that the BDS ansatz is giving the answer for the scattering amplitudes to all loops simply in terms of two functions of the coupling constant. Of course, one still would need to compute these functions, but that in principle is a much simpler task. In the next chapter we will show how to compute the four-point scattering amplitude at strong coupling and compare the result to the BDS ansatz.

The BDS ansatz satisfies a number of nontrivial tests:(i)as already mentioned, it is consistent with the well-known exponential structure of IR divergences;(ii)it agrees with explicit computations for four gluons at three loops and five gluons at two loops;(iii)it has the correct collinear factorization properties, that is, the correct limit when two adjacent external momenta become collinear.

The last point is nontrivial and very restrictive. Color-ordered amplitudes have a simple universal behavior as the momenta of two adjacent gluons, let us say that and become collinear, (very) schematically: (The correct expression contains also terms involving the point function at lower order.)with , in the collinear limit and is a universal function. The BDS ansatz turns out to be consistent with this behavior.

At the moment, there are strong reasons to believe that the BDS ansatz is correct for four and five gluons, while direct evidence against it was found by an explicit computation for the scattering of six gluons at two loops [6]. Still, even if not totally correct, such ansatz is useful as a guidance in perturbative and strong coupling computations. For instance, it is by now customary to express perturbative answers in terms of deviations from the BDS ansatz. We will further discuss this issue at the end of the next section.

Let us mention that the cusp anomalous dimension plays a very important role in recent developments regarding integrability in the AdS/CFT correspondence. It turns out that this function arises in many computations. In particular it controls the large spin limit of the anomalous dimension of the so-called twist two operators, of the form , namely, for large , to all orders in perturbation theory. Using the integrability of the gauge theory dilatation operator [7] constructed an integral equation whose solution gives the cusp anomalous dimension to all values of the coupling constant. For instance, the equation can be solved at weak and strong values of the coupling [810], giving Such values agree with previous computations, done either at weak coupling by using perturbation theory [11] or at strong coupling by using the AdS/CFT duality [12].

As for the collinear anomalous dimension , note from the expression (2.9) that it is scheme dependent. For instance, if we change , then changes by a term proportional to the cusp anomalous dimension . Hence the cusp anomalous dimension is a more "robust" quantity. This is also related to the fact that is related to the anomalous dimension of an operator, as described above, while is not.

In the next section we will explain how to compute gluon scattering amplitudes of planar MSYM at strong coupling by using the AdS/CFT duality. In particular, for the case of four gluons, the strong coupling results are in perfect agreement with (2.7).

3. Gluon Scattering Amplitudes at Strong Coupling

In order to attack the problem of computing scattering amplitudes at strong coupling (In this section, we follow closely [35], to which we refer the reader for the details.) we will make use of the AdS/CFT duality [13, 14]. This duality expresses the equivalence between four dimensional MSYM and type IIB string theory on . The duality provides us with a dictionary between the parameters on both sides of the correspondence, more precisely where is the string coupling and is the radius of the that coincides with that of the . We see that in the limit of a large number of colors and large , the strings are well described by a weakly coupled sigma model. A powerful feature of the AdS/CFT duality is that it allows to compute certain observables of MSYM at strong coupling by doing geometrical computations on AdS. A well-known example is the computation of the expectation value of supersymmetric Wilson loops, which reduces to a minimal area problem [15, 16]. In this section we will show that this is also the case for the computation of scattering amplitudes at strong coupling.

As in the gauge theory, we will need to introduce a regulator in order to define properly scattering amplitudes. In order to set up our computation we introduce a -brane as IR regulator. Actual computations, however, will be done by using the supergravity analogous of dimensional regularization, since it is simpler to proceed in this scheme, and besides we are interested in comparing our results to gauge theory expectations.

3.1. Set-Up of the Computation

As a first IR regulator we consider a -brane localized in the radial direction. In other words, we start with the metric written in Poincare coordinates and we place a -brane at some fixed large value of and extending along the coordinates. The asymptotic states are open strings that end on that -brane. We then consider the scattering of these open strings, that will have the interpretation of the gluons.

The proper momentum of the strings is , where is the momentum conjugate to , and it plays the role of gauge theory momentum and will be kept fixed as we take away the IR cut-off, . Therefore, due to the warping of the metric, the proper momentum is very large, so we are considering the scattering of strings at fixed angle and with very large momentum.

Amplitudes in such regime were studied in flat space by Gross and Mende [17]. The key feature of their computation is that the amplitude is dominated by a saddle point of the classical action. In our case we need to consider classical strings on AdS. Hence, we need to consider a world-sheet with the topology of a disk with vertex operator insertions on its boundary, which correspond to the external states (see Figure 1). A disk amplitude with a fixed ordering of the open string vertex operators corresponds to a given color-ordered amplitude.

Figure 1: World-sheet corresponding to the scattering of four open strings. In the figure on the left we see four open strings ending on the IR -brane; the world-sheet has then the topology of a disk, shown on the right, with four vertex operator insertions.

What are the boundary conditions for such world-sheet? First of all, in the vicinity of a vertex operator, the momentum of the external state fixes the form of the solution, since far away the solution looks like a free string propagating with a given momentum. Secondly, since the open strings are attached to the -brane, at the boundary.

It turns out to be hard to find classical solutions with the above properties. In order to state more simply the boundary conditions for the world-sheet, it is convenient to describe the solution in terms of -dual coordinates , defined as follows: Note that we do not -dualize along the radial direction. After defining we end up again with formally the same metric

In terms of the (3.4) coordinates, the boundary of the world-sheet is located at . Besides, the boundary conditions for the original coordinates , which are that they carry momentum , translate into the condition that has “winding” . Summarizing, the boundary of the world-sheet is located at and is a particular line constructed as follows (see Figure 3):

(i)for each particle of momentum , draw a segment joining two points separated by ;(ii)concatenate the segments according to the insertions on the disk (corresponding to a particular color ordering);(iii)as gluons are massless, the segments are light-like. Due to momentum conservation, the segments form a closed polygon.

The world-sheet, when expressed in -dual coordinates, will then end in such sequence of light-like segments (see Figure 2) located at .

Figure 2: Comparison of the world-sheet in original and -dual coordinates. The hyperplane on the picture to the right should not be interpreted as a -brane.
Figure 3: Polygon of light-like segments corresponding to the momenta of six external particles undergoing a scattering.

Before -dualizing, we can interpret the regulator -brane as several -branes, one per each external particle, at a fixed radial distance and extended in the four directions . The external states are then open strings ending in these branes, in such a way that two consecutive states end on each brane.

As we perform the four -dualities, these -branes become branes or instantons, and the external states become open strings stretching a light-like distance between two such instantons. These instantons are located at the vertices, or cusps, of the above mentioned polygons.

As we take away the IR cut-off, , the boundary of the world-sheet moves towards the boundary of the -dual metric, at . At leading order in the strong coupling expansion, the computation that we are doing is formally equivalent to the one we would do if we were to compute the expectation value of a Wilson loop given by a sequence of light-like segments.

Our prescription is then that the leading exponential behavior of the -point scattering amplitude is given by the area of the minimal surface that ends on a sequence of light-like segments on the boundary Note that at strong coupling it is more meaningful to talk about the geometric area of the world-sheet. This area is directly related to the log of the amplitude rather than to the amplitude itself. On the other hand, the BDS ansatz (2.7) is also an ansatz for the log of the amplitude. (As already mentioned in the previous section, the presence of noniterating terms in the exponent of (2.7) allows only to write the log of the amplitude, and not the amplitude itself, up to terms of order .) Hence, the strong coupling limit of the BDS ansatz can be directly compared to the geometrical area, times a proportionality factor of the form in front of the whole answer.

We should stress that the strong coupling computation is blind to the type or polarization of the external particles and hence valid for a general scattering, MHV or non-MHV, among gluons or other particles. Such information will contribute to prefactors in (3.5) and will be subleading in a . In particular, such prefactor could even vanish, as it should be the case if we were considering amplitudes of the form . The difference between different scatterings of the same number of particles should be visible if we consider quantum corrections to the classical area. This is a very important, still open, problem.

We have then reduced the problem of computing scattering amplitudes at strong coupling to the problem of finding minimal surfaces in AdS. In the next section we will use the integrability of the problem in order to give a general solution in the form of a set of integral equations. However, in the following we will show that such surface can be found for the particular case of the scattering of four gluons without using any fancy technology. To find and understand this solution in detail will also be instructive for what follows.

3.2. Scattering of Four Gluons

Consider the scattering of two particles into two particles, , and define the usual Mandelstam variables

According to our prescription we need to find the minimal surface ending in the following light-like polygon, see Figure 4.

Figure 4: Polygon corresponding to the scattering of four gluons.

As a warm-up exercise it is instructive to consider the solution close to one of the cusps.

3.2.1. The Single Cusp Solution

We start by considering the solution near the cusp where two of the light-like lines meet. So we consider two semi-infinite light-like lines meeting at a point. This case was already considered in [18]. The surface can be embedded into an subspace of the full

We are interested in computing the surface ending on a light-like Wilson loop along , ; see Figure 5.

Figure 5: Single cusp solution.

The problem has a boost and scaling symmetry that becomes explicit if we choose the following parameterization:

Boosts and scaling transformations are simply shifts of and . The Nambu-Goto action becomes

One can explicitly check that solves the equations of motion and has the correct boundary conditions. Hence the surface is given by

3.2.2. Four Cusps Solution

Now, lets come back to the case of the solution with four cusps. In order to write the Nambu-Goto action it is convenient to use Poincare coordinates , setting , and parameterize the surface by its projection to the plane. In this case we obtain an action for two fields, and , living in the space parameterized by and where , , and so forth. We then need to find a solution to the classical equations of motion with the appropriate boundary conditions. Let us consider first the case with , where the projection of the Wilson lines is a square. By scale invariance we can choose the edges of the square to be at . The boundary conditions are then given by

Near each of the cusps the form of the solution should reduce to the single cusp solution (3.10). Making educated guesses satisfying the boundary conditions and with the correct properties near the cusps we propose

Remarkably it turns out to be a solution of the equations of motion.

On dimensional grounds, the finite part of the amplitude is some function of the form . Hence, in order to capture the nontrivial dependence of the amplitude on the kinematics, we need to consider more general solutions with . In this case the projection of the surface to the plane will not be a square but a rhombus, with and given by the square of the distances between opposite vertices.

In order to find the solution for this more general case, it is instructive to study the surface (3.13) in terms of embedding coordinates. These are coordinates where we view as the following surface embedded in :

The relation between embedding and Poincare coordinates is

The surface (3.13) is then given by

Once we have written our solution in embedding coordinates, we notice that we can apply transformations, that are linearly realized in this coordinates, in order to obtain new solutions. This symmetry is sometimes referred to as “dual conformal symmetry” and should not be confused with the original symmetry associated to the original AdS space. Later we will comment on the implications of this symmetry.

Solutions with can be obtained by starting from (3.13) and performing a boost in the direction Hence, we end up with a two-parameter solution, the parameter , related to the size of the initial square and the parameter , related to the boost parameter. The solution can be conveniently written as where we have written the surface as a solution to the equations of motion in conformal gauge and encode the kinematical information of the scattering as follows:

According to the prescription, we should now plug the classical solution into the classical action in order to obtain the four-point scattering amplitude at strong coupling. However, in doing so, we obtain a divergent answer. That is of course the case, since we have taken the IR regulator away. In order to obtain a finite answer we need to reintroduce a regulator. There are two basic possibilities. (i)The first one, suggested by the strong coupling prescription, would be to set the boundary conditions at some fixed value of the radial direction , namely, in (3.12) , and then take ; however, it seems pretty hard to find classical solutions with such boundary conditions. (Something that one could do is to plug the classical solution into the action and cut the integral giving the area at a finite radial distance, as done in [21], this will give the correct finite piece for the amplitude up to terms independent of the kinematics. However, It will not give the correct collinear anomalous dimension at strong coupling [34].)(ii)The second one is to try to mimic what was done in perturbative computations and consider the amplitudes in dimensions. This route is quite convenient and besides it will allow a more direct comparison with the perturbative results.

3.2.3. Dimensional Regularization

Gauge theory amplitudes are regularized by considering the theory in dimensions. More precisely, one starts with in ten dimensions and then dimensionally reduce to dimensions. For integer this is precisely the low energy theory living on a -brane, where . We regularize the amplitudes at strong coupling by considering the gravity dual of these theories and then analytically continuing in . The string frame metric is Following the steps described above, we are led to the following action: where is the Lagrangian density for . The presence of the factor will have two important effects. On one hand, previously divergent integrals will now converge (if ). On the other hand, the equations of motion will now depend on , and we were not able to compute the full solution for arbitrary . However, we are interested in computing the amplitude only up to finite terms as we take . In that case, it turns out to be sufficient to plug the original solution into the -deformed action (up to a contribution from the regions close to the cusps that add an unimportant additional constant term). After performing the integrals we obtain Expanding in powers of we get the final answer

This should be compared with the field theory expectations (2.7). We notice that the general structure is in perfect agreement with the BDS ansatz. Once we use the strong coupling behavior for the cusp anomalous dimension [12], , we see that the leading divergence, as well as the finite piece, has the correct kinematical dependence with the correct overall coefficient.

3.3. A Conformal Ward Identity

An important ingredient of the previous computation was the existence of a dual symmetry (Actually, this symmetry was first noticed in perturbative computations [36] and then independently in the strong coupling computation described here.), associated to the isometry group of the dual space. This symmetry allowed the construction of new solutions and fixed somehow the finite piece of the scattering amplitude. Naively, this conformal symmetry would imply that the amplitude is independent of and , since they can be sent to arbitrary values by a dual conformal transformation. The whole dependence on and arises due to the necessity of introducing an IR regulator. However after keeping track of the dependence on the IR regulator, the amplitude is still determined by the dual conformal symmetry. Hence, this regulator breaks the dual conformal symmetry, but in a controlled way:

We can write a generic regularized area as follows As already mentioned, the structure of IR divergences of scattering amplitudes is known to all loops. As a result, the divergent piece of the area is well known. Here the denotes the position of the cusps. To a symmetry we associate a Ward identity, that will impose certain constraints on scattering amplitudes. The dual conformal symmetry under consideration will impose the following equation for the finite part of the area [1921] where is the value of the cusp anomalous dimension at strong coupling (properly normalized). This conformal Ward identity (3.27) was first written down in [19]. For the particular case of and , the solution to this equation is actually unique; hence, the finite piece of the scattering of four and five gluons is fixed in case dual conformal symmetry persists to all values of the coupling constant.

What are the consequences of this symmetry at strong coupling? The Ward identity (3.27) implies that the strong coupling answer can be written as a solution to the above equation plus a solution of the homogeneous equation, namely, an arbitrary function of the cross-ratios. On the other hand, the one loop answer multiplied by the appropriate factor is a solution of the conformal Ward identities (as it should be, since the symmetry is present at one loop). In the literature, this term is referred to as the BDS factor, hence where is the so-called remainder function. The BDS ansatz, if true, would imply that is a constant. On the other hand, note that this is indeed the case for , since in this case we cannot construct cross-ratios of the form . On the other hand, we can construct such cross-ratios for , so in this case we can have a nontrivial remainder function. By explicit computations it was actually checked that such remainder function is indeed nontrivial at two loops [6, 2224] and also at strong coupling [25, 26].

4. Minimal Surfaces on AdS

In the previous section we have shown how the problem of computing gluon scattering amplitudes at strong coupling reduces to the computation of the area of certain minimal surfaces in AdS. In this section we show how the integrability of the system can be used in order to give a solution to the problem, in the form of a set of integral equations. We will follow closely [2629], to which we refer the readers for the details.

4.1. Mathematical Problem

The mathematical problem is to find the area of the minimal surface ending on the boundary of AdS at a polygon of light-like edges. Such polygon is parameterized by the position of its cusps , which satisfy . See Figure 6

Figure 6: The light-like polygon in which the minimal surface ends is parameterized by the position of the cusps, , with .

As will be clear below, we will focus on certain regularized area that is invariant under conformal transformations. As such, it will depend only on cross-ratios, of the form . The mathematical problem is then to compute the minimal area as a function of such cross-ratios. As already mentioned, the full problem involves minimal surfaces on ; in which case there are cross-ratios if we are dealing with the scattering of gluons. However, for most of the exposition presented in this paper, we will restrict to special kinematical configurations in which the minimal surfaces involved are embedded in ; in this case, we have cross-ratios (For the general scattering in four dimensions we have coordinates, minus , since the distance between consecutive points has to be light-like, minus 15, that is the dimension of the conformal group . In the case of , we have minus , which is the dimension of .), and the polygon is a zig-zaged polygon living in one-plus-one dimensions, see Figure 7

Figure 7: A zig-zaged Wilson loop in dimensions is parameterized by coordinates and coordinates. If you want a closed polygon, you can fold the figure in a cylinder.

If we wanted a closed polygon we could consider the above contour on a cylinder. Note that we can consider only polygons with an even number of sides, hence . As one can see in Figure 7, the contour is parameterized by coordinates and coordinates. With each set we can form invariant cross-ratios, of the form , and all together we have the expected cross-ratios.

In order to consider minimal surfaces in we need to consider the world-sheet of classical strings on . This is the subject of the following subsection.

4.2. Strings on AdS3

Classical strings on can be described in terms of embedding coordinates, where we see as the following surface embedded in : We take the world-sheet to be the whole complex plane and . Since we are interested in classical solutions, the fields have to satisfy the conformal gauge equations of motion Furthermore, the equations of motion are supplemented by the Virasoro constraints An efficient way to focus only on the physical degrees of freedom is by performing a so-called Pohlmeyer kind of reduction; see, for instance, [3032] and consider the "reduced" fields As a consequence of the equations of motion and Virasoro constraints, can be seen to be a holomorphic function, and and can be seen to satisfy a generalized version of the Sinh-Gordon equation Hence, we have reduced the initial set of equations to a holomorphic function and to a field which satisfies a generalized Sinh-Gordon equation. Eventually, we will be interested in the area of the classical world-sheet. This can be simply written as the conformal gauge action in terms of the reduced fields An important property of the reduced fields is that they are invariant under space-time conformal transformation. In other words, they contain only the “nontrivial” information about the space-time solution.

Summarizing, a solution of classical strings on is mapped to a holomorphic function and a field satisfying (4.5). We can also wonder whether given a holomorphic function and a field satisfying (4.5) it is possible to construct a solution of classical strings on . The answer is affirmative, and there is a well-know procedure to construct such solution. What we need to do is to solve two auxiliary linear problems, which we denote as left and right where the flat connections are simple two-by-two matrices constructed from and , for instance we will denote the different components of the connections by and . On the other hand, the indices and in (4.7) denote independent solutions of the auxiliary linear problems. Each or is then a doublet. We denote the components of this doublet by and so forth.

Given then the solutions of these two auxiliary linear problems, one can show that the space-time solution is simply given by for some constant matrix (which can be taken, for instance, to be the identity matrix). One can show that constructed this way satisfies all the required properties. In order to write down the answer, we have written the space-time coordinates in an notation.

An important observation will be very relevant for what we will do next: it is possible to introduce a spectral parameter in such a way that the flat connection is still flat, more precisely One can check that the connections are still flat, namely, they satisfy for any value of . Another important fact is that the left and right connections are just particular cases of the above flat connection, more precisely

Before proceeding, let us make the following remark. Since is a holomorphic function, it is possible to make a change of coordinates from the -plane to the -plane, where . In the -plane, the generalized Sinh-Gordon equation takes the usual form (after a simple field redefinition of ) Naively, it would seem that we got rid of all the information on . However, this is not the case, since the -plane will have in general a complicated structure (for instance, it will have branch cuts, etc., depending on ). So, we can choose between a complicated equation on the complex plane or a simple equation on a more complicated space. Depending on which questions we want to answer, one description may be more convenient than another. In the -plane, the area is simply given by

4.3. Classical Solutions Corresponding to Scattering Amplitudes

A natural question to ask is what are the properties of the holomorphic function and for solutions corresponding to scattering amplitudes. In order to gain some intuition we can consider the four cusps solution found in Section 3 and perform the Pohlmeyer reduction. We find Hence, the four cusps solution simply corresponds to the vacuum solution of the Sinh-Gordon equation. What about solutions with a higher number of cusps? First of all we propose that the field is regular everywhere, since we are looking for smooth space-like solutions. Second, we expect a general solution to be similar to the four cusps solution when approaching the boundary, which corresponds to very large . So we expect that as becomes large. Finally, if we are interested in a minimal surface ending on a polygon with cusps, we propose to be a polynomial of degree : we have used rescalings and translations in order to set the coefficients of and to one and zero respectively. We see that such polynomial contains complex coefficients, or real coefficients, which exactly agrees with the amount of expected independent cross-ratios. In conclusion, minimal surfaces ending on a light-like polygon with cusps correspond to a holomorphic polynomial of degree and a field satisfying the Sinh-Gordon equations, with boundary conditions such that it decays at infinity and diverges logarithmically at the zeroes of .

Note that since decays at infinity, the integral defining the area (4.14) diverges. We define a regularized/finite area by subtracting the behavior at infinity The computation of this regularized area will be the main focus of the remaining of this paper. As already mentioned, all the fields entering into the expression for the regularized area are invariant under space-time conformal transformations. Hence, the regularized area will be a function of the cross-ratios only. (The full answer would include also the integral of the one we have subtracted. In order to compute it one would need to introduce a physical regulator and this part of the answer will not be conformal invariant. Anyway, its explicitly form can be worked out and turns out to be quite universal. In this paper we will focus on the “interesting” part of the answer .)

4.3.1. Behavior of the Solution at Infinity

In the following, we want to understand the behavior of the solutions of the linear auxiliary problems for very large values of , or , since this will tell us the behavior of the world-sheet near its boundary. Suppose, by simplicity, that , hence . As a result, as we go once around the -plane, we go around the -plane times.

Since for large values of vanishes, the flat connections drastically are simplified and we can solve for the auxiliary linear problems. A general solution will be of the form We see that the -plane divides naturally into quadrants; see Figure 8. In each quadrant one of the two solutions of each problem (left and right) dominates. For instance, if we focus on the upper right quadrant, the solution proportional to dominates in the left problem while the solution proportional to dominates in the right problem. For instance, for the first quadrant we have at large values of . This means that for large values of , the whole quadrant corresponds to a single point in the boundary, which corresponds to a cusp.

Figure 8: When looking at the left problem, (each sheet of) the -plane is naturally divided into two parts, according to the sign of . In the same way, when looking at the right problem, the -plane is naturally divided into two parts, according to the sign of . Hence, the -plane is naturally divided into four quadrants. Large values of in each of these angular sectors correspond to a cusp.

As we change quadrant, one (and only one) of the two dominant solution changes and we jump a light-like distance to the next cusp. Note that in each cusp we can write where is given by the leading contribution to the left problem and by the leading contribution to the right problem. As we change quadrant, one of the two solutions, or , changes. As we go around the -plane times, we get the expected cusps.

The above solution, in which is a homogeneous polynomial, turns out to correspond to a regular polygon with an additional symmetry. In the general case in which the polynomial is not homogeneous, we get a polygon of different shapes but with the same number of cusps. Hence the degree of the polynomial determines the number of cusps while the coefficients on the polynomial determine the shape of the polygon.

4.4. -System for Minimal Surfaces

Let us focus on the left problem. We see that each sheet on the -plane is naturally divided into two sectors, one with and the other with . Note the following simple but important fact: in each sector the small solution is well defined (up to a normalization constant), on the other hand, the large solution is not, as we can add to it a part of the small solution. Let us then introduce the following terminology.(i)The -plane is divided into sectors, since each sheet contains two sectors. We label this sectors by .(ii)We call the small solution at the th sector. By definition, this is the solution with the fastest decay along the line in the center of the th sector, for large . For instance, given solution (4.18), we see that the term proportional to has the fastest decay as we move along the positive real axis in the -plane. Such small solution is well defined, up to a normalization constant.

In order to understand why these small solutions are important, we need to introduce a new element. Given that our connections are matrices, we can introduce an invariant product Given a solution to the linear problem, also called a flat section, one can explicitly check that the above product is independent of the world-sheet coordinate . Hence, we can normalize our solutions such that As already seen, the location of the cusps is determined by the large solutions. The large component of a solution, on a given sector, can be extracted by using the small solution on such sector and the invariant product just introduced, more precisely How do we construct space-time cross-ratios? We have seen that the location of the cusps is given by . The space time cross-ratios involve distances like When seeing as a two-by-two matrix, the normalization condition (4.21) simply means that the determinant of such matrix is one. Given this, one can easily show that Which means that space-time cross-ratios can be constructed from inner products of the small solutions in the corresponding sectors, more precisely We mentioned that small solutions are defined up to a normalization constant. Note that such normalization constant would cancel out when computing cross-ratios.

The strategy we will follow is to introduce the spectral parameter as shown in the previous section and to study the small flat sections of the corresponding connection then, we can consider the cross-ratios as a function of such spectral parameter Then, if we set we obtain the physical plus/left cross-ratios while setting we obtain the physical minus/right cross-ratios: A property of the flat connection that will be used in what follows is that it possesses a symmetry where is the usual Pauli matrix. This symmetry allows to relate small sections at different values of the spectral parameter, for instance, , and in particular, it implies Besides, in order to simplify the subsequent derivation, we will assume that .

Now we have all the elements to derive the so-called Hirota equations and subsequently the -system equations. The trick is to choose and as a complete basis of small flat sections and express two arbitrary consecutive flat sections and in terms of them: Next, we use (4.30) in order to express every wedge as a wedge involving and consider ; then we obtain where we have introduced the notations , , and so forth. This equation can be written in a more suggestive way by introducing . We obtain which has the form of the so-called Hirota equations. From the definition of , we see that it is nontrivial for . The -system equations can be simply obtained by introducing ; we then obtain is nontrivial for . Notice that this agrees with the amount of (complex) cross-ratios of our scattering problem. These are the -system equations. Note that they have simply followed from a chain of rather trivial facts. These are functional equations for and are valid for any value of . One could reintroduce the normalized factors and check that the -functions are given by the usual cross-ratios introduced above. The physical cross-ratios are then obtained by evaluating at and .

Such equations were obtained by rather trivial considerations, and we do not expect them to be the whole story. In particular, note that the information of the holomorphic function does not enter in such equations. There are many solutions to such equations. The correct solution is then picked by specifying the analytic properties and boundary conditions of as we move on the -plane. This will be the subject of the following section.

4.5. Integral Equations

In order to pick the appropriate solution to the -system equations (4.34) we need to specify the analytic properties of . By analyzing the auxiliary linear problem and the definition of one can show that(i) are analytically away from ,(ii)as , the flat connection is simplified and the inverse problem can be solved by using a WKB approximation, where the role of is played by or .

More precisely, by calling , one can show that for large where is given by the periods of along the cycles , more precisely . This is how the information of the polynomial enters into the problem. These periods are usually complex, and there are of them, which exactly agrees with the quantity of expected cross-ratios.

The strategy now is well known from the study of integrable systems. We can combine the -system equations with the analytic properties and boundary conditions for the -functions, in order to write a system of integral equations for them, whose solutions will automatically satisfy the -system equations and have the required boundary conditions. The system of integral equations is given by where the convolution operation is defined by The system of equations (4.36) has the form of TBA equations, that arise when studying integrable models in finite volume; see, for example, [33]. Even though, for the sake of clarity, some coefficients have been suppressed in the derivation of these equations, the equations given have all the correct coefficients. Note that from the TBA point of view, the parameters enter as masses. Once the masses are given, the solution of the above system is unique. The physical cross-ratios can be read off by looking at for appropriate values of .

For simplicity, in (4.36) we have assumed that all the masses are real and positive, and the same will be assumed for the rest of the paper. This can be achieved if all the zeroes of are along the real line, such that is positive for real, positive, and sufficiently large.

The last question is then how we compute the regularized area once we solve the above system of integral equations. It turns out that the area can be written in terms of the -functions and has the simple expression This expression has exactly the form of the free energy of the TBA system.

The strategy to solve the full problem is then clear. The initial data (rather than the coefficients of the polynomial) are the masses . For a given choice of these masses, we find a unique solution of the TBA integral equations. From this solution we can read off the space-time cross ratios and we can also compute the regularized area as the free energy of the TBA system. Hence, we obtain the regularized area for given space-time cross-ratios.

The above equations can be easily solved numerically, by an iterative method. See the appendix for more details.

4.6. Generalization to AdS5

In the first part of this section we have seen how a solution to the problem of computing certain minimal surfaces in can be given by using the integrability of the system. In the following, we highlight the results for the general case of without deriving them. It turns out that, up to some subtleties, much of the previous analysis goes through, and the final results are a simple generalization of the results obtained for the case. We get a slightly more complicated -system for functions , with and . The system reads As before, in order to get a system of integral equations, we need to supplement these equations by the analytic properties of . As for the case, the large behavior can be understood by a WKB analysis. The boundary conditions for the -functions are depending on the signature of space-time, the reality conditions of and are fixed. For instance, for or signature, the are real while the are purely imaginary. As for the case, these boundary conditions together with regularity imply the following set of integral equations: where and the kernels are given by The data of the system of integral equations are the (in general complex) masses and the purely imaginary constants . Since for the scattering of gluons , we get degrees of freedom, which exactly agrees with the expected amount of cross ratios.

Again, as for the case of the equations, these TBA equations can be solved iteratively. Once the solution is found, the space-time cross-ratios can be read off from the solution at particular values of . Furthermore, quite remarkably, the expression for the regularized area exactly agrees with the expression for the free energy of such TBA system This is the final solution for the general problem in (where for simplicity we have assumed real masses).

5. Conclusions

In the first part of this paper we have described how the computation of scattering amplitudes of planar MSYM at strong coupling reduces to a minimal area problem in . More precisely, one is to look for the area of minimal surfaces that is at the boundary on a light-like polygon whose edges are given by the momenta of the particles undergoing the scattering.

In the second part of this paper, we have focused on the above minimal area problem. The integrability of the system plays a crucial role in finding such solution. We have seen that it is convenient to introduce a one-parameter deformation (the spectral parameter ) and study such deformed problem. One can then write down a system of functional equations, or -system, valid for any value of . One can combine these functional equations with the knowledge of the analytic behavior of the -functions in the -plane, in order to write a set of integral equations which can be solved iteratively and give the desired answer. There are many directions one could try to follow.(i)One possible direction would be to take the continuous limit of our construction, in order to study smooth Wilson loops.(ii)It would be nice to find a physical connection between the integrable system that the TBA equations describe and the original integrable system.(iii)Of course, the most interesting problem is the extension of this construction to the full quantum theory.(iv)An interesting direction would also be to look for similar structures in perturbative computations.


Numeric Implementation of TBA Equations for AdS3

In this appendix we describe how to solve the TBA equations (4.36) numerically in Mathematica in a very simple way. We encourage the reader to try the code below by him/herself and to extend it to other cases, in order to understand how much better integral equations are compared to differential equations.

The algorithm is almost trivial: we simply iterate the integral equation plugging the -functions at iteration in the right-hand side of (4.36) and reading from the left-hand side their values at the th iteration, starting with the inhomogeneous terms in the TBA equations as the initial value (or zero in the previous step). We start by defining the kernel appearing in the integral equations, K[x_]=1/(2 Pi Cosh[x]),

and specify how many gluons we want to consider. We do that by introducing a list of masses which appear in the integral equations. For example, let us focus on a polygon with sides, so , and set the following numerical values of the masses: m=1.,2.;M=Length[m].

We also introduce a cut-off for the several integrals and set the number of iterations. (The cut-off is chosen so that for the smallest mass we have .)cut = ArcCosh[8Log[]/Min[m]]; ni = 8.

At each iteration step we compute the new values of the -functions at a discrete set of points and construct the new -function as the function which interpolates through these points. For that we need an interpolating function and a function to perform the several integrals (The last command in the integration function ensures that we are really integrating something which is numerical once the integration variable takes some random numerical value. We chose but anything would work.)F[S_]:=FunctionInterpolation[S,x,-cut,cut,InterpolationPoints->100]int[S_]:= NIntegrate[S,y,-cut,cut]/;NumberQ[S/.y-.1].

This is all the machinery we need. The iterations are simply implemented by Y[a_,k_]: (0#&)/;a0aM+1,Y[i_,1]: Exp[m[[i]]Cosh[#]]&Y[i_,k_]: Y[i,k] F[Exp[m[[i]]Cosh[x]int[K[xy]Log[(1Y[i1,k1],[y])(1Y[i1,k1],[y])]]]].

In the first line we set ; in the second line we set the seed values for the iteration procedure to be simply the asymptotic expressions for the -functions and in the last command we compute the -functions at the iteration step using the right-hand side of the integral equations with the -functions at step . We can now compute the regularized area from the free energy (4.38) of the TBA system and see how it converges to a particular value as we iterate (which gives a good idea of how good the numerics are) energyTable[int[Sum[Cosh[y]m[[i]]Log[1+Y[i,s][y]]/Pi/2,i,M]],s,ni] ListPlot[energy].

For the masses chosen above we see that the free energy converges to . Finally, if we want to read the cross-ratios of the polygon whose area we have just computed, we simply need to evaluate the -functions at , Table[Y[i,ni][],i,M].

To compute the the free energy again for a different choice of masses we should first run Clear[Y].


The authors would like to thank Davide Gaiotto, Amit Sever, Pedro Vieira, and specially Juan Maldacena for collaboration on the material exposed in these lectures notes. This work was supported by U.S. Department of Energy Grant no. DE-FG02-90ER40542.


  1. L. J. Dixon, “Gluon scattering in N=4 Super-Yang-Mills theory: from weak to strong coupling,” in Proceedings of the 8th International Symposium on Radiative Corrections (RAD COR '07), Florence, Italy, October 2007.
  2. L. F. Alday and R. Roiban, “Scattering amplitudes, Wilson loops and the string/gauge theory correspondence,” Physics Reports, vol. 468, no. 5, pp. 153–211, 2008. View at Publisher · View at Google Scholar · View at MathSciNet
  3. Z. Bern, L. J. Dixon, and V. A. Smirnov, “Iteration of planar amplitudes in maximally supersymmetric Yang-Mills theory at three loops and beyond,” Physical Review D, vol. 72, no. 8, Article ID 085001, p. 27, 2005. View at Publisher · View at Google Scholar · View at MathSciNet
  4. C. Anastasiou, L. Dixon, Z. Bern, and D. A. Kosower, “Planar amplitudes in maximally supersymmetric Yang-Mills theory,” Physical Review Letters, vol. 91, no. 25, Article ID 251602, p. 4, 2003. View at Publisher · View at Google Scholar · View at MathSciNet
  5. G. Sterman and M. E. Tejeda-Yeomans, “Multi-loop amplitudes and resummation,” Physics Letters B, vol. 552, no. 1, pp. 48–59, 2003. View at Publisher · View at Google Scholar
  6. Z. Bern, L. J. Dixon, D. A. Kosower et al., “Two-loop six-gluon maximally helicity violating amplitude in maximally supersymmetric Yang-Mills theory,” Physical Review D, vol. 78, no. 4, Article ID 045007, p. 25, 2008. View at Google Scholar · View at MathSciNet
  7. N. Beisert, B. Eden, and M. Staudacher, “Transcendentality and crossing,” Journal of Statistical Mechanics, vol. 2007, article P01021, 2007. View at Publisher · View at Google Scholar
  8. M. K. Benna, S. Benvenuti, I. R. Klebanov, and A. Scardicchio, “A test of the AdS/CFT correspondence using high-spin operators,” Physical Review Letters, vol. 98, no. 13, Article ID 131603, 4 pages, 2007. View at Publisher · View at Google Scholar
  9. L. F. Alday, G. Arutyunov, B. Eden, M. K. Benna, and I. R. Klebanov, “On the strong coupling scaling dimension of high spin operators,” Journal of High Energy Physics, no. 4, article 082, 2007. View at Publisher · View at Google Scholar · View at MathSciNet
  10. B. Basso, G. P. Korchemsky, and J. Kotanski, “Cusp anomalous dimension in maximally supersymmetric Yang-Mills theory at strong coupling,” Physical Review Letters, vol. 100, no. 9, Article ID 091601, p. 4, 2008. View at Publisher · View at Google Scholar · View at MathSciNet
  11. D. J. Gross and F. Wilczek, “Asymptotically free gauge theories. 2,” Physical Review D, vol. 9, p. 980, 1974. View at Google Scholar
  12. S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, “A semi-classical limit of the gauge/string correspondence,” Nuclear Physics. B, vol. 636, no. 1-2, pp. 99–114, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  13. J. Maldacena, “The large N limit of superconformal field theories and supergravity,” Advances in Theoretical and Mathematical Physics, vol. 2, no. 2, pp. 231–252, 1998. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  14. J. Maldacena, “The large N limit of superconformal field theories and supergravity,” International Journal of Theoretical Physics, vol. 38, p. 1113, 1999. View at Google Scholar
  15. S.-J. Rey and J.-T. Yee, “Macroscopic strings as heavy quarks: large-N gauge theory and anti-de Sitter supergravity,” The European Physical Journal C., vol. 22, no. 2, pp. 379–394, 2001. View at Publisher · View at Google Scholar · View at MathSciNet
  16. J. Maldacena, “Wilson loops in large N field theories,” Physical Review Letters, vol. 80, no. 22, pp. 4859–4862, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  17. D. J. Gross and P. F. Mende, “String theory beyond the Planck scale,” Nuclear Physics. B, vol. 303, no. 3, pp. 407–454, 1988. View at Publisher · View at Google Scholar · View at MathSciNet
  18. M. Kruczenski, “A note on twist two operators in N=4 SYM and Wilson loops in Minkowski signature,” Journal of High Energy Physics, no. 12, article 024, 2002. View at Publisher · View at Google Scholar · View at MathSciNet
  19. J. M. Drummond, J. Henn, G. P. Korchemsky, and E. Sokatchev, “Conformal Ward identities for Wilson loops and a test of the duality with gluon amplitudes,” Nuclear Physics B, vol. 826, no. 1-2, pp. 337–364, 2009. View at Publisher · View at Google Scholar
  20. Z. Komargodski, “On collinear factorization of Wilson loops and MHV amplitudes in N=4 SYM,” Journal of High Energy Physics, no. 5, article 019, 2008. View at Publisher · View at Google Scholar · View at MathSciNet
  21. L. F. Alday, “Lectures on scattering amplitudes via AdS/CFT,” Fortschritte der Physik, vol. 56, no. 7–9, pp. 816–823, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  22. J. M. Drummond, J. Henn, G. P. Korchemsky, and E. Sokatchev, “Hexagon Wilson loop = six-gluon MHV amplitude,” Nuclear Physics B, vol. 815, no. 1-2, pp. 142–173, 2009. View at Publisher · View at Google Scholar · View at MathSciNet
  23. C. Anastasiou, A. Brandhuber, P. Heslop, V. V. Khoze, B. Spence, and G. Travaglini, “Two-loop polygon Wilson loops in N=4 SYM,” Journal of High Energy Physics, no. 5, article 115, 2009. View at Publisher · View at Google Scholar · View at MathSciNet
  24. V. Del Duca, C. Duhr, and V. A. Smirnov, “An analytic result for the two-loop hexagon Wilson loop in N=4 SYM,” Journal of High Energy Physics, vol. 2010, no. 3, pp. 1–17, 2010. View at Publisher · View at Google Scholar
  25. L. F. Alday and J. Maldacena, “Comments on gluon scattering amplitudes via AdS/CFT,” Journal of High Energy Physics, no. 11, article 068, 2007. View at Publisher · View at Google Scholar · View at MathSciNet
  26. L. F. Alday, D. Gaiotto, and J. Maldacena, “Thermodynamic bubble ansatz,”
  27. L. F. Alday and J. Maldacena,
  28. L. F. Alday and J. Maldacena, “Null polygonal Wilson loops and minimal surfaces in Anti-de-Sitter space,” Journal of High Energy Physics, vol. 2009, no. 11, article 082, 2009. View at Publisher · View at Google Scholar
  29. L. F. Alday, J. Maldacena, A. Sever, and P. Vieira, “Y-system for scattering ampli-tudes,”
  30. H. J. de Vega and N. Sanchez, “Exact integrability of strings in D-dimensional de Sitter spacetime,” Physical Review D, vol. 47, no. 8, pp. 3394–3404, 1993. View at Publisher · View at Google Scholar · View at MathSciNet
  31. A. Jevicki, K. Jin, C. Kalousios, and A. Volovich, “Generating AdS string solutions,” Journal of High Energy Physics, no. 3, article 032, 2008. View at Publisher · View at Google Scholar · View at MathSciNet
  32. H. Dorn, “Some comments on spacelike minimal surfaces with null polygonal boundaries in AdSm,” Journal of High Energy Physics, vol. 2010, no. 2, article 013, pp. 1–14, 2010. View at Publisher · View at Google Scholar
  33. A. B. Zamolodchikov, “Thermodynamic Bethe ansatz in relativistic models. Scaling 3-state Potts and Lee-Yang models,” Nuclear Physics B, vol. 342, no. 3, pp. 695–720, 1990. View at Publisher · View at Google Scholar · View at MathSciNet
  34. L. F. Alday, J. M. Henn, J. Plefka, and T. Schuster, “Scattering into the fifth dimension of N=4 super Yang-Mills,” Journal of High Energy Physics. In press,
  35. L. F. Alday and J. M. Maldacena, “Gluon scattering amplitudes at strong coupling,” Journal of High Energy Physics, no. 6, article 064, 2007. View at Publisher · View at Google Scholar · View at MathSciNet
  36. J. M. Drummond, J. Henn, V. A. Smirnov, and E. Sokatchev, “Magic identities for conformal four-point integrals,” Journal of High Energy Physics, no. 1, article 064, 2007. View at Publisher · View at Google Scholar · View at MathSciNet