I'm currently studying a scalar quantum field theory with a $\lambda\phi^4$ interaction (commonly referred to as $\phi^4$ theory. To study this theory non perturbatively I've written a program in python that simulates the theory in the euclidean (imaginary time) formalism. I'm trying to numerically look at the potential triviality of the theory in 3 + 1 dimensions. The problem I'm having is that I can't find a good definition of the renormalized coupling constant that I can compute with relative ease on the lattice.
On Scholarpedia they have a formula for the renormalized coupling constant in the symmetric phase is given by $$\lambda_r=(Lm_r)^4(3\left\langle M^2\right\rangle ^2 - \left\langle M^4 \right\rangle)/\left\langle M^2 \right\rangle ^2$$ where $$ M = a^4\sum_{i}\phi(x_i),$$ and where $\left\langle X \right\rangle$ deontes the expectation value or mean of said quantity. Also note, $m_r$ is the renormalized mass of the scalar particle.
I'm having a few problems with this coupling constant as it's defined. The first is that it is unreasonably large even for highly coarse lattices with a lattice spacing $a = 1/2$ and with a lattice size $L = 8$. The bigger problem however is that if we're in the the symmetric phase where $\langle\phi\rangle = 0$, it follows that $ M = a^4\sum_{i}\phi(x_i)$ should be zero as well.
So my question is the following; have I messed up somewhere with my definition of the renormalized coupling constant? Has Scholarpedia messed up with their definition? Does anyone have an equally easy to calculate definition of the coupling constant that doesn't have these (apparent) problems?
Thanks