Polyakov loop

In quantum field theory, the Polyakov loop is the thermal analogue of the Wilson loop, acting as an order parameter for confinement in pure gauge theories at nonzero temperatures. In particular, it is a Wilson loop that winds around the compactified Euclidean temporal direction of a thermal quantum field theory. It indicates confinement because its vacuum expectation value must vanish in the confined phase due to its non-invariance under center gauge transformations. This also follows from the fact that the expectation value is related to the free energy of individual quarks, which diverges in this phase. Introduced by Alexander M. Polyakov in 1975,^[1] they can also be used to study the potential between pairs of quarks at nonzero temperatures.

Thermal quantum field theory is formulated in Euclidean spacetime with a compactified imaginary temporal direction of length β {\displaystyle \beta } . This length corresponds to the inverse temperature of the field β ∝ 1 / T {\displaystyle \beta \propto 1/T} . Compactification leads to a special class of topologically nontrivial Wilson loops that wind around the compact direction known as Polyakov loops.^[2] In SU ( N ) {\displaystyle {\text{SU}}(N)} theories a straight Polyakov loop on a spatial coordinate x {\displaystyle {\boldsymbol {x}}} is given by

where P {\displaystyle {\mathcal {P}}} is the path-ordering operator and A 4 {\displaystyle A_{4}} is the Euclidean temporal component of the gauge field. In lattice field theory this operator is reformulated in terms of temporal link fields U 4 ( m , j ) {\displaystyle U_{4}({\boldsymbol {m}},j)} at a spatial position m {\displaystyle {\boldsymbol {m}}} as^[3]

Φ ( m ) = 1 N tr [ ∏ j = 0 N T − 1 U 4 ( m , j ) ] . {\displaystyle \Phi ({\boldsymbol {m}})={\frac {1}{N}}{\text{tr}}{\bigg [}\prod _{j=0}^{N_{T}-1}U_{4}({\boldsymbol {m}},j){\bigg ]}.}

The continuum limit of the lattice must be taken carefully to ensure that the compact direction has fixed extent. This is done by ensuring that the finite number of temporal lattice points N T {\displaystyle N_{T}} is such that β = N T a {\displaystyle \beta =N_{T}a} is constant as the lattice spacing a {\displaystyle a} goes to zero.

Gauge fields need to satisfy the periodicity condition A μ ( x , x 4 + β ) = A μ ( x , x 4 ) {\displaystyle A_{\mu }({\boldsymbol {x}},x_{4}+\beta )=A_{\mu }({\boldsymbol {x}},x_{4})} in the compactified direction. Meanwhile, gauge transformations only need to satisfy this up to a group center term h {\displaystyle h} as Ω ( x , x 4 + β ) = h Ω ( x , x 4 ) {\displaystyle \Omega ({\boldsymbol {x}},x_{4}+\beta )=h\Omega ({\boldsymbol {x}},x_{4})} . A change of basis can always diagonalize this so that h = z I {\displaystyle h=zI} for a complex number z {\displaystyle z} . The Polyakov loop is topologically nontrivial in the temporal direction so unlike other Wilson loops it transforms as Φ ( x ) → z Φ ( x ) {\displaystyle \Phi ({\boldsymbol {x}})\rightarrow z\Phi ({\boldsymbol {x}})} under these transformations.^[5] Since this makes the loop gauge dependent for z ≠ 1 {\displaystyle z\neq 1} , by Elitzur's theorem non-zero expectation values of ⟨ Φ ⟩ {\displaystyle \langle \Phi \rangle } 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 ⟨ Φ ⟩ = 0 {\displaystyle \langle \Phi \rangle =0} and deconfining phase when ⟨ Φ ⟩ ≠ 0 {\displaystyle \langle \Phi \rangle \neq 0} .^[6] 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 270 {\displaystyle 270} MeV.^[7] 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 G 2 {\displaystyle G_{2}} .^[8]

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.^[9][10][11]

The free energy F {\displaystyle F} of N {\displaystyle N} quarks and N ¯ {\displaystyle {\bar {N}}} antiquarks, subtracting out the vacuum energy, is given in terms of the correlation functions of Polyakov loops^[12]

e − β F = ⟨ Φ ( x 1 ) … Φ ( x N q ) Φ † ( x 1 ′ ) ⋯ Φ † ( x N ¯ ′ ) ⟩ . {\displaystyle e^{-\beta F}=\langle \Phi ({\boldsymbol {x}}_{1})\dots \Phi ({\boldsymbol {x}}_{N_{q}})\Phi ^{\dagger }({\boldsymbol {x}}'_{1})\cdots \Phi ^{\dagger }({\boldsymbol {x}}_{\bar {N}}')\rangle .}

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 e − β Δ F = ⟨ Φ ( x ) ⟩ {\displaystyle e^{-\beta \Delta F}=\langle \Phi ({\boldsymbol {x}})\rangle } .^[13] 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 r = | x 1 − x 2 | {\displaystyle r=|{\boldsymbol {x}}_{1}-{\boldsymbol {x}}_{2}|} . Here the potential V ( r ) {\displaystyle V(r)} is the first term in the free energy, so that the correlation function of two Polyakov loops is

⟨ Φ ( x 1 ) Φ ( x 2 ) ⟩ ∝ e − β V ( r ) ( 1 + O ( e − β Δ E ( r ) ) ) , {\displaystyle \langle \Phi ({\boldsymbol {x}}_{1})\Phi ({\boldsymbol {x}}_{2})\rangle \propto e^{-\beta V(r)}(1+{\mathcal {O}}(e^{-\beta \Delta E(r)})),}

where Δ E {\displaystyle \Delta E} is the energy difference between the potential and the first excited state. In the confining phase the potential is linear V ( r ) = σ r {\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.^[14]

• Quark–gluon plasma
• 't Hooft loop