Baxter Q-operator for the open $\tau_2$-model with arbitrary boundary fields

The $\tau_2$-model (Baxter–Bazhanov–Stroganov model) is the integrable lattice model built on $p$-dimensional cyclic representations of the Weyl algebra at roots of unity; it underlies the chiral Potts and lattice sine-Gordon models. For the open chain with generic site-dependent inhomogeneity and arbitrary non-diagonal integrable boundary fields, the transfer-matrix eigenvalues are known through an inhomogeneous $T$-$Q$ relation (off-diagonal Bethe ansatz), and the full spectrum — eigenvalues and eigenstates — has since been characterized by separation of variables, including an explicit inhomogeneous Baxter-type functional equation. These are eigenvalue-level or SoV-basis statements, however: the Baxter $Q$-operator itself — an explicit operator family commuting with the open-chain transfer matrix and satisfying an operator-level $T$-$Q$ relation from which the known functional relations follow — has not been constructed for this model. Recent work has built such $Q$-operators at roots of unity for the periodic 6-vertex and periodic $\tau_2$-models and explicitly identifies the open-chain case with arbitrary boundary conditions as an open problem. The problem posed here is to construct the Baxter $Q$-operator for the open $\tau_2$-model with generic site-dependent inhomogeneity and arbitrary non-diagonal boundary fields.

Background

Integrable lattice models and quantum spin chains are diagonalized through their transfer matrix $t(u)$: a one-parameter commuting family built from local Lax operators $L_n(u)$ that obey the Yang–Baxter algebra (the quadratic exchange relations imposed by the six-vertex $R$-matrix — the fundamental trigonometric solution of the Yang–Baxter equation — on the monodromy matrix, the ordered product of Lax operators along the chain). The $\tau_2$-model (also called the Baxter–Bazhanov–Stroganov model) arises when each local quantum space carries a $p$-dimensional cyclic representation of the Weyl algebra — generators $X_n, Z_n$ with $X_n^p = Z_n^p = 1$ and $X_n Z_n = q\,Z_n X_n$ at the root of unity $q = e^{\eta}$, $\eta = 2\mathrm{i}\pi/p$ with $p$ odd. The associated Bazhanov–Stroganov Lax operators, whose entries depend on site-dependent inhomogeneity parameters constrained only by the cyclic-representation condition, solve the Yang–Baxter algebra, and the resulting transfer matrix defines the $\tau_2$-model. Under further parameter constraints this model reduces to the chiral Potts model and to the lattice sine-Gordon model, and its spectral theory is the algebraic engine behind the exact solution of the chiral Potts model.

Open integrable boundaries are imposed through reflection equations (the boundary analogues of the Yang–Baxter equation) for two boundary matrices $K^\pm(u)$ (the Sklyanin framework), and the double-row transfer matrix is $t(u) = \operatorname{tr}\{K^+(u)\,T(u)\,K^-(u)\,\hat{T}(u)\}$, where $T(u)$ is the monodromy matrix built from the Lax operators and $\hat{T}(u)$ is its reflected product. The most general integrable boundary fields are non-diagonal: each $K$-matrix carries three boundary parameters, so the usual diagonal reference state of the algebraic Bethe ansatz is unavailable.

For this open $\tau_2$-model with generic site-dependent inhomogeneity and arbitrary boundary fields, Xu et al. obtained the eigenvalues by the off-diagonal Bethe ansatz (ODBA), a method in which eigenvalues are reconstructed from operator product identities, asymptotics, and analyticity instead of a reference state. The result is an inhomogeneous $T$-$Q$ relation: the eigenvalue $\Lambda(u)$ equals $a(u)\,Q(u-\eta)/Q(u) + d(u)\,Q(u+\eta)/Q(u)$ plus a third term proportional to a Laurent polynomial $F(u)$ built from averages of monodromy-matrix entries, with $Q(u)$ a trigonometric polynomial parameterized by $(p-1)N+2p$ Bethe roots fixed by Bethe ansatz equations. This is an eigenvalue-level functional ansatz, not an operator identity; the inhomogeneous term vanishes (and the relation collapses to the conventional Baxter form) only on parameter loci for which the paper gives sufficient constraints, and the average-value functions entering $F(u)$ are computed explicitly only for one- and two-site chains. Subsequently, Maillet, Niccoli and Pezelier applied Sklyanin's separation of variables (SoV) — the representation-theoretic method that diagonalizes the transfer matrix in the eigenbasis of a chosen operator, reducing the spectral problem to polynomial equations in one variable — to the same model and completely characterized the spectrum, eigenvalues and eigenstates included, for the most general boundary K-matrices, proving that an inhomogeneous Baxter-type functional equation with an explicitly computed inhomogeneous term captures the full spectrum.

What is still absent is Baxter's $Q$-operator: an explicit family of operators $\mathcal{Q}(u)$ commuting with $t(u)$ and satisfying an operator-level $T$-$Q$ relation from which these functional relations follow, with zeros of $Q$-eigenvalues reproducing the Bethe roots. Historically, the $\tau_2$ $Q$-operator is tied to the half-monodromy matrix of the chiral Potts model (an observation of Bazhanov and Stroganov) and lives over the high-genus chiral Potts spectral curve (the algebraic curve carrying the chiral Potts spectral parameters), which has made explicit constructions difficult. The periodic case is now understood — the Q-operator was reconstructed from SoV for cyclic representations and recently constructed algebraically at roots of unity from cyclic representations of the Borel subalgebra (the upper-triangular half) of the quantum affine algebra $U_q(\widehat{\mathfrak{sl}}_2)$ (the $q$-deformed enveloping algebra of the affine $\mathfrak{sl}_2$ loop algebra), for both the 6-vertex and the $\tau_2$ transfer matrices — while for open chains Q-operators are known only at generic $q$ and for restricted boundary classes, with universal $T$-$Q$ relations available only for diagonal boundaries. For the open $\tau_2$-model at roots of unity with the most general non-diagonal boundary fields — precisely the setting of the 2016 ODBA solution — no Q-operator exists in the literature, and current work (Weston, SciPost Phys. Core 9, 016 (2026)) identifies its construction as an open goal, the barrier being the solution of the reflection equations with the relevant auxiliary representations.

Problem Statement

Construct an explicit Baxter $Q$-operator for the open $\tau_2$-model at root-of-unity anisotropy $\eta = 2\mathrm{i}\pi/p$ ($p$ odd) with generic site-dependent inhomogeneity and arbitrary (non-diagonal) integrable boundary fields — the most general boundary K-matrices solving the reflection equations.

A complete answer must deliver a family of operators $\mathcal{Q}(u)$, given by a concrete formula (for example as a trace or double-row monodromy over a finite-dimensional auxiliary representation, such as the cyclic representations of the Borel subalgebra of $U_q(\widehat{\mathfrak{sl}}_2)$ at roots of unity, combined with boundary K-operators solving the corresponding reflection equations), such that: (1) $[\mathcal{Q}(u),\,t(u)] = 0$ for the double-row transfer matrix $t(u)$ of the open $\tau_2$-model; (2) the family satisfies a proven operator-level $T$-$Q$ relation from which the established eigenvalue-level relations follow as eigenvalue consequences — either the inhomogeneous $T$-$Q$ relation and Bethe ansatz equations of the off-diagonal Bethe ansatz, or the explicit inhomogeneous Baxter equation of the separation-of-variables characterization, including its inhomogeneous term; (3) the analytic structure is explicit enough that the zeros of the $Q$-eigenvalues reproduce, for generic parameters, the Bethe roots of whichever established functional relation the construction realizes — the $(p-1)N+2p$ Bethe roots of the off-diagonal Bethe ansatz formulation, or the zeros of the degree-$(p-1)N$ $Q$-polynomial of the separation-of-variables characterization — with the degenerate parameter loci where the inhomogeneous term vanishes either covered by the construction or explicitly delimited.

References

1. X. Xu, K. Hao, T. Yang, J. Cao, W.-L. Yang, K. Shi, 'Bethe ansatz solutions of the $\tau_2$-model with arbitrary boundary fields', JHEP 11 (2016) 080, DOI: 10.1007/JHEP11(2016)080, arXiv:1608.06209, https://arxiv.org/abs/1608.06209
2. V. V. Bazhanov, Yu. G. Stroganov, 'Chiral Potts model as a descendant of the six vertex model', J. Stat. Phys. 59 (1990) 799, DOI: 10.1007/BF01025851, https://doi.org/10.1007/BF01025851
3. E. K. Sklyanin, 'Boundary Conditions for Integrable Quantum Systems', J. Phys. A: Math. Gen. 21 (1988) 2375, DOI: 10.1088/0305-4470/21/10/015, https://doi.org/10.1088/0305-4470/21/10/015
4. G. Niccoli, 'Reconstruction of Baxter Q-operator from Sklyanin SOV for cyclic representations of integrable quantum models', Nucl. Phys. B 835 (2010) 263, DOI: 10.1016/j.nuclphysb.2010.03.009, arXiv:1001.0035, https://arxiv.org/abs/1001.0035
5. N. Grosjean, G. Niccoli, 'The $\tau_2$-model and the chiral Potts model revisited: completeness of Bethe equations from Sklyanin's SOV method', J. Stat. Mech. (2012) P11005, DOI: 10.1088/1742-5468/2012/11/P11005, arXiv:1205.4614, https://arxiv.org/abs/1205.4614
6. J.-M. Maillet, G. Niccoli, B. Pezelier, 'Transfer matrix spectrum for cyclic representations of the 6-vertex reflection algebra I', SciPost Phys. 2, 009 (2017), DOI: 10.21468/SciPostPhys.2.1.009, arXiv:1607.02983, https://arxiv.org/abs/1607.02983
7. J.-M. Maillet, G. Niccoli, B. Pezelier, 'Transfer matrix spectrum for cyclic representations of the 6-vertex reflection algebra II', SciPost Phys. 5, 026 (2018), DOI: 10.21468/SciPostPhys.5.3.026, arXiv:1802.08853, https://arxiv.org/abs/1802.08853
8. Y. Miao, J. Lamers, V. Pasquier, 'On the Q operator and the spectrum of the XXZ model at root of unity', SciPost Phys. 11, 067 (2021), DOI: 10.21468/SciPostPhys.11.3.067, arXiv:2012.10224, https://arxiv.org/abs/2012.10224
9. R. Weston, 'Cyclic representations of $U_q(\widehat{\mathfrak{sl}}_2)$ and its Borel subalgebras at roots of unity and Q-operators', SciPost Phys. Core 9, 016 (2026), DOI: 10.21468/SciPostPhysCore.9.1.016, arXiv:2412.14811, https://arxiv.org/abs/2412.14811
10. B. Vlaar, R. Weston, 'A Q-operator for open spin chains I. Baxter's TQ relation', J. Phys. A 53 (2020) 245205, DOI: 10.1088/1751-8121/ab8854, arXiv:2001.10760, https://arxiv.org/abs/2001.10760
11. A. Cooper, B. Vlaar, R. Weston, 'A Q-operator for open spin chains II: boundary factorization', Commun. Math. Phys. 405 (2024) 110, DOI: 10.1007/s00220-024-04973-0, arXiv:2301.03997, https://arxiv.org/abs/2301.03997
12. Z. Tsuboi, 'Universal Baxter TQ-relations for open boundary quantum integrable systems', Nucl. Phys. B 963 (2021) 115286, DOI: 10.1016/j.nuclphysb.2020.115286, arXiv:2010.09675, https://arxiv.org/abs/2010.09675
13. J.-M. Maillet, G. Niccoli, 'On Separation of Variables for Reflection Algebras', J. Stat. Mech. (2019) 094020, DOI: 10.1088/1742-5468/ab357a, arXiv:1904.00852, https://arxiv.org/abs/1904.00852