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 with unprecedented control over systemat
ic errors. This is largely a consequence of the ability to perform continuum extrapolations with physical quark masses. Here we review recent progress in lattice QCD thermodynamics, focussing 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 place a special emphasis on calculations that are directly relevant to the study of relativistic heavy ion collisions at RHIC and the LHC.
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_{tau}=6,~8,~10$, and $12$. We show that these data can be reliably extrapolated to the continuum li
mit 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~ {rm MeV} leq T leq 163$ MeV, defined by the chiral transition, is $epsilon_c=(0.18-0.5)~{rm GeV}/{rm fm}^3$, $i.e.$, $(1.2-3.1) epsilon_{rm 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 $r_0$ from the static quark potential and $w_0$ from the gradient flow.
Using numerical simulations of lattice QCD we calculate the effect of an external magnetic field on the equation of state of the quark-gluon plasma. The results are obtained using a Taylor expansion of the pressure with respect to the magnetic field
for the first time. The coefficients of the expansion are computed to second order in the magnetic field. Our setup for the external magnetic field avoids complications arising from toroidal boundary conditions, making a Taylor series expansion straightforward. This study is exploratory and is meant to serve as a proof of principle.
A calculation of the ratio of leptonic decay constants f_{K^+}/f_{pi^+} makes possible a precise determination of the ratio of CKM matrix elements |V_{us}|/|V_{ud}| in the Standard Model, and places a stringent constraint on the scale of new physics
that would lead to deviations from unitarity in the first row of the CKM matrix. We compute f_{K^+}/f_{pi^+} numerically in unquenched lattice QCD using gauge-field ensembles recently generated that include four flavors of dynamical quarks: up, down, strange, and charm. We analyze data at four lattice spacings a ~ 0.06, 0.09, 0.12, and 0.15 fm with simulated pion masses down to the physical value 135 MeV. We obtain f_{K^+}/f_{pi^+} = 1.1947(26)(37), where the errors are statistical and total systematic, respectively. This is our first physics result from our N_f = 2+1+1 ensembles, and the first calculation of f_{K^+}/f_{pi^+} from lattice-QCD simulations at the physical point. Our result is the most precise lattice-QCD determination of f_{K^+}/f_{pi^+}, with an error comparable to the current world average. When combined with experimental measurements of the leptonic branching fractions, it leads to a precise determination of |V_{us}|/|V_{ud}| = 0.2309(9)(4) where the errors are theoretical and experimental, respectively.
We present results from our simulations of quantum chromodynamics (QCD) with four flavors of quarks: u, d, s, and c. These simulations are performed with a one-loop Symanzik improved gauge action, and the highly improved staggered quark (HISQ) action
. We are generating gauge configurations with four values of the lattice spacing ranging from 0.06 fm to 0.15 fm, and three values of the light quark mass, including the value for which the Goldstone pion mass is equal to the physical pion mass. We discuss simulation algorithms, scale setting, taste symmetry breaking, and the autocorrelations of various quantities. We also present results for the topological susceptibility which demonstrate the improvement of the HISQ configurations relative to those generated earlier with the asqtad improved staggered action.
We fit lattice-QCD data for light-pseudoscalar masses and decay constants, from HISQ configurations generated by MILC, to SU(3) staggered chiral perturbation theory. At present such fits have rather high values of chi^2/d.o.f., possibly due to the la
ck of ensembles with lighter-than-physical sea strange-quark masses. We propose solutions to this problem for future work. We also perform simple linear interpolations near the physical point on two ensembles with different lattice spacings, and obtain the preliminary result (f_K / f_pi)^phys = 1.1872(41) in the continuum limit.
In calculations of the hyperfine splitting in charmonium, the contributions of the disconnected diagrams are considered small and are typically ignored. We aim to estimate nonperturbatively the size of the resulting correction, which may eventually b
e needed in high precision calculations of the charmonium spectrum. We study this problem in the quenched and unquenched QCD cases. On dynamical ensembles the disconnected charmonium propagators contain light modes which complicate the extraction of the signal at large distances. In the fully quenched case, where there are no such light modes, the interpretation of the signal is simplified. We present results from lattices with $aapprox 0.09$ fm and $aapprox 0.06$ fm.
Recently, the Fermilab heavy-quark action was extended to include dimension-six and -seven operators in order to reduce the discretization errors. In this talk, we present results of the first numerical simulations with this action (the OK action), w
here we study the masses of the quarkonium and heavy-light systems. We calculate combinations of masses designed to test improvement and compare results obtained with the OK action to their counterparts obtained with the clover action. Our preliminary results show a clear improvement.
We study the lattice spacing dependence, or scaling, of physical quantities using the highly improved staggered quark (HISQ) action introduced by the HPQCD/UKQCD collaboration, comparing our results to similar simulations with the asqtad fermion acti
on. Results are based on calculations with lattice spacings approximately 0.15, 0.12 and 0.09 fm, using four flavors of dynamical HISQ quarks. The strange and charm quark masses are near their physical values, and the light-quark mass is set to 0.2 times the strange-quark mass. We look at the lattice spacing dependence of hadron masses, pseudoscalar meson decay constants, and the topological susceptibility. In addition to the commonly used determination of the lattice spacing through the static quark potential, we examine a determination proposed by the HPQCD collaboration that uses the decay constant of a fictitious unmixed s bar s pseudoscalar meson. We find that the lattice artifacts in the HISQ simulations are much smaller than those in the asqtad simulations at the same lattice spacings and quark masses.
We extend our work on QCD thermodynamics with 2+1 quark flavors at nonzero chemical potential to finer lattices with $N_t=6$. We study the equation of state and other thermodynamic quantities, such as quark number densities and susceptibilities, and
compare them with our previous results at $N_t=4$. We also calculate the effects of the addition of the charm and bottom quarks on the equation of state at zero and nonzero chemical potential. These effects are important for cosmological studies of the early Universe.
