No Arabic abstract
We provide the most accurate results for the QCD transition line so far. We optimize the definition of the crossover temperature $T_c$, allowing for its very precise determination, and extrapolate from imaginary chemical potential up to real $mu_B approx 300$ MeV. The definition of $T_c$ adopted in this work is based on the observation that the chiral susceptibility as a function of the condensate is an almost universal curve at zero and imaganiary $mu_B$. We obtain the parameters $kappa_2=0.0153(18)$ and $kappa_4=0.00032(67)$ as a continuum extrapolation based on $N_t=10,12$ and $16$ lattices with physical quark masses. We also extrapolate the peak value of the chiral susceptibility and the width of the chiral transition along the crossover line. In fact, both of these are consistent with a constant function of $mu_B$. We see no sign of criticality in the explored range.
We study the density of states method as well as reweighting to explore the low temperature phase diagram of QCD at finite baryon chemical potential. We use four flavors of staggered quarks, a tree-level Symanzik improved gauge action and four stout smearing steps on lattices with $N_s=4,6,8$ and $N_t=6 - 16$. We compare our results to that of the phase quenched ensemble and also determine the pion and nucleon masses. In the density of states approach we applied pion condensate or gauge action density fixing. We found that the density of states method performs similarly to reweighting. At $T approx 100$ MeV, we found an indication of the onset of the quark number density at around $mu/m_N sim 0.16 - 0.18$ on $6^4$ lattices at $beta=2.9$.
Within the reweighting approach, one has the freedom to choose the Monte Carlo action so that it provides a good overlap with the finite-mu measure but remains simple to simulate. We explore several choices of action in the regime of small mu. Simulating with a finite isospin chemical potential mu_I=mu gives a better overlap than the standard choice mu=0, with no computational overhead.
We investigate the properties of QCD at finite isospin chemical potential at zero and non-zero temperatures. This theory is not affected by the sign problem and can be simulated using Monte-Carlo techniques. With increasing isospin chemical potential and temperatures below the deconfinement transition the system changes into a phase where charged pions condense, accompanied by an accumulation of low modes of the Dirac operator. The simulations are enabled by the introduction of a pionic source into the action, acting as an infrared regulator for the theory, and physical results are obtained by removing the regulator via an extrapolation. We present an update of our study concerning the associated phase diagram using 2+1 flavours of staggered fermions with physical quark masses and the comparison to Taylor expansion. We also present first results for our determination of the equation of state at finite isospin chemical potential and give an example for a cosmological application. The results can also be used to gain information about QCD at small baryon chemical potentials using reweighting with respect to the pionic source parameter and the chemical potential and we present first steps in this direction.
Taylor expansion of the equation of state of QCD suffers from shortcomings at chemical potentials $mu_B geq (2-2.5)T$. First, one faces difficulties inherent in performing such an expansion with a limited number of coefficients; second, higher order coefficients determined from lattice calculations suffer from a poor signal-to-noise ratio. In this work, we present a novel scheme for extrapolating the equation of state of QCD to finite, real chemical potential that can extend its reach further than previous methods. We present continuum extrapolated lattice results for the new expansion coefficients and show the thermodynamic observables up to $mu_B/Tle3.5$.
We compare higher moments of baryon numbers measured at the RHIC heavy ion collision experiments with those by the lattice QCD calculations. We employ the canonical approach, in which we can access the real chemical potential regions avoiding the sign problem. In the lattice QCD simulations, we study several fits of the number density in the pure imaginary chemical potential, and analyze how these fits affects behaviors at the real chemical potential. In the energy regions between $sqrt{s}_{NN}$=19.6 and 200 GeV, the susceptibility calculated at $T/T_c=0.93$ is consistent with experimental data at $0 le mu_B/T < 1.5$, while the kurtosis shows similar behavior with that of the experimental data in the small $mu_B/T$ regions $0 le mu_B/T < 0.3$. The experimental data at $sqrt{s}_{NN}=$ 11.5 shows quite different behavior. The lattice result in the deconfinement region,$T/T_c=1.35$, is far from experimental data.