Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
This paper introduces a novel framework leveraging a hierarchy of semidefinite programs to derive rigorous lower bounds on the ground-state energy of lattice gauge theories, overcoming limitations of Monte Carlo methods (sign problem) and variational approaches (upper bounds only). The methodology extends to providing certified intervals for extended observables like Wilson loops and mesonic strings, demonstrated on $\mathbb{Z}_2$ gauge theories with fermionic matter and pure gauge. The approach yields highly precise certified intervals, with relative spreads as low as $0.04\%$ for a $(1+1)$-dimensional system, offering a computational tool for exploring many-body physics and validating quantum simulations.
Rigorous lower bounds on lattice gauge theory ground states and observables are now certifiably obtainable, even in regimes intractable for Monte Carlo.
Gauge theories describe fundamental interactions in the standard model of particle physics and build effective theories in condensed matter physics. Being notoriously hard to simulate, these theories are commonly regularized as lattice gauge theories which can be evaluated with high precision through Monte Carlo algorithms or variational methods. However, Monte Carlo algorithms are not applicable in all regimes due to the sign problem; while variational methods, quantum and classical alike, only give upper-bounds on the ground-state energy. Their precision crucially depends on the chosen ansatz. Here, we present a framework based on a hierarchy of semidefinite programs to yield increasingly good lower-bounds on the ground-state energy and propagate these bounds to extended observables like Wilson loops and mesonic strings. We demonstrate the numerical capabilities by applying the algorithm to a (1+1)-dimensional $\mathbb{Z}_2$ theory with dynamic fermionic matter and its pure-gauge version in (2+1)-dimensions. In both cases, we obtain certified intervals for the ground-state energy. For a one-dimensional system with 129 sites, the certified interval has a relative spread of $0.04\%$. For a two-dimensional system with 97 sites (corresponding to a $7 \times 8$ lattice), the relative spread is $2.08\%$. Beyond energy estimates, the proposed method provides rigorous bounds on both short- and long-range observables such as Wilson loops over multiple plaquettes and the magnitude of the mesonic string. The method can be both used as a computational microscope into many-body physics and a way to certify results of quantum simulators.