Non-trivial Center Dominance in High Temperature QCD

We investigate the properties of quarks and gluons above the chiral phase transition temperature $T_c,$ using the RG improved gauge action and the Wilson quark action with two degenerate quarks mainly on a $32^3\times 16$ lattice. In the one-loop perturbation theory, the thermal ensemble is dominated by the gauge configurations with effectively $Z(3)$ center twisted boundary conditions, making the thermal expectation value of the spatial Polyakov loop take a non-trivial $Z(3)$ center. This is in agreement with our lattice simulation of high temperature QCD. We further observe that the temporal propagator of massless quarks at extremely high temperature $\beta=100.0 \, (T \simeq10^{58} T_c)$ remarkably agrees with the temporal propagator of free quarks with the $Z(3)$ twisted boundary condition for $t/L_t \geq 0.2$, but differs from that with the $Z(3)$ trivial boundary condition. As we increase the mass of quarks $m_q$, we find that the thermal ensemble continues to be dominated by the $Z(3)$ twisted gauge field configurations as long as $m_q \le 3.0 \, T$ and above that the $Z(3)$ trivial configurations come in. The transition is essentially identical to what we found in the departure from the conformal region in the zero-temperature many-flavor conformal QCD on a finite lattice by increasing the mass of quarks. We argue that the behavior is consistent with the renormalization group analysis at finite temperature.