$\beta$ function for the Gross-Neveu model

Question

In the Peskin & Schroeder textbook, the $\beta$-function for the Gross-Neveu model is discussed in problem 12.2. After computing it, I have tried checking my results with some solutions found online. My problem is that they all disagree among each other (something quite recurrent for this book actually). I have basically followed the original paper from Gross and Neveu:

• We start from the Lagrangian $\mathcal{L} = \bar{\Psi}_{i}(i\not{\partial} )\Psi_{i}-g\sigma \bar{\Psi}_i\Psi_i-\frac{1}{2}\sigma^2$.
• We compute the field strength renormalisation for $\Psi$.

For $\Psi: \quad$ $\delta_{Z_{\Psi}}=(-ig)^2\int \frac{d^{d}k}{(2\pi)^d} Tr[\frac{i\not{k}}{k^2}\frac{i}{(k+p)^2}]=0 $

For $\sigma: \quad$ $\delta_{Z_{\sigma}}=-N(-ig)^2\int \frac{d^{d}k}{(2\pi)^d}\frac{i\not{k}}{k^2} \frac{i(\not{k}+\not{p})}{(k+p)^2} =-Ng^2 \int \frac{d^{d}k}{(2\pi)^d}\frac{1}{(k+p)^2} =-Ng^2\frac{i}{4\pi\epsilon}+finite$

• We compute the coupling renormalisation (though in the Gross-Neveu original paper, the authors argue the vertex need not be renormalised. That is a point I don't get)

$\delta_{g}=(-ig)^3 \int \frac{d^{d}k}{(2\pi)^d} Tr[\frac{i(\not{k}+\not{p_1})}{(k+p_1)^2} \frac{i(\not{k}+\not{p_2})}{(k+p_2)^2}] = -ig^3 \int \frac{d^{d}k}{(2\pi)^d}\frac{-2k^2}{(k+p_1)^2 (k+p_2)^2}+...=2ig^3\frac{i}{4\pi\epsilon} $

• We use relation (12.53)

$\beta (g) = M\frac{\partial}{\partial M}\left(-\delta_g +\frac{1}{2}g\sum_i \delta_{Z_i} \right)=\frac{2g^3}{4\pi}+g\frac{-Ng^2}{4\pi}=-\frac{g^3}{4\pi}(N-2)$

It is asymptotically free, but I do fully trust my result. Any comment welcome!

Also, as a second question, the $\beta$ function derived by Gross & Neveu is $\beta(g)=-\frac{g^2N}{2\pi}$. The dependence on g is only $g^2$ and not $g^3$. Do they use a different definition of the $\beta$ function? If so, what is the advantage of this other definition?

Answer 1

I guess you have been comparing Hitoshi Murayama's solution vs. Zhong-Zhi Xianyu's. Both solutions have the right (and essentially the same) approach. Unfortunately, both contain some critical errors here and there.

Xianyu's solution first covers the two-point function before treating the four-point function. This is great. We first get $\gamma$ out of the way.

As for the Feynman rule for the vertex, those $\epsilon$'s should really be $\delta$'s. They are nothing but identity matrices used to track which propagators should be connected to which vertices.

If you see something like $\delta_{\alpha\beta} (\frac{i}{\not{k}})_{\beta\gamma} \delta_{\gamma \delta} (\frac{i}{\not{k}})_{\delta\alpha} $, it is equivalent to $-\mathrm{Tr}(\frac{i}{\not{k}} \frac{i}{\not{k}}) = -2 \frac{1}{k^2}$. (The negative sign is because now we have a Fermion loop, while the factor of 2 is from the trace in two dimensions.)

If you see something like $\delta_{\alpha\beta} (\frac{i}{\not{k}})_{\beta\gamma} \delta_{\gamma \delta} (\frac{i}{\not{k}})_{\delta\epsilon} $, it is equivalent to $\frac{i}{\not{k}} \frac{i}{\not{k}}$, which is then meant to be sandwiched between $\bar{u}(p_1)$ and $u(p_2)$. Compared to the previous situation, there are two important differences. First, the sign is the opposite. Second, we don't have the extra factor of 2.

With this in mind, you should be able to see how to interpret (12.8) in Xianyu's solution. The integrand $(A+B)\frac{i}{\not{k}}(C+D)\frac{i}{\not{k}}$ can be decomposed into four parts.

• $A\frac{i}{\not{k}}C\frac{i}{\not{k}}$: This is basically (58) in itoshi Murayama's solution.
• $A\frac{i}{\not{k}}D\frac{i}{\not{k}}$: This is basically (59) in itoshi Murayama's solution.
• $B\frac{i}{\not{k}}C\frac{i}{\not{k}}$: This is basically (59) in itoshi Murayama's solution.
• $B\frac{i}{\not{k}}D\frac{i}{\not{k}}$: This is basically (61) in itoshi Murayama's solution (t-channel).

N.B.: I am not claiming (58) is error-free. In fact, the factor 3/2 should be 2 instead in the numerator of the second term in the integrand. (See (10.33) in Peskin and Schroeder for reference.) This error propagates to (62) and a bunch of other equations. Once this is fixed, you will get the desired $(N-1)$ instead of $(N-2)$ in the final result.

The vertex Feynmann rule given in Xianyu's solution is correct. As for how to derive it, I would refer to Section 4.7 in Peskin and Schroeder. You can write something like $\langle 1(p') 2(k') | \bar{\psi}_1(x)\psi_1(x)\bar{\psi}_1(x)\psi_1(x) \bar{\psi}_1(y)\psi_1(y) \bar{\psi}_2(y)\psi_2(y) | 1(p) 2(k) \rangle$ and draw lines to connect each of them. (1 and 2 are particle types, while $p$, $k$, etc. are momenta.) Depending on how you draw those lines, this particular one here can either contain a loop or not.

When you expand $(\bar{\psi}_i\psi_i)^2$ in the Langrangian, you get $N + N(N-1)/2$ terms for each vertex. Although a term like $2\bar{\psi}_1\psi_1\bar{\psi}_2\psi_2$ has an extra term of 2 compared to $\bar{\psi}_1\psi_1\bar{\psi}_1\psi_1$. The latter is in fact twice as large because it can be connected to propagators in four different ways.

A few extra notes about Xianyu's solution.

• If you wonder what those $\gamma^\mu \gamma_\mu$ terms are about in e.g. (12.8), it has nothing to do with summation over $\mu$. It represents in fact $(\not{k}+\not{q})(\not{k}+\not{p})$ that needs to be sandwiched between $\bar{u}(q)$ and $u(p)$. Murayama's solution is a lot clearer in this regard.
• The factor $1/2$ in (12.8), (12.9), and (12.10) really should not be there. In fact, a nice feature about this theory, just as Yukawa theory, is that you don't need to worry about symmetry factors.