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:2501.00579 · Nature Communications

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

Anton T. Than, Y. Y. Atas, A. Chakraborty, Jinglei Zhang, Matthew T. Diaz, Keru Wen, Xingxin Liu, Randy Lewis, Alaina Green, Christine A. Muschik, Norbert M. Linke

Published Nov 21, 2025• 12 citations
arXivHTMLPDF

Abstract

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms. Quantum simulations of the phase diagram of quantum chromodynamics faces hard challenges, such as having to prepare mixed states and enforcing the non-Abelian gauge symmetry constraints. Here, the authors show how to solve the two above problems in a trapped-ion device using motional ancillae and charge-singlet measurements.

DOI
10.1038/s41467-025-65198-w

Related papers

SU(2) hadrons on a quantum computer via a variational approach

Y. Y. Atas, Jinglei Zhang, Randy Lewis +3 more

Quantum computers have the potential to create important new opportunities for ongoing essential research on gauge theories. They can provide simulations that are unattainable on classical computers such as sign-problem afflicted models or time evolutions. In this work, we variationally prepare the low-lying eigenstates of a non-Abelian gauge theory with dynamically coupled matter on a quantum computer. This enables the observation of hadrons and the calculation of their associated masses. The SU(2) gauge group considered here represents an important first step towards ultimately studying quantum chromodynamics, the theory that describes the properties of protons, neutrons and other hadrons. Our calculations on an IBM superconducting platform utilize a variational quantum eigensolver to study both meson and baryon states, hadrons which have never been seen in a non-Abelian simulation on a quantum computer. We develop a hybrid resource-efficient approach by combining classical and quantum computing, that not only allows the study of an SU(2) gauge theory with dynamical matter fields on present-day quantum hardware, but further lays out the premises for future quantum simulations that will address currently unanswered questions in particle and nuclear physics.

arXiv:2102.08920•Nature Communications•Nov 11, 2021•134 citations

Variational thermal quantum simulation of the lattice Schwinger model

Xu-Dan Xie, Xingyu Guo, Hongxi Xing +3 more

Confinement of quarks due to the strong interaction and the deconfinement at high temperatures and high densities are a basic paradigm for understanding the nuclear matter. Their simulation, however, is very challenging for classical computers due to the sign problem of solving equilibrium states of finite-temperature quantum chromodynamical systems at finite density. In this paper, we propose a variational approach, using the lattice Schwinger model, to simulate the confinement or deconfinement by investigating the string tension. We adopt an ansatz that the string tension can be evaluated without referring to quantum protocols for measuring the entropy in the free energy. Results of numeral simulation show that the string tension decreases both along the increasing of the temperature and the chemical potential, which can be an analog of the phase diagram of QCD. Through numerical simulations on the classical computer, we demonstrate the potential of exploiting near-term quantum computers for investigating the phase diagram of finite-temperature and finite-density nuclear matters.

arXiv:2205.12767•Physical review. D/Physical review. D.•Sep 21, 2022•27 citations

Simulating one-dimensional quantum chromodynamics on a quantum computer: Real-time evolutions of tetra- and pentaquarks

Y. Y. Atas, Jan F. Haase, Jinglei Zhang +4 more

Quantum chromodynamics (QCD)---the theory of quarks and gluons---has been studied for decades, but it is yet to be fully understood. A recent example is the prediction and experimental discovery of tetraquarks, which opened a new research field. Crucially, numerous unsolved questions regarding the standard model can exclusively be addressed by nonperturbative calculations. Quantum computers can solve problems for which well-established QCD methods are inapplicable, such as real-time evolution. We take a key step in exploring this possibility by designing a real-time evolution of tetraquark and pentaquark physics in one-dimensional SU(3) gauge theory. We also perform an experiment on a superconducting quantum computer demonstrating an elementary cell consisting of two staggered sites containing quarks and antiquarks with all three colors. This experiment represents an exciting step in quantum computation involving quarks with the gauge group of QCD.

arXiv:2207.03473•Physical Review Research•Sep 13, 2023•58 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

Toward Exploring Phase Diagrams of Gauge Theories on Quantum Computers with Thermal Pure Quantum States

Connor Powers, Zohreh Davoudi, Niklas Mueller

Aiming at evading the notorious sign problem in classical Monte-Carlo approaches to lattice quantum chromodynamics, we present an approach for quantum computing finite-temperature lattice gauge theories at non-zero density. Based on the thermal pure-quantum-state formalism of statistical mechanics when extended to gauge-theory systems, our approach allows for sign-problem-free quantum computations of thermal expectation values and non-equal time correlation functions. By taking a simple lattice gauge theory for which classical benchmarks are possible, namely $\mathbb{Z}_2$ lattice gauge theory in 1+1 dimensions at finite chemical potential, we discuss resource requirements and robustness to algorithmic and hardware imperfections for near-term quantum-hardware realizations.

arXiv:2212.11388•International Symposium on Lattice Field Theory (LATTICE)•Jan 29, 2023•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
The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms. Quantum simulations of the phase diagram of quantum chromodynamics faces hard challenges, such as having to prepare mixed states and enforcing the non-Abelian gauge symmetry constraints. Here, the authors show how to solve the two above problems in a trapped-ion device using motional ancillae and charge-singlet measurements.
Quantum computers have the potential to create important new opportunities for ongoing essential research on gauge theories. They can provide simulations that are unattainable on classical computers such as sign-problem afflicted models or time evolutions. In this work, we variationally prepare the low-lying eigenstates of a non-Abelian gauge theory with dynamically coupled matter on a quantum computer. This enables the observation of hadrons and the calculation of their associated masses. The SU(2) gauge group considered here represents an important first step towards ultimately studying quantum chromodynamics, the theory that describes the properties of protons, neutrons and other hadrons. Our calculations on an IBM superconducting platform utilize a variational quantum eigensolver to study both meson and baryon states, hadrons which have never been seen in a non-Abelian simulation on a quantum computer. We develop a hybrid resource-efficient approach by combining classical and quantum computing, that not only allows the study of an SU(2) gauge theory with dynamical matter fields on present-day quantum hardware, but further lays out the premises for future quantum simulations that will address currently unanswered questions in particle and nuclear physics.
Confinement of quarks due to the strong interaction and the deconfinement at high temperatures and high densities are a basic paradigm for understanding the nuclear matter. Their simulation, however, is very challenging for classical computers due to the sign problem of solving equilibrium states of finite-temperature quantum chromodynamical systems at finite density. In this paper, we propose a variational approach, using the lattice Schwinger model, to simulate the confinement or deconfinement by investigating the string tension. We adopt an ansatz that the string tension can be evaluated without referring to quantum protocols for measuring the entropy in the free energy. Results of numeral simulation show that the string tension decreases both along the increasing of the temperature and the chemical potential, which can be an analog of the phase diagram of QCD. Through numerical simulations on the classical computer, we demonstrate the potential of exploiting near-term quantum computers for investigating the phase diagram of finite-temperature and finite-density nuclear matters.
Quantum chromodynamics (QCD)---the theory of quarks and gluons---has been studied for decades, but it is yet to be fully understood. A recent example is the prediction and experimental discovery of tetraquarks, which opened a new research field. Crucially, numerous unsolved questions regarding the standard model can exclusively be addressed by nonperturbative calculations. Quantum computers can solve problems for which well-established QCD methods are inapplicable, such as real-time evolution. We take a key step in exploring this possibility by designing a real-time evolution of tetraquark and pentaquark physics in one-dimensional SU(3) gauge theory. We also perform an experiment on a superconducting quantum computer demonstrating an elementary cell consisting of two staggered sites containing quarks and antiquarks with all three colors. This experiment represents an exciting step in quantum computation involving quarks with the gauge group of QCD.
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.
Aiming at evading the notorious sign problem in classical Monte-Carlo approaches to lattice quantum chromodynamics, we present an approach for quantum computing finite-temperature lattice gauge theories at non-zero density. Based on the thermal pure-quantum-state formalism of statistical mechanics when extended to gauge-theory systems, our approach allows for sign-problem-free quantum computations of thermal expectation values and non-equal time correlation functions. By taking a simple lattice gauge theory for which classical benchmarks are possible, namely Z2\mathbb{Z}_2Z2​ lattice gauge theory in 1+1 dimensions at finite chemical potential, we discuss resource requirements and robustness to algorithmic and hardware imperfections for near-term quantum-hardware realizations.
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.