The full Quantum Spectral Curve for $AdS_4/CFT_3$

The spectrum of planar N=6 superconformal Chern-Simons theory, dual to type IIA superstring theory on $AdS_4 \times CP^3$, is accessible at finite coupling using integrability. Starting from the results of [arXiv:1403.1859], we study in depth the basic integrability structure underlying the spectral problem, the Quantum Spectral Curve. The new results presented in this paper open the way to the quantitative study of the spectrum for arbitrary operators at finite coupling. Besides, we show that the Quantum Spectral Curve is embedded into a novel kind of Q-system, which reflects the OSp(4|6) symmetry of the theory and leads to exact Bethe Ansatz equations. The discovery of this algebraic structure, more intricate than the one appearing in the $AdS_5/CFT_4$ case, could be a first step towards the extension of the method to $AdS_3/CFT_2$.


Introduction
The idea of a duality between gauge and string theory was put forward many years ago by 't Hooft [1], who noticed that the perturbative expansion in SU (N c ) Yang-Mills theory in the large N c limit naturally organizes in terms of the topology of Feynman diagrams, mimicking the genus expansion of string theory.
The first concrete realization of the duality [2][3][4] conjectures the exact equivalence of N = 4 super Yang-Mills (SYM) theory and type IIB string theory on AdS 5 × S 5 . The precise identification of observables and parameters in the two theories relates the perturbative region of each model to the deep non perturbative regime of the other. For this reason, the correspondence makes powerful predictions, but is also very difficult to test.
An important turning point in this field was the discovery of fingerprints of integrability, at both weak and strong coupling [5,6], in the planar limit of this duality. At least in this limit, it is hoped that the theory will be exactly solved adapting integrable model tools, and remarkable progress has been made on the study of various observables, including Wilson loops and correlation functions.
In particular, the problem of computing the conformal spectrum of the theory was tackled by tailoring integrable QFT techniques to this new setting, in particular the Bethe Ansatz [5,7,8], the TBA, the Y and T-systems [9][10][11][12][13][14][15], leading to the discovery of the very effective Quantum Spectral Curve (QSC) formulation [16,17]. The latter is a very satisfactory simplification and probably the most elementary formulation of the problem. Thanks to the mathematical simplicity of the QSC, it appears that, in the near future, the spectral problem may be completely solved also in a practical/computational sense. Already, the QSC method allows to compute the spectrum numerically with high precision [18,19] and to inspect analytically interesting regimes such as the BFKL limit [20,21] or the weak coupling expansion [22][23][24]. It has also been generalized to so-called γ deformations [25] and to the quark-antiquark potential [26,27].
Another remarkable example of AdS/CFT correspondence was introduced by Aharony, Bergman, Jafferis and Maldacena (ABJM) in [28]. The gauge side of the duality corresponds to the N = 6 superconformal Chern-Simons theory with gauge group U (N ) × U (N ), with opposite Chern-Simons levels, k and −k, for the two U (N ) factors. We will be concerned with the planar limit, where k, N → ∞ with the 't Hooft coupling λ = k N kept finite and the dual gravity theory becomes type IIA superstring theory on AdS 4 × CP 3 . In this regime, integrability emerges, making the ABJM model the only known example of 3d quantum field theory which can be exactly solved [29][30][31][32][33] (see also the review [34]).
The spectral problem in ABJM theory was approached exploiting the experience gained in AdS 5 /CF T 4 . Anomalous dimensions of single trace operators with asymptotically large quantum numbers are described at all loop by the so-called Asymptotic Bethe Ansatz equations, conjectured in [35] and derived from the exact worldsheet S-matrix of [36]. The exact result, including all finite-size corrections for short operators, is formally described by an infinite set of TBA equations, proposed in [37,38]. These equations were solved numerically for a particular operator in [39]. However, solving excited states TBA equations with high precision is a challenging task already for very simple models [40][41][42]. Besides, the form of the TBA equations depends on the state and possibly also on the range of the coupling considered, so that they can be studied only on a case-by-case basis.
It is important to look for a simpler formulation which overcomes these problems. Starting from a precise knowledge of the analytic properties of the TBA solutions [43], the basic equations characterizing the Quantum Spectral Curve of the ABJM model were obtained in [44]. These results were used to compute the so-called slope function in a near-BPS finite coupling regime [45] and to develop a generic algorithm for the weak coupling expansion in the SL(2)-like sector [46].
Although we stress that, as proved by the applications discussed above, the results of [44] contain all the analytic information necessary to solve the spectral problem, several important aspects of the full picture were still missing. First of all, the concrete recipe to describe states within the QSC framework was discussed in [44] only for the SL(2)-like sector. Secondly, the set of equations obtained in [44], the Pµ/Pν-system, can be associated, in the classical limit, to degrees of freedom related to the CP 3 part of the whole AdS 4 × CP 3 target space. A dual system of equations, only briefly mentioned in [44], may be instead associated to AdS 4 classical degrees of freedom. The interplay between the two systems is important for the development of the state-of-the-art solution algorithm at finite coupling [18], as well as at weak coupling for generic states [21,23]. Furthermore, the full algebraic structure was still not transparent, and for example the link between the formulation of [44] and the Asymptotic Bethe Ansatz of [35] was difficult to see. In this paper we will fill these gaps and present the necessary elements for the quantitative solution of the spectral problem for an arbitrary -2 -operator at finite coupling. Besides, we reveal an interesting underlying representation theory structure, which could allow for generalisations and may in particular help in the solution of the spectral problem for AdS 3 /CF T 2 dualities (see [47] for a recent review).
To conclude this introduction, let us review an important fact. In contrast with N =4 SYM, in ABJM theory integrability leaves unfixed the so-called interpolating function h(λ) [30,48], which parametrizes the dispersion relation of elementary spin chain/worldsheet excitations and enters as an effective coupling constant in the integrability-based approach, in particular in the QSC equations. An important conjecture for the exact form of this function, passing several tests at weak and strong coupling [49], was made in [45] by a comparison with the structure of localization results. This conjecture was extendend in [50] to encompass the ABJ model [51], which is based on a more general gauge group U (N ) × U (M ) and possesses two 't Hooft couplings λ 1 , λ 2 in the planar limit. According to the proposal of [50] (based on important observations of [52][53][54][55]), at the level of the spectrum the only difference between the ABJM and ABJ theories lies in the replacement of h(λ) with an explicitly defined h ABJ (λ 1 , λ 2 ) (see [50]). In the following we will simply denote the ABJM/ABJ interpolating function as h.
The contents of this paper are presented in detail below. In Section 2, we discuss the bosonic symmetry underlying the problem, namely SO(3, 2) × SO (6), the isometry group of AdS 4 × CP 3 . We will introduce important vector and spinor notation used in the rest of the paper. Besides, we comment on the interesting fact that the isometry group of CP 3 effectively appears in the Quantum Spectral Curve as SO (3,3), rather than SO (6). In Section 3, we review the results of [44] and discuss how they reflect the CP 3 symmetry. We discuss a subtle modification of the analytic properties (initially overlooked in [44]), which is needed for the study of certain non-symmetric sectors of the theory. The modified equations contain an extra nontrivial function of the coupling, which can be interpreted at weak coupling as the momentum of a single species of magnons. In Section 4, we present an explicit construction of new variables, the functions Q I , Q • and τ i , which satisfy a dual system of Riemann-Hilbert equations reflecting the symmetry of AdS 4 . In Section 5, we treat in full generality the boundary conditions which need be imposed on the solutions of the QSC at large value of the spectral parameter in order to describe a physical state. This is the place where the quantum numbers of the state make an appearance. We also discuss the correspondence between the functions P and Q and quasi-momenta of the spectral curve in the classical limit. In Section 6, based on results obtained in [21,56], we discuss a set of exact relations which are perhaps the most convenient way to repack the analytic properties discussed in Sections 3,4. It is also shown how these equations encode the quantization of the spin. In Section 7, we embed the previous results into a larger set of functional relations which may be considered as (part of) a Q-system. Q-systems are familiar in the theory of integrable models [57,58] and in the ODE/IM framework [59]: they are powerful sets of functional -3 -relations that, supplemented by simple analytic requirements, become equivalent to exact Bethe equations. The structure of Q-systems is completely fixed by symmetry: for example, the QQ relations relations appearing in the N =4 SYM case are the same as the ones for SU (4|4) spin chains. For the OSp(4|6) superalgebra relevant to ABJM theory, however, this algebraic construction was not known in the literature. While we do not treat in full generality the representation theory aspects, we construct explicitly an enlarged set of Q functions, and prove that they satisfy exact Bethe equations reflecting the full supergroup structure. Generalizing arguments of [17], we will show that, in the limit of large volume, some of these exact Bethe equations reduce to the Asymptotic Bethe Ansatz. The paper also contains four Appendices: In Appendix A, we discuss the details of the derivation (already summarized in [44]) of the QSC from the analytic properties of the T-system [43]. In Appendix B, we list some useful algebraic identities used in the derivation of the Q-system relations. In Appendix C, we deduce some of the constraints on the asymptotics of P and Q functions. Finally, in Appendix D we review the dictionary between OSp(4|6) quantum numbers and number of Bethe roots appearing in various versions of the (Asymptotic) Bethe Ansatz, which could be useful for the reader wanting to apply the prescription of Section 5 to concrete states.

Symmetries and conventions
ABJM theory is invariant under the supergroup OSp(4|6), whose bosonic subgroups are associated to the isometries of AdS 4 and CP 3 . We will see that the Quantum Spectral Curve equations encode elegantly this symmetry structure. Let us briefly introduce the main group-theoretic constructions related to the bosonic symmetries.
• CP 3 : the isometry group of CP 3 is the orthogonal group SO(6) SU (4). The invariant 6 × 6 symmetric tensor naturally associated to this symmetry is the metric. This tensor enters the QSC equations 1 , and will be denoted in this paper as η AB . Peculiarly, we will see that it appears in the QSC with a (+ + + − −−) signature. The concrete form of η AB to be used in the rest of this paper is where η AB is the inverse matrix, i.e. η AB η BC = δ C A . This particular choice for η AB emerged naturally from the derivation of the QSC, summarized in Appendix A. As explained there, the specific form of η AB in (2.1) is partly conventional, but its signature cannot be modified without spoiling the reality properties of the system. The fact that the CP 3 symmetry appears effectively as SO (3,3) can be understood heuristically considering the classical limit, where the basic variables of the QSC are related to the quasi-momenta of the algebraic curve (see Section 5.1). The quasi-momenta describing a string moving in CP 3 are defined through the diagonalization of a SO(6) block of the classical monodromy matrix. An SO(2n) orthogonal matrix in general cannot be diagonalized with a real transformation, so that the signature of the metric is not preserved in the eigenvectors basis; moreover, the signature changes precisely to the one typical of SO(n, n).
Let us introduce some conventions. We will use different index labels for objects with different symmetry properties. The indices A, B, C = 1, . . . , 6 will be assumed to carry the vector representation of SO (3,3), and will always be lowered and raised with the metric η AB and its inverse η AB , respectively. It will be useful to consider also spinor representations of SO(3, 3). The relevant 8 × 8 gamma matrices are defined by In even dimension, gamma matrices can always be written in a chiral form: where the matrices σ A ab and (σ A ) ab satisfy While all our equations will be covariant, it is convenient to specify a concrete basis. The matrices σ A ab and (σ A ) ab are defined in our conventions by for an arbitrary vector (V 1 , . . . , V 6 ). Lower-case indices a, b, c will always be taken to run over 1, . . . , 4 and will be reserved for the spinor representations. Note that there is a distinction between upper and lower spinor indices, as they belong to the chiral and antichiral spinor representations, respectively, which are equivalent to the representations 4 and4 of SU (4) SO (6). Another natural tensor that will make an appearance in the equations is the anti-symmetrized product of gamma matrices, • AdS 4 : the isometry group of AdS 4 is SO(3, 2) Sp(4). We will denote the metric of this orthogonal group as ρ IJ , and our concrete choice will be: In the following, we shall always reserve the indices I, J, K, running over 1, . . . , 5, for the vector representation of SO (3,2).
Let us remind the reader of the isomorphism between SO(3, 2) and Sp(4), the group of linear maps preserving a 4 × 4 anti-symmetric two-form. One way to see this is to view SO(3, 2) as obtained from SO(3, 3) by reducing to the subspace orthogonal to a preferred vector v, with v · v = −1.
Then we see that an anti-symmetric two-form naturally emerges: Let us denote a projection of the σ,σ matrices on the subspace orthogonal to v as Σ I ,Σ I , respectively, with I = 1, . . . , 5. By construction, they satisfy the intertwining relations Σ ij showing that there are in fact only five independent matrices Σ I . The latter give a four dimensional representation of Clifford algebra: In the following, we will use indices i, j, k, l, running over 1, . . . , 4, to refer to the fourdimensional representation of SO (3,2). Finally, one can introduce the anti-symmetric combinations which play the role of generators of SO (3,2). By construction, these generators leave invariant the two-form κ ij : therefore the spinor representation of SO(3, 2) is identified with the fundamental representation of Sp(4).
In our concrete case, we see that the metric (2.7) is obtained from (2.1) by restricting to the subspace orthogonal to v = (0, 0, 0, 0, −1, 1). Our choice for the Σ matrices will be 11) and the two-form κ ij reads

Formulation of the QSC from the TBA
In this Section, we recall the first version of the QSC equations proposed in [44]. These equations were obtained through a reduction of the T-system underlying the TBA formulation (see Appendix A), and ultimately take the form of a nonlinear Riemann-Hilbert problem defined on the complex domain of the spectral parameter u. In the u-plane, the Q functions have a characteristic pattern of branch points, whose positions depends on the coupling constant h as specified below. These branch points will all be of square-root type. This peculiar kind of analytic structure for the Q functions, beside AdS 5 /CF T 4 , is also characteristic of some non-relativistic integrable systems such as the Hubbard model [60].

Equations in vector form and analyticity conditions
In the first version of the equations derived from TBA, the basic variables are: six functions {P A (u)} 6 A=1 , and a 6 × 6 anti-symmetric matrix {µ AB (u) = −µ BA (u)} 6 A,B=1 . They are constrained by the following quadratic conditions: where η AB is defined in (2.1). All these functions live on an infinite-sheet cover of the uplane, which, however, is built out of a simple set of rules. On what we will consider the first Riemann sheet, the functions P A (u) have a single branch cut, running from −2h to +2h, see Figure 1. We assume that they have power-like asymptotics at large u, which means that they can be written as a Laurent series in the Zhukovsky variable x(u): The functions µ AB (u) instead display an infinite ladder of branch cuts, at u ∈ (−2h , +2 h) + iZ. They however have the following analyticity property (mirror periodicity 2 ): where the symbol tilde is used throughout the paper to denote analytic continuation around any of the branch points at ±2h (see Figure 1), while the shift on the rhs is evaluated avoiding all branch cuts. Finally, the discontinuities of P A and µ AB across the cut on the real u-axis are related by In addition, as common for the Q functions in integrable models, we should impose a regularity condition for the basic variables P A and µ AB . The precise statement of this condition, however, cannot be formulated in terms of the matrix entries µ AB , but of more fundamental building blocks which we introduce below. Figure 1. Cut structure of the P A functions, with a single cut on the first sheet. We denote with P A the analytic continuation to the next sheet, through the cut on the real axis. Figure 2. The quasi-periodicity property of ν a functions on a sheet with long cuts corresponds to ν a (u + i) = e −iP ν a (u) on the defining sheet with short cuts.

Equations in spinor form
As already discussed in [44], the matrix µ AB can be decomposed in terms of 4 + 4 functions ν a , ν a , as 3 (3.5) which, using the sigma matrices introduced in Section 2, can be compactly written as The constraint (µη) 2 = 0 is now equivalent to the condition ν a ν a = 0. (3.7) As suggested by the analysis of the weak coupling limit [44,46], we will require that the functions ν a , ν a are entire functions on any sheet of the Riemann surface, bounded as u approaches any of the branch points at ±2h + iZ, and with power-like asymptotics at infinity. Under these conditions, the splitting (3.6) contains nontrivial analytic information, and may be argued to be essentially unique 4 . The new functions ν a and ν a should therefore be regarded as more fundamental objects than µ AB . Indeed, at weak coupling, ν 1 and ν 4 are proportional to the Baxter polynomials containing the two types of momentum-carrying roots entering the 2-loop Bethe Ansatz of [29].
The weak coupling analysis also reveals that the periodicity of µ AB on the mirror sheet, equation (3.3), in general translates into quasi-periodicity for the basic functions ν a , ν a (see Figure 2). In the subsector considered in [46], these functions could be either periodic or antiperiodic, and this is a general characteristic of a large sector of states discussed in Section 4.4. For a completely generic state, however, we have 5 where the phase P depends on the state under consideration and may be, in general, a nontrivial function of the coupling constant h. We will make more comments on this quantity in Section 3.3 below.
It is now convenient to pack the six P functions into an anti-symmetric 4 × 4 tensor P ab , defined as while the inverse matrix reads (3.10) The constraint (3.1) can now be rewritten as the condition that P ab has unit Pfaffian: Besides, it is possible to verify that the discontinuity equations (3.4) can be split nicely as It is unique apart for trivial rescalings νa → νa z, ν a → ν a /z, where z is a constant independent of u.
This freedom is however removed by the choice of the normalization of equations (3.12),(3.13) below. 5 Notice that P has to be the same for all the components of νa, due to the fact that in (3.5) all combinations of νaν b are present, for every a, b.
-9 -As discussed in [44], in this form the equations are, from a purely algebraic point of view, exactly the same as the Pµ-system of N = 4 SYM [16,17], with the redefinitions ν a → (P a ) SYM , ν a → (P a ) SYM , P ab → (µ ab ) SYM . (3.14) The analytic properties characterizing the AdS 5 /CF T 4 case are however completely different: the map between the two models in (3.14) requires to change all periodic functions into singlecut functions, and viceversa 6 . Equations (3.7),(3.11), (3.12) and (3.13) should be supplemented with the requirement that all functions are bounded and free of singularities on every sheet of the Riemann surface, and with some information on their large-u asymptotics, see Section 5. This set of conditions is in principle already constraining enough to determine the spectrum, but it is difficult if not impossible to solve in practice at finite coupling. For this purpose it is necessary to embed them in the wider set of equations derived in Sections 4 and 6.

Interpretation of the phase P at weak coupling
The phase P appearing in (3.8) has an interesting interpretation at weak coupling. Recall that the ABJM spin chain admits two types of momentum-carrying excitations [28,36], also known as A and B particles and corresponding to excitations of type 4 and4 in our notations. These pseudoparticles satisfy collectively the zero momentum condition: The total momentum of a single type of excitations is instead in general a nontrivial function of the coupling: it can be defined in the regime of validity of the Asymptotic Bethe Ansatz as P (4) We will show that, at the first two weak coupling orders, In particular, since A and B particles are decoupled at weak coupling, this shows that P is quantized in units of the spin chain length L at leading order: P + O(h 2 ) ∈ 2πZ L . At order O(h 0 ), the identification (3.17) follows from the condition that the ν a (u) functions have a finite limit as u approaches one of the branch points. As discussed in [22], at the leading weak coupling order this requirement boils down to Using the weak coupling limits where Q 4 and Q4 are the Baxter polynomials storing 4 and4-type Bethe roots, respectively, equations (3.8), (3.18) lead to: Further, in Section 7.3, we derive an explicit expression for P for finite h in the large volume limit -equation (7.88) -which proves the identification (3.17) up to the next order at weak coupling. For a generic short operator at finite coupling, the above mentioned large-volume result is not applicable, and therefore P is in principle an undetermined, state-dependent function of the coupling. This could raise some questions on the completeness of the system of QSC equations. It is part of our proposal that P should not be seen as an input, but is rather fully fixed, for every state, from the self-consistency of the QSC. In particular, we expect that it is not necessary to know P in advance in order to compute the anomalous dimension numerically using the method of [21] 7 . On the contrary, this phase can be computed as an output from the solution of the QSC (for instance, one method to reconstruct P is presented in Appendix E). It would be interesting to clarify whether this quantity has any meaningful physical interpretation at finite h.

Construction of the AdS -related Q functions
As we will discuss in Section 5.1, the equations presented above are associated, in the classical limit, to the CP 3 degrees of freedom, and in particular the P A functions are quantum versions of the classical quasi-momenta living in this part of the target space. We shall now show how to construct an equivalent version of the QSC which is more appropriate to the description of AdS 4 degrees of freedom, and contains, in the classical limit, the four quasimomenta parametrizing the motion of a classical string solution in AdS 4 . As in the case of AdS 5 /CF T 4 considered in [17], this entails a swap between the physical and the mirror section of the Riemann surface. In addition, we will see that this alternative system naturally encodes the relevant symmetry group SO(3, 2), which was not explicitly visible in the previous formulation.

The Q a|i and Q ij functions
It is convenient to introduce the standard notation for shifts of the rapidity variable u: where we will always assume that shifts are performed on the section of the Riemann surface where all cuts are short. 7 We plan to return on this issue shortly [56].
-11 - The first step of our construction is the definition of a 4 × 4 matrix Q a|i , through the 4th order finite difference equation Notice that exactly the same equation is satisfied by ν + a , as can be verified by combining (3.8) and (3.13): and that the index i in (4.2) does not enter the matrix structure of the equation. We will take this index to run from 1 to 4, labeling a set of independent solutions of this fourth-order equation, distinguished by different asymptotic behaviours at large u (see Section 5). Despite the fact that they satisfy the same finite-difference relation, the analytic properties of ν a and Q a|i will be different: we shall require that Q a|i (u) has no singularities in the whole region Im(u) > 0. Notice that, because of the cut of P ab on the real axis, (4.3) implies that Q a|i has an infinite ladder of short branch cuts in the lower half plane, starting at Im(u) = −1/2. It will be convenient to define Q a |i ≡ (P ab ) − (Q b|i ) [−2] , so that (4.2) can be split as Now, let us construct the tensor Using (4.4), it is simple to see that k ij is invariant under a shift u → u + 2i, and, since by construction it is free of cuts in the upper half plane and has power-like asymptotics, it must be a constant matrix. In addition, notice that (4.4) implies more precisely that k + ij = −k − ji , so that k ij is an anti-symmetric matrix, i.e. a symplectic form. This shows that the space of the i-indices should be thought as carrying the fundamental representation of Sp(4) SO (3,2), the isometry group of AdS 4 . It is suggestive that this symmetry was completely hidden at the level of the equations discussed in Section 3, and the mechanism by which it has appeared evokes somehow a spontaneous symmetry breaking.
From (4.5) we see that the specific form of k ij can be adjusted by taking different linear combinations of the columns of the matrix Q a|i (we are allowed to do this since the defining relation (4.2) is linear). We use this freedom to impose that k ij = κ ij as defined in (2.12). Note in particular that 8 Pf(κ ij ) = −1.
Using (4.5), we can relate Q a |i to the inverse transposed matrix of Q a|i : Another simple consequence of (4.2) is that the determinant det(Q a|i ) is invariant under shifts of +2i; by the same arguments as above, it also must be a constant independent of u. Considering the Pfaffian of equation (4.5) and using the property Pf(A t BA) = det(A) Pf(B), we see that We proceed now to construct an object whose indices live in the product of two Sp(4) representations, as Let us discuss the algebraic properties of this tensor. First, from (4.8), we see immediately that Being a 4 × 4 anti-symmetric matrix, Q ij has six independent components. It will be convenient to decompose it into 5+1-dimensional irreducible representations of SO(3, 2) using the invariant tensor κ: the trivial representation is given by the trace while the five dimensional vector representation is the traceless part: Finally, the inverse matrix Q ij , satisfying Q ij Q jk = δ k j , can be computed as and it is simple to show (see Appendix B.3) that the following identity holds (4.14) The following relations constitute a natural counterpart of (4.4) involving the Sp(4)-invariant indices: Shortly, we will show that the elements Q ij have very simple analytic properties: starting from the upper half plane, they can be analytically continued to a Riemann section with the only branch cuts being the semi-infinite segments (−∞, −2h) and (2h, ∞).

The τ i functions
We now construct a new set of four functions, denoted as τ i and defined as Manifestly, these quantities exhibit an infinite series of short branch cuts. Applying (4.4) and (3.8), we see that, under a shift u → u + i, they transform as 17) and shifting this expression once more we find that τ i are 2i-periodic on the Riemann section with short cuts: τ The τ i functions may be seen as counterpart of the ν a functions. Their analytic properties are very similar, with a characteristic swap of short and long cuts. However, notice that, while the functions ν a and ν a are distinct objects, carrying different irreps of SO (3,3), there are only four independent functions τ i , corresponding to the spinor representation of SO(3, 2).

The Qτ -system
The functions Q ij (u) introduced above have, by their very definition, no singularities in the upper half plane, with two branch points at u = ±2h and an infinite ladder of short cuts further down in the lower half plane. Let us study the analytic continuation of Q ij and τ i through the branch cut on the real axis. Combining (4.18) and (4.17), we have and, since Q a |i has no cuts in the upper half plane, we find where we used (4.15) in the last step. By comparison with (4.19), we see that (4.20) can be rewritten as Let us now consider the discontinuity of Q ij : we find   All in all, we see that the discontinuities (4.20) and (4.22) take the form The second relation in (4.23) shows how the phase P appears in the Qτ -system, through (4.21). Finally, contracting (4.16) and (4.17) with κ ij , we find the constraint Equations (4.23), with the constraints (4.24), (4.9) may be considered as a counterpart of the Pν-system (3.7),(3.11)-(3.13). While the equations take a very similar form, they are not identical from the algebraic point of view, due to the fact that the functions τ i and τ i are simply related, for a generic state, by a shift in the spectral parameter, as expressed by (4.21). This distinction reflects the representation theory, as there is only one four-dimensional representation of Sp(4). The difference can be fully appreciated by projecting the Qτ equations on irreducible representations; this is discussed below in Section 4.3.2.

Q ij on the mirror sheet
Let us now prove that, when analytically continued from the upper to the lower half plane passing through the cut (−2h, 2h), the matrix Q ij is analytic in the whole lower half plane -15 -(see Figure 3). Therefore, on an appropriate Riemann section, it has only a pair of long cuts stretching from ±2h to infinity (see Figure 4). This is a very strong analogy with the AdS 5 /CF T 4 case considered in [16]. We start by observing that, using (4.24) and the second equation in (4.23), the discontinuity relation (4.22) can be put in the form This relation can be recast as We will now show that Q LHPA a|i has no branch cuts in the lower half plane (hence the superscript LHPA -Lower Half Plane Analytic). Therefore, the representation (4.26) manifestly shows that the same is true for Q ij , and therefore Q has a single long cut on the mirror Riemann sheet.
To prove that Q LHPA a|i has no cuts in the lower half plane, we can exploit the fact that, due to the periodicity of f j i (u), it satisfies the same fourth order difference equation (4.2) fulfilled by Q a|i . Therefore, it is sufficient to check that it has no cut on the lines Im(u) = −1/2, −3/2: the difference equation (4.2) will then automatically imply that it is analytic everywhere in the lower half plane. This leaves us with just two conditions to check. The first discontinuity to study is where we are using the notation (G) = G − G. From the first relation in (4.15), we find where we used (4.23) in the last step. We may now to use the following identities, found by inverting (4.19), (4.20): to transform (4.30) into The last equality shows the vanishing of the discontinuity (4.28). A completely analogous calculation would show that therefore also the next discontinuity is trivial which concludes the proof.

Vector form of the Qτ -system
We may rewrite the discontinuity equations (4.23) in an alternative form, more similar to the Pµ-system. To do this, let us rearrange the components of Q 5 ij into a five-vector: where we are using the matrices Σ I and the metric ρ IJ defined in Section 2. In components, this definition reads It is also convenient to define or explicitly: From (4.18),(4.21), it is simple to prove that the components of ω IJ (u) are i-periodic functions, while the components of ψ I are anti-periodic under the same shift: In terms of these new variables, the nonlinear constraints (4.9),(4.24) take the form while the discontinuity equations (4.23) can be rewritten as

Reduction to 4 ↔4 symmetric states
In this Section we consider the reduction of the QSC equations to a large subsector characterized by perfect symmetry between the contributions of A-and B-type excitations. In terms of the ABA, this subsector is characterized by the equality of the sets of momentum-carrying Bethe roots, k=1 . As discussed in Appendix A, this case is selected by the conditions: In this case we have the relation P ab = κ al P lm κ mb and we see that necessarily, e iP is either 1 or −1. By studying the large-u asymptotics of equation (4.2), we find that, in this case, the elements of the matrices Q a|i , Q a|i may be chosen as related by the symmetry: This means also that The symmetry imposes the following condition: Taking (4.44),(4.45) into account in (4.17), we see that in this subsector the periodicity of τ i is enhanced to τ which means that τ 1 and τ 4 are i-periodic, while τ 2 , τ 3 are i-anti-periodic. Since we expect all these functions to have power-like asymptotics for physical operators, we see, from the condition of anti-periodicity, that This resut will be important in the following. Finally, in terms of the variables of Section 4.3.2, the reduction to the symmetric subsector can be obtained setting Q 5 = ψ 5 = ω 5I = ω I5 = 0.

Asymptotics and global charges
The Riemann-Hilbert type equations described in Sections 3 and 4 have to be supplemented with appropriate constraints on the large-u behaviour of the functions entering the QSC. We will assume, in analogy with [17], that all the functions we have described scale as powers of u for large values of the spectral parameter, in particular An important observation is that, since the P functions have a single short cut on the first Riemann sheet, they must have trivial monodromy around infinity, which forces M A ∈ Z. These integer parameters are related to the three SO(6) R-charges J 1 , J 2 , J 3 , corresponding to three angular momenta parametrizing the motion of the string in CP 3 , through the prescription The AdS 4 charges ∆ and S, corresponding to the conformal dimension and spin of the gauge theory operator, respectively, enter the QSC through the asymptotics of ν a , described below. Equivalently, they can be read off the coefficients A A in (5.1), through the constraints (with no summation implied on the index B), where the 5-vectorM is defined aŝ The charges (∆, S, J 1 , J 2 , J 3 ) used above are defined relatively to the Dynkin diagram of Figure 5. We remind the reader that, for supersymmetric algebras, the definition of the charges depends on a choice of grading of the Dynkin diagram; if a different grading were chosen, relations (5.2) and (5.4) would be slightly different. However, we stress that the parameters M A andM I appearing in the asymptotics of the QSC are invariant under these changes, and unambiguously associated to a given multiplet (see [17] for a detailed discussion). Concretely, we may read the charges from the Asymptotic Bethe Ansatz description of the state: where L is the length parameter and K i denotes the number of Bethe roots of type i in the so-called η = +1 version of the ABA [35], while γ is the anomalous dimension. For more details and a dictionary between different forms of the ABA, see Appendix D.
The large-u asymptotics of the matrix Q a|i (u) may be determined by studying (4.2). There are four possible asymptotic behaviours where Q a|i scales as a power of u, parametrized -19 -in terms of the charges M A ,M I entering the equation through (5.1),(5.3). By choosing a suitable linear combination of solutions, we shall impose that different columns of Q a|i have distinct leading asymptotics, ordered in such a way that |Q a|i | > |Q a|j | for i < j for large u.
To describe the possible scaling behaviours, it is convenient to introduce: With these definitions, we have while ν a and ν a have the same leading asymptotic behaviour as Q a|1 , Q a |1 , namely: The asymptotics of Q ij can be computed from the definition (4.8), and turn out to be, for the vector components, where the coefficients B I are constrained by consistency conditions similar to (5.3): , (I = 1, . . . , 5), (5.11) (with no summation on the index I in (5.11)). The trace part satisfies where the constant C coincides with the value of the OSp(4|6) Casimir: A derivation of these constraints is discussed in Appendix C. Finally, let us comment on the asymptotics of the four functions τ i (u). Since the latter are 2i-periodic, and by construction grow less than exponentially for large u, they must approach a vector of constants at infinity.
-20 -There is a certain amount of freedom in normalizing these constants, but we expect that for any physical state the components of τ i with i = 2, 3 always vanish at large u: In Section 4.4 we established (5.15) for the class of 4 ↔4-symmetric operators. While we do not have a fully rigorous argument, we postulate that (5.15) is true in general even for nonsymmetric states. As we discuss in Section 6, the asymptotics (5.15) implies the quantization of the spin and is the main ingredient for deriving the so-called gluing conditions, a powerful set of constraints encoding the main analytic properties of the system. While here we presented these prescriptions on the asymptotics somehow axiomatically, we stress that they are fixed to a great extent by the internal consistency of the QSC equations (see Appendix C). Indeed, one can derive the form of the constraints presented in this Section assuming only that all functions have power-like asymptotics, and that the parameters M A are paired up as in (5.2) The latter condition was fixed by studying the weak coupling limit. The identifications (5.2),(5.4) between the parameters and the charges are strongly motivated by the classical limit discussed in Section 5.1 below, and by several tests both at weak and finite coupling [46,56].

Classical limit
The algebraic curve describing IIA string solutions on AdS 4 ×CP 3 in the classical limit where ∆, S, J i = O(h), h → ∞ 9 was proposed in [33]. In particular, a monodromy matrix was built on the basis of the Lax connection found in [31,32] and its eigenvalues λ a ≡ e iqa were shown to define a ten-sheeted Riemann surface covering the domain of the relevant strong coupling spectral parameter, the Zhukovsky variable x. It is convenient to consider the logarithm of the eigenvalues, the so-called quasi-momenta, naturally grouped as {q 3 , q 4 , q 5 , −q 3 , −q 4 , −q 5 } and {q 1 , q 2 , −q 1 , −q 2 }, corresponding respectively to the SO(6) invariant CP 3 and the Sp(4) invariant AdS 4 sectors of the monodromy matrix. The quasi-momenta are connected by logarithmic cuts 10 , which may be viewed as condensates of Bethe roots. Classical string solutions can be studied by listing algebraic curves satisfying appropriate analytic properties (see [33] for full details), and in particular the charges can be read off the asymptotics of the curve at large values of the spectral parameter: In this limit, h(λ) ∼ λ/2. 10 These cuts exist only in the classical limit and of course they should not be confused with the square-root branch cuts at u = ±2h + iZ considered in the rest of the paper for the QSC.
In the classical limit, we expect that some of the P and Q functions of the QSC are related to the quasi-momenta as follows: Given these identifications, (5.16) justifies the expressions (5.2)-(5.4) for the large-u asymptotics of P and Q functions in terms of the charges 11 . Moreover, at least some of the limits in (5.17)-(5.21), particularly the ones for P 1 , P 2 , Q 1 , Q 2 , can be directly confirmed from the large volume limit of the QSC, see Section 7.3 below.
Another important property of the classical curve is the so-called inversion symmetry, which reads [33]  The symmetry (5.22) is inherited by the transformation property of the monodromy matrix under the Z 4 automorphism of OSp(4|6) [31,32], and is related to the Riemann-Hilbert type equations valid for P and Q at finite coupling. Let us illustrate this point and discuss some consistency checks of the quasi-classical identifications (5.17)-(5.21). Consider the case of P functions. In the classical limit, one can argue (see [17]) that, due to its mirror-periodicity, the matrix µ AB (u) connecting P and P freezes to a constant value independent of u. Moreover, we see from (5.17)-(5.21) that P 1 and P 2 are exponentially suppressed as h → ∞. Therefore, from the QSC equation (3.4) we find On the other hand, starting from the classical expressions (5.17),(5.18) for P 1 and P 2 , analytically continuing to the second sheet and using the inversion symmetry (5.22), one finds (see [17] for details) 24) and the comparison between (5.24) and (5.23) gives the classical identification for P 3 and P 4 . This analysis cannot be straightforwardly repeated for the Q functions, since the τ 's are not constants in the classical limit. However, the inversion symmetry has a quantum analogue in the gluing conditions discussed in Section 6, which connect Q ij and the complex conjugate functions Q ij . From the analytic continuation of (5.17)-(5.21), combined with the inversion symmetry, we may infer that in the classical limit This is indeed consistent with the results of Section 6. As a last comment, notice that there is no classical analogue for two of the components of the matrix Q ij , namely the functions Q 5 and Q • , which enter the basic Riemann-Hilbert constraints at finite coupling, but appear to completely decouple from the dynamics in the classical limit. This is a peculiar feature, as compared with the case of AdS 5 /CF T 4 , and it would be important to find a proper interpretation. One may also speculate that there is a connection with the fact that part of the classical string solutions in ABJM theory are not captured by the classical spectral curve [61].

Unitarity conditions
The structure of the QSC also appears to automatically implement the unitarity bounds satisfied by the charges of a physical state. The discussion here will be very similar to the argument of Section C.2 of [17], so we will only sketch the main points. From the perspective of the QSC, the unitarity bounds arise from the requirement that the powers appearing in the asymptotics of P and Q functions are all distinct. This condition is very natural, since otherwise expressions like (5.3) and (5.11) for the coefficients A A , B I would become singular.
A further condition appears to be needed, namely that, for all consistent solutions of the QSC, the powers entering the asymptotics of Q functions are greater than the ones entering the asymptotics of P functions: precisely, |M A | < |M I |, I = 5. While it is more difficult to motivate this bound from first principles, it can be verified that it holds at weak coupling or in the large volume limit. Assuming a (purely conventional) ordering of magnitude for the components of P A and Q I , we can therefore argue that all non-singular solutions of the QSC can be found restricting our attention tô With the identification (5.2),(5.4), we find that these conditions coincide with the unitarity bounds 27) or equivalently, in terms of excitation numbers (see [53] 12 ): As a final comment, notice that, in principle, some of the inequalities (5.26) could be saturated exactly in the weak coupling limit, where γ → 0. Since the parameters M A , as well asM 2 −M 1 (see Section 6) are quantized, this is possible only for the conditionM 2 > M 2 . The saturation of this bound for γ → 0 is equivalent to the multiplet shortening condition: where ∆ (0) is the classical conformal dimension, or equivalently K 2 = 0 in terms of excitation numbers. The states satisfying (5.30) have a peculiar characteristic in the QSC, namely they are the ones for which one of the P functions vanishes at weak coupling. This is shown by the fact that for these operators

Gluing conditions and spin quantization
We shall now derive an exact relation (valid for real values of the charges) connecting the values of Q ij on the second sheet to the values of the complex conjugate function Q ij . A similar result was first found in the AdS 5 /CF T 4 context and exploited to solve the QSC in various regimes [18,21]. In particular the equations presented below 13 may be used to solve the QSC numerically at finite coupling [56]. For the derivation, we need an important technical assumption: we require that the matrix elements Q a|i can be expanded at large-u as In words, (6.1) means that there is no mixing among the powers occurring in the asymptotics of different columns of Q a|i . This condition was dubbed "pure asymptotics" in [21], and can always be enforced using the freedom to take linear combinations of the columns of Q a|i . We also assume that, for real values of the charges and the coupling, P ab can be chosen to be real 14 . Under these conditions, the conjugate matrix elements Q a|i satisfy the same difference equation (4.2) as Q a|i . This implies that the two matrices are related through where Ω j i (u) is a 2i-periodic function of u: Ω i j (u + 2i) = Ω i j (u). The condition of pure asymptotics (6.1) implies that, as u → ∞, the matrix Ω i j becomes diagonal. Now, we recall the discontinuity relation (4.25): 3) 13 The results presented in this Section were also obtained independently by Riccardo Conti using a slightly different argument [62]. 14 Throughout this section, reality and complex conjugation will be defined on the Riemann section with short cuts. Concretely, the reality of PA means that all coefficients cA,n in (3.2) are real.
-24 - , which, combined with (6.2), gives The crucial observation is now that L i j (u) must be a constant independent of u. In fact, the definition (6.5) can be rewritten as and the last equality shows manifestly that L j i has no cuts in the upper half plane, since this property is true for both Q LHPA a|i and Q a|j . Because of its 2i-periodicity, L i j is then entire in u, and, since it does not grow exponentially, it must be a constant.
To determine the form of L j i , we can study its definition at large u, where Ω i j becomes diagonal and many of the matrix elements of f i j vanish due to the fact that τ 2 , τ 3 → 0. The structure is further specified by several consistency conditions. For instance, since L does not depend on u, we should certainly impose the equality of the following limits: To exploit this constraint, notice that the constant limits of Ω at ±∞ are related as follows: This condition can be obtained studying the definition (6.2) as u → ±∞, using the fact that the asymptotic behaviour of Q a|i (u) (Q a|i (u), respectively) as u → −∞ must be connected to the one for u → +∞ through analytic continuation along a large semicircle in the upper (lower) half plane, where this function is free of singularities. Considering relation (6.6) for j = 2, 3, and using (6.7), we find e 2πi(Na+N i ) = 1, (6.8) for i = 2, 3, ∀a. This equation implies thatM 2 −M 1 = 2S + 1 ∈ Z, namely the spin is integer or half-integer. The other conditions in (6.6) constrain the asymptotics of the nonzero components of τ . Denoting t i,± ≡ lim u→±∞ τ i , we have in particular t 1,± t 4,± = ±i e iP tan(πM 1 ). (6.9) Finally, evaluating L at large u and using (6.9), relation (6.4) leads to the gluing conditions: 10) 11) -25 -where we are using the vector notation defined in Section 4.3.2, δ 1 = e −iP t 2 1,+ /(y 1 y 2 ), δ 2 = −e −iP t 2 4,+ /(y 3 y 4 ), and y i ≡ lim u→+∞ Ω i i (u). For completeness we point out that the constants y i , δ i may in general depend on the coupling and on various normalization choices. For the implementation of the numerical method, it is only needed to know explicitly the value of y i . These constants, which satisfy the consistency conditions y 1 = 1/y 4 = 1/(y * 1 ), y 2 = 1/y 3 = 1/(y * 2 ), are simply related 15 to the choice of normalization of the Q a|i (u) functions, and can be determined as: The relations (6.10)-(6.12) are similar to the ones obtained in [18,21], but slightly more complicated. Indeed, in the AdS 5 /CF T 4 context a single Q function appears on the rhs of the gluing conditions, which are an almost direct lift of the inversion symmetry connecting pairs of quasi-momenta in the classical limit. In the present case, the quantum version is a bit more intricate. In particular, the explicit parametric dependence of the gluing conditions on the chargeM 1 needs to be taken into account in order to develop a numerical algorithm [56]. As a last comment, we observe that the quantization of the spin is a direct consequence of the choice of vanishing asymptotics for two of the components of τ . As shown in [18], it should be possible to relax this condition and consider continuous values of S by admitting exponentially growing asymptotics in τ 2 and τ 3 .

The Q-system
In this Section we show how to embed the previous results into a larger set of functional equations reflecting the OSp(4|6) symmetry. It is important to mention that, while the form of Q-systems associated to GL(M |N )-type superalgebras is known (see e.g. [14,25,63]), there appears to be no comprehensive understanding of this mathematical structure for orthosymplectic superalgebras. Here we take a bottom-up approach to the problem and try to construct the Q-system starting from the Q functions already introduced 16 : P A , Q I , Q a|i , Q a |i , together with the relations linking them, equations (4.4),(4.5),(4.8). We will explicitly define new Q functions and prove the validity of a set of functional relations which is rich enough to contain various forms of exact Bethe Ansatz equations (equivalent to the absence of poles for the Q functions) related to the OSp(4|6) symmetry.
Before starting the construction, let us describe some of its main characteristics. Various types of Q functions will be assigned to particular nodes of the Dynkin diagram. We will almost exclusively consider the two versions of the diagram shown in Figures 5, 6, which are 15 Actually, for real values of the coupling it is always possible to choose a normalization where B (a|i),0 ∈ R, so that yi = 1. 16 Starting from these functions, we will define a Q-system where the Q functions are free of cuts in the upper half plane. An analogous construction, analytic in the lower half plane, could be performed starting from the Q functions PA, QI , and Q LHPA a|i defined in (4.27). Notice that the two systems are connected through the ν or τ functions, which therefore play the role of a symmetry transformation of the Q-system (for an interesting discussion see [17]).  the ones associated to the two known forms of Asymptotic Bethe Ansatz. The Q functions will have the general index structure 17 Q •| * , where • and * are (vector or spinor) multi-indices carrying representations of SO(3, 3) and SO(3, 2), respectively, see Section 2 for notations. Various arguments, and in particular the weak coupling analysis, suggest that Q functions of types P A and Q I carry Bethe roots associated to the first node of the two diagrams, while the Q functions Q a|i , Q a |i should be linked to the nodes corresponding to the spinorial representations, see Figures 5, 6. The main task of this Section is to complete the picture by constructing Q functions and functional relations associated to the remaining nodes. In analogy with the Q-system of [17], and in contrast to the case of standard Lie algebras, for every node of the diagram one may define equations of two basic types -fermionic or bosonic. This feature of supersymmetric Q-systems is known to be related to the existence of different gradings of the Dynkin diagram. Choosing different chains of Q functions, we will recover different sets of exact Bethe equations. Finally, as a non-trivial check of the construction, we will recover the two forms of the ABA equations in the large volume limit.

Construction of the Q-system
First step: identifying Q A|I We start the construction by some guesswork. From the form of the Bethe Ansatz, and taking inspiration from [17], it is natural to expect that one of the functional relations should read: We have marked this equation with the symbol F 1 to point out that it is a fermionic-type Q-system relation, based at the first node of the Dynkin diagram. This equation might be 17 Notice that also the Q functions PA and QI fit this pattern and we could identify them with PA ≡ Q A|∅ , QI ≡ Q ∅|I , where ∅ denotes the trivial representation.
-27 -taken as a non-local definition of the 6 × 5 matrix 18 Q A|I . However, this new type of Q functions can also be expressed as an explicit, local combination of the building blocks Q a|i , Q a |i , through the following quadratic combinations: namely, the 2 × 2 minors of the 4 × 4 matrix Q a|i . Notice that Q ab|ij is antisymmetric in both (ab) and (ij), and therefore has 6 × 6 independent components. To match the 6 × 5 components of Q A|I we need of course to project the (ij) indices on the vector representation. The correct identification, which will be important for the derivation of the rest of the Qsystem, is simply: We will show below that this definition implies the validity of (7.1). One could also consider the complementary projection on the singlet representation for the (ij) indices, and define: However, it turns out that all Q functions carrying the singlet representation of SO(3, 2), such as Q A|• and Q • , drop out of the functional relations needed for the derivation of exact Bethe equations. It would be interesting to understand from the algebraic point of view whether they should be considered as part of the Q-system. To prove the validity of (7.1), we start by rewriting the constraint Pf(P ab ) = 1 as: P ab P cd − P cb P ad − P ac P bd = abcd , (7.5) where abcd denotes the completely antisymmetric Levi-Civita tensor. Using this identity, it is immediate to prove that 19 and, inserting (4.14), we obtain 18 Notice that we are denoting Q functions carrying capital indices such as A ∈ {1, . . . , 6} or I ∈ {1, . . . , 5} with the calligraphic font Q in order to avoid possible confusion with Q a|i when the indices take some concrete value. So, for example, notice that Q 1|2 = Q 1|2 ! 19 We are using the standard notation [ , ] for the antisymmetrization of indices, e.g.
Projecting on vector indices as in (7.3) and taking into account simple algebraic identities (see (B.24)), (7.7) yields precisely the fermionic equation (7.1): For completeness, we report also the identity obtained by tracing over (ij): As anticipated, (7.9) is apparently decoupled from the rest of the Q-system and will not play a role in the following considerations. Bosonic-type Q-system relations for the first node can be introduced straightforwardly. They take the standard form: 10) which can be interpreted as definitions of the new two-index objects Q AB|∅ and Q ∅|IJ . These Q functions do not sit on the diagrams in Figures 5, 6, but appear in other choices of gradings, such as the distinguished one (see discussion below). The construction of functional relations for the second and third nodes is standard and follows the usual fusion rules, cf [17]. In particular, associated to the third node we define the Q functions which satisfy bosonic-type relations for the second node: 14) Using equation F 1 (7.8), we can also straightforwardly establish the following fermionic-type functional relations for the second node: Now let us derive the relations centered around the third node. Using (7.12)-(7.13), it is simple to obtain the bosonic-type equations -29 -while the definitions (7.12),(7.13) and relation (7.8), imply the validity of the fermionic identity As we may expect from the Dynkin diagram, the newly defined object in (7.21) represents the fusion of the spinorial Q functions Q a|i and Q a |i . Indeed, let us prove that it can be rewritten as: where Σ ij IJ ≡ 1 2 (Σ I κΣ J −Σ J κΣ I ) ij . This equation will be crucial for the derivation of closed sets of exact Bethe equations. To derive (7.22), start from the definition of Q A|I in (7.3) and rewrite (7.21) as Using formula (B.8) for the commutator of sigma matrices appearing in (7.23), we find where, in the last step, we have used the anti-symmetry in (IJ) of the whole expression by definition of Q AB|IJ .

Q-system relations for the nodes 4 and4
Let us now derive the functional relations centered at the spinor nodes. The two bosonic Q-system equations (centered at nodes 4 and4, respectively) are: while the fermionic-type relations, which cross the two spinor nodes, read To prove (7.25), start from the combination -30 -Using (4.15), (4.14), (4.11), we can eliminate all positive shifts through and we find 20 : where we have used identity (B.11) to simplify the product of Σ matrices in (7.31). Contracting with (σ A ) ab and comparing with (7.12) yields (7.25). Similarly, to prove (7.27), we consider and replace all Q functions with positive shifts using (Q a |i ) + = P aa 1 Q − a 1 |i :

Exact Bethe equations
Let us now show how to obtain exact Bethe equations for the zeros of Q functions. We will obtain equations formally identical to the various versions of 2-loop Bethe Ansatz proposed in [29], based on the underlying OSp(4|6) symmetry, with the important difference that, at finite coupling, Q functions are nontrivial functions of the spectral parameter living on infinitely many sheets (and, in general, with infinitely many zeros). In the weak coupling limit, the branch cuts shrink to zero size, usually being replaced by poles. However, some of the Q functions become polynomials, and in particular the equations presented below reduce to the weak coupling BA of [29].
To derive a version of the Bethe Ansatz related to the η = 1 grading of the Dynkin diagram, we need to consider a chain of functional relations made of equations of type F 1 (7.8), B 2 (7.14) and F 3 (7.20) for the first, second and third nodes respectively, and B 4 (7.25) 20 Notice that the terms proportional to Q• cancel out of the equation due to the symmetry (ΣIJ ) ij = (ΣIJ ) ji , see Appendix B.
-31 -and B4 (7.26) for the nodes at the bifurcation. For concreteness, let us make a specific choice of indices, and consider the following sequence of Q-system relations where we used (7.22) to evaluate Taking the ratio of (7.38) evaluated at points u 4,k + i/2 and u 4,k − i/2, where u 4,k is a generic zero of Q 1|1 , gives the massive node Bethe equation with Q 1|1 (u 4,k ) = 0, (7.42) and similarly from (7.39) one gets Auxiliary equations for the fermionic nodes are obtained simply by evaluating (7.35) and (7.37) at the respective zeros u 1,k and u 3,k of their rhs: while the Bethe equation for the second node is obtained by taking the ratio of (7.49) computed at u 2,k + i/2 and u 2,k − i/2: with Q 2|2 (u 2,k ) = 0. (7.46) In Section 7.3, we will show that in the large volume limit these equations reduce to the η = 1 form of the ABA [35]. We can describe an alternative grading by using relation B 2 * (7.15) instead of B 2 for the second node and the fermionic-type equations (7.27),(7.28) for the nodes 4 and4. Consider for example the chain of Q functions connected by the Q-system relations Using the pole-free condition, they straightforwardly lead to exact BA equations corresponding to the Dynkin diagram of Figure 6: The main difference with respect to the derivation in the η = +1 case concerns the equations for the momentum-carrying nodes: for instance, (7.53) is obtained by taking the ratio of equation (7.51) evaluated at u 4,k + i/2 and equation (7.52) at u 4,k − i/2. As shown in the next Section 7.3, equations (7.53)-(7.57) reduce to the η = −1 version of the ABA of [35] in the large-L limit. We may also consider subsets of Q functions whose zeros satisfy exact Bethe equations related to the so-called "distinguished" grading of the Dynkin diagram. An example of such a chain is: The Bethe equations associated to the momentum-carrying nodes are (7.42), (7.43). To constrain the remaining Q functions, we may use B 1 * (7.11), F 2 * (7.17) and B 3 * (7.19) with indices A, I = 1; B, J = 2. Employing standard arguments, we find the Bethe equations: At the leading weak coupling order these equations reduce to one of the variants of the 2loop Bethe Ansatz of [29]. However, it is well known that this grading is impractical when considering the large-volume limit and does not lead to simple Asymptotic Bethe equations.

The ABA limit
Let us now argue that in the large volume limit a subset of Q functions -in particular, the ones appearing in the chains (7.41) and (7.47) -reduces to a simple explicit form parametrized by a finite set of Bethe roots living on two sheets only. The exact BA equations (7.42)-(7.46) and (7.53)-(7.57) will then be shown to reproduce the Asymptotic Bethe Ansatz of [35]. The following argument is very similar to the one presented in [17]. The main origin of the simplification occurring in the large volume limit is that some of the Q functions vanish at an exponential rate at large L. To keep track of the scaling of different quantities with L, we can rely heuristically on the asymptotics discussed in Section 5. From (5.5),(5.6), we see that the charges scale as ∆, J 1 , J 2 ∼ L, while S, J 3 ∼ O(1) at large L, from which we get for example that ν a ∼ (1, 1/ε, 1/ε, 1/ε 2 ), ν a ∼ (1/ε 2 , 1/ε, 1/ε, 1), (7.62) where ε ∼ u −L represents a quantity exponentially suppressed in L. Similarly, we have Moreover, since the functions τ i approach constants with no explicit L dependence at large u, we deduce that they scale as O(1) in the large volume limit. Using this information, we -34 -obtain some simplified relations. Let us list the ones most relevant for the derivation of the ABA. First, from the scaling (7.63) we find that (4.31) reduces to: Second, from (3.12) we find, for α = 1, 2, where we used also the identity (7.25) in the last step, and we recall that ω 12 = τ 1 τ 4 . Finally, it will be useful to consider the relation between the Q functions analytic in the upper/lower half plane, which simplifies in the large volume limit. In particular, we have from which we see that equation (7.67) can be rewritten as Computing µ 12 , ω 12 and Q 12|12 The first part of the argument is essentialy the same as in [17]. We shall assume that ν 1 and ν 4 have each a finite number of real zeros on the first sheet in physical kinematics, which we denote as {u 4,j } K 4 j=1 , u4 ,j K4 j=1 respectively. We start by defining where we remind the reader that µ 12 = ν 1 ν 4 and F is manifestly free of poles on the first sheet. Using (7.66), we can rewrite this quantity as where the contribution of ω 12 = τ 1 τ 4 cancels due to its i-periodicity and we used (7.40) in the last equality. The expression (7.72) shows that, within this approximation, F 2 is built out of quantities that have manifestly no cuts in the upper half plane. On the other hand, using (7.68) we see that F 2 could equivalently be rewritten in terms of LHPA Q functions -35 -only. We therefore conclude that it must have no branch cuts apart from a short cut running on the real axis. The discontinuity across the latter is described by the condition which is a simple consequence of (7.70). These analyticity requirements completely fix F (but for a sign) as: and similarly we find Already at this stage, we can prove that the zero momentum condition (3.15) is contained in the QSC equations. Indeed, from (7.78) we have: in the ABA limit. Due to the mirror i-periodicity of µ 12 , this ratio should approach 1 at large u. Expanding the rhs of (7.80), and taking into account the dispersion relation The next order in the large-u expansion can be compared with the asymptotics (5.7)-(5.8), and fixes the ABA limit of the anomalous dimension: Let us now show that the ratio between Q 1|1 and Q 4 |1 must be, in the large-L limit, a meromorphic function without branch cuts (and therefore, due to the power-like asymptotics, it must be a rational function of u). Indeed, equation (7.68) shows that The analyticity strips of the two sides of (7.83) overlap nontrivially, showing that this ratio is indeed a ratio of polynomials. The correct way to split (7.79) is then The solution of the constraints (7.86) may be found in terms of an integral representation 21 : 2i .
(7.87) We should also impose that log F(u) has the correct bounded asymptotic behaviour as u → +∞, which leads to the condition As already discussed in Section 3.3, the expression for P in (7.88) is expected to hold only in the large-L limit, or at the first ∼ L orders at weak coupling. A general exact integral formula for P, expressed in terms of quantities computable form the numerical solution of the QSC, can be found in Appendix E.
Computing P α , Q α|12 and Q α|β Let us now derive the ABA limit of P α , with α = 1, 2 (again, we follow [17] closely). We define where σ has a single short cut on the real axis on its defining sheet. Since (7.90) is simply one of the crossing equations, it follows that σ is related to the dressing factor as in (D.14). Let us consider the quantity P α /σ, which by construction has a single cut on the first sheet. Using (7.67),(7.79), we see that, on the second sheet, it may be written as s ), (7.91) which has no cuts in the upper half plane, or alternatively from (7.69) as ), (7.92) which has no cuts in the lower half plane. Hence, P α / σ must have a single cut on the second sheet as well, so that it may be written as a rational function in the Zhukovsky variable x(u). Therefore, we have P α ∝ x −L B α|12 R α|∅ σ, α = 1, 2, (7.93) 22 Notice that the ABA expression for the total momentum of a single excitation species is given by: which agrees with the rhs of (7.88) at the first two orders at weak coupling.
-38 - where the x −L prefactor is fixed by imposing the large-u asymptotics (5.1), and we have introduced the notation B α|12 (u) (R α|∅ (u)) to indicate generic polynomials in x(u) (1/x(u), respectively), see Appendix D for a precise definition. By consistency with (7.67), we then find: where R α|12 (u) = B α|12 (u) and B α|∅ (u) = R α|∅ (u) are obtained through analytic continuation, which sends x(u) → 1/x(u). At this stage, we have computed four of the functions entering the chain (7.41); to complete the picture we still need to compute the Q functions corresponding to the second node. We start from relation which is a consequence of (7.68), and implies that ratios of the form 23 have no cuts and are therefore ratios of polynomials. We have therefore a parametrization where Q α|β is a polynomial function of u, and the f 4 f4 factor was fixed by comparison with (7.79).
Asymptotic Bethe Ansatz in η = +1 grading Generalizing the arguments of Section 7.2, we see that the Q functions for any choice of α, β ∈ {1, 2}, satisfy exact Bethe equations of the form (7.42)-(7.45). Using (7.84), (7.93), (7.94), (7.97), it is straightforward to verify that, in the large volume limit, these Bethe equations reduce precisely to the ABA of [35] in η = +1 grading (see Appendix D). In each of these four equivalent sets of ABA equations, the role of roots of types 1,2,3, is played by the zeros of the following polynomials in u: Computing Q 1 and Q 2 The large volume limit of Q β with β = 1, 2, may be computed from the Q-system relation F 1 , namely: 99) 23 Notice the restriction of the indices to the set {1, 2}. This ensures that the ratios in (7.96) are of order O(1) for large L, which is a prerequisite condition for obtaining nontrivial information in the asymptotic limit.
-39 -for α, β ∈ {1, 2}. Similarly, Q 12|β may be determined from the F 3 equation: Using the large-L expressions (7.93), (7.94) and (7.97), these relations yield where the functions R ∅|α and R 12|α (B ∅|α and B 12|α , respectively) are polynomials in x(u) (1/x(u)) defined by Notice that the fact that the newly defined R and B functions have no poles is a consequence of the ABA. Equations (7.103)-(7.103) are the well-known fermionic duality relations, which allow to switch between the η = ±1 versions of the ABA, see Section D.2. Using (7.84), (7.97), (7.93), (7.101), we may indeed check that the exact Bethe Ansatz satisfied by the chains of Q functions which in particular involves the fermionic form of the massive node equations, (7.53),(7.54), reduce precisely to the η = −1 ABA equations. As a last comment, we point out that, from the large volume limit, we can get a further confirmation of the semi-classical identifications (5.17)-(5.21). To this end, we exploit the well-known fact that the classical spectral curve can be obtained as a scaling limit of the Asymptotic Bethe Ansatz, where the roots condense to form the cuts connecting various pairs of quasi-momenta [64]. Consider for instance, the following large volume expression derived from (7.101): In the limit where h ∼ L 1, it is meaningful to concentrate on the region u > h, where the lhs of (7.105) becomes approximately exp(∂ u log Q 2 ). On the other hand, a standard calculation shows that, in the limit where the Bethe roots scale like u s,j ∼ h and condense to form a set of cuts, the rhs reconstructs precisely exp(iq 2 ), where the quasi-momentum q 2 is defined in terms of Bethe root densities as in [35,38] (see also Section 6 of [17]). Therefore we recover the identification (5.21). Similarly one could derive the classical limits of the remaining P and Q functions which we have determined in the large volume limit.
-40 -In this paper, besides a detailed derivation of the equations proposed in [44], we presented several new results on the Quantum Spectral Curve associated to the AdS 4 /CF T 3 duality, deepening our understanding of the basic integrable structures underlying this theory.
There are many directions for future work. First of all, the results of this paper make it possible to develop a high-precision numerical algorithm for the computation of anomalous dimensions at finite coupling, inspired by [18]. We already have partial results [56,65] confirming the TBA data of [39]. The QSC method however allows us to move deeper in the strong coupling region, and therefore to test more accurately the AdS/CFT predictions.
Secondly, we expect from the example of AdS 5 /CF T 4 [26,66,67] that the QSC may be used, with minimal modifications, to describe also various open string configurations. In particular, it would be very interesting to find an integrable description of some kind of generalized cusp anomalous dimension, such as the one described in [68]. This would give a direct way to test the proposals of [44,45] for the ABJM/ABJ interpolating functions, by comparison with localization results for the Brehmsstrahlung function [69][70][71][72].
Third, these results should allow to extend the weak coupling algorithm of [46] to a generic operator.
It would be very interesting to gain a complete understanding of the algebraic structures underlying our results. Especially, it would be desirable to understand the interpretation of the Q-system described in Section 7 in terms of representation theory of the full supergroup OSp(4|6).
We hope that the results presented in this paper, which exhibit some interesting differences from the AdS 5 /CF T 4 case, will also help to extend the QSC method to the integrable examples of AdS 3 /CF T 2 and AdS 2 /CF T 1 , see e.g. [47,[73][74][75]. These cases are less supersymmetric, and the construction may be expected to be even more complicated. It is important to stress that, since a TBA formulation for these models is at present still missing (and even the structure of the Asymptotic Bethe Ansatz is quite intricate and fully known only in one case, see [76]), there is presently no way to rigorously derive the QSC for these theories. However, the two examples at hand, AdS 5 /CF T 4 and AdS 4 /CF T 3 , show that the structure of the QSC is, in the end, quite universal and rigidly constrained by the symmetry. It would be very nice if these examples could help to develop a classification of several types of QSC corresponding to different gauge and string theories.
In particular, we thank Riccardo Conti for collaboration on the project [56], during which many aspects of the present work were elucidated.
This project was partially supported by the INFN (I.S. FTECP and GAST), UniTo-SanPaolo Nr TO-Call3-2012-0088 "Modern Applications of String Theory" (MAST), ESF Network HoloGrav (09-RNP-092 (PESC)), MPNS-COST Action MP1210 "The String Theory Universe", and the EU network GATIS. AC thanks King's College London for kind hospitality during two visits in 2016, where part of this work was done.

A Derivation of the QSC from the analytic properties of T functions
In this Appendix we present in detail the derivation of the QSC equations from the TBA/Tsystem framework, which was already outlined in [44]. In particular, we will obtain the QSC equations in the "Pµ" vector form presented in Section 3.1.

A.1 Summary on the properties of T functions
Let us briefly summarize the starting point of the derivation (see [43] for more details). The discrete Hirota equation, or T-system, is the following difference equation for a set of T functions defined on the nodes of the "T-hook" diagram shown in Figure 7: a,−1 = T α a+1,0 T β a−1,0 + T a,1 T β a,−1 , α, β ∈ {I, II} , α = β, where T functions with indices outside the diagram are taken to be zero and the products in (A.1) are over horizontal (↔) and vertical ( ) neighbouring nodes, with the subtlety that, for s = 0, −1, the two wings of the diagrams need to be crossed 24 . Notice that T (n) = T (u + i 2 n) denotes shifts on a specific section of the u domain where all cuts are long, connecting ±2h+iZ to infinity. This is called the mirror section and is the one where the Y-system and Tsystem are naturally defined [12]. Throughout this Appendix we will use the special notation f (n) (u) ≡ f (u + in/2) to denote a function shifted on this particular sheet. T functions are related to Y functions, the objects appearing in the TBA formulation, by This parametrization is not unique: there is a vast "gauge" freedom (which we will exploit) in choosing a set of T functions corresponding to a given solution of the TBA. In order to furnish a complete formulation of the spectral problem, the T-system must be supplemented by some information on its analytic dependence on the spectral parameter. As learnt in the AdS 5 /CF T 4 case, this extra input can be expressed in terms of discontinuity relations for the Y (u) functions across their branch cuts in the u-plane [12], but can be simplified and much better understood in the T-system framework [15]. In the case of AdS 4 /CF T 3 , similar analytic constraints on the T functions were identified in [43]. They are expressed in terms of two special gauges, denoted as T and T. The properties of the T gauge needed in the following derivation are: (i) Analyticity strips: denoting as A n the class of functions free of branch point singularities in the strip |Im(u)| < n 2 , we have (T α n,0 ) ∈ A n+1 , (T n,1 ) ∈ A n , (T n,2 ) ∈ A n−1 , n ∈ N, α ∈ {I, II} .
(ii) The two functions T I 0,0 , T II 0,0 are equal, and periodic on the mirror section: The functionμ 12 defined above will eventually be identified with an element of the µ AB matrix appearing in the QSC equations. The notationμ 12 signals that, throughout this Appendix, we will considerμ 12 (u) as a function defined on the mirror Riemann section with long cuts, where it is i-periodic. This function agrees with µ 12 (u) used in the rest of the paper in the strip 0 < Im(u) < 1/2, and elsewhere is obtained by analytic continuation keeping all cuts long. Notice that the mirror i-periodicity ofμ(u) is equivalent to the property (3.3).
(iii) Finally, the T functions enjoy the following group-theoretical properties: We expect that the T gauge defined by these properties is essentially unique (apart from rescalings by constants independent of u). The T gauge may be defined by a transformation: T n,s (u) = T n,s (u) μ where s I = −s II = +1, and d(u) = d (+2) (u) is a mirror i-periodic function, representing an additional degree of freedom in the definition which will be practically irrelevant for our derivation 25 . It is simple to check that (A.9)-(A.11) leave invariant the form of the T-system due to the mirror periodicity ofμ 12 and d(u).
In general, we expect both the T and T functions to exhibit an infinite ladder of branch points for u ∈ ±2h + iZ/2. From the TBA analysis, we know that these singularities are all of square-root type and that analytic continuation around branch points symmetric with respect to the imaginary axis leads to the same sheets. This structure is further specified by the property (i) above: some of the potential branch points in the T functions fall inside the analyticity strips and therefore they must have trivial monodromy.
Besides, the T functions enjoy some special properties when continued to the short-cut section of the Riemann surface (also known as the physical sheet). We will denote their values on this section asT: in analogy with the case ofμ and µ, the convention is that T andT are the same in the analyticity strip immediately above the real axis, while in the rest of the complex plane, they are defined by analytic continuation keeping long cuts for T and short cuts forT. TheT a,s functions have the following nontrivial properties: The goal of the following derivation is to obtain the Riemann-Hilbert type equations characterizing the QSC. We will see that the whole structure can be derived by imposing the consistency of the conditions (i), (ii) , (iii) and (a), (b). Let us make an additional comment. Here, we do not aim to derive the regularity properties of the QSC, namely the statement that P(u) and ν(u) functions are entire on the Riemann surface defined by the branch points at u ∈ ±2h ± iZ. However, it is natural to expect that this condition is equivalent to the requirement that the T functions are regular in appropriate gauges, and indeed one can verify a posteriori that, picking appropriately the function d(u) in (A.9) and assuming the regularity of the QSC, all the T and T functions can be chosen to be regular. For instance, it is possible to identify 26 T I 1,0 (u) = ν 1 (u) ν 1 (u), T II 1,0 (u) = ν 4 (u) ν 4 (u). Therefore, choosing d(u) ∝ (ν 1 (u)/ν 4 (u)) 1 2 , one can set T I 1,0 (u) ∝ ν 1 (u), T II 1,0 (u) ∝ ν 4 (u), from which we have a clear indication that the regularity properties of the ν and T functions are equivalent. This example also illustrates the fact that a requirement of regularity for the T gauge specifies the function d(u) uniquely, apart for an overall constant 27 . However, we remark that, for the purposes of the following derivation, the precise form d(u) is irrelevant: this function cancels out of all the equations reported below.

A.2 Strategy of the derivation
The main tactic of the derivation is to choose a parametrization of the T functions that makes (a), (b) explicit; we will then reconstruct the T functions through (A.9) and impose the validity of (i), (ii), (iii).
To start, we notice that the properties (a), (b) presented above can be encapsulated by the following parametrization [15]: where P 1 , P 2 are functions with a single short cut. Notice that this parametrization covers only the right tail of the T-hook diagram. To reach the rest of the diagram using the T-system relation (A.1), we need one more constraint involving at least one node outside this domain.
For this purpose we may use T 3,2 /T 2,3 =μ 12 , (A.14) which follows from the transformation (A.9) combined with the property (iii). We then see that, applying Hirota equation starting from any point in the right band, we may parametrize any of the T functions in terms of only three building blocks, the functions P 1 , P 2 , µ 12 , which as we will see will be evaluated on various Riemann sheets. The T functions, defined through (A.9), can be expressed in terms of the same data, and one can check that they satisfy the constraints (ii), (iii) by construction. However, it is not obvious that they have the correct analyticity strips described by condition (i); we still need to impose an infinite ladder of relations: where we use the symbol for the discontinuity f ≡ f −f expressing the monodromy around any of the branch points at ±2h on the real axis. The conditions (A.15) place further constraints on P 1 , P 2 andμ 12 and will lead us to the QSC equations. 26 These expressions for T α 1,0 follow from the comparison between equation (A.28) below and the Pν-system. 27 In fact, from (A.9) it is evident that this function must be chosen in such a way that it cancels the extra singularities in T α a,0 introduced by the square root factors √μ 12 in (A.9). For states with 4 ↔4-symmetry, we can simply set d(u) = 1, since in that case µ12 has only double zeros.
-45 -As a convenient notation, we will introduce a splitting function g(u), defined through 16) In particular, in the 4 ↔4-symmetric subsector in which T I n,0 = T II n,0 , one has simply g(u) = 1.

A.3 Details
Before discussing the derivation in detail, let us mention a technical point. In the following paragraphs, we will find relations between functions which are defined, by default, on different sections of the Riemann surface covering the u plane. To remove possible ambiguities, we specify that all the equations below are valid for u in a strip slightly above the real axis. With this understanding, we will use interchangeablyμ 12 and µ 12 in the following equations.
First level n = 0 The first constraint coming from (A.15) is that T 2,1 = T 2,1 has no cut on the real axis. The consequences of this requirement were already discussed in [44]. Using Hirota equation and carefully continuing the expressions (A.12) to the mirror sheet, we find Imposing the absence of a cut on the real axis, we obtain (T 2,1 ) = T 1,2 µ 12 − µ 12 − P 1 P 2 + P 2 P 1 = 0, (A. 19) and, since T 1,2 cannot be zero everywhere, we get a first relation of the Pµ-system (3.4): Using the Hirota equation centered at the node (1, 1), we can also compute which means we can parametrize The requirement that T α 1,0 have no cuts on the real axis then imposes (g(u)) = 0.

General level
As illustrated in the previous example, the functions T a,s with a > s, computed using the T-system relations, will depend not only on the values of P 1 , P 2 and µ 12 on their defining sheet, but also on their shifted values on the second sheet: this is due to the fact that Hirota equation is defined on the mirror section, while the cut in the definition of P 1 and P 2 is short. In general, the constraints (A.15) can be translated as conditions on the monodromies of P 1 , P 2 and g around the branch points lying further and further from the real axis. Remarkably, the content of (A.15) can be recast in a very simple form: the constraints on the cuts in the upper half plane yield 28 together with the condition that One obtains very similar but not identical equations describing the discontinuities in the lower half plane. For conciseness, we will only refer to (A.23) in the following arguments.
Remarkably, the form of these relations contains already the full structure of the QSC.

Constructing the Pµ-system
The equations in (A.23) can be rewritten as for A ∈ {1, 2}. Considering the discontinuity of these relations on the real axis, we see that (η n ) = 0. Inspecting expression (A.24), we then see that − P (+2n) 1 with g I = g, g II = 1/g. We will exploit (A.27) to construct two new functions with a single short cut, which we denote as P 5 and P 6 in anticipation of their role in the QSC equations. They are defined through where the functions φ A,α , with indices A ∈ {1, 2}, α ∈ {I, II}, are defined from the requirement that they are periodic on the mirror section, with power-like asymptotics and with discontinuities 29 ) = 0, ∀n ∈ N + . A simple extension of this analysis shows that (A.28) defines a function with only a single short cut on the real axis.
Let us point out that, when g = 1 (which is appropriate for 4 ↔4-symmetric states), by definition we have P 5 = P 6 , in agreement with the rules described in Section 4.4. As another side remark, notice that the definitions (A.28) can be recognized as two equations of the Pν-system (3.12) provided the mirror-periodic functions φ A,α are identified as ratios of ν functions, and g(u) is identified as In the rest of this Appendix, for simplicity we will concentrate solely on obtaining the "Pµ" vector form of the equations. Using (A.28),(A.24), equations (A.26) becomes AB , for A, B ∈ {1, 2}, whose (periodically repeated) discontinuities are Then, defining the functions P 3 and P 4 as −P 3 ≡ P 1 µ 12 + Φ 12 P 1 − Φ 11 P 2 + φ 1,II P 5 + φ 1,I P 6 , (A.33) −P 4 ≡ P 2 µ 12 + Φ 22 P 1 − Φ 21 P 2 + φ 2,II P 5 + φ 2,I P 6 , (A.34) 29 These requirements specify ωA,α uniquely apart from an additive constant independent of u. therefore P 3 and P 4 are free of branch points in the upper half plane and, by a small additional effort, we can show that they have just a single short cut on the real axis. Let us summarize the situation: by a scrutiny of the equations, we have so far found six functions with a single short cut, and eight mirror-periodic functions Φ AB , φ A,α . It remains only to check that the relations between their monodromies can be written in a closed form.
The fifteen components of the antisymmetric matrix µ AB can be defined in terms of the periodic functions introduced above. Indeed, setting we immediately recognize that (A.33),(A.34) are two equations of the Pµ-system. Besides, the form of these relations implies the existence of three quadratic constraints among the matrix elements defined in (A.36); let us discuss in detail how these conditions emerge. Consider the following equation: which can be derived from (A.33) and its analytic continuation to the second sheet using the monodromy rules (A.20),(A.29),(A.32). The form of (A.38) implies that the combination of mirror-periodic functions Φ 11 − φ 1,I φ 1,II must be free of cuts. Due to its power-like asymptotics, it must be a constant independent of u and, using the freedom to redefine the Φ's by a constant shift, we will assume that Φ 11 − φ 1,I φ 1,II = 0. Therefore, (A.38) can be recognized as another equation of the Pµ-system. Moreover, the quadratic constraint we have just found can be rewritten as µ 12 µ 13 − µ 15 µ 16 = 0, which is one of the components of the matrix equation (µη) 2 = 0. By similar reasoning, we can impose two more constraints and all in all we can set The rest of the derivation goes along the same lines. The remaining independent entries of µ AB are defined as: and it is possible to verify that all equations of the Pµ system, including the quadratic constraints on the P and µ functions, follow from the relations listed above (and their analytic continuation around the branch cut on the real axis).
-49 - The specific form of the matrix η AB entering the Pµ-system equations does depend on the normalization of our definitions (A.28),(A.33),(A.34),(A.36)-(A.42), and could be changed by rescaling some of the µ or P functions, or by a more general linear change of basis, However, η AB is clearly always a symmetric tensor, and besides its signature (+ + + − −−) is invariant under all linear transformations with H ∈ R 6×6 . This reality restriction is meaningful since it preserves the following property: for real values of the coupling, all the functions P A (u) can be chosen to be real 30 on the Riemann section with short cuts. This property is verified with our choice of conventions, and follows from the reality of the solutions of the TBA.

B.1 Identities for gamma matrices
In this Appendix we collect some useful algebraic identities, descending from the properties of gamma and sigma matrices for SO(3, 3) and SO (3,2). The defining relation for the SO(3, 3) sigma matrices is and we recall that (σ AB ) b a is defined through A useful property, specific to orthogonal groups in six and five dimensions, is the fact that gamma matrices are anti-symmetric: (σ A ) ab = −(σ A ) ba . This allows us to prove the following very useful relation: and its consequence Tr Another identity that is specific to this dimension is which implies in particular that (σ AB ) is traceless: (σ AB ) a a = 0, and moreover that, for any anti-symmetric matrix 4 × 4 matrix G ab : where the corresponding vector {G A } 6 A=1 is defined by G ab = G A (σ A ) ab . Another useful formula is: All the properties listed above are independent on for any choice of chiral representation of the gamma matrices. The situation is analogous for the representations of SO (3,2). In that case we recall that we use the symbols (Σ I ) ij , (Σ IJ ) j i , and denote the metric as ρ IJ ≡ 1 2 Tr(Σ IΣJ ). In particular the defining relation for the matrices Σ I and Σ IJ is: with Σ IJ = −Σ JI . On top of these properties, in the SO(3, 2) case the matrices Σ andΣ are related by a similarity transformation: where κ ij is an anti-symmetric 4 × 4 matrix. Equation (B.10) can be used to prove the additional symmetry property (Σ IJ ) ij = +(Σ IJ ) ji . Finally, the analogue of (B.5), (B.4) are Finally, we report below some useful identities for a generic antisymmetric 4 × 4 matrix G ab = −G ba : where we recall that the Pfaffian is defined as In particular: However, a simple linear algebra identity relates the minors of a matrix and its inverse, and shows that the two definitions are algebraically related: From (B.20), we see that and using (B.18) we find Let us define the projections: where Q ab|(ij) denotes the traceless part and satisfies Q ab|(ij) κ ij = 0. Identity (B.22) then splits as

B.3 Relation between Q ij and its inverse
From (B.16), we have and, using (B.13), we immediately find Contracting (B.26) with κ ij , we find that in fact Q • = Q • = Q ij κ ij , so that (B.26) reduces to equation (4.14) presented in the main text.

C Derivation of constraints on large-u asymptotics
Here we derive the constraints (5.3), (5.11) on the asymptotics of P and Q functions using the QQ-relations derived in Section 7. In order to find (5.3), we start from relation (7.8). At large u, its rhs is given by which constrains the asymptotic behaviour of Q A|I to be We may now use the following relation, which is a consequence of the Q-system: and gives, using the aymptotics (C.2), the constraint These relations, together with the constraint Pf (P ij ) = 1, may be solved for the terms A A A A , giving precisely (5.3). To derive (5.11), it will be convenient to use the following equaton, which can be obtained with simple manipulations from the Q-system relations: The large-u asymptotics of Q • can be fixed using the first constraint in (4.43), which yields We will also need and, from (7.9), Expanding (C.5) at NLO, we find, using (C.6), (C.7), (C.8), The solution of these equations finally yields (5.11) and fixes the coefficient C as in (5.14).

D State/charges dictionary
The purpose of this Appendix is to provide a dictionary to express the charges M A ,M I in terms of the spin chain length and excitation numbers appearing in the Asymptotic Bethe Ansatz description of a generic state.
D.2 Fermionic duality: from η = +1 to η = −1 It is expected that every state (or, more precisely, every multiplet) can be represented by a regular solution of the Asymptotic Bethe Ansatz, where regular means that for every type of root x i we have x i = 0, x i = ∞. Let us now review (see Appendix A in [35]) how to switch from a regular solution of the η = +1 ABA, characterized by the roots to a regular solution of the η = −1 ABA. This type of duality transformations is well known from the N =4 SYM case [7]. Following the standard argument, we consider the polynomial in x(u): Due to the ABA equations (D.1),(D.3), we see that this polynomial has zeros at all roots of type x = x(u 3,j ) and x = 1/x(u 1,j ); besides, due to the zero momentum condition, it vanishes for x = 0. One may then write where x 3,j label the extra zeros of P (x) outside/inside the unite circle, respectively. By considering the weak coupling limit of P (x), and considering that x s,j ∼ h −1 , -55 -one may count the two new types of roots: We have then found the fermionic duality equation 31 : with an inessential proportionality factor independent of u. It is now standard to verify that the set of roots

D.3 Asymptotics of the QSC and excitation numbers
The charges entering the asymptotics of the QSC are, in terms of the number of Bethe roots in η = +1 grading: where we have denoted L = L η=−1 .

D.4 Important subsectors
In what follows we list a set of special cases corresponding to different subsectors of the theory, described by different values of excitation numbers and subsets of BA equations in η = ±1 gradings.
SL(2|1) sector: This sector can be represented by operators made of scalars Y 1 Y † 4 , covariant derivatives and fermions ψ 4+ , ψ 1 † + . The corresponding large-volume spectrum is described by the solutions of the ABA equations (D.6)-(D.10) in η = −1 grading without any auxiliary root, namely K 3 = K 1 = K 2 = 0. The classical dimensions of these operators as realized in 31 Notice that the prefactor x δ K 2 ,0 appears here due to the fact that we insisted on enumerating only regular Bethe roots in both gradings. In the grading η = +1, the description of this sector involves some of the auxiliary roots: SL(2)-like sector: Rather than a sector, this is a subset of states belonging to the SL(2|1) sector, which satisfy the condition K 4 = K4 and {u 4,j } = u4 ,j (see [77] and [34] for a detailed discussion). In this case M 5 = 0 and the ABA equations reduce to the following single equation: This set of states were studied at weak coupling using the QSC in [46].

SU(4) sector:
The operators belonging to this sector are made of all the complex scalars of the theory: Y a , Y † b , a, b = 1, . . .

D.5 Distinguished grading
Finally, a further very common form of the Bethe Ansatz equations is the one related to the distinguished Dynkin diagram. This is the form in which the 2-loop BA was originally written in [29]; it is known that it does not admit an all-loop generalization in terms of explicit functions of the Bethe roots. At two loops, one can relate the roots appearing in this version of the BA to the ones featuring in the other two versions by a chain of fermionic dualities (see [77], Appendix A). The relation between the excitation numbers in the distinguished-grading Bethe Ansatz, denoted as K d s for s = 1, 2, 3, 4,4, and the excitation numbers in the η = −1 grading, is Finally, let us make contact with the Dynkin labels [∆, j; p 1 , q, p 2 ] defined in relation to the distinguished diagram, which are widely used in the literature, e.g. [77]. In terms of these charges, the parameters entering the asymptotics of the QSC are given by

E An integral formula for P
In this Appendix we prove an exact integral formula for P, which could be useful for computing this quantity from the numerical solution of the QSC. The expression is dz e πz log −τ 4 (z) dz e πz log τ 4 (z) τ 4 (z) , where E(h) is an elementary function of h defined in (7.89). To prove (E.1), we use (4.21) to write log −τ 4 (z) τ 1 (z) = iP + A(z), A(z) = log τ 4 (z) τ 4 (z + i) , where A(z + i) = −A(z). Assuming that A(z) has no singularities on the first sheet, we can open up the integration contour circling the cut to a couple of infinite horizontal lines lying at Im(z) = ±i/2. Thus we see that the integral over A(z) exactly cancels due to the periodicity of the integrand, leading to (E.1). Notice that the ABA expression (7.88) for P is just a an application of this formula where τ 4 /τ 1 takes its large volume value, which can be read from (7.66),(7.84),(7.85),(7.87).