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

Authors: Luis F. Alday
Journal: Advances in High Energy Physics
 (2010)

---

## Abstract

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.


---

## Body

## 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 inD=4-2ϵ 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 Section3 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 D-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 Section4 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 Section5 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.Disclaimer
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|𝒢〉=|h,pμ,a〉 are characterized by their helicity, h=±1, four momenta, pμ, with p2=0, and color indices a in the adjoint representation of the gauge group SU(Nc). 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(2.1)𝒜N(L)~gN-2(g2Nc)L∑ρTr(Taρ(1)⋯Taρ(N))AN(L)(ρ(1),…,ρ(N)),
where we denoted by 𝒜N(L) the L-loop, N-point amplitude. The sum runs over noncyclic permutations, Nc denotes the number of colors and g the gauge theory coupling constant. Furthermore, we denoted by Ta the gauge group generators in the adjoint representation. In the nonplanar case, there are also contributions with multi-traces, suppressed in the large Nc 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 AN(L), 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(2.2)A4,D=4(1)~∫d4pp2(p-k1)2(p-k1-k2)2(p+k4)2,
where p is the momentum running along the loop and ki2=0 are the momenta of the external particles. We recognize two kinds of divergences. First, from the region pμ~0, 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 pμ~αkiμ 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 pμ~αkiμ with α~0, 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 inD=4-2ϵ 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) (2.3)A4,D=4(1)→A4,D=4-2ϵ(1).(In practice, we do the computation in four dimensions and once we arrive to the final result (2.2), we replace d4p→d4-2ϵp.)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/ϵ2. In the general case of L loops, we expect the same phenomenon to occur for each subloop, hence we expect a divergence of the form(2.4)AN(L)~1ϵ2L+⋯.
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(2.5)𝒜(++⋯+)=𝒜(-+⋯+)=0.
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(2.6)MN(L)(ϵ)=AN(L)(ϵ)AN(0),
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:(2.7)ℳ≡1+∑ℓ=1αℓMn(ℓ)=exp[∑ℓ=1(f(ℓ)(ϵ)Mn1(ℓϵ)+Cℓ+𝒪(ϵ))],
where(2.8)f(ℓ)(ϵ)=f0(ℓ)+ϵf1(ℓ)+ϵ2f2(ℓ),α≈λμ2ϵα 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 amplitudeMn1(ℓϵ). Note that this relation is highly nontrivial, since the constants fi(ℓ) and Cℓ 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(2.9)A4=Atree(Adiv,s)2(Adiv,t)2exp(f(λ)8(logst)2+const),(2.10)Adivs=exp(-18ϵ2f(-2)(λμ2ϵsϵ)-14ϵg(-1)(λμ2ϵsϵ)),
where s and t 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 f(λ)=(λ∂λ)2f(-2)(λ) controls both the leading divergent piece and the finite piece while the so-called collinear anomalous dimension, g(λ), 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 amplitudesAN(k1,k2,…,kN) have a simple universal behavior as the momenta of two adjacent gluons, let us say that k1 and k2 become collinear, (very) schematically: (2.11)AN(L)(k1,k2,k3,…,kN)→R(z)AN-1(L)(kt,k3,…,kN)(The correct expression contains also terms involving the N-1 point function at lower order.)with kt=k1+k2, k1→zkt in the collinear limit and R(z) 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 dimensionf(λ) 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 TrϕDSϕ, namely,(2.12)Δ-S=f(λ)logS+⋯
for large S, 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 [8–10], giving(2.13)f(λ)=λ2π2+⋯,λ≪1,f(λ)=λπ+⋯,λ≫1.
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 dimensiong(λ), note from the expression (2.9) that it is scheme dependent. For instance, if we change μ→αμ, then g(λ) changes by a term proportional to the cusp anomalous dimension f(λ). Hence the cusp anomalous dimension is a more "robust" quantity. This is also related to the fact that f(λ) is related to the anomalous dimension of an operator, as described above, while g(λ) 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).


## 2.1. Basics

Being massless particles of spin one in four dimensions, gluon states|𝒢〉=|h,pμ,a〉 are characterized by their helicity, h=±1, four momenta, pμ, with p2=0, and color indices a in the adjoint representation of the gauge group SU(Nc). 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(2.1)𝒜N(L)~gN-2(g2Nc)L∑ρTr(Taρ(1)⋯Taρ(N))AN(L)(ρ(1),…,ρ(N)),
where we denoted by 𝒜N(L) the L-loop, N-point amplitude. The sum runs over noncyclic permutations, Nc denotes the number of colors and g the gauge theory coupling constant. Furthermore, we denoted by Ta the gauge group generators in the adjoint representation. In the nonplanar case, there are also contributions with multi-traces, suppressed in the large Nc 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 AN(L), 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(2.2)A4,D=4(1)~∫d4pp2(p-k1)2(p-k1-k2)2(p+k4)2,
where p is the momentum running along the loop and ki2=0 are the momenta of the external particles. We recognize two kinds of divergences. First, from the region pμ~0, 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 pμ~αkiμ 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 pμ~αkiμ with α~0, 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 inD=4-2ϵ 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) (2.3)A4,D=4(1)→A4,D=4-2ϵ(1).(In practice, we do the computation in four dimensions and once we arrive to the final result (2.2), we replace d4p→d4-2ϵp.)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/ϵ2. In the general case of L loops, we expect the same phenomenon to occur for each subloop, hence we expect a divergence of the form(2.4)AN(L)~1ϵ2L+⋯.
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(2.5)𝒜(++⋯+)=𝒜(-+⋯+)=0.
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(2.6)MN(L)(ϵ)=AN(L)(ϵ)AN(0),
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:(2.7)ℳ≡1+∑ℓ=1αℓMn(ℓ)=exp[∑ℓ=1(f(ℓ)(ϵ)Mn1(ℓϵ)+Cℓ+𝒪(ϵ))],
where(2.8)f(ℓ)(ϵ)=f0(ℓ)+ϵf1(ℓ)+ϵ2f2(ℓ),α≈λμ2ϵα 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 amplitudeMn1(ℓϵ). Note that this relation is highly nontrivial, since the constants fi(ℓ) and Cℓ 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(2.9)A4=Atree(Adiv,s)2(Adiv,t)2exp(f(λ)8(logst)2+const),(2.10)Adivs=exp(-18ϵ2f(-2)(λμ2ϵsϵ)-14ϵg(-1)(λμ2ϵsϵ)),
where s and t 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 f(λ)=(λ∂λ)2f(-2)(λ) controls both the leading divergent piece and the finite piece while the so-called collinear anomalous dimension, g(λ), 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 amplitudesAN(k1,k2,…,kN) have a simple universal behavior as the momenta of two adjacent gluons, let us say that k1 and k2 become collinear, (very) schematically: (2.11)AN(L)(k1,k2,k3,…,kN)→R(z)AN-1(L)(kt,k3,…,kN)(The correct expression contains also terms involving the N-1 point function at lower order.)with kt=k1+k2, k1→zkt in the collinear limit and R(z) 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 dimensionf(λ) 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 TrϕDSϕ, namely,(2.12)Δ-S=f(λ)logS+⋯
for large S, 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 [8–10], giving(2.13)f(λ)=λ2π2+⋯,λ≪1,f(λ)=λπ+⋯,λ≫1.
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 dimensiong(λ), note from the expression (2.9) that it is scheme dependent. For instance, if we change μ→αμ, then g(λ) changes by a term proportional to the cusp anomalous dimension f(λ). Hence the cusp anomalous dimension is a more "robust" quantity. This is also related to the fact that f(λ) is related to the anomalous dimension of an operator, as described above, while g(λ) 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 AdS5×S5. The duality provides us with a dictionary between the parameters on both sides of the correspondence, more precisely(3.1)1Nc≈gs,λ≡g2Nc=R2α′,
where gs is the string coupling and R is the radius of the AdS5 that coincides with that of the S5. 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 aD-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 aD-brane localized in the radial direction. In other words, we start with the AdS5 metric written in Poincare coordinates(3.2)ds2=R2dx3+12+dz2z2,
and we place a D-brane at some fixed large value of z=zIR and extending along the x3+1 coordinates. The asymptotic states are open strings that end on that D-brane. We then consider the scattering of these open strings, that will have the interpretation of the gluons.The proper momentum of the strings iskpr=kzIR/R, where k is the momentum conjugate to x3+1, and it plays the role of gauge theory momentum and will be kept fixed as we take away the IR cut-off, zIR→∞. 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 IRD-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 theD-brane, z=zIR 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 ofT-dual coordinates yμ, defined as follows:(3.3)ds2=w2(z)dxμdxμ→∂αyμ=iw2(z)ϵαβ∂βxμ.
Note that we do not T-dualize along the radial direction. After defining r=R2/z we end up again with formally the same AdS5 metric(3.4)ds2=R2dyμdyμ+dr2r2.In terms of the (3.4) coordinates, the boundary of the world-sheet is located at r=R2/zIR. Besides, the boundary conditions for the original coordinates xμ, which are that they carry momentum kμ, translate into the condition that yμ has “winding” Δyμ=2πkμ. Summarizing, the boundary of the world-sheet is located at r=R2/zIR and is a particular line constructed as follows (see Figure 3):(i)
for each particle of momentumkμ, draw a segment joining two points separated by Δyμ=2πkμ;(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 inT-dual coordinates, will then end in such sequence of light-like segments (see Figure 2) located at r=R2/zIR.Figure 2
Comparison of the world-sheet in original andT-dual coordinates. The hyperplane on the picture to the right should not be interpreted as a D-brane.Figure 3
Polygon of light-like segments corresponding to the momenta of six external particles undergoing a scattering.BeforeT-dualizing, we can interpret the regulator D-brane as several D3-branes, one per each external particle, at a fixed radial distance and extended in the four directions x1,3. 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 fourT-dualities, these D3-branes become D(-1) branes or D(-1) instantons, and the external states become open strings stretching a light-like distance Δyμ=2πkμ between two such D(-1) instantons. These instantons are located at the vertices, or cusps, of the above mentioned polygons.As we take away the IR cut-off,zIR→∞, the boundary of the world-sheet moves towards the boundary of the T-dual metric, at r=0. 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 theN-point scattering amplitude is given by the area A of the minimal surface that ends on a sequence of light-like segments on the boundary(3.5)𝒜N~e-(λ/2π)A(k1,…,kN).
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 1/λ. 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,k1+k3→k2+k4, and define the usual Mandelstam variables(3.6)s=-(k1+k2)2,t=-(k2+k3)2.According to our prescription we need to find the minimal surface ending in the following light-like polygon, see Figure4.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 AdS3 subspace of the full AdS5(3.7)ds2=-dy02+dy12+dr2r2.We are interested in computing the surface ending on a light-like Wilson loop alongy1=±y0, y0>0; 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:(3.8)y0=eτcoshσ,y1=eτsinhσ,r=eτw.Boosts and scaling transformations are simply shifts ofσ and τ. The Nambu-Goto action becomes(3.9)SNG=R22π∫dσdτ1-(w(τ)+w′(τ))2w(τ)2.One can explicitly check thatw(τ)=2 solves the equations of motion and has the correct boundary conditions. Hence the surface is given by(3.10)r=2y02-y12.


#### 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(r,y0,y1,y2), setting y3=0, and parameterize the surface by its projection to the (y1,y2) plane. In this case we obtain an action for two fields, r and t, living in the space parameterized by y1 and y2(3.11)S=R22π∫dy1dy21+(∂ir)2-(∂iy0)2-(∂1r∂2y0-∂2r∂1y0)2r2,
where r=r(y1,y2), ∂ir=∂yir, 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 s=t, where the projection of the Wilson lines is a square. By scale invariance we can choose the edges of the square to be at y1,y2=±1. The boundary conditions are then given by(3.12)r(±1,y2)=r(y1,±1)=0,y0(±1,y2)=±y2,y0(y1,±1)=±y1.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(3.13)y0(y1,y2)=y1y2,r(y1,y2)=(1-y12)(1-y22).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 formf(s/t). Hence, in order to capture the nontrivial dependence of the amplitude on the kinematics, we need to consider more general solutions with s≠t. In this case the projection of the surface to the (y1,y2) plane will not be a square but a rhombus, with s and t 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 AdS5 as the following surface embedded in R2,4:(3.14)-Y-12-Y02+Y12+Y22+Y32+Y42=-1.The relation between embedding and Poincare coordinates is(3.15)Yμ=yμr,μ=0,…,3,Y-1+Y4=1r,Y1+Y4=r2+yμyμr.The surface (3.13) is then given by(3.16)Y0Y-1=Y1Y2,Y3=Y4=0.Once we have written our solution in embedding coordinates, we notice that we can applySO(2,4) transformations, that are linearly realized in this coordinates, in order to obtain new solutions. This SO(2,4) symmetry is sometimes referred to as “dual conformal symmetry” and should not be confused with the original SO(2,4) symmetry associated to the original AdS space. Later we will comment on the implications of this symmetry.Solutions withs≠t can be obtained by starting from (3.13) and performing a boost in the 0-4 direction(3.17)Y0Y-1=Y1Y2,Y4=0→Y4-vY0=0,1-v2Y0Y-1=Y1,Y2.
Hence, we end up with a two-parameter solution, the parameter a, related to the size of the initial square and the parameter b, related to the boost parameter. The solution can be conveniently written as(3.18)r=acoshu1coshu2+bsinhu1sinhu2,y0=a1+b2sinhu1sinhu2coshu1coshu2+bsinhu1sinhu2,(3.19)y1=asinhu1coshu2coshu1coshu2+bsinhu1sinhu2,y2=acoshu1sinhu2coshu1coshu2+bsinhu1sinhu2,
where we have written the surface as a solution to the equations of motion in conformal gauge(3.20)iS=-R22π∫ℒ=-R22π∫du1du212(∂r∂r+∂yμ∂yμ)r2;a and b encode the kinematical information of the scattering as follows:(3.21)-s(2π)2=8a2(1-b)2,-t(2π)2=8a2(1+b)2,st=(1+b)2(1-b)2.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 directionr=r0, namely, in (3.12) r(±1,y2)=r(y1,±1)=r0, and then take r0→0; 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 inD=4-2ϵ 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 inD=4-2ϵ dimensions. More precisely, one starts with 𝒩=1 in ten dimensions and then dimensionally reduce to 4-2ϵ dimensions. For integer 2ϵ this is precisely the low energy theory living on a Dp-brane, where p=3-2ϵ. 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(3.22)ds2=f-1/2dx4-2ϵ2+f1/2[dr2+r2dΩ5+2ϵ2],f=(4π2eγ)ϵΓ(2+ϵ)μ2ϵλr4+2ϵ.
Following the steps described above, we are led to the following action:(3.23)S=cϵλμϵ2π∫ℒϵ=0rϵ,
where ℒϵ=0 is the Lagrangian density for AdS5. The presence of the factor rϵ will have two important effects. On one hand, previously divergent integrals will now converge (if ϵ<0). 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 ϵ→0. 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(3.24)S≈λμϵaϵF21(12,-ϵ2,1-ϵ2;b2).
Expanding in powers of ϵ we get the final answer(3.25)𝒜=e-(λ/2π)A,-λ2πA=iSdiv+λ8π(logst)2+C̃,Sdiv=2Sdiv,s+2Sdiv,t,iSdiv,s=-1ϵ212πλμ2ϵ(-s)ϵ-1ϵ14π(1-log2)λμ2ϵ(-s)ϵ.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], f(λ)=λ/π+⋯, 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 dualSO(2,4) 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 AdS5 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 s and t, since they can be sent to arbitrary values by a dual conformal transformation. The whole dependence on s and t 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(3.26)Anreg=Div+Fin(xi).
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 xi 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 [19–21](3.27)∑i(2xiμ(xi·∂xi)-xi2∂xiμ)Fin(xi)=f(λ)4∑i(xi-1μ+xi+1μ-2xiμ)log(xi-1,i+12),
where f(λ) 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 n=4 and n=5, 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 f(λ)≈λ 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(3.28)Astrong=ABDS+F(xij2xkl2xik2xjl2),
where F is the so-called remainder function. The BDS ansatz, if true, would imply that F is a constant. On the other hand, note that this is indeed the case for n=4,5, since in this case we cannot construct cross-ratios of the form xij2xkl2/xik2xjl2. On the other hand, we can construct such cross-ratios for n=6, 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, 22–24] and also at strong coupling [25, 26].


## 3.1. Set-Up of the Computation

As a first IR regulator we consider aD-brane localized in the radial direction. In other words, we start with the AdS5 metric written in Poincare coordinates(3.2)ds2=R2dx3+12+dz2z2,
and we place a D-brane at some fixed large value of z=zIR and extending along the x3+1 coordinates. The asymptotic states are open strings that end on that D-brane. We then consider the scattering of these open strings, that will have the interpretation of the gluons.The proper momentum of the strings iskpr=kzIR/R, where k is the momentum conjugate to x3+1, and it plays the role of gauge theory momentum and will be kept fixed as we take away the IR cut-off, zIR→∞. 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 IRD-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 theD-brane, z=zIR 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 ofT-dual coordinates yμ, defined as follows:(3.3)ds2=w2(z)dxμdxμ→∂αyμ=iw2(z)ϵαβ∂βxμ.
Note that we do not T-dualize along the radial direction. After defining r=R2/z we end up again with formally the same AdS5 metric(3.4)ds2=R2dyμdyμ+dr2r2.In terms of the (3.4) coordinates, the boundary of the world-sheet is located at r=R2/zIR. Besides, the boundary conditions for the original coordinates xμ, which are that they carry momentum kμ, translate into the condition that yμ has “winding” Δyμ=2πkμ. Summarizing, the boundary of the world-sheet is located at r=R2/zIR and is a particular line constructed as follows (see Figure 3):(i)
for each particle of momentumkμ, draw a segment joining two points separated by Δyμ=2πkμ;(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 inT-dual coordinates, will then end in such sequence of light-like segments (see Figure 2) located at r=R2/zIR.Figure 2
Comparison of the world-sheet in original andT-dual coordinates. The hyperplane on the picture to the right should not be interpreted as a D-brane.Figure 3
Polygon of light-like segments corresponding to the momenta of six external particles undergoing a scattering.BeforeT-dualizing, we can interpret the regulator D-brane as several D3-branes, one per each external particle, at a fixed radial distance and extended in the four directions x1,3. 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 fourT-dualities, these D3-branes become D(-1) branes or D(-1) instantons, and the external states become open strings stretching a light-like distance Δyμ=2πkμ between two such D(-1) instantons. These instantons are located at the vertices, or cusps, of the above mentioned polygons.As we take away the IR cut-off,zIR→∞, the boundary of the world-sheet moves towards the boundary of the T-dual metric, at r=0. 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 theN-point scattering amplitude is given by the area A of the minimal surface that ends on a sequence of light-like segments on the boundary(3.5)𝒜N~e-(λ/2π)A(k1,…,kN).
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 1/λ. 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,k1+k3→k2+k4, and define the usual Mandelstam variables(3.6)s=-(k1+k2)2,t=-(k2+k3)2.According to our prescription we need to find the minimal surface ending in the following light-like polygon, see Figure4.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 AdS3 subspace of the full AdS5(3.7)ds2=-dy02+dy12+dr2r2.We are interested in computing the surface ending on a light-like Wilson loop alongy1=±y0, y0>0; 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:(3.8)y0=eτcoshσ,y1=eτsinhσ,r=eτw.Boosts and scaling transformations are simply shifts ofσ and τ. The Nambu-Goto action becomes(3.9)SNG=R22π∫dσdτ1-(w(τ)+w′(τ))2w(τ)2.One can explicitly check thatw(τ)=2 solves the equations of motion and has the correct boundary conditions. Hence the surface is given by(3.10)r=2y02-y12.


### 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(r,y0,y1,y2), setting y3=0, and parameterize the surface by its projection to the (y1,y2) plane. In this case we obtain an action for two fields, r and t, living in the space parameterized by y1 and y2(3.11)S=R22π∫dy1dy21+(∂ir)2-(∂iy0)2-(∂1r∂2y0-∂2r∂1y0)2r2,
where r=r(y1,y2), ∂ir=∂yir, 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 s=t, where the projection of the Wilson lines is a square. By scale invariance we can choose the edges of the square to be at y1,y2=±1. The boundary conditions are then given by(3.12)r(±1,y2)=r(y1,±1)=0,y0(±1,y2)=±y2,y0(y1,±1)=±y1.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(3.13)y0(y1,y2)=y1y2,r(y1,y2)=(1-y12)(1-y22).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 formf(s/t). Hence, in order to capture the nontrivial dependence of the amplitude on the kinematics, we need to consider more general solutions with s≠t. In this case the projection of the surface to the (y1,y2) plane will not be a square but a rhombus, with s and t 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 AdS5 as the following surface embedded in R2,4:(3.14)-Y-12-Y02+Y12+Y22+Y32+Y42=-1.The relation between embedding and Poincare coordinates is(3.15)Yμ=yμr,μ=0,…,3,Y-1+Y4=1r,Y1+Y4=r2+yμyμr.The surface (3.13) is then given by(3.16)Y0Y-1=Y1Y2,Y3=Y4=0.Once we have written our solution in embedding coordinates, we notice that we can applySO(2,4) transformations, that are linearly realized in this coordinates, in order to obtain new solutions. This SO(2,4) symmetry is sometimes referred to as “dual conformal symmetry” and should not be confused with the original SO(2,4) symmetry associated to the original AdS space. Later we will comment on the implications of this symmetry.Solutions withs≠t can be obtained by starting from (3.13) and performing a boost in the 0-4 direction(3.17)Y0Y-1=Y1Y2,Y4=0→Y4-vY0=0,1-v2Y0Y-1=Y1,Y2.
Hence, we end up with a two-parameter solution, the parameter a, related to the size of the initial square and the parameter b, related to the boost parameter. The solution can be conveniently written as(3.18)r=acoshu1coshu2+bsinhu1sinhu2,y0=a1+b2sinhu1sinhu2coshu1coshu2+bsinhu1sinhu2,(3.19)y1=asinhu1coshu2coshu1coshu2+bsinhu1sinhu2,y2=acoshu1sinhu2coshu1coshu2+bsinhu1sinhu2,
where we have written the surface as a solution to the equations of motion in conformal gauge(3.20)iS=-R22π∫ℒ=-R22π∫du1du212(∂r∂r+∂yμ∂yμ)r2;a and b encode the kinematical information of the scattering as follows:(3.21)-s(2π)2=8a2(1-b)2,-t(2π)2=8a2(1+b)2,st=(1+b)2(1-b)2.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 directionr=r0, namely, in (3.12) r(±1,y2)=r(y1,±1)=r0, and then take r0→0; 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 inD=4-2ϵ 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 inD=4-2ϵ dimensions. More precisely, one starts with 𝒩=1 in ten dimensions and then dimensionally reduce to 4-2ϵ dimensions. For integer 2ϵ this is precisely the low energy theory living on a Dp-brane, where p=3-2ϵ. 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(3.22)ds2=f-1/2dx4-2ϵ2+f1/2[dr2+r2dΩ5+2ϵ2],f=(4π2eγ)ϵΓ(2+ϵ)μ2ϵλr4+2ϵ.
Following the steps described above, we are led to the following action:(3.23)S=cϵλμϵ2π∫ℒϵ=0rϵ,
where ℒϵ=0 is the Lagrangian density for AdS5. The presence of the factor rϵ will have two important effects. On one hand, previously divergent integrals will now converge (if ϵ<0). 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 ϵ→0. 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(3.24)S≈λμϵaϵF21(12,-ϵ2,1-ϵ2;b2).
Expanding in powers of ϵ we get the final answer(3.25)𝒜=e-(λ/2π)A,-λ2πA=iSdiv+λ8π(logst)2+C̃,Sdiv=2Sdiv,s+2Sdiv,t,iSdiv,s=-1ϵ212πλμ2ϵ(-s)ϵ-1ϵ14π(1-log2)λμ2ϵ(-s)ϵ.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], f(λ)=λ/π+⋯, we see that the leading divergence, as well as the finite piece, has the correct kinematical dependence with the correct overall coefficient.


## 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 AdS3 subspace of the full AdS5(3.7)ds2=-dy02+dy12+dr2r2.We are interested in computing the surface ending on a light-like Wilson loop alongy1=±y0, y0>0; 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:(3.8)y0=eτcoshσ,y1=eτsinhσ,r=eτw.Boosts and scaling transformations are simply shifts ofσ and τ. The Nambu-Goto action becomes(3.9)SNG=R22π∫dσdτ1-(w(τ)+w′(τ))2w(τ)2.One can explicitly check thatw(τ)=2 solves the equations of motion and has the correct boundary conditions. Hence the surface is given by(3.10)r=2y02-y12.


## 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(r,y0,y1,y2), setting y3=0, and parameterize the surface by its projection to the (y1,y2) plane. In this case we obtain an action for two fields, r and t, living in the space parameterized by y1 and y2(3.11)S=R22π∫dy1dy21+(∂ir)2-(∂iy0)2-(∂1r∂2y0-∂2r∂1y0)2r2,
where r=r(y1,y2), ∂ir=∂yir, 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 s=t, where the projection of the Wilson lines is a square. By scale invariance we can choose the edges of the square to be at y1,y2=±1. The boundary conditions are then given by(3.12)r(±1,y2)=r(y1,±1)=0,y0(±1,y2)=±y2,y0(y1,±1)=±y1.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(3.13)y0(y1,y2)=y1y2,r(y1,y2)=(1-y12)(1-y22).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 formf(s/t). Hence, in order to capture the nontrivial dependence of the amplitude on the kinematics, we need to consider more general solutions with s≠t. In this case the projection of the surface to the (y1,y2) plane will not be a square but a rhombus, with s and t 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 AdS5 as the following surface embedded in R2,4:(3.14)-Y-12-Y02+Y12+Y22+Y32+Y42=-1.The relation between embedding and Poincare coordinates is(3.15)Yμ=yμr,μ=0,…,3,Y-1+Y4=1r,Y1+Y4=r2+yμyμr.The surface (3.13) is then given by(3.16)Y0Y-1=Y1Y2,Y3=Y4=0.Once we have written our solution in embedding coordinates, we notice that we can applySO(2,4) transformations, that are linearly realized in this coordinates, in order to obtain new solutions. This SO(2,4) symmetry is sometimes referred to as “dual conformal symmetry” and should not be confused with the original SO(2,4) symmetry associated to the original AdS space. Later we will comment on the implications of this symmetry.Solutions withs≠t can be obtained by starting from (3.13) and performing a boost in the 0-4 direction(3.17)Y0Y-1=Y1Y2,Y4=0→Y4-vY0=0,1-v2Y0Y-1=Y1,Y2.
Hence, we end up with a two-parameter solution, the parameter a, related to the size of the initial square and the parameter b, related to the boost parameter. The solution can be conveniently written as(3.18)r=acoshu1coshu2+bsinhu1sinhu2,y0=a1+b2sinhu1sinhu2coshu1coshu2+bsinhu1sinhu2,(3.19)y1=asinhu1coshu2coshu1coshu2+bsinhu1sinhu2,y2=acoshu1sinhu2coshu1coshu2+bsinhu1sinhu2,
where we have written the surface as a solution to the equations of motion in conformal gauge(3.20)iS=-R22π∫ℒ=-R22π∫du1du212(∂r∂r+∂yμ∂yμ)r2;a and b encode the kinematical information of the scattering as follows:(3.21)-s(2π)2=8a2(1-b)2,-t(2π)2=8a2(1+b)2,st=(1+b)2(1-b)2.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 directionr=r0, namely, in (3.12) r(±1,y2)=r(y1,±1)=r0, and then take r0→0; 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 inD=4-2ϵ 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 inD=4-2ϵ dimensions. More precisely, one starts with 𝒩=1 in ten dimensions and then dimensionally reduce to 4-2ϵ dimensions. For integer 2ϵ this is precisely the low energy theory living on a Dp-brane, where p=3-2ϵ. 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(3.22)ds2=f-1/2dx4-2ϵ2+f1/2[dr2+r2dΩ5+2ϵ2],f=(4π2eγ)ϵΓ(2+ϵ)μ2ϵλr4+2ϵ.
Following the steps described above, we are led to the following action:(3.23)S=cϵλμϵ2π∫ℒϵ=0rϵ,
where ℒϵ=0 is the Lagrangian density for AdS5. The presence of the factor rϵ will have two important effects. On one hand, previously divergent integrals will now converge (if ϵ<0). 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 ϵ→0. 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(3.24)S≈λμϵaϵF21(12,-ϵ2,1-ϵ2;b2).
Expanding in powers of ϵ we get the final answer(3.25)𝒜=e-(λ/2π)A,-λ2πA=iSdiv+λ8π(logst)2+C̃,Sdiv=2Sdiv,s+2Sdiv,t,iSdiv,s=-1ϵ212πλμ2ϵ(-s)ϵ-1ϵ14π(1-log2)λμ2ϵ(-s)ϵ.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], f(λ)=λ/π+⋯, 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 dualSO(2,4) 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 AdS5 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 s and t, since they can be sent to arbitrary values by a dual conformal transformation. The whole dependence on s and t 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(3.26)Anreg=Div+Fin(xi).
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 xi 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 [19–21](3.27)∑i(2xiμ(xi·∂xi)-xi2∂xiμ)Fin(xi)=f(λ)4∑i(xi-1μ+xi+1μ-2xiμ)log(xi-1,i+12),
where f(λ) 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 n=4 and n=5, 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 f(λ)≈λ 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(3.28)Astrong=ABDS+F(xij2xkl2xik2xjl2),
where F is the so-called remainder function. The BDS ansatz, if true, would imply that F is a constant. On the other hand, note that this is indeed the case for n=4,5, since in this case we cannot construct cross-ratios of the form xij2xkl2/xik2xjl2. On the other hand, we can construct such cross-ratios for n=6, 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, 22–24] 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 [26–29], 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 cuspsxi, which satisfy xi,i+12=0. See Figure 6Figure 6
The light-like polygon in which the minimal surface ends is parameterized by the position of the cusps,xi, with xi,i+12=0.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χijkl=xij2xkl2/xik2xjl2. 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 AdS5; in which case there are 3N-15 cross-ratios if we are dealing with the scattering of N 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 AdS3; in this case, we have N-6 cross-ratios (For the general scattering in four dimensions we have 4N coordinates, minus N, since the distance between consecutive points has to be light-like, minus 15, that is the dimension of the conformal group SO(2,4). In the case of AdS3, we have 2N-N minus 6, which is the dimension of SO(2,2).), and the polygon is a zig-zaged polygon living in one-plus-one dimensions, see Figure 7Figure 7
A zig-zaged Wilson loop in1+1 dimensions is parameterized by nxi+ coordinates and nxi- 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, henceN=2n. As one can see in Figure 7, the contour is parameterized by nxi+ coordinates and nxi- coordinates. With each set we can form n-3 invariant cross-ratios, of the form χijkl±=xij±xkl±/xik±xjl±, and all together we have the expected 2n-6 cross-ratios.In order to consider minimal surfaces inAdS3 we need to consider the world-sheet of classical strings on AdS3. This is the subject of the following subsection.


### 4.2. Strings on AdS3

Classical strings onAdS3 can be described in terms of embedding coordinates, where we see AdS3 as the following surface embedded in R2,2:(4.1)Y·Y≡-Y-12-Y02+Y12+Y22=-1.
We take the world-sheet to be the whole complex plane and Y=Y(z,z¯). Since we are interested in classical solutions, the fields have to satisfy the conformal gauge equations of motion(4.2)∂∂¯Y-(∂Y·∂¯Y)Y=0.
Furthermore, the equations of motion are supplemented by the Virasoro constraints(4.3)∂Y·∂Y=∂¯Y·∂¯Y=0.
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, [30–32] and consider the "reduced" fields(4.4)α=log∂Y·∂¯Y,p2=∂2Y·∂2Y.
As a consequence of the equations of motion and Virasoro constraints, p can be seen to be a holomorphic function, and p=p(z) and α(z,z¯) can be seen to satisfy a generalized version of the Sinh-Gordon equation(4.5)∂∂¯α-eα+p(z)p¯(z¯)e-α=0.
Hence, we have reduced the initial set of equations to a holomorphic function p(z) and to a field α(z,z¯) 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(4.6)𝒜=∫eαd2z.
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 onAdS3 is mapped to a holomorphic function p(z) and a field α satisfying (4.5). We can also wonder whether given a holomorphic function p(z) and a field α satisfying (4.5) it is possible to construct a solution of classical strings on AdS3. 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(4.7)(d+BL)ψaL=0,(d+BR)ψȧR=0,
where the flat connections BL,R are simple two-by-two matrices constructed from p(z) and α(z,z¯), for instance(4.8)BzL=(∂α412eα/212pe-α/2-∂α4);
we will denote the different components of the connections by BαβL and Bα̇β̇R. On the other hand, the indices a and ȧ in (4.7) denote independent solutions of the auxiliary linear problems. Each ψaL or ψȧR is then a doublet. We denote the components of this doublet by ψα,aL 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(4.9)Ya,ȧ=(Y-1+Y2Y1-Y0Y1+Y0Y-1+Y2)α,α̇=ψα,aLMαβ̇ψβ̇,ȧR
for some constant matrix M (which can be taken, for instance, to be the identity matrix). One can show that Y constructed this way satisfies all the required properties. In order to write down the answer, we have written the space-time coordinates in an SL(2)×SL(2)≈SO(2,2) 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(4.10)Bz→Bz(ζ)=14(∂α00-∂α)+1ζ12(0eα/2pe-α/20),Bz¯→Bz¯(ζ)=14(-∂¯α00∂¯α)+ζ12(0p¯e-α/2eα/20).
One can check that the connections are still flat, namely, they satisfy(4.11)dB+B∧B=∂Bz¯-∂¯Bz+[Bz,Bz¯]=0,
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(4.12)B(ζ=1)=BL,B(ζ=i)=BR.Before proceeding, let us make the following remark. Sincep(z) is a holomorphic function, it is possible to make a change of coordinates from the z-plane to the w-plane, where dw=p(z)dz. In the w-plane, the generalized Sinh-Gordon equation takes the usual form (after a simple field redefinition of α)(4.13)α=α̂+14logpp¯→∂w∂¯wα̂=sinh2α̂.
Naively, it would seem that we got rid of all the information on p(z). However, this is not the case, since the w-plane will have in general a complicated structure (for instance, it will have branch cuts, etc., depending on p(z)). 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 w-plane, the area is simply given by(4.14)A=∫eα̂d2w.


### 4.3. Classical Solutions Corresponding to Scattering Amplitudes

A natural question to ask is what are the properties of the holomorphic functionp(z) and α(z,z¯) 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(4.15)p(z)=1,α=α̂=0.
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 |z|. So we expect that α̂→0 as |z| becomes large. Finally, if we are interested in a minimal surface ending on a polygon with 2n cusps, we propose p(z) to be a polynomial of degree n-2:(4.16)p(z)=zn-2+cn-4zn-4+⋯+c0;
we have used rescalings and translations in order to set the coefficients of zn-2 and zn-3 to one and zero respectively. We see that such polynomial contains n-3 complex coefficients, or 2n-6 real coefficients, which exactly agrees with the amount of expected independent cross-ratios. In conclusion, minimal surfaces ending on a light-like polygon with 2n cusps correspond to a holomorphic polynomial of degree n-2 and a field α̂ satisfying the Sinh-Gordon equations, with boundary conditions such that it decays at infinity and diverges logarithmically at the zeroes of p(z).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(4.17)Areg=∫(eα̂-1)d2w.
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 Areg.)


#### 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|z|, or |w|, since this will tell us the behavior of the world-sheet near its boundary. Suppose, by simplicity, that p(z)=zn-2, hence w≈zn/2. As a result, as we go once around the z-plane, we go around the w-plane n/2 times.Since for large values of|w|α̂ vanishes, the flat connections BL,R drastically are simplified and we can solve for the auxiliary linear problems. A general solution will be of the form(4.18)ψaL≈ca+(10)ew+w¯+ca-(01)e-(w+w¯),ψȧR≈dȧ+(10)e(w-w¯)/i+dȧ-(01)e-(w-w¯)/i.
We see that the w-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 ca+ dominates in the left problem while the solution proportional to dȧ+ dominates in the right problem. For instance, for the first quadrant we have Yaȧ≈(Large)×ca+dȧ+ at large values of |w|. This means that for large values of |w|, 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) thew-plane is naturally divided into two parts, according to the sign of Re(w). In the same way, when looking at the right problem, the w-plane is naturally divided into two parts, according to the sign of Im(w). Hence, the w-plane is naturally divided into four quadrants. Large values of w 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(4.19)Ya,ȧ≈λaλ̃ȧ,
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 w-plane n/2 times, we get the expected 2n cusps.The above solution, in whichp(z) is a homogeneous polynomial, turns out to correspond to a regular polygon with an additional Z2n 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.Y-System for Minimal Surfaces

Let us focus on the left problem. We see that each sheet on thew-plane is naturally divided into two sectors, one with Re(w)>0 and the other with Re(w)<0. 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)
Thew-plane is divided into n sectors, since each sheet contains two sectors. We label this sectors by i=0,…,n-1.(ii)
We callsiL the small solution at the ith sector. By definition, this is the solution with the fastest decay along the line in the center of the ith sector, for large |w|. For instance, given solution (4.18), we see that the term proportional to e-(w+w¯) has the fastest decay as we move along the positive real axis in the w-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 areSL(2) matrices, we can introduce an SL(2) invariant product(4.20)ψaL∧ψbL≡ϵαβψα,aLψβ,bL.
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 z. Hence, we can normalize our solutions such that(4.21)ψaL∧ψbL=ϵab.
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 SL(2) invariant product just introduced, more precisely(4.22)ψaL∧siL≈λai.
How do we construct space-time cross-ratios? We have seen that the location of the cusps is given by Yaȧi=λaiλ̃ȧi. The space time cross-ratios involve distances like(4.23)Yi·Yj=ϵabϵȧḃYaȧiYbḃj=〈λiλj〉〈λ̃iλ̃j〉,〈λiλj〉=ϵabλaiλbj.
When seeing ψaL 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(4.24)〈λiλj〉≈〈(ψaL∧siL)(ψaL∧sjL)〉=siL∧sjL.
Which means that space-time cross-ratios can be constructed from inner products of the small solutions in the corresponding sectors, more precisely(4.25)xij+xkl+xik+xjl+=(siL∧sjL)(skL∧slL)(siL∧skL)(sjL∧slL).
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(4.26)(d+B(ζ))si(ζ)=0;
then, we can consider the cross-ratios as a function of such spectral parameter(4.27)χijkl(ζ)=(si∧sj)(sk∧sl)(si∧sk)(sj∧sl).
Then, if we set ζ=1 we obtain the physical plus/left cross-ratios while setting ζ=i we obtain the physical minus/right cross-ratios:(4.28)χijkl(ζ=1)=χijkl+,χijkl(ζ=i)=χijkl-.
A property of the flat connection B(ζ) that will be used in what follows is that it possesses a Z2 symmetry(4.29)B(eiπζ)=σ3B(ζ)σ3,
where σ3 is the usual Pauli matrix. This symmetry allows to relate small sections at different values of the spectral parameter, for instance, si+1(ζ)=σ3si(eiπζ), and in particular, it implies(4.30)si∧sj(eiπζ)=si+1∧sj+1(ζ).
Besides, in order to simplify the subsequent derivation, we will assume that si∧si+1=1.Now we have all the elements to derive the so-called Hirota equations and subsequently theY-system equations. The trick is to choose s0 and s1 as a complete basis of small flat sections and express two arbitrary consecutive flat sections sk and sk+1 in terms of them:(4.31)sk=(sk∧s1)s0-(sk∧s0)s1,sk+1=(sk+1∧s1)s0-(sk+1∧s0)s1.
Next, we use (4.30) in order to express every wedge as a wedge involving s0 and consider 1=sk∧sk+1; then we obtain(4.32)-(sk-1∧s0)++(sk+1∧s0)+(sk∧s0)++(sk∧s0)=1,
where we have introduced the notations f±=f(e±iπ/2ζ), f++=f(eiπζ), and so forth. This equation can be written in a more suggestive way by introducing Tk=s0∧sk+1(e-i(k+1)π/2ζ). We obtain(4.33)Ts+Ts-=Ts+1Ts-1+1
which has the form of the so-called Hirota equations. From the definition of Ts, we see that it is nontrivial for s=0,…,n-2. The Y-system equations can be simply obtained by introducing Ys≡Ts-1Ts+1; we then obtain(4.34)Ys+Ys-=(1+Ys+1)(1+Ys-1).Ys is nontrivial for s=1,…,n-3. Notice that this agrees with the amount of (complex) cross-ratios of our scattering problem. These are the Y-system equations. Note that they have simply followed from a chain of rather trivial facts. These are functional equations for Ys(ζ) and are valid for any value of ζ. One could reintroduce the normalized factors si∧si+1 and check that the Y-functions are given by the usual cross-ratios introduced above. The physical cross-ratios are then obtained by evaluating Ys(ζ) at ζ=1 and ζ=i.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 functionp(z) 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 Ys(ζ) 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 theY-system equations (4.34) we need to specify the analytic properties of Ys(ζ). By analyzing the auxiliary linear problem and the definition of Ys(ζ) one can show that(i)
Ys(ζ) are analytically away from ζ=0,∞,(ii)
asζ→0,∞, 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 1/ζ.More precisely, by callingζ=eθ, one can show that for large θ(4.35)logYs≈-mscoshθ+⋯,
where ma is given by the periods of p(z)1/2 along the cycles γa, more precisely ma≈-∮γap(z)dz. This is how the information of the polynomial p(z) enters into the problem. These periods are usually complex, and there are n-3 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 theY-system equations with the analytic properties and boundary conditions for the Y-functions, in order to write a system of integral equations for them, whose solutions will automatically satisfy the Y-system equations and have the required boundary conditions. The system of integral equations is given by(4.36)logYs=-mscoshθ+K⋆log(1+Ys+1)(1+Ys-1),
where the convolution operation is defined by(4.37)K⋆f(θ)=∫-∞∞dθ′K(θ-θ′)f(θ′),K(θ)=12πcoshθ.
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 ma 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 Ys(θ) 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 p(z) are along the real line, such that p(z) is positive for z 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 theY-functions and has the simple expression(4.38)Areg=∑s∫dθms2πcoshθlog(1+Ys(θ)).
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 massesma. 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 inAdS3 can be given by using the integrability of the system. In the following, we highlight the results for the general case of AdS5 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 AdS3 case. We get a slightly more complicated Y-system for functions Ya,s, with a=1,2,3 and s=1,…,N-5. The system reads(4.39)Y2,m-Y2,m+Y1,mY3,m=(1+Y2,m+1)(1+Y2,m-1)(1+Y1,m)(1+Y3,m),Y3,m-Y1,m+Y2,m=(1+Y3,m+1)(1+Y1,m-1)1+Y2,m,Y1,m-Y3,m+Y2,m=(1+Y1,m+1)(1+Y3,m-1)1+Y2,m.
As before, in order to get a system of integral equations, we need to supplement these equations by the analytic properties of Ya,s(ζ). As for the AdS3 case, the large θ=logζ behavior can be understood by a WKB analysis. The boundary conditions for the Y-functions are(4.40)logY1,s→-mscoshθ-Cs±Ds+⋯,θ→±∞,logY3,s→-mscoshθ-Cs∓Ds+⋯,θ→±∞,logY2,s→-2mscoshθ+⋯,θ→±∞;
depending on the signature of space-time, the reality conditions of Cs and Ds are fixed. For instance, for (3,1) or (1,3) signature, the Ds are real while the Cs are purely imaginary. As for the AdS3 case, these boundary conditions together with regularity imply the following set of integral equations:(4.41)logY2,s=-ms2coshθ-K2⋆αs-K1⋆βs,logY1,s=-mscoshθ-Cs-12K2⋆β2-K1⋆αs-12K2⋆γs,logY3,s=-mscoshθ+Cs-12K2⋆β2-K1⋆αs+12K2⋆γs,
where(4.42)αs=log(1+Y1,s)(1+Y3,s)(1+Y2,s-1)(1+Y2,s+1),γs=log(1+Y1,s-1)(1+Y3,s+1)(1+Y1,s+1)(1+Y3,s-1),βs=log(1+Y2,s2)(1+Y1,s-1)(1+Y1,s+1)(1+Y3,s-1)(1+Y3,s+1),
and the kernels are given by(4.43)K1=12πcoshθ,K2=2coshθπcosh2θ,K3=iπtanh2θ.
The data of the system of integral equations are the (in general complex) masses ms and the purely imaginary constants Cs. Since for the scattering of N gluons s=1,…,N-5, we get 2(N-5)+N-5=3(N-5) degrees of freedom, which exactly agrees with the expected amount of cross ratios.Again, as for the case of theAdS3 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(4.44)Areg=∑s∫dθms2πcoshθlog((1+Y1,s)(1+Y3,s)(1+Y2,s)2).
This is the final solution for the general problem in AdS5 (where for simplicity we have assumed real masses).


## 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 cuspsxi, which satisfy xi,i+12=0. See Figure 6Figure 6
The light-like polygon in which the minimal surface ends is parameterized by the position of the cusps,xi, with xi,i+12=0.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χijkl=xij2xkl2/xik2xjl2. 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 AdS5; in which case there are 3N-15 cross-ratios if we are dealing with the scattering of N 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 AdS3; in this case, we have N-6 cross-ratios (For the general scattering in four dimensions we have 4N coordinates, minus N, since the distance between consecutive points has to be light-like, minus 15, that is the dimension of the conformal group SO(2,4). In the case of AdS3, we have 2N-N minus 6, which is the dimension of SO(2,2).), and the polygon is a zig-zaged polygon living in one-plus-one dimensions, see Figure 7Figure 7
A zig-zaged Wilson loop in1+1 dimensions is parameterized by nxi+ coordinates and nxi- 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, henceN=2n. As one can see in Figure 7, the contour is parameterized by nxi+ coordinates and nxi- coordinates. With each set we can form n-3 invariant cross-ratios, of the form χijkl±=xij±xkl±/xik±xjl±, and all together we have the expected 2n-6 cross-ratios.In order to consider minimal surfaces inAdS3 we need to consider the world-sheet of classical strings on AdS3. This is the subject of the following subsection.


## 4.2. Strings on AdS3

Classical strings onAdS3 can be described in terms of embedding coordinates, where we see AdS3 as the following surface embedded in R2,2:(4.1)Y·Y≡-Y-12-Y02+Y12+Y22=-1.
We take the world-sheet to be the whole complex plane and Y=Y(z,z¯). Since we are interested in classical solutions, the fields have to satisfy the conformal gauge equations of motion(4.2)∂∂¯Y-(∂Y·∂¯Y)Y=0.
Furthermore, the equations of motion are supplemented by the Virasoro constraints(4.3)∂Y·∂Y=∂¯Y·∂¯Y=0.
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, [30–32] and consider the "reduced" fields(4.4)α=log∂Y·∂¯Y,p2=∂2Y·∂2Y.
As a consequence of the equations of motion and Virasoro constraints, p can be seen to be a holomorphic function, and p=p(z) and α(z,z¯) can be seen to satisfy a generalized version of the Sinh-Gordon equation(4.5)∂∂¯α-eα+p(z)p¯(z¯)e-α=0.
Hence, we have reduced the initial set of equations to a holomorphic function p(z) and to a field α(z,z¯) 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(4.6)𝒜=∫eαd2z.
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 onAdS3 is mapped to a holomorphic function p(z) and a field α satisfying (4.5). We can also wonder whether given a holomorphic function p(z) and a field α satisfying (4.5) it is possible to construct a solution of classical strings on AdS3. 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(4.7)(d+BL)ψaL=0,(d+BR)ψȧR=0,
where the flat connections BL,R are simple two-by-two matrices constructed from p(z) and α(z,z¯), for instance(4.8)BzL=(∂α412eα/212pe-α/2-∂α4);
we will denote the different components of the connections by BαβL and Bα̇β̇R. On the other hand, the indices a and ȧ in (4.7) denote independent solutions of the auxiliary linear problems. Each ψaL or ψȧR is then a doublet. We denote the components of this doublet by ψα,aL 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(4.9)Ya,ȧ=(Y-1+Y2Y1-Y0Y1+Y0Y-1+Y2)α,α̇=ψα,aLMαβ̇ψβ̇,ȧR
for some constant matrix M (which can be taken, for instance, to be the identity matrix). One can show that Y constructed this way satisfies all the required properties. In order to write down the answer, we have written the space-time coordinates in an SL(2)×SL(2)≈SO(2,2) 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(4.10)Bz→Bz(ζ)=14(∂α00-∂α)+1ζ12(0eα/2pe-α/20),Bz¯→Bz¯(ζ)=14(-∂¯α00∂¯α)+ζ12(0p¯e-α/2eα/20).
One can check that the connections are still flat, namely, they satisfy(4.11)dB+B∧B=∂Bz¯-∂¯Bz+[Bz,Bz¯]=0,
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(4.12)B(ζ=1)=BL,B(ζ=i)=BR.Before proceeding, let us make the following remark. Sincep(z) is a holomorphic function, it is possible to make a change of coordinates from the z-plane to the w-plane, where dw=p(z)dz. In the w-plane, the generalized Sinh-Gordon equation takes the usual form (after a simple field redefinition of α)(4.13)α=α̂+14logpp¯→∂w∂¯wα̂=sinh2α̂.
Naively, it would seem that we got rid of all the information on p(z). However, this is not the case, since the w-plane will have in general a complicated structure (for instance, it will have branch cuts, etc., depending on p(z)). 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 w-plane, the area is simply given by(4.14)A=∫eα̂d2w.


## 4.3. Classical Solutions Corresponding to Scattering Amplitudes

A natural question to ask is what are the properties of the holomorphic functionp(z) and α(z,z¯) 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(4.15)p(z)=1,α=α̂=0.
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 |z|. So we expect that α̂→0 as |z| becomes large. Finally, if we are interested in a minimal surface ending on a polygon with 2n cusps, we propose p(z) to be a polynomial of degree n-2:(4.16)p(z)=zn-2+cn-4zn-4+⋯+c0;
we have used rescalings and translations in order to set the coefficients of zn-2 and zn-3 to one and zero respectively. We see that such polynomial contains n-3 complex coefficients, or 2n-6 real coefficients, which exactly agrees with the amount of expected independent cross-ratios. In conclusion, minimal surfaces ending on a light-like polygon with 2n cusps correspond to a holomorphic polynomial of degree n-2 and a field α̂ satisfying the Sinh-Gordon equations, with boundary conditions such that it decays at infinity and diverges logarithmically at the zeroes of p(z).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(4.17)Areg=∫(eα̂-1)d2w.
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 Areg.)


### 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|z|, or |w|, since this will tell us the behavior of the world-sheet near its boundary. Suppose, by simplicity, that p(z)=zn-2, hence w≈zn/2. As a result, as we go once around the z-plane, we go around the w-plane n/2 times.Since for large values of|w|α̂ vanishes, the flat connections BL,R drastically are simplified and we can solve for the auxiliary linear problems. A general solution will be of the form(4.18)ψaL≈ca+(10)ew+w¯+ca-(01)e-(w+w¯),ψȧR≈dȧ+(10)e(w-w¯)/i+dȧ-(01)e-(w-w¯)/i.
We see that the w-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 ca+ dominates in the left problem while the solution proportional to dȧ+ dominates in the right problem. For instance, for the first quadrant we have Yaȧ≈(Large)×ca+dȧ+ at large values of |w|. This means that for large values of |w|, 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) thew-plane is naturally divided into two parts, according to the sign of Re(w). In the same way, when looking at the right problem, the w-plane is naturally divided into two parts, according to the sign of Im(w). Hence, the w-plane is naturally divided into four quadrants. Large values of w 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(4.19)Ya,ȧ≈λaλ̃ȧ,
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 w-plane n/2 times, we get the expected 2n cusps.The above solution, in whichp(z) is a homogeneous polynomial, turns out to correspond to a regular polygon with an additional Z2n 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.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|z|, or |w|, since this will tell us the behavior of the world-sheet near its boundary. Suppose, by simplicity, that p(z)=zn-2, hence w≈zn/2. As a result, as we go once around the z-plane, we go around the w-plane n/2 times.Since for large values of|w|α̂ vanishes, the flat connections BL,R drastically are simplified and we can solve for the auxiliary linear problems. A general solution will be of the form(4.18)ψaL≈ca+(10)ew+w¯+ca-(01)e-(w+w¯),ψȧR≈dȧ+(10)e(w-w¯)/i+dȧ-(01)e-(w-w¯)/i.
We see that the w-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 ca+ dominates in the left problem while the solution proportional to dȧ+ dominates in the right problem. For instance, for the first quadrant we have Yaȧ≈(Large)×ca+dȧ+ at large values of |w|. This means that for large values of |w|, 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) thew-plane is naturally divided into two parts, according to the sign of Re(w). In the same way, when looking at the right problem, the w-plane is naturally divided into two parts, according to the sign of Im(w). Hence, the w-plane is naturally divided into four quadrants. Large values of w 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(4.19)Ya,ȧ≈λaλ̃ȧ,
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 w-plane n/2 times, we get the expected 2n cusps.The above solution, in whichp(z) is a homogeneous polynomial, turns out to correspond to a regular polygon with an additional Z2n 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.Y-System for Minimal Surfaces

Let us focus on the left problem. We see that each sheet on thew-plane is naturally divided into two sectors, one with Re(w)>0 and the other with Re(w)<0. 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)
Thew-plane is divided into n sectors, since each sheet contains two sectors. We label this sectors by i=0,…,n-1.(ii)
We callsiL the small solution at the ith sector. By definition, this is the solution with the fastest decay along the line in the center of the ith sector, for large |w|. For instance, given solution (4.18), we see that the term proportional to e-(w+w¯) has the fastest decay as we move along the positive real axis in the w-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 areSL(2) matrices, we can introduce an SL(2) invariant product(4.20)ψaL∧ψbL≡ϵαβψα,aLψβ,bL.
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 z. Hence, we can normalize our solutions such that(4.21)ψaL∧ψbL=ϵab.
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 SL(2) invariant product just introduced, more precisely(4.22)ψaL∧siL≈λai.
How do we construct space-time cross-ratios? We have seen that the location of the cusps is given by Yaȧi=λaiλ̃ȧi. The space time cross-ratios involve distances like(4.23)Yi·Yj=ϵabϵȧḃYaȧiYbḃj=〈λiλj〉〈λ̃iλ̃j〉,〈λiλj〉=ϵabλaiλbj.
When seeing ψaL 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(4.24)〈λiλj〉≈〈(ψaL∧siL)(ψaL∧sjL)〉=siL∧sjL.
Which means that space-time cross-ratios can be constructed from inner products of the small solutions in the corresponding sectors, more precisely(4.25)xij+xkl+xik+xjl+=(siL∧sjL)(skL∧slL)(siL∧skL)(sjL∧slL).
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(4.26)(d+B(ζ))si(ζ)=0;
then, we can consider the cross-ratios as a function of such spectral parameter(4.27)χijkl(ζ)=(si∧sj)(sk∧sl)(si∧sk)(sj∧sl).
Then, if we set ζ=1 we obtain the physical plus/left cross-ratios while setting ζ=i we obtain the physical minus/right cross-ratios:(4.28)χijkl(ζ=1)=χijkl+,χijkl(ζ=i)=χijkl-.
A property of the flat connection B(ζ) that will be used in what follows is that it possesses a Z2 symmetry(4.29)B(eiπζ)=σ3B(ζ)σ3,
where σ3 is the usual Pauli matrix. This symmetry allows to relate small sections at different values of the spectral parameter, for instance, si+1(ζ)=σ3si(eiπζ), and in particular, it implies(4.30)si∧sj(eiπζ)=si+1∧sj+1(ζ).
Besides, in order to simplify the subsequent derivation, we will assume that si∧si+1=1.Now we have all the elements to derive the so-called Hirota equations and subsequently theY-system equations. The trick is to choose s0 and s1 as a complete basis of small flat sections and express two arbitrary consecutive flat sections sk and sk+1 in terms of them:(4.31)sk=(sk∧s1)s0-(sk∧s0)s1,sk+1=(sk+1∧s1)s0-(sk+1∧s0)s1.
Next, we use (4.30) in order to express every wedge as a wedge involving s0 and consider 1=sk∧sk+1; then we obtain(4.32)-(sk-1∧s0)++(sk+1∧s0)+(sk∧s0)++(sk∧s0)=1,
where we have introduced the notations f±=f(e±iπ/2ζ), f++=f(eiπζ), and so forth. This equation can be written in a more suggestive way by introducing Tk=s0∧sk+1(e-i(k+1)π/2ζ). We obtain(4.33)Ts+Ts-=Ts+1Ts-1+1
which has the form of the so-called Hirota equations. From the definition of Ts, we see that it is nontrivial for s=0,…,n-2. The Y-system equations can be simply obtained by introducing Ys≡Ts-1Ts+1; we then obtain(4.34)Ys+Ys-=(1+Ys+1)(1+Ys-1).Ys is nontrivial for s=1,…,n-3. Notice that this agrees with the amount of (complex) cross-ratios of our scattering problem. These are the Y-system equations. Note that they have simply followed from a chain of rather trivial facts. These are functional equations for Ys(ζ) and are valid for any value of ζ. One could reintroduce the normalized factors si∧si+1 and check that the Y-functions are given by the usual cross-ratios introduced above. The physical cross-ratios are then obtained by evaluating Ys(ζ) at ζ=1 and ζ=i.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 functionp(z) 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 Ys(ζ) 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 theY-system equations (4.34) we need to specify the analytic properties of Ys(ζ). By analyzing the auxiliary linear problem and the definition of Ys(ζ) one can show that(i)
Ys(ζ) are analytically away from ζ=0,∞,(ii)
asζ→0,∞, 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 1/ζ.More precisely, by callingζ=eθ, one can show that for large θ(4.35)logYs≈-mscoshθ+⋯,
where ma is given by the periods of p(z)1/2 along the cycles γa, more precisely ma≈-∮γap(z)dz. This is how the information of the polynomial p(z) enters into the problem. These periods are usually complex, and there are n-3 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 theY-system equations with the analytic properties and boundary conditions for the Y-functions, in order to write a system of integral equations for them, whose solutions will automatically satisfy the Y-system equations and have the required boundary conditions. The system of integral equations is given by(4.36)logYs=-mscoshθ+K⋆log(1+Ys+1)(1+Ys-1),
where the convolution operation is defined by(4.37)K⋆f(θ)=∫-∞∞dθ′K(θ-θ′)f(θ′),K(θ)=12πcoshθ.
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 ma 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 Ys(θ) 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 p(z) are along the real line, such that p(z) is positive for z 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 theY-functions and has the simple expression(4.38)Areg=∑s∫dθms2πcoshθlog(1+Ys(θ)).
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 massesma. 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 inAdS3 can be given by using the integrability of the system. In the following, we highlight the results for the general case of AdS5 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 AdS3 case. We get a slightly more complicated Y-system for functions Ya,s, with a=1,2,3 and s=1,…,N-5. The system reads(4.39)Y2,m-Y2,m+Y1,mY3,m=(1+Y2,m+1)(1+Y2,m-1)(1+Y1,m)(1+Y3,m),Y3,m-Y1,m+Y2,m=(1+Y3,m+1)(1+Y1,m-1)1+Y2,m,Y1,m-Y3,m+Y2,m=(1+Y1,m+1)(1+Y3,m-1)1+Y2,m.
As before, in order to get a system of integral equations, we need to supplement these equations by the analytic properties of Ya,s(ζ). As for the AdS3 case, the large θ=logζ behavior can be understood by a WKB analysis. The boundary conditions for the Y-functions are(4.40)logY1,s→-mscoshθ-Cs±Ds+⋯,θ→±∞,logY3,s→-mscoshθ-Cs∓Ds+⋯,θ→±∞,logY2,s→-2mscoshθ+⋯,θ→±∞;
depending on the signature of space-time, the reality conditions of Cs and Ds are fixed. For instance, for (3,1) or (1,3) signature, the Ds are real while the Cs are purely imaginary. As for the AdS3 case, these boundary conditions together with regularity imply the following set of integral equations:(4.41)logY2,s=-ms2coshθ-K2⋆αs-K1⋆βs,logY1,s=-mscoshθ-Cs-12K2⋆β2-K1⋆αs-12K2⋆γs,logY3,s=-mscoshθ+Cs-12K2⋆β2-K1⋆αs+12K2⋆γs,
where(4.42)αs=log(1+Y1,s)(1+Y3,s)(1+Y2,s-1)(1+Y2,s+1),γs=log(1+Y1,s-1)(1+Y3,s+1)(1+Y1,s+1)(1+Y3,s-1),βs=log(1+Y2,s2)(1+Y1,s-1)(1+Y1,s+1)(1+Y3,s-1)(1+Y3,s+1),
and the kernels are given by(4.43)K1=12πcoshθ,K2=2coshθπcosh2θ,K3=iπtanh2θ.
The data of the system of integral equations are the (in general complex) masses ms and the purely imaginary constants Cs. Since for the scattering of N gluons s=1,…,N-5, we get 2(N-5)+N-5=3(N-5) degrees of freedom, which exactly agrees with the expected amount of cross ratios.Again, as for the case of theAdS3 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(4.44)Areg=∑s∫dθms2πcoshθlog((1+Y1,s)(1+Y3,s)(1+Y2,s)2).
This is the final solution for the general problem in AdS5 (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 inAdS. 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 Y-system, valid for any value of ζ. One can combine these functional equations with the knowledge of the analytic behavior of the Y-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.



---

*Source: 703064-2010-06-27.xml*