About two years ago I have got published a paper of mine on PRD (see here). Generally this is a reason of great accomplishment for us working in the field of high-energy and nuclear physics as well. What I would like to show you in this post is that the referee, the editor and myself were all right and I will do this with a presentation at an undergraduate level, the one of American Journal of Physics if you want. The relevance of this paper is that I modified the Bender et al. method (see here) to make it work properly and finally obtain sensible results out of a strongly coupled quantum field theory.
A quantum field theory can be solved if we are able to solve the corresponding Heisenberg equations of motion. These equations, in the classical limit, are generally solvable in special cases (e.g. free particle) but the classical limit may be used to extract non perturbative results as happens for the WKB approximation. So, let us write down the equation for a massless quartic scalar field theory:
As strange as may seem, this equation admits an exact solution that we can write down as
when the following dispersion relation does hold
Here 
![{\it sn}(u,i) = \frac{2\pi}{K(i)}\sum_{n=0}^\infty\frac{(-1)^ne^{-(n+\frac{1}{2})\pi}}{1+e^{-(2n+1)\pi}}\sin\left[(2n+1)\frac{\pi u}{2K(i)}\right]](https://s0.wp.com/latex.php?latex=%7B%5Cit+sn%7D%28u%2Ci%29+%3D+%5Cfrac%7B2%5Cpi%7D%7BK%28i%29%7D%5Csum_%7Bn%3D0%7D%5E%5Cinfty%5Cfrac%7B%28-1%29%5Ene%5E%7B-%28n%2B%5Cfrac%7B1%7D%7B2%7D%29%5Cpi%7D%7D%7B1%2Be%5E%7B-%282n%2B1%29%5Cpi%7D%7D%5Csin%5Cleft%5B%282n%2B1%29%5Cfrac%7B%5Cpi+u%7D%7B2K%28i%29%7D%5Cright%5D&bg=ffffff&fg=333333&s=0&c=20201002)
being 
that is a superposition of plane waves, i.e. we are describing a spectrum of free particles. But what is the mass spectrum of these particles? From the above solution we have to put our particles to rest. We put 
that is the one of an harmonic oscillator. The theory is trivial. Very easy, isn’t it?
Now, we put 
and we get the following table
to be compared with the one of Meyer for Yang-Mills in D=3+1
Still better are the Teper et al. values
This is really shocking. We have a satisfactory correspondence between our quartic scalar field theory and the spectrum of Yang-Mills theory in D=3+1. Indeed, we have already seen that lattice results about propagators seems to indicate that Yang-Mills theory is trivial. But this connection with the quartic scalar field turns out to be really unexpected. Indeed we will see that the solution we have found is also a solution to the classical Yang-Mills field.
The Gauge Connection 
[…] equations of motion of Yang-Mills theory. Indeed, as already seen for the quartic scalar field (see here), one can show, with undergraduate level arguments, that Yang-Mills theory has a massive exact […]
[…] I devised to solve Dirac equation using the massive solution of the massless scalar field given here and […]