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:2602.18080

Observation of Robust and Coherent Non-Abelian Hadron Dynamics on Noisy Quantum Processors

Fran Ilčić, Ritajit Majumdar, Emil Mathew, Md. Osama Ali, Nathan Earnest-Noble, Indrakshi Raychowdhury

hep-latcond-mat.str-elhep-thnucl-thquant-phSubmitted Feb 20, 2026• Announced Feb 23, 2026• 0 citations• Revised Mar 11, 2026 (v2)
v1 · Feb 20, 2026v2 · Mar 11, 2026v3 · May 29, 2026
arXivHTMLPDF

Abstract

The real-time evolution of strongly interacting matter remains a frontier of fundamental physics, as classical simulations are hampered by exponential Hilbert space growth and rapid, unmanageable growth of quantum entanglement. This study reports the quantum simulation of hadron dynamics within a (1 + 1)-dimensional SU(2) lattice gauge theory using a 156-qubit IBM superconducting processor. Leveraging a hardware-efficient Loop-String-Hadron (LSH) encoding, we simulate the dynamics of the physical degrees of freedom on a 60-site lattice in the weak-coupling regime, as a crucial step toward the continuum limit. We successfully observe the light-cone propagation of a confined meson and internal oscillations indicative of early-time hadronic breathing modes. Notably, these high-fidelity results were obtained directly from the quantum data via a differential measurement protocol, together with measurement error mitigation, demonstrating a robust pathway for large-scale simulations even on noisy hardware. To validate the results, we benchmarked the quantum algorithm and its outcomes from the quantum processor against state-of-the-art approximated classical algorithms using tensor network methods and Pauli propagation, respectively. Furthermore, we provide a quantitative comparison demonstrating that as the system approaches the weak-coupling or the continuum limit, the quantum processor maintains a consistent structural robustness where classical tensor networks and Pauli propagation methods encounter an onset of exponential complexity or symmetry violations as an artifact of approximation in the algorithm. These results establish a scalable pathway for simulating non-Abelian dynamics on near-term quantum hardware and mark a critical step toward achieving a practical quantum advantage in high-energy physics.

Comments
21 pages, 13 figures including supplementary information and figures. v2 contains one additional classical benchmark and error analysis, and a few changes in figures. The main conclusion remains unchanged

Related papers

Cold-atom quantum simulator for string and hadron dynamics in non-Abelian lattice gauge theory

Raka Dasgupta, Indrakshi Raychowdhury

We propose an analog quantum simulator for simulating real-time dynamics of $(1+1)$-dimensional non-Abelian gauge theory well within the existing capacity of ultracold-atom experiments. The scheme calls for the realization of a two-state ultracold fermionic system in a one-dimensional bipartite lattice, and the observation of subsequent tunneling dynamics. Being based on the loop string hadron formalism of $SU(2)$ lattice gauge theory, this simulation technique is completely $SU(2)$ invariant and simulates accurate dynamics of physical phenomena such as string breaking and/or pair production. The scheme is scalable and particularly effective in simulating the theory in the weak-coupling regime, and also a bulk limit of the theory in the strong-coupling regime up to certain approximations. This paper also presents a numerical benchmark comparison of the exact spectrum and real-time dynamics of lattice gauge theory to that of the atomic Hamiltonian with an experimentally realizable range of parameters.

arXiv:2009.13969•Physical Review A•Feb 22, 2022•37 citations

Dynamics of loops for SU(2) gauge theory in (2+1)-dimension

Aahiri Naskar

Pure gauge theories in more than one spatial dimension are extremely interesting as they allow pure gauge observables such as glueballs and Wilson loops. Simulating the real-time dynamics of pure non-Abelian lattice gauge theories poses a significant computational challenge due to the large number of spurious gauge degrees of freedom. This work presents a brute-force (exact diagonalization) study of 2+1-dimensional SU(2) lattice gauge theory, leveraging the gauge-invariant Loop-String-Hadron (LSH) framework. Using only CPU-based computation, we have successfully demonstrated classical simulation on systems of up to 9 plaquettes and aim to push it further by leveraging combinations of CPUs and GPUs. This work provides crucial benchmarks for GPU implementations, near-term tensor network simulators, and quantum simulations, opening a new avenue for studying fundamental phenomena like thermalization in a non-Abelian gauge theory directly from the Hamiltonian perspective in higher dimensions.

International Symposium on Lattice Field Theory (LATTICE)•Oct 11, 2026•0 citations

Digital Quantum Simulation of Dynamics in Interacting Fields

Zhiyao Li

Matter under extreme conditions behaves as a non-equilibrium, strongly-interacting quantum system. Such systems are in the non-perturbative regime of quantum chromodynamics (QCD), where their time evolution is classically intractable. Quantum computing offers a promising alternative, but present-day devices are limited in qubit count and circuit depth, so efficient formulations and algorithms are key to scaling toward experimentally relevant QCD simulations. A scalable framework is established for simulating energy loss and hadronization in 1+1-dimensional SU(2) lattice gauge theory, in which a heavy quark traverses a lattice of light quarks. Classical lattice QCD methods are translated into physics-informed circuits to prepare the interacting ground state. The evolution is run on IBM's quantum processor ibm_pittsburgh, where the charge distribution and energy loss are measured. A further study develops the quantum simulation of the fully gauge-fixed SU(2) Hamiltonian in a mixed basis, a formulation that is expected to remain efficient at all gauge couplings. Digitizing the continuous gauge fields with as few as three qubits per plaquette reproduces low-energy observables with $10^{-3}$ precision. Scaling to larger systems is limited by resource overhead, which motivates truncations that preserve the relevant observables. Sequency Hierarchy Truncation (SeqHT) exploits the hierarchy in basis-operator contributions. Demonstrated on the adiabatic preparation of the $λφ^4$ ground state on ibm_sherbrooke, SeqHT reduces circuit depth by ~30%, suppressing hardware noise in the extracted observables. All simulations in this dissertation are validated on quantum hardware against classical benchmarks. These studies advance quantum simulations of real-time QCD dynamics.

quant-pharXiv:2610.28993•Oct 22, 2026•0 citations

Tailoring Non-Abelian Gauge Theory for Digital Quantum Simulation

Jesse R. Stryker, Indrakshi Raychowdhury

The question of how to efficiently formulate Hamiltonian gauge theories is experiencing renewed interest due to advances in building quantum simulation platforms. We introduce a reformulation of an SU(2) Hamiltonian lattice gauge theory—a loop-string-hadron (LSH) formulation—in which the dynamical degrees of freedom are localized pieces of flux loops, meson strings, and hadrons. LSH operators are first derived from Schwinger bosons and used to construct a Hilbert space with the non-Abelian Gauss law built into it. They are subsequently factored into products of “normalized” ladder operators and diagonal matrices, priming them for classical or quantum information processing. The LSH formalism alleviates several disadvantages of quantum simulating the Kogut-Susskind formulation and makes little use of structures specific to SU(2). Its conceptual clarity makes it an attractive approach to apply to other non-Abelian groups like SU(3).

International Symposium on Lattice Field Theory (LATTICE)•May 27, 2020•1 citation

Trapped-ion quantum simulation of (1+1)D SU(2) lattice gauge theory via loop-string-hadron formulation

Toshiaki Kaji, Jesse R. Stryker, Yutaro Iiyama +4 more

We present a quantum simulation of the SU(2) lattice gauge theory in 1+1 dimensions with one flavor of fermions, performed using the trapped-ion Quantinuum System Model H2. This work is based on the loop-string-hadron (LSH) formulation of the gauge theory, and is distinguished by using qubits to represent both bosonic-gauge and fermionic-matter degrees of freedom. Furthermore, we access the $j=1$ electric-flux sector beyond the $j_{\max}=1/2$ truncation commonly used in SU(2) quantum simulations. To that end, we digitize an electric-field degree of freedom using a multiqubit register for representing the bosonic quantum number. Our simulation therefore synthesizes several key ingredients relevant to quantum simulations in higher spatial dimensions, in particular explicit local gauge-field degrees of freedom represented here by a multiqubit bosonic register, while relying solely on circuit-construction optimizations that are expected to scale polynomially with lattice size. This work demonstrates that modern quantum hardware is increasingly able to accommodate the features of complicated but physically important lattice gauge theories while successfully reproducing the expected Trotterized dynamics.

hep-latCross-listedarXiv:2610.24091•Oct 20, 2026•0 citations

A Framework for Quantum Simulations of Energy-Loss and Hadronization in Non-Abelian Gauge Theories: SU(2) Lattice Gauge Theory in 1+1D

Zhiyao Li, Marc Illa, Martin J. Savage

Simulations of energy loss and hadronization are essential for understanding a range of phenomena in non-equilibrium strongly-interacting matter. We establish a framework for performing such simulations on a quantum computer and apply it to a heavy quark moving across a modest-sized 1+1D SU(2) lattice of light quarks. Conceptual advances with regard to simulations of non-Abelian versus Abelian theories are developed, allowing for the evolution of the energy in light quarks, of their local non-Abelian charge densities, and of their multi-partite entanglement to be computed. The non-trivial action of non-Abelian charge operators on arbitrary states suggests mapping the heavy quarks to qubits alongside the light quarks, and limits the heavy-quark motion to discrete steps among spatial lattice sites. Further, the color entanglement among the heavy quarks and light quarks is implemented using hadronic operators, and Domain Decomposition is shown to be effective in quantum state preparation. Scalable quantum circuits that account for the heterogeneity of non-Abelian charge sectors across the lattice are used to prepare the interacting ground-state wavefunction in the presence of heavy quarks. The discrete motion of heavy quarks between adjacent spatial sites is implemented using fermionic SWAP operations. Quantum simulations of the dynamics of a system on $L=3$ spatial sites are performed using IBM's ${\tt ibm\_pittsburgh}$ quantum computer using 18 qubits, for which the circuits for state preparation, motion, and one second-order Trotter step of time evolution have a two-qubit depth of 398. A suite of error mitigation techniques are used to extract the observables from the simulations, providing results that are in good agreement with classical simulations. The framework presented here generalizes straightforwardly to other non-Abelian groups, including SU(3) for quantum chromodynamics.

quant-pharXiv:2512.05210•Dec 8, 2025•0 citations
The real-time evolution of strongly interacting matter remains a frontier of fundamental physics, as classical simulations are hampered by exponential Hilbert space growth and rapid, unmanageable growth of quantum entanglement. This study reports the quantum simulation of hadron dynamics within a (1 + 1)-dimensional SU(2) lattice gauge theory using a 156-qubit IBM superconducting processor. Leveraging a hardware-efficient Loop-String-Hadron (LSH) encoding, we simulate the dynamics of the physical degrees of freedom on a 60-site lattice in the weak-coupling regime, as a crucial step toward the continuum limit. We successfully observe the light-cone propagation of a confined meson and internal oscillations indicative of early-time hadronic breathing modes. Notably, these high-fidelity results were obtained directly from the quantum data via a differential measurement protocol, together with measurement error mitigation, demonstrating a robust pathway for large-scale simulations even on noisy hardware. To validate the results, we benchmarked the quantum algorithm and its outcomes from the quantum processor against state-of-the-art approximated classical algorithms using tensor network methods and Pauli propagation, respectively. Furthermore, we provide a quantitative comparison demonstrating that as the system approaches the weak-coupling or the continuum limit, the quantum processor maintains a consistent structural robustness where classical tensor networks and Pauli propagation methods encounter an onset of exponential complexity or symmetry violations as an artifact of approximation in the algorithm. These results establish a scalable pathway for simulating non-Abelian dynamics on near-term quantum hardware and mark a critical step toward achieving a practical quantum advantage in high-energy physics.
We propose an analog quantum simulator for simulating real-time dynamics of (1+1)(1+1)(1+1)-dimensional non-Abelian gauge theory well within the existing capacity of ultracold-atom experiments. The scheme calls for the realization of a two-state ultracold fermionic system in a one-dimensional bipartite lattice, and the observation of subsequent tunneling dynamics. Being based on the loop string hadron formalism of SU(2)SU(2)SU(2) lattice gauge theory, this simulation technique is completely SU(2)SU(2)SU(2) invariant and simulates accurate dynamics of physical phenomena such as string breaking and/or pair production. The scheme is scalable and particularly effective in simulating the theory in the weak-coupling regime, and also a bulk limit of the theory in the strong-coupling regime up to certain approximations. This paper also presents a numerical benchmark comparison of the exact spectrum and real-time dynamics of lattice gauge theory to that of the atomic Hamiltonian with an experimentally realizable range of parameters.
Pure gauge theories in more than one spatial dimension are extremely interesting as they allow pure gauge observables such as glueballs and Wilson loops. Simulating the real-time dynamics of pure non-Abelian lattice gauge theories poses a significant computational challenge due to the large number of spurious gauge degrees of freedom. This work presents a brute-force (exact diagonalization) study of 2+1-dimensional SU(2) lattice gauge theory, leveraging the gauge-invariant Loop-String-Hadron (LSH) framework. Using only CPU-based computation, we have successfully demonstrated classical simulation on systems of up to 9 plaquettes and aim to push it further by leveraging combinations of CPUs and GPUs. This work provides crucial benchmarks for GPU implementations, near-term tensor network simulators, and quantum simulations, opening a new avenue for studying fundamental phenomena like thermalization in a non-Abelian gauge theory directly from the Hamiltonian perspective in higher dimensions.
Matter under extreme conditions behaves as a non-equilibrium, strongly-interacting quantum system. Such systems are in the non-perturbative regime of quantum chromodynamics (QCD), where their time evolution is classically intractable. Quantum computing offers a promising alternative, but present-day devices are limited in qubit count and circuit depth, so efficient formulations and algorithms are key to scaling toward experimentally relevant QCD simulations. A scalable framework is established for simulating energy loss and hadronization in 1+1-dimensional SU(2) lattice gauge theory, in which a heavy quark traverses a lattice of light quarks. Classical lattice QCD methods are translated into physics-informed circuits to prepare the interacting ground state. The evolution is run on IBM's quantum processor ibm_pittsburgh, where the charge distribution and energy loss are measured. A further study develops the quantum simulation of the fully gauge-fixed SU(2) Hamiltonian in a mixed basis, a formulation that is expected to remain efficient at all gauge couplings. Digitizing the continuous gauge fields with as few as three qubits per plaquette reproduces low-energy observables with 10−310^{-3}10−3 precision. Scaling to larger systems is limited by resource overhead, which motivates truncations that preserve the relevant observables. Sequency Hierarchy Truncation (SeqHT) exploits the hierarchy in basis-operator contributions. Demonstrated on the adiabatic preparation of the λφ4λφ^4λφ4 ground state on ibm_sherbrooke, SeqHT reduces circuit depth by ~30%, suppressing hardware noise in the extracted observables. All simulations in this dissertation are validated on quantum hardware against classical benchmarks. These studies advance quantum simulations of real-time QCD dynamics.
The question of how to efficiently formulate Hamiltonian gauge theories is experiencing renewed interest due to advances in building quantum simulation platforms. We introduce a reformulation of an SU(2) Hamiltonian lattice gauge theory—a loop-string-hadron (LSH) formulation—in which the dynamical degrees of freedom are localized pieces of flux loops, meson strings, and hadrons. LSH operators are first derived from Schwinger bosons and used to construct a Hilbert space with the non-Abelian Gauss law built into it. They are subsequently factored into products of “normalized” ladder operators and diagonal matrices, priming them for classical or quantum information processing. The LSH formalism alleviates several disadvantages of quantum simulating the Kogut-Susskind formulation and makes little use of structures specific to SU(2). Its conceptual clarity makes it an attractive approach to apply to other non-Abelian groups like SU(3).
We present a quantum simulation of the SU(2) lattice gauge theory in 1+1 dimensions with one flavor of fermions, performed using the trapped-ion Quantinuum System Model H2. This work is based on the loop-string-hadron (LSH) formulation of the gauge theory, and is distinguished by using qubits to represent both bosonic-gauge and fermionic-matter degrees of freedom. Furthermore, we access the j=1j=1j=1 electric-flux sector beyond the jmax⁡=1/2j_{\max}=1/2jmax​=1/2 truncation commonly used in SU(2) quantum simulations. To that end, we digitize an electric-field degree of freedom using a multiqubit register for representing the bosonic quantum number. Our simulation therefore synthesizes several key ingredients relevant to quantum simulations in higher spatial dimensions, in particular explicit local gauge-field degrees of freedom represented here by a multiqubit bosonic register, while relying solely on circuit-construction optimizations that are expected to scale polynomially with lattice size. This work demonstrates that modern quantum hardware is increasingly able to accommodate the features of complicated but physically important lattice gauge theories while successfully reproducing the expected Trotterized dynamics.
Simulations of energy loss and hadronization are essential for understanding a range of phenomena in non-equilibrium strongly-interacting matter. We establish a framework for performing such simulations on a quantum computer and apply it to a heavy quark moving across a modest-sized 1+1D SU(2) lattice of light quarks. Conceptual advances with regard to simulations of non-Abelian versus Abelian theories are developed, allowing for the evolution of the energy in light quarks, of their local non-Abelian charge densities, and of their multi-partite entanglement to be computed. The non-trivial action of non-Abelian charge operators on arbitrary states suggests mapping the heavy quarks to qubits alongside the light quarks, and limits the heavy-quark motion to discrete steps among spatial lattice sites. Further, the color entanglement among the heavy quarks and light quarks is implemented using hadronic operators, and Domain Decomposition is shown to be effective in quantum state preparation. Scalable quantum circuits that account for the heterogeneity of non-Abelian charge sectors across the lattice are used to prepare the interacting ground-state wavefunction in the presence of heavy quarks. The discrete motion of heavy quarks between adjacent spatial sites is implemented using fermionic SWAP operations. Quantum simulations of the dynamics of a system on L=3L=3L=3 spatial sites are performed using IBM's ibm_pittsburgh{\tt ibm\_pittsburgh}ibm_pittsburgh quantum computer using 18 qubits, for which the circuits for state preparation, motion, and one second-order Trotter step of time evolution have a two-qubit depth of 398. A suite of error mitigation techniques are used to extract the observables from the simulations, providing results that are in good agreement with classical simulations. The framework presented here generalizes straightforwardly to other non-Abelian groups, including SU(3) for quantum chromodynamics.