© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2609.38848

Quantum quenches of scar states in the Affleck-Kennedy-Lieb-Tasaki model via Clifford augmented tensor network simulation

Azar C. Nakhl, Ben Harper, Muhammad Usman, Thomas Quella

cond-mat.str-elquant-phSubmitted Sep 30, 2026• Announced Oct 1, 2026• 0 citations• Revised Oct 28, 2026 (v2)

v1 · Sep 30, 2026 v2 · Oct 28, 2026

Abstract

Clifford augmented methods have emerged as a powerful technique to simulate quantum circuits and many-body systems by exploiting the stabiliser structure found in said systems, whilst additionally enabling one to readily quantify the non-stabilizerness (i.e.\ magic) present. In this work, we utilise a qudit Clifford augmented simulation method to probe the presence of magic in the Affleck-Kennedy-Lieb-Tasaki (AKLT) model, finding that its ground state and scar states are highly magical as quantified by their Stabiliser Rényi Entropy (SRE). Nevertheless, owing to its low entanglement structure, these states are still readily expressible as a Clifford augmented Matrix Product State with low bond-dimension and moreover may be time-evolved using the Clifford augmented time-dependent variational principle (TDVP) algorithm. We utilise Clifford augmented TDVP to study the string order dynamics in the scar states of the AKLT model after a quench within the Haldane phase to the anti-ferromagnetic Heisenberg Hamiltonian. We observe that there is a decay of string order following the quench, which is more pronounced with increasing bimagnon number. This finding highlights the sensitivity of string order in the AKLT scar states to symmetry-preserving quenches that break the underlying restricted spectrum generating algebra, even when the initial and final Hamiltonians belong to the same topological phase.

Comments: 11+3 pages, 8 figures

Related papers

Disentangling Haldane Phase by Generalized Clifford Circuits

Minsoo Kim, Changhun Oh, Donghoon Kim

Disentangling transformations play a central role in the classical simulation of quantum many-body systems, yet their analytic structure and underlying mechanism remain largely unexplored. Here, we study the structure of the disentangler in the Haldane phase of spin-1 systems using generalized Clifford circuits. To this end, we extend the Clifford-circuit-augmented matrix product states (CAMPS)-based density-matrix renormalization group (DMRG) method to spin-1 systems. Within this framework, we find that the local disentanglers optimized for the Haldane phase implement the generalized Kramers--Wannier (KW) transformation, and we analytically verify its optimality for the Affleck--Kennedy--Lieb--Tasaki (AKLT) state. Beyond reducing entanglement, the KW transformation maps the Haldane phase to a phase with spontaneously broken $\mathbb{Z}_{2}$ symmetry. This mapping is distinct from the Kennedy--Tasaki transformation and provides a new unitary route from symmetry-protected topological order to symmetry breaking.

quant-pharXiv:2607.03939•Jul 7, 2026•0 citations

Connection between quantum-many-body scars and the Affleck–Kennedy–Lieb–Tasaki model from the viewpoint of embedded Hamiltonians

Naoto Shiraishi

Abstract We elucidate the deep connection between the PXP model, which is a standard model of quantum many-body scars, and the Affleck–Kennedy–Lieb–Tasaki (AKLT) Hamiltonian. Using the framework of embedded Hamiltonians, we establish the connection between the PXP Hamiltonian and the AKLT Hamiltonian, which clarifies the reason why the PXP Hamiltonian has nonthermal energy eigenstates similar to the AKLT state. Through this analysis, we find that the presence of such nonthermal energy eigenstates reflects the symmetry in the AKLT Hamiltonian.

arXiv:1904.05182•Journal of Statistical Mechanics Theory and Experiment•Aug 1, 2019•87 citations

Destruction of string order after a quantum quench

Marcello Calvanese Strinati, Leonardo Mazza, Manuel Endres +2 more

We investigate the evolution of string order in a spin-1 chain following a quantum quench. After initializing the chain in the Affleck-Kennedy-Lieb-Tasaki state, we analyze in detail how string order evolves as a function of time at different length scales. The Hamiltonian after the quench is chosen either to preserve or to suddenly break the symmetry which ensures the presence of string order. Depending on which of these two situations arises, string order is either preserved or lost even at infinitesimal times in the thermodynamic limit. The fact that nonlocal order may be abruptly destroyed, what we call string-order melting, makes it qualitatively different from typical order parameters in the manner of Landau. This situation is thoroughly characterized by means of numerical simulations based on matrix product states algorithms and analytical studies based on a short-time expansion for several simplified models.

arXiv:1604.02823•Physical review. B./Physical review. B•Jul 5, 2016•17 citations

Clifford Circuits Augmented Time-Dependent Variational Principle

Xiangjian Qian, Jiale Huang, Mingpu Qin

The recently proposed Clifford circuits augmented matrix product states (CA-MPSs) [X. Qian, J. Huang, and M. Qin, Phys. Rev. Lett. 133, 190402 (2024)PRLTAO0031-900710.1103/PhysRevLett.133.190402] seamlessly augment density matrix renormalization groups with Clifford circuits. In CA-MPSs, the entanglement from stabilizers is transferred to the Clifford circuits, which can be easily handled according to the Gottesman-Knill theorem. As a result, an MPS needs only to deal with the nonstabilizer entanglement, which largely reduces the bond dimension and the resource required for the accurate simulation of many-body systems. In this Letter, we generalize CA-MPSs to the framework of the time-dependent variational principle (TDVP) for time evolution simulations. In this method, we apply Clifford circuits to the resulting MPS in each TDVP step with a two-site sweeping process similar as in density matrix renormalization groups, aiming at reducing the entanglement entropy in the MPS, and the Hamiltonian is transformed accordingly using the chosen Clifford circuits. Similar as in CA-MPSs, the Clifford circuits does not increase the number of terms in the Hamiltonian, which makes the overhead very small in the new method. We test this method in XXZ chain, 2D Heisenberg model, and Kitaev honeycomb model. The results show that the Clifford circuits augmented TDVP method can reduce the entanglement entropy in the time evolution process and hence makes the simulation reliable for a longer time. The Clifford circuits augmented time-dependent variational principle provides a useful tool for the simulation of the time evolution process of many-body systems in the future.

arXiv:2407.03202•Physical Review Letters•Apr 15, 2025•22 citations

Affleck-Kennedy-Lieb-Tasaki State on a Honeycomb Lattice from t2g Orbitals

Maciej Koch-Janusz, Daniel I. Khomskii, Eran Sela

The two-dimensional Affleck-Kennedy-Lieb-Tasaki (AKLT) model on a honeycomb lattice has been shown to be a universal resource for quantum computation. In this valence bond solid, however, the spin interactions involve higher powers of the Heisenberg coupling (S[over →](I)·S[over →](j))(n), making these states seemingly unrealistic on bipartite lattices, where one expects a simple antiferromagnetic order. We show that those interactions can be generated by orbital physics in multiorbital Mott insulators. We focus on t(2g) electrons on the honeycomb lattice and propose a physical realization of the spin-3/2 AKLT state. We find a phase transition from the AKLT to the Néel state on increasing Hund's rule coupling, which is confirmed by density matrix renormalization group simulations. An experimental signature of the AKLT state consists of protected, free S=1/2 spins on lattice vacancies, which may be detected in the spin susceptibility.

arXiv:1501.03512•Physical Review Letters•Jun 18, 2015•22 citations

Symmetry-protected topological scar subspaces

Chihiro Matsui, Thomas Quella, Naoto Tsuji

We propose a framework that extends the notion of symmetry-protected topological properties beyond the ground-state paradigm to dynamically isolated subspaces formed by exceptional non-thermal energy eigenstates of non-integrable systems, known as quantum many-body scars (QMBS). We introduce the concept of a symmetry-protected topological (SPT) scar subspace -- a Hilbert subspace stabilized by a restricted spectrum-generating algebra (rSGA) while being protected by on-site, inversion, and time-reversal symmetries. QMBS often admit a non-interacting quasiparticle description, which enables matrix-product representations with small bond dimension. Although individual QMBS do not necessarily retain the protecting symmetries of the Hamiltonian, we show that the subspace formed by the symmetry-connected QMBS does retain them, giving rise to consistently emerging topological properties across the entire scar subspace. Using the spin-$1$ Affleck--Kennedy--Lieb--Tasaki (AKLT) model, we demonstrate that its bimagnon scar subspace reflects the topological properties of the SPT ground state, as evidenced by the appropriate bond-space symmetry representations, the expected topological response, and the numerically verified long-range string order. Our findings indicate that scar subspaces can inherit -- and in inhomogeneous cases systematically modify -- the topological character of the SPT ground state, offering a new and experimentally accessible platform for probing symmetry-protected topology beyond the ground-state regime.

cond-mat.str-elCross-listedarXiv:2512.11216•Dec 15, 2025•0 citations

Clifford augmented methods have emerged as a powerful technique to simulate quantum circuits and many-body systems by exploiting the stabiliser structure found in said systems, whilst additionally enabling one to readily quantify the non-stabilizerness (i.e.\ magic) present. In this work, we utilise a qudit Clifford augmented simulation method to probe the presence of magic in the Affleck-Kennedy-Lieb-Tasaki (AKLT) model, finding that its ground state and scar states are highly magical as quantified by their Stabiliser Rényi Entropy (SRE). Nevertheless, owing to its low entanglement structure, these states are still readily expressible as a Clifford augmented Matrix Product State with low bond-dimension and moreover may be time-evolved using the Clifford augmented time-dependent variational principle (TDVP) algorithm. We utilise Clifford augmented TDVP to study the string order dynamics in the scar states of the AKLT model after a quench within the Haldane phase to the anti-ferromagnetic Heisenberg Hamiltonian. We observe that there is a decay of string order following the quench, which is more pronounced with increasing bimagnon number. This finding highlights the sensitivity of string order in the AKLT scar states to symmetry-preserving quenches that break the underlying restricted spectrum generating algebra, even when the initial and final Hamiltonians belong to the same topological phase.

Disentangling transformations play a central role in the classical simulation of quantum many-body systems, yet their analytic structure and underlying mechanism remain largely unexplored. Here, we study the structure of the disentangler in the Haldane phase of spin-1 systems using generalized Clifford circuits. To this end, we extend the Clifford-circuit-augmented matrix product states (CAMPS)-based density-matrix renormalization group (DMRG) method to spin-1 systems. Within this framework, we find that the local disentanglers optimized for the Haldane phase implement the generalized Kramers--Wannier (KW) transformation, and we analytically verify its optimality for the Affleck--Kennedy--Lieb--Tasaki (AKLT) state. Beyond reducing entanglement, the KW transformation maps the Haldane phase to a phase with spontaneously broken

\mathbb{Z}_{2}

symmetry. This mapping is distinct from the Kennedy--Tasaki transformation and provides a new unitary route from symmetry-protected topological order to symmetry breaking.

Abstract We elucidate the deep connection between the PXP model, which is a standard model of quantum many-body scars, and the Affleck–Kennedy–Lieb–Tasaki (AKLT) Hamiltonian. Using the framework of embedded Hamiltonians, we establish the connection between the PXP Hamiltonian and the AKLT Hamiltonian, which clarifies the reason why the PXP Hamiltonian has nonthermal energy eigenstates similar to the AKLT state. Through this analysis, we find that the presence of such nonthermal energy eigenstates reflects the symmetry in the AKLT Hamiltonian.

We investigate the evolution of string order in a spin-1 chain following a quantum quench. After initializing the chain in the Affleck-Kennedy-Lieb-Tasaki state, we analyze in detail how string order evolves as a function of time at different length scales. The Hamiltonian after the quench is chosen either to preserve or to suddenly break the symmetry which ensures the presence of string order. Depending on which of these two situations arises, string order is either preserved or lost even at infinitesimal times in the thermodynamic limit. The fact that nonlocal order may be abruptly destroyed, what we call string-order melting, makes it qualitatively different from typical order parameters in the manner of Landau. This situation is thoroughly characterized by means of numerical simulations based on matrix product states algorithms and analytical studies based on a short-time expansion for several simplified models.

The recently proposed Clifford circuits augmented matrix product states (CA-MPSs) [X. Qian, J. Huang, and M. Qin, Phys. Rev. Lett. 133, 190402 (2024)PRLTAO0031-900710.1103/PhysRevLett.133.190402] seamlessly augment density matrix renormalization groups with Clifford circuits. In CA-MPSs, the entanglement from stabilizers is transferred to the Clifford circuits, which can be easily handled according to the Gottesman-Knill theorem. As a result, an MPS needs only to deal with the nonstabilizer entanglement, which largely reduces the bond dimension and the resource required for the accurate simulation of many-body systems. In this Letter, we generalize CA-MPSs to the framework of the time-dependent variational principle (TDVP) for time evolution simulations. In this method, we apply Clifford circuits to the resulting MPS in each TDVP step with a two-site sweeping process similar as in density matrix renormalization groups, aiming at reducing the entanglement entropy in the MPS, and the Hamiltonian is transformed accordingly using the chosen Clifford circuits. Similar as in CA-MPSs, the Clifford circuits does not increase the number of terms in the Hamiltonian, which makes the overhead very small in the new method. We test this method in XXZ chain, 2D Heisenberg model, and Kitaev honeycomb model. The results show that the Clifford circuits augmented TDVP method can reduce the entanglement entropy in the time evolution process and hence makes the simulation reliable for a longer time. The Clifford circuits augmented time-dependent variational principle provides a useful tool for the simulation of the time evolution process of many-body systems in the future.

The two-dimensional Affleck-Kennedy-Lieb-Tasaki (AKLT) model on a honeycomb lattice has been shown to be a universal resource for quantum computation. In this valence bond solid, however, the spin interactions involve higher powers of the Heisenberg coupling (S[over →](I)·S[over →](j))(n), making these states seemingly unrealistic on bipartite lattices, where one expects a simple antiferromagnetic order. We show that those interactions can be generated by orbital physics in multiorbital Mott insulators. We focus on t(2g) electrons on the honeycomb lattice and propose a physical realization of the spin-3/2 AKLT state. We find a phase transition from the AKLT to the Néel state on increasing Hund's rule coupling, which is confirmed by density matrix renormalization group simulations. An experimental signature of the AKLT state consists of protected, free S=1/2 spins on lattice vacancies, which may be detected in the spin susceptibility.

We propose a framework that extends the notion of symmetry-protected topological properties beyond the ground-state paradigm to dynamically isolated subspaces formed by exceptional non-thermal energy eigenstates of non-integrable systems, known as quantum many-body scars (QMBS). We introduce the concept of a symmetry-protected topological (SPT) scar subspace -- a Hilbert subspace stabilized by a restricted spectrum-generating algebra (rSGA) while being protected by on-site, inversion, and time-reversal symmetries. QMBS often admit a non-interacting quasiparticle description, which enables matrix-product representations with small bond dimension. Although individual QMBS do not necessarily retain the protecting symmetries of the Hamiltonian, we show that the subspace formed by the symmetry-connected QMBS does retain them, giving rise to consistently emerging topological properties across the entire scar subspace. Using the spin-

1

Affleck--Kennedy--Lieb--Tasaki (AKLT) model, we demonstrate that its bimagnon scar subspace reflects the topological properties of the SPT ground state, as evidenced by the appropriate bond-space symmetry representations, the expected topological response, and the numerically verified long-range string order. Our findings indicate that scar subspaces can inherit -- and in inhomogeneous cases systematically modify -- the topological character of the SPT ground state, offering a new and experimentally accessible platform for probing symmetry-protected topology beyond the ground-state regime.