We investigate three-colour QCD thermodynamics at finite quark chemical potential. Lattice QCD results are compared with a generalized Nambu JonaLasinio model in which quarks couple simultaneously to the chiral condensate and to a background temporal gauge field representing Polyakov loop dynamics. This so-called PNJL model thus includes features of both deconfinement and chiral symmetry restoration. The parameters of the Polyakov loop effective potential are fixed in the pure gauge sector. The chiral condensate and the Polyakov loop as functions of temperature and quark chemical potential are calculated by minimizing the thermodynamic potential of the system. The resulting equation of state, (scaled) pressure difference and quark number density at finite quark chemical potential are then confronted with corresponding Lattice QCD data.

The present paper concludes our investigation on the QCD equation of state with 2+1 staggered flavors and one-link stout improvement. We extend our previous study [JHEP 0601:089 (2006)] by choosing even finer lattices. Lattices with $N_t=6,8$ and 10 are used, and the continuum limit is approached by checking the results at $N_t=12$. A Symanzik improved gauge and a stout-link improved staggered fermion action is utilized. We use physical quark masses, that is, for the lightest staggered pions and kaons we fix the $m_\pi/f_K$ and $m_K/f_K$ ratios to their experimental values. The pressure, the interaction measure, the energy and entropy density and the speed of sound are presented as functions of the temperature in the range $100 ...1000 \textmd{MeV}$. We give estimates for the pion mass dependence and for the contribution of the charm quark. We compare our data to the equation of state obtained by the "hotQCD" collaboration.Comment: Minor changes: Figure 1 added; Figure 15, Figure 17 and Table 5 changed. Accepted for publication in JHE

The present paper concludes our investigations on the QCD cross-over transition temperatures with 2+1 staggered flavours and one-link stout improvement. We extend our previous two studies [Phys. Lett. B643 (2006) 46, JHEP 0906:088 (2009)] by choosing even finer lattices (N t =16) and we work again with physical quark masses. The new results on this broad cross-over are in complete agreement with our earlier ones. We compare our findings with the published results of the hotQCD collaboration. All these results are confronted with the predictions of the Hadron Resonance Gas model and Chiral Perturbation Theory for temperatures below the transition region. Our results can be reproduced by using the physical spectrum in these analytic calculations. The findings of the hotQCD collaboration can be recovered by using a distorted spectrum which takes into account lattice discretization artifacts and heavier than physical quark masses. This analysis provides a simple explanation for the observed discrepancy in the transition temperatures between our and the hotQCD collaborations.

No abstract