© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2203.07443 · Physical review. B./Physical review. B

Locality optimization for parent Hamiltonians of tensor networks

Giuliano Giudici, J. I. Cirac, Norbert Schuch

Published Jul 7, 2022• 8 citations

Abstract

Tensor network states form a powerful framework for both the analytical and numerical study of strongly correlated phases. Vital to their analytical utility is that they appear as the exact ground states of associated parent Hamiltonians, where canonical proof techniques guarantee a controlled ground space structure. Yet, while those Hamiltonians are local by construction, the known techniques often yield complex Hamiltonians which act on a rather large number of spins. In this paper, we present an algorithm to systematically simplify parent Hamiltonians, breaking them down into any given basis of elementary interaction terms. The underlying optimization problem is a semidefinite program, and thus the optimal solution can be found efficiently. Our method exploits a degree of freedom in the construction of parent Hamiltonians---the excitation spectrum of the local terms---over which it optimizes such as to obtain the best possible approximation. We benchmark our method on the AKLT model and the toric code model, where we show that the canonical parent Hamiltonians (acting on 3 or 4 and 12 sites, respectively) can be broken down to the known optimal two-body and four-body terms. We then apply our method to the paradigmatic resonating valence bond (RVB) model on the kagome lattice. Here, the simplest previously known parent Hamiltonian acts on all the 12 spins on one kagome star. With our optimization algorithm, we obtain a vastly simpler Hamiltonian: we find that the RVB model is the ground state of a parent Hamiltonian whose terms are all products of at most four Heisenberg interactions, and whose range can be further constrained, providing a major improvement over the previously known 12-body Hamiltonian.

DOI: 10.1103/physrevb.106.035109

Related papers

Lower Bounds on Spectral Gaps of Parent Hamiltonians via Tensor Networks

Milán Ádám Rozmán, András Molnár, Norbert Schuch

Proving spectral gaps above the ground space is a central problem in quantum many-body physics. Yet, finding lower bounds on the gap is notoriously difficult. We revisit the martingale method, originally developed by Fannes, Nachtergaele, and Werner [Fannes '92, Nachtergaele '96] to prove the existence of spectral gaps of parent Hamiltonians of Matrix Product States (MPS), and provide improvements to the different steps of the method. Most importantly, we devise a new technique which allows to compute the key quantity in the martingale method -- the overlap of local ground spaces -- exactly and efficiently. This enables a clear improvement of the method, allowing it to outperform other existing techniques to lower bound gaps, which we demonstrate by benchmarking on several models. Remarkably, our -- numerically motivated -- approach at the same time also yields a significantly simplified proof of the fact that parent Hamiltonians of any (well-behaved) MPS are always gapped.

quant-pharXiv:2607.19078•Jul 22, 2026•0 citations

Simple Hamiltonians for Matrix Product State models

Norbert Schuch, András Molnár, David Pérez-Garcı́a

Matrix Product States (MPS) and Tensor Networks provide a general framework for the construction of solvable models. The best-known example is the Affleck-Kennedy-Lieb-Tasaki (AKLT) model, which is the ground state of a 2-body nearest-neighbor parent Hamiltonian. We show that such simple parent Hamiltonians for MPS models are, in fact, much more prevalent than hitherto known: The existence of a single example with a simple Hamiltonian for a given choice of dimensions already implies that any generic MPS with those dimensions possesses an equally simple Hamiltonian. We illustrate our finding by discussing a number of models with nearest-neighbor parent Hamiltonians, which generalize the AKLT model on various levels.

arXiv:2503.10767•Mar 13, 2025•0 citations

Microscopic Spin-1 Parent Hamiltonians for Emergent Valence-Bond Loop Manifolds

Hari Borutta, Yasir Iqbal, Kirill Shtengel

We construct an exact spin-$1$ parent Hamiltonian for constrained valence-bond loop manifolds on the checkerboard and pyrochlore lattices. The Hamiltonian is local, SU(2)- and time-reversal-invariant, and built from positive-semidefinite projectors acting on triangular faces. Each projector removes only the maximally polarized state of a triangle, so the model is frustration-free. Its zero-energy states are generated by an AKLT-like construction in which each spin-$1$ moment is resolved into two virtual spin-$\tfrac12$ degrees of freedom, singlets are formed inside every crossed plaquette or tetrahedron, and the physical spin-$1$ Hilbert space is recovered by projection. The resulting ground states are fully packed singlet-loop states on the corner-sharing lattice. Thus, a loop or dimer constraint, usually introduced as part of an effective Rokhsar--Kivelson description, appears here as the exact zero-energy manifold of a microscopic spin Hamiltonian. We analyze spin correlations within this manifold and show that, for a fixed loop covering, they are determined by loop connectivity. We also project symmetry-allowed perturbations into the ground-state manifold and derive the resulting low-energy pseudospin dynamics. The checkerboard and pyrochlore cases differ sharply. On the pyrochlore lattice, tetrahedral symmetry removes simple local bias terms, and the leading nontrivial next-nearest-neighbour Heisenberg perturbation gives an emergent spin-$\tfrac12$ XY model on the diamond lattice of tetrahedron centers. These results give an exact spin-$1$ microscopic starting point for constrained valence-bond physics in two and three dimensions, and show how loop, dimer, and gauge-theoretic descriptions can be approached from a small-spin, SU(2)-invariant frustrated magnet.

arXiv:2607.12925•Jul 14, 2026•0 citations

Non-Abelian chiral spin liquid on a spin-1 kagome lattice: Truncation of an exact Hamiltonian and numerical optimization

Błażej Jaworowski, Anne E. B. Nielsen

We search for short-range Hamiltonians of finite spin-1 kagome systems, maximizing the overlaps with lattice Moore-Read states. Our starting point is an exact, long-range parent Hamiltonian for such a state on a finite plane, obtained from conformal field theory. A truncation procedure is applied to it, which retains only short-range terms and makes it easy to define the Hamiltonian on a torus. Finally, the remaining coefficients are optimized, to yield maximized overlaps between exact diagonalization results and model ground states. In the best cases, the squared overlaps exceed 0.9 and 0.8 for the three lowest states of 12- and 18-site systems, respectively, suggesting that the obtained Hamiltonians are good parent Hamiltonians for a non-Abelian topological order.

arXiv:2202.09193•Physical review. B./Physical review. B•Sep 19, 2022•4 citations

Optimization at the boundary of the tensor network variety

Matthias Christandl, Fulvio Gesmundo, Daniel Stilck França +1 more

Tensor network states form a variational ansatz class widely used, both analytically and numerically, in the study of quantum many-body systems. It is known that if the underlying graph contains a cycle, e.g., as in projected entangled pair states, then the set of tensor network states of given bond dimension is not closed. Its closure is the tensor network variety. Recent work has shown that states on the boundary of this variety can yield more efficient representations for states of physical interest, but it remained unclear how to systematically find and optimize over such representations. We address this issue by defining an ansatz class of states that includes states at the boundary of the tensor network variety of given bond dimension. We show how to optimize over this class in order to find ground states of local Hamiltonians by only slightly modifying standard algorithms and code for tensor networks. We apply this method to different models and observe favorable energies and runtimes when compared with standard tensor network methods.

arXiv:2006.16963•Physical review. B./Physical review. B•May 18, 2021•14 citations

Parent Hamiltonians for stabilizer quantum many-body scars

Shane Dooley

Quantum many-body scars (QMBS) have attracted considerable interest due to their role in weak ergodicity breaking in many-body systems. We present a general construction that embeds stabilizer states as QMBS of local Hamiltonians. The method relies on a notion of factorizability of Pauli strings on a lattice, which is used to convert stabilizer elements into local, few-body operators that annihilate the stabilizer state. This enables the systematic construction of parent Hamiltonians with zero-energy stabilizer QMBS typically near the middle of the spectrum. The method reproduces several known results in a unified framework, including recent examples of volume-law entangled QMBS, such as the ``rainbow'' QMBS and the entangled antipodal Bell pair state. We also apply the framework to construct examples of stabilizer QMBS with a more complex entanglement structure, such as the cluster state, the toric code state, and a volume-law entangled state we dub the antipodal toric code (ATC) state. Exact diagonalization confirms our results and reveal the stabilizer states as exact eigenstates of their parent Hamiltonian.

quant-pharXiv:2601.10805•Physical Review Letters•Jan 19, 2026•0 citations

Tensor network states form a powerful framework for both the analytical and numerical study of strongly correlated phases. Vital to their analytical utility is that they appear as the exact ground states of associated parent Hamiltonians, where canonical proof techniques guarantee a controlled ground space structure. Yet, while those Hamiltonians are local by construction, the known techniques often yield complex Hamiltonians which act on a rather large number of spins. In this paper, we present an algorithm to systematically simplify parent Hamiltonians, breaking them down into any given basis of elementary interaction terms. The underlying optimization problem is a semidefinite program, and thus the optimal solution can be found efficiently. Our method exploits a degree of freedom in the construction of parent Hamiltonians---the excitation spectrum of the local terms---over which it optimizes such as to obtain the best possible approximation. We benchmark our method on the AKLT model and the toric code model, where we show that the canonical parent Hamiltonians (acting on 3 or 4 and 12 sites, respectively) can be broken down to the known optimal two-body and four-body terms. We then apply our method to the paradigmatic resonating valence bond (RVB) model on the kagome lattice. Here, the simplest previously known parent Hamiltonian acts on all the 12 spins on one kagome star. With our optimization algorithm, we obtain a vastly simpler Hamiltonian: we find that the RVB model is the ground state of a parent Hamiltonian whose terms are all products of at most four Heisenberg interactions, and whose range can be further constrained, providing a major improvement over the previously known 12-body Hamiltonian.

Proving spectral gaps above the ground space is a central problem in quantum many-body physics. Yet, finding lower bounds on the gap is notoriously difficult. We revisit the martingale method, originally developed by Fannes, Nachtergaele, and Werner [Fannes '92, Nachtergaele '96] to prove the existence of spectral gaps of parent Hamiltonians of Matrix Product States (MPS), and provide improvements to the different steps of the method. Most importantly, we devise a new technique which allows to compute the key quantity in the martingale method -- the overlap of local ground spaces -- exactly and efficiently. This enables a clear improvement of the method, allowing it to outperform other existing techniques to lower bound gaps, which we demonstrate by benchmarking on several models. Remarkably, our -- numerically motivated -- approach at the same time also yields a significantly simplified proof of the fact that parent Hamiltonians of any (well-behaved) MPS are always gapped.

Matrix Product States (MPS) and Tensor Networks provide a general framework for the construction of solvable models. The best-known example is the Affleck-Kennedy-Lieb-Tasaki (AKLT) model, which is the ground state of a 2-body nearest-neighbor parent Hamiltonian. We show that such simple parent Hamiltonians for MPS models are, in fact, much more prevalent than hitherto known: The existence of a single example with a simple Hamiltonian for a given choice of dimensions already implies that any generic MPS with those dimensions possesses an equally simple Hamiltonian. We illustrate our finding by discussing a number of models with nearest-neighbor parent Hamiltonians, which generalize the AKLT model on various levels.

We construct an exact spin-

1

parent Hamiltonian for constrained valence-bond loop manifolds on the checkerboard and pyrochlore lattices. The Hamiltonian is local, SU(2)- and time-reversal-invariant, and built from positive-semidefinite projectors acting on triangular faces. Each projector removes only the maximally polarized state of a triangle, so the model is frustration-free. Its zero-energy states are generated by an AKLT-like construction in which each spin-

1

moment is resolved into two virtual spin-

\tfrac12

degrees of freedom, singlets are formed inside every crossed plaquette or tetrahedron, and the physical spin-

1

Hilbert space is recovered by projection. The resulting ground states are fully packed singlet-loop states on the corner-sharing lattice. Thus, a loop or dimer constraint, usually introduced as part of an effective Rokhsar--Kivelson description, appears here as the exact zero-energy manifold of a microscopic spin Hamiltonian. We analyze spin correlations within this manifold and show that, for a fixed loop covering, they are determined by loop connectivity. We also project symmetry-allowed perturbations into the ground-state manifold and derive the resulting low-energy pseudospin dynamics. The checkerboard and pyrochlore cases differ sharply. On the pyrochlore lattice, tetrahedral symmetry removes simple local bias terms, and the leading nontrivial next-nearest-neighbour Heisenberg perturbation gives an emergent spin-

\tfrac12

XY model on the diamond lattice of tetrahedron centers. These results give an exact spin-

1

microscopic starting point for constrained valence-bond physics in two and three dimensions, and show how loop, dimer, and gauge-theoretic descriptions can be approached from a small-spin, SU(2)-invariant frustrated magnet.

We search for short-range Hamiltonians of finite spin-1 kagome systems, maximizing the overlaps with lattice Moore-Read states. Our starting point is an exact, long-range parent Hamiltonian for such a state on a finite plane, obtained from conformal field theory. A truncation procedure is applied to it, which retains only short-range terms and makes it easy to define the Hamiltonian on a torus. Finally, the remaining coefficients are optimized, to yield maximized overlaps between exact diagonalization results and model ground states. In the best cases, the squared overlaps exceed 0.9 and 0.8 for the three lowest states of 12- and 18-site systems, respectively, suggesting that the obtained Hamiltonians are good parent Hamiltonians for a non-Abelian topological order.

Tensor network states form a variational ansatz class widely used, both analytically and numerically, in the study of quantum many-body systems. It is known that if the underlying graph contains a cycle, e.g., as in projected entangled pair states, then the set of tensor network states of given bond dimension is not closed. Its closure is the tensor network variety. Recent work has shown that states on the boundary of this variety can yield more efficient representations for states of physical interest, but it remained unclear how to systematically find and optimize over such representations. We address this issue by defining an ansatz class of states that includes states at the boundary of the tensor network variety of given bond dimension. We show how to optimize over this class in order to find ground states of local Hamiltonians by only slightly modifying standard algorithms and code for tensor networks. We apply this method to different models and observe favorable energies and runtimes when compared with standard tensor network methods.

Quantum many-body scars (QMBS) have attracted considerable interest due to their role in weak ergodicity breaking in many-body systems. We present a general construction that embeds stabilizer states as QMBS of local Hamiltonians. The method relies on a notion of factorizability of Pauli strings on a lattice, which is used to convert stabilizer elements into local, few-body operators that annihilate the stabilizer state. This enables the systematic construction of parent Hamiltonians with zero-energy stabilizer QMBS typically near the middle of the spectrum. The method reproduces several known results in a unified framework, including recent examples of volume-law entangled QMBS, such as the ``rainbow'' QMBS and the entangled antipodal Bell pair state. We also apply the framework to construct examples of stabilizer QMBS with a more complex entanglement structure, such as the cluster state, the toric code state, and a volume-law entangled state we dub the antipodal toric code (ATC) state. Exact diagonalization confirms our results and reveal the stabilizer states as exact eigenstates of their parent Hamiltonian.