No Arabic abstract
We report on a continuum extrapolated result [arXiv:1309.5258] for the equation of state (EoS) of QCD with $N_f=2+1$ dynamical quark flavors. In this study, all systematics are controlled, quark masses are set to their physical values, and the continuum limit is taken using at least three lattice spacings corresponding to temporal extents up to $N_t=16$. A Symanzik improved gauge and stout-link improved staggered fermion action is used. Our results are available online [ancillary file to arXiv:1309.5258].
We present a full result for the 2+1 flavor QCD equation of state. All the systematics are controlled, the quark masses are set to their physical values, and the continuum extrapolation is carried out. This extends our previous studies [JHEP 0601:089 (2006); 1011:077 (2010)] to even finer lattices and now includes ensembles with Nt = 6,8,10,12 up to Nt = 16. We use a Symanzik improved gauge and a stout-link improved staggered fermion action. Our findings confirm our earlier results. In order to facilitate the direct use of our equation of state we make our tabulated results available for download.
We present results for the QCD equation of state, quark densities and susceptibilities at nonzero chemical potential, using 2+1 flavor asqtad ensembles with $N_t=4$. The ensembles lie on a trajectory of constant physics for which $m_{ud}approx0.1m_s$. The calculation is performed using the Taylor expansion method with terms up to sixth order in $mu/T$.
We study the critical point for finite temperature Nf=3 QCD using several temporal lattice sizes up to 10. In the study, the Iwasaki gauge action and non-perturbatively O(a) improved Wilson fermions are employed. We estimate the critical temperature and the upper bound of the critical pseudo-scalar meson mass.
The energy-momentum tensor plays an important role in QCD thermodynamics. Its expectation value contains information of the pressure and the energy density as its diagonal part. Further properties like viscosity and specific heat can be extracted from its correlation function. Recently a new method based on the gradient flow was introduced to calculate the energy-momentum tensor on the lattice, and has been successfully applied to quenched QCD. In this paper, we apply the gradient flow method to calculate the energy-momentum tensor in (2+1)-flavor QCD. As the first application of the method with dynamical quarks, we study at a single but fine lattice spacing a=0.07 fm with heavy u and d quarks ($m_pi/m_rho=0.63$) and approximately physical s quark. Performing simulations on lattices with Nt=16 to 4, the temperature range of T=174-697 MeV is covered. We find that the results of the pressure and the energy density by the gradient flow method are consistent with the previous results using the T-integration method at T<280 MeV, while the results show disagreement at T>350 MeV (Nt<8), presumably due to the small-Nt lattice artifact of $O((aT)^2)=O(1/N_t^2)$. We also apply the gradient flow method to evaluate the chiral condensate taking advantage of the gradient flow method that renormalized quantities can be directly computed avoiding the difficulty of explicit chiral violation with lattice quarks. We compute the renormalized chiral condensate in the MS-bar scheme at renormalization scale $mu=2$ GeV with a high precision to study the temperature dependence of the chiral condensate and its disconnected susceptibility. Even with the Wilson-type quark action, we obtain the chiral condensate and its disconnected susceptibility showing a clear signal of pseudocritical temperature at T~190 MeV related to the chiral restoration crossover.
We report results for the interaction measure, pressure and energy density for nonzero temperature QCD with 2+1 flavors of improved staggered quarks. In our simulations we use a Symanzik improved gauge action and the Asqtad $O(a^2)$ improved staggered quark action for lattices with temporal extent $N_t=4$ and 6. The heavy quark mass $m_s$ is fixed at approximately the physical strange quark mass and the two degenerate light quarks have masses $m_{ud}approx0.1 m_s$ or $0.2 m_s$. The calculation of the thermodynamic observables employs the integral method where energy density and pressure are obtained by integration over the interaction measure.