On the Higher Spin Spectrum of Chern-Simons Theory coupled to Fermions in the Large Flavour Limit

In this note, we compute the higher spin spectrum of $U(M)_k$ Chern-Simons theory coupled to $N$ flavours of fundamental fermions, in the limit $N\gg M$ with the 't Hooft coupling $\lambda_M =\frac{N}{k_m}$ held fixed, to order $M/N$. This theory possesses a slightly broken higher spin symmetry, and may be of interest from the perspective of higher-spin and non-supersymmetric holography. We find that anomalous dimensions of the higher spin currents achieve a finite value at strong coupling $\lambda_M \rightarrow \infty$, which grows with spin as $\log s$ for large $s$, as expected for gauge theories.


Introduction and Summary of Results
U (M ) k Chern Simons theory coupled to N flavours of fundamental fermions arises as a simple limit of U (N ) k N × U (M ) k M Chern-Simons theory coupled to bifundamental fermions. This theory is a generalization of the theory studied [1], and is of interest from the perspective of holography and higher-spin gauge theory, as we describe below. The theory, and its cousin obtained by replacing the fermions with critical bosons has also been used in condensed matter physics as a calculable model to estimate critical exponents in fractional quantum Hall transitions (e.g., [2][3][4]).

Motivation and Bifundamental Chern-Simons Theories
As emphasized in [1], Chern-Simons theories coupled to matter are known to provide a large class of lines of non-supersymmetric conformal fixed points, in the large N limit. Such lines are rare (or non-existent) in higher dimensions (e.g., [34]) and are therefore of great interest.
The specific example that we will focus our attention on in this paper is U (N ) k 1 × U (M ) k 2 Chern-Simons theory coupled to bifundamental fermions [35]. This family of interacting theories is conformal in perturbation theory because the Chern-Simons levels k 1 and k 2 must be integers (or half-integers, depending on the number of fermion flavours), and therefore cannot run under RG flow. Because the scaling dimension of the fermion in three dimensions is 1, there are simply no other marginal or relevant interactions one can write down (other than a mass term, which can easily be tuned to zero order-by-order in perturbation theory in any given renormalization scheme).
Let us assume N ≥ M . Then a simple 't Hooft-like large N limit that one can study is one in which α = M N , λ N = N k N and λ M = N k M are held fixed. We will be using a dimensional regularization scheme to define the Chern-Simons level k, so |k N | > N and In the 't Hooft limit described above, the parameters λ N , λ M and α are effectively continuous, so we have a three-parameter family of non-supersymmetric interacting conformal field theories. Let us study the theory in an M/N expansion, as discussed in [36] for the N = 6 ABJ theory [37,38]. (See also [39] for a recent generalization of [40] to N = 5 supersymmetry.) When α = M/N 1, the theory is essentially the vector model studied in [1], which is exactly solvable in the large N limit. (See [41][42][43][44][45][46][47][48][48][49][50][51][52][53] for some all-orders results in Chern-Simons vector models.) The single-trace primary operators of this theory consist of only a single tower of spin s operators, one for each spin s. As argued in [1,54], in this limit the anomalous dimensions of all single trace higher-spin primary operators vanish, i.e. τ s = 1 + O M N . A bulk dual description must must therefore be a higher-spin gauge theory, and these operators would correspond to a tower of massless higher-spin gauge fields. Indeed, this higher-spin/vector-model duality [55,56] has been well tested when M = 1, [57][58][59][60][61][62][63][64] and we expect these results can be straightforwardly generalized to the bifundamental case in the limit M/N 1, if one considers higher-spin gauge fields with U (M ) Chan Paton factors, as discussed for ABJ theory in [40].
When M/N = 1, the theory is not a vector model, and instead can be thought of as a large N theory with matrix degrees of freedom, completely analogous to ABJM theory [37] or N = 4 SYM [65]. Indeed, the theory can be thought of as a simple non-supersymmetric generalization of the ABJ(M) [37,38] theory, that retains conformal invariance despite the lack of supersymmetry.
It is natural to ask whether the theory has a holographic dual at strong coupling when M/N ∼ 1. We do not, at present, have a concrete way of testing any answer (see [35] for suggestions) to this question since we are unable to perform any strong-coupling calculations when M/N ∼ 1. However, as a preliminary step, we can calculate corrections to the scaling dimensions of higher-spin currents at strong coupling as a power series in M/N . In this note we present some results for the first-order M/N corrections to the scaling dimensions of the single-trace, higher spin primaries which are twist one at zero coupling. The scaling dimensions depend on two parameters, λ M and λ N , and take the following form, expressed in terms of the twist τ = ∆ − s,: If we set λ M = 0, then δ s (0, λ N ) is simply M times the 1/N anomalous dimension of j s computed in a Chern-Simons vector model in [66], which we reproduce here:  (2)) + −38s 4 +24s 3 +34s 2 −24s−32 , for even s ,  The form (1.2) follows from the analysis based on slightly-broken higher spin symmetry given in [67,68], and the spin-dependent numerical coefficients were determined using techniques developed in [9,33] as well as a direct two-loop Feynman diagram computation. When s is large, a s ∼ 8 3π 2 and b s 2 π 2 log s. Note that, as mentioned above, |λ N | < 1, so the logarthmic growth of anomalous dimension on spin is only present for intermediate values of the coupling.
In this paper, we consider the higher-spin spectrum in a different limit, namely λ N = 0 and calculate the first-order M/N spectrum of higher spin operators to all orders in λ M . This limit can be essentially thought of as a large flavour limit, with N playing the role of a large number of flavours. (However, the flavour symmetry is gauged, so the only gaugeinvariant operators are singlets of the flavour symmetry). Below, we refer to the theory with λ N = 0 as the "large flavour" theory, and the theory with λ M = 0 as the "large colour" theory.

Summary of Results
Our results are as follows: , for even s, Here H s are harmonic numbers. In section 3, we argue that the λ M -dependence of these anomalous dimensions follows from the planar three-point functions and the higher spin symmetry, up to spin dependent coefficients. We then fix these spin dependent coefficients in the direct diagrammatic calculation presented in section 4.
The anomalous dimensions of higher-spin currents in this limit acquire finite values as λ M → ∞, which increase with spin as log s for large spin: The logarithmic dependence on spin is a characteristic feature of gauge theory that is expected from general arguments [69]. In section 5 we present the large spin expansion of this spectrum in a bit more detail, and discuss to what extent these results follow from a more general analysis of expectations for the large-spin spectrum of conformal field theories given in [30,32].

M/N Gauge-Propagator
See [35] for detailed description of the theory we study and our conventions.
Here we are working in the limit λ N = 0, and calculating the M/N correction to the anomalous dimension to all orders in λ M . In this limit, the gauge propagator is suppressed by a factor of 1/N , but it receives an infinite series of self-energy corrections to first order in M/N , which are given by: The one-loop self energy of the gauge field, Σ µν is given by Let D µν be the free gauge propagator and let G µν be the gauge propagator with an infinite series of self-energy corrections. Then, the matrices G and D are related by In light-cone gauge, which we will use throughout the paper, the corrected gauge propagator is Explicitly, 3 Analysis based on slightly-broken higher-spin symmetry

Higher Spin Currents
The higher-spin currents j s (x, z) are given by the following explicit expressions [1,66]: Here, and in what follows, z is a null polarization vector satisfying z 2 = 0 (see e.g., [9,70]), and we are suppressing color-indices. These expressions are defined for the free theory. In the interacting theory, derivatives are promoted to covariant derivatives. It was argued in [1] that these operators, along with the scalar primary j 0 =ψψ, are the only single-trace primary operators in the Chern-Simons vector model. In the bifundamental theory, this remains true, if we define "single-trace" operators to mean bilinear operators constructed from a single contraction of the U (N ) color index, which is natural in the limit N large and N M [40]. In momentum space, these operators can be written as: We define the free vertex as: When q = 0, this simplifies to The anomalous dimension, δ s = τ s − 1, of j s is related to the logarithmic divergence of

General Form of the Anomalous Dimensions
Let us first derive expressions for the general form of the anomalous dimensions using the slightly-broken higher spin symmetry of the theory. As argued in [41,66], we can determine the 1/N anomalous dimension of j s from the leading order (planar) parity-violating threepoint functions j s 1 j s 2 j s , outside the 'triangle inequality', (i.e., with s 1 + s 2 < s) [41,71]. Let us review this briefly. For simplicity, let us take M = 1 in what follows (although all these expressions can be generalized to M = 0 as long as M N ).
where [j s 1 ][j s 2 ] denotes the particular combination of descendents of j s 1 and j s 2 given in [66]. Note that the exact form of [j s 1 ][j s 2 ] is uniquely determined by the requirement that ∂ · j s is a conformal primary. Also, note that the only spins that can appear on the right hand side of (3.7) are those satisfying s 1 + s 2 < s. Then, via the state operator correspondence, (see, e.g., [1,41]) This implies that, where α s 1 ,s 2 ,s is a purely numerical coefficient that can be computed using the particular combination of descendents of j s 1 and j s 2 represented by [j s 1 ][j s 2 ].
We can determine the C s 1 ,s 2 ,s (λ M ) from the planar three-point functions as follows which implies Here ∼ denotes equality up to a numerical coefficient that again depends on the details of the form of [j Although this is a theory with a slightly-broken higher spin symmetry its planar threepoint functions have not been computed in [41], which applies only to theories containing only even spins (e.g., O(N ) theories with Majorana fermions). However, it is straightforward to compute them in momentum space directly. Loops of the gauge field are suppressed by 1/N in the large flavour limit, so we only consider tree-level diagrams. It is easy to see that nearly all three point functions are those of the theory of free fermions -the only exceptions are three-point functions involving the spin 1 current j 1 , pictured in figure 1 and 2. The U (M ) Chern-Simons gauge field effectively acts as a double trace deformation of the schematic form j 1 j 1 in the action, once integrated out.  The non-vanishing parity-odd three point functions, outside the triangle inequality, are only j 1 j s j s and j 1 j 1 j s , where s, s = 1, pictured in Figures 1 and 2. The second three-point function is only non-zero when s is even [1]. Using the the corrected gauge propagator, we can directly read off the λ M -dependence of these three point functions: We normalize two-point functions as j s j s ∼ N , but the two-point function of j 1 is given by as can be seen from Figure 3. This implies that  The anomalous dimensions are then given by: where c (1) s and c (2) s are spin-dependent numerical coefficients, and c (2) s is non-zero only if s is even because j 1 j 1 j s vanishes when s is odd. When M = 1, one can show that equation (3.17) is multiplied by M , as long as M N . (However, this is not true at higher orders, as M 2 /N 2 corrections will differ from 1/N 2 corrections.) In the section below we directly evaluate the logarithmic divergences of the relevant Feynman diagrams to compute the anomalous dimension and verify that it is of the form (3.17). From the calculation we are also able to read off the coefficients c (1) and c (2) , and find c (2) = 0 when s is odd, as expected. The results are: We conclude by noting that one could also use the classical equations of motion to determine C s 1 ,s 2 ,s to order λ 2 M . The coefficients C s 1 ,s 2 ,s in this large-flavour theory are a priori expected to be different from those in the large-colour theory [66], since the theories are different. For example, for s 1 , s 2 = 1, C s 1 ,s 2 ,s are zero for the large flavour theory but non-zero for the large colour theory. Then c (2) would be related to C 2 1,1,s , and c (1) would be related to the

Direct Feynman Diagram Calculation
The calculation of the anomalous dimension of the scalarψψ was carried out in [35]. Here, we consider the anomalous dimensions for currents with s > 0 which follows [35]. The anomalous dimension is given by where δ (1) s is the contribution from the self-energy diagram, δ (2) s is the contribution from the rainbow diagram, δ s is the contribution from the two three-point functions.

Self-Energy
where We find the contribution of self-energy correction to anomalous dimensions is: The rainbow correction is determined by the following integral

Rainbow Correction
which contributes to the logarithmic divergence via and we find: where, g(s) is (4.10) The sum of the two "three point-function" diagrams is given by the following integral, when s is even,

Other Contributions
The two diagrams cancel if s is odd. Evaluating the diagrams for even s, we obtain: The contribution to the anomalous dimension is read off from the divergence via: (4.14) From here, The anomalous dimension for even s is: This is always positive. One can observe that it vanishes for s = 2.
For odd s, the third diagram does not contribute and the anomalous dimension is: One can observe that it vanishes for s = 1.
For large spin, we have:  where, H s− 1 2 is a harmonic number. At large coupling the anomalous dimensions approach a constant value, which increases logarithmically with s. The strong coupling limit for both even and odd spin is plotted in Figure 8.

Large Spin Expansion
In this section, we briefly comment on the large spin expansion of our results and the results of [66].
In the large spin limit, our results can be written as: for odd s,   We presented the anomalous dimension in this way to draw an analogy with equations (61) and (63) of [30], which give the large spin expansion of the large N anomalous dimensions of non-singlet and singlet operators in the critical O(N ) model [7,8].
The first term in each equation is the s → ∞ limit of the anomalous dimension, which contains a log s, unlike the critical O(N ) model. Hence the anomalous dimension of ψ is undefined, which is expected as it is not a gauge invariant operator.
Looking at the subsequent terms in equations (5.1) and (5.2), it appears as if, just as in the critical O(N ) model, for odd-spins (which resemble non-singlet operators in the critical O(N ) model), the anomalous dimensions are due to an intermediate operator of twist two (i.e.,ψψ), while for even spins (which resemble singlet operators in the critical O(N ) model) we have two series -the first series corresponding to the twist-two operator and the second series corresponding to twist-one operators (such as j s ). There may be some ambiguity in separating these two series, and these observations are only suggestive at present because the the analysis of [30] does not immediately apply here.

Conformal Spin
It is also interesting to express these results in terms of the conformal spin J [30], defined as a function of spin s and twist τ = ∆ − s as: In [30], it was argued that the expansion of anomalous dimensions around large conformal spin should contain only even powers of 1/J, possibly with extra log J coefficients. Due to the fact that the field ψ is not a gauge invariant operator, and the dimensionality of space-time is odd, the results of [30] do not immediately apply here. However, as will be described elsewhere, applying the ideas of [32,72,73], it is possible to refine the arguments of [30] to show that this result holds to order λ 2 M , but should not be expected to hold at higher orders. 1 Expanded in powers of M/N , the conformal spin is given by Because the leading order anomalous dimensions vanish in our theory, τ = 1 + O( M N ) and J 0 is given by the free value:  contains only even powers of 1/J in its expansion, apart from the initial log J. Hence the entire expression has an expansion in even powers of 1/J, to all orders in λ M . For even spins, δ s (J) = M N 8λ 2 M (−11 + 6γ + 12 log(2)) 3 π 2 λ 2 M + 64 Here, we see that the last piece would contain even as well as odd powers of 1/J. This corresponds to an order λ 4 M correction, so at order λ 2 M , we again get an expansion with only even powers of 1/J. The odd powers of 1/J originate from a term proportional to √ 4J 2 + 1, and could be interpreted as originating from the exchange of a tower of twist-one (almost-conserved) operators. Such a term is also present in the singlet large N anomalous dimensions of the critical O(N ) model in d = 3, given in equation (64) of [30].

Discussion
In this paper, we studied the M/N anomalous dimensions in a bifundamental U (N )×U (M ) Chern-Simons theory, δ s (λ M , λ N ), in the special case λ N = 0. These results compliment the results of [66], which effectively determine the anomalous dimensions in the case λ M = 0. It should also be possible to determine the M/N anomalous dimensions in the bifundamental theory to all orders in both λ M and λ N by extending the arguments given in section 3, however, the calculation would be substantially more involved. We hope to address this problem in the near future.
Chern-Simons theories coupled to matter have also played an important role in the study of fractional quantum Hall transitions, e.g., [2][3][4]. In particular, [3] includes a study of U (1)-Chern Simons theory coupled to N flavours of fermions as a model of a Mott Insulator -quantum Hall transition, which, at the level of 1/N corrections also has the same spectrum as the U (M ) theory we consider. The anomalous dimension of the scalar operator ψψ effectively determines the critical exponent ν associated with the transition, and so is of experimental interest. To date, the anomalous dimension ofψψ in the large colour theory is not known, and it would be interesting to compare the value of ν from this calculation to the existing large flavour calculation.
We are not aware if the anomalous dimension of the higher-spin currents have any similar physical interpretation or relevance to Quantum Hall physics, but it may be interesting to pursue this further. It might also be possible to study the spectrum at order 1/N 2 , which may perhaps provide better agreement with experiment (in which N = 1). The techniques involved would be similar to the recent computations of 1/N 2 spectrum of the critical O(N ) and Gross-Neveu models [13,15]. Correlation functions involving the stress tensor and conserved U (1) current could also be calculated in this theory possibly at finite temperature/chemical potential (e.g., [74]), and may correspond to physical observables. These should obey the bounds determined in [75][76][77][78].
It would also be of interest to consider scaling dimensions in the theory of U (M ) k Chern-Simons coupled to N critical bosons, when N is large. In [2], the anomalous dimension of the scalar operator j 0 was calculated in the U (1) version of this theory. It would be interesting to also compute the higher spin spectrum and compare it is to the spectrum of the fermionic theory studied here.
Chern-Simons theories coupled to matter also exhibit a beautiful bosonization duality [79] relating theories with fermionic matter to critical bosonic matter. This duality was discovered and has been supported via various large N computations [41][42][43][44][45][46][47][48][48][49][50], but it is now conjectured to hold at finite N as well [80][81][82][83][84], and has been of interest from the point of view of condensed matter (e.g., [85][86][87][88]). These dualities can be thought of as generalizations of the supersymmetric Giveon-Kutasov duality [89][90][91][92]. While these dualities cannot hold in the large flavour limit, they may perhaps be generalized to bifundamental Chern-Simons theories when both couplings are nonzero, as a strong/weak duality when λ N ↔ (1 − |λ N |). One form for these dualites was conjectured in [35], based on analogy with ABJ dualities [38], and could be tested by explicit computation when M N in the fermionic and critical bosonic theory. (Some perturbative calculations in the non-critical bifundamental bosonic theory appear in [93].) We hope to perform these tests in the near future.
Finally, we found the spectrum of anomalous dimensions appears to take a form similar to those of the large N critical O(N ) model, when expanded in a large spin expansion, and also expressed in terms of conformal spin, despite the additional logarithmic dependence on spin. It would be interesting if this behaviour could be better understood from the point of view of large spin perturbation theory, e.g., [30,32,72,73].