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:2508.00061 · Quantum

Truncation uncertainties for accurate quantum simulations of lattice gauge theories

Anthony N. Ciavarella, Siddharth Hariprakash, Jad C. Halimeh, C. Bauer

Published Sep 24, 2026• 0 citations
arXivHTMLPDF

Abstract

The encoding of lattice gauge theories onto quantum computers requires a discretization of the gauge field's Hilbert space on each link, which presents errors with respect to the Kogut–Susskind limit. In the electric basis, Hilbert space fragmentation has recently been shown to limit the excitation of large electric fields. Here, we leverage this to develop a formalism for estimating the size of truncation errors in the electric basis. Generically, the truncation error falls off as a factorial of the field truncation. Examples of this formalism are applied to the Schwinger model and a pure U(1) lattice gauge theory. For reasonable choices of parameters, we improve on previous error estimates by a factor of 10 306 .

DOI
10.22331/q-2026-09-24-2216

Related papers

Towards Precision Quantum Simulation of Lattice Gauge Theories

Siddharth Hariprakash

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.

•May 1, 2026•0 citations

Impact of Irrep Truncations on Quantum Simulations of QCD

Patrick Draper, Drishti Gupta

Quantum simulations of lattice gauge theories with continuous gauge groups require a choice of basis (how states of the gauge theory map to the qubit computational basis) and truncation scheme (which of the infinitely many gauge theory basis states should be omitted in the map to the finite dimensional qubit Hilbert space.) The most common choices are adapted to an electric basis for the gauge theory, where the basis states for each link are associated with irreducible representations of the gauge group. Electric bases and truncations are generally advantageous for preserving gauge symmetries and simulating the strong coupling regime, but their effectiveness in the weak coupling regimes of physical interest are less well understood. In this work, we compute relative probabilities $p(R')/p(R)$ for finding a link in the irreps $R,R'$ in the ground state of SU(3) lattice gauge theory using Euclidean Monte Carlo methods and weak-coupling lattice perturbation theory. The observables that compute these relative probabilities are constructed by introducing seams into an anisotropic lattice, and by a similar construction we are also able to measure relative probabilities of different site singlets (ways that electric flux braids together into a gauge singlet) at each vertex. With these methods we can quantify, for example, the importance of higher electric flux states in the weak coupling ground state. These results directly inform what truncations are needed to accurately capture non-abelian gauge dynamics near the continuum limit.

hep-latCross-listedarXiv:2610.23681•Oct 19, 2026•0 citations

Provably accurate simulation of gauge theories and bosonic systems

Yu Tong, Victor V. Albert, Jarrod R. McClean +2 more

Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit Λ on each local quantum number, and if the initial state has low local quantum numbers, then an error at most ϵ can be achieved by choosing Λ to scale polylogarithmically with ϵ−1 , an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on Λ that achieves accuracy ϵ , obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy ϵ by truncating local quantum numbers at Λ=polylog(ϵ−1) .

arXiv:2110.06942•Quantum•Sep 22, 2022•72 citations

A new basis for Hamiltonian SU(2) simulations

C. Bauer, Irian D’Andrea, Marat Freytsis +1 more

Due to rapidly improving quantum computing hardware, Hamiltonian simulations of relativistic lattice field theories have seen a resurgence of attention. This computational tool requires turning the formally infinite-dimensional Hilbert space of the full theory into a finite-dimensional one. For gauge theories, a widely-used basis for the Hilbert space relies on the representations induced by the underlying gauge group, with a truncation that keeps only a set of the lowest dimensional representations. This works well at large bare gauge coupling, but becomes less efficient at small coupling, which is required for the continuum limit of the lattice theory. In this work, we develop a new basis suitable for the simulation of an SU(2) lattice gauge theory in the maximal tree gauge. In particular, we show how to perform a Hamiltonian truncation so that the eigenvalues of both the magnetic and electric gauge-fixed Hamiltonian are mostly preserved, which allows for this basis to be used at all values of the coupling. Little prior knowledge is assumed, so this may also be used as an introduction to the subject of Hamiltonian formulations of lattice gauge theories.

arXiv:2307.11829•Jul 21, 2023•0 citations

Efficient truncations of SU ( N c ) lattice gauge theory for quantum simulation

Anthony N. Ciavarella, Iván Mauricio Burbano, C. Bauer

Quantum simulations of lattice gauge theories offer the potential to directly study the nonperturbative dynamics of quantum chromodynamics, but naive analyses suggest that they require large computational resources. Large N c expansions are performed to order 1 / N c to simplify the Hamiltonian of pure SU ( N c ) lattice gauge theories. A reformulation of the electric basis is introduced with a truncation strategy based on the construction of local Krylov subspaces with plaquette operators. Numerical simulations show that these truncated Hamiltonians are consistent with traditional lattice calculations at relatively small couplings. It is shown that the computational resources required for quantum simulation of time evolution generated by these Hamiltonians is 17–19 orders of magnitude smaller than previous approaches, provided that the truncations in this work can reach lattice spacings in three-dimensional simulations comparable to the two-dimensional simulations performed.

arXiv:2503.11888•Physical review. D/Physical review. D.•Aug 26, 2025•8 citations

Generic Hilbert space fragmentation in Kogut-Susskind lattice gauge theories

Anthony N. Ciavarella, C. Bauer, Jad C. Halimeh

At the heart of quantum many-body physics lies the understanding of mechanisms that avoid quantum thermalization in an isolated system quenched far from equilibrium. A prominent example is Hilbert space fragmentation, which has recently emerged as an ergodicity-breaking mechanism in constrained spin models. Here, we show that Kogut-Susskind formulations of lattice gauge theories in $d+1\mathrm{D}$ ($d$ spatial and one temporal dimensions) give rise to Hilbert space fragmentation, and discuss possible implications for understanding continuum physics. Our findings not only prove that lattice gauge theories are a natural platform for Hilbert space fragmentation, they also serve as a guide to the conditions under which these models can be faithfully used to infer the thermalization properties of quantum chromodynamics.

arXiv:2502.03533•Physical review. D/Physical review. D.•Nov 12, 2025•8 citations
The encoding of lattice gauge theories onto quantum computers requires a discretization of the gauge field's Hilbert space on each link, which presents errors with respect to the Kogut–Susskind limit. In the electric basis, Hilbert space fragmentation has recently been shown to limit the excitation of large electric fields. Here, we leverage this to develop a formalism for estimating the size of truncation errors in the electric basis. Generically, the truncation error falls off as a factorial of the field truncation. Examples of this formalism are applied to the Schwinger model and a pure U(1) lattice gauge theory. For reasonable choices of parameters, we improve on previous error estimates by a factor of 10 306 .
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.
Quantum simulations of lattice gauge theories with continuous gauge groups require a choice of basis (how states of the gauge theory map to the qubit computational basis) and truncation scheme (which of the infinitely many gauge theory basis states should be omitted in the map to the finite dimensional qubit Hilbert space.) The most common choices are adapted to an electric basis for the gauge theory, where the basis states for each link are associated with irreducible representations of the gauge group. Electric bases and truncations are generally advantageous for preserving gauge symmetries and simulating the strong coupling regime, but their effectiveness in the weak coupling regimes of physical interest are less well understood. In this work, we compute relative probabilities p(R′)/p(R)p(R')/p(R)p(R′)/p(R) for finding a link in the irreps R,R′R,R'R,R′ in the ground state of SU(3) lattice gauge theory using Euclidean Monte Carlo methods and weak-coupling lattice perturbation theory. The observables that compute these relative probabilities are constructed by introducing seams into an anisotropic lattice, and by a similar construction we are also able to measure relative probabilities of different site singlets (ways that electric flux braids together into a gauge singlet) at each vertex. With these methods we can quantify, for example, the importance of higher electric flux states in the weak coupling ground state. These results directly inform what truncations are needed to accurately capture non-abelian gauge dynamics near the continuum limit.
Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit Λ on each local quantum number, and if the initial state has low local quantum numbers, then an error at most ϵ can be achieved by choosing Λ to scale polylogarithmically with ϵ−1 , an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on Λ that achieves accuracy ϵ , obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy ϵ by truncating local quantum numbers at Λ=polylog(ϵ−1) .
Due to rapidly improving quantum computing hardware, Hamiltonian simulations of relativistic lattice field theories have seen a resurgence of attention. This computational tool requires turning the formally infinite-dimensional Hilbert space of the full theory into a finite-dimensional one. For gauge theories, a widely-used basis for the Hilbert space relies on the representations induced by the underlying gauge group, with a truncation that keeps only a set of the lowest dimensional representations. This works well at large bare gauge coupling, but becomes less efficient at small coupling, which is required for the continuum limit of the lattice theory. In this work, we develop a new basis suitable for the simulation of an SU(2) lattice gauge theory in the maximal tree gauge. In particular, we show how to perform a Hamiltonian truncation so that the eigenvalues of both the magnetic and electric gauge-fixed Hamiltonian are mostly preserved, which allows for this basis to be used at all values of the coupling. Little prior knowledge is assumed, so this may also be used as an introduction to the subject of Hamiltonian formulations of lattice gauge theories.
Quantum simulations of lattice gauge theories offer the potential to directly study the nonperturbative dynamics of quantum chromodynamics, but naive analyses suggest that they require large computational resources. Large N c expansions are performed to order 1 / N c to simplify the Hamiltonian of pure SU ( N c ) lattice gauge theories. A reformulation of the electric basis is introduced with a truncation strategy based on the construction of local Krylov subspaces with plaquette operators. Numerical simulations show that these truncated Hamiltonians are consistent with traditional lattice calculations at relatively small couplings. It is shown that the computational resources required for quantum simulation of time evolution generated by these Hamiltonians is 17–19 orders of magnitude smaller than previous approaches, provided that the truncations in this work can reach lattice spacings in three-dimensional simulations comparable to the two-dimensional simulations performed.
At the heart of quantum many-body physics lies the understanding of mechanisms that avoid quantum thermalization in an isolated system quenched far from equilibrium. A prominent example is Hilbert space fragmentation, which has recently emerged as an ergodicity-breaking mechanism in constrained spin models. Here, we show that Kogut-Susskind formulations of lattice gauge theories in d+1Dd+1\mathrm{D}d+1D (ddd spatial and one temporal dimensions) give rise to Hilbert space fragmentation, and discuss possible implications for understanding continuum physics. Our findings not only prove that lattice gauge theories are a natural platform for Hilbert space fragmentation, they also serve as a guide to the conditions under which these models can be faithfully used to infer the thermalization properties of quantum chromodynamics.