Generalized fishnets and exact four-point correlators in chiral CFT4

We study the Feynman graph structure and compute certain exact four-point correlation functions in chiral CFT4 proposed by Ö. Gürdŏgan and one of the authors as a double scaling limit of γ-deformed NUnknowncharacter0xF700from"EuclidMathOne"0x3D=4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{N}\mathrm{Unknown}\ \mathrm{character}\ \left(0\mathrm{x}\mathrm{F}700\right)\ \mathrm{from}\ "\mathrm{Euclid}\ \mathrm{Math}\ \mathrm{One}"\ \left(0\mathrm{x}3\mathrm{D}\right) = 4 $$\end{document} SYM theory. We give full description of bulk behavior of large Feynman graphs: it shows a generalized “dynamical fishnet” structure, with a dynamical exchange of bosonic and Yukawa couplings. We compute certain four-point correlators in the full chiral CFT4, generalizing recent results for a particular one-coupling version of this theory — the bi-scalar “fishnet” CFT. We sum up exactly thecorresponding Feynman diagrams, including both bosonic and fermionic loops, by Bethe-Salpeter method. This provides explicit OPE data for various twist-2 operators with spin, showing a rich analytic structure, both in coordinate and coupling spaces.


Introduction
Quantum conformal field theories in various space-time dimensions attracted recently a considerable attention, due to their phenomenological importance in physics, for subjects ranging from the description of critical phenomena to the fundamental interactions beyond the Standard Model, but also due to their beautiful mathematical structure allowing to get a deep insight into the basic features of Quantum Field Theory and, via AdS/CFT duality, of Quantum Gravity. In spite of the considerable simplifications in the properties of CFTs w.r.t. the massive QFTs, the non-perturbative structure of strongly interacting CFTs in d > 2 dimension is very complicated and in general not very well studied analytically. A considerable progress in this direction has been achieved due to the conformal bootstrap methods [1,2] based on the basic properties of CFTs following from the conformal symmetry, such as crossing symmetry in various channels for the four-point correlation functions. But this approach stays to a great extent "experimental", based on heavy numerical computations rather than on explicit analytic formulation of the final results. A great progress in the understanding of analytic structure of CFTs in d > 2 dimensions has been achieved for various superconformal QFTs, often due to the AdS/CFT correspondence. In a special case -the N = 4 SYM-the analytic study of OPE data was greatly advanced due to the planar integrability of the theory [3]. In particular, the spectral problem -exact, all-loop calculation of anomalous dimensions of local operators -found its ultimate formulation in terms of the Quantum Spectral Curve (QSC) [4,5] -a system of algebraic relations on Baxter-type Q-functions, supplied by analyticity properties and Riemann-Hilbert monodromy conditions (see recent reviews [6,7]).
The integrability appears to persist for a class of 3-parameter γ-deformations of the R-symmetry of N = 4 SYM [8][9][10] if one tunes the so-called double-trace terms, generated by the RG of the model, to their critical, in general complex values [11,12]. The γdeformed N = 4 SYMappears to be a family of non-sypersymmetric and non-unitary four-dimensional CFTs labeled by 't Hooft coupling g and three γ-deformation angles γ j , j = 1, 2, 3. The OPE data of N = 4 SYMhas been studied in numerous papers, using the integrability properties, as well as AdS/CFT correspondence for the strong coupling regime g → ∞, or a direct Feynman graph calculus at weak coupling g → 0. Apart from the spectral problem, an impressive progress has been achieved in a more difficult problem of computation of structure constants and correlation functions [13][14][15][16], as well as of 1/N 2 c corrections [17]. However, the efficient all-loop solution of these problems is still hindered by outstanding technical complexity. We also have to admit that integrability of N = 4 SYM is still a somewhat mysterious phenomenon, not very well understood, especially on the CFT side of this AdS 5 / CFT 4 duality.
In 2015, one of the authors andÖ. Gürdȏgan proposed a family of non-unitary, nonsupersymmetric CFTs [18], based on a special double scaling limit of γ-deformed N = 4 SYM combining weak coupling limit of small 't Hooft coupling, g → 0, and strong imaginary twist, γ j → i∞, with three finite effective couplings ξ j = ge −iγ j /2 . The gauge fields and the gaugino decouple in this limit and one is left with three complex scalars and three complex fermions with certain chiral structure of interactions (see the Lagrangian JHEP06(2019)078 JHEP06(2019)078 undertaken, at least until very recently, to understand deeper the physical properties and the Feynman graph structure of the full χCFT. It is worth noticing that, unlike the biscalar fishnet CFT, the reasons for the integrability of this model remain mysterious.
A few days before the completion of the current paper, a very interesting study of the one loop perturbative properties of this χCFT was undertaken [30], and especially of its two reductions: the bi-scalar fishnet CFT, as well as β-deformed N = 1 supersymmetric case, when all three couplings are equal [19]. The paper explores an interesting subject of study of non-unitary spin chains, having a rich and complicated structure of the spectrum, including the Jordan cells as specific multiplets of states. The Jordan multiplets leading to the logarithmic behavior in non-unitary CFT's [31], are noticed and studied perturbatively in the fishnet CFT [32,33]. The notion of one-loop integrability in χCFT appears to be quite different from the one-loop integrability of its mother theory -the N = 4 SYM [34,35].
The non-unitarity of the studied χCFT represents an obvious drawback from the point of view of the physical interpretations: the presence of complex OPE data violating various basic quantum-mechanical axiomes and usual analyticity constraints. On the other hand, the non-unitary theories are curious objects in themselves, having interesting OPE properties, such as a logarithmic behaviour of certain correlators χCFT is an example of logarithmic CFTs). In addition, they share many basic common features with unitary CFTs and help to understand their general features.
We attempt in this paper to answer some of the questions posed above about the χCFT. First of all, we will give the complete description of the bulk structure of Feynman graphs (far from their boundaries defined by the particular underlying physical quantities). It appears to be much richer than in the fishnet CFT, though much simpler than in the full N = 4 SYMconserving a certain lattice regularity. A pictorious way to describe these graphs is to introduce the regular triangular lattice and the to do all possible Baxter moves of all three types of lines, as shown on figure 2. These lines should represent sequences of bosonic and fermionic propagators and the mixed intersections (where both bosonic and fermionic propagators meet) should be disentangled, in a unique way, into pairs of Yukawa vertices). These configurations should be summed up, so that the collection of such graphs could be called the "dynamical fishnet". The integrability of these graphs, or the sum of them, remains to be proved, though we demonstrate it in this paper in a simpler case of the two-coupling reduction of χCFT (see eq. (2.4)), with two bosonic and one Yukawa coupling.
Then we will compute exactly the 4-point correlation functions of certain short, protected scalar operators, similar those obtained in fishnet CFT [7,12,24]. For that we identify all the graphs contributing these quantities and sum them up using the Bethe-Salpeter approach helped by the conformal invariance. In comparison to the fishnet CFT, the two-coupling dependence of these correlators in the full double-scaled CFT reveals a rich phase structure in the coupling space. We study in detail the related perturbative expansions of these correlators, as well as their strong coupling limits.

JHEP06(2019)078
2 Feynman graphs and correlators of χCFT -the strongly γ-deformed N = 4 SYM theory In this section, we will study the generic structure of planar Feynman graphs and discuss their integrability properties, in the full three-coupling chiral CFT (χCFT) proposed in [18] (see also [19] for more details). This CFT was obtained as a double scaling limit of γ-twisted N = 4 SYMdescribed above. It is defined by the Lagrangian for three complex scalars and three complex fermions transforming in the adjoint representation of SU(N c ): where the sum is taken with respect to all doubly repeated indices, including j = 1, 2, 3, and the interaction part is

(2.2)
We suppressed in the last equation the spinorial indices assuming the scalar product of both fermions in each term. We will refer to this theory as χCFT theory. The double scaling procedure and the derivation of this action from γ-deformed N = 4 SYM can be found in papers [19,36]. In the next sections, we will study the fourpoint functions obtained by point splitting of fields in coinciding points, in the two-point correlation functions of local operators of three types: j (x)] (j = 1, 2, 3), Tr[φ j φ k (x)] (i > j), Since the Lagrangian (2.2) depends on three arbitrary couplings, For some particular values of these couplings, interesting reductions of this χCFT emerge. For example, in the limit ξ 1 → 0, one fermion decouples and we obtain the following action [19] L int = N c Tr ξ 2 3 We will refer to this theory as χ 0 CFT theory. Another interesting case of (2.2) occurs when all three couplings are equal ξ 1 = ξ 2 = ξ 3 = ξ and corresponds to the doubly-scaled β-deformed SYM [37,38]. It has the following interaction Lagrangian [19] (2.2). The graphs of the first line represent the quartic scalar interactions and the ones in the second line are the Yukawa interactions. Tick solid lines and dashed lines represent scalar and fermionic propagators respectively. Colors stand for the various "flavour" of the particles φ i and ψ i : black for i = 1, red for i = 2 and green for i = 3. Arrows symbolize the fixed orientation (chirality) of the vertices and, according to our notation, it points always to the fields with bars or daggers. The second chirality of Yukawa interactions i.e. the one withψ i → ψ i and φ † i → φ i with i = 1, 2, 3, can be represented as the second line of vertices with flipped arrows.
Most of the papers on this relatively young subject were devoted to the abovementioned single coupling reduction of this model: ξ 1 = ξ 2 = 0, ξ 3 ≡ ξ = 0, i.e. the bi-scalar, fishnet CFT defined by the action [18]: Our paper is devoted to generalization of some of these results of [23,24] to the full chiral model -χCFT (2.2) and to its various limits presented above. This represents a step forward, w.r.t. the bi-scalar model (2.6), in understanding the non-perturbative structure of physical quantities of the full N = 4 SYM. We will also describe the general bulk structure of the underlying planar graphs. Indeed, one interesting feature of those models is the drastic simplification of their weak coupling expansions in terms of Feynman diagrams in the planar limit. In general, any diagram of χCFT can be built as a collection of the vertices in figure 1, connected by scalar and fermionic propagators. 1 The arrows indicate the fixed orientation (chirality) of the interactions, i.e. in a propagator it is directed from a field to its hermitian conjugate. An essential feature of (2.2) is the absence of the hermitian conjugate of every interaction vertex, or of the vertices obeying the reality condition. The chirality of this theory makes it non-unitary and plays a crucial role for the underlying conformality and integrability in the 't Hooft limit. In fact, in the absence of the hermitian conjugate vertices, all the graphs which could renormalize the couplings and the mass are non-planar. As a consequence, we will see in section 2.2 that the planar weak coupling expansion of physical quantities w.r.t.

JHEP06(2019)078
interactions (2.2) in χCFT is dominated, at least in the bulk and for high enough perturbative order, by a specific class of planar diagrams having a kind of a lattice structure, much more rigid than the structure of graphs in the original N = 4 SYM. This lattice structure is richer, and more "dynamical" than in the bi-scalar theory where the unique regular square fishnet structure dominates at any order in perturbation theory. In the full χCFT, due to the presence of Yukawa interactions and quartic scalar vertices, there are more planar graphs contributing at each perturbative order, but the chirality still dramatically reduces their number. We can dubb the structure of full χCFT graphs as "dynamical fishnet".

Double-trace interactions and conformal symmetry
The γ-deformed N = 4 SYM theory and its doubly-scaled version are not conformal in a strict sense, not even in the planar limit [36]. Indeed, the renormalization group calculations show [41] that the new, scalar double-trace interactions are generated (2.7) where in our notation j + = j + 1 with the constraint 3 + = 1. The double-trace couplings α k,j generically flow with the scale. They are needed to renormalise the 2-point correlators of the local operators Tr[φ j φ j ], Tr[φ j φ † j + ] and Tr[φ j φ j + ] respectively. For any of these planar correlators only one double-trace term contributes, that is the β-function of each α k,j depends only on couplings {ξ 1 , ξ 2 , ξ 3 } and α k,j itself. Due to permutation symmetry of flavour indices j = 1, 2, 3 in the Lagrangian (2.1), the functions β α k,j show the same symmetry in the coupling dependance, namely thus will drop in what follows the specification of subscript j in double-trace couplings. The double-trace terms (2.7) appear in the theory already at one-loop renormalization and the β-functions associated to the couplings α 2 k are not zero. In γ-deformed N = 4 SYMthe one-loop β-function associated to the double-trace interaction where γ ± k are linear combinations of the deformation parameters γ j of the theory defined in (B.3). Let us turn to the theory (2.2) with the double-trace terms (2.7). In contrast to the bi-scalar theories, where the invariance under exchange allows to identify α 2 and α 3 , the presence of Yukawa interactions in χCFT specifically breaks this symmetry, and operators Tr[φ j φ † j + ] and Tr[φ j φ j + ] show different behaviour.

JHEP06(2019)078
When only one α k coupling is running, the corresponding β-function has the following form where a, b, c are functions of the couplings ξ = {ξ 1 , ξ 2 , ξ 3 }. This quadratic behavior of β as a function of α 2 k was encounter for the first time in [40] as an example of nonsupersymmetric orbifold theories with double-trace interactions and established in [42] for a generic deformed theory in the 't Hooft limit. If the running coupling α k is associated to the double trace interaction TrOTrO † of length-two scalar operators O, the functions a, b and c are related to the normalization coefficient of the two-point function of O, the contribution of the single-traces to the anomalous dimension of O and the coefficient of the induced double-trace terms.
To make the theory conformal at the quantum level, one needs to tune the double-trace couplings to a fixed point. In the original γ-deformed N = 4 SYM, the 't Hooft coupling g 2 is not running, so the critical (conformal) point for double-trace couplings can be computed imposing the vanishing of their β-functions. In the case of a single running coupling, (2.9) has the following fixed points Similarly, the coupling constants ξ i of the theory (2.2) are not running in the 't Hooft limit and one can fine-tune the double-trace couplings α 2 i to critical values in terms of their ξ i dependence, imposing the vanishing of the underlying β-function (2.11) as follows At the two fixed points (2.13), it is possible to write the anomalous dimension γ of the operator O in terms of the discriminant of β α k = 0 [42] (2.14) At the fixed points for all double-trace couplings (2.7) of γ-deformed N = 4 SYM, the theory becomes a genuine non-supersymmetric CFT. This conformal theory appears also to be integrable [9,10,12] and its spectrum of anomalous dimensions can be treated by such a powerful tool as quantum spectral curve (QSC) [4,5,10]. The same statements hold for the double-scaling limit of the 4D χCFT theory (2.2), to which we have to add the doubletrace Lagrangian (2.7). Integrability of the full χCFT is still a conjecture, as it is for the full γ-deformed N = 4 SYM. It was demonstrated explicitly only for the simplest reduction of χCFT -the bi-scalar CFT (2.6), where the fishnet planar graphs have an iterative regular lattice structure [18], shown to be integrable long ago by A.Zamolodchikov [20] (see also [33]). We extended the proof of integrability to a larger, two-coupling sector of χCFT in section 2.4, by methods of conformal SU(2, 2) quantum spin chain. In the case of χCFT we also have good chances to prove full integrability on the level of planar Feynman graphs since, as we show below, these graphs preserve a certain rigid lattice structure.

JHEP06(2019)078
The obvious physical defect of such CFTs is the loss of unitarity. Indeed, as it will be clear with the explicit example below, the discriminant of the equation β α i = 0 is negative, inducing complex values for the fixed points (2.13) and anomalous dimension (2.14). Moreover in the AdS/CFT context, this fact can be interpreted as the presence of true tachyons in the bulk on the string theory side [42].
The one-loop anomalous dimension of the length-two operator Tr[φ j φ j ] in γ-deformed N = 4 SYMat the fixed point is [41] Notice that both the fixed points (2.12) and the anomalous dimensions (2.15) are complex conjugate, as expected. Those relations are actually valid in the full γ-deformed N = 4 SYMtheory, but in the double-scaling limit under analysis it is simple to obtain some predictions for the one-loop β, the associated critical points and the anomalous dimensions. In particular we have In section 4.3 and section 6.2 we will verify these results computing the exact spectrum of the operator Tr[φ j φ j ] with the Bethe-Salpeter method, and the first order of the fixed point α 2 1 using Feynman diagrams. Non-unitary CFTs are usually logarithmic [31], i.e. with an interesting, logarithmic behavior of certain correlators. The γ-deformed N = 4 SYM and its double-scaled version -the χCFT (2.2) (and its reductions mentioned above) are not exceptions: they show the same logaritmic properties due to the non-hermiticity of their dilatation operators [32,33].

The bulk structure of large planar graphs
Let us try to describe the general structure of an arbitrarily big Feynman graph in the bulk, far from the boundaries. The generic picture is illustrated on figure 2. The theory (2.2) contains 3 complex scalars φ i and 3 complex fermions ψ i labelled by i = 1, 2, 3. We chose to represent scalar propagators with thick solid lines and fermionic propagators with dashed lines (see figure 1), while the label denoting their U(1) ⊗3 flavour (see appendix 4) is mapped into colours: (1, 2, 3) ≡ (black, red, green). In figure 2, coloured dotted lines in a particular direction represent a generic propagator, both scalar or fermionic. In this framework, a set of parallel lines represents any combination of fermionic and scalar propagators of a given flavour.
This system of three dotted lines forms a lattice which combines the features of both regularity and irregularity. Any such lattice can be obtained from the regular triangular lattice (or a more general Kagomé lattice) by arbitrary Baxter moves of all lines: displacements in the direction orthogonal to the line, i.e. conserving its direction.
The links of the resulting lattice are propagators while nodes are quartic effective interactions. These interactions are of three kinds, depending on which lines are crossing and which propagators enter the corresponding crossing (effective vertex). They can represent JHEP06(2019)078 Figure 2. General "dynamical" fishnet bulk structure of a planar graph for 3-coupling chiral CFT (2.2). Dotted lines represent scalar or fermionic propagators (the rules for the choice of propagators will be explained below and demonstrated in figure 3) The colors and directions of the lines stand for the three "flavours" of the particles i = 1, 2, 3 with the same notation as we used in figure 1. The intersections correspond to six different effective vertices that can be written in terms of the usual ones following the map given in table 1. a set of φ 4 or various Yukawa vertices, according to the rules listed in table 1. Indeed in this framework, a quartic vertex involving fermions can be though of as a couple of Yukawa vertices, or similarly, as a split quartic vertex in which we have added a propagator in the remaining direction, according to the rules in table 1. The quartic interaction can involve four scalars, four fermions or two of each. Moreover, we chose the directions of the arrows to be consistent with the Feynman rules in figure 1. Depending on the orientation of the mixed interactions we will refer to them as crossing or scattering interactions as in table 1.
Given three sets of parallel lines crossing each other with quartic interactions, the resulting irregular lattice is formed by a finite set of convex polygons. The smallest possible n-gon is a triangle and the largest one is a hexagon. Those convex polygons can be constructed locally by the abovementioned moves of lines in two or three different directions: • 2 directions (colors): we can discard the lines in one of the directions. The local interaction of lines with only two directions (colors) forms a square lattice as in [18,20]. Since we are considering three colors, we can have three different squares depending on their directions.
• 3 directions (colors): in this case there are more possibilities to build convex polygons. Indeed let's start with the crossing of three lines with three different directions. Locally, they form a triangle that can have two different orientations. Adding another line, parallel to one of the previous three, and cutting the triangle, we will end up with a square. Since we can add a line of any color and there are two possible triangle orientations, we can draw 6 different squares. Iterating this cutting procedure by adding one and two lines we obtain pentagons and the hexagon.
In the following  Table 1. Substitution rules for the effective vertices appearing in the fishnet bulk structure of figure 2 in terms of the Feynman rules of figure 1. Any effective vertices is associated with a combination of the coupling constants ξ i with i = 1, 2, 3 of order ξ 2 . same polygon can appear in the graph.

JHEP06(2019)078
n-gon Multiplicity 2 9 6 1 It follows that for a given set of lines, the resulting lattice can be seen as a tiling of the plane with 18 different tiles. The structure of the fishnet bulk is very rich, indeed once the topology of the lattice is defined as in figure 2, some information is lost, as any quartic dotted-vertex can be associated to six different physical vertices, as listed in table 1. The number of possible Feynman diagrams N d which can be associated to a given close n-gon, defined by n quartic dotted-vertices, can be computed considering first all possible combinations of fermionic and scalar propagators for the edges of the polygon and then cancel out those vertices which does not fit in any configuration. After this tedious combinatorics we obtain the following table   n-gon  N d  28 82 244 730 This result can be written in the following compact formula Now we can estimate the number of Feynman diagrams for a given topology of the dottedfishnet bulk. This number has the sum of all the N d 's for all the polygons as an upper bound and we can estimate its order of magnitude. Then the number of possible Feynman diagrams for the topology of the fishnet bulk given in figure 2 is around 1.5×10 4 . Moreover, since any vertex is associated with a combination of the couplings ξ i with i = 1, 2, 3 of order 2, we know that the diagram in figure 2 is of order ξ 234 . One of those configurations is represented in figure 3.

JHEP06(2019)078
, and the graph is of order ξ 42 . As it results from table 1, each effective vertex can be replaced in a univocal way in terms of structure made of real vertices.

Single-trace correlation functions
We can realize the above mentioned bulk graphs (with fixed coordinates of external legs) as a single-trace operator of the form: i.e. each χ(x) under the trace is one of 18 fields of the χCFT model (2.1)-(2.2). Of course (2.18) must have zero overall R-charge, to have a non-zero answer. This implies a condition on the elementary fields under trace, namely if we define n j and m j as the differences between the number of φ j , respectively ψ j and the conjugated fields, the mentioned condition reads To describe the Feynman graph content of this quantity, let us remind that a similar singletrace correlator in bi-scalar fishnet CFT [25,26], consisting only of scalar fields, was given by a single fishnet graph of the disc topology where the disc was cut out across the edges of a regular square lattice. The ends of the cut edges represented external fixed coordinates and the integrals were taken over all vertices inside the disc. Similarly, for each of the quantities (2.18) there exist a collection of graphs of the disc shape cut out of the lattice of the type drawn on figure 3. The types of external legs -the cut edges along the boundary -define the species of fields from the set χ following in the same order under the trace in (2.18). We present an example in figure 4, where the disc is drawn on the concrete realization of the lattice as given in figure 3. A big difference w.r.t. the bi-scalar singletrace correlators is that in the full χCFT such a quantity is defined by the sum of all graphs with the same order of fields on the boundary (same sequence of external legs) which are related to each other by the orthogonal moves of three types of parallel lines described in the JHEP06(2019)078 previous subsection (as example, see figure 5 (right)). Furthermore, even at fixed topology, one can change the interaction vertices inside the graph, namely switching some dashed (fermionic) lines to solid (scalar) lines and vice-versa (figure 5 (left)). This corresponds to different realizations of a disc segment of the dotted-lattice in figure 2 with boundary conditions fixed by the external legs. The number of possible graphs can be estimated by considerations of the previous subsection. This single-trace correlator can be used to define the scattering amplitudes via Lehmann-Symanzik-Zimmerman procedure, by going to the dual momentum space and taking on-shell external momenta, in the spirit of the papers [25,26]. It is worth noticing that not all the planar single-trace correlators are obtained out of this procedure. Indeed certain external states can be cut out only drawing a circle on the actual Feynman graph (see table 1) where all propagators are explicitly drawn. Moreover, for a given correlation function, there are lower order graphs in the coupling which cannot be cut out of the planar lattice, but from two or more sheets of such lattice as explained in [25]. It would be interesting to show the Yangian invariance of these single-trace correlators, in the same spirit as it was done in [25,26] for bi-scalar case. Namely, to define the monodromy ("lasso") around the boundary for which each of these graphs, or sum of all graphs, is an eigenfunction. This is one of the ways to show the integrability of the full χCFT. In this paper we will limit ourselves by the proof of integrability of the model (2.4) which will be given in the next subsection.

Integrability of Wheel graphs in χCFT
A statement of integrability, milder than the lattice integrability of the bulk of large planar graphs, can be made for the scaling dimension of Tr[φ L j ] operators at any L. These operators, protected in the original N = 4 SYM due to supersymmetry, are described in the planar limit of bi-scalar theory by a perturbative expansion in globe-like fishnet graphs [18] with an integrable square-lattice bulk [20]. These graphs can be built up by the action of JHEP06(2019)078 an integral "graph-building" kernelĤ (2.21) It represents one of the conserved charges generated by the transfer matrix of the integrable quantum SU(2, 2) spin chain of L sites in the scalar (∆, J 1 , J 2 ) = (1, 0, 0) representation [33]. Similarly, in the two-coupling version (2.4) of χCFT the perturbative expansion can be described by graphs which, in spite of more complicated structure (see figure 6), can be still built by integrals of motion of the conformal spin chain. Namely, every planar graph in the ξ j expansion is a certain permutation of multiple action of operatorsĤ (L) B and H (L) F , where the latter operator is responsible for fermionic loops contribution. As we will see, the order in the permutation doesn't matter, since any fermionic loop can be moved through scalar wrappings, due to their commutativity, and this fact lays at the basis of integrability of these graphs. The action ofĤ and it builds up an integrable "brick-wall" domain [26]. Its commutation withĤ B can be proven directly by star-triangle relation (C.2), as shown in figure 7.
In order to show thatĤ F is a conserved charge of the conformal scalar spin chain, we should prove its commutation with the transfer matrix at any value of the spectral JHEP06(2019)078 x 1 x 2 x 3 y 1 y 2 y 3 x 1 x 2 x 3 y 1 y 2 y 3 y 1 y 2 y 3 x 1 x 2 x 3 parameter u, (2.23) For this purpose, we rewrite the kernel integrating out z k variables and we recall the definition of T (L) (u) where R ij (u) is the R-operator of the scalar conformal chain. It satisfies the Yang-Baxter equation [43] R Then operator (2.21) coincides with 4 −2L T (L) in the limit u → −1, as pointed out in [33], since the first propagator under the integral in (2.24) disappears and the last one effectively becomes a δ-function. Now we introduce the transfer matrix for the brick-wall domain 2 T (L) and we check, similarly to the above scalar case, that lim u→−1 T  (2.28). Both sides can be transformed in the hexagonal object in the middle. First the triangle is opened into a star integral using (C.1) (left side) or (C.2) (right side). Doing so, each of the three black dots will become the end of only three lines. Then integration can be performed by (C.1), (C.2) and leads to the hexagon. which will be done by means of a Yang-Baxter type relatioñ and Ω then iterating (2.28) we can write and we finally trace over space L 2 (x 0 ) ⊗ L 2 (x 0 ) and over spinorial indices getting which is equivalent to (2.27). Our derivation straightforwardly shows that from the point of view of integrability the regular square lattice and the brick-wall lattice built by Yukawa vertices can be combined into the same integrable structure and form a mixed lattice. This concludes the demonstration of integrability of the two-coupling model (2.4). The proof of integrability of the full χCFT (2.2) is a more tricky exercise and we leave it for the future.

Bethe-Salpeter equation for four-point correlators and conformal data
In this and the next sections of this paper, we will exploit conformal symmetry and the Bethe-Salpeter method to obtain the exact 4-point correlations functions

JHEP06(2019)078
where the operators O i are protected operators in the double-scaled γ-deformed N = 4 SYM theory, the χCFT. Then we will extract from it the OPE data, anomalous dimensions and structure constants, for length-2 unprotected operators exchanged in the s-channel of (3.1). In the current section, we present the generalities of conformal Bethe-Salpeter approach, generalizing the one applied in [12,21,24] to the bi-scalar fishnet CFT, to sum up the Feynman graphs for these quantities in χCFT.
At the fixed point (2.13) and in the planar limit, the correlation function (3.1) is a finite function of the couplings ξ i with i = 1, 2, 3. The correlation functions can be remarkably written as a geometric sum of primitive divergencies in the perturbative expansion. For this reason, we will study those diagrams with the help of the Bethe-Salpeter equation. In the following we will review the Bethe-Salpeter method pointing out how to extract the spectrum and the OPE data from the four-point functions (3.1). In section 2.2, we presented the bulk fishnet structure of large planar diagrams in the general double-scaled γ deformed N = 4 SYM theory. In this section we will focus on the correlation functions defined by (3.1) for matrix (untraced) operators with bare dimensions ∆ O 1 and ∆ O 2 . Since to preserve the renormalizability of the theory we have to supplement it with double-trace counter-terms (2.7), diagrams in the perturbative expansion of (3.1) will take the following chain structure where the black dots are insertions of the double-trace operator 3 and the links of the chain are periodically repeating configurations of propagators (a special case of the topologies presented in section 2.2) generated by the kernel of integral operators. We will refer to this set of operators as Hamiltonian graph-building operatorsĤ i . In the family of theories we are considering,Ĥ i Hamiltonians can be of three different kind: the double trace operator V, the bosonic operatorĤ B and the fermionic operatorĤ F . The operatorsV,Ĥ B and H F separately produce divergent integrals. However, at the fixed point, their combination is finite due to conformal symmetry (see section 6). These integral operators commute among themselves and they are diagonalized by the same basis of conformal trianglesthe 3-point correlators of the protected operators with un protected operator with spin, described below. The correlation function (3.1) can be written in general as a geometric series of a linear combination of the Hamiltonian graph-building operators as followŝ

2)
3 Such insertions should always split a graphs, and its color structure, into two disconnected pieces.

JHEP06(2019)078
where c B = 1/(4π 2 ) is the normalization factor of the free scalar propagator, χ V , χ B an χ F are combinations of the couplings α i and ξ i with i = 1, 2, 3, which will be introduced later (see section 4). The spacetime dimension D in this paper is always taken to be D = 4. 4 The correlation function (3.1) can be obtained from the operatorĜ O 1 O 2 as follows where the Hamiltonian operators are represented by the corresponding integration kernels such that (3.4) In order to compute the correlators G O 1 O 2 , given the set of Hamiltonian graph-building operatorsĤ i , we need to compute their eigenvalues and decomposeĜ O 1 O 2 over a complete basis of their eigenfunctions.
To compute the eigenvalues ofĤ i , we can use the fact that these integral operators transform covariantly with respect to the (1, 0, 0) ⊗ (1, 0, 0) conformal spin chain generators. 5 This property completely fixes their eigenstates to be the conformal triangle Φ ∆,S,x 0 (x 1 , x 2 ), the three-point function of two (scalar) operators in x 1 and where p = ∆−S 2 and n µ an auxiliary light-cone vector. In the case S = 0, the conformal triangle is composed by simply three scalar propagators that we can graphically represents as follows Finally, given the eigenstate (3.5), we can compute the spectrum of the Hamiltonian oper-atorsĤ i as follows 4 It is possible to generalize the bi-scalar fishnet theory to any integer dimension D, as in [21], at the cost of losing locality. It is not evident that such a generalization is possible for the full χCFT. 5 In particular, defining the inversion . We can check that for every integral operator:

JHEP06(2019)078
where h i∆,S is the eigenvalue. More specifically, given the Hamiltonians operators defined in (3.2), we have In section 4 and section 5, we will verify that (3.5) diagonalizes these Hamiltonians and perform a direct computation of the eigenvalues. The scaling dimension appearing in (3.5) is defined as [44] ∆ = 2 + 2iν , with ν a non-negative real number. For such values of ∆, the state Φ ∆,S,x 0 belongs to the principal series of type-I irreducible representations (∆, S, 0) of the conformal group labelled by ν and the discrete compact spin S and satisfies the orthogonality condition [44,45] where the 4-dimensional coefficients c 1 and c 2 are given by . (3.13) The eigenfunction Φ ∆,S,x 0 forms an orthonormal basis for ν ≥ 0 implying the following representation for the identity (3.14) that, together with the definition (3.7), leads to the diagonalized representation

JHEP06(2019)078
Plugging the representation (3.15) for the graph-building operators H i into (3.2), we obtain the representation of the 4-point function in terms of their eigenvalues h i∆,S The integral over the auxiliary point x 0 can be expressed in terms of the four-dimensional conformal blocks g ∆,S [44,46,47] where the cross-ratios are u = zz = x 2 12 ) and we recall from [46] that where we defined Notice that we extended the integral over ν on the full real axis with the change of variable ν → −ν in the second term of (3.17). This is allowed by the symmetry of eigenvalues appearing in the spectral equation and can be interpreted as the fact that, for a given spin S, states with dimension ∆ and 4−∆ belong to a unitary equivalent representation of the conformal group. This symmetry is indeed satisfied for every studied case, (4.14), (4.25) and (5.31). Before studying the integral in (3.20), we want to focus on the role of the doubletrace Hamiltonian and its eigenvalues in the perturbative and Bethe-Salpeter approaches. To find the correlation function (3.1), we have to sum up diagrams of the kind shown at the beginning of this section. These diagrams contain an involved scalar and fermionic JHEP06(2019)078 structure generated by the operatorsĤ B andĤ F interspersed with the contributions of the double-trace vertices introduced in (2.7). Since in general the integrals over the positions of the single-trace vertices develop ultraviolet (UV) divergencies at short distances, one needs the double-trace interactions to produce other UV divergent contributions which cancels against them. Therefore, the weak coupling expansion of the four-point correlation function remains UV finite at any order as expected for protected O 1 and O 2 .
In the context of the Bethe-Salpeter equation the story is slightly different. Indeed consider the Hamiltonian operatorV associated to the double-trace kernel defined as follows where Φ(y 1 , y 2 ) is a test function. We have to compute its spectrum by means of (3.8) that, when applied to (3.5), reads First of all, due to the form of the eigenvalue, the double-trace term can affect only the contribution to the sum in (3.20) with spin S = 0. Then we expect that the contribution to (3.20) given by the Hamiltonian operatorsĤ B andĤ F are well-defined for S = 0 but in principle we have to take into account the double-trace term for S = 0.
Since we want to write G O 1 O 2 in the standard OPE form, we will consider the limit in which two of the external points are approaching, i.e. |x 12 | → 0 (or u → 0 and v → 1). Since the conformal block scales as u p (1 − v) S decaying exponentially for Re(iν) → ∞, one can close the contour in the integral over ν in lower-half plane and then compute it by residues.
At short distances, the eigenstate (3.5) scales as Φ ∆,S, and thus it vanishes in the lower-half plane (which is true in our case, since ∆ O 1 = ∆ O 2 = 1 and (∆) > 2). In this case, the bosonic and fermionic operators do not develop UV divergencies (one can verify it in the two special cases that we study in detail in section 4 and section 5). Moreover, given the definition (3.22) and the formula (3.23), the doubletrace operatorV annihilates Φ ∆,S,x 0 with Im(ν) < 0 for any S and therefore, it does not contribute. With this argument, we are able to neglect the double-trace contributions when we compute the four-point function G O 1 O 2 with the Bethe-Salpeter method. Then we can rewrite (3.20) as follows where C − is the close path in the lower-half plane.
In order to compute the integral over ν in (3.24) with residues, we have to identify the poles of the integrand. The physical poles are given by the zeros of denominator under the integral, i.e. by solutions of the equation

JHEP06(2019)078
We will refer to (3.25) as spectral equation: indeed given the eigenvalues h i∆,S and the constants χ i , solving the equation we will obtain the scaling dimensions ∆ as functions of the couplings ξ i with i = 1, 2, 3 and the spin S. In the integrand of (3.24), two series of spurious poles are generated by the measure c 2 and the conformal block g ∆,S . In appendix E we will prove that the contribution of those poles cancel under the condition which happens to be satisfied. Finally G O 1 O 2 is given by the sum of only the residues at the physical pole (3.25). Then we can rewrite the correlation function in the standard form of a conformal partial wave expansion as follows with the OPE coefficients C ∆,S defined as The sum over ∆ in (3.27) runs over the solutions of the spectral equation for scaling dimensions of exchanges operators with Re(∆) > 2. 6 In the following sections, we will focus on the computation of the point-split four-point correlation functions of the operators introduced in (2.3), establishing the Hamiltonian op-eratorsĤ i and the constants χ i appearing in (3.1) from their Feynman diagram expansion. We closely follow in our analysis the logic of [24], but in contrast to this paper which treats the bi-scalar fishnet CFT, we have to introduce new types of diagrams into the Bethe-Salpeter procedure, reflecting a richer structure of the full three-coupling χCFT. To write the correlation function (3.1) in the standard OPE representation requires, as the only dynamical input, the knowledge of eigenvalues h i∆,S of the Hamiltonian operators. We will diagonalizeĤ i to extract the conformal data, i.e. the scaling dimensions of the exchanged operators and the OPE coefficients. In what follows we consider only single scalar fields as protected external operators and then we should set ∆ O 1 = ∆ O 2 = 1. Since the four-point correlator constructed from the second operator of (2.3) is trivial (see section 6.1), in the following two sections we will analyze the remaining two.

Exact four-point correlations function for
In this section we consider the four-point correlators associated to the first operator of (2.3), Since the computation of the correlators is the same for any j, we will consider the case j = 1 and then the four-point function we want to study takes the following form This correlation function was extensively studied in [12] in the simplest case of the family of theories we are inspecting, i.e. the bi-scalar theory (2.6).
In the planar limit N c → ∞, once chosen j = 1, the weak coupling expansion of (4.1) in terms of Feynman diagrams is given by a combination of the following bosonic vertices and the following Yukawa vertices In the following we will study the correlation function (4.1) with the Bethe-Salpeter method.

The Bethe-Salpeter method for the correlator
The perturbative expansion of (4.1) can be written in the following form where G ( ) φ 1 φ 1 at any perturbative order contains contributions from the bosonic and fermionic integrals with different coupling dependencies. In figure 9, we present an example of an arbitrary Feynman diagram contributing to G ( ) φ 1 φ 1 . The black dots represents insertions of the double-trace vertex in the last line of (4.2) that in the Bethe-Salpeter picture are associated with the operatorV defined in (3.22). Then it is straightforward to fix the normalization of its coupling constant in (3.2) as follows In section 3, we discussed the role of the double-trace terms in the computation of the fourpoint function, discovering that they are not contributing to the spectral equation. Then, similarly to observations of [12,24], as far as we consider the perturbative expansion (4.4) JHEP06(2019)078 in the point-splitting x 1 = x 2 and x 3 = x 4 , we need only to sum over the single trace contributions, namely the diagrams inside the chain link of figure 9. In section 6 we will present in detail how the relation between single-and double-trace terms is crucial for the setting of the fixed point (2.13).
The first two orders of the perturbative expansion are given by the diagrams represented in figure 10 and they can be written as follows where each scalar propagator brings in the factor c B /x 2 ij and each fermionic propagator the factor Since the fermionic propagator can also be written as . These functions can be expressed in terms of a combination of the Hamiltonian graph-building operatorŝ H i . Indeed defining the following kernels represented in figure 11, we can rewrite (4.6) as follows The kernels (4.7) transform covariantly under conformal transformations, 7 then the corresponding Hamiltonian integral operators commute with the generators of the conformal group.
JHEP06(2019)078 When carrying on the perturbative expansion, it becomes clear that an arbitrary diagram at order is given bŷ Then the correlator (4.4) can be presented in the following operatorial form Comparing it with the definition (3.2) we fix the values of the constants χ i (in this caseV is not contributing)

Eigenvalues of the Hamiltonian graph-building operators
Writing the four-point correlation function in the standard OPE form, as presented in detail in section 3, involves the computation of the spectrum of the graph-building operators (4.7). The eigenstate that diagonalize the Hamiltonians is defined in (3.5) for ∆ O 1 = ∆ O 2 = 1 and the eigenvalues are defined by means of equations (3.9) and (3.10). Substituting in the latter the kernels (4.7) and using the definition (3.7), we will end up with a set of integrals that can be computed with the help of the star-triangle relations presented in appendix C (also known as uniqueness method ). The fact that all the integrals that we have to compute can be computed by means of the star-triangle relations is a consequence of the underlying conformal symmetry.
Bosonic eigenvalue. The bosonic eigenvalue h B ∆,S is defined in (3.9). Using the bosonic Hamiltonian (4.7), this relation can be written in the following integral form (4.12)

JHEP06(2019)078
In the case of S = 0, the function Φ ∆,S,x 0 reduces to (3.6) and the computation is straightforward. Indeed, one needs to apply the star-triangle relations two times as follows 8 The eigenvalue at S = 0 can be computed in the same way, using the generalization of startriangle relation to any-spin case (C.3) derived in [50]. The computation can be otherwise done in a more tedious and explicit way as presented in detail in [24]. The result reads [12] h The eigenvalue is invariant under ∆ → 4 − ∆, as expected from (3.21).
Fermionic eigenvalue. The fermionic eigenvalue h F ∆,S is defined in (3.10). This is a new object, absent in the similar correlator of bi-scalar model treated in [12]. First of all we can simplify the fermionic Hamiltonian in (4.7) integrating the primed variables by means of the Yukawa star-triangle identity (C.2) as follows (red lines are spin-1/2 fermionic propagators) where the computation and figures are made with the STR package (see footnote 8). We obtain the following kernel Using the formula for the trace of four σ-matrices (A.11) and simplifying the scalar products by means of (A.4), we can rewrite the fermionic hamiltonian in the following form where we used the symmetry ) of the bosonic hamiltonian studied in the previous paragraph, andH F is defined bỹ (4.17)

JHEP06(2019)078
Then the fermionic eigenvalue h F ∆,S consists of the bosonic eigenvalue (4.14) and the eigenvalue ofH F defined as follows Let's focus on the relation (4.18). It can be written in the following integral form In order to simplify the computation, we consider the limit in which x 0 → ∞ on both sides of (4.20). In this limit the eigenvalueh F ∆,S is given by the following integral where p = ∆−S 2 and we put x 2 12 = (nx 12 ) = 1 for convenience. Notice that the integrand is antisymmetric in the exchange y 1 ↔ y 2 for odd S, then the eigenvalueh F ∆,S is non-zero only for even S.
In the S = 0 case, the integral (4.21) is known as a massless two-loop self-energy Feynman integral, or kite. Its value is known for any power of the propagator 1/y 2 12 in terms of an hypergeometric function [51], theñ where ∆ = 2 + 2iν. Expanding (4.22) around ν = 0, one can notice that the cotangent cancels all the odd terms of the hypergeometric functions. The analytic properties of (4.22) are more clear when writing it in the following equivalent form where ψ (1) (x) = dψ(x)/dx and ψ(x) is the digamma function. When S = 0, we can appeal to a similar computation made in [24]. In fact, the same integral of (4.21) appears in the study of the spectrum of the graph-building operator associated to the 2-magnon correlation function. The 2-magnon Hamiltonian is H 2-magnon = x 2 34 /x 2 12H F but, when applied to the eigenstate Φ ∆,S,x 0 , that has in this case ∆ O 1 = ∆ O 2 = 2, it leads to an eigenvalue with the same integral representation as (4.21). The strategy to compute the eigenvalue is to write the following recursion relation for the integrals (4.24) Solving the recurrence with the eigenvalueh F ∆,S=0 , given by (4.23), as initial condition, we obtainh We can conclude that the eigenvalue (4.19) is manifestly invariant under ∆ → 4 − ∆, as expected from (3.21).

Spectrum of exchanged operators of
In this section we will use the eigenvalues (4.14) and (4.19) to compute the scaling dimensions of the operators contributing to the correlation function (3.27) for The spectrum of the exchanged operators is defined by the solutions of the equation for the physical poles (3.25). Substituting in (3.25) the definition of bosonic and fermionic eigenvalues (4.14) and (4.19) and the constants χ i computed in (4.11), we can rearrange the spectral equation in the following form where we defined the new couplings and Plugging (4.14) and (4.25) into (4.26), we obtain the following spectral equation with the additional constraint Im ν < 0. This equation can be studied perturbatively, for each individual anomalous dimension, expanding in ν around the value ν 0 corresponding to a bare dimension ∆ 0 = 2 + 2iν 0 at weak coupling, and in 1/ν at strong coupling.
Weak coupling expansion. The small coupling limit suggests that the equation has solutions with bare dimensions 2 + S and 4 + S, in analogy with the same quantity in the bi-scalar theory [12,24]. There are six such solutions, but only half satisfies the physical requirement Re ∆ ≥ 2: one of them corresponds to the scaling dimensions of exchanged operator with bare dimension ∆ 0 = 2+ S and two -to the scaling dimensions of operators with bare dimensions ∆ 0 = 4 + S. The remaining three solutions are related to the first ones b the transformation ∆ → 4 − ∆ and describe shadow operators, with Re ∆ < 2. In addition to that, there is an infinite series of physical solutions around the bare dimensions ∆ 0 = t + S with t = 6, 8, . . ., due to the non algebraic eigenvalue h F ∆,S , similarly to the two-magnon case studied in [24]. For each value of the twist t there are two solutions; they describe the exchange of an infinite tower of local primary operators in the OPE of 4.5.
Writing ν as a function of the two couplings (4.27) and expanding around the physical pole ν = −iS/2 at weak coupling κ, ω → 0, we obtain the following expansion for the twist-two operator ∆ (4) = 4 + S + 4κ 2 (S + 1)(S + 2) + (S + 2)ω 4 − 8κ 4 (S + 1)(S + 2) 2 + S+2 −H  where H (2) k are Harmonic numbers of order 2. Remarkably, the expressions in square brackets in (4.29), (4.30) are in fact rational numbers. In both cases, we present only the first few terms since the following ones are cumbersome. We notice that the weak coupling expansions of ∆ (4) , ∆ (4 ) are divergent but, as it will be pointed out later in the analysis of section 4.5, the sum of the two corresponding OPE contributions has a well defined expansion. 9 Similar considerations can be made for the solutions at higher twist t = ∆ 0 − S = 6, 8, . . .
The twist-2 solution corresponds to the operator namely the traceless symmetric S-tensor obtained by insertion of light cone derivatives ∂ = n µ ∂ µ , n 2 = 0 into the operator Tr[φ 2 1 ]. At twist-4 the matter content of the theory allows to find several S-tensor operators satisfying the condition ∆ 0 − S = 4 and having the right U(1) ⊗3 quantum numbers (e.g. for i = 1, j = 2: (2, 0, 0)). Twist-4 operators start mixing with each other. We perform an introductory analysis of this phenomenon for the simple scalar case S = 0 in appendix F.1. At this stage the log-CFT effects [31] due to chiral interaction vertices in (2.2) show up. The analysis suggests the presence at twist-4 of only two non-protected physical operators, which should be identified with the two solutions ∆ (4) and ∆ (4 ) at S = 0, in contrast to the bi-scalar fishnet CFT where only one type of twist-4 operators appears [12].Similar considerations apply to the higher twist operators ∆ (t) ± . Indeed also for value of twist t > 4 it is possible to find several S-tensor primary operators with the correct set of Cartan's charges. The detailed study of these operators and their mixing would be an interesting insight in the structure of operator algebra of χCFT. We will restrict from here on most of our analysis to solution of twist two and four, whose contribution to the OPE expansion appears to be enough for complete description of the first non-trivial order of the weak coupling expansion, confirmed by direct computations in terms of Feynman diagrams.

JHEP06(2019)078
Recalling the definition (4.27), the weak coupling expansion (4.29) for the twist-two operator goes in powers of ξ 4 of the original couplings which is exactly the expected behavior considering that the perturbative expansion in figure 10 alternates bosonic and fermionic wheels attached to the diagrams with two quartic or four Yukawa single-trace vertices. On the contrary, the weak coupling expansions (4.30) for the twist-four operators goes in power of ξ 2 of the original coupling. This fact can be understood looking at the expansion of (4.28) around the physical pole located at ν = −i(S + 2)/2. Indeed this expansion starts from κ 4 /(ν + i(S + 2)/2) 2 and as a consequence the four-point correlation function (4.10) is convergent when ν → −i(S + 2)/2 if κ is finite while it produces a divergence when we consider the weak coupling limit κ, ω → 0 such that G φ 1 φ 1 ∼ ±κ 2 .
The zero-spin case presents some peculiar behaviours. Indeed, expanding (4.28) for S = 0 around the physical poles ν = 0, −i at weak coupling, we obtain the following expansions for the solutions of (4.26) where the one-loop order of the scaling dimension ∆ (2) S=0 is in agreement with the prediction (2.16). This twist-2 solution is the scaling dimension of the operator Tr[φ 1 φ † 2 ], while the two solutions of twist-4 arise from the operatorial mixing in a similar way as to S = 0 case analysed in appendix F.1.
Notice that the weak-coupling expansion of the twist-two operator is drastically different as compared to the S = 0 case, indeed it goes in powers of ξ 2 . The same behavior was noticed in [24] and the reason is similar to the one explained above. We observe that, expanding around the physical pole ν = 0, the spectral equation (4.28) goes as ω 4 /ν 2 . Then, when ν → 0, the correlation function (4.10) is convergent if µ is finite, but it produces a square-root divergence when we expand at weak coupling, as in the previous case. The fact that the weak-coupling and S → 0 limits are not commutative is related to this divergence.
The divergence in the expansion of the scaling dimension of the twist-two operator is not a surprise. In fact, as noticed also in some different contexts in [12], in order to write the correlation function in the OPE form as in (3.27), we assumed that in the integral (3.24) no physical poles are located on the real ν-axes. However the poles that at weak coupling and when S = 0 are situated at ν = ∓iS/2 pinch the integration contour at the origin when S = 0, thus producing a divergence. Hence, the contribution of the double-traces is needed in this case to produce a non-vanishing term that cancels this divergence at weak coupling. Again, we stress that at finite couplings the solutions of (4.28) are well-defined even at zero spin.

JHEP06(2019)078
Strong coupling expansion. At strong coupling, κ , ω → ∞, we consider the four solutions of eq. (4.28) of lowest twist. The solutions are related to the physical poles of the spectral equation located at ν = e iπ k 2 4 √ ω 4 + 2κ 4 + . . . with k = 0, 1, 2, 3 but only two of them satisfy the condition Im ν < 0, the remaining solutions being associated to the shadow operators. However we stress that we are neglecting all the infinite non-algebraic solutions of higher twist, purely generated by h F ∆,S . Then we have Notice that, if all couplings scale as ξ j ∼ ξ 1, the strong coupling expansion (4.36) is growing linearly with ξ. This becomes clear if one expands the eigenvalues appearing in (4.26). Indeed both of them decay at large ν as h B ∆,S ,h F ∆,S ∼ 1/ν 4 then, since in the spectral equation the couplings appear in power of ξ 4 , it is evident that the expansion will contain terms linear in ξ. The S → 0 limit is not singular at strong coupling and one can compute ∆ ∞ S=0 directly from (4.36).
The spectrum of exchanged operators in reductions of χCFT. In section 2.2, we presented the γ-deformed N = 4 SYM theory in the double-scaling limit as a family of theories. In fact, playing with the three couplings ξ j with j = 1, 2, 3 it is possible to describe different Lagrangians with different matter contents and symmetries. Thus we want to obtain the spectrum of exchanged operators for each theory of this family simply taking the limit on the couplings in the spectral equation (4.28) of the most general doubly-scaled theory. First of all, we recall the well-known result for the spectrum for the simplified Lagrangian (2.6) also known as 4D bi-scalar fishnet CFT. In this theory the only non-trivial four-point correlation function is G φ 1 φ 1 , and it can be written in the same OPE form as the one we are considering as (3.27). By the Bethe-Salpeter method it is possible to compute the correlator at all-loops, since its perturbative expansion is generated only by a bosonic graph-building operator H B of (4.7), then we can extract the non-perturbative scaling dimension of the exchanged operators in the OPE s-channel. The corresponding spectral equation is the same of (4.28) with ω 4 = ξ 4 and κ 4 = 0 (indeed the bi-scalar theory has only one coupling ξ 2 ) and it has two solutions corresponding to the twist-two and -four operators with the following scaling dimensions bi (ξ 4 ) = 2 + (S + 1) 2 + 1 + 2 (S + 1) 2 + 4ξ 4 , (4.37) together with two shadow solutions with ∆ = 4 − ∆ for Re ∆ < 2. The analytic properties of those solution and their weak-and strong-coupling expansions have been studied in detail in [24].
The scaling dimensions of the exchanged operators in the correlation function G φ 1 φ 1 for theories defined as a reduction of χCFT as in section 2.2, can be computed as solutions of the spectral equation (4.28) in which we are applying some limits on the couplings, or JHEP06(2019)078 even directly on the weak-and strong-coupling expansions. In the table 2 we present the summary of our results.
• χ 0 CFT: since the spectrum of the exchanged operators for the four-point function G φ 1 φ 1 doesn't depend on ξ 1 , the limit in which one of the couplings of the full χCFT is going to zero (reducing the theory to the χ 0 CFT) is not unique. Indeed if we set ξ 1 = 0, the scaling dimensions of the exchanged operators in the χ 0 CFT are the same of the full χCFT. On the contrary, if we set ξ 2 or ξ 3 to zero, the spectrum reduces to that of the bi-scalar theory (4.37) depending on a single coupling. Notice that in this case the number of solution of twist-four operators reduces to a single one, while the higher-twist operators get protected.
• bi-scalar theory: the reduction to bi-scalar theory corresponds to the limit in which two couplings of χCFT vanish. If one of the vanishing couplings is ξ 1 , the spectrum is the usual one of the bi-scalar theory (4.37) while if ξ 2 = ξ 3 = 0 the operators are protected because the only remaining interaction vertex is not contributing.
• β-deformed theory: when all the couplings are equal we reduce the full theory to its βdeformation. In this case, due of the restoration of one supersymmetry, the operator of twist-two is protected as pointed out in [19] and confirmed by our computation (this reduction in terms of the new couplings κ and ω corresponds to κ → ξ and ω → 0). The symmetry doesn't constrain the operators of twist-four to be protected, as well as for higher twist t > 4. Indeed, their spectrum can be easily read applying the limit, for example at weak coupling, to the expansions (4.30), (4.31).

The structure constant of the exchanged operators
Once the spectrum of the exchanged operators is computed, in order to obtain the full set of conformal data for the four point function G φ 1 φ 1 , one has to compute the OPE coefficients. From their definition (3.27), we get

JHEP06(2019)078
where Here c 2 is given in (3.13) and one puts ∆ O 1 = ∆ O 2 = 1. The eigenvalues h B ∆,S andh F ∆,S are presented in (4.14) and (4.25), and the constants c F = −2c B = −1/(2π 2 ). Plugging these eigenvalues into (4.39) and performing the derivative, we obtain a rather cumbersome result that we will present in the next paragraph.
Weak coupling expansion. Performing the derivative in (4.39) and substituting the weak coupling expansions of the scaling dimensions computed in (4.29) and (4.30), we obtain the following expansions for the structure constants associated to the exchanged operators for S = 0 k are harmonic numbers. Again, the expressions in square brackets are in fact rational numbers. Similarly to the expansion of the scaling dimension, the OPE coefficient of the twist-two operator is singular for S = 0. Indeed, as discussed in section 4.3, due to the singularity arising at zero spin, the weak coupling and S → 0 limits don't commute. In order to obtain the correct weak coupling expansion for the twist-two operator, one has to set S = 0 in (4.38) and then expand it in the coupling. The zero spin expansion the OPE coefficients of exchanged operators reads  In analogy with the spectrum analysis, the power counting shows that the twist-two operator goes in power of ξ 4 as expected if S = 0. In the S = 0 case it is going in powers of ξ 2 , suggesting that the weak coupling expansion is sensitive to the double trace counterterms.

JHEP06(2019)078
Moreover in both cases (4.40) and (4.41), the twist-four OPE coefficients are suppressed by a factor of order ξ 2 as compared to those of the twist-2.
Strong coupling expansion. Since we know from the expansion at strong coupling of the scaling dimension (4.36) that the scaling dimension becomes large, we can expand (4.38) in the limit ∆ → ∞ and obtain 42) where we have to substitute ∆ from the strong coupling spectrum ∆ ∞ computed in (4.36) for low-twist operators. Naively, the expansion (4.42) looks the same as the one of the structure constant of the bi-scalar model [24], but actually it is not. Indeed, one can notice from the definition (4.38) that the OPE coefficient in our model depends explicitly on the coupling. Then in the expansion (4.42) some coefficients at higher order will start to depend on κ 4 . The first contribution different from the bi-scalar expansion appear as κ 4 /∆ 6 which, after the substitution ∆ ∞ , contributes at order O(1/ξ 2 ) in the inverse coupling expansion. Hence, it is convenient to write (4.42) as follows where k = 0, 1, 2, 3 labels the four solutions of the spectral equation (4.28) and dots stand for higher orders in 1/κ and 1/ω. Thus, given the scaling dimension ∆ ∞ (4.36), the OPE coefficient is exponentially small at strong coupling due to the factor 1 2 2∆∞ . The S → 0 limit is not singular at strong coupling and one can compute C ∆∞,0 directly from (4.43).

The four-point correlation function
Once the conformal data in sections 4.3 and 4.4 is computed, one can determine the fourpoint function (4.1) by means of (3.19). In the case O 1 = O 2 = φ 1 we obtain with the cross-ratios defined as u = x 2 12 x 2 34 /(x 2 13 x 2 24 ) and v = x 2 14 x 2 23 /(x 2 13 x 2 24 ) and ∆ φ 1 = 1. The function G φ 1 φ 1 (u, v) can be written in terms of the OPE representation (3.27) as a sum over the non-negative integer Lorentz spin S and the states with scaling dimensions ∆. From the study of the spectrum of exchanged operators in section 4.3 it turns out that infinitely many operators are exchanged in the OPE channel. Then we have

JHEP06(2019)078
where the scaling dimensions ∆ (i) are defined by the spectral equation (4.28) and computed at weak coupling in (4.29), (4.30) and (4.31), and for low twist t = 2, 4 at strong coupling in (4.36). The structure constants C ∆ (i) ,S associated to the exchanged operators are defined by (4.38) are computed at weak coupling in (4.40) and for low twist and strong coupling in (4.43). The four-dimensional conformal blocks g ∆,S are defined in (3.18). The proper definition of the four-point correlation function G φ 1 φ 1 takes into account the symmetrization x 3 ↔ x 4 . Under this symmetry, the cross-ratios transform as u → u/v and v → 1/v. Correspondingly, from the definition (3.18) the conformal blocks obey the symmetry g ∆,S (u/v, 1/v) = (−1) S g ∆,S (u, v). Combining together this relation with (4.45), it's easy to see that, imposing the symmetry x 3 ↔ x 4 , the terms in (4.45) with odd S cancel out whereas those with even S get doubled.
Despite of the presence of singularities in the weak-coupling expansions of scaling dimensions and OPE coefficients, their combination in (4.45) is well-defined. Indeed, plugging the conformal data into (4.45), we obtain an expansion in powers of the couplings that is compatible with the interpretation of the correlation function as a sum of Feynman diagrams in perturbation theory (see section 6.2 for an explicit example). In particular, since the first non-trivial order is fixed by the S = 0 conformal data, it is easy to write the very first contributions to G φ 1 φ 1 in terms of the known functions, as follows where Φ (L) is the ladder three-point function [52] that in the case L = 1 is given by the Bloch-Wigner dilogarithm function (4.48)

Exact four-point correlations function for
In this section we consider the four-point correlators associated to the last operator of (2.3), namely when O 1 (x) = φ j (x) and O 2 (x) = φ † k =j (x) with j, k = 1, 2, 3. Since the computation of the correlators is the same for any j and k, we will consider the case j = 1 and k = 2 and then the four-point function we want to study takes the following form This correlation function was trivial in the bi-scalar model [12] but in the general doublescaled theory (2.2) it has a rich diagrammatic structure. Indeed a generic Feynman diagram in the weak coupling expansion of (5.1) in the planar limit N c → ∞ is given by a combination of all the single-trace vertices in (2.2) and the following double-trace vertex coming from the counterterm Lagrangian (2.7). In the following we will compute the conformal data of (5.1) with the Bethe-Salpeter method.

The Bethe-Salpeter method for the correlator
The perturbative expansion of (5.1) can be written in the following form at any perturbative order contains contributions from the bosonic and fermionic integrals, with different dependence on couplings. In figure 12, we present an example of a generic Feynman diagram contributing to (5.3). As in the previous case, the Feynman diagrams of this four-point correlation function have a cylindric topology and, at arbitrary order , they take an iterative form allowing us to write the full expansion as an infinite geometric sum of the primitive divergencies, as in (3.2). In contrast to the previous case, the nodes of the chain diagrams in the expansion of G φ 1 φ † 2 are not only insertions of double-trace vertices but also of the single-trace vertex In the Bethe-Salpeter procedure, both vertices enter only as insertions of the operatorV defined in (3.22). Then it's easy to conclude, as it was done in [24] for the biscalar fishnet model, that the coefficient of this operator in (3.2) is In section 3, we discussed the role of the operatorV in the computation of the four-point function, arguing that it is not contributing to the spectral equation for finite coupling or S. Then, as far as we consider the perturbative expansion (5.3) in the point-splitting The first few orders of the perturbative expansion are given by the diagrams represented in figure 13. They can be written as follows where each scalar propagator brings in the factor c B /x 2 ij and each fermionic propagatorthe factor c F / x ij /x 4 ij , where / x can be σ µ x µ orσ µ x µ and c F = −2c B = −1/(2π 2 ). These diagrams can be expressed in terms of a combination of the Hamiltonian graphbuilding operatorsĤ i . Indeed, considering the bosonic part of (5.6), the bosonic kernel is that is clearly the same as studied in the previous case (see section 4). The fermionic kernel is more involved. Considering the diagrams in figure 13(c) and 13(d) and their integral representation (5.6), it is clear that they are not generated by the same repeated Hamiltonian operator. In fact, going on with the perturbative expansion of G φ 1 φ † 2 , one can notice that at any order for > 1, at least one fermionic diagram with the following ladder topology appears Since any of them carries a ξ 2 power of the coupling, the maximum number of bosonic rungs in the ladder depends on the perturbative order, in particular is − 2. However for > 3, also the superpositions of ladders with less rungs contribute. For these reasons it is JHEP06(2019)078 H r (y 2j+1 , y 2j+2 |y 2j−1 , y 2j )H t (y 1 , y 2 |x 3 , x 4 ), (5.8) where the fermionic sub-kernels H b , H r and H t are contracted in the spin indices in order to recompose the trace of σ-matrices. In our convention, when n = 0 the rung-building operator H r is not contributing to H F . These Hamiltonians, graphically represented in figure 14, are defined as follows and they can be used to rewrite the expansion (5.6) obtaining F (x 1 , x 2 |y 1 , y 2 ) H B (y 1 , y 2 |x 3 , x 4 ) , where, using the definition (5.8), we have d 4 y i Hβ β b (x 1 , x 2 |y 3 , y 4 )H r βα;αβ (y 3 , y 4 |y 1 , y 2 )Hα α t (y 1 , y 2 |x 3 , x 4 ) .
The kernels (5.7) and (5.8) transform covariantly under conformal transformations, then the corresponding Hamiltonian integral operators commute with the generators of the conformal group. The fermionic sub-kernels (5.9) have spinorial indices carried by the σ-matrices and transform as two-components spinors. Following the conventions we are using for the raising and lowering of spin indices, as explained in appendix A, we have the following transformations (5.12) which corresponds to the exchange σ ↔σ. The rung-building operator H r contains a couple of un-contracted σ-matrices, then it appears with two pairs of indices. In order to build the general fermionic diagram, one has to contract the fermionic sub-kernels with the only constraint to obtain the trace of all the σ-matrices around the fermionic loop alternating σ's withσ's. Once chosen if the top sub-kernel H t contains the combination σσσ orσσσ, the first rung sub-kernel H r has to have the right combination of indices to be contracted, in particularσσ and σσ respectively. Then the other H r kernels have to alternate upper and lower indices. Depending on parity of the number of rungs of the ladder n, the bottom sub-kernel H b can carry σ orσ. Indeed, we can distinguish two different index structures, for odd or even number, of repeated applications of H r as follows Carrying on in the perturbative expansion, one can find that for example the perturbative order = 3 is given by the sum of same combinations of kernels appearing at order F and so on. For this motivation, the -th perturbative order G cannot be written as the contribution for = 1 to the power as in the case studied in section 4, but its sum takes the following form Since the operatorial form of the fermionic Hamiltonian (5.8) in terms of the subkernels (5.9) isĤ

JHEP06(2019)078
one can sum the two geometric series in (5.14). Then the correlator (5.3) can be written as followsĜ Finally, comparing (4.10) with the definition (3.2), we can fix the value of the remaining constants χ i

Eigenvalues of the Hamiltonian graph-building operators
In order to compute the four-point correlation function with the operator method presented in section 3, one has to obtain the spectrum of the graph-building operators (5.7) and (5.8) as was done in section 4.2 for another four-point correlator. The eigenstate that diagonalizes all these Hamiltonians is defined in (3.5) for ∆ O 1 = ∆ O 2 = 1 and the eigenvalues are defined by means of equations (3.9) and (3.10). Substituting in the latter the kernels (5.7) and (5.9) and using the definition (3.7), we will end up with a set of integrals that can be computed with the help of the star-triangle relations (C.2). The fact that all the integrals that we have to compute can be solved by means of the star-triangle relations is a strong evidence of the underlying conformal symmetry.
Bosonic eigenvalue. Since the bosonic Hamiltonian (5.7) is the same as studied in the previous case, its eigenvalue is already computed: it is given by (4.14) for any S and by (4.13) in the case of S = 0.
Fermionic eigenvalue. The fermionic eigenvalue is defined in (3.10). In the case we are studying, the Hamiltonian operator (5.8) depends on the number of rungs n, thus its eigenvalue will be also a function of n, as follows The computation of the fermionic spectrum is simpler if we consider the Hamiltonian H (n) F in terms of the sub-kernels H b , H r and H t , as in (5.15). While the full fermionic Hamiltonian is diagonalized by the state Φ ∆,S,x 0 , the sub-kernels are not. Indeed, due to their fermionic structure, individually they will turn the scalar conformal eigenfunction into a fermionic state. However, compatibly with (5.19), their combination (5.8) will leave the state Φ ∆,S,x 0 unchanged. For instance, at S = 0 when the state is given by the conformal where is a non-negative integer, p = ∆/2, black and white dots are positions with and without integrations over R 4 , and we defined the state (5.21) In (5.20) we consider the action of an even and odd number of rung-building operators separately. Indeed, while the top sub-kernel is changing the conformal triangle into a fermionic object and the bottom sub-kernel is turning it back to the original state, the n copies of the operatorĤ r are exchanging σ ↔σ in the states (5.21) depending on the parity of n, according to (5.13).
These arguments hold also in the case of S = 0, then we have to solve the following equations 22) and the same for the transformed kernels by means of (5.12). Considering the first two equations of (5.22) and the definitions (5.9), we have Moreover, focusing on the last equation of (5.22) and the definition of bottom sub-kernel given in (5.9), we have where we used equation (A.10) to simplify the fermionic structure. The integral in the right-hand side of the equation is the same as appearing in the computation of the bosonic eigenvalue related to the operatorĤ B given by (4.12). Then eq. (5.25) together with (5.22) and (4.12) leads to where h B ∆,S corresponds to (4.14) for any S and (4.13) in the S = 0 case. With these arguments it is clear that, in order to compute the fermionic eigenvalue of the operator (5.15), we need only to compute the eigenvalue h t∆,S defined by the first equation in (5.22). In the S = 0 case the computation of h t∆,S can be performed using only star-triangle relation (C.2). Indeed starting from the integral (5.23) together with the first equation of (5.22) and performing two fermionic star-triangle integrations, as follows (see footnote 8) we obtain . (5.28) The computation of the eigenvalue h t∆,S for S = 0 is more involved. Without any loss of generality, we can consider the limit x 0 → 0 in both sides of the first equation of (5.22) in order to simplify the spin structure. In this limit we are able to compute the resulting integral going to momentum space. We leave the details of the computation in appendix D. Going through the calculation we obtain (D.7). Notice that also this eigenvalue can be written in terms of the bosonic one as follows  (5.30) where in the last line we used the definition (4.14). Notice that setting S = 0 we obtain (5.28) as expected. Finally, using the definition (5.17), we can resum the fermionic eigenvalue obtaining We can conclude that the eigenvalue (5.31) is manifestly invariant under ∆ → 4 − ∆, as expected from (3.21).

Spectrum of exchanged operators of
In this section we will use the eigenvalues (4.14) and (5.31) to compute the scaling dimensions of the operators contributing to the correlation function (3.27) for The spectrum of the exchanged operators is defined by the solutions of the equation for the physical poles (3.25). Substituting in (3.25) the definition of bosonic and fermionic eigenvalues (4.14) and (5.31) and the constants χ i computed in (5.18), we can rearrange the spectral equation in the following form where we defined the new couplings Plugging (4.14) into (5.32), we obtain the following equation but instead of these bulky formulas we prefer to present their perturbative expansions at weak and strong coupling.
The equation (5.34) can be viewed as defining an algebraic curve of 8th degree in terms of ν, or of the corresponding dimension ∆ = 2 + 2iν as a complex variable. In terms of ν 2 variable it is a 4th order algebraic curve whose branches -the four sheets of the related Riemann surface -describe directly the four physical dimensions as functions of couplings. Changing the couplings we can pass from one sheet to another, observing the transitions between various dimensions. The branch points correspond to the collisions of physical dimensions. In contrast to this case, the spectral equation for the previous four-point function (4.28) is not algebraic and its Riemann surface contains infinitely many sheets.
Weak coupling expansion. Expanding around the physical pole ν = −iS/2 at weak coupling λ, µ → 0, we obtain the following expansions of dimensions for the two twist-two operators where we introduced the following short-hand notation ζ = S(S + 1)λ 2 − 16µ 2 and τ = (S + 1)(S + 2)λ 2 − 16µ 2 . (5.37) In both cases we presented only the first few terms of the expansions since the following ones are quite cumbersome. Moreover, we notice that the weak coupling expansions of the JHEP06(2019)078 solutions are divergent but, as we point out later, by the analysis of section 5.5 the sum of the corresponding OPE contributions has a well defined expansion (see footnote 9). Taking into account the quantum numbers of the external state Tr[φ 1 φ † 2 ] we can list the operators of twist-2 and -4 which mix, obtaining the expression for the exchanged operators in the OPE from the diagonalization of the mixing matrix. At ∆ 0 − S = 2 we recognize that Tr[φ 1 ∂ S + φ † 2 ] and Tr[ψ 1 ∂ S−1 + (n · σ)ψ 2 ], where n · n = 0 and ∂ + = (n · ∂), operators can mix if S > 0. Interestingly, due to trace ciclicity, at twist-2, the mixing transitions are symmetric and then there are no logarithmic operators despite chiral interactions of (2.2) (see appendix F.2). On the other hand, at ∆ 0 − S = 4 many more operators can realise the U(1) ⊗3 quantum numbers (1, −1, 0), building up a large mixing matrix. As pointed out in appendix F.2, already for S = 0 the anomalous dimension matrix is not diagonalizable. We conjecture there that its canonical Jordan form presents several logarithmic multiplets (Jordan blocks) of various ranks, together with two non-zero diagonal elements, corresponding to solutions (5.36) (or (5.39) and S = 0).
The zero-spin case presents some peculiar behaviors. Indeed, expanding (5.34) for S = 0 around the physical poles ν = 0 at weak coupling, we obtain the following expansion for the twist-two operator and around the pole ν = −i -the following twist-four operators The fact that at S = 0 there is only one solution corresponding to operators of length-two is not surprising: the only scalar operator which has twist-2 and the correct quantum numbers is indeed Tr[φ 1 φ † 2 ], that is the one exchanged in the OPE channel (see appendix F.2). Moreover, similarly to the case of the spectrum of the length-two operator in the correlation function computed in section 4.3, the limit S → 0 and the weak-coupling limit are not commutative, thus the scaling dimension ∆ (2) for S = 0 presents a different expansion in power of the coupling w.r.t. the case S > 0. The explanation is the same of the previous case. Indeed, since the physical pole is situated at ν = −iS/2, and the mirror pole at ν = +iS/2, at weak coupling, when S → 0, they pinch the integration contour of (3.24) at the origin producing a divergence. Hence, the contribution of the double-traces is needed in this case to produce a non-vanishing term that cancels this divergence at weak coupling. Again we stress that at finite couplings the solutions (5.34) are well-defined even at zero spin.
Strong coupling expansion. Equation (5.34) has 6 solutions at strong coupling i.e. for λ, µ → ∞. The explanation of the fact that we have a different number of solutions JHEP06(2019)078 Table 3. In this table we summarize the operator and dimension content of exchange operators in three reductions of χCFT for the correlator G φ1φ † 2 . In our notation and ∆ (i) bi is defined in (4.37).
w.r.t. the weak coupling case can be identified by the analysis of the exchanged operators in appendix F.2. In particular when λ, µ → ∞ the two operators O ± defined in (F.10) collapse to a single operator. Then the resulting physical solutions are The spectrum of exchanged operators in reductions of χCFT. The scaling dimensions of the exchanged operators in the correlation function G φ 1 φ † 2 for the reductions of χCFT introduced in section 2.2, can be computed as solutions of the spectral equation (5.34) in which we are applying the limits on the couplings, or even directly on the weak-and strong-coupling expansions. In table 3 we summarize our results.
• χ 0 CFT: sending one of the couplings to zero we perform the reduction from the full χCFT to the χ 0 CFT. In this case, the reduction doesn't give a unique spectrum. Indeed, setting ξ 1 or ξ 2 to zero, or equivalently λ → 0, we obtain two protected operators. This is clear when considering the spectral equation (5.34) in which only the first term on the left-hand side will contribute and clearly it describes two protected solutions (plus the shadow operators associated to them). If we consider ξ 3 to vanish, or equivalently µ → 0, and again referring to the spectral equation (5.34), it is clear that the second term in the left-hand side vanishes but, in contrast to the previous case, in the right-hand side of the equation the dependence on the coupling is still present. The resulting equation is the same as describing the spectrum of exchanged

JHEP06(2019)078
operators for the four-point function G φ 1 φ 1 in the bi-scalar theory. Then we obtain 2 solutions with the spectrum given by (4.37) where ξ = λ.
• bi-scalar theory: the reduction to bi-scalar theory consists of the limit in which two couplings of χCFT vanish. For the correlator we are considering, any choice of the vanishing couplings leads to protected solutions, as expected. Indeed, in the bi-scalar theory the lack of interaction vertices doesn't develop an anomalous dimension for exchanged operators with external states given by φ 1 and φ † 2 as previously noticed in [12,24].
• β-deformed theory (all three couplings are equal): in this case, the restoration of one supersymmetry is not sufficient to constrain any of the solutions. Then their spectrum can be easily read off applying the equal couplings limit, for example at weak coupling, to the expansions (5.35) and (5.36).

The structure constant of the exchanged operators
Once the spectrum of the exchanged operators is computed, in order to obtain the full set of conformal data for the four point function G φ 1 φ †

2
, one has to compute the OPE coefficients. Their definition is given by (3.27) or, equivalently, by (4.38) where where the eigenvalues h B ∆,S is defined in (4.14) and the constants c F = −2c B = −1/(2π 2 ). Plugging the eigenvalues into (5.41) and performing the derivative, we obtain explicit but rather cumbersome result that we will not present here. In the following paragraphs we will consider its weak-and strong-coupling expansions.
Weak coupling expansion. Substituting the weak coupling expansions of the scaling dimensions computed in (5.35) into (4.38), (5.41), we obtain the following weak coupling expansions for the structure constants associated to the twist-two operators for S = 0 and a similar substitution of (5.36) for the twist-four operators at any S gives as discussed in section 5.3, due to the singularity arising at zero spin, the weak coupling and S → 0 limits don't commute. In order to obtain the correct weak coupling expansion for the twist-two operators, one has to set S = 0 in (5.41) and then expand in the couplings.
In the S = 0 case we found that the spectral equation (5.32) has 6 solutions instead of 8, as expected from the operator mixing analysis in appendix F.2 for twist-2 solutions.
In particular the zero spin expansion of the structure constant related to the twist-two operator is Since the S → 0 singularity is not arising for the twist four operators, one can read off the spinless OPE coefficients for the twist-four operators directly setting S = 0 in (5.43), thus obtaining

(5.45)
In analogy with the analysis of spectrum, the divergence arising at S = 0 for the twisttwo operators suggests that the weak coupling expansion is sensitive to the double trace counterterms. Moreover, for any S the twist-four OPE coefficients are suppressed by a factor of order ξ 3/2 as compared to the twist-two ones.
Strong coupling expansion. In section 5.3 we studied the strong coupling spectrum of the exchanged operator for the four-point function Since the scaling dimensions are growing linearly with the coupling λ, the OPE coefficients can be computed expanding (4.38) and (5.41) in ∆ → ∞ and then plugging the strong coupling asymptotics of scaling dimensions in the resulting expansion. Since in this limit the coefficient c 2 (∆, S) in (4.38) is dominant w.r.t. the one of (5.41), the first terms of the expansion are exactly the same as in (4.42). Obviously at higher orders in 1/∆, (5.41) starts contributing and then we obtain Thus, given the linear growth of the scaling dimensions ∆ ∞ with coupling, the OPE coefficients are exponentially small at strong coupling. The S → 0 limit is not singular at strong coupling, thus one can compute C ∆∞,0 directly from expansion (5.46), setting there S = 0.

The four-point correlation function
Once we computed the conformal data in section 5.3 and section 5.4, we can determine the four-point function (5.1) by means of (3.19). In the case O 1 = φ 1 and O 2 = φ † 2 we obtain The cross-ratios are defined as u = x 2 12 x 2 34 /(x 2 13 x 2 24 ) and v = x 2 14 x 2 23 /(x 2 13 x 2 24 ) and ∆ φ i = 1. The function G φ 1 φ 1 (u, v) can be written in terms of the OPE representation (3.27) as a sum over the non-negative integer Lorentz spin S and the states with scaling dimensions ∆. From the study of the spectrum of exchanged operators in section 5.3 it turns out that in the OPE channel two operators of length-two for S = 0 are exchanged, but only one when S = 0. Thus, in order to write the four point function in a compact form, we set C ∆ (2 ) ,S to vanish at S = 0, or equal to (5.42) otherwise, then we have As already pointed out for the four-point correlator G φ 1 φ 1 , the proper definition of the four-point correlation function G φ 1 φ † 2 , takes into account the symmetrization x 3 ↔ x 4 . Under this symmetry, the cross-ratios transform as u → u/v and v → 1/v and correspondingly, from the definition (3.18), the conformal blocks g ∆,S (u/v, 1/v) = (−1) S g ∆,S (u, v). Combining together this relation with (5.48), it's easy to see that including the symmetry x 3 ↔ x 4 means that the terms in (5.48) with odd S cancel out whereas those with even S get doubled.
Despite the presence of singularities in the weak-coupling expansions of scaling dimensions and OPE coefficients, their sum in (5.48) is well-defined and non-singular. Indeed plugging the conformal data into (5.48), we obtain an expansion in powers of the couplings that is compatible with the interpretation of the correlation function as a sum of Feynman diagrams in perturbation theory (see section 6.3 for an explicit example). In particular, since the first non-trivial order is fixed by the S = 0 conformal data, it is easy to write the two leading contributions to G φ 1 φ † 2 in terms of known functions as follows where Φ (1) is the Bloch-Wigner dilogarithm function defined in (4.47). Expanding (5.48) at higher order in the couplings we obtain a cumbersome result. However, we notice that the maximum transcendental weight of the involved functions grows linearly with the order.

JHEP06(2019)078
This suggest, together with the presence of the Bloch-Wigner dilogarithm at the first nontrivial order, the possibility to interpret the four-point function (5.48) as a combination of a class of special iterated integrals, the so-called harmonic polylogarithms [53] following the same idea as of [24].

Correlation functions at weak coupling from Feynman diagrams
In section 4 and section 5, we have analyzed two different four-point functions computing their conformal data by means of the Bethe-Salpeter method. With this procedure we were able to diagonalize the graph-building operators and write exact equations for the spectrum of exchanged operators even though we ignored on the way the contribution of the doubletrace interactions (2.7). The double-trace counterterms are necessary in the action to have a consistent description of the double-scaled theory (2.2) in the perturbative regime and in particular for the restoration of conformal symmetry.
In this section, we will study the weak coupling expansions of the four-point functions related to the operators (2.3) and clarify the role of the double-trace terms for this expansion. As we already mentioned, bosonic and fermionic wrappings in the related Feynman graphs develop UV divergencies at short distances. Adding the double-trace vertices in the perturbative expansion we will be able to determine the conformal fixed points (2.13) canceling divergencies generated by the single trace terms. However, they will not affect the finite coupling solutions computed in the section above.
The double-trace counterterms are given by the Lagrangian (2.7). In general, this action contains 9 terms but, due of the cylindric topology of Feynman diagrams for the observables we are computing. For any four-point function only one double-trace term is contributing generating a new local four-scalar vertex (see for instance figure 9). This fact is crucial to ensure that conformal symmetry is restored. Indeed, in this case we know that the β-function can be written as (2.11) and it admits two fixed points α 2 j as in (2.13). If we focus on the Feynman diagram expansion of the four point-functions we have to deal with divergent integrals. Then we have to introduce dimensional regularization setting D = 4 − 2 . One important observation is that the diagrams containing fermionic contributions produce the same divergence as the bosonic ones. In other words the fermionic kernels contains the divergent part of the bosonic one plus a remainder function that is finite in D = 4 which therefore does not require regularization. Then the divergent operator can be written as where Φ(x 1 , x 2 ) is a test function.

JHEP06(2019)078
Given the Hamiltonians (6.1) and the definition (3.2), the four-point correlation function is defined as follows it is taken at the fixed point α 2 j = α 2 j . It is clear that for = 0 conformal symmetry is broken. However, expanding (6.2) at weak-coupling in terms of Feynman diagrams, one can demonstrate order-by-order how conformal symmetry is restored. In the following sections, we present some examples of this mechanism for the four-point correlation functions associated to the operators (2.3).
Let's start from this simplistic example. Consider O 1 (x) = φ j (x) and O 2 (x) = φ k =j (x) with j, k = 1, 2, 3. Since the correlators are the same for any j and k, we choose j = 1 and k = 2. Then the four-point function we want to study is as follows In the planar limit N c → ∞, the weak coupling expansion of (6.3) in terms of Feynman diagrams is given by a combination of the following vertices and it can be written as follows where represents the perturbative order. It is straightforward to see that Feynman diagrams at any order form a chain structure alternating groups of single-and doubletrace vertices. These vertices insert into the graphs identical primitive divergencies. Then at the conformal critical point α 2 3 = −ξ 2 3 the graphs cancel each other except for the = 0 term, namely Restoring the two point function in the limit x 1 → x 2 and x 3 → x 4 we notice that the spectrum of the operator Tr[φ 1 φ 2 ] is not affected by quantum corrections. Indeed, its scaling dimension is protected and equal to the bare one In section 4 we studied the contribution to the four-point function G φ 1 φ 1 of diagrams generated by the bosonic and fermionic HamiltoniansĤ B andĤ F but ignored the doubletrace vertices, which works well for finite couplings. Since any diagram in figure  for the trivial leading order diagram) is UV divergent, we will re-introduce in this section the double-trace counterterms in the perturbative expansion in order to make the weakcoupling expansion UV finite and to restore conformal symmetry.
Let's compute the first few orders of the weak coupling expansion of G φ 1 φ 1 . In terms of Feynman diagrams, we have to compute the graphs given in figures 10 and 15. Defining a function G The first correction is given by the diagram in figure 15(a). This contribution is finite and it can be written as follows where Φ (L) is the ladder three-point function [52] that in the case L = 1 is given by the Bloch-Wigner dilogarithm function The cross-ratios are u = x 2 12 x 2 34 /(x 2 13 x 2 24 ) and v = x 2 14 x 2 23 /(x 2 13 x 2 24 ) and the constant c B = 1/(4π 2 ).
The bosonic part of the two-loop correction given by G where we defined the short-hand notation

JHEP06(2019)078
This integral is UV divergent at short distances y 2 12 → 0. Using the identity 1/(y 2 12 ) 2−2 = π 2 δ 4−2 (y 12 )/ + O( 0 ), one can compute the divergent part of the integral (6.13) that is proportional to the same one-loop function found in (6.9), as follows where for the purpose of this section we are not interested in the finite part. Let us finally consider the fermionic contribution G (0,0,4) φ 1 φ 1 . This term corresponds to the Feynman diagram in figure 10(d) and its integral representation in four dimensions is given in the last line of (4.6). Since we are only interested in the computation of the UV divergent part of the diagram, we can avoid computing the whole integral (4.6) in dimensional regularization and proceed in a more naive way. Indeed, representing the integral as in the last line of (4.8) and considering that the fermionic Hamiltonian can be written as a combination of the bosonic one and some finite reminder function as in (4.16), we know that all the UV divergence is arising from H B . Then we can write Combining this result with (6.14) and (6.9), we obtain that the expansion (6.8) takes the expected form (3.19) with the function G φ 1 φ 1 (u, v) given by (6.16) where the new couplings κ and ω are defined in (4.27) and finite(κ 4 , ω 4 ) stands for the finite part of G φ 1 φ 1 at two-loop. Finally, imposing the UV finiteness of the correlation function we obtain the first order of the fixed point as follows 17) and notice that it matches exactly the prediction given in (2.16). Replacing the doubletrace coupling in (6.16) with its value 10 α 2 1 = − i ω 2 2 we obtain the one-loop expansion of the correlation function as follows This perfectly matches the same quantity computed via OPE, with conformal data fixed by the Bethe-Salpeter method (4.46).
the weak coupling perturbative expansion of the correlator using Feynman diagrams, the graphs generated by this operator are needed in order to cancel the UV divergencies arising from the diagrams in figure 13. There is here a substantial difference w.r.t. the case of G φ 1 φ 1 studied in section 4: indeed the diagrams generated byV are not only produced by the double-trace vertex . For this reason, we redefine the double-trace coupling as (5.5) taking into account both contributions.
Let's compute the first few orders of the weak coupling expansion of G φ 1 φ †

2
. In terms of Feynman diagrams, we have to compute the first three graphs given in figure 13 and the ones in figure 16. Defining the function G where a counts the number of double-and single-trace vertices contributing with couplingα 2 2 and b and c the number of bosonic and fermionic vertices, respectively, we have The first correction is given by the diagram in figure 16(a) that is a half 11 of the one computed in the previous section in (6.9). The bosonic part of the two-loop correction given by G Their integral representation is exactly the same as in the previous case, namely (6.12), and then we can write their sum as follows Let's finally consider the fermionic contribution G This term corresponds to the Feynman diagram in figure 13(c) and its integral representation in four dimensions is given in the next-to-the-last line of (5.6). In this case, our goal is also to identify the UV divergent part of the diagram. Thus we will not compute the whole integral in (5.6) in dimensional regularization but we will rather proceed following the method of the previous section. Performing two integrations by means of the star-triangle relation (C.2) and then

JHEP06(2019)078
simplifying the spin structure with the help of (A.11) and (A.4), one can identify a divergent integral of the same kind as in (6.13), together with a finite remainder integral. Computing the integral I for small distances y 12 → 0, we can extract the following pole Combining all this results, we obtain that the expansion (6.19) takes the expected form (3.19) with the function G φ 1 φ † 2 (u, v) given by where finite(λ 4 , µ 4 ) stands for the finite part of G φ 1 φ † 2 at two-loop. Finally, imposing the UV finiteness of the correlation function and recalling the definition (5.5), we obtain the first order of the fixed point as follows Replacing the effective coupling in (6.22) with its valueα 2 2 = iλ λ 2 + 2µ 2 , 12 we obtain the one-loop expansion of the correlation function in the following form that perfectly matches the same quantity computed via OPE with conformal data fixed by the Bethe-Salpeter method (5.49).

Conclusion and discussion
This paper represents an attempt of a deeper understanding of physical properties and analytic structure of the four-dimensional, three-coupling chiral CFT -the χCFT 13proposed byÖ. Gürdȏgan and one of the authors in [18] as a double scaling limit of γdeformed N = 4 SYMtheory, combining the weak coupling with the strong imaginary γ-twist. We study here two aspects of this χCFT with three effective couplings, given by the Lagrangian (2.2): i) the explicit description of the Feynman graph content of the perturbative expansion, partially uncovering their integrability properties; ii) the exact computation, via conformal symmetry, of two four-point correlation functions of shortest protected scalar operators of the theory. As concerns the planar Feynman graphs content of the theory, we found here the complete description of possible Feynman graphs in the "bulk" -inside a generic big Feynman graph, far from its boundaries defined by a particular studied physical quantity. These graphs can be dubbed as "dynamical fishnet", since, unlike the usual regular fishnet of the bi-scalar model (2.6) they have a certain dynamics (summation over many of such graphs) preserving at the same time a kind of irregular fishnet structure shown on figures 2, 3.

JHEP06(2019)078
Interestingly, this bulk structure is neatly realized as Feynman graphs describing arbitrary single-trace correlation functions of all elementary fields, as shown on figures 4, 5. It would be very interesting to find the realisation of the Yangian symmetry of these correlators, and of the related planar amplitudes (with disc topology), generalizing the results of [25,26] for the bi-scalar CFT. It would be the neatest demonstration of the integrability of the full model. In section 2 we demonstrate such integrability in the two-coupling reduction of the full χCFT, having a much simpler fishnet structure (combination of regular "brick wall" graphs with Yukawa vertices and regular square lattice fishnets). A considerably more involved analysis of the integrability of the full dynamical fishnet of χCFT, in particular, via the Yangian symmetry of single-trace correlators, is underway. We believe that it will be another important step to the understanding of integrability of the mother theorythe N = 4 SYM. It is worth noticing here that γ-deformation represents a rather mild, "topological" modification of the planar graphs of original N = 4 SYM, altering only the boundaries of these graphs, and not the bulk.
In the second part of our paper, we managed to compute two non-trivial four point correlation functions of elementary fields of the full three coupling χCFT, generalizing the bi-scalar fishnet CFT results of [12,24]. As in these papers, we employed the Bethe-Salpeter method and the conformal symmetry to do the computations, but the procedure is more complicated and the corresponding analytic structures, both in coordinate and in the coupling spaces, are considerably richer, due to a more "dynamical" nature of the involved Feynman graphs. A new phenomenon presented in the correlators of the full theory is the non-perturbative behavior of certain individual OPE data -anomalous dimensions and structure constants of exchange operators, in the weak coupling limit. But the perturbative behavior of the four point correlator is restored in the sum over all OPE terms. The equations on the anomalous dimensions, obtained from the pole structure of integrands in spectral decomposition of these correlators, appear to have a few interesting singularities in the coupling space, whose physical significance for the theory is left to understand. We also demonstrate the relevance of the double-trace terms for the correct Feynman graph interpretation of our results obtained via Bethe-Salpeter conformally symmetric procedure.
To get a further insight to these intereting chiral CFTs, we have to compute more complicated correlation functions, involving the exchanged operators of higher R-charges, such as trφ L 1 , or even more complicated multi-magnon operators. For the moment, only the exact anomalous dimensions of L = 3 case of such operators and of some related operators with the same R-charge have been computed for bi-scalar fishnet CFT in [33] 14 via the double-scaling limit of the QSC equations. Similar results on L = 4, 5 and magnon operators will be reported in [23]. Not much is yet done in this direction for the full χCFT, apart from the ABA approach of [19] to long operators L 1 and the one-loop study of [30], as well as the results of the current paper on the shortest exchange operators. As concerns the study of the structure constants, the first all-loop results for multi-magnon operators in bi-scalar fishnet CFT have been obtained in the very recent paper [54]. The generalization

JHEP06(2019)078
to four-point functions and to more complicated operators, and to the full χCFT, will necessitate a considerable new insight into integrability properties of these models.
The generalization of bi-scalar fishnet CFT to any dimension D [21] poses a natural question whether the D-dimensional generalization of the full 3-coupling χCFT exists. A related question: can we generalize the Basso-Dixon type fishnet integrals -the four-point single-trace correlators of scalar fields in bi-scalar CFT, explicitly computed in D = 4 [55] and D = 2 [56], to the case of dynamical fishnets of χCFT?
It would be also interesting to understand the behavior of large Feynman graphs in the full χCFT, in-line with the early results of [20] and the recent observations of [57] for the fishnet reduction of the χCFT. In particular, if the σ-model interpretation of the latter paper can be generalized to the full χCFT, it could be a big step in the explicit construction of the AdS dual of this chiral CFT, if such one exists at all after the double scaling limit of γ-deformed N = 4 SYM theory.
As a final comment: it would be interesting to find a realisation of these non-unitary theories in physical systems, if not in the fundamental quantum field theory (at least as an effective theory) than my be for certain statistical-mechanical and, presumably nonequilibrium, condensed matter models. The beautiful mathematical structures behind the χCFT promises more of such applications in the future.

JHEP06(2019)078
Remarkably, the null vector n µ in 4 dimensions can be realized in terms of Pauli matrices and two 2-spinors as n µ = λ σ µλ or n µ =η σ µ η. In such a context two S-tensor propagators can be contracted to form a tensorial loop as follows: where A µν ≡ (λσ µσν η) .
Under the condition α + β + γ = 4 + S we can write 3) which for S = 0 reduces to (C.1) and for S = 1 is equivalent to (C.2). In the same formalism we can encode the tensor structure of (3.6), that is n µ , allowing the application of formula (C.3) to the direct computation of (4.14).

D Eigenvalue of the fermionic graph-building operatorĤ t
In this section, we compute the relevant part of the fermionic eigenvalueĤ (n) F defined in (5.15) which is involved in the computation of the four-point function G φ 1 φ † 2 (5.1). The eigenfunction of the fermionic operator is determined, as usual, by the conformal representation of external fields, thus it corresponds to (3.5). As shown in section 5.2, the action of the sub-kernelsĤ t ,Ĥ r andĤ b transforms the conformal triangle to a new one, adding a factor / x/x 2 , with chirality depending on the number n of sub-kernels we apply. Moreover, since the spectrum ofĤ b is the same as of the bosonic operator (5.26) and the spectrum ofĤ r is expressed through the one forĤ t (5.24), the latter is the only operator we have to study.

(E.1)
This results in the following extra contribution to (3.27): (E.2) In addition to that, the measure factor 1/c 2 (∆, S) develops poles at ∆ = S + 3 + k, k = 0, 1, 2, . . . . The corresponding contribution can be expressed as 3) The overall contribution of these terms is the sum over all non-negative integers S, k of the generic term A possible vanishing condition for the full contribution is then It is easy to check that plugging (E.6) into (E.5), one is left with the condition (E.7), which means that (E.6) together with (E.7) are a sufficient condition for (E.5). To prove these equations hold, we notice first of all that at ∆ 1 = ∆ 2 = 1 (E.1) vanishes at odd n, so it would be sufficient to prove (E.6), (E.7) at even k ∈ 2N. Moreover, equation (E.6) has been checked in [24], where it was enough to state the cancelation of spurious poles in Tr[φ 2 1 ] sector. Let us verify the second condition (E.7) at even integer k. Starting from the sector Tr[φ 2 1 ] it is equivalent toh f 3+S+k,S −h f 3+S,S+k = 0 , (E.8)

F Operator mixing and logarithmic multiplet
In both sectors Tr[φ 2 1 ] and Tr[φ 1 φ † 2 ] of our theory the exchanged physical operators in the OPE s-channel of the 4-point correlators under analysis present mixing. Namely, due to the wide matter content of the theory, the renormalized operators are not just rescaled and normal-ordered monomials of elementary fields and derivatives, but linear couplingdependant combination of several such terms which share the same symmetries. Concretely, in our theory we deal with single trace primary operators as hermitian conjugation of (2.2) prevent to come to similar conclusions for χCFT theory. In particular, performing the planar limit can lead to "one-way" transitions and the correspondent mixing matrix can be only brought into Jordan canonical form, e.g. for the mixing of four primaries: The matrix (F.4) contains a 2 × 2 Jordan block, together with two diagonal terms γ 3 and γ 4 , corresponding to two renormalized operators with such anomalous dimensions. The physical interpretation of Jordan blocks leads to the formulation of logarithmic CFT (see [31,61] where ∆ 0 is the bare dimension of O i operators, and µ the energy scale. This phenomenon, the presence of log-multiplets in χCFT has first be noticed by J.Caetano [32] for its bi-scalar reduction (2.6) and some examples of its occurance in the context of Fishnet CFT have been presented in [33]. Despite such logarithmic operators appear in our theory, we are mostly interested in selecting the non-logarithmic ones: indeed these are the only exchanged in the OPE of the correlators under study, as the solutions of spectral equations (4.26) and (5.32) correspond to non-protected operators ∆(ξ) = ∆ 0 .
F.1 tr[φ 2 1 ] sector This first sector is characterized by the Cartan R-charge of two φ 1 fields (2, 0, 0). The equation (4.26) shows physical solutions for every even twist. In particular there is only one solution at twist-2 (4. 29), and two at twist-4 (4.30) and higher (4.31) both for spin S = 0 and S > 0. The twist-2 solution is easily interpreted as the scaling dimension of tr[φ 1 (n · ∂) S φ 1 ] + permutations S = 0, 2, . . . , (F.5) indeed for any S there is no other twist-2 conformal primary with charge (2, 0, 0). On the other hand for ∆ 0 − S = t ≥ 4 we can list several primaries with the right set of Cartan's charges. Let us concentrate on the scalar case S = 0 of twist four; we find 9 scalar conformal primaries which have the right set of charges

JHEP06(2019)078
As also the structure of (4.1) shows, this sector is fully described in terms of the χ 0 CFT, thus the mixing transitions are realized by the vertices of (2.4). At any coupling O 1 shows no planar transitions, and we deal with a set of 8 conformal primaries which at non-zero couplings ξ 2 , ξ 3 mix among themselves. This fact is apparently in contrast with the presence of only two twist-4 exchanged operators in the OPE expansion of section 4.5, and can be explained with the arising of logarithmic multiplets of operators with ∆ = ∆ 0 = 4, not being solutions of (4.26). Indeed, for instance, the following planar transitions can happen respectively starting from order ξ 2 and ξ 3 , while they lack the hermitian conjugate due to chirality of (2.4). One can actually check that there is no conjugate transition to (F.6) at any order. This suggest that matrix γ ij won't be diagonalizable and presents Jordan blocks in its canonical form, i.e. logarithmic operators.
In this sector the Cartan's U(1) ⊗3 charge is the difference between those of φ 1 and φ 2 , that is (1, −1, 0). The solutions of equation (5.32) are three in the scalar case S = 0, one of twist-2 (5.38) and two of twist-4 (5.39). Interestingly, for tensors S > 0 there is an additional solution of twist-2 (5.35). The explanation of such difference is that there is a unique way to realize a scalar ∆ 0 = 2 operator with (1, −1, 0) charge, that is tr[φ 1 φ † 2 ], while for higher rank tensors S > 0 two operators can mix up. Starting from S = 1 we find that in addition to Tr[φ 1 (n · ∂)φ † 2 ], the operator Tr[ψ 1 (n · σ)ψ 2 ] = −Tr[ψ 2 (n ·σ)ψ 1 ] satisfy the constraint of quantum numbers and behave as an SO(4) vector. Therefore operator mixing happens among two conformal primaries and O 2 = tr[ψ 1 (n · σ)(n · ∂) S−1ψ 2 ] , (F.7) at any S ∈ 2N + . Indeed, being interested in symmetric traceless tensors (type I [44]), no other bilinear covariant made of fermions can enter the game at higher S. The allowed planar transitions are: together with the "diagonal" terms O i (x)O † i (0) , for i = 1, 2. Dealing with single trace operators made of two elementary fields the effect of chiral interaction is cancelled by the trace cyclicity. The resulting transitions are symmetric and the anomalous dimension matrix is diagonalizable. Recalling the couplings redefinition (5.33) and (5.37), its first perturbative order ξ 2 should coincide with: (F.9)