# 1 Introduction

## 1 Introduction

Despite the success of quantum field theory in the description of the behaviour of elementary particles in the perturbative regime of interactions, it still remains a challenge to understand the non-perturbative domain satisfactorily. One of the methods which has gained attention in this regard in recent years is the study of Schwinger-Dyson equations (SDEs) . Despite the difficulties involved in finding a non-perturbative truncation of these equations, this approach has been very successful in addressing issues like dynamical mass generation for fundamental fermions when they are involved in sufficiently strong interactions . Moreover, recent attempts, e.g.,  to improve the reliability of the approximations used have increased the credibility of the results obtained through such studies.

Application of Schwinger-Dyson formalism to supersymmetric (SUSY) models has been less extensive. In supersymmetric Quantum Electrodynamics (SQED), based upon the arguments of non-renormalization theorem and gauge invariance, it is expected to be impossible to obtain dynamical mass generation for fermions  though some studies  argue that it is probably possible to break chiral symmetry dynamically in SQED. We postpone the study of SQED for a future work. In this paper, we take the simplest SUSY model, i.e the Wess-Zumino model and attempt to solve the corresponding Schwinger-Dyson equations for the fermion and scalar propagators. We believe this exercise will provide us with a deeper insight into how the role of supersymmetry in the context of dynamical mass generation translates into the language of Schwinger-Dyson equations.

Such a study should provide us with a better starting point for more complicated SUSY theories such as SQED and SQCD. In the latter theory, a need also exists to further explore connections between the Holomorphic approach and that of the Schwinger-Dyson equations .

We first study the Yukawa model with one real scalar and one Majorana fermion, which can be considered as a truncated Wess-Zumino model. We discuss this model in some detail for two reasons, first being that it is interesting in its own right because, after all, it is Yukawa interactions which are responsible for giving masses to fermions in the Standard Model (SM). Secondly, extending the Yukawa model by doubling the scalar degrees of freedom provides us with a clear understanding of how supersymmetry works. We use the quenched approximation. Keeping in mind the perturbative expansion of the 3-point vertex beyond the lowest order and its transformation under charge conjugation symmetry, we propose an ansatz for the full effective vertex. One of the advantages of using this vertex is that the equations for the mass function $`(p^2)`$ and the wavefunction renormalization $`F(p^2)`$ decouple completely in the neighbourhood of the critical coupling, $`\alpha _c`$ above which mass is generated for the fermions, and partly above it. We solve both the equations to find analytical expressions for $`F(p^2)`$ and the anomalous mass dimensions in the neighbourhood of $`\alpha _c`$. The results show that non-perturbative interaction of fermions with fundamental scalars can give masses to fermions in a dynamical way provided the interaction is strong enough. We use numerical calculation to draw Euclidean mass of the fermions as a function of the coupling, and confirm that it obeys Miransky scaling. We also evaluate $`F(p^2)`$ numerically. We then extend the particle spectrum by doubling the number of scalars and imposing relations for the couplings that define the Wess-Zumino model. Due to the presence of the additional symmetry, we are able to extract useful information beyond the rainbow approximation.

## 2 The Yukawa Model

Consider a massless Lagrangian with one Majorana fermion and one real scalar interacting with each other through a $`\gamma _5`$-type interaction:

$`={\displaystyle \frac{1}{2}}(_\mu A)^2+{\displaystyle \frac{1}{2}}i(\overline{\psi }\gamma ^\mu _\mu \psi ){\displaystyle \frac{1}{2}}g^2A^4+ig\overline{\psi }\gamma _5\psi A.`$ (1)   
The corresponding Schwinger-Dyson equation for the fermion propagator, $`S_F(p)`$, is displayed in Fig. 1. Motivated by the success of the quenched approximation in QED and QCD, we neglect the fermion loops. Moreover, as a first step towards truncating the infinite set of Schwinger-Dyson equations, we drop all 4-point functions. The full scalar propagator can then be replaced by its bare counterpart. Using Feynman rules, the Schwinger-Dyson equation can be written as:

$`iS_F^1(p)=iS_F^{0^1}(p){\displaystyle \frac{d^4k}{(2\pi )^4}(2g\gamma _5)(iS_F(k))(2g\gamma _5\mathrm{\Gamma }_A(k,p))(\frac{i}{q^2})},`$ (2)   
where $`q=kp`$ and $`S_F(p)`$ can be expressed in terms of two Lorentz scalar functions, $`F(p^2)`$, the wavefunction renormalization and $`(p^2)`$, the mass function, so that

$`S_F(p)`$ $`=`$ $`{\displaystyle \frac{F(p^2)}{\overline{)}p(p^2)}}.`$ (3)   
The bare propagator $`S_F^0(k)=1/(\text{/}p)`$, where the bare mass has been taken to be zero. We can project out equations for $`F(p^2)`$ and $`(p^2)`$ by taking the trace of Eq. (2) having multiplied by $`\overline{)}p`$ and 1 in turn. On Wick rotating to Euclidean space,

$`F(p^2)`$ $`=`$ $`1{\displaystyle \frac{\alpha }{\pi ^3}}{\displaystyle \frac{1}{p^2}}{\displaystyle d^4k\frac{F(k^2)F(p^2)}{k^2+^2(k^2)}\frac{kp}{q^2}\mathrm{\Gamma }_A(k,p)}`$ (4)   
$`(p^2)`$ $`=`$ $`{\displaystyle \frac{\alpha }{\pi ^3}}{\displaystyle d^4k\frac{(k^2)}{k^2+^2(k^2)}\frac{F(k^2)F(p^2)}{q^2}\mathrm{\Gamma }_A(k,p)}`$ (5)   
where $`\alpha =g^2/4\pi `$.

It is here that we cannot proceed any further without making an ansatz for $`\mathrm{\Gamma }_A(k,p)`$. Any ansatz for the 3-point vertex must fulfill at least the following requirements:

* Perturbatively, we must have $`\mathrm{\Gamma }_A(k,p)=1+𝒪(g^2)`$.
* It must be symmetric in $`k`$ and $`p`$.

Moreover, as the SDEs relate the 2-point function with the 3-point function, the expression for the full vertex is expected to involve functions $`F(p^2)`$ or/and $`(p^2)`$.

The commonly used ansatz in the non-perturbative study of the SDEs is the bare vertex ansatz. For the Yukawa Model under discussion, it implies $`\mathrm{\Gamma }_A(k,p)=1`$. It agrees with the lowest order perturbation theory. The only truncation of the complete set of Schwinger-Dyson equations known so far that avoids any assumptions other than the smallness of the coupling at every level of this approximation is the perturbation theory. Therefore, it is natural to assume that physically meaningful solutions of the Schwinger-Dyson equations must agree with perturbative results in the weak coupling regime. It requires, e.g., that every non-perturbative ansatz chosen for the full vertex must reduce to its perturbative counterpart when the interactions are weak. Bare vertex fulfills this requirement to the lowest order in perturbation theory. Any other vertex which fulfills this condition and does not violate other requirements is at least as good as the bare vertex. One of the simplest non-perturbative vertices can be constructed by realizing that Eq.(4) yields the following expansion of $`F(p^2)`$ in perturbation theory:

$`F(p^2)=1+𝒪(\alpha )`$ (6)   
Therefore, a simple candidate for the 3-point vertex can be written as:

$`\mathrm{\Gamma }_A(k,p)={\displaystyle \frac{1}{F(k^2)F(p^2)}}.`$ (7)   
Perturbatively, it gives

$`\mathrm{\Gamma }_A(k,p)={\displaystyle \frac{1}{[1+𝒪(\alpha )][1+𝒪(\alpha )]}}=1+𝒪(\alpha )`$ (8)   
Therefore, to the lowest order in perturbation theory, our vertex ansatz reduces to the bare vertex.

It is exceedingly complicated to ensure that at higher orders, the non-perturbative vertex reduces to its perturbative counterpart in the weak coupling regime. We do not aim at it in this paper. However, we demonstrate that even up to next to lowest order, i.e, to $`𝒪(\alpha )`$, our ansatz is correct to the extent that both the ansatz and the real vertex have the logarithmically divergent behaviour in the ultraviolet regime.

Fig. 2 represents the perturbative expansion of the 3-point fermion-scalar vertex to $`𝒪(\alpha )`$. Using Feynman rules, we can write it as follows:

$`2g\gamma _5\mathrm{\Gamma }_A=2g\gamma _5+{\displaystyle \frac{d^4w}{(2\pi )^4}(2g\gamma _5)iS_F(pw)(2g\gamma _5)iS_F(kw)(2g\gamma _5)\frac{i}{w^2}}`$ (9)   
which can be simplified to:

$`\mathrm{\Gamma }_A`$ $`=`$ $`16\pi i\alpha ^2\left[\overline{)}k\overline{)}pJ^{(0)}(\overline{)}k\gamma ^\nu +\gamma ^\nu \overline{)}p)J_\nu ^{(1)}+K^{(0)}\right]`$ (10)   
where

$`J^{(0)}`$ $`=`$ $`{\displaystyle d^4w\frac{1}{w^2(pw)^2(kw)^2}}`$ (11)   
$`J_\mu ^{(1)}`$ $`=`$ $`{\displaystyle d^4w\frac{w_\mu }{w^2(pw)^2(kw)^2}}`$ (12)   
$`K^{(0)}`$ $`=`$ $`{\displaystyle d^4w\frac{1}{(pw)^2(kw)^2}}`$ (13)   
The exact analytical expressions for these three integrals are known . They involve basic functions of momenta $`k`$ and $`p`$ and a spence function. We believe that it is highly non-trivial to construct a non-perturbative vertex which reduces to this complicated form in the weak coupling regime. However, asymptotic behaviour can be reproduced to some extent. Simple power counting reveals that the integrals $`J^{(0)}`$ and $`J_\mu ^{(1)}`$ are perfectly well-behaved in the ultraviolet regime. However, $`K^{(0)}`$ is logarithmically divergent. We now show that our vertex ansatz also exhibits this behaviour. Using the fact that perturbatively $`(p^2)=0`$, we can re-write Eq. (4) as follows:

$`F(p^2)`$ $`=`$ $`1{\displaystyle \frac{\alpha }{\pi ^3}}{\displaystyle \frac{1}{p^2}}{\displaystyle d^4k\frac{F(k^2)F(p^2)}{k^2}\frac{kp}{q^2}}`$ (14)   
where we have employed Eq. (6), and the Feynman rule for the vertex. Carrying out angular and radial integrations respectively and retaining the leading log terms, we get

$`F(p^2)`$ $`=`$ $`1+{\displaystyle \frac{\alpha }{2\pi }}\mathrm{ln}{\displaystyle \frac{p^2}{\mathrm{\Lambda }^2}}`$ (15)   
Therefore, the proposed vertex ansatz can be written perturbatively as follows

$`\mathrm{\Gamma }_A(k,p)={\displaystyle \frac{1}{F(k^2)F(p^2)}}=1+{\displaystyle \frac{\alpha }{\pi }}\mathrm{ln}{\displaystyle \frac{k^2p^2}{\mathrm{\Lambda }^2}}+𝒪(\alpha ^2)`$ (16)   
which is logarithmically divergent in the ultraviolet regime just as the real vertex to $`𝒪(\alpha )`$. Therefore, perturbatively our vertex ansatz is more realistic than the bare vertex.

An added advantage of using the proposed vertex ansatz is that Eq.(5) can be solved independently of Eq.(4). Therefore, it serves a purpose similar to that of Mandlestam’s choice  of the 3-gluon vertex in studying the Schwinger-Dyson equation of the gluon propagator.

As the unknown functions $`F`$ and $``$ do not depend upon the angle between $`k`$ and $`p`$, we can perform angular integration to arrive at

$`F(p^2)`$ $`=`$ $`1{\displaystyle \frac{\alpha }{2\pi }}{\displaystyle 𝑑k^2\frac{1}{k^2+^2(k^2)}\left[\frac{k^4}{p^4}\theta (p^2k^2)+\theta (k^2p^2)\right]}`$ (17)   
$`(p^2)`$ $`=`$ $`{\displaystyle \frac{\alpha }{\pi }}{\displaystyle 𝑑k^2\frac{(k^2)}{k^2+^2(k^2)}\left[\frac{k^2}{p^2}\theta (p^2k^2)+\theta (k^2p^2)\right]}.`$ (18)   
Such equations are known to have a non-trivial solution for the mass function above a critical value of the coupling $`\alpha =\alpha _c`$. In the neighbourhood of the critical coupling, when the generated mass is still small, we can put $`^2=0`$. Then Eqs. (17,18) decouple from each other completely. The leading log solution for $`F(p^2)`$ is then:

$`F(p^2)`$ $`=`$ $`1+{\displaystyle \frac{\alpha }{2\pi }}\mathrm{ln}{\displaystyle \frac{p^2}{\mathrm{\Lambda }^2}}.`$ (19)   
As for the mass function, multiplicative renormalizability demands a solution of the type $`M(p^2)(p^2)^s`$. Substituting this in Eq. (18) and performing radial integration, we find

$`s`$ $`=`$ $`{\displaystyle \frac{1}{2}}\pm {\displaystyle \frac{1}{2}}\sqrt{1{\displaystyle \frac{\alpha }{\alpha _c}}}`$ (20)   
where $`\alpha _c=\pi /4`$. For $`\alpha >\alpha _c`$ the solution of the mass function enters the complex plane indicating that a phase transition has taken place from perturbative to non-perturbative solution corresponding to the dynamical generation of mass. Numerically, above $`\alpha _c`$, we solve Eq. (18) in a two-step process. We first use the iterative method to get close to the solution and then refine the answer by converting the integral equation into a set of simultaneous nonlinear equations to be solved by Newton-Raphson method. In Fig. 3., we have drawn the Euclidean mass $`M`$(which can be taken to be $`(0)`$) as a function of the coupling $`\alpha `$. We see that it obeys Miransky scaling law and can be fitted to the form

$`{\displaystyle \frac{M}{\mathrm{\Lambda }}}`$ $`=`$ $`\mathrm{exp}\left[{\displaystyle \frac{A}{\sqrt{\frac{\alpha }{\alpha _c}1}}}+B\right]`$ (21)   
very well by the choice $`A=0.97\pi `$ and $`B=1.45`$. These numbers are close to the ones found in  although the value of the critical coupling is of course different. The slight mismatch in the values of $`A`$ and $`B`$ is due to the fact that in the logarithmic grid of momenta, we choose 30 points per decade and do not extrapolate the result to an infinite number of points, an exercise carried out in . We also compute $`F(p^2)`$ for various values of $`\alpha `$ and find that the closer we approach $`\alpha _c`$, where the generated mass is still small, starting from a larger value of $`\alpha `$, the numerical result gets closer and closer to the analytical result as expected, Fig. 4.

## 3 The Wess-Zumino Model

We now extend the particle spectrum by doubling the number of scalars to discuss the massless Wess-Zumino model, characterized by the following Lagrangian:

$`={\displaystyle \frac{1}{2}}(_\mu A)^2+{\displaystyle \frac{1}{2}}(_\mu B)^2+{\displaystyle \frac{1}{2}}(i\overline{\psi }\gamma ^\mu _\mu \psi ){\displaystyle \frac{1}{2}}g^2(A^2+B^2)^2g\overline{\psi }(BiA\gamma _5)\psi .`$ (22)   
The Schwinger-Dyson equation for the fermion propagator in this model is depicted in Fig. 5. Before we embark on solving this equation, we must re-address the validity of the ansatz for the 3-point vertex which we made for the Yukawa case. To the one loop level in perturbation theory, one now has contributions from both the scalars as depicted in Fig. 6. Therefore, e.g., for scalar $`A`$, we can write

$`2g\gamma _5\mathrm{\Gamma }_A`$ $`=`$ $`2g\gamma _5+{\displaystyle \frac{d^4w}{(2\pi )^4}(2g\gamma _5)iS_F(pw)(2g\gamma _5)iS_F(kw)(2g\gamma _5)\frac{i}{w^2}}`$ (23)   
$`+`$ $`{\displaystyle \frac{d^4w}{(2\pi )^4}(2ig)iS_F(pw)(2g\gamma _5)iS_F(kw)(2ig)\frac{i}{w^2}}=0.`$   
The same is the case for $`\mathrm{\Gamma }_B`$, i.e., in the presence of both the scalars $`A`$ and $`B`$, none of the 3-point vertices gets modified at $`𝒪(\alpha )`$ in perturbation theory. Therefore, in the Wess-Zumino case, it is more reasonable to use the bare vertex instead of the ansatz earlier made. Though one would now expect to solve coupled integral equations for $`F(p^2)`$ and $`(p^2)`$, a miraculous cancellation of terms takes place as evident from the following Schwinger-Dyson equation for the fermion propagator:

$`{\displaystyle \frac{\overline{)}p(p^2)}{iF(p^2)}}`$ $`=`$ $`{\displaystyle \frac{\overline{)}p}{i}}{\displaystyle \frac{\alpha }{\pi ^3}}{\displaystyle d^4k\left[\frac{F(k^2)}{\overline{)}k(p^2)}\frac{1}{q^2}\right]}+{\displaystyle \frac{\alpha }{\pi ^3}}{\displaystyle d^4k\left[\gamma _5\frac{F(k^2)}{\overline{)}k(p^2)}\gamma _5\frac{1}{q^2}\right]}.`$   
Taking the trace of this equation, we get

$`(p^2)`$ $`=`$ $`0.`$   
As the cancellation of terms takes place at the very beginning, it is easy to see that dynamical mass generation will remain an impossibility for the full vertex and the full scalar propagator as long as they are identical for both the scalars. The vertex corrections for $`A`$ and $`B`$ have been proven to be equal up to $`𝒪(\alpha ^2)`$ . We shall shortly see that the same is true for the full scalar propagator at least up to $`𝒪(\alpha )`$. This is in accordance with the arguments based on non-renormalization theorem. SUSY plays a role in providing same number of bosonic and fermionic degrees of freedom. We have seen from the case of the Yukawa Model that without this equality, it will not be possible to prevent dynamical mass generation. Secondly, SUSY imposes relations on couplings of the two scalars with the fermions. This relationship is crucial in preventing dynamical mass generation.

As far as wavefunction renormalization $`F(p^2)`$ is concerned, its leading log behaviour gets modified slightly, by the inclusion of the other scalar, to:

$`F(p^2)`$ $`=`$ $`1+{\displaystyle \frac{\alpha }{\pi }}\mathrm{ln}{\displaystyle \frac{p^2}{\mathrm{\Lambda }^2}}.`$ (25)   
Although it is an interesting conclusion in its own right that supersymmetry prevents dynamical mass generation for fermions in the Wess-Zumino model, another important issue to probe will be whether supersymmetry itself remains intact, i.e., whether the scalars can also be kept massless. This is what we discuss now. The Schwinger-Dyson equation for the scalar (for example $`A`$) has been depicted in Fig. 7. A scalar propagator, unlike a fermion, needs only one unknown function to describe it. But we shall prefer to split it into two parts and write the full scalar propagator as follows:

$`S_A(p)`$ $`=`$ $`{\displaystyle \frac{F_A(p^2)}{p^2_A^2(p^2)}}.`$ (26)   
The non-zero value of the mass function $`_A(p^2)`$ will be responsible for shifting the pole from $`p^2=0`$ to some finite value, generating the mass for the scalar dynamically. $`F_A(p^2)`$ on the other hand is the scalar wavefunction renormalization. The SD-equation for the scalar propagator in Euclidean space can now be written as:

$`{\displaystyle \frac{p^2+_A^2(p^2)}{F_A(p^2)}}`$ $`=`$ $`p^2+{\displaystyle \frac{3\alpha }{2\pi ^3}}{\displaystyle d^4k\frac{F_A(k^2)}{k^2+_A^2(k^2)}}+{\displaystyle \frac{\alpha }{2\pi ^3}}{\displaystyle d^4k\frac{F_B(k^2)}{k^2+_B^2(k^2)}}`$ (27)   
$`{\displaystyle \frac{2\alpha }{\pi ^3}}{\displaystyle d^4k\frac{F(k^2)F(q^2)}{k^2q^2}\mathrm{\Gamma }_A(k,p)kq}`$   
where we have used the fact that the fermions do not acquire mass. If we want to preserve supersymmetry and do not want the scalars to acquire mass, we must have:

$`_A(p^2)=_B(p^2)=0.`$ (28)   
We are then left with:

$`{\displaystyle \frac{1}{F_A(p^2)}}`$ $`=`$ $`1+{\displaystyle \frac{\alpha }{2\pi ^3p^2}}{\displaystyle \frac{d^4k}{k^2}\left[3F_A(k^2)+F_B(k^2)4\frac{kq}{q^2}F(k^2)F(q^2)\mathrm{\Gamma }_A(k,p)\right]}`$ (29)   
and there is a similar equation for the scalar $`B`$:

$`{\displaystyle \frac{1}{F_B(p^2)}}`$ $`=`$ $`1+{\displaystyle \frac{\alpha }{2\pi ^3p^2}}{\displaystyle \frac{d^4k}{k^2}\left[3F_B(k^2)+F_A(k^2)4\frac{kq}{q^2}F(k^2)F(q^2)\mathrm{\Gamma }_B(k,p)\right]}.`$ (30)   
These equations should yield a solution for $`F_A(p^2)`$ and $`F_B(p^2)`$ such that it does not change the position of the pole for the scalar propagator and that the quadratic divergences cancel. It is well-known that it does happen in perturbation theory to $`𝒪(\alpha )`$. In fact, one can evaluate $`F_A(p^2)`$ and $`F_B(p^2)`$. The leading log expression for these functions to $`𝒪(\alpha )`$ is

$`F_A(p^2)=F_B(p^2)=1+{\displaystyle \frac{\alpha }{\pi }}\mathrm{ln}{\displaystyle \frac{p^2}{\mathrm{\Lambda }^2}}`$ (31)   
which is exactly the same expression as that for $`F(p^2)`$ for the fermion propagator. This result indicates that supersymmetry need not be broken.

## 4 Conclusions

We have studied the Schwinger-Dyson equations for the Yukawa (a scalar interacting with a fermion with a $`\gamma _5`$ type interaction) and the Wess-Zumino models. In the simple Yukawa model, we propose a vertex $`\mathrm{𝑎𝑛𝑠𝑎𝑡𝑧}`$ which we argue should perform better than the bare vertex. In the quenched approximation, we find dynamical mass generation for fermions above a critical value of the coupling $`\alpha _c=\pi /4`$. The generated Euclidean mass obeys Miransky scaling. When we extend this Yukawa model to equate the scalar and fermionic degrees of freedom (Wess-Zumino model), we find that a neat cancellation of terms occurs and there is no mass generation for the fermions. This fact remains true beyond the rainbow approximation and is supported by perturbative calculations available for the 3-point vertex to $`𝒪(\alpha ^2)`$ and of the scalar propagator. This result was expected on the basis of non-renormalization theorem. The two approaches will remain in agreement provided the full vertex and the full scalar propagator as long as they are identical for both the scalars to higher orders in perturbation theory as well.

If supersymmetry has to be preserved, the scalars should also acquire no mass dynamically. Studying the Schwinger-Dyson equations for the scalars, we observe that such a solution is allowed and in fact leads to the wavefunction renormalization function for the scalars which is exactly the same as that for the fermion.

It is more interesting to see the role of supersymmetry in more complicated theories such as SQED. The studies so far carried out in superfield and component formalism seem to arrive at different conclusions. We plan to present our work in this context in a future publication.

Acknowledgements

This work was partly supported by a TWAS-AIC award and CONACYT-SNI (México). AB is also grateful for the hospitality of Instituto de Física, Benemérita Universidad Autónoma de Puebla (BUAP) during his stay there, where a part of the work was done.

Figure Captions

Fig. 1. Schwinger-Dyson equation for the fermion propagator in the Yukawa model. The solid lines represent fermions and the dashed the scalar. The solid dots indicate full, as opposed to bare, quantities.

Fig. 2. One loop perturbative expansion for the vertex in the Yukawa model.

Fig. 3. The dynamically generated mass $`(M/\mathrm{\Lambda })`$ versus the 3-point coupling $`\alpha `$ in the Yukawa model. The critical coupling is $`\alpha _c=\pi /4`$, above which the mass can be seen to be bifurcating away from the chirally symmetric solution. $``$s represent the numerical result and $`+`$s the numerical fit to the form $`M=\mathrm{\Lambda }\mathrm{exp}\left[A/\sqrt{\alpha /\alpha _c1}+B\right]`$ with $`A=0.97\pi `$ and $`B=1.45`$.

Fig. 4. The wavefunction renormalization function $`F(p^2)`$ in the Yukawa model for various values of the coupling $`\alpha `$. The solid line corresponds to the analytical expression $`1\alpha /4\pi +(\alpha /2\pi )\mathrm{ln}(p^2/\mathrm{\Lambda }^2)`$ in case of no mass generation for $`\alpha =0.78`$.

Fig. 5. Schwinger-Dyson equation for the fermion propagator in the Wess-Zumino model.

Fig. 6. One loop perturbative expansion for the vertex in the Wess-Zumino model. We have shown the case for scalar $`A`$. A similar diagram exists for scalar $`B`$.

Fig. 7. Schwinger-Dyson equation for the scalar propagator in the Wess-Zumino model.