QCD equation of state to O(μB6) from lattice QCD

We calculated the QCD equation of state using Taylor expansions that include contributions from up to sixth order in the baryon strangeness and electric charge chemical potentials. Calculations have been performed with the Highly Improved Staggered Quark action in the temperature range T epsilon [135 MeV 330 MeV] using up to four different sets of lattice cutoffs corresponding to lattices of size N sigma 3x N tau with aspect ratio N sigma/N tau = 4 and N tau=6-16. The strange quark mass is tuned to its physical value and we use two strange to light quark mass ratios ms/ml = 20 and 27 which in the continuum limit correspond to a pion mass of about 160 and 140 MeV respectively. Sixth order results for Taylor expansion coefficients are used to estimate truncation errors of the fourth order expansion. We show that truncation errors are small for baryon chemical potentials less then twice the temperature (mu(B) <= 2T). The fourth order equation of state thus is suitable for the modeling of dense matter created in heavy ion collisions with center of mass energies down to root SNN similar to 12 GeV. We provide a parametrization of basic thermodynamic quantities that can be readily used in hydrodynamic simulation codes. The results on up to sixth order expansion coefficients of bulk thermodynamics are used for the calculation of lines of constant pressure energy and entropy densities in the T mu(B) plane and are compared with the crossover line for the QCD chiral transition as well as ith xperimental results on freeze out parameters in heavy ion collisions. These coefficients also provide estimates for the location of a possible critical point. We argue that results on sixth order expansion coefficients disfavor the existence of a critical point in the QCD phase diagram for B= T <= 2 and T/Tc(mu(B) = 0) > 0.9.