Efficient construction of mathbb{Z}₂ gauge-invariant bases for the Quantum Minimally Entangled Typical Thermal States algorithm

In quantum computations of gauge theories at finite temperature and finite density, it is challenging to enforce Gauss's law for all states contributing to the thermal ensemble. While various techniques for implementing gauge constraints have been proposed, they often involve practical trade-offs. In this work, we adopt the Quantum Minimally Entangled Typical Thermal States (QMETTS) algorithm for mathbb{Z}₂ gauge-constrained systems, which allows us to capture thermal equilibrium states with chemical potential while mitigating these trade-offs. To ensure that gauge invariance is preserved throughout the procedure while maintaining computational efficiency, we derive the specific measurement bases within the algorithm. Furthermore, since the estimation of expectation values on quantum hardware is inherently noisy, we rigorously account for shot noise in estimating expectation values, and propose a sampling method that is more efficient than those in previous works. We validate our approach numerically by studying a (1+1)-dimensional mathbb{Z}₂ gauge theory coupled to staggered fermions. Our proposed algorithm reproduces the correct equilibrium states at finite temperature and finite density.