Correlation functions in N = 2 Supersymmetric vector matter Chern-Simons theory

We compute the two, three point function of the opearators in the spin zero multiplet of N = 2 Supersymmetric vector matter Chern-Simons theory at large N and at all orders of ’t Hooft coupling by solving the Schwinger-Dyson equation. Schwinger-Dyson method to compute four point function becomes extremely complicated and hence we use bootstrap method to solve for four point function of scaler operator J 0 = ψ̄ψ and J b 0 = φ̄φ. Interestingly, due to the fact that 〈J 0 J f 0 J b 0〉 is a contact term, the four point function of J 0 operator looks like that of free theory up to overall coupling constant dependent factors and up to some bulk AdS contact terms. On the other hand the Jb 0 four-point function receives an additional contribution compared to the free theory expression due to the J 0 exchange. Interestingly, double discontinuity of this single trace operator J 0 vanishes and hence it only contributes to AdS-contact term. E-mail: karthikin@tauex.tau.ac.il, sachin.jain@iiserpune.ac.in, vinaymm@iiserpune.ac.in, abhishek.mehta@students.iiserpune.ac.in, pranjal.nayak@uky.edu, tarun sharma2@brown.edu, tarun.sharma@niser.ac.in ar X iv :1 90 7. 11 72 2v 1 [ he pth ] 2 6 Ju l 2 01 9


Abstract
We compute the two, three point function of the opearators in the spin zero multiplet of N = 2 Supersymmetric vector matter Chern-Simons theory at large N and at all orders of 't Hooft coupling by solving the Schwinger-Dyson equation. Schwinger-Dyson method to compute four point function becomes extremely complicated and hence we use bootstrap method to solve for four point function of scaler operator J f 0 =ψψ and J b 0 =φφ. Interestingly, due to the fact that J f 0 J f 0 J b 0 is a contact term, the four point function of J f 0 operator looks like that of free theory up to overall coupling constant dependent factors and up to some bulk AdS contact terms. On the other hand the J b 0 four-point function receives an additional contribution compared to the free theory expression due to the J f 0 exchange. Interestingly, double discontinuity of this single trace operator J f 0 vanishes and hence it only contributes to AdS-contact term.

Introduction
The perturbative technique to compute observables in quantum field theories involving Feynman diagrams is effective only when the coupling is weak and breaks down in the strong coupling regime. In the past few decades, various strong-weak dualities have been discovered proving to be extremely useful in understanding some of the most interesting -1 -non-perturbative properties of strongly coupled quantum field theories. One such class of dualities which have been studied extensively in recent times are the Bosonization dualities in Chern Simons gauge theories coupled to fundamental matter at large N [1][2][3][4][5][6]. Though, one of the main indications for these dualities initially came from their holographic duality with Vasiliev theories in AdS 4 [7][8][9][10][11][12], by now there exists a plethora of evidence for these dualities coming from exact computations of correlation functions, thermal partition functions, scattering amplitudes and RG flow analysis relating these theories to known supersymmetric dualities [3][4][5].
Primary example of these bosonization dualities are those among the quasi-fermionic (critical bosonic and regular fermionic theory) and the quasi-bosonic theories (regular bosonic and critical fermionic theory) 1 . A particularly interesting case of these dualities is present in the N = 2 supersymmetric U (N ) Chern Simons theory coupled to a single fundamental chiral multiplet. This theory exhibits a strong-weak self duality [34][35][36] generalizing the well known Giveon-Kutosov duality [34,37]. The self duality of this supersymmetric theory serves as a parent duality for the non-supersymmetric bosonization dualities mentioned above since they can obtained from the supersymmetric theory via RG flows seeded by mass deformations [38,39]. Taking hints from the supersymmetric dualities and the Level-Rank duality of pure Chern Simons theory, finite N extensions for the non supersymmetric dualities have also been proposed [40][41][42][43][44][45][46][47][48][49]. These theories were also investigated recently in presence of background magnetic field [50].
In this article, we will focus our attention on the N = 2 theory. Various large N computations in this theory show remarkable features which are absent in the non supersymmetric couterparts. For example, the all loop 2 → 2 scattering amplitude is tree-level exact except in anyonic channel [26] and it was shown that these amplitudes are also invariant under Dual superconformal symmetry [26,29]. [28] further showed that the tree level m → n scattering amplitudes in this theory can be constructed using the BCFW recursions relations.
Although many interesting non supersymmetric physical observables, as mentioned above, are amenable to direct exact computations by solving corresponding Dyson-Schwinger equations, the computation of 4-point correlation function of even the simplest of single trace operators, namely the scalar operatorsφφ andψψ, appears prohibitively difficult 2 to compute via this direct approach. Given the remarkable simplicity of the results for other known observables one expects the 4-point functions in this theory to also have a simple structure. In the present article, our main goal will be to determine the exact 2, 3 and leading connected 4-point correlation functions of scalar operators in this supersymmetric theory.
For the quasi-bosonic and quasi-fermionic theories mentioned above, [51,52] recently determined the 4-point function of scalar operator using recently developed ideas from conformal Bootstrap. In particular, one of the central of the central objects used in [51] is the double discontinuity of the 4-point function which determines the coefficients in the 1 In the terminology of [3,4] 2 See appendix (D) for a discussion of our attempt.
-2 -OPE expansion of external operators via the Lorenzian inversion formula (LIF) discovered by Caron-Huot in [53]. The authors of [51] first demonstrated that for large-N CFTs the double discontinuity of the 4-point function of identical scalars determines the full 4-point function up to three AdS 4 contact Witten diagrams. The authors further showed that for the quasi-bosonic and quasi-fermionic theories the coefficients of these contact terms vanish. In the present work, we apply some of these ideas in conjunction with the self duality, to the case of scalar 4-point functions in our N = 2 theory.
Our article is structured as follows. In section 2, we describe the N = 2 theory of interest in this paper its operators spectrum in some detail. In section 3, we determine the scalar multiplet 2 and 3-point functions via a directly computation. In section 4, we determine the 4-point function of the bosonic and the fermionic scalar operators in this theory using the double discontinuity technique developed in [51]. Finally, in the section 5, we summarize our results and outline related open questions and future directions. In various appendix, we collect our notation and conventions, some technical details of the results in main text of the paper and briefly summarize our attempt at the direct computation of 4-point function.
Note added in the proof : While we were in the process of finishing up our article, we were informed about the related work [54] by the authors which has overlap with the results of our section 3.

N = theory and its Operator Spectrum
In this paper, we are interested in N = 2 U (N ) Chern-Simons theory coupled to single chiral multiplet, Φ ≡ (φ, ψ), in the fundamental representation of the gauge group. The position space Lagrangian for the theory is The theory above has two parameters : the rank of the gauge group, N , and the Chern-Simons level, κ, which is quantized to take only integer values [55]. κ −1 controls the strength of gauge interactions and the theory is perturbative for large values of κ at any finite N .
This theory is conjectured to be self-dual under a strong-weak type duality, [34]. In the 't Hooft like large N limit of interest in this paper, the duality transformation is Apart from the matching of many of the supersymmetric observables which can be computed at finite N and κ using supersymmetric localization techniques, recent exact computation of many non-supersymmetric observables, e.g. the thermal partition function, in the large N limit [1,2,5,25,30,38] has provided ample evidence for this conjectured duality.
The theory is quantum mechanically (super) conformal for all values of κ and N . In the 't Hooft limit, one can focus on the single trace superconformal primary operator spectrum of the theory. Though our theory has N = 2 superconformal symmetry, in this paper we will work in the N = 1 superspace formulation to allow us to use the relevant results of [26] for our computations. In the N = 1 language, the operators spectrum of the theory consists of a set of supercurrent operators [56] which are written in terms of the superfields, using the auxiliary commuting polarisation spinors, λ α i , which keep track of the spin; and ∇ α i are the standard supersymmetry invariant gauge-covariant derivatives. Their action on the matter superfields of our theory is given by, The explicit expressions for the first few spin-s currents are, In the free limit of the theory i.e. λ → 0, each of these supercurrents, J (s) with s = 0, satisfies the conservation equation and constitutes two component conserved current operators {J (s) , J (s+ 1 2 ) } in its θ expansion [56]. At finite λ, the conservation equation (2.8) is violated at order 1 N by double trace operators [1,56].
In this article, we are interested in the scalar operator J (0) (θ, x). There is no conservation equation associated with this operator and it constitutes 2 scalar and 1 spin half operator as follows In the subsequent sections, we compute the 2 and 3-point functions of the J (0) operator and two component of the 4-point function.

Correlation functions
In this section, we compute the two and three point correlation function of the J 0 (θ, p) operator. Two of the main ingredients for these computations are the exact propagator (3.2) and the renormalized four point vertex for the fundamental superfield Φ(θ, p) (ν 4 in 3.3). These were computed in [26] for a more general class of theories with N = 1 supersymmetry which can be thought of as one parameter 3 deformation of the N = 2 theory of interest in this paper. Below, we list these results for our N = 2 theory, conveniently stated in term of the exact quantum effective action The quadratic part of the effective action receives no quantum corrections at large N in the N = 2 theory. The propagator is thus tree level exact and given by The quartic superspace vertex, ν 4 , does receive quantum corrections and takes the following form ν 4 (θ 1 , θ 2 , θ 3 , θ 4 ; p, q, k) = e 1 4 X.(p.X 12 +q.X 13 +k.X 43 ) F 4 (X 12 , X 13 , X 43 ; p, q, k), 3 Quartic superpotential term : − π(ω−1) κ d 3 x d 2 θ (ΦΦ) 2 . ω = 1 is the N = 2 point.
-5 - The overall exponential factor is determined by supersymmetric Ward identity (3.7), while the coefficient functions A, C and D require explicit computation and given by [26] A(p, q, k) = − 2πi κ e 2iλ tan −1 2ks Note that the vertex ν 4 was computed in a special momentum configuration, namely while the momenta p and k are arbitrary 4 . For this reason, our computation of correlation functions will also be restricted configuration in which the momentum of J 0 operators are restricted to lie only in the 3-direction. Diagrammatically, the exact four point vertex will be represented as in Figure 1.

Constraints on correlation functions from supersymmetry
To begin with, let us study the constraints on an arbitrary correlation function due to supersymmetry. As stated earlier, although our theory has N = 2 supersymmetry, we will be working in N = 1 superspace following [26]. A general n-point correlation function of N = 1 scalar superfield is constrained by supersymmetry and translation invariance to take the following form [26] where θ ij = θ i − θ j . The δ 3 ( i p i ) follows from translation invariance while the overall Grassmann exponential factor follows from invariance under N = 1 supersymmetry. Note 4 We refer the reader to appendix A for conventions for labelling momenta.
-6 -that the function F n above only depend on the differences of the Grassmann coordinates. Following [26], the form is easily derived as follows In the last line above, we used the momentum conservation to replace p n with n−1 i=1 (−p i ). Defining the sum and the difference of Grassmann variables as The factorised form in (3.6) follows as the solution to last equation in (3.7).

J 0 -vertex
Before proceeding to the computation of correlation functions, it would be useful to compute an intermediate quantity, the J 0 -vertex. It is defined by stripping of the propagators from J 0 ΦΦ as follows and satisfies the same Ward identity as a three point function (3.6). The vertex receives contribution both from the free propagation of the fundamental field as well as from the interaction vertices in the theory. The free part vertex is simply proportional to the momentum and the Grassmannian δ-functions while the interacting part of the vertex can be computed from the exact ν 4 vertex. Figure (2) shows the relevant diagrams. Explicit computation of the above integral, with constraint p + = p − = 0 following from the (3.5), leads to the following result for the full J 0 -vertex factor The J 0 -vertex computed above will be useful in further computations of 2 and 3 point functions of the J 0 operator.

3.3
J 0 J 0 correlation function The 2 point function can be straightforwardly computed from the J 0 -vertex determined in the previous section by combining the exact vertex on one side with the free vertex on the other side. Figure (3) shows the relevant diagram which leads to the following integral for the two point function Again, the collinear constraint (3.5) restricts the momenta p and s to lie in 3-direction.
Computing the integrals with this constraint leads to the following result The result can be straightforwardly generalized for arbitrary external momenta to give The non vanishing component correlators can easily be read off to give Let us compare the above two-point functions with the corresponding two-point functions in the regular fermionic and regular bosonic theories studied in [5] and [30] respectively.
Note that as opposed to the regular bosonic and regular fermionic theories studied in [5] and [30], the λ dependence of the two-point function of J b 0 and J f 0 operators is the same as that of the higher spin currents in the non-supersymmetric cases. Further, using the double trace factorization argument of [20] relating the two-point function of current operators in the supersymmetric and the above mentioned non-supersymmetric theories, we know that the two-point function of all the current operators in our supersymmetric theory is exactly the same as those of the corresponding regular boson/fermion theory. Thus, we see that in our theory the two-point function of scalar operators is the same as that for the higher spin current operators. The reason for this is supersymmetry. Though we are working in N = 1 superspace language, our theory has underlying N = 2 supersymmetry under which the scalar operators J b 0 , J f 0 belong to the same supersymmetry multiplet as There are two such diagrams, as shown above, which turn out to be equal.
the spin 1 conserved current and thus the two-point function of the two are thus related by supersymmetry.
The full 3-point function can be constructed by combining three J 0 vertices with exact propagators. There are two such diagrams shown in figure (4). Each of these two diagrams can easily be shown to be cyclically symmetric and related to each other by pair-exchange of any two J 0 insertions. An explicit computation of the diagram shows that each of the diagrams is completely symmetric (cyclic as well as under pair-exchange) by itself and the two diagrams are equal. The full 3 point function is then just twice the contribution of the first diagram which we write down below.
The overall factor of 2 in the above equation is from the sum over two triangle diagrams in figure (4) which turn out to be equal while the factor of N results from index contractions. Explicit computation of the above integrals in the collinear limit of the external momenta gives the following result The overall factor of 2 in the expression of F 3 above is from the sum over two triangle diagrams which turn out to be equal. The coefficients {A i } are given by The non vanishing components of the three point functions can easily be extracted from (3.17) and (3.18) to be Notice that in the above result for 3 point functions, two different functional forms of λ dependences appear, namely sin πλ πλ and sin 2 πλ πλ . The two of them differ in a crucial way. The first one has a finite λ → 0 limit and is invariant under parity under which λ is odd. The second is odd under parity and vanishes in λ → 0 limit. This result thus provides some support for the conjecture made in [56] that the three-point functions in N = 1 superconformal theories with higher spin symmetry have exactly one parity even and one parity odd structure. The results (3.14) and (3.18) for the 2 and 3-point are clearly invariant under the duality transformation (2.3).

Four point functions
In the previous section, we evaluated the 3-point functions involving the J 0 operator in the N = 2 supersymmetric theory by computing the required vertex. However, the direct computation of the four-point function of J 0 operator following the same technique has proven to be intractable in our attempt till now. We describe our attempt to evaluate this four-point function in momentum space through the required vertices in the Appendix (D).
In this section, we determine the four-point correlators of the J b 0 and J f 0 operators using a novel method developed in [51], which we briefly review below. Note that we will be evaluating the 4-point correlation function in the position space as in [51].
Consider the position space four-point correlator of the identical external operators with conformal dimensions ∆. The function A which is known as the reduced correlator is defined as follows Here, u, v are the standard cross-ratios: The conformal block expansion expressed in terms of the reduced correlator A(u, v) is given as where G ∆ k ,J k (u, v) is known as the conformal block corresponding to the operator O k with scaling dimension ∆ k and spin J k . In the supersymmetric four point functions of J 0 operators, the relevant exchanges are schematically shown below

Review of the double discontinuity technique
In [51], the authors determine the four-point correlation functions of the scalar operator in the non-supersymmetric scalar/fermion coupled to Chern Simons gauge field i.e. quasibosonic and quasi-fermionic theory respectively. In order to obtain the required four-point -12 -  functions, the authors utilize the inversion formula which relates the double discontinuity to the OPE coefficients [53]. The authors first prove an interesting theorem that in the large-N limit of a CF T d , the double discontinuity constrains the four-point correlator up to three contact terms in AdS d+1 . Suppose there are two solutions G 1 and G 2 to the crossing equation with the same double discontinuity then they are related by the contact interactions in the AdS as follows Furthermore, the authors showed 5 that for the four-point function of single trace scalar operator in Chern-Simons coupled fundamental scalar/fermion theories these AdS 4 contact terms do not contribute and hence the double discontinuity completely determines the fourpoint functions.
Consider the normalized three point functions of the operators O i (i = 1, 2, 3) which are defined as follows In [3,4,51], it was noticed that the square of this normalized coefficients in the quasifermionic theories (C 2 s,qf ) are related to that of a single free Majorana fermion (C 2 s,f f ) as whereÑ is related to the the rank of the gauge group N and coupling λ qf by, Note that the normalized coefficients of quasi-fermionic theory and free fermionic theory are proportional to each other as given in eq.(4.5). Hence, the double discontinuity of the scalar four point function in the free fermionic theory is same as that of the quasi-fermionic theories up to an overall factor which depends only on N and λ qf . On the other hand, the square of the normalized coefficients of the quasi-bosonic theories (C 2 s,qb ) are related to the theory of a free real boson (C 2 s,f b ) as follows whereÑ andλ are related to N and coupling λ qb as λ qb = tan πλ qb 2 . (4.10) Note that unlike the normalized coefficients of the quasi-fermionic theories, in the quasibosonic theories, the spin s = 0 and s = 0 coefficients given above have different factors in front of their free bosonic counterparts. In order to account for the second term on the RHS of eq.(4.8) one needs to add a conformal partial wave with spin-0 exchange which is given by the well knownD-function with the correct pre-factor [51]. We now proceed to employ this technique for the supersymmetric case.

Double discontinuity and the supersymmetric correlators
Here, we utilize the technique described above to compute the four-point correlators for spin-0 operators J b 0 and J f 0 in our supersymmetric theory. Since we are considering correlators of identical external operators 6 , only even spin operators will contribute to the block expansion. 6 Although we have all the three-point correlators required, we do not compute mixed correlators such The four point function of the J b 0 operators is expressed as follows Here, disc corresponds to the disconnected part given by while F (u, v) is given by In order to determine the double discontinuity, and hence the four-point functions in the supersymmetric case using the method described above, we first obtain all the normalized three-point functions coefficients defined in eq.(4.4). Utilizing the two and three point functions obtained in the previous section, the normalized coefficients for spins s = 0 operators that contribute to the J b 0 four point function are given as follows 7 On the other hand the normalized coefficients involving one of operators with non-zero spin (s > 0) are given by where C 2 s,f b is the corresponding normalized coefficient of three point function for the free bosonic theory. The derivations for the above relations are provided in the Appendix C.
Since, the above relations were computed in momentum-space, they must be converted to position-space to make our results useful. 8 Similarly, for s > 0 are given bỹ Note that we may re-express both the spin zero coefficients given by eq.(4.18) and eq.(4.19) as followsC 2(BBB) 0,susy where ∆ φ being the conformal dimension of the external operator. Notice that for ∆ = 2∆ φ + J + 2m, the double-discontinuity vanishes. Therefore, for the double-trace exchange, the double-discontinuity vanishes. That is why the OPE of single-trace operators are sufficient to construct a function that has a double-discontinutiy equal to the four-point correlator. However, notice that the single-trace exchange J F F 0 with quantum numbers (∆, J) = (2, 0) also vanish. Coincidently, the double-trace operator [J b 0 , J b 0 ] 0,0 also has the same quantum numbers. 9 By inspection, we can see that the function below has the right double-discontinuity (v, u) + 1 where, the function f f b (u, v) is the free bosonic part given by. 10 The contact terms are explicitly provided in E.18. Note that c 1 contains contribution from both single-trace and double-trace operators which we have separated in the following equation as a 1 andc 1 O is a single-trace operator. 10 Note that we may have used two separate tree-level φ 3 exchange Witten diagrams corresponding to ∆ = 1 and ∆ = 2 bulk exchange with arbitrary coefficients instead [58]. But Witten diagrams themselves admitting an expansion in contact terms would compound the problem. TheD-functions, therefore, represents the choice with the least number of contact terms and the right double-discontinuity.

-16 -
To determine a 1 we take the OPE limit. In the OPE limit, the conformal blocks go like 11 For (∆, J) = (2, 0) i.e. for J f 0 exchange, we have G 2,0 (u, v) ≈ u in the OPE limit. Since, we are interested in the single-trace operator J f 0 , hence, we have In the OPE limit, we have for φ 4 contact term By only looking at the single-trace contributions we obtain Now, we focus our attention to double-trace operators. Coefficientc 1 can now be determined by looking at the double-trace trace operator [J b 0 , J b 0 ] 0,0 . Since, (∆, J) = (2, 0) for the double-trace is same as that of the single-trace operator J f 0 , we use the same method to obtainc 1 Hence, we have determined the first coefficient of the AdS contact terms. The results are a little cumbersome and we report it in the appendix F. We leave the explicit computation of these ope coefficients for future work.
The four point function of J f 0 is given by the following expression where, disc denotes the disconnected piece given by We now proceed to determine the four-point function J f 0 utilizing the same technique. The normalized coefficients that are required in this case are given by 12 where C 2 s,f f is the normalized three point functions for free fermionic theory. Unlike the previous section, changing the OPE coefficients to position-space is redundant here as both sides of the eqaution change by the same factor. Once again these relations are derived in Appendix C. Note that the three point functions of the spin-0 exchanges given by C 2(F F F ) 0,susy and C 2(F F B) 0,susy are contact terms in this case which, therefore, may be set to zero. This implies that the above relation is trivially satisfied for the spin s = 0 case as the free fermionic coeffcient C 2 0,f f = 0. Hence, both the s = 0 and s = 0 coefficients in this case come with the same pre-factor. This implies that the function which has the correct double discontinuity is given by where f f f (u, v) is the free fermionic part given by However, in this case, the separation of the double-trace and the single trace contributions is redundant as there are no single-trace operators with (∆, J) = (4, 0). After a faithful implementation of the methodology of the previous section, we obtain This completes our analysis of the determination of the four point functions of J b 0 and J f 0 operators in the Chern Simons matter theories in the supersymmetric scenario. The contact terms are explicitly provided in E. 19. The results are a little cumbersome and we report it in the appendix F. We leave the explicit computation of these ope coefficients for future work.

Summary and Discussion
In this article, we have focused our attention on the N = 2 U (N ) Chern Simons theory coupled with a single fundamental chiral multiplet in the 't Hooft large N limit and presented the computations for the exact 2 and 3-point functions for the scalar supermultiplet.
The result are invariant the duality transformation (2.3) and can be seen as an independent confirmation of the duality. For the case of 4-point function, though we are not able to perform the direct computation for the full scalar supermultiplet, we are able to use the a combination of techniques from conformal bootstrap, factorization of 3-point functions via double trace interactions along with the self duality of our theory to write down two of the component of the full J 0 4-point function. Though this leave room for 3 undetermined coefficient, we could formally related these to 3 point function coefficients involving specific double trace operators. We plan to report on these in near future.
Though we have focused on the N = 2 theory in this paper, the approach used to compute the four point function is can be straightforwardly applied to the one parameter deformed N = 1 theory. These differ from our N = 2 theory only via a double trace term in N = 1 superspace 13 . The 2 and 3-point functions of the two theories can thus be related via the double trace type factorization also used in this paper.
The approach used in this paper, following [51], to compute the J b 0 and J f 0 4-point functions relies crucially on the fact that the double discontinuity of the 4-point function in the interacting theory is almost the same as that of the free theory. We could thus write down the full interacting 4-point function in term of the free 4-point function. For the case of mixed 4-point functions, e.g.
this approach is not directly useful since a free theory analogue of such mixed correlator is not available since bosons and fermions decouple from each other the mixed 4 point correlators vanish in λ → 0 limit. One approach that might be useful in this regard is to first study the single trace OPE coefficients in the N = 1 deformed theory (for general w) in λ = 0 limit. We expect this limit to be significantly simpler then N = 2 theory and one can compute not only the exact 2, 3 point functions (see e.g. [54]) but perhaps even the exact 4-point function (we expect it to be non vanishing for w = 0) of J 0 operators in this limit since the only interaction term present is a double trace term. If this is indeed turn out to be the case, one can compare the double discontinuity of mixed J b 0 , J f 0 correlators in N = 2 theory with this limit and see these are closely related in a similar way as in [51] and in this paper for the identical scalar 4-point function.
As we have noticed in this paper, the coefficients {c i } can be determined in term of the normalized 3-point function coefficients of specific double trace operator. An interesting property of the AdS 4 contact Witten diagrams is that their series expansions contain Log terms. This implies that the coefficients {c i } not only contribute to the OPE coefficients of double trace terms but also to their leading anomalous dimensions as well, but in a coordinated way. The absence of Log term in the free 4-point function along with the vanishing of these coefficients for quasi-bosonic and quasi-fermionic theories [51] thus means that these double trace operators in the leading large-N order do not receive corrections to there anomalous dimensions these theories. Whether this is also the case in the supersymmetric theory studied in this paper as well requires the computation of anomalous dimensions of these double trace operators which we leave for future investigation.

B Component 3 point functions
In this appendix, we write down the component 3 functions abstractly in term of the functions {A i } appearing in form of full superspace 3 point function (3.17) determined by supersymmetric Ward identity.

C Normalized higher spin three point functions
In this section, we aim to derive the normalized coefficients of the three point functions involving two J 0 operators and one higher spin operator J s (s > 1) leading to eq.(4.16), (4.17), (4.35) and (4.36). The classical action for N = 2 supersymmetric matter coupled to Chern Simons gauge field is given by Note that we will be using the following notation for correlation functions in the superymmetric and the non-supersymmetric theories · · · T k,N Correlator in SUSY theory · · · B k,N Correlator in regular bosonic theory · · · F k,N Correlator in regular fermionic theory · · · f b Correlator in free bosonic theory · · · f f Correlator in free fermionic theory We begin with the derivation of relation between the normalized three point functions of the J b s J b 0 J b 0 in supersymmetric theory and the corresponding three point function in the free bosonic theory. We will employ the path integral technique and utilize some of the relations derived in [20]. The correlation function J b s J b 0 J b 0 T k,N in the path integral representation may be expressed as follows In the first line, we have implemented factorization of planar correlators through the multitrace interactions of the SUSY theory by splitting the action into non-SUSY components [20]. Now, by usual arguments of connectedness and by Wick's theorem, The various integrations and the position labels in the above are implicit. The implicit notation and the momentum space representation are look identical and therefore, the implicit format may be viewed in momentum-space representation. In eq.(20) of [20], an interesting relation was derived between the two point function which is given as Hence eq.(C.4) simplifies upon utilizing the above relation as Note that in order to arrive at the last line in the above relation we have utilized the explicit form of the two point function Thus, from eq.(C.9) we have obtained the required relation A similar calculation allows one to write Upon using the two point functions leads to the normalized coefficient Hence we obtain the required coefficient as 16) In order to derive the other two coefficients in eq.(4.17) and eq.(4.36), we make use of the Giveon-Kutasov duality which relates the magnetic theory to the electric counterpart through the following exchanges [20,34,37] where now the labels (k e , N e ) and (k m , N m ) represent the electric (e) and the magnetic (m) theories, respectively. Now, by the normalization convention of (7), we can rewrite the following, The first equality uses the duality described in eq.(C.17). The second equality makes use of the relation derived in eq.(C.14). The third equality utilizes the results of Maldacena and Zhiboedov derived in [3,4] and the relation between the parameters of the magnetic and electric theory namelyÑ m =Ñ e andλ m = 1 λe . This leads to the following expression for the normalized coefficient This leads to the normalized coefficient C In this appendix, we will describe the relevant diagrams, and corresponding integrals, constructed using the exact 4 point vertex which contribute to the full J (0) four point function. Figure 7 show all the relevant diagrams. For diagram one in 7, Note that the exact vertex (3.10) is a function of two internal grassmann variables. The propagators in fig (7) that emanate from/to the crossed vertices connect these internal Grassmann variables. So one has additional integration of 8 twocomponent internal Grassmann variables, these are not shown in the figure. In the exact three vertex (crossed), for the internal theta arguments, V 3 (θ 1 , θ 1 , θ 1 , q, p) the primed theta is on the direction of the leaving momenta, and the double primed momenta are on the direction of the incoming momenta. The insertion momentum is the first momentum argument, and the incoming momentum is the second momentum argument. In fig (7) the -25 -diagram a is given by There are a total of 6 additional diagrams due to permutations of the operators. and the interaction part is given by The bosonic and fermionic correlators for the diagram fig (7) are given by 4 (q, q , q , θ 1 , θ 2 , θ 3 , θ 4 ) Although we were able to successfully perform the integrals for the components p 3 , θ p and k 3 , θ k in the expression for V (B) 4 given by eq.(D.2) k s and p s integrals out be intractable analytically. Due to this difficulty we were not able to obtain a closed form expression for the four point function of the scalar operatorsJ b 0 and J f 0 in eq.(D.3).

E.2 Decomposition in terms of conformal blocks
The contact diagrams may be written as an expansion in conformal blocks [59]

F.2 Fermionic case
We compute the OPE coefficient relations for scalar double-trace operators with ∆ ≥ 4 using the methodology of the previous subsection.