Polyakov loop observables

The Polyakov loop \(P\) gives the order parameter of the deconfinement transition in the quenched limit. The imaginary part is the order parameter of the Roberge-Weiss transition. It is also occasionally of interest at finite quark mass and real \(\mu\). You can find a bunch of Polyakov loop related observables in the polyakovTools module.

Initialize a polyakovTools object with

pt = polyakovTools(Ns, Nt)

with spatial extent Ns and Euclidean time extent Nt. One can then get \(\langle |P|\rangle\) from

numpy.mean( pt.absP(RePmeas, ImPmeas) )

where RePmeas and ImPmeas are numpy arrays of measurements of real and imaginary parts of \(P\), respectively.