Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2610.06272

Two-Sided Bounds on Ground State Properties in Lattice Gauge Theories

Eloïc Vallée, Ariel Kelman, Julius Mildenberger, Jordi Tura, Patrick Emonts

quant-phhep-latSubmitted Oct 5, 2026• Announced Oct 6, 2026• 0 citations
arXivHTMLPDF

Abstract

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.

Comments
14 pages, 10 figures

Related papers

Variational Monte Carlo algorithm for lattice gauge theories with continuous gauge groups: A study of ( 2 + 1 ) -dimensional compact QED with dynamical fermions at finite density

Julian Bender, Patrick Emonts, J. I. Cirac

Lattice gauge theories coupled to fermionic matter account for many interesting phenomena in both high-energy physics and condensed-matter physics. Certain regimes, e.g., at finite fermion density, are difficult to simulate with traditional Monte Carlo algorithms due to the so-called sign problem. We present a variational, sign-problem-free Monte Carlo method for lattice gauge theories with continuous gauge groups and apply it to $(2+1)$-dimensional compact QED with dynamical fermions at finite density. The variational ansatz is formulated in the full gauge-field basis, i.e., without having to resort to truncation schemes for the $\text{U}(1)$ gauge-field Hilbert space. The ansatz consists of two parts: first, a pure gauge part based on Jastrow-type ansatz states (which can be connected to certain neural-network ansatz states) and, second, a fermionic part based on gauge-field-dependent fermionic Gaussian states. These are designed in such a way that the gauge-field integral over all fermionic Gaussian states is gauge-invariant and at the same time still efficiently tractable. To ensure the validity of the method we benchmark the pure gauge part of the ansatz against another variational method and the full ansatz against an existing Monte Carlo simulation where the sign problem is absent. Moreover, in limiting cases where the exact ground state is known we show that our ansatz is able to capture this behavior. Finally, we study a sign-problem affected regime by probing density-induced phase transitions.

arXiv:2304.05916•Physical Review Research•Nov 7, 2023•13 citations

Ground state search of interacting boson–fermion systems using variational Monte Carlo

Pranay Naredi, Stefan Kühn, Enrico Brehm

We study lattice QED in $2+1$ dimensions coupled to Wilson fermions in the Hamiltonian formulation and develop a sign-problem-free variational Monte Carlo algorithm. Using a continuous basis to describe the gauge degrees of freedom and fermionic Gaussian states for the matter fields, we construct a variational ansatz for the ground state. The parameters of the ansatz are optimized using gradient descent and/or imaginary time evolution based on the time-dependent variational principle. For a lattice with $N$ matter sites, the computational complexity for sampling a gauge configuration is $\mathcal{O}(N^2\log N)$, and the complexity for computing the energy per gauge configuration is $\mathcal{O}(N^4)$. Comparing our method to results from exact diagonalization on small lattices, we demonstrate good agreement. In addition, we present preliminary results for the ground state for two fermion flavors, $N_f=2$, in the presence of a non-zero chemical potential, a regime in which the conventional action-based Monte Carlo methods suffer from the sign problem.

International Symposium on Lattice Field Theory (LATTICE)•May 30, 2026•0 citations

Certifying Ground-State Properties of Many-Body Systems

Jie Wang, Jacopo Surace, Irénée Frérot +4 more

A ubiquitous problem in quantum physics is to understand the ground-state properties of many-body systems. Confronted with the fact that exact diagonalization quickly becomes impossible when increasing the system size, variational approaches are typically employed as a scalable alternative: Energy is minimized over a subset of all possible states and then different physical quantities are computed over the solution state. Despite remarkable success, rigorously speaking, all that variational methods offer are upper bounds on the ground-state energy. On the other hand, so-called relaxations of the ground-state problem based on semidefinite programming represent a complementary approach, providing lower bounds to the ground-state energy. However, in their current implementation, neither variational nor relaxation methods offer provable bound on other observables in the ground state beyond the energy. In this work, we show that the combination of the two classes of approaches can be used to derive certifiable bounds on the value of any observable in the ground state, such as correlation functions of arbitrary order, structure factors, or order parameters. We illustrate the power of this approach in paradigmatic examples of 1D and 2D spin- 1 / 2 Heisenberg models. To improve the scalability of the method, we exploit the symmetries and sparsity of the considered systems to reach sizes of hundreds of particles at much higher precision than previous works. Our analysis therefore shows how to obtain certifiable bounds on many-body ground-state properties beyond energy in a scalable way. Published by the American Physical Society 2024

arXiv:2310.05844•Physical Review X•Jul 11, 2024•15 citations

Finding the ground state of a lattice gauge theory with fermionic tensor networks: A 2+1D Z2 demonstration

Patrick Emonts, Ariel Kelman, Umberto Borla +3 more

Tensor network states, and in particular projected entangled pair states (PEPS) have been a strong ansatz for the variational study of complicated quantum many-body systems, thanks to their built-in entanglement entropy area law. In this work, we use a special kind of PEPS---gauged Gaussian fermionic PEPS (GGFPEPS)---to find the ground state of $2+1$ dimensional pure ${\mathbb{Z}}_{2}$ lattice gauge theories for a wide range of coupling constants. We do so by combining PEPS methods with Monte-Carlo computations, allowing for efficient contraction of the PEPS and computation of correlation functions. Previously, such numerical computations involved the calculation of the Pfaffian of a matrix scaling with the system size, forming a severe bottleneck; in this work we show how to overcome this problem. This paves the way for applying the method we propose and benchmark here to other gauge groups, higher dimensions, and models with fermionic matter, in an efficient, sign-problem-free way.

arXiv:2211.00023•Physical review. D/Physical review. D.•Jan 4, 2023•27 citations

Local manipulation and measurement of nonlocal many-body operators in lattice gauge theory quantum simulators

Erez Zohar

Lattice gauge theories form a very successful framework for studying nonperturbative gauge field physics, in particular in quantum chromodynamics. Recently, their quantum simulation on atomic and solid-state platforms has been discussed, aiming at overcoming some of the difficulties still faced by the conventional approaches (such as the sign problem and real time evolution). While the actual implementations of a lattice gauge theory on a quantum simulator may differ in terms of the simulating system and its properties, they are all directed at studying similar physical phenomena, requiring the measurement of nonlocal observables, due to the local symmetry of gauge theories. In this work, general schemes for measuring such nonlocal observables (Wilson loops and mesonic string operators) in general lattice gauge theory quantum simulators that are based merely on local operations are proposed.

arXiv:1911.11156•Physical review. D/Physical review. D.•Feb 24, 2020•12 citations

Neural quantum states for non-Abelian lattice gauge theories with dynamical fermions

Gabriel Rouxinol, Julian Bender, Michele Grossi +2 more

Determining the ground state of non-Abelian lattice gauge theories coupled to dynamical fermions is key to understanding confinement and the phase structure of gauge--matter systems. We present a variational Monte Carlo framework for the ground state of the untruncated fully-continuous SU$(2)$ lattice gauge theory coupled to dynamical staggered fermions on an $L\times L$ square lattice. We work in the magnetic basis with a neural-network representation of the gauge wavefunction. The fermions are described by a gauge-covariant Gaussian fermionic correction built on a fixed Néel reference state where, for each sampled gauge configuration $\mathbf{U}$, the correction is generated by a Hermitian operator. This operator is constructed from short Wilson lines and the eigenvectors of the mass--hopping Hamiltonian, with number of variational parameters polynomial in the system size. This Gaussian structure also gives analytical expressions for all fermionic contributions to the energy and related observables in terms of the fermion occupation matrix. The results are validated against strong-coupling perturbation theory, where they recover the expected effective antiferromagnetic spin Hamiltonian. Using this framework, we map a coarse ground state phase diagram in the plane of independent electric and magnetic couplings $(g^2, λ)$ and show that a hysteresis analysis can identify the existence of phase transitions. Restoring the physical relation $λ=4/g^2$, we characterize how increasing the system size and changing the electric coupling $g^2$ move the state away from the reference Néel state, for lattice sizes $L=4,6,8$. More broadly, the method offers a sign-problem-free variational framework for continuous non-Abelian gauge groups with dynamical matter that should extend to other matter content and higher-dimensional lattices.

hep-latCross-listedarXiv:2607.26131•Jul 30, 2026•0 citations
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 Z2\mathbb{Z}_2Z2​ 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%0.04\%0.04%. For a two-dimensional system with 97 sites (corresponding to a 7×87 \times 87×8 lattice), the relative spread is 2.08%2.08\%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.
Lattice gauge theories coupled to fermionic matter account for many interesting phenomena in both high-energy physics and condensed-matter physics. Certain regimes, e.g., at finite fermion density, are difficult to simulate with traditional Monte Carlo algorithms due to the so-called sign problem. We present a variational, sign-problem-free Monte Carlo method for lattice gauge theories with continuous gauge groups and apply it to (2+1)(2+1)(2+1)-dimensional compact QED with dynamical fermions at finite density. The variational ansatz is formulated in the full gauge-field basis, i.e., without having to resort to truncation schemes for the U(1)\text{U}(1)U(1) gauge-field Hilbert space. The ansatz consists of two parts: first, a pure gauge part based on Jastrow-type ansatz states (which can be connected to certain neural-network ansatz states) and, second, a fermionic part based on gauge-field-dependent fermionic Gaussian states. These are designed in such a way that the gauge-field integral over all fermionic Gaussian states is gauge-invariant and at the same time still efficiently tractable. To ensure the validity of the method we benchmark the pure gauge part of the ansatz against another variational method and the full ansatz against an existing Monte Carlo simulation where the sign problem is absent. Moreover, in limiting cases where the exact ground state is known we show that our ansatz is able to capture this behavior. Finally, we study a sign-problem affected regime by probing density-induced phase transitions.
We study lattice QED in 2+12+12+1 dimensions coupled to Wilson fermions in the Hamiltonian formulation and develop a sign-problem-free variational Monte Carlo algorithm. Using a continuous basis to describe the gauge degrees of freedom and fermionic Gaussian states for the matter fields, we construct a variational ansatz for the ground state. The parameters of the ansatz are optimized using gradient descent and/or imaginary time evolution based on the time-dependent variational principle. For a lattice with NNN matter sites, the computational complexity for sampling a gauge configuration is O(N2log⁡N)\mathcal{O}(N^2\log N)O(N2logN), and the complexity for computing the energy per gauge configuration is O(N4)\mathcal{O}(N^4)O(N4). Comparing our method to results from exact diagonalization on small lattices, we demonstrate good agreement. In addition, we present preliminary results for the ground state for two fermion flavors, Nf=2N_f=2Nf​=2, in the presence of a non-zero chemical potential, a regime in which the conventional action-based Monte Carlo methods suffer from the sign problem.
A ubiquitous problem in quantum physics is to understand the ground-state properties of many-body systems. Confronted with the fact that exact diagonalization quickly becomes impossible when increasing the system size, variational approaches are typically employed as a scalable alternative: Energy is minimized over a subset of all possible states and then different physical quantities are computed over the solution state. Despite remarkable success, rigorously speaking, all that variational methods offer are upper bounds on the ground-state energy. On the other hand, so-called relaxations of the ground-state problem based on semidefinite programming represent a complementary approach, providing lower bounds to the ground-state energy. However, in their current implementation, neither variational nor relaxation methods offer provable bound on other observables in the ground state beyond the energy. In this work, we show that the combination of the two classes of approaches can be used to derive certifiable bounds on the value of any observable in the ground state, such as correlation functions of arbitrary order, structure factors, or order parameters. We illustrate the power of this approach in paradigmatic examples of 1D and 2D spin- 1 / 2 Heisenberg models. To improve the scalability of the method, we exploit the symmetries and sparsity of the considered systems to reach sizes of hundreds of particles at much higher precision than previous works. Our analysis therefore shows how to obtain certifiable bounds on many-body ground-state properties beyond energy in a scalable way. Published by the American Physical Society 2024
Tensor network states, and in particular projected entangled pair states (PEPS) have been a strong ansatz for the variational study of complicated quantum many-body systems, thanks to their built-in entanglement entropy area law. In this work, we use a special kind of PEPS---gauged Gaussian fermionic PEPS (GGFPEPS)---to find the ground state of 2+12+12+1 dimensional pure Z2{\mathbb{Z}}_{2}Z2​ lattice gauge theories for a wide range of coupling constants. We do so by combining PEPS methods with Monte-Carlo computations, allowing for efficient contraction of the PEPS and computation of correlation functions. Previously, such numerical computations involved the calculation of the Pfaffian of a matrix scaling with the system size, forming a severe bottleneck; in this work we show how to overcome this problem. This paves the way for applying the method we propose and benchmark here to other gauge groups, higher dimensions, and models with fermionic matter, in an efficient, sign-problem-free way.
Lattice gauge theories form a very successful framework for studying nonperturbative gauge field physics, in particular in quantum chromodynamics. Recently, their quantum simulation on atomic and solid-state platforms has been discussed, aiming at overcoming some of the difficulties still faced by the conventional approaches (such as the sign problem and real time evolution). While the actual implementations of a lattice gauge theory on a quantum simulator may differ in terms of the simulating system and its properties, they are all directed at studying similar physical phenomena, requiring the measurement of nonlocal observables, due to the local symmetry of gauge theories. In this work, general schemes for measuring such nonlocal observables (Wilson loops and mesonic string operators) in general lattice gauge theory quantum simulators that are based merely on local operations are proposed.
Determining the ground state of non-Abelian lattice gauge theories coupled to dynamical fermions is key to understanding confinement and the phase structure of gauge--matter systems. We present a variational Monte Carlo framework for the ground state of the untruncated fully-continuous SU(2)(2)(2) lattice gauge theory coupled to dynamical staggered fermions on an L×LL\times LL×L square lattice. We work in the magnetic basis with a neural-network representation of the gauge wavefunction. The fermions are described by a gauge-covariant Gaussian fermionic correction built on a fixed Néel reference state where, for each sampled gauge configuration U\mathbf{U}U, the correction is generated by a Hermitian operator. This operator is constructed from short Wilson lines and the eigenvectors of the mass--hopping Hamiltonian, with number of variational parameters polynomial in the system size. This Gaussian structure also gives analytical expressions for all fermionic contributions to the energy and related observables in terms of the fermion occupation matrix. The results are validated against strong-coupling perturbation theory, where they recover the expected effective antiferromagnetic spin Hamiltonian. Using this framework, we map a coarse ground state phase diagram in the plane of independent electric and magnetic couplings (g2,λ)(g^2, λ)(g2,λ) and show that a hysteresis analysis can identify the existence of phase transitions. Restoring the physical relation λ=4/g2λ=4/g^2λ=4/g2, we characterize how increasing the system size and changing the electric coupling g2g^2g2 move the state away from the reference Néel state, for lattice sizes L=4,6,8L=4,6,8L=4,6,8. More broadly, the method offers a sign-problem-free variational framework for continuous non-Abelian gauge groups with dynamical matter that should extend to other matter content and higher-dimensional lattices.