Lattice QCD Thermodynamics Research Articles

Quasi particle model vs lattice QCD thermodynamics: extension to Nf=2+1+1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N_f=2+1+1$$\\end{document} flavors and momentum dependent quark masses

In the last decade a quasi-particle model (QPM) has supplied the basis for the study of heavy quark (HQ) production in ultra-relativistic collisions, allowing for a phenomenological estimate of the HQ diffusion coefficient Ds(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$D_s(T)$$\\end{document}. Using the new lattice QCD results for the equation of state (EoS) with 2+1+1 dynamical flavors, we extend the QPM from Nf=2+1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N_f=2+1$$\\end{document} to Nf=2+1+1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N_f=2+1+1$$\\end{document}, where the charm quark is included. Fixing the coupling g(T) by a fit to the lQCD energy density ϵ(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\epsilon (T)$$\\end{document}, we evaluate the impact of different temperature parametrizations of charm quark mass on EoS and susceptibilities χq(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _q(T)$$\\end{document} of light, χs(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _s(T)$$\\end{document} of strange and χc(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _c(T)$$\\end{document} of charm quarks, the last favouring a charm quark mass increasing toward Tc\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$T_c$$\\end{document}. We also explore the extension of the QPM to a more realistic approach called QPMp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_p$$\\end{document}, where quark and gluon masses explicitly depend on their momentum converging to the current quark mass at high momenta, as expected from asymptotic free dynamics. The QPMp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_p$$\\end{document} allows for a simultaneous quantitative description not only of the EoS but also of the quark susceptibilities (χq(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _q(T)$$\\end{document}, χs(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _s(T)$$\\end{document}), which instead are underestimated in the simple QPM. Furthermore, evaluating the spatial diffusion coefficient 2πTDs(T)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$2\\pi T D_s(T)$$\\end{document} in the QPMp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_p$$\\end{document}, we find it is also closer than QPM to the recent lQCD data performed including dynamical fermions. Finally, in a 1+1D expanding system, we evaluate the RAA(pT)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$R_{AA}(p_T)$$\\end{document} in the QPM and QPMp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_p$$\\end{document}, finding a significant reduction at low momenta for QPMp\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$_p$$\\end{document} which could lead in a realistic scenario to a better agreement to experimental data.

Towards a Warm Holographic Equation of State by an Einstein–Maxwell-Dilaton Model

The holographic Einstein–Maxwell-dilaton model is employed to map state-of-the-art lattice QCD thermodynamics data from the temperature (T) axis towards the baryon–chemical potential (μB) axis and aims to gain a warm equation of state (EoS) of deconfined QCD matter which can be supplemented with a cool and confined part suitable for subsequent compact (neutron) star (merger) investigations. The model exhibits a critical end point (CEP) at TCEP=O(100) MeV and μBCEP=500…700 MeV with an emerging first-order phase transition (FOPT) curve which extends to large values of μB without approaching the μB axis. We consider the impact and peculiarities of the related phase structure on the EoS for the employed dilaton potential and dynamical coupling parameterizations. These seem to prevent the design of an overall trustable EoS without recourse to hybrid constructions.

Thermodynamics of quark matter with multiquark clusters in an effective Beth-Uhlenbeck type approach

We describe multiquark clusters in quark matter within a Beth-Uhlenbeck approach in a background gluon field coupled to the underlying chiral quark dynamics using the Polyakov gauge. An effective potential for the traced Polyakov loop is used to establish the center symmetry of the SU(3) color which suppresses colored states and its dynamical breaking as an aspect of the confinement/deconfinement transition. Quark confinement is modeled by a large quark mass in vacuum which is motivated by a confining density functional approach. A multiquark cluster containing n quarks and antiquarks is described as a binary composite of smaller subclusters n1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$n_1$$\\end{document} and n2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$n_2$$\\end{document} (n1+n2=n\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$n_1+n_2=n$$\\end{document}). It has a spectrum consisting of a bound state and a scattering state continuum. For the corresponding cluster-cluster phase shifts we use simple ansätze that capture the Mott dissociation of clusters as a function of temperature and chemical potential. We go beyond the simple “step-up-step-down” model that ignores continuum correlations and introduce an improved model that includes them in a generic form. In order to explain the model, we restrict ourselves here to the cases where the cluster size is 1≤n≤6\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1 \\le n \\le 6$$\\end{document}. A striking result is the suppression of the abundance of colored multiquark clusters at low temperatures by the coupling to the Polyakov loop and their importance for a quantitative description of lattice QCD thermodynamics at non-vanishing baryochemical potentials. An important ingredient are Polyakov-loop generalized distribution functions of n-quark clusters which are derived here for the first time. Within our approach we calculate thermodynamic properties such as baryon density and entropy. We demonstrate that the limits of a hadron resonance gas at low temperatures and O(g2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {O}(g^2)$$\\end{document} perturbative QCD at high temperatures are correctly reproduced. A comparison with lattice calculations shows that our model is able to give a unified, systematic approach to describe properties of the quark-gluon-hadron system.

Holographic transport coefficients and jet energy loss for the hot and dense quark-gluon plasma

We employ an Einstein-Maxwell-dilaton model, based on the gauge/gravity correspondence, to obtain the thermodynamics and transport properties for the hot and dense quark-gluon plasma. The model, which is constrained to reproduce lattice QCD thermodynamics at zero density, predicts a critical point and a first order line at finite temperature and density, is used to quantify jet energy loss through simulations of high-energy collision events.

Transport coefficients of the quark-gluon plasma at the critical point and across the first-order line

A bottom-up Einstein-Maxwell-Dilaton holographic model is used to compute, for the first time, the behavior of several transport coefficients of the hot and baryon-rich strongly coupled quark-gluon plasma at the critical point and also across the first-order phase transition line in the phase diagram. The observables under study are of the shear and bulk viscosities, baryon diffusion, thermal conductivity, the jet quenching parameter $\hat{q}$, as well as the heavy-quark drag force and Langevin diffusion coefficients. These calculations provide a phenomenologically promising estimate for these coefficients, given that our model quantitatively reproduces lattice QCD thermodynamics results, both at zero and finite baryon density, besides naturally incorporating the nearly-perfect fluidity of the quark-gluon plasma. We find that the diffusion of baryon charge, and also the shear and bulk viscosities, are suppressed with increasing baryon density, indicating that the medium becomes even closer to perfect fluidity at large densities. On the other hand, the jet quenching parameter and the heavy-quark momentum diffusion are enhanced with increasing density. The observables display a discontinuity gap when crossing the first-order phase transition line, while developing an infinite slope at the critical point. The transition temperatures associated with different transport coefficients differ in the crossover region but are found to converge at the critical point.

The Role of Strange Quasiquarks in Transport Properties of the QGP

We study the role of dynamical quarks in the transport properties of the deconfined hot medium utilizing the quasiparticle kinetic approach. The system is then composed of the quasiparticle excitations, whose interactions are encoded in their dynamical masses with temperature dependence specified by the effective coupling extracted from lattice QCD thermodynamics. Evaluating the temperature and flavor profiles of the shear and bulk viscosities, as well as the electrical conductivity we examine how different particle species, in particular strange quarks, modify the transport parameters of the deconfined matter.

Impact of quark quasiparticles on transport coefficients in hot QCD

We study the bulk and shear viscosity and the electrical conductivity in a quasiparticle approach to Yang-Mills theory and QCD with light and strange quarks to assess the dynamical role of quarks in transport properties at finite temperature. The interactions with a hot medium are embodied in effective masses of the constituents through a temperature-dependent running coupling extracted from the lattice QCD thermodynamics. In Yang-Mills theory, the bulk viscosity to entropy density ratio exhibits a non-monotonous structure around the phase transition temperature. In QCD, this is totally dissolved because of a substantial contribution from quark quasiparticles. The bulk to shear viscosity ratio near the phase transition behaves consistently to the scaling with the speed of sound derived in the AdS/CFT approach, whereas at high temperature it obeys the same parametric dependence as in perturbation theory. Thus, the employed quasiparticle model is adequate to capture the transport properties in the weak and strong coupling regimes of the theory. This feature is not altered by including dynamical quarks which, however, retards the system from restoring conformal invariance. We also examine the individual flavor contributions to the electrical conductivity and show that the obtained behavior agrees qualitatively well with the recent results of lattice simulations and with a class of phenomenological approaches.

The possibility of twin star solutions in a model based on lattice QCD thermodynamics

The properties of compact stars and in particular the existence of twin star solutions are investigated within an effective model that is constrained by lattice QCD thermodynamics. The model is modified at large baryon densities to incorporate a large variety of scenarios of first order phase transitions to a phase of deconfined quarks. This is achieved by matching two different variants of the bag model equation of state, in order to estimate the role of the Bag model parameters on the appearance of a second family of neutron stars. The produced sequences of neutron stars are compared with modern constrains on stellar masses, radii, and tidal deformability from astrophysical observations and gravitational wave analyses. It is found that those scenarios in our analysis, in which a third family of stars appeared due to the deconfinement transition, are disfavored from astrophysical constraints.

Jet quenching parameters in strongly coupled anisotropic plasmas in the presence of magnetic fields

I use the holographic gauge/gravity duality to systematically calculate the jet quenching parameters in strongly coupled anisotropic plasmas in the presence of external magnetic fields. The magnetic field breaks down spatial rotation symmetry from $SO(3)$ to $SO(2)$, leading to the presence of multiple anisotropic jet quenching parameters, which are evaluated here in two quite different holographic settings. One of them corresponds to a top-down deformation of the strongly coupled $\mathcal{N} = 4$ Super Yang-Mills plasma triggered by an external magnetic field, while the other one is a bottom-up Einstein-Maxwell-Dilaton model of phenomenological relevance for high energy peripheral heavy ion collisions, since it is able to provide a quantitative description of $(2+1)$-flavors lattice QCD thermodynamics with physical quark masses at zero and nonzero magnetic fields. I find for both models an overall enhancement of all the anisotropic jet quenching parameters with increasing magnetic fields. Moreover, I also conclude that for both models transverse momentum broadening is larger in transverse directions than in the direction of the magnetic field. Since these conclusions are shown to hold for two rather different holographic setups at finite temperature and magnetic fields, they are suggested as fairly robust features of strongly coupled anisotropic magnetized plasmas.

QCD at high density: Equation of state for nuclear collisions and neutron stars

A unified chiral mean field approach is presented for QCD thermodynamics in a wide range of temperatures and densities. The model simultaneously gives a satisfactory description of lattice QCD thermodynamics and fulfills nuclear matter and astrophysical constraints. The resulting equation of state can be incorporated in relativistic fluid-dynamical simulations of heavy-ion collisions and neutron stars mergers. Access to different regions of the QCD phase diagram can be obtained in simulations of heavy-ion data and observations of neutron star mergers.

Lattice QCD Thermodynamics and RHIC-BES Particle Production within Generic Nonextensive Statistics

The current status of implementing Tsallis (nonextensive) statistics on high-energy physics is briefly reviewed. The remarkably low freezeout-temperature, which apparently fails to reproduce the first-principle lattice QCD thermodynamics and the measured particle ratios, etc. is discussed. The present work suggests a novel interpretation for the so-called {\it "Tsallis-temperature"}. It is proposed that the low Tsallis-temperature is due to incomplete implementation of Tsallis algebra though exponential and logarithmic functions to the high-energy particle-production. Substituting Tsallis algebra into grand-canonical partition-function of the hadron resonance gas model seems not assuring full incorporation of nonextensivity or correlations in that model. The statistics describing the phase-space volume, the number of states and the possible changes in the elementary cells should be rather modified due to interacting correlated subsystems, of which the phase-space is consisting. Alternatively, two asymptotic properties, each is associated with a scaling function, are utilized to classify a generalized entropy for such a system with large ensemble (produced particles) and strong correlations. Both scaling exponents define equivalence classes for all interacting and noninteracting systems and unambiguously characterize any statistical system in its thermodynamic limit. We conclude that the nature of lattice QCD simulations is apparently extensive and accordingly the Boltzmann-Gibbs statistics is fully fulfilled. Furthermore, we found that the ratios of various particle yields at extreme high and extreme low energies of RHIC-BES is likely nonextensive but not necessarily of Tsallis type.

Role of repulsive interactions in the interplay with missing strange resonances

The standard implementation of the HRG model has been shown to be unable to describe all the available data on QCD matter. Here we show the balance of repulsive and attractive hadronic interactions on QCD thermodynamics through observables both calculated by lattice simulations and measured in experiment. Attractive interactions are mediated by resonance formation, which are here implemented through extra states predicted by the Quark Model, while repulsive interactions are modelled by means of Excluded Volume (EV) effects. Informations on flavour dependent effective sizes are extracted. It is found that EV effects are present in lattice QCD thermodynamics, and are essential for a comprehensive description of higher order fluctuations of conserved charges.

Lattice QCD thermodynamics up to the perturbative regime

We study the thermodynamics of the quark gluon plasma with lattice simulations in the continuum limit up to 1 GeV temperature where we show that a perturbative description already applies. We calculate the effect of the presence of charm quark in the equation of state and also describe the topological features of quantum chromodynamics. This conference contribution is based on Nature 539 (2016) no. 7627, 69–71 and Phys. Rev. D 92 (11) (2015), 114505.

Momentum transport in strongly coupled anisotropic plasmas in the presence of strong magnetic fields

We present a holographic perspective on momentum transport in strongly coupled, anisotropic non-Abelian plasmas in the presence of strong magnetic fields. We compute the anisotropic heavy quark drag forces and Langevin diffusion coefficients and also the anisotropic shear viscosities for two different holographic models, namely, a top-down deformation of strongly coupled $\mathcal{N} = 4$ Super-Yang-Mills (SYM) theory triggered by an external Abelian magnetic field, and a bottom-up Einstein-Maxwell-dilaton (EMD) model which is able to provide a quantitative description of lattice QCD thermodynamics with $(2+1)$-flavors at both zero and nonzero magnetic fields. We find that, in general, energy loss and momentum diffusion through strongly coupled anisotropic plasmas are enhanced by a magnetic field being larger in transverse directions than in the direction parallel to the magnetic field. Moreover, the anisotropic shear viscosity coefficient is smaller in the direction of the magnetic field than in the plane perpendicular to the field, which indicates that strongly coupled anisotropic plasmas become closer to the perfect fluid limit along the magnetic field. We also present, in the context of the EMD model, holographic predictions for the entropy density and the crossover critical temperature in a wider region of the $(T,B)$ phase diagram that has not yet been covered by lattice simulations. Our results for the transport coefficients in the phenomenologically realistic magnetic EMD model could be readily used as inputs in numerical codes for magnetohydrodynamics.

Recent results on lattice QCD thermodynamics

I review recent results on QCD thermodynamics from lattice simulations. In particular, I will focus on the QCD equation of state at zero and finite chemical potential, the curvature of the phase diagram and fluctuations of conserved charges. The latter are compared to experimental data, to the purpose of extracting the chemical freeze-out temperature and chemical potential from first principles.

Pion and kaon in the Beth-Uhlenbeck approach

In the present work the Mott effect for pions and kaons is described within a Beth- Uhlenbeck approach on the basis of the PNJL model. The contribution of these degrees of freedom to the thermodynamics is encoded in the temperature dependence of their phase shifts. A comparison with results from Nf = 2 + 1 lattice QCD thermodynamics is performed.

Mott-hadron resonance gas and lattice QCD thermodynamics

We present an effective model for the generic behaviour of hadron masses and phase shifts at finite temperature which shares basic features with recent developments within the PNJL model for correlations in quark matter. On this basis we obtain the transition between a hadron resonance gas phase and the quark gluon plasma in the spirit of the generalized Beth-Uhlenbeck approach where the Mott dissociation of hadrons is encoded in the hadronic phase shifts. Here we restrict ourselves to low-lying hadronic channels and perform a discussion of recent lattice QCD thermodynamics results from this perspective. We find agreement in the asymptotic regions while for the description of the transition itself the inclusion of further hadronic channels as well as a selfconsistent determination of the continuum thresholds is required.

Cross-over versus first-order phase transition in holographic gravity–single-dilaton models of QCD thermodynamics

A dilaton potential is adjusted to recently confirmed lattice QCD thermodynamics data in the temperature range $(0.7 \ldots 3.5) T_c$ where $T_c = 155 \text{MeV}$ is the pseudo-critical temperature. The employed holographic model is based on a gravity--single-field dilaton dual. We discuss conditions for enforcing (for the pure gluon plasma) or avoiding (for the QCD quark-gluon plasma) a first-order phase transition, but still keeping a softest point (minimum of sound velocity).

Electric Conductivity of the QGP

The transport coefficients of strongly interacting matter have attracted a great interest in the field of Quark-Gluon Plasma (QGP). In this work we compute electric conductivity σel solving numerically the Relativistic Boltzmann Transport (RBT) equation in a uniform box with periodic boundary conditions considering 2—body scatterings. We compare numerical results obtained using two methods, Green-Kubo correlator and E-field method, with analytic formulas in Relaxation Time Approximation (RTA). We present results for the realistic case of the QGP system considering both a quasi-particle model tuned to lattice QCD thermodynamics as well as the case of a pQCD gas with a running coupling constant. Calculations based on RTA underestimate σel of about 60%.

Lattice QCD Thermodynamics with Physical Quark Masses

Over the past few years, new physics methods and algorithms as well as the latest supercomputers have enabled the study of the QCD thermodynamic phase transition using lattice gauge theory numerical simulations and allowing unprecedented control over systematic errors. This progress is largely a consequence of the ability to perform continuum extrapolations with physical quark masses. We review recent progress in lattice QCD thermodynamics, focusing mainly on results that benefit from the use of physical quark masses: the crossover temperature, the equation of state, and fluctuations of the quark number susceptibilities. In addition, we emphasize calculations that are directly relevant to the study of relativistic heavy-ion collisions at RHIC and the LHC.