Marginal deformations of 3d supersymmetric U(N) model and broken higher spin symmetry

We examine the marginal deformations of double-trace type in 3d supersymmetric U(N) model with N complex free bosons and fermions. We compute the anomalous dimensions of higher spin currents to the 1/N order but to all orders in the deformation parameters by mainly applying the conformal perturbation theory. The 3d field theory is supposed to be dual to 4d supersymmetric Vasiliev theory, and the marginal deformations are argued to correspond to modifying boundary conditions for bulk scalars and fermions. Thus the modification should break higher spin gauge symmetry and generate the masses of higher spin fields. We provide supports for the dual interpretation by relating bulk computation in terms of Witten diagrams to boundary one in conformal perturbation theory.


Introduction
Conformal field theories (CFTs) with large N structure are quite important since they can be utilized to investigate holographic dual gravity theory among others. In particular, CFTs with higher spin symmetry broken in 1/N are useful to examine symmetry breaking in dual higher spin gauge theories. In [1,2], the critical O(N ) scalar model and the Gross-Neveu model [3,4] in d dimensions were studied, and the anomalous dimensions of higher spin currents in these models were reproduced by applying the method of conformal perturbation theory. A famous example of higher spin gauge theory on AdS space is given by Vasiliev theory [5][6][7], and the two models are supposed to be holographic dual to the type A and type B Vasiliev theories [8][9][10]. The anomalous dimensions of higher spin currents correspond to the masses of higher spin fields, and the CFT computation in the conformal perturbation theory was interpreted from the bulk theory in terms of Witten diagrams. This paper addresses a continuation of these works.
Higher spin gauge theory is expected to describe superstring theory in the tensionless limit [11], and turning on string tension should correspond to breaking higher spin gauge symmetry. In [1,2], the double-trace type deformations of free O(N ) bosons and free U(N ) fermions were examined, but there are only two fixed points, where one of them is free and the other is interacting with higher spin symmetry broken in 1/N . Since the string tension is a parameter of theory, it should be mapped to a marginal deformation of CFT dual to higher spin gauge theory. Therefore it is better to work on a model which has higher spin dual and admits marginal deformations. As a simple example, we examine N = 2 supersymmetic U(N ) model with N free complex scalars and fermions in three dimensions. 1 The model is proposed to be dual to a supersymmetric Vasiliev theory on AdS 4 [9,10], and it admits deformations exactly marginal at least to the 1/N order.
We compute the anomalous dimensions of higher spin currents in the deformations of supersymmetric model to the 1/N order but to all orders in deformation parameters mainly with the help of conformal perturbation theory. Furthermore, we interpret the CFT computation in terms of bulk higher spin theory. Let us denote φ i and ψ i (i = 1, 2, . . . , N ) as the free complex bosons and fermions in the supersymmetric model. A marginal deformation of the theory is given by The scalar operators are dual to bulk scalar fields, and the deformation corresponds to mixing the boundary conditions of these scalars [16]. Higher spin symmetry is broken with non-zero λ, and the anomalous dimensions of higher spin currents are obtained to the 1/N order but to all orders in λ, see (3.22) below. There is another type of marginal deformation as with K =φ i ψ i andK =ψ i φ i . The fermionic operators are dual to bulk fermionic fields, and the deformation corresponds to mixing the boundary conditions of these spinors. The anomalous dimensions of higher spin currents are computed to the 1/N order but to all orders in κ, see (4.19) below. 2 We can include both the deformations simultaneously, and supersymmetry is preserved when 2λ = κ. 3  We adopt the conformal perturbation theory to compute the anomalous dimensions, which enables us to borrow the previous results in [1,2]. The main reason to use the method is that the CFT computation can be understood in terms of bulk theory. The relation is known between the masses of higher spin fields and the anomalous dimensions of dual higher spin currents, which may be read off from two point function of the currents. The two point function is computable from the bulk theory through Witten diagrams and the contributions to anomalous dimensions arise from loop diagrams as in figure 1. However, the bulk one-loop computation is notoriously difficult in higher spin gauge theory as is well known for gravity theory with spin two gauge field.
In [1,2] (based on previous works in [19][20][21]) the difficulty is evaded by utilizing the Witten's argument in [16] that deforming boundary theory by double-trace type operators is dual to changing the boundary condition of dual bulk fields. Once we admit the dictionary, the diagram with modified bulk-to-bulk propagators along the loop can be rewritten as a product of tree diagrams only with bulk-to-boundary propagators and boundary two point functions as in figure 2. Furthermore, we can see that the bulk computation with tree Witten diagrams is equivalent to the boundary one in the conformal perturbation theory to the 1/N order. Therefore, what we have to do is to show that the modification of bulk-tobulk propagators can be written only in terms of bulk-to-boundary operators and boundary two point functions without relying the AdS/CFT dictionary. This was already done for the bulk scalar fields in [19][20][21]. For the holography with 3d O(N ) scalars, it was applied in [1], and the same idea was already suggested in [20]. In this paper, we first repeat the analysis of [21] in a slightly different way as was done in [19,20], and then apply it to the case with bulk spinor fields.
The organization of this paper is as follows; In the next section, we introduce the supersymmetric model with free bosons and fermions, and write down our notation for higher spin currents. Then we explain how to read off the anomalous dimensions form two point functions of higher spin currents in conformal perturbation theory. In section 3, we compute the anomalous dimensions in the presence of deformation (1.1). We first obtain them to the λ 2 order and then incorporate higher order corrections. In subsection 3.3, we re-examine the model in the formulation with auxiliary fields. In section 4, we study the case with another marginal deformation (1.2). In section 5, we derive the map from the one-loop diagram as in figure 1 to a product of tree Witten diagrams as in figure 2. We first reproduce the known result with bulk scalar propagators and then extend the analysis to the case with bulk spinor propagators. Section 6 is devoted to conclusion and discussions. In appendices, the detailed computations of Feynman integrals are given.

Preliminary
We examine supersymmetric U(N ) model in three dimensions consisting of N complex scalars φ i and Dirac fermions ψ i (i = 1, 2, · · · , N ). The free action is given by where the Gamma matrices satisfy {γ µ , γ ν } = 2g µν . Two point functions are in the coordinate representation and in the momentum representation. Here and in the following we adopt the formula of Fourier transform as (2.4) In the current case, we should set Φ = φ i or ψ i . The free boson and fermion theories have conserved currents J µ 1 ···µs andJ µ 1 ···µs with traceless symmetric indices. Introducing null polarization vector with · = 0, we denote J s (x) = J µ 1 ···µs (x) µ 1 · · · µs andJ s (x) =J µ 1 ···µs (x) µ 1 · · · µs . For the explicit forms, we use the conventions in [22] as (2.6) Here we have represented as∂ With the expressions, the two point functions are computed as (see [2] and references therein) The relationsC will be used for later analysis. We deform the free action (2.1) by adding where ∆ λ S and ∆ κ S are given in (1.1) and (1.2), respectively. We treat the additional part in conformal perturbation theory as Here Φ a is an operator and the right hand side is computed in the free theory.
In the presence of deformations, the higher spin currents J s (x) andJ s (x) are not conserved anymore and possess anomalous dimensions. We read off them from the two point functions of higher spin currents by using the modified correlators as in (2.12). In the current case, there are two types of higher spin currents, and we would meet a diagonalization problem to obtain the anomalous dimensions, see, e.g., [23]. Denoting two independent linear combinations as J ,α s (x) with α = 1, 2, the two point functions are of the forms 13) which are fixed by the symmetry. Here we have usedx 12 = · x 12 . We may change the basis for J ,α s (x) such that with n αβ s , γ α s of order 1/N . Then the two point functions become (2.15) Therefore, the anomalous dimension γ α s of J ,α s to the 1/N order can be read off from the factor in front of log(x 2 12 ).

Deformation with bosonic operators
In this section we deform the supersymmetric model by (1.1) and compute the anomalous dimensions of higher spin currents. We first study them to the first non-trivial order in the deformation parameter λ in the next subsection and then systematically include higher order corrections in subsection 3.2. In subsection 3.3, we rewrite the deformed action by introducing auxiliary fields, and examine the relations to previous works on the critical O(N ) scalars and the Gross-Neveu model such as in [2,24].

Anomalous dimensions to the first non-trivial order
We deform the free theory by the marginal deformation in (1.1). The two point functions of the scalar operators O(x) =φ i φ i andÕ(x) =ψ i ψ i at λ = 0 are given by (3.1) We would like to compute the two point functions of higher spin currents in the presence of marginal deformation. The modified correlators in (2.12) to the first non-trivial orders both in λ and 1/N are Here I 1 , I 2 , and I 3 are defined by the following integrals as In order to read off the anomalous dimensions, we need to extract terms proportional to log(x 2 12 ) from the integrals. Such terms from the above integrals have been already obtained in [1,2], so we can just borrow the results there. See appendix A.1 for details.
As explained in the previous section, we need to solve the diagonalization problem to obtain the anomalous dimensions by choosing a proper linear combination of currents to satisfying the conditions (2.14). In the present case, the proper choice turns out to be The two point functions of these currents are where we have used a relation in (2.10). Using the results in appendix A.1, the anomalous dimensions to the λ 2 order and to the leading order in the 1/N -expansion are read off as 4 Thus we have for odd s and for even s. Here J ,± 1 correspond to spin 1 currents for global symmetry U (1) ⊗ U (1) ∈ U (N ) ⊗ U (N ), which are not broken by the deformation as γ ± 1 = 0. Moreover, J ,+ 2 corresponds to the energy momentum tensor, and hence we have γ + 2 = 0.

Higher order corrections
We would like to move to the higher order corrections in the λ-expansion. A large part of corrections can be incorporated simply by replacing the parameter λ 2 in (3.12) and (3.13) by, say, (3.14) However, not all of the corrections can be treated in this way, and new integrals should be evaluated as well.
A type of correction can be summarized by the replacements of Õ ( In the momentum representation, the two point functions (3.1) at λ = 0 become Therefore, at finite λ, the two point functions become and similarly Therefore, the corrections of this type simply replaceλ 2 byλ 2 /(1 +λ 2 ). For (3.7) we may rewrite as Higher order corrections modify the delta functions (or the identities in the momentum representation) as This means that we should replace eachλ byλ/(1 +λ 2 ), i.e.,λ 2 by (λ/(1 +λ 2 )) 2 . At the λ 4 order, we can see that the following new types of integral appear as The integrals have been already evaluated in [2], and we use the results, see appendix A.1.
Summarizing the above analysis, the two point functions of higher spin currents at finite λ can be given by where we can still use the linear combinations in (3.8). Referring to the results in appendix A.1, the anomalous dimensions to the 1/N order are thus 5 .
This is a main finding in this paper. Setting s = 1, we obtain γ ± 1 = 0, thus the global U(1) ⊗ U(1) ∈ U(N ) ⊗ U(N ) symmetry is not broken to the full orders inλ. For γ + s with even s, we can rewrite as In particular, we have γ + 2 = 0, and this means that the the deformation (1.1) is exactly marginal to the 1/N order.
It is useful to examine the results (3.22) in the limitλ → ∞, which are given as . (3.24) The limit of supersymmetric model should correspond to the sum of the critical U(N ) vector model and the Gross-Neveu model. This can be seen both from CFT viewpoint as in the next subsection and from the dual gravity theory as in [9,10], see also section 5. The anomalous dimensions for the critical U(N ) model can be found in [2] with a slight modification of O(N ) case in [25], and those for the Gross-Neveu model were computed in [17]. See [26][27][28] for recent related works. Comparing these results and (3.24), we find agreements.

Introducing auxiliary fields
In this subsection, we describe the deformation of supersymmetric model by (1.1) as where we have introduced auxiliary fields σ,σ. We may study the model by treating the last term in (3.25) perturbatively. It is useful to move to this formulation for the purpose to see more direct relations to previous works as in [2,24]. This is because the last term in (3.25) vanishes at λ → ∞ limit, where the theory reduces to the sum of the critical U(N ) model and the Gross-Neveu model.
Using the large N factorization of O andÕ, the effective propagators for σ,σ become (see, e.g., [24]) At the limit λ → ∞, corrections to the two point functions of higher spin currents come from the integrals I 1 , I 2 , I 4 , I 5 in (3.5), (3.6), (3.19), (3.20) but with ÕÕ 0 , OO 0 replaced by σσ 0 , σσ 0 . This can be seen from (3.21) atλ → ∞, for instance. The regularization with ∆ used in the appendix A.1 may be understood by the shifts of exponents as where we should set ∆ 1 , ∆ 2 properly. Higher order corrections in 1/λ can be understood also in this formulation. Here we set ∆ 1 + ∆ 2 = 0 for simplicity, but it is not so difficult to work with ∆ 1 + ∆ 2 = 0. Due to the last term in (3.25), the effective propagators should be corrected as Replacing σσ 0 , σσ 0 by σσ 1/λ , σσ 1/λ , we reproduce the previous result in (3.21) except for the contribution arising from the integral of I 3 -type. The integral of I 3 -type contributes from the 1/λ 2 order. This type of correction is absent for the model made from simply the sum of the critical U(N ) scalars and the Gross-Neveu model. At the order, the integral can be written as 6 and similarly for the x 6 -integral. Using these delta functions, we can see that I 3 reduces to I 3 . The higher order corrections can be included by replacing 1 by 1/(1 +λ −2 ) as shown above. Noticing the contribution from this type of integral in (3.21) can be also reproduced.

Deformation with fermionic operators
The supersymmetric model admits another marginal deformation as in (1.2) along with the one in (1.1), which was examined in the previous section. In this section, we compute the anomalous dimensions of higher spin currents in the deformed theory by (1.2). We may use both the deformations simultaneously such as the supersymmetric one with 2λ = κ. In that case we just need to add each contribution to anomalous dimensions. The deformation (1.2) consists of spinor operator K =φ i ψ i (K =ψ i φ i ), whose two point function is at κ = 0. We compute the two point functions of higher spin currents in the presence of marginal deformation. The modified correlators (2.12) to the first non-trivial orders both in κ and 1/N are HereĨ a (a = 1, . . . , 6) are integrals defined as These integrals actually coincide with those evaluated in [2], so we can utilize the results obtained there as summarized in appendix A.2. The anomalous dimensions of higher spin currents can be computed as in the previous section, and the results happen to be the same as (3.10) but λ replaced by κ.
Higher order corrections in κ can be included as in subsection 3.2. One type of corrections can be incorporated by replacing K( In the momentum representation, the two point function (4.1) becomes Thus the two point function receives corrections at finite κ as Here F (p, κ) o linearly depends on the gamma matrix but F (p, κ) e does not as Noticing that trace over the odd number of gamma matrix vanishes, corrections to the inte-gralsĨ 1 ,Ĩ 3 ,Ĩ 5 can be included by replacing F (p) by F (p, κ) o . Therefore we should replacẽ κ 2 byκ 2 /(1 +κ 2 ). Similarly, forĨ 2 ,Ĩ 4 ,Ĩ 6 , we replace 1 in the momentum representation by (4.14) This means that we should replaceκ 2 by (κ/(1 +κ 2 )) 2 .
At the order of κ 4 , there are new types of correction appearing as Further higher order corrections in κ can be included by replacingκ 4 with (κ 2 /(1 +κ 2 )) 2 . These integrals have not appeared in [2], so we have to analyze them in some way. There are several methods to evaluate these integrals, for instance, as in [24]. In appendix B, we do not directly compute these integrals but rather evaluate different ones which are closely related. Using the relation between them, we deduce the log(x 2 12 ) dependence of the integrals as in (B.9), (B.10), and (B.11).
From the above arguments, we conclude that the two point functions of higher spin current are corrected as Using the log(x 2 12 ) dependence of integralsĨ a (a = 1, . . . , 9) obtained in appendix A.2 and appendix B, we can read off the anomalous dimensions as .
This is another main result of this paper along with (3.22). We can see γ − 1 = 0 for anỹ κ but γ + 1 = 0. This means that only a diagonal global symmetry in U(1) ⊗ U(1) is left unbroken. We can also show that γ + 2 = 0 for even finiteκ, which is consistent with the unbroken conformal symmetry.

Bulk interpretations
In the previous sections, we have computed the anomalous dimensions of higher spin currents mainly in the conformal perturbation theory by following the previous works [1,2,21]. The main motivation to use the method is that the CFT computation is closely related to the evaluation of bulk Witten diagram for one-loop corrections to the masses of dual higher spin fields as in figure 1. A crucial point is to rewrite bulk-to-bulk propagators with modified boundary conditions in terms of the product of bulk-to-boundary propagators and boundary two point functions [19,20]. Using this, the one-loop corrections can be computed by using only tree Witten diagrams as in figure 2. In the next subsection, we reproduce the result in [21] by using the intrinsic coordinates instead of embedding ones as was done in [19,20]. In subsection 5.2, we examine the bulk spinor propagators by applying the method.

Bulk scalar propagators
We work on d + 1 dimensional Euclidean AdS space, and use the Poincaré coordinates as Here z 0 is the radial coordinate and z are the coordinates parallel to the boundary, where the boundary is located at z 0 = 0. We consider two bulk scalar fields Φ (1) , Φ (2) with the same mass but with alternative boundary conditions. They are dual to boundary scalar operators O ∆ and O d−∆ with the scaling dimensions ∆ and d − ∆, respectively. In the current case, d = 3 and ∆ = 1 or 2. The bulk-to-boundary propagators are given as [29,30] This leads to the normalization corresponding to the standard kinetic term of dual bulk scalars as Compared with the two point functions in (3.1), the relations are We denote G ∆ (z, w) as the bulk-to-bulk propagator, and we will need the following expression as (see, e.g., [19,31]) Here k = | k| and K ν (z) and I ν (z) are the modified Bessel functions. We deform the dual boundary CFT with where f can be identified withλ in (3.14). As shown in [16], this deformation is dual to mixing boundary conditions for Φ (1) , Φ (2) . Due to the changes of boundary condition, the scalar propagators (2) are exchanged, and this is consistent with the fact that the dual CFT reduces to the sum of the critical U(N ) scalars and the Gross-Neveu model.
Here we would like to examine the bulk Witten diagram in figure 1, where the bulk-tobulk scalar propagators with modified boundary conditions in (5.7) are used for the loop. It is a quite difficult task to directly evaluate loop diagrams with higher spin gauge fields generically. However, the AdS/CFT duality suggests that the modified scalar propagators should be written in terms of scalar propagators with f = 0 but with the insertions of boundary operators. More precisely speaking, we should have the following identities as The same should be true for G 22 f and G 21 f by exchanging ∆ with d−∆. With these identities, the loop Witten diagram in figure 1 can be evaluated only from the product of tree Witten diagrams as in figure 2. 7 Moreover, we can see that the latter description is dual to the boundary computation in conformal perturbation theory. Therefore, we can obtain the map between the bulk computation with Witten diagrams and boundary one in conformal perturbation theory, once we can confirm the identities (5.8) and (5.9). We can check that the identities (5.8) and (5.9) follow the two basic ingredients In the rest of this subsection, we shall derive them. They were already shown as (B.21) and (5.7) in [21] by adopting the embedding formulation. For the extension to bulk spinor propagators in next subsection, we re-derive the results using the intrinsic coordinates in (5.1). The derivation of (5.11) here is the same as the one in [19], see also [20]. For our purpose, it is convenient to work in the momentum representation as Using these expressions, we can show as in (5.10). Similarly, we have for (5.11) Applying the formula (ν = ∆ − d/2) we find With the expression of bulk-to-bulk propagators in (5.5), we arrive at (5.11).

Bulk spinor propagators
There is another type of marginal deformation as in (1.2), and they should be dual to modifying boundary conditions for bulk spinor fields. Thus the anomalous dimensions of higher spin currents can be computed from one-loop Witten diagram as in figure 1, but now bulk spinor fields with modified boundary condition are running along the loop. In this subsection, we show that the one-loop Witten diagram can be evaluated from the product of tree Witten diagrams as in figure 2 just like for the bulk scalar propagators. We examine d+1 dimensional Euclidean space with the coordinate system (5.1). Moreover, we consider two bulk Dirac fermions Ψ (1) , Ψ (2) with the same mass m but alternative boundary conditions. The spinor operators dual to Ψ (1) , Ψ (2) are given as K ∆ + , K ∆ − , whose conformal dimensions are ∆ ± = d/2 ± m. For application to the current case, we should set d = 3, m = 0, and consider only half of the fermions.
Spinor propagators in AdS d+1 have been investigated in [31,34,35], and we follow their analysis. We introduce the Gamma matrices Γ a (a = 0, 1, · · · , d) satisfying {Γ a , Γ b } = 2δ ab . Defining the bulk-to-boundary propagators are given as The boundary two point functions are Compared with (4.1), the normalization of spinor operators differs from K by the factor 16/N . As in [31], the bulk-to-bulk propagators should satisfy and furthermore regularity and boundary condition. Solutions were obtained as Here we have defined Since we have the solution with −m can be obtained from S m ± (z, w) by replacing Γ 0 by −Γ 0 and (z, w) by (−z, −w). This implies that which will be used for later analysis.
We deform the dual boundary CFT by Here we have two terms since now we are working with two Dirac fermions. The deformation parameter h can be identified withκ in (4.13). The deformation is dual to mixing the boundary conditions for Ψ (1) , Ψ (2) , and the modified spinor propagators Ψ (a) (z)Ψ (b) (w) = S ab f (z, w) are [33] S ab h (z, w) = .
We can see that the propagators for Ψ (1) , Ψ (2) are exchanged at the limit h → ∞. Below we shall show that Using these identities, we can rewrite the modified spinor propagators (5.30) as those with h = 0 but with the insertions of boundary operators precisely as for the bosonic case, see (5.8) and (5.9). In order to proof the identities (5.31), (5.32) and (5.33), it is convenient to express the propagators (5.18) and (5.20) in the momentum representation as ±,m (z 0 , k)P ∓ , and Then we can show for (5.31) and similarly for (5.32). For (5.33) we rewrite ±,m (z 0 , k)P ∓φ ±,m (z 0 , k) = ∓ we find (5.33) with the help of (5.28).

Conclusion and discussions
In this paper, we have studied 3d supersymmetric U(N ) model, which is supposed to be dual to 4d supersymmetric Vasiliev theory [9,10]. The model admits two types of marginal deformation as in (1.1) and (1.2), and the deformations should be dual to modifying boundary conditions for dual scalars and spinors, respectively [16]. There are two main results of this paper. One of them is on the anomalous dimensions of higher spin currents in the deformed models. Using basically conformal perturbation theory, we obtained them to the leading order in 1/N but to all orders in deformation parameters as in (3.22) and (4.19). For higher order corrections in κ in the deformed theory by (1.2), we utilized another technique as used in [17]. The other is on the dual higher spin interpretation of the computations in conformal perturbation theory, which can be done by rewriting one-loop Witten diagram as in figure 1 to the product of tree Witten diagrams as in figure 2. We have derived the ingredients essential for the rewritening in (5.10), (5.11) for bulk scalars and in (5.31), (5.32), (5.33) for bulk spinors. This was already done in [21] for bulk scalars (see [19,20] for (5.11)) but it was new for bulk spinors. We examined the 3d supersymmetric model as a simple model which has marginal deformations and higher spin holographic dual. We would like to relate deforming the CFT marginally and turning on the string tension. In order to reveal more precise relations, however, we need to extend the analysis to more complicated systems. For instance, a concrete relation (named as ABJ triality) was proposed in [18] through 3d ABJ(M) theory in [36,37]. For application to the ABJ triality, we have to work with coupling to Chern-Simons gauge fields as in [12][13][14]. See [15] for a recent work. There are also lower dimensional proposals with 2d CFTs which have N = 4 supersymmetry in [38,39] and N = 3 supersymmetry in [33,40,41]. Symmetry breaking for higher spin gauge theory on AdS 3 has been studied in [21,23,33,42], but more detailed analysis would be required to say more concretely.
There are other open problems as follows; We have evaluated the log(x 2 12 ) dependence of integralsĨ 7 ,Ĩ 8 ,Ĩ 9 in appendix B with an indirect method. However, it is desired to evaluate them explicitly, see also footnote 2. A drawback of our method is that we cannot identify the Nambu-Goldstone modes arising due to symmetry breaking as mentioned in [1,2]. A direct way to evaluate the loop corrections to the bulk higher spin propagators explicitly, but it looks a quite hard task. Along with the difficulty for evaluating loops, there is an additional problem on our limited understanding of the dual Vasiliev theory, where only classical equations of motion are available presently. An indirect way might be given by group theoretic analysis as was done for the 4d Vasiliev theory dual to 3d critical O(N ) model in [43]. It would be also possible by writing the divergence of dual higher spin current in terms of products of higher spin currents, see, e.g., [14,15,26,27]. The 3d supersymmetric U(N ) model itself is also worth studying furthermore. The deformations in (1.1) and (1.2) are exactly marginal to the 1/N order, but it is not clear what happens at the next order. Even if we could find non-trivial fixed lines (or points), it is notorious to be a hard problem to go beyond the leading order in 1/N . Hopefully, this issue would be tractable along the line in, e.g., [28,[44][45][46]. It is also important to extend our analysis to, say, the theory with both bosons and fermions in generic d = 3 as in [47,48]. In the case, the deformations (1.1) and (1.2) are not marginal anymore, so we may examine the fixed points of RG flow induced by these deformations.

Acknowledgements
The work of YH is supported by JSPS KAKENHI Grant Number 24740170, 16H02182.

A Evaluations of integral
In this appendix, we examine the integrals I a (a = 1, . . . , 5) andĨ a (a = 1, . . . , 6) defined in the main context. We denote the log(x 2 12 ) dependence of these integrals as I log a andĨ log a , and we shall obtain them explicitly by borrowing the results in [1,2].

A.1 Deformation with bosonic operators
In the case with (1.1), the integrals are almost the identical to those evaluated in [1,2], so we can directly use the results there. The integral I 1 in (3.5) is almost the same as the one in section 3.2 of [2] with d = 3, once we replace the correlator by the propagator of auxiliary field σ. Therefore, we may shift the exponent as (x 2 12 ) −2 → (x 2 12 ) −2+∆ . Then the reside at ∆ = 0 should be the factor in front of log(x 2 12 ). Since we are working with complex scalars, we need to multiply some factors as discussed in section 5.2 of [2]. Moreover, we have to replace the overall factor C σ for the auxiliary field propagator by 2N C 2 φ . We then find We can evaluate I 2 in (3.6) as for I 1 . Replacing by the propagator of auxiliary fieldσ, the computation of I 2 is almost the same as that of I 1 in subsection 4.2 of [2]. Thus we shift as (x 2 12 ) −1 → (x 2 12 ) −1+∆ and pick up the 1/∆-pole term. Replacing the overall factor Cσ for the auxiliary field propagator by N C 2 φ , we find The integral I 3 in (3.7) vanishes for s odd, so we set s even here. Using the three point functions obtained in [2], the integral I 3 can be written as The same integral arises from the O(N ) scalars as I we find For the integral I 4 in (3.19), we adopt the regularization used for I 1 . Then the integral is the same as the one computed in section 3.3 of [2]. Since we are working with U(N ) global symmetry, we need to multiply each three point function by 1/4. Moreover, we should replace C σ by 2N C φ . Thus we find For the other integral I 5 in (3.20), we adopt the regularization used for I 2 . Then the integral is the same as the one computed in subsection 4.3. Replacing Cσ by N C φ , we obtain Here we should remark that the contributions to the log(x 2 12 ) terms are not the 1/∆-pole terms but twice of them. See [24] for more details, for instance.

A.2 Deformation with fermionic operators
We move to examine the integrals appeared in section 4, where we have deformed the theory by (1.2). The integralĨ 1 in (4.5) is computed as The integral is essentially the same as I 1 computed in subsection 3.2.1 of [2], so we evaluate the integral in the same way. Namely, we replace (x 2 34 ) −5/2 by (x 2 34 ) −5/2+∆ and pick up the 1/∆-pole term. The log(x 2 12 ) dependence can be read off as The integralĨ 2 in (4.6) can be written as where B k,l was defined in (A.6). As for I 3 , we find The integralĨ 3 in (4.7) becomes The integral is almost the same as I Similarly, we obtaiñ The integral is almost the same as I The integralsĨ 5 in (4.9) andĨ 6 in (4.10) are of the similar form as Using the relation (see (4.41) and (4.42) of [2]) we can rewrite asĨ a k a l B k,l (A. 19) with B k,l in (A.6). The integralĨ 6 is twice ofĨ 2 evaluated as in (A.13), which leads tõ

B Another method for anomalous dimensions
In the previous appendix, we have obtained the log(x 2 12 ) dependence of integrals I a (a = 1, . . . , 5) andĨ a (a = 1, . . . , 6). Therefore, only the integralsĨ 7 ,Ĩ 8 ,Ĩ 9 in (4.15), (4.16), (4.17) are left. In principal, it is straightforward to evaluate the integrals as in [24], but it is quite tedious. In this section, we compute the log(x 2 12 ) dependence by developing more simple but indirect method.
For explanation, let us focus on the case of the Gross-Neveu model. In case with our main method, we compute J s (p)J s (−p) and read off the anomalous dimensions as explained in section 2. The same information can be actually read off from the three point function J s (0)ψ i (p)ψ j (−p) by solving the RG flow equation as in [17], see also [15]. It is obvious to identify which Feynman diagrams for the three point function correspond tõ I 7 ,Ĩ 8 ,Ĩ 9 , thus we can successfully deduce the terms proportional to log(x 2 12 ) in the integrals. However, in this way, we loose direct connection to the bulk computation as explained in section 5, which was actually the main purpose in previous works [1,2], see also footnote 2. Moreover, if we are interested in the overall normalization of two point function as in [24], then we have to directly evaluate the complicated integrals.
We need to extend the method in [17] to the case of the supersymmetric model deformed by (1.2). As in subsection 3.3, we rewrite the action by introducing auxiliary spinor fields η,η as The effective propagator for η,η is obtained as As in [17] we shift the exponent as which is the same as changing the scaling dimension of η,η from 3/2 to 3/2 − ∆, see, e.g., [24]. This also implies that interaction terms become at κ → ∞ with the renormalization scale µ. We read off the anomalous dimensions from the µ dependence of three point functions, such as, J s (0)ψ i (p)ψ j (−p) . The vertex operators corresponding to J s (0) andJ s (0) are We may denote the logarithmic corrections to the vertices V 0 α,s (0, p) as Here γ αβ s were given in (2.13), and they become anomalous dimensions after diagonalization. 8 There are three types of contributions as in figure 3. Two of them can be computed by left two Feynman diagrams, which correspond to the integralsĨ a (a = 1, . . . , 6). Another type of corrections come from the right most Feynman diagram. They correspond to the integralsĨ 7 ,Ĩ 8 ,Ĩ 9 , so we shall focus on the case from now on.
In the rest of this appendix, we compute the log µ corrections δV αβ s (0, p) corresponding toĨ 7 ,Ĩ 8 ,Ĩ 9 . Let us expand the value of Feynman integral (denoted as X(∆)) by the shift ∆ introduced in (B.3) as We will compute the Feynman integrals of the form in the right most diagram of figure 3, which has µ 4∆ factor from the interaction terms (B.4). Therefore, we can read off the log µ dependence from the 1/∆-pole term X (−1) as At κ → ∞, the contributions to anomalous dimensions can be identified with some factors in front of log(x 2 12 ), see section 2. Carefully treating the (16/N ) 4 factor as in (4.18) and the normalization of higher spin currents, we would obtain In order to evaluate integrals, we utilize several formulas; We introduce Feynman parameters as .
Moreover, we use the following integrals Let us start from δV 11 s corresponding to the integralĨ 7 in (4.15). The integral we need to evaluate is where F κ (k) is given in (B.3) and Let us work on A 11 s (k) first. Introducing the Feynman parameter as in (B.12), integration over l can be performed as .

(B.24)
The second term is The singular term arises as (1 − x)( / l + x/ k)(l + xk) s ((l ) 2 Collecting the two contributions, we find There are two terms in A 22 s (k), and we start from the contribution to δV 22 s from the first term. We compute For the contribution from the second term, we find The second term is twice of (B.34) and the first term becomes As mentioned in footnote 8, we need to change the relative normalization of vertices V 0 1,s (0, p) and V 0 2,s (0, p) in (B.5) so as to obtain γ 12 = γ 21 .