The utility of quantum computers for simulating lattice gauge theories is currently limited by the noisiness of the physical hardware. Various quantum error mitigation strategies exist to reduce the statistical and systematic uncertainties in quantum simulations via improved algorithms and analysis strategies. We perform quantum simulations of Z_{2} gauge theory with matter to study the efficacy and interplay of different error mitigation methods: readout error mitigation, randomized compiling, rescaling, and dynamical decoupling. We compute Minkowski correlation functions in this confining gauge theory and extract the mass of the lightest spin-1 state from fits to their time dependence. Quantum error mitigation extends the range of times over which our correlation function calculations are accurate by a factor of 6 and is therefore essential for obtaining reliable masses.
The utility of quantum computers for simulating lattice gauge theories is currently limited by the noisiness of the physical hardware. Various quantum error mitigation strategies exist to reduce the statistical and systematic uncertainties in quantum simulations via improved algorithms and analysis strategies. We perform quantum simulations of Z_{2} gauge theory with matter to study the efficacy and interplay of different error mitigation methods: readout error mitigation, randomized compiling, rescaling, and dynamical decoupling. We compute Minkowski correlation functions in this confining gauge theory and extract the mass of the lightest spin-1 state from fits to their time dependence. Quantum error mitigation extends the range of times over which our correlation function calculations are accurate by a factor of 6 and is therefore essential for obtaining reliable masses.
Simulating Z2 lattice gauge theory on a quantum computer
The utility of quantum computers for simulating lattice gauge theories is currently limited by the noisiness of the physical hardware. Various quantum error mitigation strategies exist to reduce the statistical and systematic uncertainties in quantum simulations via improved algorithms and analysis strategies. We perform quantum simulations of 1+1dZ2 gauge theory with matter to study the efficacy and interplay of different error mitigation methods: readout error mitigation, randomized compiling, rescaling, and dynamical decoupling. We compute Minkowski correlation functions in this confining gauge theory and extract the mass of the lightest spin-1 state from fits to their time dependence. Quantum error mitigation extends the range of times over which our correlation function calculations are accurate by a factor of six and is therefore essential for obtaining reliable masses.
This research paper provides a comprehensive review of quantum simulation methods for lattice gauge theories during the current era of noisy intermediate-scale quantum computing. It focuses on the fundamental framework used to understand strong interactions in high-energy physics and explores how quantum devices can overcome the limitations of traditional computational methods, such as the sign problem found in classical simulations. The manuscript is built around three primary pillars that are critical for near-term progress. The first pillar involves encoding strategies, describing the movement away from simple mapping techniques toward more advanced, compact methods that significantly decrease the number of required quantum bits by focusing only on physically relevant states and exploiting local gauge constraints. The second pillar analyzes specific variational quantum algorithms that use flexible circuit structures to find low-energy states, comparing different approaches based on their circuit depth and suitability for current hardware. The third pillar evaluates various error mitigation techniques used to extract accurate information from imperfect quantum hardware, such as methods that amplify and then remove noise or those that use physical symmetries to detect and discard errors. Additionally, the paper compares different hardware technologies—including those based on superconducting circuits, trapped atoms, and ions—to weigh their respective benefits, such as gate speed and connectivity, against their drawbacks. It also introduces emerging hybrid approaches that combine quantum circuits with artificial neural networks to improve simulation power and state representation. Ultimately, the review identifies pathways toward achieving practical advantages in the field and discusses the challenges that remain for scaling these systems to larger sizes.
Quantum simulation provides a promising framework for studying non-perturbative phenomena in lattice gauge theories that are difficult to access using classical computational methods. In particular, classical Monte Carlo approaches suffer from the sign problem, limiting their applicability to real-time dynamics and finite-density regimes. Quantum computers, by directly encoding quantum dynamics, offer a natural pathway to overcoming these limitations. In this thesis, we study quantum simulations of a truncated SU(2) lattice gauge theory using a Hamiltonian formulation, representing an important non-Abelian step toward quantum chromodynamics. We focus on three-dimensional gauge-theory simulations in the noisy intermediate-scale quantum (NISQ) era, where limited qubit resources pose a significant challenge. To address this, we introduce the triamond lattice, a three-dimensional lattice geometry that enforces local gauge invariance while providing an economical encoding of gauge degrees of freedom, thereby reducing qubit overhead compared to conventional lattice constructions. We derive the corresponding Hamiltonian, construct physical gauge-invariant states, and implement quantum simulations using quantum imaginary time evolution to prepare ground states, alongside real-time evolution methods to study dynamical behavior. Error-mitigation techniques are employed to extract reliable physical observables from noisy quantum hardware. Our results demonstrate that the triamond lattice enables resource-efficient three-dimensional quantum gauge simulations and represents a practical step toward scalable quantum simulations of non-perturbative gauge theories.
Quantum simulation of lattice gauge theories offers a path to first-principles calculations of real-time dynamics and finite-density phenomena that are inaccessible to classical Euclidean Monte Carlo methods. Realizing this potential requires a pipeline of controlled approximations---lattice discretization, truncation, digitization, state preparation, time evolution, measurement, and continuum extrapolation---each contributing errors and computational costs. General algorithmic asymptotic complexity arguments capture only part of the cost picture: the practical resources required for end-to-end simulation depend on numerical pre-factors, the cost of the subroutines that various algorithms call, and structural features of the lattice theory, all of which the asymptotic scaling leaves unspecified. This thesis develops contributions that aim to make such dependencies explicit and quantitative. We begin with a rigorous analysis of the asymptotic gate complexity of product formulas and Quantum Signal Processing (QSP) for the time evolution of general classes of lattice Hamiltonians, identifying the parameter regimes in which each method is favorable and estimating the associated pre-factors numerically. Since the cost of QSP-based simulation is dominated (in regimes of practical interest) by the block-encoding subroutine it calls, we also develop several novel block-encoding constructions tailored to bosonic lattice field theories and show that signal-processing techniques can be repurposed as highly efficient block-encoding subroutines. We then turn to the problem of truncating the Kogut--Susskind electric-field Hilbert space. Existing error bounds rely on energy-conservation arguments and require the truncation to grow with simulation time. By exploiting Hilbert space fragmentation we derive the tightest error bounds known to date: they decay factorially in the truncation parameter and are independent of simulation time. We validate them numerically for U(1) pure gauge theory and the Schwinger model. Finally, we address the continuum limit in the presence of approximate time evolution. For product formulas, we show that Trotter errors correspond to irrelevant operators that vanish in the continuum limit, which yields a renormalization trajectory independent of the Trotter step size. More generally, we introduce the Statistically-Bounded Time Evolution (SBTE) protocol, an algorithm-independent framework that treats approximation errors as a systematic uncertainty to be driven below the working statistical precision. This decouples the continuum-limit prescription from the choice of time-evolution algorithm and provides a basis for end-to-end cost comparisons across simulation methods.
This thesis investigates advances in quantum simulation methods for lattice gauge theories, with a focus on (2+1)-dimensional quantum electrodynamics (QED), which shares important features with (3+1)-dimensional quantum chromodynamics (QCD), such as confinement and asymptotic freedom. Two main projects are developed. The first project investigates the ongoing coupling in a compact pure U(1) gauge theory using a step-scaling function approach. A quantum variational method allows calculations over a wide range of couplings. Plaquette expectation value results are compared with classical Monte Carlo simulations to obtain, in principle, the physical value of the lattice spacing. This method can potentially be extended to the study of QED, which also includes fermionic matter. The second project investigates the static potential between two charges, analyzing the Coulomb, confining, and string-breaking regimes. It involves the development of a tailored variational approach that incorporates system symmetries and mutual information to represent fermionic degrees of freedom. Results from a trapped-ion quantum device agree well with classical simulations, validating the method. Furthermore, this work evaluates superconducting quantum devices, noting issues with noise but also demonstrating improvements in hardware and potential through error limitation. Overall, this work underscores the feasibility of variational quantum algorithms for lattice gauge theory simulations and demonstrates their promise for advancing theoretical physics with quantum computers. This thesis explores advancements in quantum simulation methods for lattice gauge theories, focusing on (2+1)-dimensional Quantum Electrodynamics (QED), which shares key features with (3+1)-dimensional Quantum Chromodynamics (QCD), such as confinement and asymptotic freedom. Two main projects are developed. The first investigates the running coupling in a compact U(1) pure gauge theory using a step scaling function approach. A quantum variational method enables computation across a wide range of couplings. The plaquette expectation value is matched with classical Monte Carlo simulations to obtain, in principle, the physical value of the lattice spacing. This method can be eventually extended to study QED, where fermionic matter is included. The second project examines the static potential between two charges, analyzing Coulomb, confining, and string-breaking regimes. It involves developing a custom variational Ansatz that incorporates system symmetries and mutual information to represent fermionic degrees of freedom. Results from a trapped-ion quantum device align well with exact classical simulations, validating the method. Additionally, the thesis evaluates superconducting devices, identifying noise challenges but noting improvements in hardware and the potential of error mitigation. Overall, this work highlights the viability of variational quantum algorithms for lattice gauge theory simulations and demonstrates their promise for advancing theoretical physics with quantum computing.
Quantum computers offer the possibility to implement lattice gauge theory in Minkowski rather than Euclidean spacetime, thus allowing calculations of processes that evolve in real time. In this work, calculations within SU(2) pure gauge theory are able to show the motion of an excitation traveling across a spatial lattice in real time. This is accomplished by using a simple yet powerful method for error mitigation, where the original circuit is used both forward and backward in time. For a two-plaquette lattice, meaningful results are obtained from a circuit containing hundreds of CNOT gates. The same method is used for a five-plaquette lattice, where calculations show that residual systematic effects can be reduced through follow-up mitigation.