arxiv: v3 [hep-th] 25 Jun 2009

Transcription

1 UUITP-4/09 Spiky strings in the SL Bethe Ansatz arxiv: v3 [hep-th] 5 Jun 009 L. Freyhult a,, M. Kruczenski b, and A. Tirziu b,3 a Department of Physcis and Astronomy, Uppsala University, P.O. Box 803, S-7508, Uppsala, Sweden b Department of Physics, Purdue University, W. Lafayette, IN , USA. Abstract We study spiky strings in the context of the SL Bethe ansatz equations. We find an asymmetric distribution of Bethe roots along one cut that determines the all loop anomalous dimension at leading and subleading orders in a large S expansion. At leading order in strong coupling large λ we obtain that the energy of such states is given, in terms of the spin S and the number of spikes n by E S = n λ π ln 4πS + ln 4 λ n sin π n + O lns S. This result matches perfectly the same expansion obtained from the known spiky string classical solution. We then discuss a two cut spiky string Bethe root distribution at one-loop in the SL Bethe ansatz. In this case we find a limit where n, keeping E+S, E S n n, J n fixed. This is the one loop version of a limit previously considered in the context of the string classical solutions in AdS 5 S 5. In that case it was related to a string solution in the AdS pp-wave background. lisa.freyhult@physics.uu.se markru@purdue.edu 3 atirziu@purdue.edu

2 Introduction Remarkable progress was achieved in the understanding of the AdS/CFT correspondence by studying the SL sector of the theory. Much of this progress relates to twist two gauge theory operators of the type trφd+ S Φ. On the string side they are described by the folded string solution [] as can also be verified by an alternative computation in terms of the cusp anomaly []. In the field theory a proposed all loop Bethe ansatz [3] can be used to compute their anomalous dimension at all loops in the planar approximation. Computations of quantum corrections to the folded string solution [4] at leading order in large S as well as its generalization to AdS 3 S [5, 6, 7] were recently used to successfully check such all loop asymptotic Bethe ansatz [8, 9, 0, ]. The all loop asymptotic Bethe ansatz computation was then extended to the subleading orders in large S expansion [,, 3] and precise matching with the corresponding string theory computations [4] was found. Beyond twist two, the SL sector contains operators of higher twist which have also been investigated. For example operators of the type O = tr D S n + Φ D S n + Φ...D S n + Φ. are described, on the string side, by the spiky string solutions [5]. There has been some recent progress in this area. In [6] the solutions were shown to correspond to multi-soliton solutions of a generalized sinh-gordon model. This allowed the construction of new, more general, solutions where the spikes move with respect to each other opening up the possibility to study a whole new class of operators. The relation with field theory was examined in detail in [7] where the elliptic curves associated with the classical solution where analyzed and a map was proposed to a similar structure emerging from the study of the field theory operators. However, with the current methods, full understanding of these operators from the field theory point of view requires the construction of the Bethe ansatz solution that describes them. It is the goal of this paper to fill this gap by providing a proposal to describe the spiky string solution in the all loop Bethe ansatz. It was shown in [8, 9] and extended in [7] that such operators can be described by a spin chain with a number of sites n being the same as the number of spikes. At leading order in large S the spiky string solution touches the boundary of AdS 3. It was shown [0] that at this order the spiky string solution can be mapped to the folded string solution by a conformal transformation, and thus the all loop energy at leading order in large S is given by E S = n fλ lns + OS0, n. where fλ is a universal scaling function known as the cusp anomaly. In this paper we obtain this result from the field theory side by an all loop Bethe ansatz computation. In a semiclassical analysis in string theory the spiky string solution was extended to AdS 3 S in []. It turns out that, in the limit of large number of spikes n, there is an interesting scaling limit in which E+S, E S and J remain fixed. It was shown in [] that the solution in n n n this limit is related to a solution in the AdS pp-wave background introduced in [0]. In this paper, by using the found Bethe ansatz solution, it is shown that the same pp-wave like limit can be taken at weak coupling.

3 To construct the Bethe ansatz solution, an important input is a set of integers n k which may be interpreted as the bosonic mode numbers of the waves propagating in the spin chain. Different such sets of numbers give rise to different solutions of the Bethe ansatz equations. The folded string solution in AdS 3 reduces [] to the folded string solution in flat space in S the limit λ. In flat space we can quantize the string exactly in terms of right and left moving waves propagating on the string. It is easily seen that the folded string has the same number of right and left moving excitations, i.e. n L = n R, with the same wave numbers k L = k R =. Extending this idea to AdS 3 for all values of spin S and coupling λ, suggests that in a Bethe ansatz computation the folded string solution is described by two sets of modes, one with n k = and the other with n k =. These numbers were used then as the input in the one loop SL Bethe ansatz equations for two cuts corresponding to n k = ±, and a symmetric Bethe root distribution was found to describe the folded string solution [3]. The same bosonic quantum numbers were used later [4, 3] to find the corresponding solution to the all loop Bethe ansatz. In the large S limit and finite J the two cuts in fact merge into one cut with a discontinuous wave-number distribution since we still have n k = ±. Following the same procedure, we look at the spiky string solution in flat space, find the left and right wave numbers and use them as an input to find the corresponding Bethe root equation that we then proceed to solve. The spiky string solution was already discussed in [5]. It turns out that the bosonic quantum numbers that should be used in the Bethe ansatz equations are n L = and n R = n. Using this input we find an asymmetric distribution of Bethe roots spread over one cut along the real axis at one loop in weak coupling and leading order in the large S expansion. We then extend the computation to all loops and further obtain the all loop result also for the subleading order in large S. More precisely, we obtain that the energy of the spiky string at all loops can be written as E S = n + fg n ln S + γ E + ln n sin π n + n B g + O ln S S, n, g = λ 6π.3 where B g denotes the virtual scaling function of twist operators obtained in []. The leading large S term, i.e. E S fg n ln S, matches the string theory result. at the same order. If we further expand the all loop result.3 for large coupling we obtain E S = n λ ln 4πS 4 + ln π λ n sin π.4 n + n + 6 ln 3 ln 4 π π ln n sin π 3 ln 4πS ln + O, ln S n π λ λ S Remarkably, the first line precisely matches the result known from expanding the classical string solution [4], that is, not only at order lns but also S 0. The second line is a prediction for the quantum corrections to the classical solution which will be interesting to check explicitly. Let us observe that, since the asymptotic Bethe ansatz reproduces the correct strong coupling result to order S 0, wrapping effects play no role to this order in S. Having understood the solution with one cut we proceed to consider a two cut solution to the one-loop Bethe ansatz equations. Again we use the spiky string bosonic quantum numbers 3

4 for the integers n k, i.e. one cut has n k = and the other n k = n where n is the number of spikes. The solution has four real parameters d < c < a < b which represent the position of the two cuts given by the segments [d, c] and [b, a]. The solution obtained allows us to compute the energy E, spin S, R-charge J and number of spikes n in terms of these parameters. Although there are four parameters it turns out that, for the solution to exist, a consistency condition has to be satisfied which implies that only three are independent. Equivalently this means that we can, implicitly, express the energy as a function of the other physical quantities E = ES, J, n. As a check, in the limit c, b 0, namely when the cuts merge, we recover the results of the previous one cut solution. More interestingly, in the limit d a pp-wave type scaling is obtained, i.e. n, with E+S, E S and J remaining fixed. In this limit the equations simplify since one parameter n n n is removed. To be able to compare with the string calculation of [] it would be necessary to extend this results to all loops and then consider the strong coupling limit. This should be an important test of the all-loop Bethe ansatz since these solutions contain a lot of structure and a highly non-trivial dependence of the energy with the other parameters. Another direction one can pursue is to consider the Bethe ansatz description of the more general solutions described in [6]. In fact at this point it might be more interesting to find a direct mapping between the classical action of the string and the classical action of the spin chain at strong coupling in a similar way as in [5]. Mapping the actions would determine a mapping between the configurations of the string and those of the spin chain. In particular this would allow to understand how to see the spikes and their motion from the spin chain/field theory point of view. In the approach we follow here this cannot be seen since we map a semiclassical or coherent state on the string side to an eigenstate of the Hamiltonian Bethe ansatz. A more detailed map would have to be constructed between coherent states on both sides. The paper is organized as follows. In section we review the spiky string solution in flat space which leads us to the correct mode numbers to be used in the Bethe ansatz. In section 3 we present in detail the construction of the one-loop Bethe root distribution for a one cut solution at leading order in large S. In section 4 we extend the computation to all loops using the all loop asymptotic Bethe ansatz. Section 5 deals with the extension of the all loop solution at leading order in large S to the subleading order OS 0. In section 6 we analyze the spiky string solution with two cuts using the one loop Bethe ansatz equations. We then take two limits: first a limit when the two cuts collide recovering the one cut solution, and second, we consider the pp-wave type limit. In appendix A we present some details about the solution of the Bethe ansatz equations which we used throughout the paper. In Appendix B we present a review of the OS 0 corrections to the ground state, which corresponds to the folded spinning string solution. Spiky strings in flat space As with the folded string, it is convenient to start by studying solutions in flat space. In that case we can quantize the theory exactly and it turns out that the solutions are just a 4

5 superpositions of a left and a right moving waves: x = A cos n σ + + A n cos σ. y = A sin n σ + + A n sin σ. t = A n τ = A n σ + + σ.3 where σ + = τ + σ, σ = τ σ, n is the number of spikes and A is constant that determines the size of the string. Here, τ, σ parameterize the world-sheet of a string which is moving in a Minkowski space with metric: ds = dt + dx + dy.4 The solutions are periodic in σ with period π and satisfy the equations of motion in conformal gauge, τ σx µ = 0, as well as the constraints + X = X = 0. Quantum mechanically the state has n R = A n right moving excitations of wave number k R = and n L = A n left moving excitations with wave number k L = n satisfying the level matching condition n R k R = n L k L. All excitations carry one unit of angular momentum and therefore the total angular momentum and energy are given by S = n L + n R = A nn, E = n n L k L + n R k R = An, E = n S.5 which agrees with a classical computation. For n = we recover the standard Regge trajectory E = S and for n > we get a Regge trajectory of modified slope. At fixed time, the shape of the string for different values of n takes the form depicted in figs. and. It can be seen that the string has n spikes or cusps and, analyzing the time dependence, that it rotates rigidly in such a way that the end points of the spikes move at the speed of light. More details on these solutions can be found in [5]. Here, we will just need the fact that the spiky string is obtained by superposing left and right moving waves of different momentum. This suggests that when extending to all loops we should take a distributions of Bethe roots such that the wave number can take the values n w = or n w = n where n is the number of spikes. 3 Spiky strings in the Bethe Ansatz: -cut solution We follow the same prescription used in [4] to obtain the lns dependence from the Bethe ansatz at large S. In the SL sector the Bethe roots are on the real axis and, for large S, they accumulate on various segments or cuts. In the case of the folded string with J = and S, there is only one cut and the root distribution is symmetric. In this case, namely for the lowest twist, there is only one state whereas for J > there is more then one state. Nevertheless, for the lowest energy state at higher twist, the root distribution is again symmetric as discussed In the case when S J there are two cuts. This situation was analyzed in [3] and will be studied later in this paper. 5

6 Figure : Rotating string with 3 spikes in flat space. We plot x, y in eq..3 parametrically as a function of σ = 0 π for n = 3, A = and τ = 0. in [4]. In this section we consider excited energy states which describe the n-spike solution in the large S limit for a fixed finite J. Following the discussion of the previous section, in the case of the spiky strings at leading order in large S, we assume to have one cut but with an asymmetric root distribution. The mode numbers of each excitation will be taken to be n w = or n w = n where n is the number of spikes. Let us start with the one loop SL spin chain. The one loop Bethe ansatz equations corresponding to a XXX nearest neighbor spin chain can be written as uk + i J = u k i S j k u k u j i, k =,,..., S 3. u k u j + i Taking the logarithm of 3., and further expanding in u as appropriate for large S gives [4] J u k = πn k S j=,j k u k u j 3. where u k are the Bethe roots and J is the length of the corresponding spin chain. The total 6

7 Figure : Rotating string with 0 spikes in flat space. We plot x, y in eq..3 parametrically as a function of σ = 0 π for n = 0, A = and τ = 0. momentum should vanish due to the cyclicity of the trace. This condition reads S u j + i u j i j= = 3.3 The one loop energy is given by E S J g = S k= u k At this point one should solve eq.3. for the real numbers u k subject to the condition 3.3. However, when the number of roots S is very large, the roots accumulate in cuts on the real axis and one can conveniently approximate the problem by defining a density of roots. Thus we introduce the root distribution function ρ 0 u = k δu u k. Furthermore, for large S, the cut grows as S and we can make the approximation that u. Therefore, at leading order in 7

8 Figure 3: Wave numbers n w ū = n w u for n=,3,6. By analogy with flat space, we take the S modes with ū < 0 to have n w = and those with ū > 0 to have n w = n. the large S expansion 3., we need to solve the equations du ρ 0u u u = πn wu 3.5 duρ 0 u = S 3.6 du ρ 0u u = 0, 3.7 which are, the Bethe equation, the normalization of the density, and the zero momentum condition. The function ρ 0 u has support on a finite interval that needs to be determined as part of the solution. To solve these equations we use the method described in the appendix. We first define the cut, namely the interval where ρ 0 is non-vanishing, to extend from d < x < a with d < 0 and a > 0. The function n w u is defined as n w u = {, d < u < 0 n, 0 < u < a 3.8 In figure 3 we plot such functions in terms of a rescaled ū = u. We then introduce a function S Fw = w a w d which has a cut on the real axis extending precisely from d to a. Furthermore, right above below the cut it is purely imaginary with positive negative imaginary 8

9 part. Using this function we define the resolvent as Gw = π a d Fw n w u du 3.9 Fu + i0 + u w which, following the appendix A, determines ρ 0 u to be ρ 0 u = a π Fu n w u Fu u u du 3.0 Besides, we need the behavior of Gw at infinity which is readily seen to be with G 0 G = π d Gw G 0 + G +..., w 3. w a d = a + d G 0 π n w u du = in Fu + i0 + a d + in a + d π arcsin a d 3. un w u du 3.3 Fu + i0 + As explained in the appendix we need to satisfy the consistency condition [ ] Gw Res w u, = G 0 = for ρ 0 u to solve eq.3.5. On the other hand, equation 3.6 gives a S = ρ 0 udu = Gwdw = iπres [Gw, ] = iπg 3.5 and eq.3.7 is d a ρ 0 u 0 = d u du = C C [ ] Gw Gw dw = iπres w w, = iπg which is satisfied if 3.4 is. In the previous equations we used a contour C that encircles the cut in a clockwise direction. Using 3. we then find a π d = tan 3.7 n Using 3.3 and 3.5 gives S = n ad 3.8 Before computing the energy let us compute the generating function G Gw = Fw iπ a d n w u du u w a uu d = i ± n aw d + da w π ln 3.9 da w aw d 9

10 where ± refer to w having positive or negative imaginary part. w takes values in the complex plane with a cut along the real axis from d to a. We now need to compute the energy which, within the Bethe ansatz is given by E S J a ρ 0 u = g d u + du = Gw 4 C w + dw { [ Gw = iπ Res w +, i ] [ Gw + Res w 4 +, i ]} = π G i G i 3. 4 Using the definition of Gw and through an explicit computation we find E S J g = πreg i F i = Im n w u Fu u i du 3. a a i d + i d = n ln i 3.3 a d For fixed n and large S we see from eqs.3.7 and 3.8 that a d S so we can approximate the last equation as E S J = n ln 8ad, S 3.4 g a d Therefore, through eqs.3.8, 3.4, we have solved the problem by writing E and S in terms of two parameters a and d related by the condition 3.7. Explicitly we can write E S J g a = S n tan π n d = S n cot π n = n ln 4S n sin π n Notice that up to now we have not needed the actual density ρ 0 u. However, in the following section it will be needed for extending the result to all loops and subleading terms in the large S expansion. Using eq.3.0 it can be written as ρ 0 u = n π ln au d + da u au d da u 3.8 = { n arctanh d a u π a u d, a u > 0 n u d arctanh a π d a u, d u < Plots of the root distribution are shown in Figure 4. It is clear that S can be thought as 0

11 Figure 4: Root density ρ 0 ū = ρ 0 u for n =, n = 3 and n = 6. The larger the n the more S skewed the distribution becomes. determining the scale of the variable u. For that reason it is convenient to define a rescaled variable ū = u. Defining new parameters S ã = n cot π n, b = n sin π n 3.30 such that d S = ã + b, a S = b ã we find another convenient way to write the density: ρ 0 ū = n b π ln ū + ãã + b ū + ã b ã 3.3 b ū + ãã b ū + ã b ã Summarizing, the result from this section is that the energy of the spiky string at one-loop is given by E S g = n lns + OS When n = this indeed recovers the twist two result in [4]. The leading order expression is proportional to n similarly as at strong coupling.

12 4 Spiky strings in the all loop Bethe ansatz After understanding the root distribution at one-loop we now proceed to find a solution to the all-loop Bethe equations. We will follow the procedure used in [4]. The asymptotic all loop Bethe equations in the SL sector read J arctanu k + ij ln + g /x k + g /x + k = πñ k + i S/ j= S/ j 0 S/ j= S/ j 0 arctanu k u j 4. ln g /x + k x j g /x k x+ j S/ j= S/ j 0 θu k, u j, where θu, v denotes the dressing phase and x ± u are defined through u ± i = x± u + g x ± u. For the ground state the fermionic mode numbers are given by [4] ñ k = k + J ǫk, k = ±, ±3,...,±S, 4. where k = k ǫk/. For the highest excited state there is no gap in the mode numbers [8] In the limit of large S we define the density of roots as ñ k = k, k = ±, ±3,...,±S. 4.3 ρu = S dk /S du and the all loop equations can then be written as J u + /4 + ij d du ln + g /x u d a du ρu d du d + g /x + u = πρu a i ln g /x + ux u g /x + u x u θu, u ρu u u At one-loop this reduces to J u + /4 = πρ 0u a d du ρ 0 u u u Rescaling the roots, u = Sū, we find in the large S limit J b ã S πδu = π ρ 0ū dū ρ 0 ū πδū ū + b ã S P ū ū 4.7

13 and the one-loop equation reduces to 0 = πjδū + b ã b ã dū ρ 0 ū ū ū. 4.8 For the highest excited state J is equal to the number of cusps, n, and we recognize the above equation as the derivative of the one-loop equation. To simplify the presentation of the all-loop computation we will omit the dressing phase and only restore its contribution in the end. We can now proceed with the all loop equation by splitting off the one-loop piece of the density, ρu = ρ 0 u + g σu, σu 0 = π σu u u + + a du ρ g 0 u d i ln g /x + ux u d du g /x + u x u + du σu d i ln g /x + ux u. 4.9 du g /x + u x u In the above we have extended the integral boundaries to ± in the terms containing the higher loop density, σu. For the term containing the one-loop density we find, in the large S limit, a du ρ g 0 u d i ln g /x + ux u d du g /x + u x u a = i du ρ 0 u d d du x + ux u x ux + u = π d G i du x + u G i +..., 4.0 x u Using 3.9 we find and G ± i 0 = π σu + n ln S d du x + u + + du σu d du n = ± lns π σu u u + x u i ln g /x + ux u g /x + u x u After Fourier transformation and a redefinition of the density,. 4. ˆσt = ne t/ˆ σt ln S, 4.3 We use the convention ˆft = dw e iwt fw. Details on the Fourier transforms used here can be found in [4]. 3

14 the result is ˆσt = t Kgt, 0 4g dt Kgt, gt ˆσt e t 0 with the kernel defined by 4.4 Kt, t = K 0 t, t + K t, t 4.5 K 0 t, t = n J tt n tj n t n= 4.6 K t, t = nj tt n tj n t. 4.7 Including the dressing phase in the computation replaces the kernel above with where n= Kt, t = K 0 t, t + K t, t + K d t, t 4.8 K d t, t = 8g dt K t, gt t e t K 0gt, t. 4.9 Further we find that the energy is given by 0 E S = 8g nˆσ0 lns + OS which gives E S = n fg lns + OS0 4. which matches the string theory result at leading order in large S [0]. 5 Subleading corrections for the highest excited state The computation of the first subleading corrections for the highest excited state follow the same logic as the corresponding computation for the ground state but uses an asymmetric root density. The subleading corrections for the ground state were computed in [, ] using a method different than what will be convenient here. We will here use the same method as in [4, 3] but also include subleading corrections, for convenience we include the corresponding computation for the ground state in Appendix B. 5. One loop The one-loop equation is given in 4.6, where J should be set to n for the highest excited state. An approximate solution to this equation is given by 3.8 which we expand for large S ρ 0 u = n π ln S n π lnu + n π ln n sin π + O/S. 5. n 4

15 The above density captures only the leading order in the large spin expansion correctly. To get the subleading corrections of order S 0 we split off the above solution from the density, ρ 0 u = ρ 0 u + ru, 5. and solve for ru. We are now allowed to extend the limits of integration to ±. 0 = πru + n ln S n lnu + n ln n sin π n n u + /4 n dv π ln S n π lnu + n π ln n sin π n u v + rv dv u v The equation can be solved by Fourier transformation e t / ˆrt = n e. 5.4 t t Computing the energy including the first two orders in the large S expansion we find a E S n ρ 0 u = g d u + /4 + du ru u + /4 4 = n ln S + n ln n sin π + n n = n ln S + γ E + ln n sin π n 5. All loops The all loop equation can be written as n d xi/ + u ln du xi/ u = πρu a + du ρu d du d a d e t e t / dt e t t. 5.5 ρu u u + i ln g /x + ux u g /x + u x u θu, u. 5.6 Again we will omit the dressing phase, only to restore it in the end. Splitting off the one-loop density ρu = ρ 0 u + g σu, 5.7 where ρ 0 u satisfies 4.6, we find n d xi/ + u ln g du xi/ u u + /4 + du σu d du i ln g /x + ux u g /x + u x u = π σu σu u u + + Iu 5.8 5

16 where Iu = g a d du ρ 0 u d i ln g /x + ux u. 5.9 du g /x + u x u Note that the limits of integration in the integrals containing the higher loop density can be extended to ± to this level of approximation. To obtain the all loop energy, E S n a i = ρu g x + u i, 5.0 x u d we can follow a similar strategy as above. We split off the one-loop density and then extend the limits of integration in 5.0 to ±. It is then clear that we only need the combination σu+ σ u to obtain the energy. To write the integral equation for the symmetric combination σu + σ u we need the combination Iu + I u which we can split into two parts Iu + I u = π G i G i d du x + u x u [ + du ρ 0 u d i ln g /x + ux u g du g /x + u x u d u + /4 du x + u ] 5. x u Note that in the integral in 5. we have extended the integration limits to ±. Fourier transformation of the equation 5.8 gives n J g πe t/ 0 gt = π ˆ σt + ˆ σ t e t + Ît + Î t +4πg te t/ dt ˆ σt + ˆ σ t e t / Kgt, gt, 5. 0 where the kernel is given by 4.5. For Ît we find, to first two orders in the expansion, Ît + Î t = 8πn ln S + γ E + ln n sin π te t/ Kgt, 0 n + 8πte t/ dt n e t Kgt, gt Kgt, Using the result for the one-loop energy 5.5 and introducing the density, ˆσt ˆσt = ˆ σt 8 e t/ + ˆ σ t we hence find for the all-loop equation t [ ˆσt = Kgt, 0 n ln S + γ e t E + ln n sin π n n 8g t J 0gt + dt n 0 e t Kgt, gt Kgt, 0 ] 4g dt Kgt, gt ˆσt

17 Taking the dressing function into account replaces the kernel by 4.8. The string energy can now be written as E S = n + 6g ˆσ We note that the integral equation is very similar to the equation in [] and the solution can consequently be written as E S = n + fg n ln S + γ E + ln n sin π + n n B g where B g denotes the virtual scaling function of twist operators. Using the first orders in the expansions at strong coupling [], we find fg = 4g 3 ln + O/g 5.8 π B g = γ E ln gfg 4g ln 6 ln + 3ln + O/g 5.9 π π E S = ng + n ln S 4 g + ln n sin π n + 6 ln 3 ln π π ln 4 n sin π n 3 ln π ln S + O/g 5.0 g We see that the leading strong coupling result is in agreement with the known string theory result [4]. For n = we recover the result for the folded string [], this result is in full agreement with the one-loop computation from string theory [4, 6]. For arbitrary n 5.0 is expected to be in agreement with the one-loop result computed from the sigma model [7]. Nevertheless, as we mentioned already, it would be interesting to compare the result in the second line in 5.0 with an explicit one-loop string computation of the OS 0 term. 6 Spiky strings in the Bethe Ansatz: -cuts solution We start again with equation 3. where we do not rescale u by S, i.e. we introduce the root distribution ρ 0 u = S k= δu u k du ρ 0u = πn u u w u J u πn wu 6. The condition that the total momentum should vanish gives duρ 0 u ln u + i u i = 0 6. The root density ρ 0 is normalized as duρ 0 u = S 6.3 7

18 The -loop anomalous dimension is E S J g = du ρ 0u u Let us consider two cuts d c b a. Extending the bosonic wave number distribution from the previous section we take n w u = {, d < u < c n, b < u < a 6.5 Then the root distribution function is ρ 0L u = n c u b u a dπ[ b au d π u da u a cb d b du a, r] + u ak[r] + J π u B u B u 6.6 for the left interval u [d, c], and ρ 0R u = n u cu b d aπ[ d cu a π u da u a cb d a cu d, r] + u dk[r] + J π u B u + B u 6.7 for the right one u [b, a]. Here r = a bc d a cb d and we used the following definitions of the elliptic integrals 6.8 K[m] = π 0 dα m sin α, Also, we denote B, B the integrals Π[n, m] = π 0 dα m sin α n sin α 6.9 B u = Fu A 4, B u = Fu A We compute the residue term in A. and obtain [ ] Gw πres w z, = z[na + J π A 4 A 5 ] + A + n A 3 a + b + c + d [na + J π A 4 A 5 ] 6. where A i integrals are defined in the Appendix A and can be written in terms of elliptic integrals. 8

19 As explained in Appendix A we need to take the above residue to zero in order for the root distribution 6.6,6.7 to be a solution of 6.. This gives the following relationship na + J π A 4 A 5 = 0, A + n A 3 = 0 6. Then one can write J, n as functions of a, b, c, d J = πna A 5 A 4, n = A A The normalization condition can be computed using the function G as ρ 0 udu = ires [Gw, ] 6.4 = na + J π A 4 A 5 + A + n A 3 + J π A A 3 + A 6 + n A 7 Taking into account equation 6. the normalization condition reduces to J π A A 3 + A 6 + n A 7 = S 6.5 The consistency of equation 6. gives the equation J du ρ 0u = π duρ 0 unu 6.6 u The momentum condition 6. expanded in large u gives du ρ 0u = u We use again the residues of the function G to compute the integral in 6.7 du ρ 0u = dw Gw [ ] [ ] Gw Gw u w = iπres w, + iπres w, 0 C 6.8 We already computed the first term and set it to zero in 6.. What remains is to compute the second term in 6.8. We obtain du ρu u = abcd A 4 n A 5 J π A 8 + J π A Therefore the momentum condition 6.7 gives the condition abcd[a4 + n A 5 + J π A 8 A 9 ] =

20 All A i integrals depend on the parameters a, b, c, d. Let us write the physical charges in terms of these parameters n = A A 3, J = πa A 5 A 4 A A 3 6. A S = A 6 + A A A 3 A 3 A 4 A 5 A 7, A 3 E S J = πre[g i ] 6. g where we wrote the energy in terms of the function G whose explicit expression in terms of elliptic integrals is in A.0. We also have the condition 6.0 giving a relationship among parameters a, b, c, d A 4 A A 5 + A A 8 A 9 A 3 A = A 3 A 3 A 5 A 4 We observe that using the above equations one can in principle solve for the unknown constants a, b, c, d and obtain the -loop anomalous dimension E = En, S, J. However, this is not possible to find explicitly because of the complicated elliptic integrals involved. Below we consider particular limits. 6. Particular limit: cut solution From the above general -cut solution for an arbitrary number of spikes let us recover the -cut solution again with arbitrary number of spikes. The -cut solution should be obtained in the limit b = c = 0, d < 0, a > To obtain this limit we take b = ǫ, c = ηǫ and take ǫ 0 +. The parameter η > 0 is to be determined. The integrals in this limit become A = ln ǫ ad + Oǫ 0, A = π + arcsin a + d a d, A 3 = π + arcsin a + d a d 6.5 A 8 = A 4 = cot η ad η ǫ + Olnǫ, A 6 = ad π a + d a + d arcsin a + d 4 d a A 7 = ad + π a + d a + d 4 η η cot η adη η ǫ + Olnǫ, A 9 = A 5 = tan η + Olnǫ 6.6 ad η ǫ arcsin a + d d a η + η tan η Olnǫ 6.9 adη η ǫ Then the second equation in 6. implies a + d a d = cos π n

21 which indeed matches the result from the section 3. The first equation in 6.3 becomes η J = πnǫ ln ǫ tan η + cot η The normalization equation 6.5 implies S = n a + d tan π n 6.3 which again matches the results in section 3. Finally, the momentum condition 6.0 determines η to be n η = + π n lnǫ Equations 6.30, 6.3 can be used to find the parameters a, d in terms of S, n. The Bethe root density 6.6, 6.7 in this limit reduces precisely to the -cut density expression Particular limit: pp-wave type scaling We want to consider the limit of large n. To achieve this we take the limit d while keeping a, b, c fixed. The integrals become A = K[q] +..., a c d q = a b a c 6.34 A = π + O, A 3 = [ck[q] + a ce[q]] d a c d A 4 = a a c [K[q] Π[ a a c, q]] + O d d 3/ 6.36 A 5 = a a c Π[a b a, q] +..., A 6 = π d d π a + b + c + O 6.37 d A 7 = 3 a c [a ca + b + ce[q] + ac b + cb + ck[q]] d A 8 = A 9 = a bc a c [ac ae[q] + ba + ck[q] bc + ab + acπ[ a a c, q]] d a bc b [ac ae[q] ack[q] + bc + ab + acπ[a a c a, q]] d The root distributions 6.6,6.7 in the large d limit reduce to ρ 0L u = n π u bu c a u K[q] + Π[ a b a c a u u a u, q] + O d 6.4 ρ 0R u n = π u bu c a u a c a u u K[q] Π[a a c, q] + O d 6.4

22 The number of spikes in this limit is n = π a c d ck[q] + a ce[q] In this limit the physical parameters scale as n, S n, J n, E S n 6.44 which is the same as the pp-wave limit scaling found at strong coupling in [0, ]. More precisely we obtain S = ck[q] + a ce[q] n, J = π a c The condition 6.0 gives the relationship Π[ a b a πak[q] 6.45 Π[ a b a, q] + Π[, q] K[q]n a a c ab + bc + ac, q] = K[q] 6.46 bc The above three equations are to be solved for a, b, c as functions of S, J, n and then plug them in the energy. In the large d limit the function Gw becomes Gw = i w aw bw c ww a where Hw d +... a ce[q] + ck[q]k[q] Π[ a b a, q] Π[, q] a a c 6.47 Hw = wk[q] + w Π[ a b a, q] + Π[ a a c, q] Π[ a b, q] 6.48 a w + K[q] [ a wπ[ a b a, q] + a wπ[ a a b w, q] a + wπ[, q] aπ[a a c a w a c, q]] To find the energy we need to take w = i in Gw. The end of the right cut a is still large, i.e. d a, therefore we can expand in large a and finally obtain the energy at leading order in large n E S J g = n + 4b 4 + 4c 4 sin[ aa c tan b + b tan ]K[q] Π[a c a, q] 6.49 When b, c are also large, i.e d b, c the expression simplifies E S J nb + c = K[q] Π[ a b, q] 6.50 g a c abc a We therefore obtain a set of three equations 6.45, 6.46 which are to be solved for a, b, c in terms of S, J, n. Furthermore, such solution for a, b, c should be replaced in the energy 6.49, 6.50.

23 7 Conclusions We have found the Bethe root distribution that describes the spiky strings in the context of the spin chain model that describes the field theory operators in the SL sector. At one-loop we constructed solutions where the roots condense on one cut and two cuts. The one cut solutions we were able to extend to all-loops. In particular expanding the result for large coupling we see that at order ln S and S 0 in the large S expansion, we recover exactly the result obtained from the classical solutions. The result contains a non-trivial dependence on the number of spikes and therefore is an important test of both, our proposal of describing the spiky strings by this particular Bethe ansatz solution and, more broadly, of the all-loop Bethe ansatz used in the calculation. In fact, the all-loop Bethe ansatz also provides a prediction for the one-loop quantum corrections to the spiky string. This is a doable calculation that should be interesting to perform in order to verify the prediction. In the case of the two cut solution, we did not extend the solution to all-loops. This should be interesting further work since the result at strong coupling is known from the string side. We have found that the equations simplify in a pp-wave-like limit where the number of spikes n grows to infinity keeping E+S, E S and J n n n fixed. We showed that such a limit is well defined at one-loop suggesting that it actually can be taken also at higher loops. Finding the all-loop solution in this limit would be an interesting problem for further studies. Acknowledgments We are grateful to N. Gromov, J. Minahan, A. Tseytlin and S. Zieme for useful discussions. M. K. and A.T. were supported in part by NSF under grant PHY and DOE under grant DE-FG0-9ER4068. L.F. was supported in part by the Swedish research council. Appendix A: Solving the Bethe equations In this appendix we summarize the procedure that we used to solve the Bethe Ansatz equations. Although such methods are well known we adapt them here to our particular needs. However, it remains generic enough to be applied to other situations e.g. at strong coupling.. To be concrete we have to find a function ρx defined on the union of several segments C i=...n on the real axis i.e. cuts such that π ρx C...C n x x = Nx, x C...C n A. where denotes principal part of the integral and Nx is a given real function defined on the cuts. As a first step we have to find a complex function Gz, analytic except for cuts at C i and a pole at infinity. Furthermore, we require that, on the cuts, the real part changes sign and the imaginary part is equal to Nx: Gx ± iǫ = ± ρx + inx, x C...C n, ǫ 0 +. A. 3

24 The, yet undetermined, real part we denote as ρx. It is straightforward to find such a function in general for any Nx. While the procedure following below can also be directly extended for an arbitrary number of cuts, we take for simplicity two cuts. To be specific lets consider two cuts C = [d, c], C = [b, a] with d < c < b < a and the following function: Fw = w a w b w c w d A.3 where the square roots are defined with a standard cut on the negative real axis. In figure 5 we show the phase of Fw for w close to the real axis. Now we can write Gw = Fw Nx π Fx + iǫ x w dx A.4 C C It is easy to verify that Gw so defined has cuts in C and C and satisfies eq.a. with ρx = π Fx Nx Fx x x dx A.5 Furthermore it behaves as Gw wω. Now consider an arbitrary point z away from the cuts and a small contour γ encircling it as shown in figure 6. We have Gz = Gw πi w z dw A.6 Deforming the contour we obtain Gz = Res γ [ ] Gw w z, + Gw πi C C w z dw A.7 where the integral on the right hand side is over contours encircling the cuts clockwise. From the property A. we find [ ] Gw Gz = Res w z, + ρx πi C C x z dx A.8 From this result and A. we find that, if z approaches one of the cuts from above, then: [ ] Gw ρx + inx = Gx + iǫ = Res w x, + ρx πi C C x x iǫ dx A.9 [ ] Gw = Res w x, + ρx + πi ρx C C x x dx, ǫ 0 + A.0 Therefore [ ] Gw Nx = ires w x, + π ρx, C C x x dx A. So we obtain that ρx solves the problem A. under the condition [ ] Gw Res w x, = 0 A. 4

25 + i +i d +i c b i a + Figure 5: Phase of the function Fw defined in the text, when w approaches the real axes. Along the two cuts, the imaginary part changes sign as depicted. C z γ C Figure 6: Contours of integration surrounding z and the two cuts. The residue at infinity can be computed by expanding Gw = G 0 w + G + ω G +... A.3 which results in [ ] Gw Res w x, = G 0 xg A.4 namely we need G 0 = G = 0. Since Gz is already fixed we can only solve this problem for certain positions of the cuts. Namely, G 0 = G = 0 is an equation for a,b,c,d. In that case, formula A.5 gives the solution to the problem A.. Furthermore, it is now straight forward to compute integrals of the type C C ρu u n du = Gw dw = πires C C wn [ Gw w n, 0 ] + πires [ Gw w n, ] A.5 For example ρudu = Gwdw = πires [Gw, ] = iπg C C C C ρu C C u du = Gw C C w dw = πig0 πig 0 ρu C C u z du = Gw dw = iπgz C C w z A.8 A.6 A.7 and so on. In particular the energy E and spin S can be computed by these formulas. Also, the last equation can be used as an alternative definition of Gz [4]. In particular we can 5

26 recover A. by noticing that iπgu ρu ± iǫ = C C u u iǫ du = ρu C C u u du ± iπρu = πnu ± iπρu A.9 The above discussion of two cuts is valid for any function Nu. However, in this paper Nu is simply 6. Nu = n w u J. For such Nu, the function G can be written in terms of πu elliptic integrals as Gw = i [ w aw bw cw d πw aw d n J a cb d πw b aw d a dπ[, r] + w ak[r] n J b dw a πw ] d cw a d aπ[, r] + w dk[r] a cw d + ij w aw bw cw d A π 4 A 5 A.0 w where n, n are integers. In our case of interest n =, n = n. In the remaining of this appendix let us define the integrals used in this paper. First it is convenient to define the function We can now compute A 4 = A 5 = A = A 3 = c d a b A 6 = c d a b Fx = x ax bx cx d A = a b xdx Fx = xdx Fx = dx x Fx = ad a cb d dx x Fx = ad a cb d c d dx Fx = a cb d K[r] [d aπ[ d c, r] + ak[r]] a cb d a c A.3 [a dπ[ b a, r] + dk[r]] a cb d b d A.4 c [a dπ[ad, r] + dk[r]] da c A.5 a [d aπ[db, r] + ak[r]] ab d A.6 x dx Fx = [a cb de[r] a cb d + aa + c + da ck[r] a da + b + c + dπ[ d c a c, r]] A. A. A.7 6

27 A 7 = a b x dx Fx = a cb d [a cb de[r] + ca + c + bc ak[r] c ba + b + c + dπ[ a b a c, r]] A.8 c dx A 8 = d x Fx = a d bc [ada cb de[r] a cb d A.9 + bd a ad c + cd + ac + dk[r] + a dbcd + acd + abc + abdπ[ da c, r]] a dx A 9 = b x Fx = c b ad [cba cb de[r] a cb d A.30 + bd c cb a + ab + ca + bk[r] + c bdab + cab + acd + cbdπ[ bc a, r]] Appendix B: Subleading corrections for the ground state In the main text we compute the subleading corrections in the large S expansion to the density of roots. It is convenient to remind ourselves how this is done in the case of the ground state which is described by a symmetric root distribution. We do this at one-loop and and then make us of that solution to construct the all-loop solution. One loop The one-loop equation for the ground state reads πρ 0 u + πj δu J u + /4 S/ S/ ρ 0 v dv = 0. B. u v + An approximate solution to the equation is given by the following density [8] [4] which we expand for large S ρ 0 u = π ln + 4u /S 4u /S = π ln S π lnu + O/S. B. This density captures the leading order in the large spin expansion correctly. To get the subleading corrections of order S 0 we split the density, ρ 0 u = ρ 0 u + ru, B.3 and solve for ru. We are now allowed to extend the limits of integration to ±. πru + 4 ln S lnu J + πj δu u + /4 dv π ln S π lnv u v + rv dv u v +. B.4 7

28 The equation can now be solved by Fourier transformation e t / ˆrt = e J t t + e Computing the energy we find E S J g = All loops S/ S/ The all loop equation reads ρ 0 u u + /4 + du ru u + /4 e t dt t /. B.5 e t / = 4 lns + 4 ln + dt e t t + e t / = 4 ln S + γ E J ln B.6 J d xi/ + u ln = πρu + πj δu du xi/ u + S/ S/ du ρu d du J S/ S/ i ln g /x + ux u g /x + u x u θu, u e t / ρu u u +. B.7 In the following the dressing phase will be omitted for simplicity. It is easily restored in the end. Splitting the density ρu = ρ 0 u + g σu, B.8 where ρ 0 u satisfies B. we find J d xi/ + u ln g du xi/ u u + /4 where + du σu d du = π σ i ln g /x + ux u g /x + u x u σu u u + + Iu B.9 Iu = S/ du ρ g 0 u d i ln g /x + ux u. B.0 S/ du g /x + u x u Note that the limits of integration in the integrals containing the higher loop density can extended to ± as we are considering the large S limit. Further we can split the integral, Iu, into two parts Iu = π G i G i d du x + u x u [ + du ρ 0 u d i ln g /x + ux u g du g /x + u x u d u + /4 du x + u ] B. x u 8

29 where Gu is the one-loop resolvent, defined similarly to 3.9. Here we have used that the root distribution is symmetric around the origin. Note that in the integral in B. we have extended the integration limits to ±. Fourier transformation of the equation B.9 gives J J g πe t/ 0 gt = πˆ σt e t + Ît + 8πg te t/ dt ˆ σt e t / Kgt, gt, 0 B. with the kernel given by 4.5. For Ît we find, to first two orders in the expansion, Ît = 6 ln S + γ E J ln πte t/ Kgt, 0 + 4πte t/ dt J e t / Kgt, gt Kgt, 0. e t + e t / Redefining the density as 0 ˆσt = 4 e t/ˆ σt B.3 B.4 we hence find for the all-loop equation t [ ˆσt = Kgt, 0lnS + γ e t E J ln J 8g t J 0gt + dt 0 e t J Kgt, gt Kgt, 0 e t / + ] 4g dt Kgt, gt ˆσt. B.5 The energy is then given by 0 E S = J + 6g ˆσ0. B.6 The effect of the dressing phase can now be restored by replacing the kernel above with 4.8 and we have hence obtained the same integral equation as in []. Note that this derivation explicitly made use of the fact that the density is symmetric. In the main text we do a similar calculation for the asymmetric distribution we are considering in this paper. References [] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, A semi-classical limit of the gauge/string correspondence, Nucl. Phys. B 636, [hep-th/00405]. [] M. Kruczenski, A note on twist two operators in N = 4 SYM and Wilson loops in Minkowski signature, JHEP 0, [arxiv:hep-th/005], Y. Makeenko, Light-cone Wilson loops and the string / gauge correspondence, JHEP 030, [arxiv:hep-th/0056]. 9

30 [3] N. Beisert, B. Eden and M. Staudacher, Transcendentality and crossing, J. Stat. Mech. 070, P0 007 [arxiv:hep-th/0605]. [4] S. Frolov and A. A. Tseytlin, Semiclassical quantization of rotating superstring in AdS5 x S5, JHEP 006, [hep-th/0046]. [5] S. Frolov, A. Tirziu and A. A. Tseytlin, Logarithmic corrections to higher twist scaling at strong coupling from AdS/CFT, Nucl. Phys. B 766, [arxiv:hep-th/0669]. [6] R. Roiban, A. Tirziu and A. A. Tseytlin, Two-loop world-sheet corrections in AdS5 x S5 superstring, JHEP 0707, [arxiv: ]. R. Roiban and A. A. Tseytlin, Strong-coupling expansion of cusp anomaly from quantum superstring, JHEP 07, [arxiv: ]. [7] R. Roiban and A. A. Tseytlin, Spinning superstrings at two loops: strong-coupling corrections to dimensions of large-twist SYM operators, Phys. Rev. D 77, [arxiv:07.479]. [8] M. K. Benna, S. Benvenuti, I. R. Klebanov and A. Scardicchio, A test of the AdS/CFT correspondence using high-spin operators, Phys. Rev. Lett. 98, [arxiv:hep-th/0635]. [9] B. Basso, G. P. Korchemsky and J. Kotanski, Cusp anomalous dimension in maximally supersymmetric Yang-Mills theory at strong coupling, Phys. Rev. Lett. 00, [arxiv: ]. [0] P. Y. Casteill and C. Kristjansen, The Strong Coupling Limit of the Scaling Function from the Quantum String Bethe Ansatz, Nucl. Phys. B 785, 007 [arxiv: ]. L. F. Alday and J. M. Maldacena, Comments on operators with large spin, JHEP 07, [arxiv: ]. I. Kostov, D. Serban and D. Volin, Functional BES equation, JHEP 0808, [arxiv:080.54]. D. Bombardelli, D. Fioravanti and M. Rossi, Large spin corrections in N = 4 SYM sl: still a linear integral equation, arxiv: [hep-th]. B. Basso and G. P. Korchemsky, Embedding nonlinear O6 sigma model into N=4 super-yang-mills theory, Nucl. Phys. B 807, [arxiv: ]. D. Fioravanti, P. Grinza and M. Rossi, Strong coupling for planar N = 4 SYM theory: an all-order result, arxiv: [hep-th]; The generalised scaling function: a note, arxiv: [hep-th]. F. Buccheri and D. Fioravanti, The integrable O6 model and the correspondence: checks and predictions, arxiv: [hep-th]. N. Gromov, Generalized Scaling Function at Strong Coupling, JHEP 08, [arxiv: ]. Z. Bajnok, J. Balog, B. Basso, G. P. Korchemsky and L. Palla, Scaling function in AdS/CFT from the O6 sigma model, arxiv: D. Fioravanti, G. Infusino and M. Rossi, On the high spin expansion in the sl N = 4 SYM theory, arxiv: [hep-th]. 30

31 [] L. Freyhult, A. Rej and M. Staudacher, A Generalized Scaling Function for AdS/CFT, J. Stat. Mech. 0807, P [arxiv:07.743]. [] L. Freyhult and S. Zieme, arxiv: [hep-th]. [3] D. Fioravanti, P. Grinza and M. Rossi, Beyond cusp anomalous dimension from integrability, arxiv: [hep-th]. [4] M. Beccaria, V. Forini, A. Tirziu and A. A. Tseytlin, Structure of large spin expansion of anomalous dimensions at strong coupling, Nucl. Phys. B [arxiv: [hep-th]]. [5] M. Kruczenski, Spiky strings and single trace operators in gauge theories, JHEP 0508, [arxiv:hep-th/0406] [6] A. Jevicki and K. Jin, Solitons and AdS String Solutions, Int. J. Mod. Phys. A 3, [arxiv: ]. A. Jevicki and K. Jin, Moduli Dynamics of AdS 3 Strings, arxiv: [hep-th]. [7] N. Dorey, A Spin Chain from String Theory, arxiv: N. Dorey and M. Losi, Spiky Strings and Spin Chains, arxiv: [hep-th]. [8] A. V. Belitsky, A. S. Gorsky and G. P. Korchemsky, Gauge / string duality for QCD conformal operators, Nucl. Phys. B 667, [arxiv:hep-th/030408]. A. V. Belitsky, A. S. Gorsky and G. P. Korchemsky, Logarithmic scaling in gauge / string correspondence, Nucl. Phys. B 748, [arxiv:hep-th/060]. A. V. Belitsky, G. P. Korchemsky and R. S. Pasechnik, Fine structure of anomalous dimensions in N=4 super Yang-Mills theory, Nucl. Phys. B 809, [arxiv: ]. [9] V. A. Kazakov and K. Zarembo, Classical / quantum integrability in non-compact sector of AdS/CFT, JHEP 040, [arxiv:hep-th/04005]. [0] M. Kruczenski and A. A. Tseytlin, Spiky strings, light-like Wilson loops and pp-wave anomaly, Phys. Rev. D 77, [arxiv: [hep-th]]. [] R. Ishizeki, M. Kruczenski, A. Tirziu and A. A. Tseytlin, Spiky strings in AdS 3 xs and their AdS-pp-wave limits, Phys. Rev. D 79, [arxiv:08.43 [hep-th]]. [] A. Tirziu and A. A. Tseytlin, Quantum corrections to energy of short spinning string in AdS5, Phys. Rev. D 78, [arxiv: [hep-th]]. [3] N. Beisert, S. Frolov, M. Staudacher and A. A. Tseytlin, Precision spectroscopy of AdS/CFT, JHEP 030, [arxiv:hep-th/03087]. [4] B. Eden and M. Staudacher, Integrability and transcendentality, J. Stat. Mech. 06, P [arxiv:hep-th/060357]. 3

32 [5] M. Kruczenski, Spin chains and string theory, Phys. Rev. Lett. 93, [arxiv:hep-th/0303]. [6] N. Gromov, unpublished [7] N. Gromov and P. Vieira, Complete -loop test of AdS/CFT, JHEP [arxiv: [hep-th]]. [8] G. P. Korchemsky, Quasiclassical QCD pomeron, Nucl. Phys. B [arxiv:hep-th/950805]. 3