Continuum limit of fishnet graphs and AdS sigma model

We consider the continuum limit of 4d planar fishnet diagrams using integrable spin chain methods borrowed from the $\mathcal{N}=4$ Super-Yang-Mills theory. These techniques give us control on the scaling dimensions of single-trace operators for all values of the coupling constant in the fishnet theory. We use them to study the thermodynamical limit of the BMN operator corresponding to the spin chain ferromagnetic vacuum. We find that its scaling dimension exhibits a critical behaviour when the coupling constant approaches Zamolodchikov's critical coupling. Analysis close to that point suggests that the continuum limit of the fishnet graphs is controlled by the two-dimensional $AdS_{5}$ non-linear sigma model. More generally, we present evidence that the fishnet diagrams define an integrable lattice regularization of the $AdS_{5}$ model. A system of massless TBA equations is derived for the tachyon energy by dualizing the TBA equations of the weakly coupled planar $\mathcal{N}=4$ SYM theory.


Introduction
String sigma models are believed to provide a general solution for the sum over planar diagrams in gauge theories [1]. The AdS-CFT correspondence [2] shed light on this old idea and suggested new embodiments for conformally invariant gauge theories. The most famous example is the 4d N = 4 SYM theory which is conjectured to be dual to string theory in AdS [3]. This supersymmetric gauge theory is also special in that it is believed to be integrable at large N [4]. The latter property gives us a handle on the AdS-CFT dictionary, enabling the development of new techniques for carrying out the large N resummation of the field theory diagrams, at both planar [4][5][6][7][8] and non-planar level [9][10][11]. Furthermore, these methods allow us to explore a larger chunk of the correspondence between planar diagrams and sigma models by means of partial or twisted re-summations, which are naturally associated with some integrable deformations of N = 4 SYM.

JHEP01(2019)002
L T Figure 1. A planar fishnet diagram of order L × T . Every intersection point x is integrated over spacetime and every connecting line stands for a massless scalar propagator ∼ 1/(x − y) 2 . In the fishnet theory, the horizontal and vertical lines correspond to trajectories of the φ 1 and φ 2 particles, respectively, and the diagram is weighted by L × T powers of the coupling g 2 .
In this paper we consider such an integrable daughter of N = 4 SYM, which comes with no gauge fields nor any clear-cut stringy interpretation. The theory to be studied was introduced recently by Gurdogan and Kazakov [12] and is known as the fishnet theory, see also [13] for earlier work and [14][15][16][17] for further developments. It consists of two complex matrix scalar fields interacting by means of a single quartic coupling 1 with the trace taken over the N × N matrix indices, which here are just flavour indices. It can be viewed as a truncation of weakly coupled N = 4 SYM in which gluons and gauginos are forcefully decoupled and only two of the three complex scalar fields are retained. The proper procedure goes through the extremal twisting [12] of the γ-deformed SYM theory [20][21][22][23] which involves sending the YM coupling to zero and the deformation parameter γ to i∞, while keeping the suitably rescaled coupling g 2 fixed. Owing to this "embedding", the theory is expected to be conformally invariant and integrable for any g 2 [15,24], at least in the planar regime. In fact, the integrability of the fishnet vertex was recognized by Zamolodchikov almost 40 years ago [13]. One appealing feature of the fishnet theory is that it produces many fewer graphs than N = 4 SYM. Quite often only a single graph contributes at a given order in perturbation theory, in the planar limit. The price to pay for this massive cut is the loss of unitarity, as the strict ordering of the fields in the potential clashes with the reality of the action. Another drawback is that the duality with string theory is uncertain, for the AdS radius is naively small. Still, the fishnet theory proves to be a remarkable testing ground for integrability, which, in turn, sheds light on families of conformal Feynman integrals [25,26] and suggests new ways of evaluating them [12, 14-16, 24, 27].
The planar diagrams of the theory, the fishnet graphs, are special in that they all look locally like the L × T square lattice shown in figure 1. Accordingly, every diagram can be JHEP01(2019)002 viewed as a partition function Z L,T for a 2d vertex model [13], with the bulk spacetime points acting as classical "spins", the propagators as nearest neighbour couplings and with the graph's external lines setting up the boundary conditions. Different observables of the planar fishnet theory correspond to different boundary conditions and all the graphs obeying the same boundary conditions are summed over.
An important observation concerning the large order behaviour of the fishnet diagrams was made by Zamolodchikov [13], see also [28] for a recent discussion, who computed, using integrable vertex model techniques, the free energy per site in the thermodynamical limit log Z L,T ∼ −LT log g 2 cr , (1.2) for graphs subject to periodic boundary conditions, and found that This constant determines a critical coupling for thermodynamically large observables in the fishnet theory and one might expect, in analogy with matrix models [29][30][31][32], that a "dual" continuum description is taking over at that point. 2 In this paper we examine the thermodynamical limit of fishnet graphs using integrable methods borrowed from the N = 4 SYM theory and argue that the continuum description is given by the 2d (bosonic) AdS 5 sigma model.
Our discussion will center around the scaling dimension ∆ of the BMN operator tr φ L 1 which maps to the ground-state energy of a ferromagnetic non-compact spin chain [24].
Integrability will allow us to study the thermodynamical limit of this energy for generic coupling g by means of a linear integral equation. It will confirm the existence of a nontrivial thermodynamical scaling ∆ ∼ Lf (g), for sufficiently "strong" coupling, and the emergence of a critical behaviour ∼ √ g cr − g close to the critical point, in line with the results of [15,24] for L = 2, 3. 3 A sketch of the thermodynamical behaviour of the scaling dimension is shown in figure 2. Dualizing our equation, by means of a particle-hole transformation, will reveal the nature of the critical point and suggest the interpretation of the BMN operator as describing the "tachyon" ground-state of the AdS sigma model, labelled by the global time energy ∆ of the BMN operator and implicitly by the size L of the worldsheet. Though we will actually never cross the line where the AdS mass squared turns negative, we will stick to the name of tachyon for the dual object. 4 The correspondence (1.4) is best summarized by the formula log g 2L = log g 2L cr + E 2d (∆, L) , (1.5) 2 Note that the critical coupling does not refer to a point at which the 4d theory is becoming critical; the fishnet theory is conformal for any g 2 . It is a point at which the planar diagrams become dense. 3 The location of the branch point gcr(L) is function of the length L; in particular [15], gcr(L = 2) = 0. 4 The tachyonic domain maps to g > gcr(L), where the scaling dimension has an imaginary part, ∆ = 2 + iν; see [15,19,24] for discussions. Schematic plot of the scaling dimension ∆ of the BMN operator tr φ L 1 as function of the coupling constant g 2 , at large L. The scaling dimension is nearly independent of the coupling constant, ∆ = L, until g = 1/2 where it starts developing a non-trivial thermodynamical scaling, ∆ ∼ Lf (g). At g = g cr , the scaling function f (g) vanishes and its derivative is infinite. The physics close to this point is controlled by the 2d low energy dynamics of a large fishnet graph. The overall shape of the thermodynamical curve agrees with the findings of [24] for L = 3. The value of the scaling dimension at the branch point appears to be ∆ = 2 at both small and large L.
which relates the 4d coupling g 2 to the sigma model energy E 2d (∆, L) and shows that the vicinity of the critical point maps to low energies on the worldsheet. In this regime the AdS 5 sigma model is weakly coupled and we will be able to verify our claims directly.
We shall also test the correspondence at the level of the 1/L corrections and obtain a system of TBA equations for the tachyon ground state, valid in principle for any L and ∆. Its form will support the more general conjecture that the fishnet diagrams define an integrable lattice regularization of the AdS 5 sigma model.
The plan is as follows. In section 2 we introduce the integrability set up and derive the integral equation for the large L limit of the scaling dimension ∆. With its help we reproduce Zamolodchikov's prediction for the critical coupling. In section 3 we obtain the dual system of equations, give arguments for its interpretation as describing the tachyon in the AdS sigma model and carry out some perturbative tests. In section 4 we dualize the full system of TBA equations and compute the IR central charge. We conclude in section 5. The appendices contain a detailed analysis of the dual linear equation, a discussion of spinning operators and a brief study of the thermodynamical limit of 3d triangular fishnet diagrams.

Thermodynamics of fishnet graphs
We start with a light review of the integrability methods for computing the scaling dimension of interest, emphasizing the correspondence with the Feynman diagrams. Readers familiar with these techniques may jump directly to subsections 2.2 and 2.3 where we restrict our attention to the thermodynamical limit and re-derive Zamolodchikov's critical coupling.

Grand canonical ensemble
In the planar N = 4 SYM integrable framework the BMN operator tr φ L 1 is identified with the ferromagnetic ground state of a periodic spin chain of length L and its scaling JHEP01(2019)002 Figure 3. On the left panel we show a two-wheel diagram contributing to the anomalous dimension of the BMN operator tr φ L 1 , represented in the middle by a spin chain with L sites. The operator is local but magnified here for the sake of illustration. Each wheel corresponds to the trajectory of a φ 2 particle. On the right panel a wheel diagram is wrapped on a cylinder of size R along the global time direction and length L along the spin chain one. Taking R → ∞ projects on the ground state at the boundary. Each vertical wheel gives rise to a tower of mirror magnons which are integrated over. Each magnon is labelled by a rapidity u, or momentum p = 2u along the time direction, and a Lorentz index a = 1, 2, . . ., which comes from the partial wave decomposition of the scalar field φ 2 on S 3 . dimension ∆ maps to the spin chain vacuum energy. The same can be said in the fishnet theory. The main difference is that in the fishnet theory the vacuum is not protected and its energy ∆ = ∆ L (g) is a very complicated function of the length L and coupling g. It was studied extensively in [15,24] for small values of L and perturbatively at weak coupling for any L in [12].
For any length L > 2, otherwise see [15], the diagrams contributing to the scaling dimension in the planar limit are the wheel diagrams [12] where virtual φ 2 particles loop around the operator tr φ L 1 as shown in figure 3. A general wheel diagram is obtained by pinching the end points of the L horizontal lines in figure 1 and periodically identifying the vertical ones. Obviously, each wheel costs L powers of the coupling g 2 and, consequently, the scaling dimension admits an expansion in powers of g 2L at weak coupling, ∆ = L + O(g 2L ), referred to as the wheel expansion in the following.
In the spin chain picture, the wheels map to long-range interactions mediated by virtual magnons travelling around the chain [33][34][35], also known as wrapping or Lüscher corrections. The magnons circulating around the chain are not the familiar ones parameterizing spin waves on top of the ferromagnetic vacuum [36]. Instead they live in the orthogonal, so called mirror, kinematics where the time t is interpreted as a space direction and the spin-chain length L as an inverse temperature, see figure 3. The mirror picture, which in the general context is inherited from the dual string worldsheet theory, see [37,38] for further discussion, can be motivated by considering the problem on the 4d euclidean cylinder R × S 3 , with S 3 the unit 3-sphere surrounding the operator and with the global time corresponding to the spin chain time t. The relevant 1 + 1 dimensional picture is obtained by dropping the 3-sphere, which becomes internal, and replacing it by the spin-chain circle, R × S 3 → R × S L . 5 The partition function on this geometry returns the ground state 5 ZL would be more appropriate, but the difference is immaterial here.

JHEP01(2019)002
energy, in the large volume limit R → ∞, Z L,R (g) = e −∆ L (g)R + . . . , (2.1) where the dots stand for heavier single-trace operators, with φ 1 charge L, and where R is the length of the cylinder along the t direction. The wheel expansion of the partition function, and thus of ∆, is more naturally interpreted as the decomposition in the orthogonal, open string, channel, with the sum running over a complete basis of states along the t direction, labelled by the total number of mirror magnons T . To be more precise, each wheel gives rise to a semi-infinite family of mirror magnons, stemming from the partial wave decomposition of the scalar field φ 2 on S 3 . Suppressing the Lorentz indices, each mirror magnon ∼ ∂ a−1 φ 2 is tagged with an integer a = 1, 2, . . ., for its Lorentz spins ( 1 2 (a − 1), 1 2 (a − 1)) and, after smearing over the t direction, with a momentum p = 2u, with u ∈ R the so called Bethe rapidity. The micro-canonical contribution Z L,T (R) follows then from integration and summation over the T magnon phase space, with proper thermodynamical weights and measures; see [39] for a recent discussion of the relation between the micro-canonical partition functions and the canonical one.
The Thermodynamical Bethe Ansatz allows one to take advantage of the factorized scattering between mirror magnons and compute the bulk free energy (2.1) to all orders in the wheel expansion and for any temperature 1/L, see [37,40,41] and references therein. The latter free energy takes the usual form with p = 2 in the case at hand and with the Y functions describing the thermal distribution of the energy among magnons. In the weak coupling regime g 2 1, the gas is rarefied and the TBA equations can be expanded around the Fermi-Dirac distribution, characterized by the Boltzmann weights Y a , where a 2 is the dimension of the a-th Lorentz representation and where a (u) is the mechanical energy of the associated magnon, The departure from the free distribution is controlled by the interaction among magnons.
To leading order, it takes the rather universal form where S a,b is the dynamical factor of the magnon S-matrix, which specifies the model and reads with Γ(z) the Euler Gamma function. The dots in (2.6) capture the effect of the O(4) degrees of freedom, or scattering of Lorentz indices, and are simply absent for the lightest magnons (a = 1 or b = 1), which are Lorentz singlets. (This isotopic component of the scattering is controlled by a rational R matrix, which is function of the difference of rapidities and is coupling independent.) Formula (2.6) is a particular case of the NLO Lüscher formula for the ground state of the twisted N = 4 SYM spin chain [40], obtained by truncating the N = 4 SYM magnon super-multiplet to its φ 2 component and by sending the coupling constant of the gauge theory to zero in all spectral and scattering data. Once plugged into (2.3), it yields the wheel expansion of the scaling dimension with the obvious map "Y = wheel", such that the first term is the free wheel, the next one the double wheel, etc. It was used for comparison with the direct integration of the Feynman integrals in [12].
All the information for moving to higher orders is in principle contained in the TBA equations, which are pretty simple if not for the dots, which accommodate for the matrix degrees of freedom and cannot be spelled out without introducing an auxiliary set of variables. They will be given in their full form in section 4. Nonetheless, the general term in the wheel expansion (2.6) has not been worked out explicitly. Also, the integration over the magnons' phase space is nearly impossible in general. A more powerful analytical treatment relies on the Baxter equation [24], which relates to the (twisted) N = 4 SYM Quantum Spectral Curve [42,43] and enables higher loop computations of the scaling dimension at finite length L. We will not need so much improvement in this paper however, since in the thermodynamical limit L → ∞ the TBA equations (2.8) simplify drastically. A distinguished feature of the fishnet TBA equations is that all the dependence on the coupling constant g 2 comes along with the energy, like in (2.8) or equivalently (2.4). Put differently, the 4d coupling of the fishnet theory is a fugacity for the number of mirror magnons and we can think of the scaling dimension ∆ as defining the free energy of a grand canonical ensemble at chemical potential h = log g 2 . This parallels the fact that the mirror magnons are in one-to-one correspondence with the wheels. Hence, one can easily probe fishnet graphs of arbitrarily large orders and obtain information about their continuum limit by playing with h.

Thermodynamical limit
The thermodynamical limit L → ∞ is uninteresting at weak coupling, since the wheels are heavily suppressed, see equations (2.4) and (2.6), and the scaling dimension ∆ = L JHEP01(2019)002 up to exponentially small effects. The situation changes drastically as soon the chemical potential h = log g 2 gets bigger than the "mass gap" 1 (0) = log 1/4. Above this threshold the s-wave magnons, with a = 1, start filling a Fermi sea, see figure 4. The filling is strict in the limit L → ∞ with all the modes outside the sea being unoccupied, with θ the step function and with the pseudo energy defined here in such a way that χ 1 > < 0 inside/outside the sea. The Fermi rapidity B fixes the edges of the sea and is determined by the condition χ 1 (±B) = 0. Moreover, as long as g 2 is smaller than the masses of the higher a magnons, that is, naively, g 2 < 1, the gas is mono-atomic and consists solely of fundamental magnons with a = 1. As we shall see later on we will never reach the next threshold, so in the following we drop the Lorentz index and assume that a = 1. An immediate consequence of having a Fermi sea, obtained by plugging (2.9) into (2.3), is that the scaling dimension develops a non-trivial thermodynamical scaling where f (h) is the free energy density of the gas at chemical potential h. The same linearization is observed at the level of the TBA equations (2.8), see also section 4, which can be written concisely as with (u) = log (u 2 + 1/4). Here, for convenience, we split the scattering kernel, into its "boost" invariant component, K, and the rest, k, which depends on a single rapidity. Their explicit forms follow from equation (2.7), with a = b = 1, and with ψ(z) = ∂ z log Γ(z) the digamma function. Note that both are even functions. The k part of the kernel merely renormalizes the chemical potential and was absorbed into the constant C,

JHEP01(2019)002
This system of equations, with the boundary condition χ(±B) = 0, admits a unique (even) solution, which can be constructed iteratively for finite value of B. The numerical solution for the free energy density is shown in figure 6. Finally, notice that the free energy density (2.11) can also be read out from the large u behaviour of the pseudo energy, which, according to (2.12) and (2.14), scales like This way of obtaining f will turn out to be useful later on.

Critical coupling
Let us now look for Zamoldchdikov's scaling (1.2). First we note that at non zero B the typical number of magnons T in the ensemble is always large, since the gas has a finite particle density, This is necessary for matching with the scaling (1.2), but the condition is not enough to get an actual match. The continuum limit also requires that a low-energy approximation be taken w.r.t. the fishnet Hamiltonian δ/δT . Since, heuristically, δT ∼ j, we expect that the fishnet dynamics will freeze at large magnon densities, that is when j → ∞. Another way of seeing it is that the scaling (1.2) is a statement about the micro-canonical energy density which, in our thermodynamical variables, translates into the requirement that for h → log g 2 cr . This behaviour cannot be observed at small B, that is for a dilute gas, j ∼ 0, since then ε ∼ f ∼ 1. On the contrary, the critical regime appears for B = ∞, when all the energy levels are filled, which again means that the density is infinite, j = ∞.
We can verify it explicitly using the integral equation. Denoting by χ cr , C cr , etc, the limiting values at B = ∞, we get to solve where * denotes the convolution over the full real axis,

JHEP01(2019)002
and using the Fourier integrals immediately tell us that Hence, the free energy density (2.11) vanishes 25) in accordance with our previous discussion. The transformation back to rapidity space yields the critical pseudo energy which is positive definite and decays exponentially quickly at infinity. Plugging it back into (2.20) and taking a large rapidity limit fix the constant of integration, C cr = lim u→∞ (χ cr + − K * χ cr ) = 0. The critical coupling follows from that condition, after recalling (2.15),

Dualization and AdS sigma model
Having reached the critical point, we want now to study its neighbourhood, corresponding to a large but finite Fermi rapidity B. Since "almost all" of the energy levels are filled, it is convenient to analyze this regime by means of a particle-hole transformation, which flips the notions of filled and empty states. As we shall see, the dual equation, the one for the holes, is of a totally different nature and lends itself the interpretation of a thermodynamical equation for the tachyon of the hyperbolic sigma model.

Particle-hole transformation
Formally, the particle-hole transformation amounts to introducing the dual kernel  Low energies mean small rapidities and all the energy levels inbetween ±B are occupied. The sea level rises upon increase of the chemical potential h = log g 2 . On the right, the dual Fermi sea, χ ∼ E. Low energies mean large rapidities, the dual excitations accumulate at infinity and all the energy levels beyond the Fermi rapidity are filled. There is no chemical potential here. Instead, the filling is controlled by the charge density ρ = ∆/L, which pilots the large u asymptotics χ ∼ −ρ log u 2 . When B → ∞, this density vanishes, the dual sea is empty and the fishnet freezes. This happens at the critical coupling h = log g 2 cr .
and acting on both sides of (2.12) with (1 − K * ). Straightforward algebra gives then where the convolution is now supported on the complementary support and where I dual = (1 − K * )I is the dual of the driving term I(u) = C − (u). The procedure is a bit formal since K scales logarithmically at large rapidity and thus the self-convolutions in the r.h.s. of (3.1) are not well defined. Fortunately, one reaches the same point by defining the dual kernel more implicitly, as the solution to Taking derivative of this equation, going to Fourier space and fixing the constant of integration yield coincides with the S-matrix of the O(6) sigma model in the symmetric channel [44,45]. It is written here in terms of the Bethe rapidity u, which relates to the sigma model hyperbolic rapidity by θ = πu/2. Formula (3.4) hints at a connection between the dual model and the O(6) sigma model. However, this cannot be the full story, as we shall discuss shortly. The next essential piece of information comes from I dual . Using the normalization of the dual kernel,

JHEP01(2019)002
one concludes that the constant C, and hence the chemical potential, drops out of the dual equation, (1 − K * )C = 0, leaving us with where we used (2.20) and (3.3). The dual driving term is thus simply given by the critical pseudo energy (2.26). It acquires here the meaning of a dual energy, E = χ cr . As noticed earlier, this one decays exponentially at large rapidity. Therefore, the dual energy describes a gapless particle, since one can lower arbitrarily the energy of a dual excitation, by sending it to larger and larger rapidities, where m = 4 √ 2 sets a reference energy scale. This is in line with the fact that the dual Fermi sea has a non-compact support, see figure 4. The dual low energy modes accumulate at infinity, which is typical for gapless systems, see e.g. [46][47][48][49][50][51].
It is also natural in the dual picture to exchange the roles of the thermodynamical quantities. The lack of a dual chemical potential, for instance, invites us to view the free energy density f as being part of the specification of the state. Namely, we can simply think of ρ = f = ∆/L as the charge density that pilots the large u behaviour, (3.9) One can also say that it triggers the formation of the dual Fermi sea, as ρ and B play interchangeable roles. 6 In particular, the approach to the critical point B → ∞ corresponds to the low density regime ρ → 0. The last important quantity is the 4d coupling log g 2 , which, at the moment, is buried inside the constant C. Fortunately, one can substitute to (2.15) the more transparent relation log g 2 = log g 2 cr + where is the critical micro-canonical distribution density, which is Legendre conjugated to χ cr , and solution to the integral equation

JHEP01(2019)002
Formula (3.10) is obtained by integrating both sides of the original equation (2.12) against φ cr , Here, in the first line, we used and transferred the action of K from χ to φ cr , then applied (3.12) in the second line and finally used (2.15).
Introducing then a dual momentum P , by means of a 90 o rotation θ → θ + iπ/2 of the energy, and noting that φ cr = P = ∂ u P , one concludes that in the dual description the coupling of the 4d theory is simply related to the 2d energy of the state through the formula (1.5) with, in the thermodynamical limit, As we shall see later on, formula (1.5) is quite general, and not restricted to the thermodynamical limit. Finally, let us add that putting energy and momentum together yields the dual dispersion relation, which reads after eliminating the parametric dependence on the rapidity in (3.15) and (3.8). This formula makes the square lattice and its symmetries manifest: it preserves a Z 4 subgroup of euclidean rotations, generated by E ↔ iP , and shows a maximum momentum, P = ±π, and a maximum energy, E = 2 log (1 + √ 2). Both features disappear at low momentum where one recovers the dispersion relation for a massless relativistic particle. The dispersion relation is depicted in figure 5, in the 1st Brillouin zone. It is analogous to the energy of a spinon above the anti-ferromagnetic vacuum of the XXX spin chain and, following this analogy, we could say that at the critical point the spin chain settles down in its symmetric vacuum.  Figure 5. Plot of the dual energy E(P ) for P ∈ (−π, π). At low momentum, we have a massless spectrum E ∼ |P |. The energy reaches a maximum at P = ±π.

Sigma model interpretation
Let us come to the interpretation of the dual equations. They are very similar to the thermodynamical equation for a (zero temperature) finite density gas of particles, carrying maximal U(1) charge, in the non-linear O(6) sigma model [52], where θ = πu/2. The difference only comes from the dispersion relation which, in the relativistic low momentum approximation, amounts to substituting and reversing the support of the distribution. It makes a big difference for the interpretation. In the O(6) model, the particles are massive and though the scattering kernel is repulsive the particles remain confined on a compact support. The kernel has the same effect in our case but the energy is not gapped and the potential runs away. Therefore, at finite charge density ρ, the excitations start filling the energy levels around u = ±∞ and spread in the opposite directions, towards smaller rapidities. Also, in the compact case, the charge density ρ matches with the particle density, obtained by integrating the distribution χ over its support. In our case, the support is non-compact and the distribution is not normalizable, suggesting that the gapless excitations in the condensate cannot be counted. Given the symmetries of our problem, the most natural guess is that we are dealing with the AdS 5 sigma model. This non-compact model is known for not developing a mass gap and for having a continuous spectrum in finite volume. Furthermore, the change (3.19) in the energy for a given scattering kernel K embodies the "inversion of the RG flow", which is the formal perturbative way of relating the sphere and the hyperbolic sigma model. We discuss it in more detail below.
Finally, note that there are similarities with the equations obtained for massless factorized scattering theories [46][47][48][49], see also [50,51]. In our case, since we cannot enumerate the particles in the condensate, it is not clear whether we can talk about an S matrix.

JHEP01(2019)002
Put differently, we do not think of our equation as describing the continuum limit of a dense but fundamentally discrete distribution of Bethe roots. On the contrary, the distribution is fundamentally continuous, and will remain continuous after introducing finite size corrections. It defines a one-parameter family of ground states, labelled by ρ, or better ∆.

Perturbative analysis
We can verify the interpetation of the dual equation by comparing its predictions against a direct finite density calculation in the sigma model. This is standard analysis for sigma models. It was carried out through two loops in [52,53] for the sphere. We can follow the same lines for the hyperboloid. In fact, the results for the sphere carry over to the hyperboloid, since compact and non-compact models only differ perturbatively in the sign of the coupling constant, as expected on geometrical grounds [54][55][56]. We recall how this comes about below.
We consider the non-linear sigma model in Minkowskian AdS d+1 , where d = 4 in our case. Its 2d action is given by where the embedding coordinates X A = X A (σ, τ ) take values on the hyperboloid and where X = (X 1 , . . . , X d ) are d transverse directions. The worldsheet metric is taken to be flat, with signature (−+), and we assume periodic boundary conditions in σ, σ ∼ = σ + L.
The coupling e 2 sets the curvature of the hyperboloid. In stringy notation e 2 = 2πα and the model is weakly coupled when e 2 ∼ 0. The theory has exact SO(2, d) symmetry, with associated conserved Noether currents and A, B = 0, . . . , d + 1.
The model is classically integrable and presumably quantum integrable, perturbatively, for the arguments supporting the integrability of the sphere [57] also work for the hyperboloid. The integrability of the non-compact model remains puzzling at the quantum level however. The model's spectrum has no good particle interpretation and thus cannot be handled by the conventional factorized scattering methods. On top of that the model has well known problems in the UV and might have to be completed. However, none of these complications really is a limitation here, as we do not need particles to make sense of our equation and we have a lattice to make sense of the model at short distances. Nonetheless, it is interesting to note that both difficulties somehow relate to the running of the coupling e 2 , which, owing to the negative curvature of the AdS space, is governed by a positive beta function, see [55] for a recent discussion, (3.23)

JHEP01(2019)002
It says that perturbation theory can be trusted at low momentum and this is more than enough for what we intend to do here. Getting back to our problem, we seek a state with a charge ∆, along the global time direction, which is uniformly distributed along σ, with J α the Noether current (3.22), and which corresponds to a local extremum of the sigma model energy, The most natural candidate is the "tachyon", that is, a point-like and time-like geodesics at the center of AdS. It is given classically by the σ independent solution 26) or, equivalently, t cl = Hτ , where t is the global time coordinate. Owing to the signature of the AdS space, "excitations" along the time-like direction contribute negatively to the sigma model energy. This applies to our state that has classical energy and charge density with a frequency H which is negative for ∆ positive. Eliminating H, we obtain that the energy is quadratic in ρ, E 2d /L = −e 2 ρ 2 /2 , As we shall see, the agreement gets even better when the running of the coupling and the perturbative corrections are included. Loop corrections are more efficiently computed by exploiting the thermodynamical nature of the state under consideration. Like for the sphere [52], one can access to the energy density by coupling the model to a constant electric field H and extremizing the path integral over an Euclidean worldsheet. Unlike the sphere and owing to the indefinite signature of the target space, one must be careful with the Wick rotation. The kinetic term of the global time coordinate comes out with the wrong sign. We evade the problem by rotating the global time coordinate along with the worldsheet one, 7 It brings us to a doubly euclideanized partition function, which is perturbatively well defined to any order in e 2 . It remains then to covariantize the τ derivatives, with X d+1 E = iX d+1 the Euclidean coordinate, and expand at weak coupling around the center of the space, a.k.a. Goldstone vacuum, where the fields y and z = (z 1 , . . . , z d ) are canonically normalized. It gives where L E is the Euclidean Lagrangian density and where the dots stand for cubic and higher couplings. The first term in (3.32) is the classical free energy density F (H)/L. The next one shows that the d transverse excitations z acquire a mass |H| in the tachyon background, as in the compact case [52]. In fact, everything is as for the sphere model up to e → ie. The substitution flips the signs below the square root in (3.31), resulting in the compactification of X 0 , and turns the derivatives (3.30) into which is the canonical way of boosting the system along a big circle ⊂ S d+1 , see [52,53]. Therefore, in agreement with the discussions in [54][55][56], the compact and non-compact problems are the same, perturbatively, if not for the sign of the coupling, e 2 ↔ −e 2 . At this stage, we could import the result directly from the perturbative studies of the sphere sigma model [52,53] by doing the continuation to negative coupling. In particular the one loop free energy is the same in the two cases and comes directly from the determinants for the quadratic actions in (3.32). Their evaluation using dimensional regularization gives, in the M S scheme, where e 2 = e 2 (µ) solves (3.23), with Λ the M S scale. Taking the running of the coupling into account, the free energy (3.34) is verified to be independent of the subtraction scale µ. The energy and charge of the state are obtained by a Legendre transformation, which, after plugging (3.34), eliminating the coupling (3.35) and fixing µ = 2πρ/d, yields The two-loop calculation [53] would also give us the next contribution ∼ 1/ log 3 (Λ/µ). Importantly, since the coupling has disappeared, the same formula (3.37) applies to both the sphere and the hyperboloid. The sole difference is that the expansion is valid for ρ Λ in the case of the hyperboloid and for ρ Λ in the case of the sphere. This is the "inversion of the RG flow" alluded to before.
Equation which matches with the Λ-to-m ratio of the O(6) model [52]. We recall that in our case m is not a mass gap. More generally, as shown in appendix A, our integral equation turns out to be identical at large B with the equation (3.18) for the compact sigma model, after continuing B to −B. Since the Fermi rapidity plays the role of the radius of the AdS space, changing its sign has the same effect as changing the sign of the coupling e 2 . Hence, we can bypass the direct comparison with the perturbative sigma model and, assuming the validity of the equation

TBA equations and central charge
The finite density equation only probes a diagonal subsector of the 2d theory and as such might miss some features of the model. More compelling evidence for our proposal can be found by looking at the finite size corrections. We will see that the observations made earlier at the linearized level uplift to the full set of TBA equations. We will then discuss briefly the finite size corrections to the tachyon energy level.

Massive TBA
First recall the original form of the TBA equations for the ferromagnetic vacuum tr φ L 1 . It is obtained by taking the fishnet limit of the system of TBA equations for the ground state of the twisted N = 4 SYM spin chain [40]. The relevant symmetry group is the Lorentz group O(4) and the relevant excitations describe the Lorentz harmonics of the scalar field φ 2 , introduced in section 2. These modes appear on an equal footing in the TBA with each mode mapping to a massive (non-relativistic) magnon with bare energy a , see (2.5), and thermodynamical weight Y a , see (2.6). These Y functions are subject to the equations where L = log (1 + Y ) and * denotes the convolution, defined in (2.21). The kernel K a,b controls the part of the interaction that depends on the difference of rapidity, see (2.7), and K 1,1 = K, with K the kernel used before, see (2.14). The driving term in (4.1) is given in terms of the bare energy (2.5) and of a constant C, which does not depend on the mode number a. The latter constant captures the dependence on the coupling constant of the 4d theory and absorbs the part of the kernel that depends on a single rapidity. It reads where k a = 2ψ generalizes (2.14) to a 1. The last sum in the r.h.s. of (4.1) describes the couplings to the matrix degrees of freedom, represented in the form of O(4) ∼ SU(2) + × SU(2) − dispersion-less magnons, with wave functions Y a,± . They are labelled by the dimensions a = 2, 3, . . . of the SU (2) representations. The interactions between momentum carrying magnons and isotopic ones are controlled by the scattering kernels of the XXX spin chain, (4.5) Note in particular that K 1,b = 0, meaning that there is no coupling between the s-wave magnons and the isotopic ones, as expected. The TBA equations for the isotopic magnons are entirely controlled by the kernels (4.5) and read log Y a,± = − b 2 which follows from the universal logarithmic scaling of the kernels K a,b ∼ 4 log u and energies a ∼ 2 log u. The non-universal part log m a of the sub-leading behaviour comes from the isotopic Y functions Y a,± . The latter tend at large rapidity to the constant solution of the O(4) Y system, which will be given later on, see equation (4.37). Plugging this solution into (4.1) and using the normalization of the XXX kernels, one recovers that m a = a 2 gives the dimension of the a-th Lorentz representation, see (2.4).
The relation between global quantum numbers and large rapidity asymptotics is well known in the spin chain context, see e.g. [58], and plays an essential role in the Quantum Spectral Curve [42]. Equipped with the full set of TBA equations we can verify the claims made earlier about the thermodynamical limit L → ∞ for 1/2 g g cr . In particular, we can check that the higher modes Y a>1 stay under control in the presence of the condensate log Y 1 ∼ Lχ, all the way to the critical point where the back reaction is maximal. Plugging the critical pseudo energy, log Y 1 = LE 1, inside (4.1) and using the identity (4.18) reveal that the functions Y a>1 become L independent at the critical point, see (4.38). This behaviour is indicative of a symmetry enhancement and is further discussed below. Nonetheless, it does not change the fact that the higher modes are negligible thermodynamically. Also, the isotopic Y functions do not couple to the length L, nor to Y 1 , and though their actual values depend on g they remain of order O(L 0 ) all the time, see (4.37) and (4.38). Therefore, in the thermodynamical limit one can legitimately substitute L b → δ b,1 log Y 1 θ(B 2 − u 2 ) and L b,± → 0 inside (4.1) and this way recover the linear integral equation (2.12). 8

Massless TBA
The dualization introduced earlier can be applied to the TBA equations. Let us start with the dual singlet equation, which extends (3.2) beyond the thermodynamical limit. As before, it follows from acting with (1 − K O(6) * ) on both sides of the a = 1 TBA equation (4.1). Then, combining where L 1 = log (1 + 1/Y 1 ) and where K b,2 is the XXX kernel defined in (4.5). The rationalization of the couplings between the singlet and the higher modes is the first indication of the restoration of the full symmetry. Notice also that the new TBA equation (4.11) is consistent with the logarithmic scaling (4.7), thanks to the normalization of the O(6) kernel (3.6).
To complete the picture we also verify that the higher modes become dispersion-less in the new vacuum. Indeed, as shown below, the equations for the a > 1 modes can be brought to the form which shows no coupling anymore to the length L of the system. Together with the equations (4.6) for the O(4) magnons, which stay untouched in the dual picture, these equations form an O(6) system of magnons sourced by the "vector" node Y 1 . The derivation of the dual equations (4.12) relies on the relations which are valid for a, b > 1. Together with (4.9), they allow us to re-write the first term in the r.h.s. of (4.1) in the form where we introduced It remains then to show that rest a = L a − C , (4.16) which is the statement that the driving term cancels out after shifting to the new vacuum. This can be done in two steps: (1) evaluate b>1 K 2,b * L b using (4.11), plug the expression in (4.15) and simplify the result, such as to get and (2) apply the identities K a,1 * LE = L a , (4.18)

JHEP01(2019)002
and which both hold for a > 1. This is readily seen to imply (4.16). Let us add that though formula (4.18) is straightforwardly derived, the following one (4.19) requires more attention.
It is tempting to open up the brackets in the l.h.s. of (4.19) and apply (4.9). This algebra is not correct however and misses the constant C. The problem is that K a,1 * log Y 1 and K a,1 * K O(6) * log Y 1 are not separately meaningful, because log Y 1 and K O(6) * log Y 1 scale logarithmically at large u, and so does the kernel K a,1 . However, the large asymptotic behaviours cancel out in their difference as in the l.h.s. of (4.19). To reproduce (4.19) one can first strip out the problematic part by writing where Y 1 = Y 1 e ∆ −C → 1 + O(1/u 2 ) at large rapidity, and then apply (4.9). It yields after using that K a,2 * = a and K a,2 * C = C for a > 1. Finally, adding gives (4.19) and completes the proof. The TBA energy is defined canonically by where P is the dual momentum introduced in (3.15). It is related to the coupling constant of the 4d theory by means of (1.5) as found earlier in the thermodynamical limit. The proof is similar to the one outlined below equation (3.10). Namely, one integrates both sides of the TBA equation (4.1) for log Y 1 against φ cr (u) = P (u), transfer the action of the kernels to φ cr , and use the identities as well as the definition (4.3) and the relations (3.14). Equations (4.11), (4.12) and (4.6) form the dual set of TBA equations. They take the same form as for the ground state of the O(6) sigma model if not for the driving term, which is given by the gapless energy (3.8). This feature allows us to impose the boundary condition (4.7) at large rapidity, or equivalently small momentum, which in the dual picture is part of the specification of the state. Therefore, as seen before, the TBA equations do not describe a single isolated vacuum, as in the compact case, but a one-parameter family of vacua labelled by the scaling dimension ∆. The sigma model energy (4.23) is then determined as a function of the system size L and quantum number ∆. Formula (1.5) is the only place where the coupling constant of the 4d theory enters. Figure 7. Summary of the massive and dual massless TBA equations, on the left and right panel, respectively. Each node represents a Y function. Filled nodes stand for the momentum carriers while white nodes denote dispersion-less magnons. On the left panel, the momentum carrying nodes form the array Y a , with a = 1, 2, . . . and with a = 1 being the singled out node. This array is connected on both sides to an SU(2) system of magnons, if not for a = 1, which is a Lorentz singlet. In the dual description an O(6) tail of magnons is attached to a single momentum carrying node, which is in the vector of O(6).

JHEP01(2019)002
Note finally that although our model is gapless we do not witness the splitting into left and right movers that is characteristic of massless factorized scattering theories [46,47]. This is due to the fact that we are dealing with a lattice completion of the sigma model. The 2d boost symmetry is broken at short distances, not at the level of the scattering kernels, which are as in the O(6) model, but because of the dispersion relation, see figure 5. The latter permits the left and right movers to join up at high energy. Alternative UV completions might differ in the way left and right movers are put together into the TBA equations. It would be interesting to investigate what are the possible options if one insists on having exact relativistic symmetry.

Y system
The two formulations of the TBA equations are summarized in figure 7. As a final check of their equivalence one can verify that they give rise to the same Y system. The latter refers to functional relations that the Y functions must obey. One verifies, using fusion relations among the TBA kernels, that the O(6) tail of equations takes the universal form for a > 1 and s = 0, ±1, while the equation for the singlet, a = 1, is more model dependent and reads Here, Y a,s=0 = Y a , Y a,s=±1 = Y a,± , Y ± (u) = Y (u ± i 2 ), etc. They are seen to match the Y system equations for the O(6) sigma model [59][60][61] up to change of notations.

IR central charge
Since our model is gapless, the typical energy scale for the finite size corrections is 1/L. Hence, the thermodynamical description no longer applies when the charge density ρ be-

JHEP01(2019)002
comes comparable to 1/L 1. On the sigma model side, we can get to this far-infrared domain by sending L → ∞ at fixed ∆ and then read out the energy level from the 2d scaling dimension of the tachyon vertex operator, using the operator-state correspondence [54,62]. Alternatively, we can push the thermodynamical analysis to the very low density domain ρL ∼ 1 by including the finite size corrections to the one loop determinants. We will follow the latter approach here.
Recall that we found 1 longitudinal massless boson and d transverse massive bosons, with a mass −H set by the chemical potential H (which is negative). The sole effect of the massless boson is to contribute one unit of central charge to the free energy The contribution of the d transverse bosons is more interesting and takes the form where K 1 is the 1st modified Bessel function of the second kind. It is exponentially small at large L, for |H| ∼ 1, as expected. At very low density, |H|L 1, when the bosons become effectively gapless, it can be expanded as with γ E the Euler-Mascheroni constant. Besides the Casimir energy, we find a term linear in H, which has the interesting effect of shifting the extremum value of ∆, at H = 0, The next term in (4.29) is quadratic in H and resembles the bulk term (3.34). Its role is to switch the argument of the coupling from |H| to 1/L. Adding up the bulk and finite size contributions and performing the Legendre transformation, one obtains the energy where e 2 ∼ 2π/(d log L) and c = d + 1. This formula formally agrees with the CFT prediction for the tachyon energy [54,55,62], being the sum of the Casimir energy and of the tachyon kinetic energy, which is yet another Casimir, for the quadratic invariant of the tachyon representation this time. However, as written the formula is not complete and misses the loop correction to the central charge; the latter is of the same order as the kinetic part but originates from the two loop correction to the free energy F (H) at H = 0.

JHEP01(2019)002
This addition has a minor effect though and to summarize the main points are that (a) the extremum of the sigma model energy is shifted by the finite size effects to ∆ = d/2 in agreement with the observations made in [15,24] and (b) the coupling e 2 does not run all the way down to zero and freezes at the IR energy scale µ = 1/L 1. Notice also that the extremum energy is not simply given by the would-be "vacuum energy" −πc eff (L)/6L and gets corrected by the kinetic term, which moves closer to the findings of [24] and fits remarkably well the finite L numerical points obtained from the Baxter equation. 9 The far-infared domain is significantly harder to study on the TBA side. Although it might be possible to address it using advanced TBA techniques, like the ones developed in e.g. [46,[63][64][65][66], here we will content ourselves with running the standard dilog routine for computing the central charge [64,[67][68][69][70].
The main simplification that comes about at large L is that the solution splits into decoupled left-and right-moving kinks, centered around ± log ( 1 2 mL) respectively. Each kink corresponds to a scale invariant solution of the TBA equations and, by parity, the two kinks contribute the same amount to the energy, which takes the scale invariant form E 2d = −πc T BA /6L. Remarkably, the exact shape of the kink appears to be immaterial at first order and the central charge c T BA is directly determined by the boundary values of the Y functions at the extremities of the kink, denoted by Y 0 and Y ∞ in the following. The precise relation is is a sum over all the nodes of the Y system and where L denotes Rogers dilogarithm, In our case the kink interpolates between the two phases described earlier. Each is characterized by a symmetry group and an associated constant solution to the Y system (4.25), with vanishing boundary conditions at a = ∞. At large u we are deep inside the dual sea, the symmetry is broken down to O(4) and so the relevant solution is the stationary (4.37) 9 We are very thankful to David Grabner and Kolya Gromov for sharing with us their findings for the critical coupling at L = 4, 5 and for comparing our asymptotic expression against their numerical results.

JHEP01(2019)002
On the other hand, fixing the rapidity and sending L to infinity place us far from the dual sea, where the O(6) symmetry is restored. In this case we need the constant A 3 solution, .
(4.38) (Note, in passing, that this solution is telling us how the Y functions for the higher modes behave close to the critical point; they tend to constants Y a < 1.) Evaluating the dilog sum (4.35) numerically on these two solutions, one infers that 39) and recovers that c T BA = 5, as expected. Since this analysis is not sensitive on the way the Y functions decay at infinity, it does not capture the information about the scaling dimension ∆. Note finally that the same algebra applies to the compact model, with the same central charge. Generalizing to d dimensions, one expects that a scale invariant solution interpolating between O(d + 2) and O(d) boundary conditions will reproduce the central charge of the O(d + 2) sigma model. This expectation can be verified numerically for even d using the constant solution to the D r Y system given in [59] (see also references therein).
This computation relates to earlier TBA analyses for the central charge of the sphere sigma models [59,61], although the details of the calculations are a bit different. The latter references considered the integrable deformation of the level k coset WZNW model, conjectured in [61] to approach the sphere sigma model in the limit k → ∞. This deformation translates at the level of the Y system into a truncation at a = k, where a hard wall Y k,s = ∞ is located. Although the kink solutions for the truncated system tend to the undeformed ones when k → ∞, one finds a remnant of the wall in c 0 and c ∞ even after sending k → ∞. The latter coefficients indeed receive additional contributions from the wall region a ∼ k and the excess only drops out in the physical central charge c = c 0 − c ∞ .
It would be interesting to study the truncation of the TBA equations for the noncompact model as well as its field theoretical realization along the lines of [61]. A truncated system would certainly be more tractable both numerically and analytically. One puzzling question that comes to mind is whether a truncation can be found for every ∆. Indeed it might not be possible for a truncated system to support the logarithmic behaviour (4.7) at u = ∞. In fact, looking at the expressions obtained in [59] for the compact model, one might expect all Y functions to approach non-zero constant values at u = ∞, for finite k, although some of them can be made arbitrarily small at large k, log Y a ∼ −d log k . (4.40) The similarity of this scaling with the ∆ = d/2 asymptotic behaviour (4.7) suggests that the tachyon will settle down at the extremum of the energy in the k → ∞ limit. More generally, it could be that the truncation is only possible for certain quantized values of the spectrum and comes along with a certain compactification of the target space.

Conclusion
In this paper, we applied integrability to the study of the thermodynamical limit of the 4d planar fishnet graphs. The general proposal is that the fishnet diagrams correspond to an integrable lattice regularization of the AdS 5 sigma model. We tested this correspondence perturbatively for the scaling dimension of the BMN operator tr φ L 1 which maps to the tachyon on the sigma model side. It would be very interesting to see if the correspondence can be understood within the formalisms of [13,24] for these ones relate more directly to the Feynman graphs.
Our discussion fits with the familiar correspondence between spin chains and sigma models in condensed matter physics, which relates in its emblematic form the antiferromagnetic Heisenberg magnet to the O(3) sigma model at topological angle θ = π. Our set up appears orthogonal to the traditional one however, for (1) the symmetry group is non compact, (2) the vacuum is ferromagnetic and (3) the interactions are non local. Ferromagnetic spin chains come usually with a non-relativistic spectrum [71]. What saves us from that fate is that our "magnet" loses its ferromagnetic property at the critical point, as indicated by the vanishing of the "scaling dimension per field" ∆/L in the thermodynamical limit. (The latter is an order parameter for the conformal symmetry breaking.) This loss comes along with the emergence of the symmetric phase described by the AdS sigma model. It is not clear to us whether the non locality of the interactions is essential in this respect. In any event it would be interesting to see how generic the correspondence is by considering other integrable fishnets [13,14,16,25,72], with dimension d = 4, deformed propagators, Yukawa couplings, etc. In particular, the d = 3 triangular fishnets can be embedded into the ABJM theory [14] and studied along the lines of this paper. The linearized analysis is carried out in appendix C and hints at a connection with the AdS 4 sigma model. It would also be interesting to explore the potential connections with the work of [73] which demonstrates that certain non-compact sigma models can be engineered from the continuum limit of anti-ferromagnetic spin chains with compact symmetry group but non-Hermitian dynamics.
Another important question is whether all the local operators of the theory find room in the sigma model description. Some operators are known to be protected and thus seem to evade the sigma model. The excited operators look more promising. In appendix B we give evidence that adding spin to the operator is in line with adding spin to the tachyon. A more systematic analysis would be needed for comparison with the full spectral curve of the classical AdS 5 sigma model, along the lines of [74,75].
The relation (1.5) between the 4d coupling g 2 and the sigma model energy is also worth a few comments. First one notices that this coupling decouples from the story and only enters in the latter relation to the energy. This is not too surprising given that the fishnet diagrams themselves do not know about it. The latter is inserted in the sigma model in the form of a cosmological constant and connects to the energy because of the summation over the discrete modular time, or Schwinger parameter, T 0 (g/g cr ) 2LT e −T E 2d (L,∆) = 1 1 − (g/g cr ) 2L e −E 2d (L,∆) . (5.1)

JHEP01(2019)002
It explains why the relation to the 4d coupling takes the form of a "marginality condition", for the on-shell states of the fishnet theory are associated to the poles of the "string propagator". The condition is exact, being valid for any L, but only at large L is the sigma model weakly coupled. Nonetheless, even at small L one should be able to write (part of) the spectrum in this implicit form. This is in line with the observation made in [15] for L = 2, where a four-point function of short operators was cast into the form of an integral over a continuum of scaling dimensions, with physical states sitting at the poles of the integrand. The physicality condition in this case was found to take the simple form , for spin-less operators. More generally, the marginality condition of the sigma model should relate to the eigenvalues of the graph building operator [24]. Note that similar representations were also found for the correlators of the SYK model [76], which is graphically a close relative of the fishnet theory. It remains to be seen if this on-shell condition admits a genuine stringy interpretation. String worldsheet theories with a free tunable intercept do exist in flat space, if one quantizes the Polyakov actionà la Gupta-Bleuler [77]. These non-critical worldsheets do not seem to lift to consistent interacting string theories, but they seem to be reasonably well behaved classically. An estimate of the asymptotic number of states in the fishnet theory or relatedly of the Hagedorn temperature [78,79] could shed light on this issue.
It would also be interesting to see what the sigma model can teach us about the structure constants and the higher point functions. It might naturally relate to the integral representations obtained in [15] and [76].
Finally, our analysis could also find applications in other critical corners of the N = 4 SYM theory. For instance, the critical behaviour of scaling dimensions in the fishnet theory is reminiscent of the tachyonic instabilities discussed in [80]. Also, the BFKL limit, which relates to the point in the spin plane where the derivative of the scaling dimension diverges, is another possible place where a contact could be made; see [81] for recent discussion. More generally, one might be able to find similitudes with observables that are dominated by generalized ladder dynamics, like the one recently considered in [82].

JHEP01(2019)002
useful feedback. B.B. would like to thank Nordita for hospitality during the last stage of the writing of this paper. This work was supported by the French National Agency for Research grant ANR-17-CE31-0001-02.

A Low density analysis
In this appendix we study the dual integral equation at large B, corresponding to a low charge density ρ m. We shall demonstrate that its solution matches the one for the compact sigma model at every order in the 1/B expansion, up to B → −B. Since B relates to the inverse of the sigma model coupling, switching its sign embodies the change of curvature which maps the sphere into the hyperboloid. A similar rule maps the UV and IR regimes of the finite density equations discussed in [46,47] related to marginally relevant and irrelevant (integrable) deformations of SU(2) 1 WZNW model.
We consider the equation for the derivative of χ, which is technically simpler. It reads dχ du where C = {u ∈ R : u 2 > B 2 } is the support of the dual distribution, E is the derivative of (3.8), (3.4). The equation is obtained by differentiating (3.2), performing an integration by part and using that χ(±B) = 0. Notice also that χ is odd and scales like at large u, where ρ = ∆/L is the charge density. This equation can be analyzed systematically at large B by following closely the method developed by Volin in [83,84] for the compact sigma model. The reader is also referred to [52] for an earlier study based on the Wiener-Hopf method.
To begin with, one introduces the resolvent which is an even analytical function of u if not along the contour C where it has the discontinuity with θ the step function. It is required to scale like R(u) ∼ ±iρ/u when u → ±i∞, such as to fulfill (A.3). 10 One can solve for the resolvent directly after casting (A.1) into the form with u ∈ C and where D = e i∂u is the shift operator, D k f (u) = f (u + ik). For u/B = O(1) the large B analysis of (A.6) takes the form of a gradient expansion, since then the derivatives are small, ∂ u ∼ O(1/B), and the shift operator can be expanded, D = 1 + i∂ u + . . . . The equation reduces then to a Riemann-Hilbert problem, which can be solved iteratively, order by order in 1/B. At first, one would like to get rid of the driving term, which gives rise to exponentially small corrections in this "hydrodynamical" regime. This can be done neatly, following [83], by (1) acting on both sides of (A.6) with (D +D −1 ), which explicitly removes the r.h.s. , owing to its 2i-periodicity, and (2) factorizing out (D 1/2 − D −1/2 ). Inverting the latter operator perturbatively ∼ −i∂ −1 u + O(∂ u ) results in a constant of integration which must vanish in our case, by parity. Going along these lines brings the equation which is free of inhomogeneous term, but is nonetheless equivalent to (A.6) to all orders in 1/B for u/B = O(1). This is the same equation as for the O(6) sigma model [83], since in both cases the driving term drops out after performing the (D + D −1 ) projection. The only difference resides in the domain of definition, which is non-compact in our case, u ∈ C, see figure 8.
To leading order, D → 1, and the problem reduces to a well-known singular equation, with the solution Straightforward integration yields the pseudo energy where the constant C 0 = 2ρ log ( 1 2 B) was fixed by demanding that χ 0 vanishes at ±B. This solution is the seed for the next steps and the full iteration can be taken at once by JHEP01(2019)002 means of a suitable ansatz. Drawing inspiration from [83], we write the general solution as an infinite series of terms, with increasing singularities at u = ±B, The ansatz (A.11) solves the problem in generic terms in the domain connected to the boundary conditions at u = ±i∞. To fix its free parameters we must carry out a similar analysis close to the edge of the Fermi sea, at u = B. In particular, χ should be regular at this point, since both the driving term and the kernel are smooth functions. This property is not built into the ansatz (A.11), which shows an accumulation of singular behaviours at u = B. The coefficients must thus be adjusted such that their sum has a smooth discontinuity at u = B.
One can zoom on the neighbourhood of the branch point of the resolvent by using the variable z = B − u. Then, re-expanding (A.11) at large B, with z kept fixed, indicates that the series gets re-organised into a pertubative expansion with each coefficient being a non-trivial function of z, of which we only see the large z tail here. The goal is now to determine all these functions for z = O(1). To this end, it is convenient to introduce the Laplace transform Notice also that the Laplace transform (A.14) determines the energy, which reads, in the low momentum approximation, in response to the large u behaviour of χ . By contour manipulation, we can access to the discontinuity of F (s) across s < 0. First, note that rotating the contour of integration in (A.14) into the upper or lower half plane allows us to continue F (s) to m s < > 0. Then, by taking the difference, we obtain an integral representation for the discontinuity, where b is an arbitrary real number greater than B. The integral is zero for s > 0, since then one can send b → ∞. On the other hand, for s < 0, we can evaluate the discontinuity by taking the inverse Laplace transform of the derivative of the equation (2.12). Using where to make the units apparent we introduced m/4 √ 2 = 1 in front of the driving term. On the r.h.s. we find the Laplace transform of χ on the inner support, which is compact. The latter function is strictly positive for s real, since χ is positive on this interval. Therefore, we readily derive from (A.21) that disc F (s) has simple poles at s = −2πk for k ∈ N and nowhere else. We also observe that disc F (s) is completely determined at s = s n = −π(n − 1 2 ) by the driving term. Since, according to (A.17), disc F (s) ∝ e −Bs at large B, the driving term appears exponentially small compared to the rest. Therefore, up to negligible corrections, its sole effect is to set the discontinuity at s = s 1 = − 1 2 π and one can assume that the discontinuity vanishes at s = s n>1 .
These analytical properties of F (s) together with the requirement that it should admit an expansion in powers of 1/s at large s, for the regularity of χ at u = B, fix the shapes of the functions in (A.17). The same problem was addressed in [83] for the compact model and its general solution was given in the form where Q(s) = n,m 0 Q n,m /(B m+n+1 s n+1 ) represents a general zero mode solution. Comparing the discontinuity of (A.21) at s = −π/2 determines the overall scale, We where the primed quantities refer to the expressions for the O(6) model [83] and where the numerical factors result from the u to θ conversion. In the end, since B behaves like the radius of the target space, the B → −B rule mimics the field theory recipe for connecting the sphere and the hyperboloid perturbatively. Note however that eliminating B leads to formula (A.26) in any case. The difference lies then entirely in the domain of validity of this formula, that is high densities, ρ m, for the sphere and low densities, ρ m, for the hyperboloid.

B Spinning the wheels
In this appendix we discuss a simple class of excited states obtained by adding spin to the operator, with all the derivatives pointing in the same light-like direction, ∂ = n µ ∂ µ with n 2 = 0. We expect these operators to correspond to states with transverse excitations ∼ (X 1 + iX 2 ) M in the sigma model. Below we obtain Bethe ansatz equations for them and, with their help, test the correspondence with the sigma model energy levels E 2d (L, ∆, M ). Adding derivatives to the operator brings an extra layer of difficulty since one must diagonalize a mixing problem. In the one-loop planar gauge theories this mixing results from gluon exchange and leads to the appearance of the sl(2) XXX Hamiltonian. There are no gluons, nor any short-ranged spin chain interactions, in the fishnet theory in this sub-sector of operators. So, here again, the 2d dynamics originates from virtual scalar fields traveling all around the operators. Nevertheless, the eigenstates turn out be the same, to leading order at weak coupling, as those for the Heisenberg spin chain. The point is that the twist that connects the fishnet theory to N = 4 SYM does not affect the Bethe ansatz equations in the sl(2) subsector. Hence, conformal primaries are described by Bethe states, that are linear superpositions of spin waves of derivatives on top of the BMN vacuum. They are generated by the action of the spin-chain creation operator B(v),

JHEP01(2019)002
and are parameterized by the Bethe rapidities v = {v 1 , . . . , v M }. For a periodic spin chain, the rapidities are quantized by the Bethe ansatz equations, where Φ = 0 to leading order at weak coupling. The scaling dimension of the corresponding operator, in the absence of short-ranged interactions, is given to leading order at weak coupling by the Lüscher formula [34,85]. Once reduced to the fishnet theory, one finds where S a * (u, v) is the diagonal part of the scattering matrix between a mirror scalar magnon and a derivative, and where T a =Ṫ a is a rational transfer matrix in the spin 1 2 (a − 1) representation of SU (2). The gas of wheels also induces finite size corrections to the Bethe roots, represented in (B.3) by the additional phase shift Φ = O(g 2L ). The leading-order formula can also be obtained from [34,85] and reads where the trace is taken over the O(4) indices of the S-matrix S 1 * (u, v) between a mirror magnon and a derivative. The all order formula for the scaling dimension takes the same form as before, though the Y functions are now shifted by extra source terms which accommodate for the vdependent part of the integrand in (B.4). The general formula for the phase shift Φ is not known to us, except for its s-wave component. In the latter case the scattering is abelian, S 1 * → S 1 * , and there is no need for the trace in (B.5). By the same token, there are no transfer matrices in (B.4), T 1 =Ṫ 1 = 1. In these circumstances, one can immediately write down the all order formula which is valid up to the contributions of the heavier magnons, a > 1. The only ingredient is the S matrix 11 with Q(u) = M j=1 (u − v j ) the Baxter polynomial, and the same factor should be added to the r.h.s. of the Y 1 TBA equation, which is shifted by Figure 9. The longitudinal excitations ∼ (X 0 + iX d+1 ) ∆ in the thermal bath cross the cylinder and scatter with the transverse excitations ∼ (X 1 + iX 2 ) M in the state. The amplitude for this process is controlled by σ 2 (u + i, v), where the shift of the first argument implements the rotation to the crossed channel, θ → θ + iπ/2.
Note that this shift does not alter the asymptotic behaviour of the Y function, which reads with the anomalous dimension γ = − a dp log (1 + Y a )/(2π). Indeed, the extra source term (B.8) only brings the "M " in this expression, which is the total bare dimension for the derivatives in the state.
We can now proceed with the dualization of the equations. We begin with the equation for the condensate. We just need to act with (1 − K O(6) * ) on the extra source term (B.8). Straightforward algebra gives and thus the dual TBA equation is where the dots stand for the couplings to the higher modes, see equation (4.11). The amplitude (B.11) has a nice interpretation: it is the transmission amplitude for the scattering of two orthogonal complex scalar fields, X 0 + iX d+1 and X 1 + iX 2 , in the O(6) model. Hence, the Y function differs in that the thermal excitations it describes pick an extra contribution owing to their scattering with the derivatives in the state, as shown in figure 9. A similar analysis can be done for the energy of the state which becomes (B.14)

JHEP01(2019)002
Nicely E 2d is the sum of the mechanical energies of the transverse excitations and of the tachyon energy. Furthermore, the transverse excitations contribute positively to the total energy, in agreement with the expectation based on the AdS model. Equation (B.12) captures the pseudo-energy shift induced by the derivatives. The reciprocal action can be found by dualizing equation (B.3). One first rewrites (B.6) as and then plug (B.12) in its r.h.s. The identities show that the spin chain interactions are erased and replaced by the sigma model ones.
The net result is and with the dots representing possible couplings to the Y functions that we discarded. In the following we set the tachyon component to zero, Ψ = 0. Our goal is to compare the BAEs (B.17) with the similar equations obtained in the compact sigma model. It does not hurt being slightly more general and consider the d dimensional version of the problem. In logarithmic form, after going to hyperbolic rapidities, we get to solve the equations They can only be trusted at large L, which means small momenta, P k ∼ ∓ m 2 e −|θ k | ∼ 0. The roots are thus pushed towards θ = ±∞ and split into left and right movers, θ ± k , depending on the sign of their mode number n k . Shifting the roots according to their chirality, where ξ = log ( 1 2 mL) 1, yields the equivalent equations where n ± k are positive integers, regardless of the chirality. 12 In this form it is clear that left and right movers decouple at large L, since [44,45] 12 Note also that there cannot be vanishing mode numbers, n k = 0.

JHEP01(2019)002
where e 2 ∼ 2π/(d log (mL)) is the coupling at distance L. Hence, the perturbative expansion is controlled by the strength ∼ 1/ξ of the left-right interactions. The leading order "scale invariant" or kink solutions are obtained by discarding these interactions. Since the remaining interactions, among particles of same chirality, are of order O(1), the kink solutions can only be found numerically, for generic distributions of mode numbers. Fortunately we will not need to know them explicitly.
What is important is that the massless Bethe ansatz equations (B.21) also apply to the compact model, with the twist that they hold for small length, L 1. Hence, the equations stay the same but ξ = − log (2/mL) is now large and negative. (Also, left and right movers exchange their roles, θ ± k ∼ ∓ log (2/mL) ∼ ∓∞, when going from compact to non-compact.) We could also bring the equations to a more standard form by performing ξ → −ξ, which has the same effect as flipping the sign of the coupling.
Lastly, let us discuss the direct comparison with the sigma model energy levels. A well known trick allows us to evaluate the first two perturbative terms in the energy, without actually solving the equations. Shifting the roots as in (B.20) gives and summing over the r.h.s. of (B.21) yields where we used that log S(θ) = − log S(−θ) and where N ± = M ± j=1 n ± j are the levels for left and right movers. Finally, using (B.22), one concludes that with e 2 ∼ 2π/(d log (mL)). Energy levels on a cylinder can be computed in the CFT limit using the operator-state correspondence. In our case this is applicable in the IR, at large L. The vertex operators of interest take the form where ∂ ± are worldsheet light-cone derivatives. The 2d scaling dimensions δ are known through one loop in the sphere sigma model and converting them to the AdS d+1 space yields [62,86]  fundamental fields in the sigma model. It relates to the small volume "mass gap" which entirely comes from the finite size corrections in the compact case [60,87]. Therefore, it is not surprising to find that it evades the BAEs here as well. Nonetheless, it is interesting to find that the BAEs capture the rest correctly. In the sigma model analysis, this part is the most interesting one and comes from the exchange contribution in the terminology of [62,86]. It would be nice to see if this correspondence holds for generic states-operators, including those carrying non zero ∆.
Finally, let us mention that there is a cleaner way of testing the BAEs, by going to the semi-classical regime of large quantum numbers, e 2 M ± = O(1). A systematic analysis of this regime, like the one carried out in [75] for the sphere, could reveal how much of the Hilbert space of the sigma model is realized in the fishnet theory.

C 3d fishnet graphs
In this appendix we discuss the thermodynamical limit of the triangular fishnet graphs in 3d. These graphs lie at the end point of a suitable twisting of the weakly coupled ABJM theory [14]. The corresponding 3d fishnet theory admits the Lagrangian where every complex scalar field carries U(N ) × U(N ) bi-fundamental indices, and the trace is taken over these indices. Since the gauge fields of the ABJM theory are entirely decoupled at the end of the twisting procedure, the latter group is merely a flavour group. The planar limit corresponds to sending N → ∞ at fixed g 2 and the theory is conformal for any g 2 . In particular, no double-trace interactions are needed here [72]. The simplest single-trace operator is the BMN vacuum operator tr (φ 1 φ * 2 ) L , (C.2) with engineering dimension L. Its anomalous dimension in the planar limit is induced by virtual particles of the field φ 3 that loop around the operator, as displayed in [14].
In the integrability framework, these particles map to mirror magnons and the scaling dimension computes the free energy of the gas of magnons at temperature 1/L. The mechanical energies carried by the mirror magnons take the same form (2.5) as in the 4d theory, but their momentum, at a rapidity u, is half the one in the 4d theory, p = u. In the thermodynamical limit L → ∞ the s-wave magnons, with a = 1, will condense for g 2 > 1/4 and fill a Fermi sea. The scaling dimension scales then thermodynamically and with k(u) as in (2.14). These identities follow directly from the relation between S-matrices for mirror magnons in the 3d and 4d theory [88,89] (C. 6) In the limit where all the energy levels are filled, that is for B = ∞, the solution to the integral equation is obtained by going to Fourier space and reads At this point, the scaling dimension (C.3) and effective chemical potential (C.5) both vanish. The latter condition determines the critical coupling, (C.8) or, numerically, g cr = 0.936 . . .. It agrees with Zamolodchikov's prediction [13].
There is a little twist, compared to what we had in 4d, which relates to the discrete, triangular, nature of the 3d lattice. Namely, the coupling constant of the 3d fishnet theory is not exactly the sigma model energy. The "marginality condition" still takes the form log g 2 = log g 2 cr + u 2 B 2 du 2π φ cr (u)χ(u) , (C.11) JHEP01(2019)002 but φ cr (u) = − 1 2 (1 − K * )k = 1 2i ∂ u E(u + i) + 1 2i ∂ u E(u + 2i) , (C.12) is not identical to the derivative of the dual momentum P (u) = −iE(u + 3i/2). Instead, once written in terms of the hyperbolic rapidity θ, the imaginary shifts in the r.h.s. of (C.12) are seen to implement rotations by π/3 and 2π/3, respectively. These are the angles that characterize the triangular fishnets. However, the distinction is small at low energy, φ cr (u) ∼ √ 3 2 P (u) , (C. 13) after disregarding irrelevant, exponentially small, corrections.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.