We present a high-precision lattice calculation of the equation of state in the confining phase of SU(2) Yang-Mills theory. We show that the results are described very well by a gas of massive, non-interacting glueballs, provided one assumes an exponentially growing Hagedorn spectrum. The latter can be derived within an effective bosonic closed-string model, leading to a parameter-free theoretical prediction, which is in perfect agreement with our lattice results. Furthermore, when applied to SU(3) Yang-Mills theory, this effective model accurately describes the lattice results reported by Borsányi et al. in JHEP 07 (2012) 056.

Jarzynski's theorem is a well-known equality in statistical mechanics, which relates fluctuations in the work performed during a non-equilibrium transformation of a system, to the free-energy difference between two equilibrium ensembles. In this article, we apply Jarzynski's theorem in lattice gauge theory, for two examples of challenging computational problems, namely the calculation of interface free energies and the determination of the equation of state. We conclude with a discussion of further applications of interest in QCD and in other strongly coupled gauge theories, in particular for the Schrödinger functional and for simulations at finite density using reweighting techniques.

A precise lattice determination of the equation of state in SU(3) Yang-Mills theory is carried out by means of a simulation algorithm, based on Jarzynski's theorem, that allows one to compute physical quantities in thermodynamic equilibrium, by driving the field configurations of the system out of equilibrium. The physical results and the computational efficiency of the algorithm are compared with other state-of-the-art lattice calculations, and the extension to full QCD with dynamical fermions and to other observables is discussed.

The evaluation of these diagrams is required for many phenomenologically interesting observables, but suffers from large statistical errors due to the vacuum and random-noise contributions to their variances. Motivated by a theoretical analysis of the variances, we introduce a new family of stochastic estimators of single-propagator traces built upon a frequency splitting combined with a hopping expansion of the quark propagator, and test their efficiency in two-flavour QCD with pions as light as 190 MeV. The use of these estimators reduces the cost of the computation by one to two orders of magnitude over standard estimators depending on the fermion bilinear. As a concrete application, we show the impact of these findings on the computation of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment.