QCD Equation of State at very high temperature: computational strategy, simulations and data analysis

TL;DR

This paper non-perturbatively determines the QCD Equation of State at very high temperatures using innovative computational strategies, achieving high precision and validating results against perturbative predictions, thus advancing understanding of thermal QCD.

Contribution

It introduces a novel computational framework combining finite-volume coupling and shifted boundary conditions for precise high-temperature QCD thermodynamics.

Findings

Achieved 1% accuracy in the QCD Equation of State at temperatures up to 165 GeV.

Validated non-perturbative results against perturbative predictions, highlighting the importance of non-perturbative effects.

Demonstrated the method's robustness through extensive consistency checks.

Abstract

We present a detailed account of the theoretical progress and the computational strategy that led to the non-perturbative determination of the QCD Equation of State at temperatures up to the electroweak scale reported in [Phys. Rev. Lett. 134, 201904 (2025)]. The two key ingredients that make such a calculation feasible with controlled uncertainties are: (i) the definition of lines of constant physics through the running of a non-perturbatively defined finite-volume coupling across a wide range of energy scales, and (ii) the use of shifted boundary conditions which allow a direct determination of the entropy density thus without the need for a zero-temperature subtraction. Considering the case of QCD with $N_f =3$ massless flavours in the temperature interval between 3 GeV and 165 GeV, we describe the numerical strategy based on integrating in the bare coupling and quark mass, the perturbative improvement of lattice observables, the optimization of numerical simulations, and the continuum extrapolation. Extensive consistency checks, including finite-volume and topological-freezing effects, confirm the robustness of the method. The final results have a relative accuracy of about $1\%$ or better, and the errors are dominated by the statistical fluctuations of the Monte Carlo ensembles. We also compare our non-perturbative results with predictions from standard and hard thermal loop perturbation theory showing that at the level of $\%$-precision contributions beyond those known, including non-perturbative ones due to ultrasoft modes, are relevant up to the highest temperatures explored. The methodological framework is general and readily applicable to QCD with four and five massive quark flavours and to other thermal observables, paving the way for systematic non-perturbative studies of thermal QCD at very high temperatures.