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

Towards Precision Quantum Simulation of Lattice Gauge Theories

Siddharth Hariprakash

Published May 1, 2026• 0 citations
OpenAlex

Abstract

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.

Related papers

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

Obtaining continuum physics from dynamical simulations of Hamiltonian lattice gauge theories

Kane, Christopher F., Hariprakash, Siddharth, Bauer, Christian W.

Taking the continuum limit is essential for extracting physical observables from quantum simulations of lattice gauge theories. Achieving the correct continuum limit requires careful control of all systematic uncertainties, including those arising from approximate implementations of the time evolution operator. In this work, we review existing approaches based on renormalization techniques, and point out their limitations. To overcome these limitations, we introduce a new general framework -- the Statistically-Bounded Time Evolution (SBTE) protocol -- for rigorously controlling the impact of approximate time evolution on the continuum limit. The central insight is that, since exact time evolution introduces no UV divergences, errors from approximate evolution can be treated as a source of systematic uncertainty that can be neglected if reduced below the working statistical uncertainty. We show that, using the SBTE protocol, which prescribes driving the approximate time evolution error below the working statistical uncertainty, leads to a simplified renormalization procedure. Furthermore, we show that, due to the existence of rigorous error bounds, one can guarantee a priori that such errors are negligible and do not affect the continuum limit. Ultimately, our protocol lays the foundation for performing systematic and fair comparisons between different simulation algorithms for lattice gauge theory simulations.

arXiv:2506.16559•Jun 19, 2025•0 citations

Truncation uncertainties for accurate quantum simulations of lattice gauge theories

Anthony N. Ciavarella, Siddharth Hariprakash, Jad C. Halimeh +1 more

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 .

arXiv:2508.00061•Quantum•Sep 24, 2026•0 citations

Efficient Quantum Gauge Simulations with the Triamond

Ali Hosseinzadehkavaki

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.

•Dec 5, 2025•0 citations

Strategies for simulating the time evolution of Hamiltonian lattice field theories

Siddharth Hariprakash, Neel S. Modi, Michael Kreshchuk +2 more

Simulating the time evolution of quantum field theories given some Hamiltonian $H$ requires developing algorithms for implementing the unitary operator ${\mathrm{e}}^{\ensuremath{-}iHt}$. A variety of techniques exist that accomplish this task, with the most common technique used so far being Trotterization, which is a special case of the application of a product formula. However, other techniques exist that promise better asymptotic scaling in certain parameters of the theory being simulated, the most efficient of which are based on the concept of block encoding. In this work we study the performance of such algorithms in simulating lattice field theories. We derive and compare the asymptotic gate complexities of several commonly used simulation techniques in application to Hamiltonian lattice field theories. Using the scalar ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\varphi}}}^{4}$ theory as a test, we also perform numerical studies and compare the gate costs required by product formulas and signal-processing-based techniques to simulate time evolution. For the latter, we use the linear combination of unitaries (LCU) construction augmented with the quantum Fourier transform circuit to switch between the field and momentum eigenbases, which leads to immediate order-of-magnitude improvement in the cost of preparing the block encoding. This paper also includes a pedagogical review of the techniques used, in particular product formulas, LCU, qubitization, quantum signal processing, as well as the technique for simulating geometrically-local Hamiltonians developed by Haah, Hastings, Kothari, and Low.

Physical Review A•Feb 12, 2025•8 citations

A resource efficient approach for quantum and classical simulations of gauge theories in particle physics

Jan F. Haase, Luca Dellantonio, Alessio Celi +4 more

Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present limitations of MCMC techniques may be overcome by Hamiltonian-based simulations on classical or quantum devices, which further provide the potential to address questions that lay beyond the capabilities of the current approaches. However, for continuous gauge groups, Hamiltonian-based formulations involve infinite-dimensional gauge degrees of freedom that can solely be handled by truncation. Current truncation schemes require dramatically increasing computational resources at small values of the bare couplings, where magnetic field effects become important. Such limitation precludes one from `taking the continuous limit' while working with finite resources. To overcome this limitation, we provide a resource-efficient protocol to simulate LGTs with continuous gauge groups in the Hamiltonian formulation. Our new method allows for calculations at arbitrary values of the bare coupling and lattice spacing. The approach consists of the combination of a Hilbert space truncation with a regularization of the gauge group, which permits an efficient description of the magnetically-dominated regime. We focus here on Abelian gauge theories and use 2 + 1 dimensional quantum electrodynamics as a benchmark example to demonstrate this efficient framework to achieve the continuum limit in LGTs. This possibility is a key requirement to make quantitative predictions at the field theory level and offers the long-term perspective to utilise quantum simulations to compute physically meaningful quantities in regimes that are precluded to quantum Monte Carlo.

arXiv:2006.14160•Quantum•Feb 4, 2021•129 citations
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 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) .
Taking the continuum limit is essential for extracting physical observables from quantum simulations of lattice gauge theories. Achieving the correct continuum limit requires careful control of all systematic uncertainties, including those arising from approximate implementations of the time evolution operator. In this work, we review existing approaches based on renormalization techniques, and point out their limitations. To overcome these limitations, we introduce a new general framework -- the Statistically-Bounded Time Evolution (SBTE) protocol -- for rigorously controlling the impact of approximate time evolution on the continuum limit. The central insight is that, since exact time evolution introduces no UV divergences, errors from approximate evolution can be treated as a source of systematic uncertainty that can be neglected if reduced below the working statistical uncertainty. We show that, using the SBTE protocol, which prescribes driving the approximate time evolution error below the working statistical uncertainty, leads to a simplified renormalization procedure. Furthermore, we show that, due to the existence of rigorous error bounds, one can guarantee a priori that such errors are negligible and do not affect the continuum limit. Ultimately, our protocol lays the foundation for performing systematic and fair comparisons between different simulation algorithms for lattice gauge theory simulations.
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 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.
Simulating the time evolution of quantum field theories given some Hamiltonian HHH requires developing algorithms for implementing the unitary operator ${\mathrm{e}}^{\ensuremath{-}iHt}$. A variety of techniques exist that accomplish this task, with the most common technique used so far being Trotterization, which is a special case of the application of a product formula. However, other techniques exist that promise better asymptotic scaling in certain parameters of the theory being simulated, the most efficient of which are based on the concept of block encoding. In this work we study the performance of such algorithms in simulating lattice field theories. We derive and compare the asymptotic gate complexities of several commonly used simulation techniques in application to Hamiltonian lattice field theories. Using the scalar ${\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\varphi}}}^{4}$ theory as a test, we also perform numerical studies and compare the gate costs required by product formulas and signal-processing-based techniques to simulate time evolution. For the latter, we use the linear combination of unitaries (LCU) construction augmented with the quantum Fourier transform circuit to switch between the field and momentum eigenbases, which leads to immediate order-of-magnitude improvement in the cost of preparing the block encoding. This paper also includes a pedagogical review of the techniques used, in particular product formulas, LCU, qubitization, quantum signal processing, as well as the technique for simulating geometrically-local Hamiltonians developed by Haah, Hastings, Kothari, and Low.
Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present limitations of MCMC techniques may be overcome by Hamiltonian-based simulations on classical or quantum devices, which further provide the potential to address questions that lay beyond the capabilities of the current approaches. However, for continuous gauge groups, Hamiltonian-based formulations involve infinite-dimensional gauge degrees of freedom that can solely be handled by truncation. Current truncation schemes require dramatically increasing computational resources at small values of the bare couplings, where magnetic field effects become important. Such limitation precludes one from `taking the continuous limit' while working with finite resources. To overcome this limitation, we provide a resource-efficient protocol to simulate LGTs with continuous gauge groups in the Hamiltonian formulation. Our new method allows for calculations at arbitrary values of the bare coupling and lattice spacing. The approach consists of the combination of a Hilbert space truncation with a regularization of the gauge group, which permits an efficient description of the magnetically-dominated regime. We focus here on Abelian gauge theories and use 2 + 1 dimensional quantum electrodynamics as a benchmark example to demonstrate this efficient framework to achieve the continuum limit in LGTs. This possibility is a key requirement to make quantitative predictions at the field theory level and offers the long-term perspective to utilise quantum simulations to compute physically meaningful quantities in regimes that are precluded to quantum Monte Carlo.