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

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.

💡 Research Summary

The paper presents a comprehensive account of the theoretical developments, computational strategy, and data analysis that enabled the first non‑perturbative determination of the QCD equation of state (EoS) up to temperatures of order the electroweak scale (≈165 GeV). The authors focus on three‑flavour massless QCD in the temperature interval 3 GeV ≤ T ≤ 165 GeV and achieve a relative precision of about 1 % or better across the whole range.

Two methodological breakthroughs make this possible. First, the authors define lines of constant physics (LCPs) by matching the non‑perturbatively measured Schrödinger‑functional coupling (\bar g^2_{\rm SF}(\mu)) in a finite volume to its known continuum running at the scale (\mu=1/L_0). This fixes the bare gauge coupling (g_0) as a function of the lattice spacing (a) for any chosen temperature, allowing simulations at several lattice resolutions (L0/a = 4, 6, 8, 10) and a controlled continuum extrapolation. Second, they employ shifted (or moving‑frame) boundary conditions characterized by a spatial shift vector (\xi). In this setup the derivative of the free‑energy density with respect to (\xi) yields directly the entropy density, \