Supersymmetric Tensor Model at Large $N$ and Small $\epsilon$

We study the $O(N)^3$ supersymmetric quantum field theory of a scalar superfield $\Phi_{abc}$ with a tetrahedral interaction. In the large $N$ limit the theory is dominated by the melonic diagrams. We solve the corresponding Dyson-Schwinger equations in continuous dimensions below $3$. For sufficiently large $N$ we find an IR stable fixed point and computed the $3-\epsilon$ expansion up to the second order of perturbation theory, which is in agreement with the solution of DS equations. We also describe the $1+\epsilon$ expansion of the model and discuss the possiblity of adding the Chern-Simons action to gauge the supersymmetric model.


Introduction and summary
In recent literature, there has been strong interest in theories whose dynamical fields are tensors of rank 3 or higher (for reviews, see [1][2][3]). Such theories possess a number of interesting features. For example, only the melonic diagrams dominate in the large N limit, in contrast to the vector models, where only snail diagrams dominate [3], and the matrix models, where all the planar diagrams survive in the large N limit. This fact makes the tensor models similar to the famous Sachdev-Ye-Kitaev (SYK) model [4][5][6]. The SYK model contains a disordered coupling constant, making it hard to use standard tools of quantum field theory. The SYK model is believed to describe quantum properties of the extremal charged black holes [7][8][9] and therefore may help to serve as a toy model for understanding the AdS/CFT correspondence [10][11][12]. It is already used for understanding the properties of the traversable wormholes [13][14][15][16]. While the tensor models [1] exhibit the same properties at the large N limit, they do not have disorder therefore giving us hope that they can be understood at finite N via standard techniques of quantum field theories. These techniques have already brought many interesting results [17][18][19][20][21][22][23][24][25][26][27].
We shall consider a supersymmetric analogue of such theories, which has been recently considered as a generalization of SYK model [28][29][30] or as a quantum mechanical supersymmetric tensor model [31][32][33][34]. Here we will present a similar model in continuous dimension d. We consider a minimal N = 1 supersymmetric model, where we have some number of scalar superfields Φ abc (x, θ), and indices a, b, c run from 1 to N . These fields are coupled via a "tetrahedral" superpotential 1 S = d d x d 2 θ 1 2 (D α Φ abc ) 2 + gΦ abc Φ ab c Φ a bc Φ ab c . . This model has been proposed in the paper [21] as a generalization of the scalar melonic theory. It was proved that the non-supersymmetric analogue of this theory has a so-called melonic dominance in the limit when N → ∞, g → 0 but gN 3 2 is kept fixed [36]. The proof of this peculiar fact relies on the combinatorial properties of the potential, and therefore is applicable in any dimensions and in various theories, provided that the combinatorial properties are left the same. In the case of the supersymmetric theories, the Feynman diagrams, written down in terms of the components, look quite complicated and, at first glance, do not possesses a melonic limit as in the case of scalar model or the SYK model. However, one can develop a supersymmetric version of the usual Feynman diagrams technique and work explicitly with the superfields Φ abc and see that the combinatorial and topological properties are the same as in the case of the scalar tensor models. Therefore, the proof of the dominance of melonic diagrams [20,21,[36][37][38] is applicable in this case and the theory (1.1) also possesses a melonic dominance in the large N limit. We generalize the theory (1.1) where the tethrahydral term is replaced by qvalent maximally single-trace operator to study models with different numbers of the internal propagators in each melon [37,39].
The properties of such theories in the IR limit can be investigated by solving the Dyson-Schwinger (DS) equations, which are drastically simplified if the theory is melonic. Namely, the dominance of the melonic diagrams in the large N limit can be understood as a suppression of the corrections to the vertex operators in the system of DS equations. The solution of the DS equation in the IR yields a conformal propagator, suggesting that the theory in the IR flows to the fixed point, which is described by some conformal field theory. The existence of the stress-energy tensor with the correct dimension, and the spectra of the operators confirm this hypothesis. Therefore, one can wonder whether it is possible to describe such a transition from the UV scale (where we have a bare conformal propagator determined by commutation relations) to the IR region by means of RG flow and expansion. Several attempts have been made towards this idea. For example, the melonic scalar theory in 4 1 Here we will refer to the appendix A and the paper [35] for the notations and the other helpful formulas that will be used through the paper.
dimensions [40] has been considered at the second order of the perturbation theory. For this theory, a melonic fixed point of RG flow was found, even though the corresponding couplings are complex. The complex couplings indicate that the theory is unstable. For example, the dimensions of some operators have imaginary part. One of the reasons of instability could be that the potential is unbounded from below, leading to the decay of the vacuum state.
It is quite interesting that if one drops the fermionic part of the action (1.1) and integrates out the auxiliary field, the theory still possesses the melonic dominance in the large N limit.
Such a "prismatic" theory was considered in the paper [41]. The solution of this theory was found in the large N limit and the RG properties were investigated at two loops. As opposed to the standard melonic theory [40], the fixed point is real and first order of expansion recovers the exact solution in the large N limit.
In this paper we solve the model (1.1) in the large N limit, assuming that the supersymmetry is not broken and that in the IR region the theory is described by the conformal propagator. The solution is found for general dimension d and general q-valent MST potential [37,39]. The dimension of the operators at given d and spin s can be found as a solution of the corresponding transcendental equation. It is shown that at any dimension d, there is always a stress-energy operator of dimension d and a supercurrent operator of dimension d − 1 2 , which indicates that the theory is indeed described by a conformal field theory. While the model (1.1) exists only in the fractional dimensions between one and three dimensions, the counterpart SYK model with q = 3 can work at the integer dimension d = 3 and describe a good conformal field theory with the melonic dominance in the large N limit. After that we derive a perturbation theory in 3 − dimensions of the theory (1.1) to find a fixed point that could describe the IR solution of the large N limit of the model (1.1). We find that the expansion is consistent with the exact large N solution up to the first order in . The two-loop analysis also suggests that the found melonic fixed point is IR stable.
The structure of the paper is as follows: in section two, we discuss the properties of the theory (1.1) in the large N limit. The dimensions of the operators are found and the DS equation is solved in the superspace formalism. In section three, we consider q = 3 supersymmetric SYK model and study the stability of such a theory. In section four, we study the RG properties of the quartic super theories in 3 dimensions and compare to the exact solutions in the large N limit. In section five, we discuss the possibility of introducing higher order supersymmetry and speculate about the consequences of gauging the supersymmetric tensor models. The appendix provides supplemental materials including the notations and = + + . . . useful formulas that are used throughout the paper.

Solution of the Large N Theory
In this section, we will try to find the solution of DS equations for the theory (1.1) in the large N limit. As mentioned in the introduction, the theory possesses a melonic dominance in the large N limit. This means that only specific diagrams survive in the large N limit, namely the ones generated recursively by the Dyson-Schwinger (DS) equation (schematically depicted in the fig.(1)). The resulting equation for scalar or fermion field theories was investigated analytically and numerically for many different theories [6,21,42].
For example, the DS equation can be solved in the IR limit and the solution possesses a conformal symmetry in that limit. In the case of the supersymmetric theories, one of the important differences is that one can demand the solution to respect supersymmetry. In order to do it manifestly the DS equation should be formulated in terms of the superfields.
Of course, one can do this calculation in terms of the components as in the paper [31] and check that these two approaches give the same answers. To make the discussion more general we consider the case where there are q − 1 internal propagators in the melon diagrams and suitable MST operator is considered [37]. The DS equation in the supersymmetric case reads as G(p; θ, θ ) = G 0 (p; θ, θ )+ (2.1) where G 0 (p; θ, θ ) is a bare superpropagator (A.10), G(p; θ, θ ) is an exact superpropagator and g = λN where ∆ < ∆ 0 = 1 for the solution to be valid in the IR limit [7] (namely, we can neglect by bare propagator in comparison to the exact one G −1 0 G −1 , p → 0 ). Substituting the ansatz in the DS equation (2.1) we get As soon as ∆ < 1 we can neglect the LHS of the equation by the RHS in the limit p → 0.
After that one can integrate out Grassman variables using identities for the superderivative to get This equation gives the dimension of the superfield to be ∆ = d(q−2)+2 The solution suggests that we cannot work directly in d crit (q) = 2q−2 q−2 dimensions because the bare propagator is not suppressed in the IR limit and change the solution. For example, for the case of tetrahydral potential q = 4, d crit = 3, therefore the tensorial melonic theory is not conformal in 3 dimensions. Nevertheless, we can still study the theory slightly below 3 dimensions and compare it with the expansion.
If one chooses the case of q = 3, the critical dimension is d crit = 4 and such a melonic theory should describe a conformal field theory in 3 dimensions. In the next section we will review this model in more details.
We calculated the propagator (2.2) in the momentum representation. One can carry out the calculation in the coordinate space. With the use of the relation the propagator in the coordinate representation is Another way to see that the dimension of the superfield is d−1 q is to rewrite the action in terms of the components and impose the conditions ∆ ψ = ∆ φ + 1 2 , then the action contains a term The solution (2.2) suggests that in the IR limit, the theory is described by some conformal field theory (CFT). One of the interesting questions that one may ask is, what is the spectrum of the bipartite conformal operators in this theory? The supersymmetric theory (1.1) has different types of the bipartite operators, as the prismatic one [41]. We should consider these families separetly. The most simple ones have the following structure [29] A similar set of the operators was considered in the paper [29] in 2 dimensions and [31]   for arbitrary d. We can try to put more D 2 in (2.9) to get more familes, but with the use of the identity (D 2 ) 2 = , one can descend these operators to the BB or FF series. That's why we can consider only these two families to get the whole spectrum of bipartite operators with the lowest component having spin s = 0.
As usual, the corrections to the bilinear operators in the large N limit are given by the ladder diagrams (but again, in comparison to [7,21], these diagrams should be considered to be in superspace). We assume the following ansatz in momentum space for the three-point correlation function for these families, where we have set the operators V BB , V F F to be at infinity and made a Fourier transformation with respect to the spatial coordinates, and ∆ V is the corresponding dimensions of the operator. The derivation of the equations for the dimensions ∆ V is just a straightforward generalization of the analogous calculation for the scalar model [21] or the SYK model [6].
Here we will show the derivation of such equation for the BB operators.
The addition of one step of the ladder can be considered as the action of the kernel (2.12) We act on the (2.9) by one step of the ladder, The Grassman variables can be integrated out with the use of identities from the section A.
After that we are left with a simple integral In order for the operator to be primary, the equation g B (∆ V ) = 1 must hold. An analogous equation can be written for the V F F operator, but one can see that This suggests that we might build a bigger multiplet and enhance the supersymmetry to be N = 2 (later we shall see that this does not actually happen, because there is no additional fermionic counterparts to finish supermultiplet).
From now on we shall consider the case only q = 4 to get 3 − expansion unless the other is specified. Thus, we can get the expansion in the large N limit of the Φ 2 operator The plot of the ∆ Φ 2 as a function of the dimension is depicted in the figure 3. Analogously we get the dimension of ΦD 2 Φ operator We can discuss dimensions of non-singlet operators of the form Φ abc Φ a bc . The equation for the dimension of this operator can be rewritten as where a factor q − 1 appears from the combinatorics [43], and ∆ aa is the dimension of the operator. The expansion near three dimensions for q = 4 has the following form Later, we shall show that the solution coincides with the expansion in the second level of perturbation theory.
From this, the next step would be to study the spectrum of the higher-spin operators. A generalization for the higher spin operators is with the corresponding modifications for the ansatz. For example, for higher spin spectrum of the BB operators the ansatz is (2.22) In this case we consider an arbitrarily chosen null-vector ξ µ and consider the convolution of the ansatz (2.22) with the vector ξ. After that one can integrate out the Grassman variables and carry out the integration over the real pace with the use of a relation [40]: Eventually, the equation for the dimension at given spin s reads as One would expect that there is a solution at any d and s = 2 with ∆ = d, because of the existence of the stress-energy tensor. However one cannot find this solution. The reason is quite simple. First of all, there is no stress-energy tensor in the field decomposition of the BB and FF operators. Second, the stress-energy tensor has a superpartner S α µ (corresponding to supertranslations) that has spin 3 2 , and therefore to find it we should consider a whole different family of operators, with lowest component being a Rarita-Schwinger field . Namely, let us consider a Fermi conformal primary operator where the odd number of the space-time derivatives should be inserted to get a primary operator. Indeed, if we consider a zero number of the derivatives The derivation of the equation for the spectrum of the dimensions is straightforward where the spin should be chosen to be of the form s = 2n − 1 2 . Now we can try to find the stress-energy momentum and its partner. And indeed at any d, q and s = 3 2 there is an operator with dimension ∆ = d − 1 2 that corresponds to the usual stress-energy supermultiplet. At this point one can wonder whether the current J aa , responsible for the O(N )'s transformations, is a primary operator. The supersymmetric multiplet containing the current should be also a Fermi supermultiplet with spin s = 1/2 (this operator is not a singlet operator and therefore (2.25) is not applicable). The current should satisfy the equation [43] g aa BF (d, ∆ V , s) = The relation (2.24) can be thought as a generalization of the equation for the kernel at 2 dimensions derived by Murugan et al. [29]. In this case they introduced two dimensions, h = ∆+s 2 andh = ∆−s 2 , and one can check that which coincides with the equation (7.17) in [29].
The relation (2.28) also shows that if there is a scalar bilinear multiplet with dimension h, there is no BF operator with higher spin and the dimension ∆ = ∆ + 1 2 . This shows that we cannot complete the N = 2 supermultiplet and the enhancement does not happen.
It is interesting that there is an argument in d = 1 stating that it actually must happen.
Basically, it comes from the fact that group of diffeormorphims of supertransformations in 1-dimension comprises the N = 2 superalgebra [29].
Finally we discuss the dimension of the quartic operators, because there is a fundamental relation between their dimensions and the eigenvalues of the matrix ∂β i ∂g j . We can find the dimensions of some quartic operators in the large N limit. For example, in the matrix models the anomalous dimension of a double trace operator is just the sum of the anomalous dimensions of the corresponding single trace operators. By the same analysis, we get that the anomalous dimension of the double trace operator is (2.31) Analogous analysis gives that Finally, the dimension of the tetrahedral operator can be determined as the dimension of the operator Φ abc D 2 Φ abc (namely, it follows from the equations of motion) and it gives us

The expansion near one dimension
One can try to study the behaviour of the model (1.1) near 1 dimension. The case of d = 1 supersymmetric tensor models was considered recently (see [31]). It was found that the supersymmetry is broken in the IR region. The easiest way to see this is to assume a conformal ansatz and plug it in the DS equation (  3 Supersymmetric SYK model with q = 3 in d = 3 In the previous section we mostly work with the tensor models in non-integer dimensions.
The main problem that did not allow us to work directly in 3 dimensions was that the critical dimension for such a interaction is d cr = 2q−2 q−2 = 3, meaning that directly at 3 dimensions the conformal IR solution does not work. Nevertheless, if one considers q = 3 case the critical dimension becomes d cr = 4 and therefore should work perfectly in 3 dimensions.
Unfortunately, we do not know any q = 3 tensor model and in order to somehow study this melonic model we shall consider a SYK like model with disorder, which is a special case of the models [29].
Thus, we shall try to study the following model where we consider a quenched disorder for the coupling C ijk . One might worry, that such a theory violates the causality, because the field C ijk is assumed to have the same value across the space-time and therefore the excitation of such a field changes the value of it everywhere, thus violating causality. But the procedure of quenching requires firstly to fix the value of C ijk that makes the theory casual and after that average over this field. It means that we can not excite the field C ijk and violate casuality.
This model is similar to the tensor one considered in the previous section, because again only melonic diagrams survive in the large N limit, but with two iternal propagators in each melon. Therefore the formulas derived in the previous section are applicable in this case and with the replacement of λ → J and setting q = 3 we can recover the large N solution of this model. For example, the propagator in this case is and the dimension of the field Φ i is ∆ = 2 3 . Again the spectrum of the operators could be separated into three sectors, described in the previous section. The equation for the BB operators is determined by the equation where s is the spin and should be chosen even. One can try to find the spectra of low lying states (5)  It is easy to see that the spectrum has the following asymptotic behavior at large spins On a principal line h = d 2 + iα the kernel is complex, it is connected to the fact that there is no well-defined metric in the space of two-point functions [29]. Therefore there is no problems with the complex modes, that could possibly destroy the conformal solution in the IR [16]. Thus q = 3 supersymmetric SYK model is stable at least in the BB channel. Also one can check there are no additional solutions to the equation g 3,3 BB (h, s) = 1 in the complex plane except the ones on the real line. The spectrum of the FF operators coincides with the spectrum of the BB operators but shifted with h → h + 1, therefore we don't have to worry about the instabilities of the theory in this sector.
Analogous calculations could be conducted for the BF series where the spin s should be in the form s = 2n − 1 2 . One can notice that there is a solution g 3,3 BF (5/2, 3/2) = 1 corresponding to the existence of the supercurrent and energy momentum tensor (the energy momentum is not seen directly because it belongs to the supermultiplet of the supercurrent, but if one studies the theory in terms of the components, he or she will of course find the energy momentum tensor). There is a list of some low lying operators in the FF sector (5) [∂ µ ΦD α Φ] θ=0 s =    [40,41,44]. We include all possible O(N ) 3 symmetric marginal interactions that respect the supersymmetry. Thus the superpotential has the following form where we imposed a symmetry under the exchange of the colors. In comparison to the "prismatic" theory [41], which has 8 coupling constants, the supersymmetric theory has only 3; this is a significant simplification.
Let us first consider the general renormalizable d = 3 theory of N = 1 superfields Φ i , i = 1, . . . n: where Y ijkl is a real symmetric tensor. Adapting the results from [45,46], we find that the two-loop corrections to the gamma and beta functions are These two-loop results are closely related to those in a non-supersymmetric theory with Yukawa coupling 1 4 Y ijkl ψ i ψ j φ k φ l (see [46]), except the supersymmetry requires Y ijkl to be fully symmetric.
If we choose N = 1, the couplings g 1 , g 2 , g 3 becomes degenerate because they describe the same operator. Therefore, the beta-functions should be added to get the right expression.
And indeed, if we choose N = 1 and sum up the couplings we get which is the correct beta function for the theory with superpotential (g 1 + g 2 + g 3 )Φ 4 for a single chiral superfield Φ. This special case of our theory is conformal in the entire range 2 ≤ d < 3. Indeed, in d = 2 the N = 1 supersymmetric theory with superpotential Φ m for one superfield Φ flows to the superconformal minimal model with central charge . Let us consider the large N limit where we scale the coupling constants in the following way: The scaling is taken to be the same as in the paper [40]. Applying this scaling to the formula (4), we get 10) From this one can find the fixed point in the large N limit. Namely, We may try to compute the 1/N corrections to these results to get We can study the dimension of various operators in the fixed point (4.11). One of these operators is Φ 2 abc , which belongs to the BB spectrum. We can find that the anomalous dimension of this operator is (4.13) where we have used the relation γ Φ 2 = 6γ Φ , which is true only at the second level of perturbation theory. The answer coincides with the exact solution found earlier (2.17).
As one can see, the fixed point (4.11) is IR stable, which means that the dimensions of the operators is bigger than the dimension of the space-time. Indeed, the linearized equations of RG flow near the fixed point (4.11) have the following eigenvalues but as it is known the eigenvalues of this matrix gives the dimensions of quartic operators

N = 2 supersymmetry and gauging
One can try to consider N = 2 supersymmetry and study the properties of such a model.
Here we are not going to present the solution of the corresponding DS equation , but we will just calculate the beta-functions and find the fixed point of the resulting equations. The SYK model with N = 2 supersymmetry at 2 dimensions was considered in the paper [30].
The theory is built analogously to the N = 1 case. It can be obtained by dimensional reduction from N = 1 supersymmetry in 4 dimensions. In this case, we have a set of chiral superfields Ψ abc with the action where the superpotential is taken to be the same as in the case of N = 1 supersymmetry.
The beta-function for a general quartic superpotential was considered in the paper [47]. The beta-function receives corrections only from the field renormalizations, meaning that it has the following form The fixed point is determined by demanding that the anomalous dimension of the field must , as we got for a general melonic theory in arbitrary dimensions. Apparently, for N = 2 models this fact comes not from the melonic dominance, but from the consideration of the supersymmetric algebra that fixes the dimensions to be proportional to the R charge of the corresponding operator. This condition defines a whole manifold in the space of marginal couplings. Applying the scaling (4.9), in the large N limit we get the equation It is quite interesting that this equation does not fix λ 2 , λ 3 in the large N limit. One can study the stability of these fixed points at arbitrary λ 2,3 . The RG flow near the fixed point could be linearized to get the stability matrix The given solution is marginally stable, because of the existence of two marginal operators. These two zero directions correspond to the previously discussed existence of a whole manifold of IR fixed points.
From this consideration, it would be interesting to study the large N limit of the considered N = 2 theory and corresponding DS equations. This model must have the same combinatorial properties as the N = 1 and scalar tensor model, but some cancellation happens that drastically simplifies the theory.
One can try to examine a gauged version of N = 2 theory. The gauging of the tensor models is one of the important aspects that makes them different from the SYK model. In the latter, due to the presence of the disorder in the system, the theory can possess only the global O(N ) symmetry and can not be gauged, while in the tensor models there are no such obstructions and one can add gauge field and couple to the tensor models at any dimensions.
Gauging should be important for understanding the actual AdS/CFT correspondence.
In 1 dimension, the gauging singles out from the spectrum all non-singlet states from the Hilbert states. There have been many attempts to understand of the structure of the tensorial quantum mechanics of Majorana fermions from numerical and analytical calculations [48][49][50][51]. These gave some interesting results, such as the structure of the spectrum of the matrix quantum mechanics and the importance of the discrete symmetries for explaining huge degenaracies of the spectra. Still, the general impact of gauging of the tensorial theory is not clear and demands a new approach. Here, we will give some comments of the combinatorial character and study how the gauging of N = 2 theory, studied in the previous section, changes.
In 3 dimensions one can gauge a theory by adding a Chern-Simons term instead of the usual Yang-Mills term where W (Φ abc ) is the same as in the (4.1), T a are the generators of the group O(N ) × O(N ) × O(N ), and Γ α are vector superfields that have a gauge potential A a bµ as one of the components. If one rewrites the kinetic term for the gauge field in terms of usual components, he will get a usual Chern-Simons theory. Since the theory is gauge invariant, we can choose an axial gauge to simplifty the action 4 A a b3 = 0, which eliminates the non-linear term from the theory and the Fadeev-Popov ghosts decouple from the theory. Therefore the A a b1 , A a b2 can be integrated out to get an effective potential. For example, such a term appears in the 4 I would like to thank S.Prakash for the suggested argument.
action (5.6) which can be considered as a non-local pillow operator with the wrong scaling, because the level of CS action usually scales as k = λN . Therefore some diagrams would have large N factor and diverge in the large N limit. To fix it we should consider the unusual scaling for the CS level k = λN 2 .
One can check that only specific Feynman propagators containing the non-local vertex (5.6) contribute in the large N limit [44]. Namely only snail diagrams contribute in the large N limit and usually are equal to zero by dimensional regularization for massless fields.
Therefore, one can suggest that the gauge field in the large N limit does not get any large corrections and does not change the dynamics of the theory. This argument being purely combinatorial should be applied for any theory coupled to the CS action.
We can confirm this argument by direct calculation of the dimensions of the fields in the expansion for the N = 2 supertensor model at two-loops and see whether the dimensions of the fields gets modified. The beta-functions for a general N = 2 theory coupled to a CS action was considered in the paper [47] and have the following form at finite N where γ Φ is the same as in the equation (

A Supersymmetry in 3 dimensions
In this section we will introduce the notations and useful identities for the N = 1 supersymmetric theories in 3 dimensions. We will mostly follow the lectures [35]. The Lorentz group in 3 dimensions is SL(2, R); that is a group of all unimodular real matrices of dimension 2. The gamma matrices can be chosen to be real There is no γ 5 matrix, so we can't split the spinor representation into small Weyl ones.
Because of this, the smallest spinor representation is 2 dimensional and real. It is endowed with a scalar product defined as ξη = ξ α η α = iξ α γ 0 αβ η β , θ 2 = 1 2θ θ. (A.2) Because of these facts, the N = 1 superspace, in addition to the usual space-time coordinates, will include two real Grassman variables θ ± . The fields on the superspace can be decomposed in terms of fields in the usual Minkowski space. For instance, a scalar superfield (that is our major interest) has the following decomposition Φ(x, θ α ) = φ(x) +θψ(x) + θ 2 F (x). (A.3) As usual, the algebra supersymmetry in superspace can be realized via the derivatives that act on the superfields (A.3) and mix different components where ∂ µ stands for differentiation with respect to the usual space-time variables, and ∂ α for the anticommuting ones. One can define a superderivative that anticommutes with supersymmetry generators, and therefore preserves the supersymmetry Out of these ingredients, namely (A.3),(A.5), we can build an explicit version of a supersymmetric Lagrangian. For example, we can consider the following Lagrangian where the integral over Grassman variables is defined in the usual way with the normalization d 2 θθθ = 1. Writing out the explicit form of (A.6) we get S = d 3 x 1 2 (∂ µ φ) 2 + iψ α γ µ αβ ∂ µ ψ β + F 2 + W (φ)F + W (φ)ψ 2 . (A.7) The field F does not have a kinetic term, and therefore is not dynamical and can be integrated out (that we will not do). For a further investigation we have to develop the technique of super Feynman graphs. We start with considering the partition function of the theory (A.6) The last integral is gaussian and therefore can be evaluated and is equal to From this one can recover the usual Feynman diagrammatic technique, where the vertex is taken from the superpotential W (Φ) rather than the integrated version, and the propagator is defined as Φ(x 1 , θ 1 )Φ(x 2 , θ 2 ) = 1 D 2 δ 2 (θ 1 − θ 2 ) = D 2 δ 2 (θ 1 − θ 2 ), (A. 10) which can be calculated by double differentiation of the partition function (A.8), and the operator is the usual laplacian.