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:1407.4262 · Journal of Statistical Mechanics Theory and Experiment

On integrable matrix product operators with bond dimensionD= 4

Hosho Katsura

Published Jan 6, 2015• 7 citations
arXivHTMLPDF

Abstract

We construct and study a two-parameter family of matrix product operators of bond dimension $D=4$. The operators $M(x,y)$ act on $({\mathbb C}_2)^{\otimes N}$, i.e., the space of states of a spin-$1/2$ chain of length $N$. For the particular values of the parameters: $x=1/3$ and $y=1/\sqrt{3}$, the operator turns out to be proportional to the square root of the reduced density matrix of the valence-bond-solid state on a hexagonal ladder. We show that $M(x,y)$ has several interesting properties when $(x,y)$ lies on the unit circle centered at the origin: $x^2 + y^2=1$. In this case, we find that $M(x,y)$ commutes with the Hamiltonian and all the conserved charges of the isotropic spin-$1/2$ Heisenberg chain. Moreover, $M(x_1,y_1)$ and $M(x_2,y_2)$ are mutually commuting if $x^2_i + y^2_i=1$ for both $i=1$ and $2$. These remarkable properties of $M(x,y)$ are proved as a consequence of the Yang-Baxter equation.

DOI
10.1088/1742-5468/2015/01/p01006

Related papers

Two-parameter families of matrix product operator integrals of motion in Heisenberg spin chains

Vsevolod I. Yashin

Recently, Fendley et al. (2025) [arXiv:2511.04674] revealed a new simple way to demonstrate the integrability of XYZ Heisenberg model by constructing a one-parameter family of integrals of motion in the matrix product operator (MPO) form with bond dimension 4. In this work, I report on the discovery of two-parameter families of MPOs that commute with Heisenberg spin chain Hamiltonian in case of various anisotropies (XXX, XXZ, XX, XY and XYZ). These solutions are connected by taking appropriate limits. For all cases except XYZ, I also write down Floquet charges of two-step Floquet protocols corresponding to the Trotterization. I describe a symbolic algebra approach for finding such integrals of motion and speculate about possible generalizations and applications.

cond-mat.stat-mechCross-listedarXiv:2602.19741•Feb 24, 2026•0 citations

Matrix product operator representations for the local conserved quantities of the Heisenberg chain

Kyoichi Yamada, Kohei Fukai

We present the explicit expressions for the matrix product operator (MPO) representation for the local conserved quantities of the Heisenberg chain. The bond dimension of the MPO grows linearly with the locality of the charges. The MPO has more simple form than the local charges themselves, and their Catalan tree patterns naturally emerge from the matrix products. The MPO representation of local conserved quantities is generalized to the integrable SU(N) invariant spin chain.

arXiv:2306.03431•SciPost Physics Core•Oct 19, 2023•5 citations

Matrix-product operator dualities in integrable lattice models

Yuan Miao, Andras Molnar, Nick G. Jones

Matrix-product operators (MPOs) appear throughout the study of integrable lattice models, notably as the transfer matrices. They can also be used as transformations to construct dualities between such models, both invertible (including unitary) and non-invertible (including discrete gauging). We analyse how the local Yang--Baxter integrable structures are modified under such dualities. We see that the $\check{R}$-matrix, that appears in the baxterization approach to integrability, transforms in a simple manner. We further argue for a broad class of MPO dualities that the usual Yang--Baxter $R$-matrix should be extended. This extended $R$-matrix satisfies modified algebraic relations, including a modified Yang--Baxter algebra previously identified in the unitary case, that altogether give a local integrable structure underlying the commuting transfer matrices of the dual model. We illustrate these results with two case studies, analysing an invertible unitary MPO and a non-invertible MPO both applied to the canonical XXZ spin chain. The former is the cluster entangler, arising in the study of symmetry-protected topological phases, while the latter is the Kramers--Wannier duality. We show several results for MPOs with exact MPO inverses that are of independent interest.

cond-mat.stat-mechCross-listedarXiv:2602.17436•Journal of Physics A Mathematical and Theoretical•Feb 20, 2026•0 citations

Integrable Matrix Product States from boundary integrability

Balázs Pozsgay, Piroli, L., Vernier, E.

We consider integrable Matrix Product States (MPS) in integrable spin chains and show that they correspond to "operator valued" solutions of the so-called twisted Boundary Yang-Baxter (or reflection) equation. We argue that the integrability condition is equivalent to a new linear intertwiner relation, which we call the "square root relation", because it involves half of the steps of the reflection equation. It is then shown that the square root relation leads to the full Boundary Yang-Baxter equations. We provide explicit solutions in a number of cases characterized by special symmetries. These correspond to the "symmetric pairs" $(SU(N),SO(N))$ and $(SO(N),SO(D)\otimes SO(N-D))$, where in each pair the first and second elements are the symmetry groups of the spin chain and the integrable state, respectively. These solutions can be considered as explicit representations of the corresponding twisted Yangians, that are new in a number of cases. Examples include certain concrete MPS relevant for the computation of one-point functions in defect AdS/CFT.

•May 27, 2019•50 citations

Integrable Matrix Product States from boundary integrability

Balázs Pozsgay, Lorenzo Piroli, Eric Vernier

We consider integrable Matrix Product States (MPS) in integrable spin chains and show that they correspond to “operator valued” solutions of the so-called twisted Boundary Yang-Baxter (or reflection) equation. We argue that the integrability condition is equivalent to a new linear intertwiner relation, which we call the “square root relation”, because it involves half of the steps of the reflection equation. It is then shown that the square root relation leads to the full Boundary Yang-Baxter equations. We provide explicit solutions in a number of cases characterized by special symmetries. These correspond to the “symmetric pairs” (SU(N),SO(N)) and (SO(N),SO(D) \otimes ⊗ SO(N-D)), where in each pair the first and second elements are the symmetry groups of the spin chain and the integrable state, respectively. These solutions can be considered as explicit representations of the corresponding twisted Yangians, that are new in a number of cases. Examples include certain concrete MPS relevant for the computation of one-point functions in defect AdS/CFT.

arXiv:1812.11094•SciPost Physics•May 27, 2019•49 citations

Matrix product solutions of boundary driven quantum chains

Tomaž Prosen

We review recent progress on constructing non-equilibrium steady state density operators of boundary driven locally interacting quantum chains, where driving is implemented via Markovian dissipation channels attached to the chain's ends. We discuss explicit solutions in three different classes of quantum chains, specifically, the paradigmatic (anisotropic) Heisenberg spin- chain, the Fermi–Hubbard chain, and the Lai–Sutherland spin-1 chain, and discuss universal concepts which characterize these solutions, such as matrix product ansatz and a more structured walking graph state ansatz. The central theme is the connection between the matrix product form of nonequilibrium states and the integrability structures of the bulk Hamiltonian, such as the Lax operators and the Yang–Baxter equation. However, there is a remarkable distinction with respect to the conventional quantum inverse scattering method, namely addressing nonequilibrium steady state density operators requires non-unitary irreducible representations of Yang–Baxter algebra which are typically of infinite dimensionality. Such constructions result in non-Hermitian, and often also non-diagonalisable families of commuting transfer operators which in turn result in novel conservation laws of the integrable bulk Hamiltonians. For example, in the case of the anisotropic Heisenberg model, quasi-local conserved operators which are odd under spin reversal (or spin flip) can be constructed, whereas the conserved operators stemming from orthodox Hermitian transfer operators (via logarithmic differentiation) are all even under spin reversal.

arXiv:1504.00783•Journal of Physics A Mathematical and Theoretical•Aug 19, 2015•136 citations
We construct and study a two-parameter family of matrix product operators of bond dimension D=4D=4D=4. The operators M(x,y)M(x,y)M(x,y) act on (C2)⊗N({\mathbb C}_2)^{\otimes N}(C2​)⊗N, i.e., the space of states of a spin-1/21/21/2 chain of length NNN. For the particular values of the parameters: x=1/3x=1/3x=1/3 and y=1/3y=1/\sqrt{3}y=1/3​, the operator turns out to be proportional to the square root of the reduced density matrix of the valence-bond-solid state on a hexagonal ladder. We show that M(x,y)M(x,y)M(x,y) has several interesting properties when (x,y)(x,y)(x,y) lies on the unit circle centered at the origin: x2+y2=1x^2 + y^2=1x2+y2=1. In this case, we find that M(x,y)M(x,y)M(x,y) commutes with the Hamiltonian and all the conserved charges of the isotropic spin-1/21/21/2 Heisenberg chain. Moreover, M(x1,y1)M(x_1,y_1)M(x1​,y1​) and M(x2,y2)M(x_2,y_2)M(x2​,y2​) are mutually commuting if xi2+yi2=1x^2_i + y^2_i=1xi2​+yi2​=1 for both i=1i=1i=1 and 222. These remarkable properties of M(x,y)M(x,y)M(x,y) are proved as a consequence of the Yang-Baxter equation.
Recently, Fendley et al. (2025) [arXiv:2511.04674] revealed a new simple way to demonstrate the integrability of XYZ Heisenberg model by constructing a one-parameter family of integrals of motion in the matrix product operator (MPO) form with bond dimension 4. In this work, I report on the discovery of two-parameter families of MPOs that commute with Heisenberg spin chain Hamiltonian in case of various anisotropies (XXX, XXZ, XX, XY and XYZ). These solutions are connected by taking appropriate limits. For all cases except XYZ, I also write down Floquet charges of two-step Floquet protocols corresponding to the Trotterization. I describe a symbolic algebra approach for finding such integrals of motion and speculate about possible generalizations and applications.
Matrix-product operators (MPOs) appear throughout the study of integrable lattice models, notably as the transfer matrices. They can also be used as transformations to construct dualities between such models, both invertible (including unitary) and non-invertible (including discrete gauging). We analyse how the local Yang--Baxter integrable structures are modified under such dualities. We see that the Rˇ\check{R}Rˇ-matrix, that appears in the baxterization approach to integrability, transforms in a simple manner. We further argue for a broad class of MPO dualities that the usual Yang--Baxter RRR-matrix should be extended. This extended RRR-matrix satisfies modified algebraic relations, including a modified Yang--Baxter algebra previously identified in the unitary case, that altogether give a local integrable structure underlying the commuting transfer matrices of the dual model. We illustrate these results with two case studies, analysing an invertible unitary MPO and a non-invertible MPO both applied to the canonical XXZ spin chain. The former is the cluster entangler, arising in the study of symmetry-protected topological phases, while the latter is the Kramers--Wannier duality. We show several results for MPOs with exact MPO inverses that are of independent interest.
We consider integrable Matrix Product States (MPS) in integrable spin chains and show that they correspond to "operator valued" solutions of the so-called twisted Boundary Yang-Baxter (or reflection) equation. We argue that the integrability condition is equivalent to a new linear intertwiner relation, which we call the "square root relation", because it involves half of the steps of the reflection equation. It is then shown that the square root relation leads to the full Boundary Yang-Baxter equations. We provide explicit solutions in a number of cases characterized by special symmetries. These correspond to the "symmetric pairs" (SU(N),SO(N))(SU(N),SO(N))(SU(N),SO(N)) and (SO(N),SO(D)⊗SO(N−D))(SO(N),SO(D)\otimes SO(N-D))(SO(N),SO(D)⊗SO(N−D)), where in each pair the first and second elements are the symmetry groups of the spin chain and the integrable state, respectively. These solutions can be considered as explicit representations of the corresponding twisted Yangians, that are new in a number of cases. Examples include certain concrete MPS relevant for the computation of one-point functions in defect AdS/CFT.
We consider integrable Matrix Product States (MPS) in integrable spin chains and show that they correspond to “operator valued” solutions of the so-called twisted Boundary Yang-Baxter (or reflection) equation. We argue that the integrability condition is equivalent to a new linear intertwiner relation, which we call the “square root relation”, because it involves half of the steps of the reflection equation. It is then shown that the square root relation leads to the full Boundary Yang-Baxter equations. We provide explicit solutions in a number of cases characterized by special symmetries. These correspond to the “symmetric pairs” (SU(N),SO(N)) and (SO(N),SO(D) \otimes ⊗ SO(N-D)), where in each pair the first and second elements are the symmetry groups of the spin chain and the integrable state, respectively. These solutions can be considered as explicit representations of the corresponding twisted Yangians, that are new in a number of cases. Examples include certain concrete MPS relevant for the computation of one-point functions in defect AdS/CFT.
We review recent progress on constructing non-equilibrium steady state density operators of boundary driven locally interacting quantum chains, where driving is implemented via Markovian dissipation channels attached to the chain's ends. We discuss explicit solutions in three different classes of quantum chains, specifically, the paradigmatic (anisotropic) Heisenberg spin- chain, the Fermi–Hubbard chain, and the Lai–Sutherland spin-1 chain, and discuss universal concepts which characterize these solutions, such as matrix product ansatz and a more structured walking graph state ansatz. The central theme is the connection between the matrix product form of nonequilibrium states and the integrability structures of the bulk Hamiltonian, such as the Lax operators and the Yang–Baxter equation. However, there is a remarkable distinction with respect to the conventional quantum inverse scattering method, namely addressing nonequilibrium steady state density operators requires non-unitary irreducible representations of Yang–Baxter algebra which are typically of infinite dimensionality. Such constructions result in non-Hermitian, and often also non-diagonalisable families of commuting transfer operators which in turn result in novel conservation laws of the integrable bulk Hamiltonians. For example, in the case of the anisotropic Heisenberg model, quasi-local conserved operators which are odd under spin reversal (or spin flip) can be constructed, whereas the conserved operators stemming from orthodox Hermitian transfer operators (via logarithmic differentiation) are all even under spin reversal.