A method of quantizing the Yang–Mills theory is by functional methods, i.e.
path integrals. One introduces a generating functional for -point functions as :\ Z[j] = \int [\mathrm{d}A]\ \exp\left[- \tfrac{i}{2} \int \mathrm{d}^4x\ \operatorname{tr}\left( F^{\mu \nu}\ F_{\mu \nu}\right) + i\ \int \mathrm{d}^4x\ j^a_\mu(x)\ A^{a\mu}(x) \right]\ , but this integral has no meaning as it is because the potential vector can be arbitrarily chosen due to the
gauge freedom. This problem was already known for quantum electrodynamics but here becomes more severe due to non-abelian properties of the gauge group. A way out has been given by
Ludvig Faddeev and
Victor Popov with the introduction of a
ghost field (see
Faddeev–Popov ghost) that has the property of being unphysical since, although it agrees with
Fermi–Dirac statistics, it is a complex scalar field, which violates the
spin–statistics theorem. So, we can write the generating functional as :\begin{align} Z[j,\bar\varepsilon,\varepsilon] & = \int [\mathrm{d}\ A] [\mathrm{d}\ \bar c] [\mathrm{d}\ c]\ \exp\Bigl\{ i\ S_F\ \left[\partial A, A\right] + i\ S_{gf}\left[\partial A\right] + i\ S_g\left[\partial c, \partial\bar c, c,\bar c, A \right] \Bigr\} \\ &\exp\left\{i\int \mathrm{d}^4x\ j^a_\mu(x)A^{a\mu}(x)+i\int \mathrm{d}^4x\ \left[\bar c^a(x)\ \varepsilon^a(x) + \bar\varepsilon^a(x)\ c^a(x)\right]\right\} \end{align} being :S_F=- \tfrac{1}{2} \int \mathrm{d}^4 x\ \operatorname{tr}\left( F^{\mu \nu}\ F_{\mu \nu} \right)\ for the field, :S_{gf} = -\frac{1}{2\xi} \int \mathrm{d}^4 x\ (\partial\cdot A)^2\ for the gauge fixing and :\ S_g = -\int \mathrm{d}^4 x\ \left(\bar c^a\ \partial_\mu\partial^\mu c^a + g\ \bar c^a\ f^{abc}\ \partial_\mu\ A^{b\mu}\ c^c \right)\ for the ghost. This is the expression commonly used to derive Feynman's rules (see
Feynman diagram). Here we have for the ghost field while fixes the gauge's choice for the quantization. Feynman's rules obtained from this functional are the following These rules for Feynman's diagrams can be obtained when the generating functional given above is rewritten as :\begin{align} Z[j,\bar\varepsilon,\varepsilon] &= \exp\left(-i\ g\int \mathrm{d}^4x\ \frac{\delta}{i\ \delta\ \bar\varepsilon^a(x)}\ f^{abc}\partial_\mu\ \frac{i\ \delta}{\delta\ j^b_\mu(x)}\ \frac{i\ \delta}{\delta\ \varepsilon^c(x)} \right)\\ & \qquad \times \exp\left(-i\ g\int \mathrm{d}^4x\ f^{abc}\partial_\mu\frac{i\ \delta}{\delta\ j^a_\nu(x)}\frac{i\ \delta}{\delta\ j^b_\mu(x)}\ \frac{i\ \delta}{\delta\ j^{c\nu}(x)}\right)\\ & \qquad \qquad \times \exp\left(-i\ \frac{g^2}{4}\int \mathrm{d}^4x\ f^{abc}\ f^{ars}\frac{i\ \delta}{\delta\ j^b_\mu(x)}\ \frac{i\ \delta}{\delta\ j^c_\nu(x)}\ \frac{\ i\delta}{\delta\ j^{r\mu}(x)} \frac{i\ \delta}{\delta\ j^{s\nu}(x)} \right) \\ & \qquad \qquad \qquad \times Z_0[j,\bar\varepsilon,\varepsilon] \end{align} with : Z_0[j,\bar\varepsilon,\varepsilon] = \exp \left( -\int \mathrm{d}^4x\ \mathrm{d}^4y\ \bar\varepsilon^a(x)\ C^{ab}(x-y)\ \varepsilon^b(y) \right)\exp \left( \tfrac{1}{2} \int \mathrm{d}^4x\ \mathrm{d}^4y\ j^a_\mu(x)\ D^{ab\mu\nu}(x-y)\ j^b_\nu(y) \right)\ being the generating functional of the free theory. Expanding in and computing the
functional derivatives, we are able to obtain all the -point functions with perturbation theory. Using
LSZ reduction formula we get from the -point functions the corresponding process amplitudes,
cross sections and
decay rates. The theory is renormalizable and corrections are finite at any order of perturbation theory. For quantum electrodynamics the ghost field decouples because the gauge group is abelian. This can be seen from the coupling between the gauge field and the ghost field that is \ \bar c^a\ f^{abc}\ \partial_\mu A^{b\mu}\ c^c ~. For the abelian case, all the structure constants \ f^{abc}\ are zero and so there is no coupling. In the non-abelian case, the ghost field appears as a useful way to rewrite the quantum field theory without physical consequences on the observables of the theory such as cross sections or decay rates. One of the most important results obtained for Yang–Mills theory is
asymptotic freedom. This result can be obtained by assuming that the
coupling constant is small (so small nonlinearities), as for high energies, and applying
perturbation theory. The relevance of this result is due to the fact that a Yang–Mills theory that describes strong interaction and asymptotic freedom permits proper treatment of experimental results coming from
deep inelastic scattering. To obtain the behavior of the Yang–Mills theory at high energies, and so to prove asymptotic freedom, one applies perturbation theory assuming a small coupling. This is verified
a posteriori in the
ultraviolet limit. In the opposite limit, the infrared limit, the situation is the opposite, as the coupling is too large for perturbation theory to be reliable. Most of the difficulties that research meets is just managing the theory at low energies. That is the interesting case, being inherent to the description of hadronic matter and, more generally, to all the observed bound states of gluons and quarks and their confinement (see
hadrons). The most used method to study the theory in this limit is to try to solve it on computers (see
lattice gauge theory). In this case, large computational resources are needed to be sure the correct limit of infinite volume (smaller lattice spacing) is obtained. This is the limit the results must be compared with. Smaller spacing and larger coupling are not independent of each other, and larger computational resources are needed for each. As of today, the situation appears somewhat satisfactory for the hadronic spectrum and the computation of the gluon and ghost propagators, but the
glueball and
hybrids spectra are yet a questioned matter in view of the experimental observation of such exotic states. Indeed, the resonance is not seen in any of such lattice computations and contrasting interpretations have been put forward. This is a hotly debated issue. == Open problems ==