Chiral Integrable Matrix Product States

• Chiral integrable matrix product states are exactly solvable many-body states that break spatial parity while preserving integrability.
• They leverage algebraic frameworks like the Onsager algebra and reflection equations to construct explicit ground and excited state representations.
• Their construction employs systematic fusion and dressing techniques, offering insights into spectral properties, overlaps, and nonlocal order parameters.

Chiral integrable matrix product states (MPS) are exactly solvable many-body states in quantum spin chains or vertex models, characterized by their compatibility with integrable parent Hamiltonians and their manifest breaking of spatial parity (chirality). These states constitute an essential class for understanding the interplay between integrability, topological order, and exact solvability of strongly correlated systems. Integrable MPS provide explicit, tractable representatives of ground and excited states in diverse models, ranging from the Onsager-integrable chiral clock chain to open spin-chains with twisted Yangian symmetry and alternating fundamental/anti-fundamental structures. Chirality is embedded either by spectral asymmetry, nontrivial site alternation, or boundary/defect construction, and integrability is ensured by algebraic structures such as the Onsager algebra, reflection equations, or fused transfer-matrix identities.

1. Algebraic Frameworks: Onsager Algebra, Reflection Equations, and Chiral Integrability

Integrable chiral MPS emerge from the confluence of algebraic symmetries constraining the Hamiltonians. In the $N$-state chiral clock chain, the local Hilbert space carries $N$-dimensional clock operators $X_k, Z_k$ with $X_k^N = Z_k^N = 1$ and $X_k Z_l = \omega^{\delta_{kl}} Z_l X_k$, $\omega = e^{2\pi i/N}$. The Hamiltonian is expressed as $H_A = \sum_m t_m A_m$, where the $A_m$ are fixed linear combinations of on-site and nearest-neighbor clock operator monomials. These operators generate the Onsager algebra, characterized by the Dolan–Grady relations and an infinite set of nontrivial commutator identities:

$$[A_l, A_m] = G_{l-m}, \quad [G_l, A_m] = A_{m+l} - A_{m-l}/2, \quad [G_l, G_m] = 0.$$

Integrability in open boundary spin chains and alternating chains is implemented via reflection equations (RE), both soliton-preserving (SP) and soliton-non-preserving (SNP), for the boundary $K$-matrices:

$$R_{12}(u-v) K_1(u) R_{21}(u+v) K_2(v) = K_2(v) R_{12}(u+v) K_1(u) R_{21}(u-v),$$

and their fused and twisted variants. The untwisted (chiral) integrability of a boundary MPS $|\mathcal{B}\rangle$ is checked by the transfer-matrix identity:

$$\tau(u) |\mathcal{B}\rangle = \bar\tau(-u-2) |\mathcal{B}\rangle,$$

with $\tau(u), \bar\tau(u)$ the fundamental and 'dual' transfer matrices and $u$ the spectral parameter, ensuring chiral symmetry and integrability (Camp et al., 10 Nov 2025, Liu et al., 2 Feb 2026, Liu et al., 4 Jul 2025).

2. MPS Skeletons and Density in Chiral Integrable Families

The MPS skeleton is the dense subset of parent Hamiltonians in the chiral integrable family for which exact finite-bond-dimension MPS ground states are known. In the $N$-state chiral clock chain, this subfamily is defined by requiring the Laurent polynomial $f(z) = \sum_m t_m z^m$ to factor as $f(z) = \pm z^p g(z)^2$, with $g(z)$ an ordinary polynomial and $p \in \mathbb{Z}$. For $N=2$, this condition is necessary and sufficient for an exact MPS ground state; for $N>2$, it is sufficient but may not be exhaustive. Any Hamiltonian $H_A$ in a gapped Onsager sector can be arbitrarily well approximated, in the thermodynamic limit, by skeleton Hamiltonians $H_A^{(D)}$ with $f_D(z)$ a truncation of $g(z)^2$. Thus, the exact-MPS skeleton yields a systematically improvable analytic approximation to all ground states in the gapped phase diagram (Camp et al., 10 Nov 2025).

For alternating $SU(4)$ spin chains and ABJM spin chains, the chiral integrable subspace is constructed by imposing the untwisted transfer-matrix condition and is shown—numerically and algebraically—to possess a nontrivial but constrained set of simple and fused product states (i.e., MPS with block or super-tensors) (Liu et al., 4 Jul 2025, Liu et al., 2 Feb 2026).

3. Construction of Chiral Integrable Matrix Product States

The explicit MPS construction in integrable chiral models distills the algebraic ingredients into an exact many-body ansatz. For the Onsager chiral clock chain, the skeleton condition $f(z)=\pm z^p g(z)^2$ is leveraged via the Schur–Cohn recursion to compute a set of parameters $\{\beta_k\}$, yielding (Camp et al., 10 Nov 2025):

$$|\phi\rangle = M^{(d)} M^{(d-1)} \cdots M^{(1)} |\psi_p^\pm\rangle, \quad M^{(k)} = \exp(\mp \beta_k A_{p+k}),$$

where $|\psi_p^\pm\rangle$ is the unique ground state of $\pm A_p$, and each $M^{(k)}$ is a local exponentiated MPO, so that the product is a translationally invariant MPS of bond dimension at most $\prod_k O(N^{|p+k|})$.

In alternating and ABJM spin chains, two-site and four-site chiral product states are built as translationally invariant MPS of block size $2n$, with tensors (super-tensors) constructed from fused $K$-matrices obtained by systematic fusion procedures applied to the reflection equations (Liu et al., 2 Feb 2026, Liu et al., 4 Jul 2025). These states typically exhibit chiral pairing of physical indices (alternating fundamental and anti-fundamental representations) and are supported only on certain cyclic or alternating site configurations. High-bond-dimension chiral integrable MPS arise via “dressing” with additional R- or monodromy matrices.

4. Spectral Properties: Exact Eigenstates, Excitations, and Overlaps

For chiral MPS skeleton Hamiltonians in the gapped region, the constructed MPS is an exact ground state. Outside this region, the same state remains an exact eigenstate in a specific symmetry sector but does not necessarily minimize the full spectrum. Beyond the ground state, MPS-form excitations (e.g., single-particle states with definite momentum) are constructed by applying local ladder operators and the same sequence of Onsager pivots as in the ground state construction, generalizing the Bogoliubov excitations known from the free-fermion case (Camp et al., 10 Nov 2025).

Overlaps between chiral integrable MPS and on-shell Bethe eigenstates are expressed via determinant expressions involving Gaudin matrices, imposing strict selection rules: only parity-symmetric Bethe root configurations ($\pm$-paired rapidities) yield nonvanishing overlaps, reflecting the underlying chiral symmetry. In ABJM and $SU(4)$ alternating chains, explicit determinant-type overlap formulas hold for both two-site and fused four-site chiral product states. The structure of these overlaps is dictated by the symmetry (symmetric or antisymmetric nature) of the fused $K$-matrices and additional sectorial constraints (e.g., matching parity and even-oddness of total quantum numbers) (Liu et al., 2 Feb 2026, Liu et al., 4 Jul 2025).

5. Disorder and Order Parameters, Boundary CFT Realizations, and Physical Interpretation

On certain branches of the chiral MPS skeleton (e.g., $d=1$ leg in chiral clock), correlators of nonlocal order parameters can be calculated analytically. For even $N$, the long-distance disorder parameter (half-flux string) admits a closed-form expression in terms of transfer matrices and simplifies to $\sqrt{1 - a^2}$ in the thermodynamic limit, where $a$ is a coefficient parametrizing $f(z)$ on the skeleton branch (Camp et al., 10 Nov 2025).

Open Haldane-Shastry–type spin chains exhibit chiral integrable ground states realized as infinite-dimensional MPS arising from chiral correlators of boundary CFTs of $m-1$ free bosons. These systems display polynomial spectral growth and exponentially enhanced degeneracy, reflecting enlarged Yangian-type symmetries (Basu-Mallick et al., 2015).

6. Numerical and Algebraic Characterization of the Chiral Integrable Subspace

Systematic linear-algebraic and numerical studies of low-length (e.g., $L=2,3$) ABJM and $SU(4)$ chains yield complete characterizations of the chiral integrable subspace. For product-type MPS built from tensor products of rank-1 matrices, all such solutions are fully classified by the reflection equation structure and found to exhaust only a subspace of the full chiral integrable space. Fused structures generated via the reflection equation and fusion—chiral blocks spanning more than two sites—yield additional integrable chiral states not accessible by simple product constructions (Liu et al., 2 Feb 2026).

In the ABJM model, the full chiral integrable subspace goes beyond the reflection/fusion-generated families, demonstrated by dimension counts (e.g., $142$ of $196$ for $L=2$), suggesting the existence of further, yet-unclassified chiral integrable MPS, possibly of higher bond dimension or more intricate algebraic structure (Liu et al., 2 Feb 2026).

7. Generalizations and Universality

Nearly every step in the chiral integrable MPS construction relies solely on properties of the underlying algebra (e.g., Onsager algebra, reflection algebra, Yangian symmetry) and not on the specific representation (local Hilbert space). This universality implies that the dense MPS skeleton and approximation techniques extend to arbitrary finite-range spin chain realizations of the algebra, insofar as a unique ground state for a generator exists. The explicit MPS form depends only on the realization; the spectral and overlap properties are controlled by the algebraic relations (Camp et al., 10 Nov 2025, Basu-Mallick et al., 2015).

A comparable algebraic structure appears in the analysis of one-point functions in defect CFTs (e.g., D3–D5 interface) where chiral integrable MPS encode the overlaps with Bethe eigenstates of the bulk spin chain, and the construction recurs for all $k$-dimensional representation sectors with explicit determinant formulas and recursive transfer-matrix relations (Buhl-Mortensen et al., 2015).

---

Chiral integrable matrix product states thus provide a universal and highly structured class of exactly solvable states, combining matrix product representations with the algebraic machinery of integrability and chirality, and serve as a unifying framework for understanding exact solvability, spectral statistics, and nontrivial order parameters in a broad family of quantum integrable systems (Camp et al., 10 Nov 2025, Liu et al., 2 Feb 2026, Liu et al., 4 Jul 2025, Buhl-Mortensen et al., 2015, Basu-Mallick et al., 2015).