I'm having hard times in finding the Feynman rules for the Gross-Neveu model in two dimensions: $$ \mathcal L = \bar \psi^i (i \gamma^{\mu} \partial_{\mu} ) \psi^i + g^2 (\bar \psi^i \psi^i )^2 \; ,$$ where $i = 1, \ldots, N$ and $\mu = 1, 2$. Moreover $\gamma^0 = \sigma^2$ and $\gamma^1 = i\sigma^1$ where $\sigma^i$ are the Pauli matrices.
while the 4-point diagram should be:
Question: How are the $\epsilon$'s defined? Where they come from?
Moreover, since any propagator brings two Dirac index, one at the beginning and one at the end of the propagator, there shouldn't be another Dirac index at the center of the four-point diagram?