No Arabic abstract
We calculate the equation of state in 2+1 flavor QCD at finite temperature with physical strange quark mass and almost physical light quark masses using lattices with temporal extent Nt=8. Calculations have been performed with two different improved staggered fermion actions, the asqtad and p4 actions. Overall, we find good agreement between results obtained with these two O(a^2) improved staggered fermion discretization schemes. A comparison with earlier calculations on coarser lattices is performed to quantify systematic errors in current studies of the equation of state. We also present results for observables that are sensitive to deconfining and chiral aspects of the QCD transition on Nt=6 and 8 lattices. We find that deconfinement and chiral symmetry restoration happen in the same narrow temperature interval. In an Appendix we present a simple parametrization of the equation of state that can easily be used in hydrodynamic model calculations. In this parametrization we also incorporated an estimate of current uncertainties in the lattice calculations which arise from cutoff and quark mass effects. We estimate these systematic effects to be about 10 MeV
We determine the equation of state of QCD at finite chemical potential, to order $(mu_B/T)^6$, for a system of 2+1 quark flavors. The simulations are performed at the physical mass for the light and strange quarks on several lattice spacings; the results are continuum extrapolated using lattices of up to $N_t=16$ temporal resolution. The QCD pressure and interaction measure are calculated along the isentropic trajectories in the $(T,~mu_B)$ plane corresponding to the RHIC Beam Energy Scan collision energies. Their behavior is determined through analytic continuation from imaginary chemical potentials of the baryonic density. We also determine the Taylor expansion coefficients around $mu_B=0$ from the simulations at imaginary chemical potentials. Strangeness neutrality and charge conservation are imposed, to match the experimental conditions.
The QCD equation of state at finite baryon density is studied in the framework of a Cluster Expansion Model (CEM), which is based on the fugacity expansion of the net baryon density. The CEM uses the two leading Fourier coefficients, obtained from lattice simulations at imaginary $mu_B$, as the only model input and permits a closed analytic form. Excellent description of the available lattice data at both $mu_B = 0$ and at imaginary $mu_B$ is obtained. We also demonstrate how the Fourier coefficients can be reconstructed from baryon number susceptibilities.
We study the equation of state at finite temperature and density in two-flavor QCD with the RG-improved gluon action and the clover-improved Wilson quark action on a $ 16^3 times 4$ lattice. Along the lines of constant physics at $m_{rm PS}/m_{rm V} = 0.65$ and 0.80, we compute the second and forth derivatives of the grand canonical partition function with respect to the quark chemical potential $mu_q = (mu_u+mu_d)/2$ and the isospin chemical potential $mu_I = (mu_u-mu_d)/2$ at vanishing chemical potentials, and study the behaviors of thermodynamic quantities at finite $mu_q$ using these derivatives for the case $mu_I=0$. In particular, we study density fluctuations at none-zero temperature and density by calculating the quark number and isospin susceptibilities and their derivatives with respect to $mu_q$. To suppress statistical fluctuations, we also examine new techniques applicable at low densities. We find a large enhancement in the fluctuation of quark number when the density increased near the pseudo-critical temperature, suggesting a critical point at finite $mu_q$ terminating the first order transition line between hadronic and quark gluon plasma phases. This result agrees with the previous results using staggered-type quark actions qualitatively. Furthermore, we study heavy-quark free energies and Debye screening masses at finite density by measuring the first and second derivatives of these quantities for various color channels of heavy quark-quark and quark-anti-quark pairs. The results suggest that, to the leading order of $mu_q$, the interaction between two quarks becomes stronger at finite densities, while that between quark and anti-quark becomes weaker.
We study the phase structure of lattice QCD with heavy quarks at finite temperature and density by a histogram method. We determine the location of the critical point at which the first-order deconfining transition in the heavy-quark limit turns into a crossover at intermediate quark masses through a change of the shape of the histogram under variation of coupling parameters. We estimate the effect of the complex phase factor which causes the sign problem at finite density, and show that, in heavy-quark QCD, the effect is small around the critical point. We determine the critical surface in 2+1 flavor QCD in the heavy-quark region at all values of the chemical potential mu including mu=infty.
Neither the chiral limit nor finite baryon density can be simulated directly in lattice QCD, which severely limits our understanding of the QCD phase diagram. In this review I collect results for the phase structure in an extended parameter space of QCD, with varying numbers of flavours, quark masses, colours, lattice spacings, imaginary and isospin chemical potentials. Such studies help in understanding the underlying symmetries and degrees of freedom, and are beginning to provide a consistent picture constraining the possibilities for the physical phase diagram.