Density of states method for symplectic gauge theories at finite temperature

We study the finite-temperature behavior of the <a:math xmlns:a="http://www.w3.org/1998/Math/MathML" display="inline"><a:mi>S</a:mi><a:mi>p</a:mi><a:mo stretchy="false">(</a:mo><a:mn>4</a:mn><a:mo stretchy="false">)</a:mo></a:math> Yang-Mills lattice theory in four dimensions, by applying the logarithmic linear relaxation algorithm. We demonstrate the presence of coexisting (metastable) phases, when the system is in the proximity of the transition. We measure observables such as the free energy, the expectation value of the plaquette operator and of the Polyakov loop, as well as the specific heat, and the Binder cumulant. We use these results to obtain a high-precision measurement of the critical coupling at the confinement-deconfinement transition, and assess its systematic uncertainty, for one value of the lattice extent in the time direction. Furthermore, we perform an extensive study of the finite-volume behavior of the lattice system, by repeating the measurements for fixed lattice time extent, while increasing the spatial size of the lattice. We hence characterize the first-order transition on the lattice and present the first results in the literature on this theory for the infinite volume extrapolation of lattice quantities related to latent heat and interface tension. Gauge theories with <e:math xmlns:e="http://www.w3.org/1998/Math/MathML" display="inline"><e:mi>S</e:mi><e:mi>p</e:mi><e:mo stretchy="false">(</e:mo><e:mn>4</e:mn><e:mo stretchy="false">)</e:mo></e:math> group have been proposed as new dark sectors to provide a fundamental origin for the current phenomenological evidence of dark matter. A phase transition at high temperature, in such a new dark sector, occurring in the early Universe, might have left a relic stochastic background of gravitational waves. Our results represent a milestone toward establishing whether such a new physics signal is detectable in future experiments, as they enter the calculation of the parameters, <i:math xmlns:i="http://www.w3.org/1998/Math/MathML" display="inline"><i:mi>α</i:mi></i:math> and <k:math xmlns:k="http://www.w3.org/1998/Math/MathML" display="inline"><k:mi>β</k:mi></k:math>, controlling the power spectrum of gravitational waves. We also outline the process needed in the continuum extrapolation of our measurements and test its feasibility on one additional choice of temporal extent of the lattice.