The exactly solvable interacting vertex models introduced by Alcaraz and Lazo in 2007 generalize the six-vertex model by hard-core interactions that forbid certain vertex pairs at horizontal separations s1s_1 and s2s_2, so that the vertical arrows behave as particles of hard-core size s1+s2+1s_1+s_2+1. Their row-to-row transfer matrices Ts1,s2(μ,ν)T_{s_1,s_2}(\mu,\nu) form an infinite commuting family whose conserved charges include the hard-core (constrained) XXZ chains, but the exact solution proceeds through a matrix product ansatz: because the transfer matrices contain nonlocal operators along the row, the authors could not find a generalized R-matrix placing the models within the quantum inverse scattering method, and posed this as an open problem, together with the construction of analogous interacting vertex models generating the hard-core extensions of other fundamental integrable chains (spin-1 Fateev-Zamolodchikov, Izergin-Korepin, SqU(N)S_qU(N) Sutherland, Perk-Schultz, Bariev). Yang-Baxter-type (RLL) formulations now exist for closely related but distinct objects - local hard-rod-deformed XXZ chains (2021) and the spin-1/2 constrained XXZ quantum chain at exclusion distance 2 (2025) - while the 2025 authors explicitly state that the algebraic counterpart for constrained models is still missing. No generalized R-matrix for the two-dimensional interacting vertex models themselves, and no vertex-model realization for the hard-core extensions of the other listed chains, appears in the literature; the problem remains open.

Background

The six-vertex model is the archetypal exactly solved lattice model: on each vertex of a square lattice, ice-type arrows are arranged so that two arrows point in and two point out, leaving six allowed configurations weighted by parameters a0,a1,b1,b2,c1,c2a_0,a_1,b_1,b_2,c_1,c_2. In the quantum inverse scattering method (QISM) the row-to-row transfer matrix T(u)T(u), which propagates one row of arrow configurations to the next, is written as the trace over an auxiliary space of a product of local R-matrices, matrices acting on one auxiliary and one physical site; the Yang-Baxter equation, an intertwining relation guaranteeing the consistency of factorized three-body scattering, then implies [T(u),T(v)]=0[T(u),T(v)]=0. Expanding T(u)T(u) near a regular point yields an infinite family of commuting conserved charges, one of which is the spin-1/2 XXZ Heisenberg chain - the standard explanation of XXZ integrability and the gateway to the algebraic Bethe ansatz (the QISM construction of exact eigenstates from the R-matrix and an auxiliary space) and to correlation functions via the inverse scattering method.

In the late 1990s a family of constrained, or hard-core, XXZ chains was introduced and solved by coordinate Bethe ansatz (a trial wavefunction with particle coordinates and scattering phases, fitted to the eigenvalue equation). These chains are the XXZ chain conjugated by a global projector that removes every configuration in which two up spins (particles) sit closer than a hard-core size ss: the flip-flop terms remain nearest-neighbor while the Ising coupling acts at range ss, and the model equivalently describes asymmetric diffusion of particles of finite extent on a lattice. The surprise was that the hard-core exclusion does not destroy integrability, but no Yang-Baxter explanation was known. Alcaraz and Lazo (2007) supplied one at the level of two-dimensional classical models. They introduced interacting vertex models generalizing the six-vertex model: in addition to the usual nearest-neighbor vertex weights, the transfer-matrix elements carry factors z1nz1z2nz2z_1^{n_{z_1}}z_2^{n_{z_2}} counting specific vertex pairs at horizontal separations s1s_1 and s2s_2 (integers s1,s2≥0s_1,s_2\ge 0; both parameters nonzero give interacting four-vertex models, exactly one parameter zero gives interacting five-vertex models, and s1=s2=0s_1=s_2=0 recovers the standard six-vertex model). The vertical arrows thus behave as particles of hard-core size s1+s2+1s_1+s_2+1. Diagonalizing Ts1,s2T_{s_1,s_2} with a matrix product ansatz - an exact eigenstate ansatz in which amplitudes are products of auxiliary matrices obeying quadratic algebraic relations - they found the eigenfunctions to be independent of two of the four model parameters (γ,ν,μ,h)(\gamma,\nu,\mu,h), which implies a commuting family [Ts1,s2(μ,ν),Ts1,s2(μ′,ν′)]=0[T_{s_1,s_2}(\mu,\nu),T_{s_1,s_2}(\mu',\nu')]=0, and [Ts1,s2,Ts1′,s2′]=0[T_{s_1,s_2},T_{s'_1,s'_2}]=0 whenever s1+s2=s1′+s2′s_1+s_2=s'_1+s'_2. Expanding in μ\mu produces conserved charges, among them the hard-core XXZ Hamiltonian, so the 2007 paper explains why hard-core exclusion preserves XXZ integrability; the spectra of Ts1,s2T_{s_1,s_2} reduce to those of the six-vertex model on a cylinder of reduced perimeter with a momentum-dependent seam.

The construction nevertheless sits outside the QISM: the hard-core factors make Ts1,s2T_{s_1,s_2} a nonlocal operator along the row, whereas the Yang-Baxter framework builds transfer matrices from strictly local building blocks. The authors state they were unable to find a generalized R-matrix reproducing their models and pose this explicitly as a question for the future. They pose a second, closely linked question: several other fundamental integrable chains also admit hard-core generalizations solved by coordinate Bethe ansatz - the spin-1 Fateev-Zamolodchikov chain (a three-state chain whose critical behavior is described by the Z3Z_3 parafermion theory), the Izergin-Korepin chain (a nineteen-vertex spin-1 chain), the SqU(N)S_qU(N) Sutherland and Perk-Schultz chains (higher-rank or multi-species deformations of the XXZ chain), and the Bariev model (a spin chain with two- and three-site interactions) - and the authors expect each to be generated by a larger integrable interacting vertex model of the same type, to be constructed.

Adjacent recent progress, on distinct objects, shapes what a solution could look like. The medium-range integrability theory of Gombor and Pozsgay (2021) develops multi-site Lax operators - spectral-parameter-dependent matrices acting on an auxiliary space and several physical sites - satisfying RLL intertwining relations with ordinary R-matrices, thereby extending the QISM logic to chains whose Hamiltonians act on three or more consecutive sites. Pozsgay, Gombor and Hutsalyuk (2021) used it to construct new local hard-rod-deformed XXZ chains with explicit three-site Lax operators and an explicit R-matrix; the constrained XXZ chain arises there as a projected single-rod-length sector, and the Alcaraz-Lazo paper is cited as the earlier vertex-model explanation of that chain's integrability. Corcoran, de Leeuw and Pozsgay (2025) extended the RLL formalism to constrained Hilbert spaces, classified integrable spin-1/2 chains under the Rydberg blockade constraint (exclusion distance 2), and for the constrained XXZ model constructed an explicit Lax operator and solved the projected RLL relation for an R-matrix, establishing commutativity of the transfer matrices; they simultaneously state that the algebraic Bethe ansatz counterpart for constrained models is still missing, and they treat neither higher exclusion distances nor higher-rank chains nor the two-dimensional interacting vertex models. The question posed by Alcaraz and Lazo - a generalized R-matrix for the interacting vertex models themselves, and vertex-model realizations of the other hard-core chain families - therefore remains without an answer in the literature to date.

Problem Statement

The problem, as posed in the conclusions of Alcaraz and Lazo (J. Stat. Mech. P08008 (2007)), consists of two linked components with a single scientific objective: placing hard-core-deformed integrable lattice systems within the Yang-Baxter / quantum inverse scattering framework, as the six-vertex R-matrix does for the ordinary XXZ chain. A complete solution requires both components. Within component (2) the listed chains are independent instances of the same construction, so a realization for any single one of them is a partial result rather than a complete answer.

(1) Construct a generalized R-matrix formulation of the interacting vertex models: find an R-matrix, or an equivalent local Lax-operator/intertwining structure (for example multi-site Lax operators satisfying an RLL relation, possibly with projected or enlarged auxiliary spaces), whose transfer-matrix construction reproduces the commuting family of row-to-row transfer matrices Ts1,s2(μ,ν)T_{s_1,s_2}(\mu,\nu) of the Alcaraz-Lazo interacting vertex models for arbitrary integers s1,s2≥0s_1,s_2\ge 0, with the six-vertex R-matrix recovered at s1=s2=0s_1=s_2=0, and thereby solve the models within the quantum inverse scattering method. The building blocks must be of bounded, system-size-independent dimension: re-encoding the nonlocal transfer matrix as an LL-dependent object does not answer the question.

(2) Construct, by the same mechanism, interacting vertex models (or their generalized-R-matrix equivalents) whose row-to-row transfer matrices generate the hard-core extensions of the other known integrable quantum chains: the spin-1 Fateev-Zamolodchikov chain, the Izergin-Korepin chain, the SqU(N)S_qU(N) Sutherland and Perk-Schultz chains, and the Bariev model - realizing the authors' expectation that the known nearest-neighbor vertex models are special cases of larger integrable interacting families.

Both components ask for explicit, verifiable Yang-Baxter-type objects. The spectra of all these models are already known by other methods (matrix product ansatz, coordinate Bethe ansatz); what is requested is the algebraic, R-matrix-level explanation of that solvability, not a recomputation of eigenvalues.

References

  1. F. C. Alcaraz and M. J. Lazo, Exactly solvable interacting vertex models, J. Stat. Mech. P08008 (2007). DOI: 10.1088/1742-5468/2007/08/P08008; arXiv:0706.2178. https://arxiv.org/abs/0706.2178
  2. F. C. Alcaraz and R. Z. Bariev, An Exactly Solvable Constrained XXZ Chain, in Statistical Physics on the Eve of the Twenty-First Century, Series on Advances in Statistical Mechanics Vol. 14, edited by M. T. Batchelor and L. T. Wille (World Scientific, Singapore, 1999). arXiv:cond-mat/9904042. https://arxiv.org/abs/cond-mat/9904042
  3. F. C. Alcaraz and R. Z. Bariev, Integrable models of strongly correlated particles with correlated hopping, Phys. Rev. B 59, 3373 (1999). DOI: 10.1103/PhysRevB.59.3373; preprint arXiv:cond-mat/9904129. https://arxiv.org/abs/cond-mat/9904129
  4. R. Z. Bariev, Integrable spin chain with two- and three-particle interactions, J. Phys. A 24, L549 (1991). DOI: 10.1088/0305-4470/24/10/010. https://doi.org/10.1088/0305-4470/24/10/010
  5. B. Pozsgay, T. Gombor and A. Hutsalyuk, Integrable hard-rod deformation of the Heisenberg spin chains, Phys. Rev. E 104, 064124 (2021). DOI: 10.1103/PhysRevE.104.064124; arXiv:2108.13724. https://arxiv.org/abs/2108.13724
  6. T. Gombor and B. Pozsgay, Integrable spin chains and cellular automata with medium-range interaction, Phys. Rev. E 104, 054123 (2021). DOI: 10.1103/PhysRevE.104.054123; arXiv:2108.02053. https://arxiv.org/abs/2108.02053
  7. L. Corcoran, M. de Leeuw and B. Pozsgay, Integrable models on Rydberg atom chains, SciPost Phys. 18, 139 (2025). DOI: 10.21468/SciPostPhys.18.4.139; arXiv:2405.15848. https://arxiv.org/abs/2405.15848
  8. B. Pozsgay, T. Gombor, A. Hutsalyuk, Y. Jiang, L. Pristyak and E. Vernier, Integrable spin chain with Hilbert space fragmentation and solvable real-time dynamics, Phys. Rev. E 104, 044106 (2021). DOI: 10.1103/PhysRevE.104.044106; arXiv:2105.02252. https://arxiv.org/abs/2105.02252
  9. M. Borsi, L. Pristyak and B. Pozsgay, Matrix Product Symmetries and Breakdown of Thermalization from Hard Rod Deformations, Phys. Rev. Lett. 131, 037101 (2023). DOI: 10.1103/PhysRevLett.131.037101; arXiv:2302.07219. https://arxiv.org/abs/2302.07219
  10. T. Yang, F.-K. Wen, K. Hao, L.-K. Cao and R.-H. Yue, The Effect of a Long-Range Correlated-Hopping Interaction on Bariev Spin Chains, Entropy 17, 6044 (2015). DOI: 10.3390/e17096044. https://www.mdpi.com/1099-4300/17/9/6044
  11. Y. Jiang, F. Loebbert and D. Zhong, Irrelevant deformations with boundaries and defects, J. Stat. Mech. 2022, 043102 (2022). DOI: 10.1088/1742-5468/ac6251; arXiv:2109.13180. https://arxiv.org/abs/2109.13180
  12. M. Blakeney, L. Corcoran and M. de Leeuw, Constrained integrability and anyonic chains, Int. J. Mod. Phys. A (2026). DOI: 10.1142/S0217751X26300140; arXiv:2605.28946. https://arxiv.org/abs/2605.28946