We present results on the chiral and deconfinement properties of the QCD transition at finite temperature. Calculations are performed with 2 + 1 flavors of quarks using the p4, asqtad and HISQ/tree actions. Lattices with temporal extent Nτ = 6, 8 and 12 are used to understand and control discretization errors and to reliably extrapolate estimates obtained at finite lattice spacings to the continuum limit. The chiral transition temperature is defined in terms of the phase transition in a theory with two massless flavors and analyzed using O(N ) scaling fits to the chiral condensate and susceptibility. We find consistent estimates from the HISQ/tree and asqtad actions and our main result is Tc = 154 ± 9 MeV.

We present results for the equation of state in (2+1)-flavor QCD using the highly improved staggered quark action and lattices with temporal extent Nτ = 6, 8, 10, and 12. We show that these data can be reliably extrapolated to the continuum limit and obtain a number of thermodynamic quantities and the speed of sound in the temperature range (130-400) MeV. We compare our results with previous calculations, and provide an analytic parameterization of the pressure, from which other thermodynamic quantities can be calculated, for use in phenomenology. We show that the energy density in the crossover region, 145 MeV ≤ T ≤ 163 MeV, defined by the chiral transition, is c = (0.18 − 0.5) GeV/fm 3 , i.e., (1.2 − 3.1) nuclear . At high temperatures, we compare our results with resummed and dimensionally reduced perturbation theory calculations. As a byproduct of our analyses, we obtain the values of the scale parameters r0 from the static quark potential and w0 from the gradient flow.

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 ∈ [135 MeV, 330 MeV] using up to four different sets of lattice cut-offs corresponding to lattices of size N 3 σ × Nτ with aspect ratio Nσ/Nτ = 4 and Nτ = 6 − 16. The strange quark mass is tuned to its physical value and we use two strange to light quark mass ratios ms/m l = 20 and 27, which in the continuum limit correspond to a pion mass of about 160 MeV 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 (µ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 √ sNN ∼ 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 -µB plane and are compared with the crossover line for the QCD chiral transition as well as with experimental 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(µB = 0) > 0.9.

We present results for pseudo-critical temperatures of QCD chiral crossovers at zero and non-zero values of baryon (B), strangeness (S), electric charge (Q), and isospin (I) chemical potentials µ X=B,Q,S,I . The results were obtained using lattice QCD calculations carried out with two degenerate up and down dynamical quarks and a dynamical strange quark, with quark masses corresponding to physical values of pion and kaon masses in the continuum limit. By parameterizing pseudo-critical temperatures as (0)) 4 , we determined κ X 2 and κ X 4 from Taylor expansions of chiral observables in µ X . We obtained a precise result for T c (0) = (156.5 ± 1.5) MeV. For analogous thermal conditions at the chemical freeze-out of relativistic heavy-ion collisions, i.e., µ S (T, µ B ) and µ Q (T, µ B ) fixed from strangeness-neutrality and isospin-imbalance, we found κ B 2 =0.012(4) and κ B 4 =0.000(4). For µ B 300 MeV, the chemical freeze-out takes place in the vicinity of the QCD phase boundary, which coincides with the lines of constant energy density of 0.42(6) GeV/fm 3 and constant entropy density of 3.7(5) fm −3 .

We calculate the vector current correlation function for light valence quarks in the deconfined phase of QCD. The calculations have been performed in quenched lattice QCD at T ≃ 1.45T c for four values of the lattice cut-off on lattices up to size 128 3 × 48. This allows to perform a continuum extrapolation of the correlation function in the Euclidean time interval 0.2 ≤ τ T ≤ 0.5, which extends to the largest temporal separations possible at finite temperature, to better than 1% accuracy. In this interval, at the value of the temperature investigated, we find that the vector correlation function never deviates from the free correlator for massless quarks by more than 9%. We also determine the first two non-vanishing thermal moments of the vector meson spectral function. The second thermal moment deviates by less than 7% from the free value. With these constraints, we then proceed to extract information on the spectral representation of the vector correlator and discuss resulting consequences for the electrical conductivity and the thermal dilepton rate in the plasma phase.