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:2608.16783 · Physical review. D/Physical review. D.

Quantum Simulation of QCD in Axial Gauge

Xiaojun Yao

hep-lathep-phhep-thnucl-thquant-phSubmitted Aug 17, 2026• Announced Aug 18, 2026• 0 citations
arXivHTMLPDF

Abstract

We study quantum simulation of SU(3) non-Abelian gauge theory dynamically coupled with fundamental fermions in $3+1$ dimensions by employing the lattice Hamiltonian in axial gauge that avoids Gauss's law constraints. The temporal component of the gauge field is analytically solved in terms of independent field degrees of freedom and a lattice regulated Green's function. The axial gauge condition is trivially maintained in time evolution, even under Trotterization. The gauge field degrees of freedom are expressed in the local field basis and can be efficiently transformed into the canonical conjugate momentum basis by local quantum Fourier transforms. We prove the number of qubits needed for describing all states up to an energy $E$ with an accuracy $ε$ on a lattice of volume $V$ at bare coupling $g$ is bounded as $16n_A V + 12n_f V$, where $n_A \approx \log_2 (\frac{64 E' V^{4/3}}{π^2ε} + \frac{32\sqrt{2}g n_f E'^{1/2} V^{7/6}}{\sqrt{3}π^3ε^{1/2}} ) $ is the number of qubits needed for each independent gauge field per site with a shifted energy $E'$, and $n_f$ denotes the number of fermion flavors. We then analyze a quantum algorithm for time evolution that is based on Trotterization, quantum Fourier transform, and Jordan-Wigner transformation, for which quantum circuits can be explicitly constructed under arbitrary gauge field truncation and digitization. We find the numbers of CNOT and single-qubit rotation gates both scale as $O(n_A^4 V^{4/3}) + O(V^{5/3})$ per Trotter step for fixed $n_f\leq 6$. We conclude that quantum resources needed for simulating real-time dynamics of lattice QCD scale polynomially with volume, energy, time, accuracy, and bare Hamiltonian parameters.

Comments
31 pages

Related papers

Quantum simulation of QED in Coulomb gauge

Xiaojun Yao

A recent work considered quantum simulation of Quantum Electrodynamics on a lattice in the Coulomb gauge with gauge degrees of freedom represented in the occupation basis in momentum space. Here we consider the more efficient representation of the gauge degrees of freedom in field basis in position space and develop a quantum algorithm for real-time simulation. We show that the continuum Coulomb gauge Hamiltonian is equivalent to the temporal gauge Hamiltonian when acting on physical states consisting of fermion and transverse gauge fields. The Coulomb gauge Hamiltonian is discretized by using the Green's function of the discrete Laplacian operator under the Dirichlet boundary conditions. Both the continuum Coulomb gauge Hamiltonian and the discretized one proposed here guarantee that the unphysical longitudinal gauge fields are decoupled and commute with the corresponding Hamiltonian. Thus there is no need to impose any constraint. The local gauge field basis and the canonically conjugate variable basis are swapped efficiently using the quantum Fourier transform. We prove that the qubit cost to represent physical states and the gate count for real-time simulation scale polynomially with the lattice size, energy, time, accuracy, and Hamiltonian parameters in lattice units. The gate cost here for implementing the time evolution of the gauge field is reduced at least by a factor on the order of $10^8$ for modest lattice size and accuracy level compared with the previous work.

arXiv:2507.01089•Physical review. D/Physical review. D.•Jan 27, 2026•2 citations

Preparations for quantum simulations of quantum chromodynamics in 1+1 dimensions. I. Axial gauge

Roland C. Farrell, Ivan A. Chernyshev, Sarah J. M. Powell +3 more

Tools necessary for quantum simulations of $1+1$ dimensional quantum chromodynamics are developed. When formulated in axial gauge and with two flavors of quarks, this system requires 12 qubits per spatial site with the gauge fields included via nonlocal interactions. Classical computations and $D$-wave's quantum annealer advantage are used to determine the hadronic spectrum, enabling a decomposition of the masses and a study of quark entanglement. Color ``edge'' states confined within a screening length of the end of the lattice are found. IBM's seven-qubit quantum computers, ibmq_jakarta and ibm_perth, are used to compute dynamics from the trivial vacuum in one-flavor QCD with one spatial site. More generally, the Hamiltonian and quantum circuits for time evolution of $1+1$ dimensional $SU({N}_{c})$ gauge theory with ${N}_{f}$ flavors of quarks are developed, and the resource requirements for large-scale quantum simulations are estimated.

arXiv:2207.01731•Physical review. D/Physical review. D.•Mar 30, 2023•95 citations

Efficient Quantum Simulations of Yang-Mills theory with Maximal-tree Gauge

Tianyin Li, Ying-Ying Li, Xiaoyang Wang +1 more

We develop a quantum algorithmic framework for the efficient simulation of Yang--Mills theories, including the $\mathrm{SU}(3)$ gauge theory in Quantum Chromodynamics (QCD). The framework uses maximal-tree gauge in terms of gauge field variables that removes all local gauge redundancies. In the resulting gauge-fixed formulation and digitization in the field-amplitude basis, we show that Hamiltonian time evolution admits an efficient implementation based on quantum singular value transformation (QSVT). We derive upper bounds on the total number of qubits and gate complexity, finding polynomial scaling with the inverse simulation precision $1/\varepsilon_s$, lattice volume $\mathcal{V}$, gauge coupling $g$, and target energy scale $E$. Our results provide a rigorous complexity-theoretic demonstration that non-Abelian Yang--Mills theories can be simulated efficiently on quantum computers, paving the way toward first-principles quantum simulations of non-perturbative QCD dynamics.

hep-latCross-listedarXiv:2608.27267•Aug 28, 2026•0 citations

Hamiltonian formulations of lattice gauge theories for quantum simulation

Indrakshi Raychowdhury

The development of quantum technology and the advent of quantum computation have renewed interest in Hamiltonian formulations of lattice gauge theories. In contrast to Euclidean Monte Carlo calculations, Hamiltonian simulation requires a direct representation of the physical Hilbert space, an efficient treatment of enforcing Gauss' law as constraints, and a formulation of real-time dynamics in terms of degrees of freedom that can be encoded on quantum or quantum-inspired classical platforms. The starting point corresponds to finding a suitable basis for Hamiltonian lattice gauge theories. The existing literature focuses on several complementary approaches: electric basis, magnetic basis, gauge-fixed, quantum-link, qubit-regularised, finite-group, and orbifold formulations. This article primarily focuses on the prepotential and loop--string--hadron formulations, where non-Abelian Gauss laws are solved locally and the dynamics is written directly in terms of gauge-invariant loop, string and hadron excitations. The discussion covers the SU(2) and SU(3) constructions in one spatial dimension, the role of residual Abelian Gauss laws, and the extension to higher dimensions through point splitting. We conclude by highlighting its direct applications: eliminating redundant gauge degrees of freedom makes exact diagonalisation and tensor network calculations feasible for studying real-time dynamics, provides insights into thermalisation properties, enables mapping of physical degrees of freedom onto analog and digital quantum hardware, and enables quantum simulation of dynamics.

International Symposium on Lattice Field Theory (LATTICE)•Jun 28, 2026•0 citations

Toward scalable simulations of lattice gauge theories on quantum computers

Simon V. Mathis, Guglielmo Mazzola, Ivano Tavernelli

The simulation of real-time dynamics in lattice gauge theories is particularly hard for classical computing due to the exponential scaling of the required resources. On the other hand, quantum algorithms can potentially perform the same calculation with a polynomial dependence on the number of degrees of freedom. A precise estimation is however particularly challenging for the simulation of lattice gauge theories in arbitrary dimensions, where gauge fields are dynamical variables, in addition to the particle fields. Moreover, there exist several choices for discretizing particles and gauge fields on a lattice, each of them coming at different prices in terms of qubit register size and circuit depth. Here we provide a resource counting for real-time evolution of $U(1)$ gauge theories, such as quantum electrodynamics, on arbitrary dimension using the Wilson fermion representation for the particles, and the quantum link model approach for the gauge fields. We study the phenomena of flux-string breaking up to a genuine bidimensional model using classical simulations of the quantum circuits and discuss the advantages of our discretization choice in simulation of more challenging $SU(N)$ gauge theories such as quantum chromodynamics.

arXiv:2005.10271•Physical review. D/Physical review. D.•Nov 4, 2020•68 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
We study quantum simulation of SU(3) non-Abelian gauge theory dynamically coupled with fundamental fermions in 3+13+13+1 dimensions by employing the lattice Hamiltonian in axial gauge that avoids Gauss's law constraints. The temporal component of the gauge field is analytically solved in terms of independent field degrees of freedom and a lattice regulated Green's function. The axial gauge condition is trivially maintained in time evolution, even under Trotterization. The gauge field degrees of freedom are expressed in the local field basis and can be efficiently transformed into the canonical conjugate momentum basis by local quantum Fourier transforms. We prove the number of qubits needed for describing all states up to an energy EEE with an accuracy εεε on a lattice of volume VVV at bare coupling ggg is bounded as 16nAV+12nfV16n_A V + 12n_f V16nA​V+12nf​V, where nA≈log⁡2(64E′V4/3π2ε+322gnfE′1/2V7/63π3ε1/2)n_A \approx \log_2 (\frac{64 E' V^{4/3}}{π^2ε} + \frac{32\sqrt{2}g n_f E'^{1/2} V^{7/6}}{\sqrt{3}π^3ε^{1/2}} ) nA​≈log2​(π2ε64E′V4/3​+3​π3ε1/2322​gnf​E′1/2V7/6​) is the number of qubits needed for each independent gauge field per site with a shifted energy E′E'E′, and nfn_fnf​ denotes the number of fermion flavors. We then analyze a quantum algorithm for time evolution that is based on Trotterization, quantum Fourier transform, and Jordan-Wigner transformation, for which quantum circuits can be explicitly constructed under arbitrary gauge field truncation and digitization. We find the numbers of CNOT and single-qubit rotation gates both scale as O(nA4V4/3)+O(V5/3)O(n_A^4 V^{4/3}) + O(V^{5/3})O(nA4​V4/3)+O(V5/3) per Trotter step for fixed nf≤6n_f\leq 6nf​≤6. We conclude that quantum resources needed for simulating real-time dynamics of lattice QCD scale polynomially with volume, energy, time, accuracy, and bare Hamiltonian parameters.
A recent work considered quantum simulation of Quantum Electrodynamics on a lattice in the Coulomb gauge with gauge degrees of freedom represented in the occupation basis in momentum space. Here we consider the more efficient representation of the gauge degrees of freedom in field basis in position space and develop a quantum algorithm for real-time simulation. We show that the continuum Coulomb gauge Hamiltonian is equivalent to the temporal gauge Hamiltonian when acting on physical states consisting of fermion and transverse gauge fields. The Coulomb gauge Hamiltonian is discretized by using the Green's function of the discrete Laplacian operator under the Dirichlet boundary conditions. Both the continuum Coulomb gauge Hamiltonian and the discretized one proposed here guarantee that the unphysical longitudinal gauge fields are decoupled and commute with the corresponding Hamiltonian. Thus there is no need to impose any constraint. The local gauge field basis and the canonically conjugate variable basis are swapped efficiently using the quantum Fourier transform. We prove that the qubit cost to represent physical states and the gate count for real-time simulation scale polynomially with the lattice size, energy, time, accuracy, and Hamiltonian parameters in lattice units. The gate cost here for implementing the time evolution of the gauge field is reduced at least by a factor on the order of 10810^8108 for modest lattice size and accuracy level compared with the previous work.
Tools necessary for quantum simulations of 1+11+11+1 dimensional quantum chromodynamics are developed. When formulated in axial gauge and with two flavors of quarks, this system requires 12 qubits per spatial site with the gauge fields included via nonlocal interactions. Classical computations and DDD-wave's quantum annealer advantage are used to determine the hadronic spectrum, enabling a decomposition of the masses and a study of quark entanglement. Color ``edge'' states confined within a screening length of the end of the lattice are found. IBM's seven-qubit quantum computers, ibmq_jakarta and ibm_perth, are used to compute dynamics from the trivial vacuum in one-flavor QCD with one spatial site. More generally, the Hamiltonian and quantum circuits for time evolution of 1+11+11+1 dimensional SU(Nc)SU({N}_{c})SU(Nc​) gauge theory with Nf{N}_{f}Nf​ flavors of quarks are developed, and the resource requirements for large-scale quantum simulations are estimated.
We develop a quantum algorithmic framework for the efficient simulation of Yang--Mills theories, including the SU(3)\mathrm{SU}(3)SU(3) gauge theory in Quantum Chromodynamics (QCD). The framework uses maximal-tree gauge in terms of gauge field variables that removes all local gauge redundancies. In the resulting gauge-fixed formulation and digitization in the field-amplitude basis, we show that Hamiltonian time evolution admits an efficient implementation based on quantum singular value transformation (QSVT). We derive upper bounds on the total number of qubits and gate complexity, finding polynomial scaling with the inverse simulation precision 1/εs1/\varepsilon_s1/εs​, lattice volume V\mathcal{V}V, gauge coupling ggg, and target energy scale EEE. Our results provide a rigorous complexity-theoretic demonstration that non-Abelian Yang--Mills theories can be simulated efficiently on quantum computers, paving the way toward first-principles quantum simulations of non-perturbative QCD dynamics.
The development of quantum technology and the advent of quantum computation have renewed interest in Hamiltonian formulations of lattice gauge theories. In contrast to Euclidean Monte Carlo calculations, Hamiltonian simulation requires a direct representation of the physical Hilbert space, an efficient treatment of enforcing Gauss' law as constraints, and a formulation of real-time dynamics in terms of degrees of freedom that can be encoded on quantum or quantum-inspired classical platforms. The starting point corresponds to finding a suitable basis for Hamiltonian lattice gauge theories. The existing literature focuses on several complementary approaches: electric basis, magnetic basis, gauge-fixed, quantum-link, qubit-regularised, finite-group, and orbifold formulations. This article primarily focuses on the prepotential and loop--string--hadron formulations, where non-Abelian Gauss laws are solved locally and the dynamics is written directly in terms of gauge-invariant loop, string and hadron excitations. The discussion covers the SU(2) and SU(3) constructions in one spatial dimension, the role of residual Abelian Gauss laws, and the extension to higher dimensions through point splitting. We conclude by highlighting its direct applications: eliminating redundant gauge degrees of freedom makes exact diagonalisation and tensor network calculations feasible for studying real-time dynamics, provides insights into thermalisation properties, enables mapping of physical degrees of freedom onto analog and digital quantum hardware, and enables quantum simulation of dynamics.
The simulation of real-time dynamics in lattice gauge theories is particularly hard for classical computing due to the exponential scaling of the required resources. On the other hand, quantum algorithms can potentially perform the same calculation with a polynomial dependence on the number of degrees of freedom. A precise estimation is however particularly challenging for the simulation of lattice gauge theories in arbitrary dimensions, where gauge fields are dynamical variables, in addition to the particle fields. Moreover, there exist several choices for discretizing particles and gauge fields on a lattice, each of them coming at different prices in terms of qubit register size and circuit depth. Here we provide a resource counting for real-time evolution of U(1)U(1)U(1) gauge theories, such as quantum electrodynamics, on arbitrary dimension using the Wilson fermion representation for the particles, and the quantum link model approach for the gauge fields. We study the phenomena of flux-string breaking up to a genuine bidimensional model using classical simulations of the quantum circuits and discuss the advantages of our discretization choice in simulation of more challenging SU(N)SU(N)SU(N) gauge theories such as quantum chromodynamics.
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.