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:2511.20619

Extracting conserved operators from a projected entangled pair state

Wen-Tao Xu, Miguel Frías Pérez, Mingru Yang

quant-phcond-mat.stat-mechcond-mat.str-elSubmitted Nov 25, 2025• Announced Nov 26, 2025• 0 citations
v1 · Nov 25, 2025v2 · Apr 14, 2026
arXivHTMLPDF

Abstract

Given a tensor network state, how can we determine conserved operators (including Hamiltonians) for which the state is an eigenstate? We answer this question by presenting a method to extract geometrically $k$-local conserved operators that have the given infinite projected entangled pair state (iPEPS) in 2D as an (approximate) eigenstate. The key ingredient is the evaluation of the static structure factors of multi-site operators through differentiating the generating function. Despite the approximation errors, we show that our method is still able to extract from exact or variational iPEPS to good precision both frustration-free and non-frustration-free parent Hamiltonians that are beyond the standard construction and obtain better locality. In particular, we find a 4-site-plaquette local Hamiltonian that approximately has the short-range RVB state as the ground state. Moreover, we find a Hamiltonian that has the deformed toric code state at any string tension as excited eigenstates at the same energy, which might be potential candidates for quantum many-body scars.

Comments
5+12 pages

Related papers

Variational methods for characterizing matrix product operator symmetries

Anna Francuz, Laurens Lootens, Frank Verstraete +1 more

We present a method of extracting information about topological order from the ground state of a strongly correlated two-dimensional system represented by an infinite projected entangled pair state (iPEPS). As in previous works [A. Francuz et al., Phys. Rev. B 101, 041108(R) (2020) and A. Francuz and J. Dziarmaga, Phys. Rev. B 102, 235112 (2020)] we begin by determining symmetries of the iPEPS represented by infinite matrix product operators (iMPO) that map between the different iPEPS transfer matrix fixed points, to which we apply the fundamental theorem of matrix product states to find zipper tensors between products of iMPO's that encode fusion properties of the anyons. The zippers can be combined to extract topological $F$ symbols of the underlying fusion category, which unequivocally identifies the topological order of the ground state. We bring the $F$ symbols to the canonical gauge, and also compute the Drinfeld center of this unitary fusion category to extract the topological $S$ and $T$ matrices encoding mutual statistics and self-statistics of the emergent anyons. The algorithm is applied to Abelian toric code, Kitaev model, double semion, and twisted quantum double of ${Z}_{3}$, as well as to non-Abelian double Fibonacci, double Ising, and quantum double of ${S}_{3}$ and $\mathrm{Rep}({S}_{3})$ string-net models.

arXiv:2107.05265•Physical review. B./Physical review. B•Nov 29, 2021•1 citation

Determining topological order from infinite projected entangled pair states

Anna Francuz, Jacek Dziarmaga, Guifré Vidal +1 more

We present a method of extracting information about the topological order from the ground state of a strongly correlated two-dimensional system computed with the infinite projected entangled pair state (iPEPS). For topologically ordered systems, the iPEPS wrapped on a torus becomes a superposition of degenerate, locally indistinguishable ground states. Projectors in the form of infinite matrix product operators onto states with well-defined anyon flux are used to compute topological $S$ and $T$ matrices (encoding mutual and self-statistics of emergent anyons). The algorithm is shown to be robust against a perturbation driving string-net toric code across a phase transition to a ferromagnetic phase. Our approach provides accurate results near quantum phase transition, where the correlation length is prohibitively large for other numerical methods. Moreover, we used numerically optimized iPEPS describing the ground state of the Kitaev honeycomb model in the toric code phase and obtained topological data in excellent agreement with theoretical prediction.

arXiv:1910.09661•Physical review. B./Physical review. B•Jan 21, 2020•10 citations

Towards next-generation methods to optimize two-dimensional tensor networks: Algorithmic differentiation and applications to quantum magnets

Juraj Hašík

Infinite Projected entangled-pair state (iPEPS) is a two-dimensional tensor network (TN) ansatz targeting the ground states of lattice models directly in the thermodynamic limit. It is the most successful among the true two-dimensional TN methods, challenging the more common density matrix renormalization group computations conducted on a finite-width cylinders. iPEPS, being a variational method does not suffer from the infamous sign problem, and hence could be readily applied to both fermionic problems off the half-filling and frustrated spin system. Using the imaginary-time optimization, the so called Full Update, led to the celebrated results of iPEPS such as the analysis of stripes in doped t-J and Hubbard models or the insight into the magnetization plateaus in Shastry-Sutherland model. However, the physics of the continuous phase transitions from the symmetry-broken Neel phase into the valence-bond solids or spin-liquids in frustrated magnets is still eluding iPEPS. The main difficulty lies in the optimization, since Full Update is just not precise enough. In the pursuit of the best variational iPEPS states, the gradient methods were put forward in 2016 by Corboz and Vanderstraeten. However, due to their complex nature, they have not been widely adopted especially for iPEPS ansatze with larger unit cells. We will introduce an alternative approach first applied in the context of iPEPS by Liao et al., based on the well established method of Algorithmic differentiation (AD) which allows us to evaluate gradients for iPEPS in a robust and conceptually simple way. Finally, we will present the results of our investigation into the physics of frustrated magnets putting AD to work.

•Oct 25, 2019•3 citations

The local characterization of global tensor network eigenstates

José Garre Rubio, András Molnár, Norbert Schuch +1 more

We study the conditions under which Matrix Product States (MPS) or Matrix Product Operators are exact eigenvectors of an extensive local operator, such as a Hamiltonian. By suitably choosing the local operator, this covers a wide range of settings: Exact eigenstates of Hamiltonians, including scar states, exact MPS trajectories for driven quantum systems, steady states of local Lindbladians, generalized symmetries of either Hamiltonians or density matrices, and many more. Our key result is that a local, fixed-size equation---namely, how a single term in the operator acts on a block of tensors---provides a necessary and sufficient condition for exact solutions. This local characterization allows us to obtain the full space of MPS solutions to all of the aforementioned problems; in particular, we exemplify the power of our results by recovering the quantum group symmetries of the XXZ model. We also discuss applications to numerical algorithms with MPS and the generalization of our results to 2D, i.e., projected entangled pair states (PEPS).

quant-pharXiv:2603.28349•Mar 31, 2026•0 citations

The iPEPS algorithm for ruby and triangle-honeycomb lattices

Saeed S. Jahromi, Román Orús, Mehdi Kargarian +1 more

The infinite Projected Entangled-Pair State (iPEPS) algorithm is one of the most efficient techniques for studying the ground-state properties of two-dimensional quantum lattice Hamiltonians in the thermodynamic limit. Here, we show how the algorithm can be adapted to explore nearest-neighbuor local Hamiltonians on the ruby and triangle-honeycomb lattices, using Corner Transfer Matrix (CTM) renormalization group for 2D tensor network contraction. Additionally, we show how the CTM method can be used to calculate the ground state fidelity per lattice site and the boundary density operator and entanglement entropy (EE) on an infinite cylinder. As a benchmark, we apply the iPEPS method to the ruby model with anisotropic interactions and explore the ground-state properties of the system. We further extract the phase diagram of the model in different regimes of the couplings by measuring two-point correlators, ground sate fidelity and EE on an infinite cylinder. Our phase diagram is in agreement with previous studies of the model by exact diagonalization.

•Nov 13, 2017•0 citations

Infinite projected entangled-pair state algorithm for ruby and triangle-honeycomb lattices

Saeed S. Jahromi, Román Orús, Mehdi Kargarian +1 more

The infinite projected entangled-pair state (iPEPS) algorithm is one of the most efficient techniques for studying the ground-state properties of two-dimensional quantum lattice Hamiltonians in the thermodynamic limit. Here, we show how the algorithm can be adapted to explore nearest-neighbor local Hamiltonians on the ruby and triangle-honeycomb lattices, using the corner transfer matrix (CTM) renormalization group for 2D tensor network contraction. Additionally, we show how the CTM method can be used to calculate the ground-state fidelity per lattice site and the boundary density operator and entanglement entropy (EE) on an infinite cylinder. As a benchmark, we apply the iPEPS method to the ruby model with anisotropic interactions and explore the ground-state properties of the system. We further extract the phase diagram of the model in different regimes of the couplings by measuring two-point correlators, ground-state fidelity, and EE on an infinite cylinder. Our phase diagram is in agreement with previous studies of the model by exact diagonalization.

arXiv:1711.04798•Physical review. B./Physical review. B•Mar 29, 2018•32 citations
Given a tensor network state, how can we determine conserved operators (including Hamiltonians) for which the state is an eigenstate? We answer this question by presenting a method to extract geometrically kkk-local conserved operators that have the given infinite projected entangled pair state (iPEPS) in 2D as an (approximate) eigenstate. The key ingredient is the evaluation of the static structure factors of multi-site operators through differentiating the generating function. Despite the approximation errors, we show that our method is still able to extract from exact or variational iPEPS to good precision both frustration-free and non-frustration-free parent Hamiltonians that are beyond the standard construction and obtain better locality. In particular, we find a 4-site-plaquette local Hamiltonian that approximately has the short-range RVB state as the ground state. Moreover, we find a Hamiltonian that has the deformed toric code state at any string tension as excited eigenstates at the same energy, which might be potential candidates for quantum many-body scars.
We present a method of extracting information about topological order from the ground state of a strongly correlated two-dimensional system represented by an infinite projected entangled pair state (iPEPS). As in previous works [A. Francuz et al., Phys. Rev. B 101, 041108(R) (2020) and A. Francuz and J. Dziarmaga, Phys. Rev. B 102, 235112 (2020)] we begin by determining symmetries of the iPEPS represented by infinite matrix product operators (iMPO) that map between the different iPEPS transfer matrix fixed points, to which we apply the fundamental theorem of matrix product states to find zipper tensors between products of iMPO's that encode fusion properties of the anyons. The zippers can be combined to extract topological FFF symbols of the underlying fusion category, which unequivocally identifies the topological order of the ground state. We bring the FFF symbols to the canonical gauge, and also compute the Drinfeld center of this unitary fusion category to extract the topological SSS and TTT matrices encoding mutual statistics and self-statistics of the emergent anyons. The algorithm is applied to Abelian toric code, Kitaev model, double semion, and twisted quantum double of Z3{Z}_{3}Z3​, as well as to non-Abelian double Fibonacci, double Ising, and quantum double of S3{S}_{3}S3​ and Rep(S3)\mathrm{Rep}({S}_{3})Rep(S3​) string-net models.
We present a method of extracting information about the topological order from the ground state of a strongly correlated two-dimensional system computed with the infinite projected entangled pair state (iPEPS). For topologically ordered systems, the iPEPS wrapped on a torus becomes a superposition of degenerate, locally indistinguishable ground states. Projectors in the form of infinite matrix product operators onto states with well-defined anyon flux are used to compute topological SSS and TTT matrices (encoding mutual and self-statistics of emergent anyons). The algorithm is shown to be robust against a perturbation driving string-net toric code across a phase transition to a ferromagnetic phase. Our approach provides accurate results near quantum phase transition, where the correlation length is prohibitively large for other numerical methods. Moreover, we used numerically optimized iPEPS describing the ground state of the Kitaev honeycomb model in the toric code phase and obtained topological data in excellent agreement with theoretical prediction.
Infinite Projected entangled-pair state (iPEPS) is a two-dimensional tensor network (TN) ansatz targeting the ground states of lattice models directly in the thermodynamic limit. It is the most successful among the true two-dimensional TN methods, challenging the more common density matrix renormalization group computations conducted on a finite-width cylinders. iPEPS, being a variational method does not suffer from the infamous sign problem, and hence could be readily applied to both fermionic problems off the half-filling and frustrated spin system. Using the imaginary-time optimization, the so called Full Update, led to the celebrated results of iPEPS such as the analysis of stripes in doped t-J and Hubbard models or the insight into the magnetization plateaus in Shastry-Sutherland model. However, the physics of the continuous phase transitions from the symmetry-broken Neel phase into the valence-bond solids or spin-liquids in frustrated magnets is still eluding iPEPS. The main difficulty lies in the optimization, since Full Update is just not precise enough. In the pursuit of the best variational iPEPS states, the gradient methods were put forward in 2016 by Corboz and Vanderstraeten. However, due to their complex nature, they have not been widely adopted especially for iPEPS ansatze with larger unit cells. We will introduce an alternative approach first applied in the context of iPEPS by Liao et al., based on the well established method of Algorithmic differentiation (AD) which allows us to evaluate gradients for iPEPS in a robust and conceptually simple way. Finally, we will present the results of our investigation into the physics of frustrated magnets putting AD to work.
We study the conditions under which Matrix Product States (MPS) or Matrix Product Operators are exact eigenvectors of an extensive local operator, such as a Hamiltonian. By suitably choosing the local operator, this covers a wide range of settings: Exact eigenstates of Hamiltonians, including scar states, exact MPS trajectories for driven quantum systems, steady states of local Lindbladians, generalized symmetries of either Hamiltonians or density matrices, and many more. Our key result is that a local, fixed-size equation---namely, how a single term in the operator acts on a block of tensors---provides a necessary and sufficient condition for exact solutions. This local characterization allows us to obtain the full space of MPS solutions to all of the aforementioned problems; in particular, we exemplify the power of our results by recovering the quantum group symmetries of the XXZ model. We also discuss applications to numerical algorithms with MPS and the generalization of our results to 2D, i.e., projected entangled pair states (PEPS).
The infinite Projected Entangled-Pair State (iPEPS) algorithm is one of the most efficient techniques for studying the ground-state properties of two-dimensional quantum lattice Hamiltonians in the thermodynamic limit. Here, we show how the algorithm can be adapted to explore nearest-neighbuor local Hamiltonians on the ruby and triangle-honeycomb lattices, using Corner Transfer Matrix (CTM) renormalization group for 2D tensor network contraction. Additionally, we show how the CTM method can be used to calculate the ground state fidelity per lattice site and the boundary density operator and entanglement entropy (EE) on an infinite cylinder. As a benchmark, we apply the iPEPS method to the ruby model with anisotropic interactions and explore the ground-state properties of the system. We further extract the phase diagram of the model in different regimes of the couplings by measuring two-point correlators, ground sate fidelity and EE on an infinite cylinder. Our phase diagram is in agreement with previous studies of the model by exact diagonalization.
The infinite projected entangled-pair state (iPEPS) algorithm is one of the most efficient techniques for studying the ground-state properties of two-dimensional quantum lattice Hamiltonians in the thermodynamic limit. Here, we show how the algorithm can be adapted to explore nearest-neighbor local Hamiltonians on the ruby and triangle-honeycomb lattices, using the corner transfer matrix (CTM) renormalization group for 2D tensor network contraction. Additionally, we show how the CTM method can be used to calculate the ground-state fidelity per lattice site and the boundary density operator and entanglement entropy (EE) on an infinite cylinder. As a benchmark, we apply the iPEPS method to the ruby model with anisotropic interactions and explore the ground-state properties of the system. We further extract the phase diagram of the model in different regimes of the couplings by measuring two-point correlators, ground-state fidelity, and EE on an infinite cylinder. Our phase diagram is in agreement with previous studies of the model by exact diagonalization.