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• Revised Apr 14, 2026 (v2)
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. These generating functions define a manifold of the given tensor network state deformed by some parameters, endowed with a quantum geometry, where conserved operators correspond to vanishing fidelity susceptibility. 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 for which the deformed toric code state at arbitrary string tension is an excited eigenstate with the same energy, thereby potentially realizing 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

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

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

Determining non-Abelian topological order from infinite projected entangled pair states

Anna Francuz, Jacek Dziarmaga

The iPEPS tensor network, which is free of the sign problem notorious for the quantum Monte Carlo, is the method of choice to search for topological order in two-dimensional strongly correlated systems. However, finding a ground state represented by iPEPS is one problem, while extracting from this representation the desired information about topological order is quite another. The algorithm, proposed by the authors, does precisely this in a wide class of non-Abelian string net models.

arXiv:2008.06391•Physical review. B./Physical review. B•Dec 4, 2020•10 citations

Normal projected entangled pair states generating the same state

Molnár, A., Garre-Rubio, J., Pérez-García, D. +2 more

Tensor networks (TNs) are generated by a set of small rank tensors and define many-body quantum states in a succinct form. The corresponding map is not one-to-one: different sets of tensors may generate the very same state. A fundamental question in the study of TNs naturally arises: what is then the relation between those sets? The answer to this question in one-dimensional setups has found several applications, like the characterization of local and global symmetries, the classification of phases of matter and unitary evolutions, or the determination of the fixed points of renormalization procedures. Here we answer this question for projected entangled pair states in any dimension and lattice geometry (including, for example, the Kagome lattice, hyperbolic lattices, or tree tensor networks), as long as the tensors generating the states are normal, which constitute an important and generic class.

•Nov 12, 2018•31 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. These generating functions define a manifold of the given tensor network state deformed by some parameters, endowed with a quantum geometry, where conserved operators correspond to vanishing fidelity susceptibility. 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 for which the deformed toric code state at arbitrary string tension is an excited eigenstate with the same energy, thereby potentially realizing 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 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).
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.
Tensor networks (TNs) are generated by a set of small rank tensors and define many-body quantum states in a succinct form. The corresponding map is not one-to-one: different sets of tensors may generate the very same state. A fundamental question in the study of TNs naturally arises: what is then the relation between those sets? The answer to this question in one-dimensional setups has found several applications, like the characterization of local and global symmetries, the classification of phases of matter and unitary evolutions, or the determination of the fixed points of renormalization procedures. Here we answer this question for projected entangled pair states in any dimension and lattice geometry (including, for example, the Kagome lattice, hyperbolic lattices, or tree tensor networks), as long as the tensors generating the states are normal, which constitute an important and generic class.