Polyakov loop

Order parameter

Gauge fields need to satisfy the periodicity condition ${\displaystyle A_{\mu }(\mathit{x},x_4+\beta )=A_{\mu }(\mathit{x},x_4)}$ in the compactified direction. Meanwhile, gauge transformations only need to satisfy this up to a group center term ${\displaystyle h}$ as ${\displaystyle \mathrm{\Omega }(\mathit{x},x_4+\beta )=h\mathrm{\Omega }(\mathit{x},x_4)}$. A change of basis can always diagonalize this so that ${\displaystyle h=zI}$ for a complex number ${\displaystyle z}$. The Polyakov loop is topologically nontrivial in the temporal direction so unlike other Wilson loops it transforms as ${\displaystyle \mathrm{\Phi }(\mathit{x})\to z\mathrm{\Phi }(\mathit{x})}$ under these transformations. Since this makes the loop gauge dependent for ${\displaystyle z\ne 1}$, by Elitzur's theorem non-zero expectation values of ${\displaystyle ⟨\mathrm{\Phi }⟩}$ imply that the center group must be spontaneously broken, implying confinement in pure gauge theory. This makes the Polyakov loop an order parameter for confinement in thermal pure gauge theory, with a confining phase occurring when ${\displaystyle ⟨\mathrm{\Phi }⟩=0}$ and deconfining phase when ${\displaystyle ⟨\mathrm{\Phi }⟩\ne 0}$. For example, lattice calculations of quantum chromodynamics with infinitely heavy quarks that decouple from the theory shows that the deconfinement phase transition occurs at around a temperature of ${\displaystyle 270}$ MeV. Meanwhile, in a gauge theory with quarks, these break the center group and so confinement must instead be deduced from the spectrum of asymptotic states, the color neutral hadrons.

For gauge theories that lack a nontrivial group center that could be broken in the confining phase, the Polyakov loop expectation values are nonzero even in this phase. They are however still a good indicator of confinement since they generally experience a sharp jump at the phase transition. This is the case for example in the Higgs model with the exceptional gauge group ${\displaystyle G_2}$.

The Nambu–Jona-Lasinio model lacks local color symmetry and thus cannot capture the effects of confinement. However, Polyakov loops can be used to construct the Polyakov-loop-extended Nambu–Jona-Lasinio model which treats both the chiral condensate and the Polyakov loops as classical homogeneous fields that couple to quarks according to the symmetries and symmetry breaking patters of quantum chromodynamics.

Quark free energy

The free energy ${\displaystyle F}$ of ${\displaystyle N}$ quarks and ${\displaystyle \overline{N}}$ antiquarks, subtracting out the vacuum energy, is given in terms of the correlation functions of Polyakov loops

${\displaystyle e^{-\beta F}=⟨\mathrm{\Phi }({\mathit{x}}_1)\dots \mathrm{\Phi }({\mathit{x}}_{N_q}){\mathrm{\Phi }}^{†}({\mathit{x}}_1^{\prime })\cdots {\mathrm{\Phi }}^{†}({\mathit{x}}_{\overline{N}}^{\prime })⟩.}$

This free energy is another way to see that the Polyakov loop acts as an order parameter for confinement since the free energy of a single quark is given by ${\displaystyle e^{-\beta \mathrm{\Delta }F}=⟨\mathrm{\Phi }(\mathit{x})⟩}$. Confinement of quarks means that it would take an infinite amount of energy to create a configuration with a single free quark, therefore its free energy must be infinite and so the Polyakov loop expectation value must vanish in this phase, in agreement with the center symmetry breaking argument.

The formula for the free energy can also be used to calculate the potential between a pair of infinitely massive quarks spatially separated by ${\displaystyle r=|{\mathit{x}}_1-{\mathit{x}}_2|}$. Here the potential ${\displaystyle V(r)}$ is the first term in the free energy, so that the correlation function of two Polyakov loops is

${\displaystyle ⟨\mathrm{\Phi }({\mathit{x}}_1)\mathrm{\Phi }({\mathit{x}}_2)⟩\propto e^{-\beta V(r)}(1+\mathcal{O}(e^{-\beta \mathrm{\Delta }E(r)})),}$

where ${\displaystyle \mathrm{\Delta }E}$ is the energy difference between the potential and the first excited state. In the confining phase the potential is linear ${\displaystyle V(r)=\sigma r}$, where the constant of proportionality is known as the string tension. The string tension acquired from the Polyakov loop is always bounded from above by the string tension acquired from the Wilson loop.