Structural open problems for the general two-parameter higher-spin six-vertex K-matrix: compact form, positivity, q-oscillator matrix-product construction, and non-integer spin weights

For the higher-spin stochastic six-vertex model (R-matrices of $U_q(\widehat{sl_2})$ in symmetric tensor representations), Mangazeev and Lu (2019) solved the reflection equation for arbitrary positive integer spin weight J, obtaining the general non-diagonal boundary K-matrix: its elements depend on the spectral parameter and two free parameters (nu, t) and are given by an explicit double-sum formula obtained from a terminating balanced 4φ3 generating function. They posed a bundle of structural questions about this object: whether the double sum can be rewritten as a single basic-hypergeometric sum; whether the matrix elements admit a positive (probabilistic) regime for all J; whether a q-oscillator matrix-product construction reproducing the full two-parameter family exists and extends to higher rank; and whether the construction extends to non-integer spin weights, where crossing-unitarity fails. The questions below are the open issues stated in the primary source. One further conjecture from the same paper, the 4th-degree Phi-function identity behind the degenerate reflection equation, has since been proven by Kolyaskin and Mangazeev (2024) via the star-star relation. The four remaining questions, all about the same K-matrix, remain open; none of a single-sum form, a positivity theorem for the general family, a two-parameter matrix-product realization, or a non-integer-weight extension is known.

Background

The six-vertex model assigns Boltzmann weights to the vertices of a square lattice under an ice-type conservation rule (the number of arrows entering each vertex equals the number leaving). Its higher-spin (higher-weight) generalization allows each edge to carry 0,1,...,I arrows, and is equivalent to evaluating the quantum-group R-matrix of the affine algebra $U_q(\widehat{sl_2})$ in the symmetric tensor representation of weight $I$ (a vector space of dimension $I+1$ , corresponding to spin $I/2$ ). After a gauge transformation these weights satisfy a stochasticity (column-sum) rule, so the same object defines a discrete-time Markov chain of exclusion type - the higher-spin stochastic six-vertex model, a relative of the asymmetric simple exclusion process (ASEP) in the Kardar-Parisi-Zhang universality class - whose systematic stochastic interpretation was given by Kuniba, Mangazeev, Maruyama and Okado. The R-matrix elements are built from a function $\Phi_{q^2}(i|j;x,y)$ - a closed product of $q$ -Pochhammer symbols $(a;q)_n=\prod_{k=0}^{n-1}(1-aq^k)$ and $q$ -binomial coefficients - summed over internal indices fixed by the conservation law.

To put an integrable lattice model on an interval with open ends one uses Sklyanin's extension of the quantum inverse scattering method, which requires, besides an R-matrix satisfying the Yang-Baxter equation (the cubic consistency relation that makes row-to-row transfer matrices commute), a boundary K-matrix solving the reflection equation

S_{I,J}(x/y)\,K_I(x)\,S_{J,I}(y/\bar{x})\,K_J(y)=K_J(y)\,S_{I,J}(x/\bar{y})\,K_I(x)\,S_{J,I}(\bar{y}/\bar{x}),

together with a 'dual' reflection equation built on a crossed R-matrix. For weight $I=J=1$ the general non-diagonal $2\times2$ solution with four complex parameters has been known since de Vega and Gonzalez Ruiz; imposing the stochasticity condition selects the branch relevant for Markov chains coupled to reservoirs, parametrized by $t_+$ , $t_-$ and $\nu$ .

Mangazeev and Lu (2019) solved the reflection equation for the higher-spin R-matrix at arbitrary positive integer weight $J$ . Writing the matrix elements through coefficients $N_{j,l}$ ( $0\le j,l\le J$ ) that depend on the spectral parameter $y$ and two free parameters $(\nu,t)$ (with $t^2=t_+/t_-$ , in terms of the two parameters $t_\pm$ of the spin-1/2 boundary solution), they reduced the problem to a pair of linear recurrence relations in the indices (their eqs. (6.5)-(6.6)), solved the difference equations for a generating function in terms of a terminating balanced $_4\phi_3$ series, and expanded the result into an explicit double-sum formula for $N_{j,l}$ , and hence for $K_J$ (their eq. (7.24)). In the lower- and upper-triangular cases ( $t_+=0$ or $t_-=0$ ) the elements collapse to single $\Phi$ -functions, i.e. closed products of $q$ -Pochhammer symbols and $q$ -binomial coefficients; in the generic case the double sum is the most compact known form. The derivation rests on integer weights: crossing symmetry (and crossing unitarity) of the higher-spin R-matrix, needed for the dual reflection equation, is established only for integer weights, and a numerical check in the source shows that crossing unitarity fails for generic complex weights.

A parallel line, often called 3D integrability, constructs K-matrices as matrix products of local operators acting in an auxiliary $q$ -oscillator algebra (a $q$ -deformed bosonic Fock space whose creation and annihilation operators obey $q$ -commutation relations), tied to solutions of tetrahedron-type equations (the three-dimensional analogues of the Yang-Baxter equation). Kuniba, Pasquier and Kuniba, Okado, Yoneyama obtained such matrix-product solutions associated with coideal subalgebras (subalgebras stable under the coproduct); however, their solutions depend essentially only on the spectral parameter, with no analogue of the pair $(\nu,t)$ . The absence of a $q$ -oscillator construction reproducing the two-parameter family, the lack of a single-sum closed form, the unproven positivity of the general matrix elements, and the obstruction at non-integer weights are the subject of the problem below.

Problem Statement

All four questions below concern one object: the general non-diagonal boundary K-matrix $K_J(y)$ of the higher-spin stochastic six-vertex model found by Mangazeev and Lu - the solution of the reflection equation for the $U_q(\widehat{sl_2})$ higher-spin R-matrix at arbitrary positive integer weight $J$ , whose matrix elements are expressed through coefficients $N_{j,l}$ ( $0\le j,l\le J$ ) depending on the spectral parameter $y$ and two free parameters $\nu$ and $t$ , and are given explicitly by the double-sum formula (7.24) of the cited paper (an expansion of a terminating balanced $_4\phi_3$ generating function; in the triangular cases $t_+=0$ or $t_-=0$ the elements reduce to single $\Phi$ -functions, i.e. closed products of $q$ -Pochhammer symbols and $q$ -binomial coefficients). The common objective is to complete the structural theory of this specific K-matrix. Determine, with proof:

(a) Compact form. Can the double-sum expression for $N_{j,l}$ , and hence the matrix elements of $K_J$ , be rewritten as a single basic hypergeometric sum for general $J$ and generic values of the parameters?

(b) Positivity. For every positive integer $J$ , does there exist an explicit regime of the parameters ( $y$ , $\nu$ , $t$ ; e.g. $0<q<1$ with suitable real ranges) in which all matrix elements of the general two-parameter $K_J$ are non-negative, so that in the stochastic normalization they define legitimate boundary transition probabilities? This is known for $J=1$ and expected for all $J$ on fusion grounds (the fusion procedure builds higher-weight K-matrices by projecting tensor products of fundamental ones), but has not been proved for the general family.

(c) Matrix-product (3D) realization. Find an explicit matrix product of local operators acting in an auxiliary $q$ -oscillator algebra that solves the defining recurrence relations (6.5)-(6.6) of Mangazeev and Lu with arbitrary $\nu$ and $t$ and is equivalent to the double-sum formula (7.24); and determine whether this construction extends to higher rank, i.e. to general (non-triangular) stochastic K-matrices of $U_q(\widehat{sl_n})$ with one extra parameter beyond the spectral one. Known matrix-product constructions cover only the one-parameter (spectral) case and, at higher rank, only triangular solutions.

(d) Non-integer weights. Determine whether the K-matrix construction admits a meaningful extension to non-integer (complex) spin weights $J$ . The recurrence relations and their solutions are well defined for complex $J$ (the triangular solutions of the cited paper already are), but crossing symmetry and crossing-unitarity of the higher-spin R-matrix are established only for integer weights and fail generically for complex weights, so the dual-reflection route underlying the construction of the general family is not available there; the source expects that the solution of the reduced system does not satisfy the full reflection equation for complex weights, and leaves this negative expectation unproved. Determine whether the general two-parameter $K_J$ nevertheless solves the reflection equation at non-integer weights (or prove that it cannot).

Each facet is settled either affirmatively (explicit formula, regime, construction, or extension, with proof) or negatively (proof of impossibility or of failure: for facet (b) a proof that for some fixed $J$ no regime of the parameters yields all non-negative matrix elements, and for facet (d) an explicit non-integer weight at which the reflection equation provably fails). The problem as a whole is completed when all four facets are settled for the general two-parameter family, without restriction to triangular specializations, fixed small $J$ , or particular answer forms; facet (c) is complete only when both the two-parameter $U_q(\widehat{sl_2})$ realization and the higher-rank $U_q(\widehat{sl_n})$ extension are settled.

References

V. V. Mangazeev and X. Lu, 'Boundary matrices for the higher spin six vertex model', Nuclear Physics B 945 (2019) 114665. DOI: 10.1016/j.nuclphysb.2019.114665. arXiv:1903.00274. https://arxiv.org/abs/1903.00274
D. Kolyaskin and V. V. Mangazeev, 'Triangular solutions to the reflection equation for $U_q(\widehat{sl_n})$ ', Journal of Physics A: Mathematical and Theoretical 57 (2024) 245201. DOI: 10.1088/1751-8121/ad4d2f. arXiv:2402.05442. https://arxiv.org/abs/2402.05442
Z. Yang, 'Stationary measures for higher spin vertex models on a strip', arXiv:2309.04897 (2023). https://arxiv.org/abs/2309.04897
J. He, 'Shift invariance of half space integrable models', Probability Theory and Related Fields 194 (2025) 541-611. DOI: 10.1007/s00440-025-01358-9. arXiv:2205.13029. https://arxiv.org/abs/2205.13029
R. Frassek, 'Integrable boundaries for the q-Hahn process', Journal of Physics A: Mathematical and Theoretical 55 (2022) 404008. DOI: 10.1088/1751-8121/ac901b. arXiv:2205.10512. https://arxiv.org/abs/2205.10512
A. Kuniba, M. Okado, and A. Yoneyama, 'Reflection $K$ matrices associated with an Onsager coideal of $U_p(A^{(1)}_{n-1})$ , $U_p(B^{(1)}_n)$ , $U_p(D^{(1)}_n)$ and $U_p(D^{(2)}_{n+1})$ ', Journal of Physics A: Mathematical and Theoretical 52 (2019) 375202. DOI: 10.1088/1751-8121/ab3715. arXiv:1904.05653. https://arxiv.org/abs/1904.05653
A. Kuniba, M. Okado, and A. Yoneyama, 'Matrix product solution to the reflection equation associated with a coideal subalgebra of $U_q(A^{(1)}_{n-1})$ ', Letters in Mathematical Physics 109 (2019) 2049-2067. DOI: 10.1007/s11005-019-01175-x. arXiv:1812.03767. https://arxiv.org/abs/1812.03767
A. Kuniba and V. Pasquier, 'Matrix product solutions to the reflection equation from three dimensional integrability', Journal of Physics A: Mathematical and Theoretical 51 (2018) 255204. DOI: 10.1088/1751-8121/aac3b4. arXiv:1802.09164. https://arxiv.org/abs/1802.09164
R. Frassek, C. Giardina, and J. Kurchan, 'Non-compact quantum spin chains as integrable stochastic particle processes', Journal of Statistical Physics 180 (2020) 135-171. DOI: 10.1007/s10955-019-02375-4. arXiv:1904.01048. https://arxiv.org/abs/1904.01048
A. Kuniba, V. V. Mangazeev, S. Maruyama, and M. Okado, 'Stochastic $R$ matrix for $U_q(A^{(1)}_n)$ ', Nuclear Physics B 913 (2016) 248-277. DOI: 10.1016/j.nuclphysb.2016.09.016. arXiv:1604.08304. https://arxiv.org/abs/1604.08304
E. K. Sklyanin, 'Boundary conditions for integrable quantum systems', Journal of Physics A: Mathematical and General 21 (1988) 2375-2389. DOI: 10.1088/0305-4470/21/10/015. https://doi.org/10.1088/0305-4470/21/10/015
H. J. de Vega and A. Gonzalez Ruiz, 'Boundary K-matrices for the six vertex and the n(2n-1)A_{n-1} vertex models', Journal of Physics A: Mathematical and General 26 (1993) L519-L524. DOI: 10.1088/0305-4470/26/12/007. https://doi.org/10.1088/0305-4470/26/12/007
V. V. Bazhanov and S. M. Sergeev, 'An Ising-type formulation of the six-vertex model', Nuclear Physics B 986 (2023) 116055. DOI: 10.1016/j.nuclphysb.2022.116055. https://doi.org/10.1016/j.nuclphysb.2022.116055
Z. Tsuboi, 'Generic triangular solutions of the reflection equation: $U_{q}(\widehat{sl_2})$ case', Journal of Physics A: Mathematical and Theoretical 53 (2020) 225202. DOI: 10.1088/1751-8121/ab8853. arXiv:1912.12808. https://arxiv.org/abs/1912.12808
Z. Tsuboi, 'Universal Baxter TQ-relations for open boundary quantum integrable systems', Nuclear Physics B 963 (2021) 115286. DOI: 10.1016/j.nuclphysb.2020.115286. arXiv:2010.09675. https://arxiv.org/abs/2010.09675
D. Chernyak, A. M. Gainutdinov, J. L. Jacobsen, and H. Saleur, 'Algebraic Bethe Ansatz for the open XXZ spin chain with non-diagonal boundary terms via U_q(sl_2) symmetry', SIGMA 19 (2023) 046. DOI: 10.3842/SIGMA.2023.046. arXiv:2212.09696. https://arxiv.org/abs/2212.09696
G. Barraquand and M. Rychnovsky, 'Random walk on nonnegative integers in beta distributed random environment', Communications in Mathematical Physics 398 (2022) 823-875. DOI: 10.1007/s00220-022-04536-1. arXiv:2201.07270. https://arxiv.org/abs/2201.07270