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

Quantum simulation of massive Thirring and Gross--Neveu models for arbitrary number of flavors

Bojko N. Bakalov, Joao C. Getelina, Raghav G. Jha, Alexander F. Kemper, Yuan Liu

quant-phhep-lathep-thSubmitted Feb 25, 2026• Announced Feb 27, 2026• 0 citations• Revised May 9, 2026 (v2)
v1 · Feb 25, 2026v2 · May 9, 2026
arXivHTMLPDF

Abstract

The study of fermionic quantum field theories is an important problem for realizing the standard model of particle physics on a quantum computer. As a step towards this goal, we consider the massive Thirring and Gross--Neveu models with arbitrary number of fermion flavors, $N_f$, discretized on a spatial one-dimensional lattice of size $L$ in the Hamiltonian formulation. We compute the gate complexity using the higher-order product formula and using block-encoding/qubitization and quantum singular value transformations in the limit of large $N_f$ and $L$. We also prepare the ground states of both models with excellent fidelity for system sizes up to 20 qubits with $N_f = 1,2,3,4$ using the adaptive-variational quantum imaginary time algorithm. In addition, we also classify the dynamical Lie algebras of these relativistic fermionic models and show that they belong to the same isomorphism class. Our work is a concrete step towards the quantum simulation of real-time dynamics of large $N_f$ fermionic quantum field theories models relevant for chiral symmetry breaking, understanding dimensional transmutation, and exploring the conformal window of field theories on near-term and early fault-tolerant quantum computers.

Journal reference
Phys. Rev. D 113, 114508 (2026)
DOI
10.1103/h7mq-cbz4
Comments
v2: Version consistent with one accepted for publication in Phys. Rev. D. v1: 1+29 pages, 4+1 figures

Related papers

Quantum Simulation of the Real-time Dynamics in the multi-flavor Gross-Neveu Model at the utility scale using Superconducting Quantum Computers

Talal Ahmed Chowdhury, Seokwon Choi, Kyoungchul Kong +1 more

We present a scalable quantum simulation framework for real-time dynamics of the multi-flavor Gross-Neveu model in 1+1 dimensions. Using superconducting quantum processors at utility scale, we develop a hardware-efficient Trotterization whose per-step circuit depth scales with fermion flavor number rather than total system size, enabling simulations beyond 100 qubits. A central contribution of this work is the Localized Diagonal Operator Approximation (LDOA), which systematically reduces the overhead associated with quartic interactions. We formulate diagonal unitary synthesis as a structured least-squares problem in phase space and obtain analytic solutions via the Moore-Penrose pseudoinverse. This formulation provides a principled and quantitatively controlled approximation: in the small Trotter-step regime, the unitary error is directly linked to the phase reconstruction error and vanishes asymptotically as the Trotter step size decreases. This establishes a clear mathematical foundation for the LDOA while significantly reducing two-qubit gate counts and circuit depth, and is broadly applicable to diagonal quantum operators with long-range structure, making it particularly well suited for quantum hardware with limited qubit connectivity. Using these techniques, we run large-scale simulations on IBM superconducting processors and study real-time observables, including density-density correlators. We benchmark against exact diagonalization and tensor network-based methods, finding strong agreement across system sizes. These results show that combining hardware-aware circuit design with rigorous approximations enables practical near-term simulation of interacting fermionic field theories and provides a scalable pathway toward more complex quantum field theory simulations.

quant-pharXiv:2605.05479•May 8, 2026•0 citations

Faster quantum algorithm to simulate fermionic quantum field theory

Ali Hamed Moosavian, Stephen P. Jordan

In quantum algorithms discovered so far for simulating scattering processes in quantum field theories, state preparation is the slowest step. We present an algorithm for preparing particle states to use in simulation of fermionic quantum field theory (QFT) on a quantum computer, which is based on the matrix product state ansatz. We apply this to the massive Gross-Neveu model in one spatial dimension to illustrate the algorithm, but we believe the same algorithm with slight modifications can be used to simulate any one-dimensional massive fermionic QFT. In the case where the number of particle species is 1, our algorithm can prepare particle states using $O\left({\ensuremath{\epsilon}}^{\ensuremath{-}3.23...}\right)$ gates, which is much faster than previous known results, namely, $O\left({\ensuremath{\epsilon}}^{\ensuremath{-}8\ensuremath{-}o\left(1\right)}\right)$. Furthermore, unlike previous methods which were based on adiabatic state preparation, the method given here should be able to simulate quantum phases unconnected to the free theory.

arXiv:1711.04006•Physical Review A•Jul 30, 2018•40 citations

Quantum Algorithms for Fermionic Quantum Field Theories

Stephen P. Jordan, Keith S. M. Lee, John P. Preskill

Extending previous work on scalar field theories, we develop a quantum algorithm to compute relativistic scattering amplitudes in fermionic field theories, exemplified by the massive Gross-Neveu model, a theory in two spacetime dimensions with quartic interactions. The algorithm introduces new techniques to meet the additional challenges posed by the characteristics of fermionic fields, and its run time is polynomial in the desired precision and the energy. Thus, it constitutes further progress towards an efficient quantum algorithm for simulating the Standard Model of particle physics.

arXiv:1404.7115•Apr 28, 2014•52 citations

Gross–Neveu–Wilson model and correlated symmetry-protected topological phases

A. Bermúdez, Emanuele Tirrito, Matteo Rizzi +2 more

We show that a Wilson-type discretization of the Gross–Neveu model, a fermionic N-flavor quantum field theory displaying asymptotic freedom and chiral symmetry breaking, can serve as a playground to explore correlated symmetry-protected phases of matter using techniques borrowed from high-energy physics. A large-N study, both in the Hamiltonian and Euclidean formalisms, yields a phase diagram with trivial, topological, and symmetry-broken phases separated by critical lines that meet at a tri-critical point. We benchmark these predictions using tools from condensed matter and quantum information science, which show that the large-N method captures the essence of the phase diagram even at N=1. Moreover, we describe a cold-atom scheme for the quantum simulation of this lattice model, which would allow to explore the single-flavor phase diagram.

arXiv:1807.03202•Annals of Physics•Oct 18, 2018•46 citations

Quantum simulation of the N -flavor Gross-Neveu model

Muhammad Asaduzzaman, Goksu Can Toga, Simon M. Catterall +2 more

We discuss the use of quantum simulation to study an $N$-flavor theory of interacting relativistic fermions in ($1+1$) dimensions on noisy intermediate-scale quantum (NISQ)-era machines. The case of two flavors is particularly interesting as it can be mapped to the Hubbard model. We derive the appropriate qubit Hamiltonians and associated quantum circuits. We compare classical simulation and density matrix renormalization group/time-evolving blocked decimation calculations with the results of quantum simulation on various platforms for $N=2$ and 4-flavors. We demonstrate that the four steps needed for calculations of real-time scattering can be implemented using current NISQ devices.

arXiv:2208.05906•Physical review. D/Physical review. D.•Dec 30, 2022•14 citations

Simulation of state evolutions in Gross-Neveu model by matrix product state representation

Desheng Li, Hao Wang, Chu Yu Guo +2 more

A quantum algorithm to simulate the real time dynamics of two-flavor massive Gross-Neveu model is presented in Schrodinger picture. We implement the simulation on a classic computer by applying the matrix product state representation. The real time evolutions of up to four particles on a site in initial state are figured out in space-time coordinate. The state evolutions are effectively affected by fermion mass and coupling constant of the model. Especially when the mass of fermion is small enough and the coupling is strong enough, the fundamental fermions evolve synchronistically in space from the two-fermion and four-fermion initial states. These are also the conditions on which the bound states made up of fundamental fermion pairs were found to arise automatically in the literatures.

arXiv:2011.07538•Nov 15, 2020•1 citation
The study of fermionic quantum field theories is an important problem for realizing the standard model of particle physics on a quantum computer. As a step towards this goal, we consider the massive Thirring and Gross--Neveu models with arbitrary number of fermion flavors, NfN_fNf​, discretized on a spatial one-dimensional lattice of size LLL in the Hamiltonian formulation. We compute the gate complexity using the higher-order product formula and using block-encoding/qubitization and quantum singular value transformations in the limit of large NfN_fNf​ and LLL. We also prepare the ground states of both models with excellent fidelity for system sizes up to 20 qubits with Nf=1,2,3,4N_f = 1,2,3,4Nf​=1,2,3,4 using the adaptive-variational quantum imaginary time algorithm. In addition, we also classify the dynamical Lie algebras of these relativistic fermionic models and show that they belong to the same isomorphism class. Our work is a concrete step towards the quantum simulation of real-time dynamics of large NfN_fNf​ fermionic quantum field theories models relevant for chiral symmetry breaking, understanding dimensional transmutation, and exploring the conformal window of field theories on near-term and early fault-tolerant quantum computers.
We present a scalable quantum simulation framework for real-time dynamics of the multi-flavor Gross-Neveu model in 1+1 dimensions. Using superconducting quantum processors at utility scale, we develop a hardware-efficient Trotterization whose per-step circuit depth scales with fermion flavor number rather than total system size, enabling simulations beyond 100 qubits. A central contribution of this work is the Localized Diagonal Operator Approximation (LDOA), which systematically reduces the overhead associated with quartic interactions. We formulate diagonal unitary synthesis as a structured least-squares problem in phase space and obtain analytic solutions via the Moore-Penrose pseudoinverse. This formulation provides a principled and quantitatively controlled approximation: in the small Trotter-step regime, the unitary error is directly linked to the phase reconstruction error and vanishes asymptotically as the Trotter step size decreases. This establishes a clear mathematical foundation for the LDOA while significantly reducing two-qubit gate counts and circuit depth, and is broadly applicable to diagonal quantum operators with long-range structure, making it particularly well suited for quantum hardware with limited qubit connectivity. Using these techniques, we run large-scale simulations on IBM superconducting processors and study real-time observables, including density-density correlators. We benchmark against exact diagonalization and tensor network-based methods, finding strong agreement across system sizes. These results show that combining hardware-aware circuit design with rigorous approximations enables practical near-term simulation of interacting fermionic field theories and provides a scalable pathway toward more complex quantum field theory simulations.
In quantum algorithms discovered so far for simulating scattering processes in quantum field theories, state preparation is the slowest step. We present an algorithm for preparing particle states to use in simulation of fermionic quantum field theory (QFT) on a quantum computer, which is based on the matrix product state ansatz. We apply this to the massive Gross-Neveu model in one spatial dimension to illustrate the algorithm, but we believe the same algorithm with slight modifications can be used to simulate any one-dimensional massive fermionic QFT. In the case where the number of particle species is 1, our algorithm can prepare particle states using $O\left({\ensuremath{\epsilon}}^{\ensuremath{-}3.23...}\right)$ gates, which is much faster than previous known results, namely, $O\left({\ensuremath{\epsilon}}^{\ensuremath{-}8\ensuremath{-}o\left(1\right)}\right)$. Furthermore, unlike previous methods which were based on adiabatic state preparation, the method given here should be able to simulate quantum phases unconnected to the free theory.
Extending previous work on scalar field theories, we develop a quantum algorithm to compute relativistic scattering amplitudes in fermionic field theories, exemplified by the massive Gross-Neveu model, a theory in two spacetime dimensions with quartic interactions. The algorithm introduces new techniques to meet the additional challenges posed by the characteristics of fermionic fields, and its run time is polynomial in the desired precision and the energy. Thus, it constitutes further progress towards an efficient quantum algorithm for simulating the Standard Model of particle physics.
We show that a Wilson-type discretization of the Gross–Neveu model, a fermionic N-flavor quantum field theory displaying asymptotic freedom and chiral symmetry breaking, can serve as a playground to explore correlated symmetry-protected phases of matter using techniques borrowed from high-energy physics. A large-N study, both in the Hamiltonian and Euclidean formalisms, yields a phase diagram with trivial, topological, and symmetry-broken phases separated by critical lines that meet at a tri-critical point. We benchmark these predictions using tools from condensed matter and quantum information science, which show that the large-N method captures the essence of the phase diagram even at N=1. Moreover, we describe a cold-atom scheme for the quantum simulation of this lattice model, which would allow to explore the single-flavor phase diagram.
We discuss the use of quantum simulation to study an NNN-flavor theory of interacting relativistic fermions in (1+11+11+1) dimensions on noisy intermediate-scale quantum (NISQ)-era machines. The case of two flavors is particularly interesting as it can be mapped to the Hubbard model. We derive the appropriate qubit Hamiltonians and associated quantum circuits. We compare classical simulation and density matrix renormalization group/time-evolving blocked decimation calculations with the results of quantum simulation on various platforms for N=2N=2N=2 and 4-flavors. We demonstrate that the four steps needed for calculations of real-time scattering can be implemented using current NISQ devices.
A quantum algorithm to simulate the real time dynamics of two-flavor massive Gross-Neveu model is presented in Schrodinger picture. We implement the simulation on a classic computer by applying the matrix product state representation. The real time evolutions of up to four particles on a site in initial state are figured out in space-time coordinate. The state evolutions are effectively affected by fermion mass and coupling constant of the model. Especially when the mass of fermion is small enough and the coupling is strong enough, the fundamental fermions evolve synchronistically in space from the two-fermion and four-fermion initial states. These are also the conditions on which the bound states made up of fundamental fermion pairs were found to arise automatically in the literatures.