Yang–Mills theory

Open problem in physics: Yang–Mills theory in the non-perturbative regime: The equations of Yang–Mills remain unsolved at energy scales relevant for describing atomic nuclei. How does Yang–Mills theory give rise to the physics of nuclei and nuclear constituents? (more open problems in physics)

Yang–Mills theory is a gauge theory based on the SU(N) group, or more generally any compact, semi-simple Lie group. Yang–Mills theory seeks to describe the behavior of elementary particles using these non-Abelian Lie groups and is at the core of the unification of the electromagnetic and weak forces (i.e. U(1) × SU(2)) as well as quantum chromodynamics, the theory of the strong force (based on SU(3)). Thus it forms the basis of our understanding of particle physics, the Standard Model.

History and theoretical description

In a private correspondence, Wolfgang Pauli formulated in 1953 a six-dimensional theory of Einstein's field equations of general relativity, extending the five-dimensional theory of Kaluza, Klein, Fock and others to a higher-dimensional internal space. However, there is no evidence that Pauli developed the Lagrangian of a gauge field or the quantization of it. Because Pauli found that his theory "leads to some rather unphysical shadow particles”, he refrained from publishing his results formally. Although Pauli did not publish his six-dimensional theory, he gave two talks about it in Zürich. Recent research shows that an extended Kaluza–Klein theory is in general not equivalent to Yang–Mills theory, as the former contains additional terms.

In early 1954, Chen Ning Yang and Robert Mills extended the concept of gauge theory for abelian groups, e.g. quantum electrodynamics, to nonabelian groups to provide an explanation for strong interactions. The idea by Yang–Mills was criticized by Pauli, as the quanta of the Yang–Mills field must be massless in order to maintain gauge invariance. The idea was set aside until 1960, when the concept of particles acquiring mass through symmetry breaking in massless theories was put forward, initially by Jeffrey Goldstone, Yoichiro Nambu, and Giovanni Jona-Lasinio.

This prompted a significant restart of Yang–Mills theory studies that proved successful in the formulation of both electroweak unification and quantum chromodynamics (QCD). The electroweak interaction is described by SU(2) × U(1) group while QCD is an SU(3) Yang–Mills theory. The electroweak theory is obtained by combining SU(2) with U(1), where quantum electrodynamics (QED) is described by a U(1) group, and is replaced in the unified electroweak theory by a U(1) group representing a weak hypercharge rather than electric charge. The massless bosons from the SU(2) × U(1) theory mix after spontaneous symmetry breaking to produce the 3 massive weak bosons, and the photon field. The Standard Model combines the strong interaction, with the unified electroweak interaction (unifying the weak and electromagnetic interaction) through the symmetry group SU(2) × U(1) × SU(3). In the current epoch the strong interaction is not unified with the electroweak interaction, but from the observed running of the coupling constants it is believed they all converge to a single value at very high energies.

Phenomenology at lower energies in quantum chromodynamics is not completely understood due to the difficulties of managing such a theory with a strong coupling. This may be the reason why confinement has not been theoretically proven, though it is a consistent experimental observation. Proof that QCD confines at low energy is a mathematical problem of great relevance, and an award has been proposed by the Clay Mathematics Institute for whoever is also able to show that the Yang–Mills theory has a mass gap and its existence.

Mathematical overview

Yang–Mills theories are a special example of gauge theory with a non-abelian symmetry group given by the Lagrangian

$\mathcal{L}_\mathrm{gf} = -\frac{1}{2}\operatorname{Tr}(F^2)=- \frac{1}{4}F^{a\mu \nu} F_{\mu \nu}^a$

with the generators of the Lie algebra corresponding to the F-quantities (the curvature or field-strength form) satisfying

$\operatorname{Tr}(T^aT^b)=\frac{1}{2}\delta^{ab},\quad [T^a,T^b]=if^{abc}T^c$

and the covariant derivative defined as

$D_\mu=I\partial_\mu-igT^aA^a_\mu$

where I is the identity for the group generators, $A^a_\mu$ is the vector potential, and g is the coupling constant. In four dimensions, the coupling constant g is a pure number and for a SU(N) group one has $a,b,c=1\ldots N^2-1.$

The relation

$F_{\mu \nu}^a = \partial_\mu A_\nu^a-\partial_\nu A_\mu^a+gf^{abc}A_\mu^bA_\nu^c$

can be derived by the commutator

$[D_\mu, D_\nu] = -igT^aF_{\mu\nu}^a.$

The field has the property of being self-interacting and equations of motion that one obtains are said to be semilinear, as nonlinearities are both with and without derivatives. This means that one can manage this theory only by perturbation theory, with small nonlinearities.

Note that the transition between "upper" ("contravariant") and "lower" ("covariant") vector or tensor components is trivial for a indices (e.g. $f^{abc}=f_{abc}$), whereas for μ and ν it is nontrivial, corresponding e.g. to the usual Lorentz signature, $\eta_{\mu \nu }={\rm diag}(+---)$.

From the given Lagrangian one can derive the equations of motion given by

$\partial^\mu F_{\mu\nu}^a+gf^{abc}A^{\mu b}F_{\mu\nu}^c=0.$

Putting $F_{\mu\nu}=T^aF^a_{\mu\nu}$, these can be rewritten as

$(D^\mu F_{\mu\nu})^a=0.$

A Bianchi identity holds

$(D_\mu F_{\nu \kappa})^a+(D_\kappa F_{\mu \nu})^a+(D_\nu F_{\kappa \mu})^a=0$

which is equivalent to the Jacobi identity

$[D_{\mu}, [D_{\nu},D_{\kappa}]]+[D_{\kappa},[D_{\mu},D_{\nu}]]+[D_{\nu},[D_{\kappa},D_{\mu}]]=0$

since $[D_{\mu},F^a_{\nu\kappa}]=D_{\mu}F^a_{\nu\kappa}$. Define the dual strength tensor $\tilde{F}^{\mu\nu}=\frac{1}{2}\varepsilon^{\mu\nu\rho\sigma}F_{\rho\sigma}$, then the Bianchi identity can be rewritten as

$D_{\mu}\tilde{F}^{\mu\nu}=0.$

A source $J_\mu^a$ enters into the equations of motion as

$\partial^\mu F_{\mu\nu}^a+gf^{abc}A^{b\mu}F_{\mu\nu}^c=-J_\nu^a.$

Note that the currents must properly change under gauge group transformations.

We give here some comments about the physical dimensions of the coupling. We note that, in D dimensions, the field scales as $[A]=[L^\frac{2-D}{2}]$ and so the coupling must scale as $[g^2]=[L^{D-4}]$. This implies that Yang–Mills theory is not renormalizable for dimensions greater than four. Further, we note that, for D = 4, the coupling is dimensionless and both the field and the square of the coupling have the same dimensions of the field and the coupling of a massless quartic scalar field theory. So, these theories share the scale invariance at the classical level.

Quantization of Yang–Mills theory

A method of quantizing the Yang–Mills theory is by functional methods, i.e. path integrals. One introduces a generating functional for n-point functions as

$Z[j]=\int [dA]\exp\left[- \frac{i}{2} \int d^4x\operatorname{Tr}(F^{\mu \nu} F_{\mu \nu})+i\int 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 [dA] [d\bar c] [dc] \exp\left\{iS_F[\partial A,A]+iS_{gf}[\partial A]+iS_g[\partial c,\partial\bar c,c,\bar c,A]\right\} \\ &\exp\left\{i\int d^4x j^a_\mu(x)A^{a\mu}(x)+i\int d^4x[\bar c^a(x)\varepsilon^a(x)+\bar\varepsilon^a(x) c^a(x)]\right\} \end{align}$

being

$S_F=- \frac{1}{2}\operatorname{Tr}(F^{\mu \nu} F_{\mu \nu})$

for the field,

$S_{gf}=-\frac{1}{2\xi}(\partial\cdot A)^2$

for the gauge fixing and

$S_g=-(\bar c^a\partial_\mu\partial^\mu c^a+g\bar c^a f^{abc}\partial_\mu A^{b\mu}c^c)$

for the ghost. This is the expression commonly used to derive Feynman's rules (see Feynman diagram). Here we have c^a 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 diagrams can be obtained when the generating functional given above is rewritten as

$\begin{align} Z[j,\bar\varepsilon,\varepsilon] &= \exp\left(-ig\int 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(-ig\int d^4xf^{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 d^4xf^{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 d^4xd^4y\bar\varepsilon^a(x)C^{ab}(x-y)\varepsilon^b(y)\right)\exp\left(\tfrac{1}{2}\int d^4xd^4yj^a_\mu(x)D^{ab\mu\nu}(x-y)j^b_\nu(y)\right)$

being the generating functional of the free theory. Expanding in g and computing the functional derivatives, we are able to obtain all the n-point functions with perturbation theory. Using LSZ reduction formula we get from the n-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 g 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.

Propagators

In order to understand the behavior of the theory at large and small momenta, a key quantity is the propagator. For a Yang–Mills theory we have to consider both the gluon and the ghost propagators. At large momenta (ultraviolet limit), the question was completely settled with the discovery of the asymptotic freedom. In this case it is seen that the theory becomes free (trivial ultraviolet fixed point for renormalization group) and both the gluon and ghost propagators are those of a free massless particle. The asymptotic states of the theory are represented by massless gluons that carry the interaction. The coupling runs to zero as we will see in the next section.

At low momenta (infrared limit) the question has been more involved to settle. The reason is that the theory becomes strongly coupled in this case and perturbation theory cannot be applied. The only reliable approach to get an answer is performing lattice computation on a computer powerful enough to afford large volumes. An answer to this question is a fundamental one as it would provide an understanding to the problem of confinement. On the other side, it should not be forgotten that propagators are gauge-dependent quantities and so, they must be managed carefully when one wants to get meaningful physical results.

On the other side, theoretical approaches were conceived to get an understanding of the theory in this case. Pioneering works were due to Vladimir Gribov and Daniel Zwanziger. Gribov uncovered the question of the gauge-fixing in a Yang–Mills theory: He showed that, even once a gauge is fixed, a freedom is left yet (Gribov ambiguity). Besides, he was able to provide a functional form for the gluon propagator in the Landau gauge

$D_{\mu\nu}^{ab}(p)=\delta^{ab}\left(\eta_{\mu\nu}-\frac{p_\mu p_\nu}{p^2}\right)\frac{p^2}{p^4+M^4}.$

This propagator cannot be correct in this way as it would violate causality. On the other side, it provides a linear rising potential, $V(r)\propto r$, that would give reason to quark confinement. An important aspect of this functional form is that the gluon propagator appears to go to zero with momenta. This will become a crucial point in the following. From these studies by Gribov, Zwanziger extended his approach. The inescapable conclusion was that the gluon propagator should go to zero with momenta while the ghost propagator should be enhanced with respect to the free case running to infinity. This became known in literature as the Gribov-Zwanziger scenario. When this scenario was proposed, computational resources were insufficient to decide whether it was correct or not. Rather, people pursued a different approach using the Dyson-Schwinger equations. This is a set of coupled equations for the n-point functions of the theory forming a hierarchy. This means that the equation for the n-point function will depend on the (n+1)-point function. So, to solve them one needs a proper truncation. On the other side, these equations are non-perturbative and could permit to obtain the behavior of the n-point functions in any regime. A solution to this hierarchy through truncation was proposed by Reinhard Alkofer, Andreas Hauck and Lorenz von Smekal. This paper and the following publications from this group, the German group, set the agenda for the determination of the behavior of the propagators in the Landau gauge in the subsequent years. The main conclusions these authors arrived to were to confirm the Gribov-Zwanziger scenario and that the running coupling should reach a finite non-null fixed point when momenta runs to zero. This paper represents the birth of the so-called scaling solution as the propagators are seen to obey scaling laws with given exponents. A proposal in the eighties by John Cornwall was in contrast with this scenario rather showing that the gluons get massive when momenta goes to zero and the propagator should be finite and non-null there but went ignored at that time because the theoretical evidence appeared overwhelming for the Gribov-Zwanziger scenario. Attempts to solve the Dyson-Schwinger equations numerically seemed to provide a different scenario but this could have been due to the way truncation and approximations were applied.

The significant improvement in the computational resources made possible to unveil the proper behavior of the propagators in the Landau gauge. These results where firstly announced in Regensburg at the Lattice 2007 Conference. The results were somewhat unexpected and an example is given in the following figure for the gluon propagator

that was obtained for the SU(2) case with a lattice of $128^4\sim10^8$ points reaching momenta in the very deep infrared. This result from a huge lattice shows that the gluon propagator never goes to zero with momenta but rather reaches a plateau with a finite value at zero momenta. This went called the decoupling solution in literature. Similarly, the ghost propagator is seen to behave as that of a free particle. The ghost field just decouples from the gauge field and becomes free in the deep infrared. Other groups at the same conference confirmed similar results.

The decoupling scenario is consistent with a Yukawa-like propagator in the very deep infrared

$D_{\mu\nu}^{ab}(p)\stackrel{p\rightarrow 0}{=}\delta^{ab}\left(\eta_{\mu\nu}-\frac{p_\mu p_\nu}{p^2}\right)\frac{Z}{p^2-M^2+i0},$

with $Z$ a constant. The gluon field develops a mass gap parametrized by $M$ in the above formula, while the BRST symmetry appears to be dynamically broken. These results hold in dimensions greater than 2 while for two dimensions the scaling solution holds. Today, this scenario is generally accepted as the correct one for Yang-Mills theories in the infrared limit having such a strong support from lattice computations. Researches are ongoing for a deeper theoretical understanding of these results and eventual phenomenological applications.

Beta function and running coupling

One of the key properties of a quantum field theory is the behavior over all the energy range of the running coupling. Such a behavior can be obtained from a theory once its beta function is known. Our ability to extract results from a quantum field theory relies on perturbation theory. Once the beta function is known, the behavior at all energy scales of the running coupling is obtained through the equation

$\mu^2\frac{d\alpha_s}{d\mu^2}=\beta(\alpha_s).$

being $\alpha_s=g^2/4\pi$. Yang–Mills theory has the property of being asymptotically free in the large energy limit (ultraviolet limit). This means that, in this limit, the beta function has a minus sign driving the behavior of the running coupling toward even smaller values as the energy increases. Perturbation theory permits to evaluate beta function in this limit producing the following result for SU(N)

$\beta(\alpha_s)=-\frac{11N}{12\pi}\alpha_s^2-\frac{17N^2}{24\pi^2}\alpha_s^3+O\left(\alpha_s^4\right).$

In the opposite limit of low energies (infrared limit), the beta function is not known. It is note the exact one for a supersymmetric Yang–Mills theory. This has been obtained by Novikov, Shifman, Vainshtein and Zakharov and can be written as

$\beta(\alpha_s)=-\frac{\alpha_s^2}{4\pi}\frac{3N}{1-\frac{N\alpha_s}{2\pi}}.$

With this starting point, Thomas Ryttov and Francesco Sannino were able to postulate a non-supersymmetric version of it writing down

$\beta(\alpha_s)=-\alpha_s^2\frac{11N}{12\pi}\frac{1}{1-\frac{17N}{11}\frac{\alpha_s}{2\pi}}.$

As can be seen from the beta function of the supersymmetric theory, the limit of a large coupling (infrared limit) implies

$\beta(\alpha_s)\approx\frac{3}{2}\alpha_s.$

and so the running coupling in the deep infrared limit goes to zero making this theory trivial. This implies that the coupling reaches a maximum at some value of the energy turning again to zero as the energy is lowered. Then, if Ryttov and Sannino hypothesis is correct, the same should be true for ordinary Yang–Mills theory. This would be in agreement with recent lattice computations.

Open problems

Yang–Mills theories met with general acceptance in the physics community after Gerard 't Hooft, in 1972, worked out their renormalization, relying on a formulation of the problem worked out by his advisor Martinus Veltman. (Their work was recognized by the 1999 Nobel prize in physics.) Renormalizability is obtained even if the gauge bosons described by this theory are massive, as in the electroweak theory, provided the mass is only an "acquired" one, generated by the Higgs mechanism.

Concerning the mathematics, it should be noted that presently, i.e. in 2014, the Yang–Mills theory is a very active field of research, yielding e.g. invariants of differentiable structures on four-dimensional manifolds via work of Simon Donaldson. Furthermore, the field of Yang–Mills theories was included in the Clay Mathematics Institute's list of "Millennium Prize Problems". Here the prize-problem consists, especially, in a proof of the conjecture that the lowest excitations of a pure Yang–Mills theory (i.e. without matter fields) have a finite mass-gap with regard to the vacuum state. Another open problem, connected with this conjecture, is a proof of the confinement property in the presence of additional Fermion particles.

In physics the survey of Yang–Mills theories does not usually start from perturbation analysis or analytical methods, but more recently from systematic application of numerical methods to lattice gauge theories.

See also

• Aharonov–Bohm effect
• Coulomb gauge
• Electroweak theory
• Gauge covariant derivative
• Kaluza–Klein theory
• Lattice gauge theory
• Lorenz gauge
• N = 4 supersymmetric Yang–Mills theory
• Propagator
• Quantum chromodynamics
• Quantum gauge theory
• Field theoretical formulation of the standard model
• Symmetry in physics
• Weyl gauge
• Yang–Mills existence and mass gap
• Yang–Mills–Higgs equations

Further reading

Books

Articles

External links

• Yang–Mills theory on DispersiveWiki
• The Clay Mathematics Institute
• The Millennium Prize Problems