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:2411.15132 · PRX Quantum

Efficient Eigenstate Preparation in an Integrable Model with Hilbert Space Fragmentation

Roberto Ruiz, Alejandro Sopena, Balázs Pozsgay, Esperanza López

Published Jun 10, 2025• 2 citations
arXivHTMLPDF

Abstract

We demonstrate that eigenstates of selected interacting spin chains can be efficiently prepared using quantum circuits with polynomial depth in system size and particle number. While this is well established for free-fermionic spin chains, we extend this result to the folded X X Z model, an integrable rigid-rod deformation of the X X model with simple interactions that exhibits Hilbert space fragmentation. We construct explicit quantum circuits that efficiently prepare arbitrary eigenstates of this model on an open chain. Error-mitigated noisy simulations with up to 13 qubits and various qubit connectivities achieve a relative error below 5%. As a byproduct, we extend a recent reformulation of the Bethe ansatz as a quantum circuit to open boundary conditions.

DOI
10.1103/g9f9-p8ks

Related papers

Adiabatic quantum state preparation in integrable models

Maximilian Lutz, Lorenzo Piroli, Georgios Styliaris +1 more

We propose applying the adiabatic algorithm to prepare high-energy eigenstates of integrable models on a quantum computer. We first review the standard adiabatic algorithm to prepare ground states in each magnetization sector of the prototypical XXZ Heisenberg chain. Based on the thermodynamic Bethe ansatz, we show that the algorithm circuit depth is polynomial in the number of qubits N , outperforming previous methods explicitly relying on integrability. Next, we propose a protocol to prepare arbitrary eigenstates of integrable models that satisfy certain conditions. For a given target eigenstate, we construct a suitable parent Hamiltonian written in terms of a complete set of local conserved quantities. We propose using such Hamiltonian as an input for an adiabatic algorithm. After benchmarking this construction in the case of the non-interacting XY spin chain, where we can rigorously prove its efficiency, we apply it to prepare arbitrary eigenstates of the Richardson-Gaudin models. In this case, we provide numerical evidence that the circuit depth of our algorithm is polynomial in N for all eigenstates, despite the models being interacting.

arXiv:2503.21741•Quantum•Mar 18, 2026•0 citations

Integrable spin chain with Hilbert space fragmentation and solvable real-time dynamics

Balázs Pozsgay, Tamás Gombor, Arthur Hutsalyuk +3 more

We revisit the so-called folded XXZ model, which was treated earlier by two independent research groups. We argue that this spin-1/2 chain is one of the simplest quantum integrable models, yet it has quite remarkable physical properties. The particles have constant scattering lengths, which leads to a simple treatment of the exact spectrum and the dynamics of the system. The Hilbert space of the model is fragmented, leading to exponentially large degeneracies in the spectrum, such that the exponent depends on the particle content of a given state. We provide an alternative derivation of the Hamiltonian and the conserved charges of the model, including an alternative interpretation of the so-called "dual model" considered earlier. We also construct a nonlocal map that connects the model with the Maassarani-Mathieu spin chain, also known as the SU(3) XX model. We consider the exact solution of the model with periodic and open boundary conditions, and also derive multiple descriptions of the exact thermodynamics of the model. We consider quantum quenches of different types. In one class of problems the dynamics can be treated relatively easily: we compute an example for the real-time dependence of a local observable. In another class of quenches the degeneracies of the model lead to the breakdown of equilibration, and we argue that they can lead to persistent oscillations. We also discuss connections with the TT[over ¯] and hard rod deformations known from quantum field theories.

arXiv:2105.02252•Physical review. E•Oct 5, 2021•87 citations

Deterministic Bethe state preparation

David Raveh, Rafael I. Nepomechie

We present an explicit quantum circuit that prepares an arbitrary U ( 1 ) -eigenstate on a quantum computer, including the exact eigenstates of the spin- 1 / 2 X X Z quantum spin chain with either open or closed boundary conditions. The algorithm is deterministic, does not require ancillary qubits, and does not require QR decompositions. The circuit prepares such an L -qubit state with M down-spins using ( L M ) − 1 multi-controlled rotation gates and 2 M ( L − M ) CNOT-gates.

arXiv:2403.03283•Quantum•Oct 24, 2024•14 citations

Local integrals of motion encoded in a few eigenstates

J. Pawłowski, P. Łydżba, M. Mierzejewski

Many properties of a quantum system can be obtained from just a single eigenstate of its Hamiltonian. For example, a single eigenstate can be used to determine whether a system is integrable or chaotic and, in the latter case, to establish its thermal properties. Focusing on the XXZ model, we show that the local integrals of motion, which lie at the heart of integrability, can also be estimated from a small number of eigenstates. Moreover, as the system size increases, fewer eigenstates are required, so that in the thermodynamic limit, the integrals of motion can be obtained from a vanishingly small fraction of all eigenstates. Interestingly, this property does not extend to integrals of motion arising solely from Hilbert space fragmentation, as found in the folded XXZ model, where the majority of eigenstates has to be used. This represents one of the few fundamental differences known between integrability and Hilbert space fragmentation.

cond-mat.str-elCross-listedarXiv:2603.01859•Mar 3, 2026•0 citations

Preparing Bethe Ansatz Eigenstates on a Quantum Computer

John S. Van Dyke, George S. Barron, Nicholas J. Mayhall +2 more

Several quantum many-body models in one dimension possess exact solutions via the Bethe ansatz method, which has been highly successful for understanding their behavior.Nevertheless, there remain physical properties of such models for which analytic results are unavailable and which are also not well described by approximate numerical methods.Preparing Bethe ansatz eigenstates directly on a quantum computer would allow straightforward extraction of these quantities via measurement.We present a quantum algorithm for preparing Bethe ansatz eigenstates of the spin-1/2 XXZ spin chain that correspond to real-valued solutions of the Bethe equations.The algorithm is polynomial in the number of T gates and the circuit depth, with modest constant prefactors.Although the algorithm is probabilistic, with a success rate that decreases with increasing eigenstate energy, we employ amplitude amplification to boost the success probability.The resource requirements for our approach are lower than for other state-of-the-art quantum simulation algorithms for small error-corrected devices and thus may offer an alternative and computationally less demanding demonstration of quantum advantage for physically relevant problems.

arXiv:2103.13388•PRX Quantum•Nov 9, 2021•41 citations

Preparing exact eigenstates of the open XXZ chain on a quantum computer

John S. Van Dyke, Edwin Barnes, Sophia E. Economou +1 more

Abstract The open spin-1/2 XXZ spin chain with diagonal boundary magnetic fields is the paradigmatic example of a quantum integrable model with open boundary conditions. We formulate a quantum algorithm for preparing Bethe states of this model, corresponding to real solutions of the Bethe equations. The algorithm is probabilistic, with a success probability that decreases with the number of down spins. For a Bethe state of L spins with M down spins, which contains a total of L M 2 M M ! terms, the algorithm requires L + M 2 + 2 M qubits.

arXiv:2109.05607•Journal of Physics A Mathematical and Theoretical•Dec 23, 2021•27 citations
We demonstrate that eigenstates of selected interacting spin chains can be efficiently prepared using quantum circuits with polynomial depth in system size and particle number. While this is well established for free-fermionic spin chains, we extend this result to the folded X X Z model, an integrable rigid-rod deformation of the X X model with simple interactions that exhibits Hilbert space fragmentation. We construct explicit quantum circuits that efficiently prepare arbitrary eigenstates of this model on an open chain. Error-mitigated noisy simulations with up to 13 qubits and various qubit connectivities achieve a relative error below 5%. As a byproduct, we extend a recent reformulation of the Bethe ansatz as a quantum circuit to open boundary conditions.
We propose applying the adiabatic algorithm to prepare high-energy eigenstates of integrable models on a quantum computer. We first review the standard adiabatic algorithm to prepare ground states in each magnetization sector of the prototypical XXZ Heisenberg chain. Based on the thermodynamic Bethe ansatz, we show that the algorithm circuit depth is polynomial in the number of qubits N , outperforming previous methods explicitly relying on integrability. Next, we propose a protocol to prepare arbitrary eigenstates of integrable models that satisfy certain conditions. For a given target eigenstate, we construct a suitable parent Hamiltonian written in terms of a complete set of local conserved quantities. We propose using such Hamiltonian as an input for an adiabatic algorithm. After benchmarking this construction in the case of the non-interacting XY spin chain, where we can rigorously prove its efficiency, we apply it to prepare arbitrary eigenstates of the Richardson-Gaudin models. In this case, we provide numerical evidence that the circuit depth of our algorithm is polynomial in N for all eigenstates, despite the models being interacting.
We revisit the so-called folded XXZ model, which was treated earlier by two independent research groups. We argue that this spin-1/2 chain is one of the simplest quantum integrable models, yet it has quite remarkable physical properties. The particles have constant scattering lengths, which leads to a simple treatment of the exact spectrum and the dynamics of the system. The Hilbert space of the model is fragmented, leading to exponentially large degeneracies in the spectrum, such that the exponent depends on the particle content of a given state. We provide an alternative derivation of the Hamiltonian and the conserved charges of the model, including an alternative interpretation of the so-called "dual model" considered earlier. We also construct a nonlocal map that connects the model with the Maassarani-Mathieu spin chain, also known as the SU(3) XX model. We consider the exact solution of the model with periodic and open boundary conditions, and also derive multiple descriptions of the exact thermodynamics of the model. We consider quantum quenches of different types. In one class of problems the dynamics can be treated relatively easily: we compute an example for the real-time dependence of a local observable. In another class of quenches the degeneracies of the model lead to the breakdown of equilibration, and we argue that they can lead to persistent oscillations. We also discuss connections with the TT[over ¯] and hard rod deformations known from quantum field theories.
Many properties of a quantum system can be obtained from just a single eigenstate of its Hamiltonian. For example, a single eigenstate can be used to determine whether a system is integrable or chaotic and, in the latter case, to establish its thermal properties. Focusing on the XXZ model, we show that the local integrals of motion, which lie at the heart of integrability, can also be estimated from a small number of eigenstates. Moreover, as the system size increases, fewer eigenstates are required, so that in the thermodynamic limit, the integrals of motion can be obtained from a vanishingly small fraction of all eigenstates. Interestingly, this property does not extend to integrals of motion arising solely from Hilbert space fragmentation, as found in the folded XXZ model, where the majority of eigenstates has to be used. This represents one of the few fundamental differences known between integrability and Hilbert space fragmentation.
Several quantum many-body models in one dimension possess exact solutions via the Bethe ansatz method, which has been highly successful for understanding their behavior.Nevertheless, there remain physical properties of such models for which analytic results are unavailable and which are also not well described by approximate numerical methods.Preparing Bethe ansatz eigenstates directly on a quantum computer would allow straightforward extraction of these quantities via measurement.We present a quantum algorithm for preparing Bethe ansatz eigenstates of the spin-1/2 XXZ spin chain that correspond to real-valued solutions of the Bethe equations.The algorithm is polynomial in the number of T gates and the circuit depth, with modest constant prefactors.Although the algorithm is probabilistic, with a success rate that decreases with increasing eigenstate energy, we employ amplitude amplification to boost the success probability.The resource requirements for our approach are lower than for other state-of-the-art quantum simulation algorithms for small error-corrected devices and thus may offer an alternative and computationally less demanding demonstration of quantum advantage for physically relevant problems.
Abstract The open spin-1/2 XXZ spin chain with diagonal boundary magnetic fields is the paradigmatic example of a quantum integrable model with open boundary conditions. We formulate a quantum algorithm for preparing Bethe states of this model, corresponding to real solutions of the Bethe equations. The algorithm is probabilistic, with a success probability that decreases with the number of down spins. For a Bethe state of L spins with M down spins, which contains a total of L M 2 M M ! terms, the algorithm requires L + M 2 + 2 M qubits.