No Arabic abstract
We propose a novel Bayesian method to analytically continue observables to real baryochemical potential $mu_B$ in finite density QCD. Taylor coefficients at $mu_B=0$ and data at imaginary chemical potential $mu_B^I$ are treated on equal footing. We consider two different constructions for the Pade approximants, the classical multipoint Pade approximation and a mixed approximation that is a slight generalization of a recent idea in Pade approximation theory. Approximants with spurious poles are excluded from the analysis. As an application, we perform a joint analysis of the available continuum extrapolated lattice data for both pseudocritical temperature $T_c$ at $mu_B^I$ from the Wuppertal-Budapest Collaboration and Taylor coefficients $kappa_2$ and $kappa_4$ from the HotQCD Collaboration. An apparent convergence of $[p/p]$ and $[p/p+1]$ sequences of rational functions is observed with increasing $p.$ We present our extrapolation up to $mu_Bapprox 600$ MeV.
We study the phase diagram of QCD at finite isospin density using two flavors of staggered quarks. We investigate the low temperature region of the phase diagram where we find a pion condensation phase at high chemical potential. We started a basic analysis of the spectrum at finite isospin density. In particular, we measured pion, rho and nucleon masses inside and outside of the pion condensation phase. In agreement with previous studies in two-color QCD at finite baryon density we find that the Polyakov loop does not depend on the density in the staggered formulation.
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.
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.
We study different estimators of the radius of convergence of the Taylor series of the pressure in finite density QCD. We adopt the approach in which the radius of convergence is estimated first in a finite volume, and the infinite-volume limit is taken later. This requires an estimator for the radius of convergence that is reliable in a finite volume. Based on general arguments about the analytic structure of the partition function in a finite volume, we demonstrate that the ratio estimator cannot work in this approach, and propose three new estimators, capable of extracting reliably the radius of convergence, which coincides with the distance from the origin of the closest Lee-Yang zero. We also provide an estimator for the phase of the closest Lee-Yang zero, necessary to assess whether the leading singularity is a true critical point. We demonstrate the usage of these estimators on a toy model, namely 4 flavors of unimproved staggered fermions on a small $6^3 times 4$ lattice, where both the radius of convergence and the Taylor coefficients to any order can be obtained by a direct determination of the Lee-Yang zeros. Interestingly, while the relative statistical error of the Taylor expansion coefficients steadily grows with order, that of our estimators stabilizes, allowing for an accurate estimate of the radius of convergence. In particular, we show that despite the more than 100% error bars on high-order Taylor coefficients, the given ensemble contains enough information about the radius of convergence.