Bootstrap for lattice Yang-Mills theory

We study the SU(infinity) lattice Yang-Mills theory at the dimensions D = 2, 3, 4 via the numerical bootstrap method. It combines the loop equations, with a cutoff Lmax on the maximal length of loops, and positivity conditions on certain matrices of Wilson loop averages. Our algorithm is inspired by the pioneering paper of P. D. Anderson and M. Kruczenski [Nucl. Phys. B921, 702 (2017)] but it is significantly more efficient, as it takes into account the symmetries of the lattice theory and uses the relaxation procedure in line with our previous work on matrix bootstrap. We thus obtain rigorous upper and lower bounds on the plaquette average at various couplings and dimensions. For D = 4; Lmax = 16 the lower bound data appear to be close to the Monte Carlo data in the strong coupling phase and the upper bound data in the weak coupling phase reproduce well the 3-loop perturbation theory. Our results suggest that this bootstrap approach can provide a tangible alternative to the, so far uncontested, Monte Carlo approach.