Kalman Szabo

The order of the quantum chromodynamics transition predicted by the standard model of particle physics

Quantum chromodynamics (QCD) is the theory of the strong interaction, explaining (for example) the binding of three almost massless quarks into a much heavier proton or neutron—and thus most of the mass of the visible Universe. The standard model of particle physics predicts a QCD-related transition that is relevant for the evolution of the early Universe. At low temperatures, the dominant degrees of freedom are colourless bound states of hadrons (such as protons and pions). However, QCD is asymptotically free, meaning that at high energies or temperatures the interaction gets weaker and weaker, causing hadrons to break up. This behaviour underlies the predicted cosmological transition between the low-temperature hadronic phase and a high-temperature quark–gluon plasma phase (for simplicity, we use the word ‘phase’ to characterize regions with different dominant degrees of freedom). Despite enormous theoretical effort, the nature of this finite-temperature QCD transition (that is, first-order, second-order or analytic crossover) remains ambiguous. Here we determine the nature of the QCD transition using computationally demanding lattice calculations for physical quark masses. Susceptibilities are extrapolated to vanishing lattice spacing for three physical volumes, the smallest and largest of which differ by a factor of five. This ensures that a true transition should result in a dramatic increase of the susceptibilities. No such behaviour is observed: our finite-size scaling analysis shows that the finite-temperature QCD transition in the hot early Universe was not a real phase transition, but an analytic crossover (involving a rapid change, as opposed to a jump, as the temperature varied). As such, it will be difficult to find experimental evidence of this transition from astronomical observations.

The QCD equation of state with dynamical quarks

The present paper concludes our investigation on the QCD equation of state with 2 + 1 staggered flavors and one-link stout improvement. We extend our previous study [JHEP01 (2006) 089] by choosing even finer lattices. Lattices with Nt=6, 8 and 10 are used, and the continuum limit is approached by checking the results at Nt= 12. A Symanzik improved gauge and a stout-link improved staggered fermion action is utilized. We use physical quark masses, that is, for the lightest staggered pions and kaons we fix the mπ/fK and mK/fK ratios to their experimental values. The pressure, the interaction measure, the energy and entropy density and the speed of sound are presented as functions of the temperature in the range 100 ... 1000MeV. We give estimates for the pion mass dependence and for the contribution of the charm quark. We compare our data to the equation of state obtained by the “hotQCD” collaboration.

Is there still any T c mystery in lattice QCD? Results with physical masses in the continuum limit III

The present paper concludes our investigations on the QCD cross-over transition temperatures with 2+1 staggered flavours and one-link stout improvement. We extend our previous two studies [Phys. Lett. B643 (2006) 46, JHEP 0906:088 (2009)] by choosing even finer lattices (Nt = 16) and we work again with physical quark masses. The new results on this broad cross-over are in complete agreement with our earlier ones. We compare our findings with the published results of the hotQCD collaboration. All these results are confronted with the predictions of the Hadron Resonance Gas model and Chiral Perturbation Theory for temperatures below the transition region. Our results can be reproduced by using the physical spectrum in these analytic calculations. The findings of the hotQCD collaboration can be recovered by using a distorted spectrum which takes into account lattice discretization artifacts and heavier than physical quark masses. This analysis provides a simple explanation for the observed discrepancy in the transition temperatures between our and the hotQCD collaborations.

Ab Initio Determination of Light Hadron Masses

More than 99% of the mass of the visible universe is made up of protons and neutrons. Both particles are much heavier than their quark and gluon constituents, and the Standard Model of particle physics should explain this difference. We present a full ab initio calculation of the masses of protons, neutrons, and other light hadrons, using lattice quantum chromodynamics. Pion masses down to 190 mega–electron volts are used to extrapolate to the physical point, with lattice sizes of approximately four times the inverse pion mass. Three lattice spacings are used for a continuum extrapolation. Our results completely agree with experimental observations and represent a quantitative confirmation of this aspect of the Standard Model with fully controlled uncertainties.

The equation of state in lattice QCD: with physical quark masses towards the continuum limit

The equation of state of QCD at vanishing chemical potential as a function of temperature is determined for two sets of lattice spacings. Coarser lattices with temporal extension of Nt=4 and finer lattices of Nt=6 are used. Symanzik improved gauge and stout-link improved staggered fermionic actions are applied. The results are given for physical quark masses both for the light quarks and for the strange quark. Pressure, energy density, entropy density, quark number susceptibilities and the speed of sound are presented.

High-precision scale setting in lattice QCD

A bstractScale setting is of central importance in lattice QCD. It is required to predict dimensional quantities in physical units. Moreover, it determines the relative lattice spacings of computations performed at different values of the bare coupling, and this is needed for extrapolating results into the continuum. Thus, we calculate a new quantity, w0, for setting the scale in lattice QCD, which is based on the Wilson flow like the scale t0 (M. Luscher, JHEP 08 (2010) 071). It is cheap and straightforward to implement and compute. In particular, it does not involve the delicate fitting of correlation functions at asymptotic times. It typically can be determined on the few per-mil level. We compute its continuum extrapolated value in 2 + 1-flavor QCD for physical and non-physical pion and kaon masses, to allow for mass-independent scale setting even away from the physical mass point. We demonstrate its robustness by computing it with two very different actions (one of them with staggered, the other with Wilson fermions) and by showing that the results agree for physical quark masses in the continuum limit.

Lattice QCD as a video game

Abstract The speed, bandwidth and cost characteristics of todays PC graphics cards make them an attractive target as general purpose computational platforms. High performance can be achieved also for lattice simulations but the actual implementation can be cumbersome. This paper outlines the architecture and programming model of modern graphics cards for the lattice practitioner with the goal of exploiting these chips for Monte Carlo simulations. Sample code is also given.

QCD equation of state at nonzero chemical potential: continuum results with physical quark masses at order μ 2

A bstractWe determine the equation of state of QCD for nonzero chemical potentials via a Taylor expansion of the pressure. The results are obtained for Nf = 2 + 1 flavors of quarks with physical masses, on various lattice spacings. We present results for the pressure, interaction measure, energy density, entropy density, and the speed of sound for small chemical potentials. At low temperatures we compare our results with the Hadron Resonance Gas model. We also express our observables along trajectories of constant entropy over particle number. A simple parameterization is given (the Matlab/Octave script parameterization.m, submitted to the arXiv along with the paper), which can be used to reconstruct the observables as functions of T and μ, or as functions of T and S/N.

Precision SU(3) lattice thermodynamics for a large temperature range

A bstractWe present the equation of state (pressure, trace anomaly, energy density and entropy density) of the SU(3) gauge theory from lattice field theory in an unprecedented precision and temperature range. We control both finite size and cut-off effects. The studied temperature window (0.7…1000Tc) stretches from the glueball dominated system into the perturbative regime, which allows us to discuss the range of validity of these approaches. We also determine the preferred renormalization scale of the Hard Thermal Loop scheme and we fit the unknown g6 order perturbative coefficient at extreme high temperatures T > 100Tc. We furthermore quantify the nonperturbative contribution to the trace anomaly using a simple functional form. Our high precision data allows one to have a complete theoretical description of the equation of state from T = 0 all the way to the phase transition, through the transition region into the perturbative regime up to the Stefan-Boltzmann limit. We will discuss this description, too.

Lattice QCD at the physical point: light quark masses

Abstract Ordinary matter is described by six fundamental parameters: three couplings (gravitational, electromagnetic and strong) and three masses: the electronʼs ( m e ) and those of the up ( m u ) and down ( m d ) quarks. An additional mass enters through quantum fluctuations: the strange quark mass ( m s ). The three couplings and m e are known with an accuracy of better than a few per mil. Despite their importance, m u , m d (their average m u d ) and m s are far less accurately known. Here we determine them with a precision below 2% by performing ab initio lattice quantum chromodynamics (QCD) calculations, in which all systematics are controlled. We use pion and quark masses down to (and even below) their physical values, lattice sizes of up to 6 fm, and five lattice spacings to extrapolate to continuum spacetime. All necessary renormalizations are performed nonperturbatively.