[Paper Review] Two-parameter families of matrix product operator integrals of motion in Heisenberg spin chains

V. I. Yashin|arXiv (Cornell University)|Feb 23, 2026

Quantum many-body systems0 citations

TL;DR

The author discovers two-parameter families of matrix product operator (MPO) integrals of motion that commute with the XXX, XXZ, and XYZ Heisenberg spin chains, using a symbolic-algebra approach and bond-dimension-4 MPOs. Expanding these MPOs yields local charges and suggests broader integrability insights.

ABSTRACT

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.

Motivation & Objective

• Motivate and extend the study of integrals of motion in one-dimensional Heisenberg spin chains.
• Provide explicit two-parameter MPO charges that commute with XXX, XXZ, and XYZ Hamiltonians.
• Describe a symbolic-algebra approach to constructing such MPOs and explore potential applications.

Proposed method

• Construct two-parameter families of MPOs with bond dimension 4 that commute with the Heisenberg Hamiltonians under a simple sufficient condition for MPO charges.
• Use an error-term framework (E_j) to ensure [H, M] = 0 across sites.
• Represent MPOs with four-by-four operator matrices A_j whose entries depend on local spin operators.
• Expand the MPO charges as a series in the parameters to recover local charges of the models.
• Highlight a geometry where the parameters lie on a sphere for XXX and XXZ cases.
• Discuss how this approach avoids reliance on special functions and relates to known integrability constructions.

Experimental results

Research questions

• RQ1 Do two-parameter MPOs commuting with H exist for XXX, XXZ, and XYZ spin chains?
• RQ2 What is the structural form (bond dimension) and parameter geometry of these MPO charges?
• RQ3 Can expanding the two-parameter MPOs generate all local charges of the models?
• RQ4 What potential applications and extensions does this MPO framework enable beyond prior one-parameter results?

Key findings

• Existence of two-parameter families of MPO integrals of motion commuting with XXX, XXZ, and XYZ Hamiltonians with bond dimension 4.
• Expanding the MPO charges in the parameters yields local charges of the spin chains.
• The parameter space has convenient geometry, e.g., XXX and XXZ solutions depend on points on a sphere.
• The MPO construction provides a straightforward, function-free description compared to Baxter’s eight-vertex solutions, and connects to transfer-matrix ideas.
• The approach uses symbolic algebra to discover these charges and may extend to broader integrability contexts.