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:2211.05140 · Physical Review Letters

Extracting Quantum Many-Body Scarred Eigenstates with Matrix Product States

Shunyao Zhang, Yuan Dong, Thomas Iadecola, Shenglong Xu, Dong-Ling Deng

Published Jul 11, 2023• 17 citations
arXivHTMLPDF

Abstract

Quantum many-body scarred systems host nonthermal excited eigenstates immersed in a sea of thermal ones. In cases where exact expressions for these special eigenstates are not known, it is computationally demanding to distinguish them from their exponentially many thermal neighbors. We propose a matrix-product-state (MPS) algorithm, dubbed DMRG-S, to extract such states at system sizes far beyond the scope of exact diagonalization. Using this technique, we obtain scarred eigenstates in Rydberg-blockaded chains of up to 80 sites and perform a finite-size scaling study to address the lingering question of the stability for the Néel state revivals in the thermodynamic limit. Our method also provides a systematic way to obtain exact MPS representations for scarred eigenstates near the target energy without a priori knowledge. In particular, we find several new scarred eigenstates with exact MPS representations in kinetically constrained spin and clock models. The combination of numerical and analytical investigations in our work provides a new methodology for future studies of quantum many-body scars.

DOI
10.1103/physrevlett.131.020402

Related papers

Exact many-body scars and their stability in constrained quantum chains

Federica Maria Surace, Matteo Votto, Eduardo González Lazo +3 more

Quantum scars are nonthermal eigenstates characterized by low entanglement entropy, initially detected in systems subject to nearest-neighbor Rydberg blockade, the so-called PXP model. While most of these special eigenstates elude an analytical description and seem to hybridize with nearby thermal eigenstates for large systems, some of them can be written as matrix product states with size-independent bond dimension. We study the response of these exact quantum scars to perturbations by analyzing the scaling of the fidelity susceptibility with system size. We find that some of them are anomalously stable at first order in perturbation theory, in sharp contrast to the eigenstate thermalization hypothesis. However, this stability seems to break down when all orders are taken into account. We further investigate models with larger blockade radius and find a set of exact quantum scars that we write down analytically and compare with the PXP exact eigenstates. We show that they exhibit the same robustness against perturbations at first order.

arXiv:2011.08218•Physical review. B./Physical review. B•Mar 16, 2021•63 citations

Quantum many-body scars through the lens of correlation matrix

Zhiyuan Yao, Pengfei Zhang

Quantum many-body scars (QMBS)---rare eigenstates that evade thermalization---are typically characterized by their low entanglement entropies compared to surrounding thermal eigenstates. However, due to finite-size effects in systems accessible via exact diagonalization, this measure can be ambiguous. To address this limitation, we propose using the correlation matrix spectrum as an alternative probe to identify QMBS. In cases of exact QMBS that either have known analytic expressions or are captured by various frameworks of QMBS, we find that the dimensionality of the null space of the correlation matrix---an integer value, and thus immune to finite-size effects---can qualitatively identify QMBS. Beyond serving as a diagnostic tool, the correlation matrix method enables the manipulation of the QMBS subspace. For approximate QMBS, such as those in the PXP model, we observe that the correlation matrix spectrum features numerous approximate zero eigenvalues, thereby distinguishing these states. We demonstrate the effectiveness and utility of this method with several paradigmatic QMBS examples.

arXiv:2410.21812•Physical review. B./Physical review. B•Sep 15, 2025•4 citations

Thermalization of exact quantum many-body scars in spin-1 XY chain under perturbation

Himadri Halder

Quantum many-body scars are special eigenstates that violate the eigenstate thermalization hypothesis while residing at finite energy density along with thermalizing eigenstates. The spin-1 XY model is known to host a family of such exceptional states originating from long-lived quasiparticle excitations that exhibit anomalously low entanglement entropy and long-time periodic revivals, resulting in weak ergodicity breaking. We study the stability of these scarred states against typical U(1) symmetry preserving perturbation in the XY chain. While perturbation theory can describe the deformed scar states at small system sizes, finite-size scaling of the perturbation matrix elements indicate that the scars ultimately thermalize in larger chains. Nonetheless, we demonstrate that the long-range order associated with the scars decays under the perturbation, and we estimate the relaxation timescale of oscillatory dynamics in certain local observables to be of order $λ^{-2}$, where $λ$ is the perturbation strength.

cond-mat.str-elCross-listedarXiv:2511.17476•Nov 24, 2025•0 citations

Gaussian state approximation of quantum many-body scars

Wouter Buijsman, Yevgeny Bar Lev

Quantum many-body scars are atypical, highly nonthermal eigenstates embedded in a sea of thermal eigenstates that have been observed in, for example, kinetically constrained quantum many-body models. These special eigenstates are characterized by a bipartite entanglement entropy that scales as most logarithmically with the subsystem size. We use numerical optimization techniques to investigate if quantum many-body scars of the experimentally relevant PXP model can be well approximated by Gaussian states. Gaussian states are described by a number of parameters that scales quadratically with system size, thereby having a much lower complexity than generic quantum many-body states, for which this number scales exponentially. We find that while quantum many-body scars can typically be well approximated by (symmetrized) Gaussian states, this is not the case for ergodic (thermal) eigenstates. This observation suggests that the non-ergodic part of the PXP Hamiltonian is related to certain quadratic parent Hamiltonians, thereby hinting on the origin of the quantum many-body scars.

arXiv:2307.11816•SciPost Physics•Aug 15, 2024•0 citations

Exact quantum many-body scars in higher-spin kinetically constrained models

Yuan Dong, Shunyao Zhang, Dong-Ling Deng

We discover a variety of exact quantum many-body scars in higher-spin kinetically constrained models through the recently developed DMRG-S algorithm [Zhang et al., Phys. Rev. Lett. 131, 020402 (2023)]. Specifically, for the higher-spin PXP model on arbitrary bipartite lattices of any spatial dimension, we find exact many-body scars that are equidistantly spaced in the energy spectrum and exhibit similar structures to the ground state of the Affleck-Kennedy-Lieb-Tasaki model. For the one-dimensional Fermi-Hubbard model with a tilted potential in a certain parameter regime, whose effective model is equivalent to a kinetically constrained spin model with four degrees of freedom on each site, we find several many-body scars at energy $E=0$ and $E=\ifmmode\pm\else\textpm\fi{}\sqrt{2}$ that can be exactly represented as matrix product states with finite bond dimensions. Our results demonstrate that larger local degrees of freedom in the kinetically constrained models provide a much broader space for the emergence of quantum many-body scars and weak ergodicity breaking.

arXiv:2307.06357•Physical review. B./Physical review. B•Nov 16, 2023•5 citations

Preparing quantum many-body scar states on quantum computers

Erik Gustafson, Andy C. Y. Li, Abid Khan +6 more

Quantum many-body scar states are highly excited eigenstates of many-body systems that exhibit atypical entanglement and correlation properties relative to typical eigenstates at the same energy density. Scar states also give rise to infinitely long-lived coherent dynamics when the system is prepared in a special initial state having finite overlap with them. Many models with exact scar states have been constructed, but the fate of scarred eigenstates and dynamics when these models are perturbed is difficult to study with classical computational techniques. In this work, we propose state preparation protocols that enable the use of quantum computers to study this question. We present protocols both for individual scar states in a particular model, as well as superpositions of them that give rise to coherent dynamics. For superpositions of scar states, we present both a system-size-linear depth unitary and a finite-depth nonunitary state preparation protocol, the latter of which uses measurement and postselection to reduce the circuit depth. For individual scarred eigenstates, we formulate an exact state preparation approach based on matrix product states that yields quasipolynomial-depth circuits, as well as a variational approach with a polynomial-depth ansatz circuit. We also provide proof of principle state-preparation demonstrations on superconducting quantum hardware.

arXiv:2301.08226•Quantum•Nov 7, 2023•21 citations
Quantum many-body scarred systems host nonthermal excited eigenstates immersed in a sea of thermal ones. In cases where exact expressions for these special eigenstates are not known, it is computationally demanding to distinguish them from their exponentially many thermal neighbors. We propose a matrix-product-state (MPS) algorithm, dubbed DMRG-S, to extract such states at system sizes far beyond the scope of exact diagonalization. Using this technique, we obtain scarred eigenstates in Rydberg-blockaded chains of up to 80 sites and perform a finite-size scaling study to address the lingering question of the stability for the Néel state revivals in the thermodynamic limit. Our method also provides a systematic way to obtain exact MPS representations for scarred eigenstates near the target energy without a priori knowledge. In particular, we find several new scarred eigenstates with exact MPS representations in kinetically constrained spin and clock models. The combination of numerical and analytical investigations in our work provides a new methodology for future studies of quantum many-body scars.
Quantum scars are nonthermal eigenstates characterized by low entanglement entropy, initially detected in systems subject to nearest-neighbor Rydberg blockade, the so-called PXP model. While most of these special eigenstates elude an analytical description and seem to hybridize with nearby thermal eigenstates for large systems, some of them can be written as matrix product states with size-independent bond dimension. We study the response of these exact quantum scars to perturbations by analyzing the scaling of the fidelity susceptibility with system size. We find that some of them are anomalously stable at first order in perturbation theory, in sharp contrast to the eigenstate thermalization hypothesis. However, this stability seems to break down when all orders are taken into account. We further investigate models with larger blockade radius and find a set of exact quantum scars that we write down analytically and compare with the PXP exact eigenstates. We show that they exhibit the same robustness against perturbations at first order.
Quantum many-body scars (QMBS)---rare eigenstates that evade thermalization---are typically characterized by their low entanglement entropies compared to surrounding thermal eigenstates. However, due to finite-size effects in systems accessible via exact diagonalization, this measure can be ambiguous. To address this limitation, we propose using the correlation matrix spectrum as an alternative probe to identify QMBS. In cases of exact QMBS that either have known analytic expressions or are captured by various frameworks of QMBS, we find that the dimensionality of the null space of the correlation matrix---an integer value, and thus immune to finite-size effects---can qualitatively identify QMBS. Beyond serving as a diagnostic tool, the correlation matrix method enables the manipulation of the QMBS subspace. For approximate QMBS, such as those in the PXP model, we observe that the correlation matrix spectrum features numerous approximate zero eigenvalues, thereby distinguishing these states. We demonstrate the effectiveness and utility of this method with several paradigmatic QMBS examples.
Quantum many-body scars are special eigenstates that violate the eigenstate thermalization hypothesis while residing at finite energy density along with thermalizing eigenstates. The spin-1 XY model is known to host a family of such exceptional states originating from long-lived quasiparticle excitations that exhibit anomalously low entanglement entropy and long-time periodic revivals, resulting in weak ergodicity breaking. We study the stability of these scarred states against typical U(1) symmetry preserving perturbation in the XY chain. While perturbation theory can describe the deformed scar states at small system sizes, finite-size scaling of the perturbation matrix elements indicate that the scars ultimately thermalize in larger chains. Nonetheless, we demonstrate that the long-range order associated with the scars decays under the perturbation, and we estimate the relaxation timescale of oscillatory dynamics in certain local observables to be of order λ−2λ^{-2}λ−2, where λλλ is the perturbation strength.
Quantum many-body scars are atypical, highly nonthermal eigenstates embedded in a sea of thermal eigenstates that have been observed in, for example, kinetically constrained quantum many-body models. These special eigenstates are characterized by a bipartite entanglement entropy that scales as most logarithmically with the subsystem size. We use numerical optimization techniques to investigate if quantum many-body scars of the experimentally relevant PXP model can be well approximated by Gaussian states. Gaussian states are described by a number of parameters that scales quadratically with system size, thereby having a much lower complexity than generic quantum many-body states, for which this number scales exponentially. We find that while quantum many-body scars can typically be well approximated by (symmetrized) Gaussian states, this is not the case for ergodic (thermal) eigenstates. This observation suggests that the non-ergodic part of the PXP Hamiltonian is related to certain quadratic parent Hamiltonians, thereby hinting on the origin of the quantum many-body scars.
We discover a variety of exact quantum many-body scars in higher-spin kinetically constrained models through the recently developed DMRG-S algorithm [Zhang et al., Phys. Rev. Lett. 131, 020402 (2023)]. Specifically, for the higher-spin PXP model on arbitrary bipartite lattices of any spatial dimension, we find exact many-body scars that are equidistantly spaced in the energy spectrum and exhibit similar structures to the ground state of the Affleck-Kennedy-Lieb-Tasaki model. For the one-dimensional Fermi-Hubbard model with a tilted potential in a certain parameter regime, whose effective model is equivalent to a kinetically constrained spin model with four degrees of freedom on each site, we find several many-body scars at energy E=0E=0E=0 and $E=\ifmmode\pm\else\textpm\fi{}\sqrt{2}$ that can be exactly represented as matrix product states with finite bond dimensions. Our results demonstrate that larger local degrees of freedom in the kinetically constrained models provide a much broader space for the emergence of quantum many-body scars and weak ergodicity breaking.
Quantum many-body scar states are highly excited eigenstates of many-body systems that exhibit atypical entanglement and correlation properties relative to typical eigenstates at the same energy density. Scar states also give rise to infinitely long-lived coherent dynamics when the system is prepared in a special initial state having finite overlap with them. Many models with exact scar states have been constructed, but the fate of scarred eigenstates and dynamics when these models are perturbed is difficult to study with classical computational techniques. In this work, we propose state preparation protocols that enable the use of quantum computers to study this question. We present protocols both for individual scar states in a particular model, as well as superpositions of them that give rise to coherent dynamics. For superpositions of scar states, we present both a system-size-linear depth unitary and a finite-depth nonunitary state preparation protocol, the latter of which uses measurement and postselection to reduce the circuit depth. For individual scarred eigenstates, we formulate an exact state preparation approach based on matrix product states that yields quasipolynomial-depth circuits, as well as a variational approach with a polynomial-depth ansatz circuit. We also provide proof of principle state-preparation demonstrations on superconducting quantum hardware.