Confinement and Z₃ symmetry in three-flavor QCD

We investigate the confinement mechanism in three-flavor QCD with imaginary isospin chemical potentials (μᵤ,μd,μₛ)=(iθ T,-iθ T,0), using the Polyakov-loop extended Nambu--Jona-Lasinio (PNJL) model, where T is temperature. As for three degenerate flavors, the system has Z₃ symmetry at θ=2π/3 and hence the Polyakov loop Φ vanishes there for small T. As for 2+1 flavors, the symmetry is not preserved for any θ, but Φ becomes zero at θ=θconf < 2π/3 for small T. The confinement phase defined by Φ=0 is realized, even if the system does not have Z₃ symmetry exactly. In the θ-T plane, there is a critical endpoint of deconfinement transition. The deconfinement crossover at zero chemical potential is a remnant of the first-order deconfinement transition at θ=θconf. The relation between the non-diagonal element χᵤₛ of quark number susceptibilities and the deconfinement transition is studied. The present results can be checked by lattice QCD simulations directly, since the simulations are free from the sign problem for any θ.