The Mangazeev-Bazhanov-Sergeev (MBS) model is the first solvable three-dimensional (3D) lattice model with non-negative Boltzmann weights: its local weights, built from a q-oscillator solution of the tetrahedron equation, generate a two-parameter commuting family of layer-to-layer transfer matrices. The 2013 construction left the model formally integrable but not solved: the authors explicitly flagged the convergence of the full partition function Z=∑I≥0uIZI\mathcal{Z}=\sum_{\mathcal{I}\ge0}u^{\mathcal{I}}\mathcal{Z}_{\mathcal{I}} as requiring additional study, computed no spectrum of the layer-to-layer transfer matrices, deferring 'further properties of the proposed model' to elsewhere, and called a detailed analysis of the rank-size duality's implications for the Bethe ansatz extremely interesting. This problem asks for the exact solution of the 3D model: a rigorous convergence domain for the I\mathcal{I}-sum, the spectral solution of the layer-to-layer transfer matrix, and an account of how the rank-size duality reorganizes the spectral equations. The two-dimensional descendants were developed extensively over 2013-2026, while the genuinely three-dimensional questions posed here remain open.

Background

An exactly solvable lattice model in statistical mechanics is one whose partition function Z\mathcal{Z}, a sum over configurations of products of local Boltzmann weights, can be analyzed exactly; in the strongest cases this yields the free energy per site in the thermodynamic limit (the limit of an infinite lattice, in which the free energy per site becomes well defined and governs the phase diagram) and the full spectrum of the transfer matrix, the operator that builds the lattice slice by slice and whose eigenvalues encode the thermodynamics. In two dimensions the key consistency condition on the weights is the Yang-Baxter equation, which makes transfer matrices with different spectral parameters commute and underlies the Bethe ansatz, the technique that parameterizes eigenvectors by auxiliary quasimomenta obeying algebraic Bethe equations. The three-dimensional analogue is the tetrahedron equation, introduced by A. B. Zamolodchikov in 1980-81 as a factorization condition for the scattering of straight strings in (2+1)-dimensional space-time. R. J. Baxter verified Zamolodchikov's solution and computed the partition function per site of the associated 3D lattice model. Later 3D solvable models (Bazhanov-Baxter 1992 and successors) had weights of varying sign or were solved only in layered form; a long-standing gap was a solvable 3D model whose weights are all non-negative, the regime in which the partition function has a direct probabilistic interpretation. Mangazeev, Bazhanov and Sergeev (2013) filled this gap. Spins sit on the edges of an L×M×NL\times M\times N cubic lattice with periodic boundary conditions and take infinitely many values 0,1,2,…0,1,2,\ldots; each spin is a basis vector of a qq-oscillator Fock space (the occupation-number space of the qq-deformed harmonic oscillator). The local weight is a matrix element of a constant solution R\mathsf{R} of the vertex-form tetrahedron equation acting on three such Fock spaces (Eq. (31) of the source), dressed with eight continuous parameters (Eq. (38)); for 0<q<10<q<1 all weights are non-negative, the first such example. The partition function on the finite lattice is organized into restricted components ZI\mathcal{Z}_{\mathcal{I}}, I=0,1,2,…\mathcal{I}=0,1,2,\ldots (Eq. (55)), where I\mathcal{I} fixes the total spin carried by one family of edges in each horizontal layer, and the full partition function is

Z=∑I≥0uIZI,u<1,\mathcal{Z}=\sum_{\mathcal{I}\ge0}u^{\mathcal{I}}\mathcal{Z}_{\mathcal{I}},\qquad u<1,

where uu is related to the vertical edge weights and a regularizing parameter μ\mu and fugacities v,wv,w enter through Eq. (55). Each ZI\mathcal{Z}_{\mathcal{I}} converges under the sufficient conditions μ>1\mu>1, or μ=1\mu=1 with v<1,w<1v<1,w<1 (Eq. (57)), but the convergence of the I\mathcal{I}-sum itself, in the authors' words, 'requires an additional study' and may need an I\mathcal{I}-dependent damping factor. The model's integrability rests on a two-parameter commuting family of layer-to-layer transfer matrices T(v,w)\mathsf{T}(v,w) (Eq. (60)). Composing weights along a row produces a composite R-matrix (Eq. (66)) that satisfies the Yang-Baxter equation and decomposes as an infinite direct sum of R-matrices of the quantum affine algebra Uq(sl^N)U_q(\widehat{sl}_N), the Hopf-algebraic qq-deformation of the loop algebra of slNsl_N that is the standard algebraic origin of two-dimensional solvable vertex models, intertwining symmetric tensor representations πI\pi_{\mathcal{I}}, the rank-I\mathcal{I} symmetric powers of the fundamental representation (Eq. (74)). Consequently the 3D transfer matrix is an infinite sum of fused (higher-spin) two-dimensional Uq(sl^N)U_q(\widehat{sl}_N) transfer matrices (Eqs. (84)-(85)). A remarkable corollary, the rank-size duality, rewrites the same transfer matrix with the rank NN of the algebra exchanged with the length MM of the spin chain, and the spectral parameter exchanged with the horizontal fugacity; special instances of this duality were found earlier by Bazhanov-Baxter (1992) and Bazhanov-Sergeev (2006). At the time of the source paper, no genuinely 3D transfer-matrix spectrum had ever been computed: Bethe ansatze existed only for two-dimensional layered reductions of 3D models, such as the three-layer Zamolodchikov model (Boos-Mangazeev 1999). The 2013 paper leaves exactly this frontier open: it explicitly flags the convergence of the I\mathcal{I}-sum and calls for a study of the implications of the rank-size duality for the Bethe ansatz, while the spectrum of the layer-to-layer transfer matrix is left among the 'further properties of the proposed model' to be 'considered elsewhere'.

Problem Statement

Complete the exact solution of the Mangazeev-Bazhanov-Sergeev integrable 3D lattice model, taking the model exactly as defined in the source paper (J. Phys. A 46 (2013) 465206): spins 0,1,2,…0,1,2,\ldots on the edges of the cubic lattice, weights of Eqs. (31) and (38) with 0<q<10<q<1, regularizer μ\mu, fugacities v,w,uv,w,u, periodic boundary conditions - without restricting the parameter regime, lattice geometry, or model class. Three linked results are required, jointly constituting the exact solution of the model. (1) Convergence of the partition function. Prove or disprove convergence of the series Z=∑I≥0uIZI\mathcal{Z}=\sum_{\mathcal{I}\ge0}u^{\mathcal{I}}\mathcal{Z}_{\mathcal{I}} (Eq. (58) of the source) on finite lattices; determine the exact domain of the parameters (q,μ,v,w,u)(q,\mu,v,w,u) in which it converges, and whether an I\mathcal{I}-dependent damping factor, as anticipated by the authors, is necessary. Where the series diverges, characterize the large-I\mathcal{I} asymptotics of ZI\mathcal{Z}_{\mathcal{I}}. (2) Transfer-matrix spectrum. Solve the spectral problem for the two-parameter commuting family of layer-to-layer transfer matrices T(v,w)\mathsf{T}(v,w) (Eq. (60)): eigenvalues and eigenvectors in each sector labelled by the vertical spin configuration, through a Bethe-ansatz or equivalent treatment of the infinite sums of fused Uq(sl^N)U_q(\widehat{sl}_N) transfer matrices in the decomposition (74), including control of the convergence of those infinite sums. (3) Rank-size duality. Explain how the rank-size duality of Eqs. (84)-(85) constrains and reorganizes the Bethe-ansatz structure and the spectral equations - specifically, how the exchange of the algebra rank NN with the chain length MM, together with the exchange of the spectral parameter with the horizontal fugacity, acts on the Bethe equations and the eigenvalue formulae. These three components form one scientific objective: item (1) is the rigorous foundation without which item (2) has no well-defined object, item (2) is the exact solution of the model's spectral structure, and item (3) is a structural symmetry that the spectral equations of item (2) must realize. The problem is fully solved when all three items are established for the model as defined in the source; partial resolutions of individual items constitute progress but not completion.

References

  1. V. V. Mangazeev, V. V. Bazhanov and S. M. Sergeev, 'An integrable 3D lattice model with positive Boltzmann weights', J. Phys. A: Math. Theor. 46 (2013) 465206. DOI: 10.1088/1751-8113/46/46/465206, arXiv:1308.4773, https://arxiv.org/abs/1308.4773
  2. A. B. Zamolodchikov, 'Tetrahedron equations and the relativistic S-matrix of straight strings in (2+1)-dimensions', Commun. Math. Phys. 79 (1981) 489-505. DOI: 10.1007/BF01209309, https://doi.org/10.1007/bf01209309
  3. R. J. Baxter, 'On Zamolodchikov's solution of the tetrahedron equations', Commun. Math. Phys. 88 (1983) 185-205. DOI: 10.1007/BF01209476, https://doi.org/10.1007/bf01209476
  4. V. V. Bazhanov and R. J. Baxter, 'New solvable lattice models in three dimensions', J. Stat. Phys. 69 (1992) 453-485. DOI: 10.1007/BF01050423, https://doi.org/10.1007/bf01050423
  5. V. V. Bazhanov and S. M. Sergeev, 'Zamolodchikov's tetrahedron equation and hidden structure of quantum groups', J. Phys. A: Math. Gen. 39 (2006) 3295-3310. DOI: 10.1088/0305-4470/39/13/009, arXiv:hep-th/0509181, https://arxiv.org/abs/hep-th/0509181
  6. H. E. Boos and V. V. Mangazeev, 'Bethe ansatz for the three-layer Zamolodchikov model', J. Phys. A: Math. Gen. 32 (1999) 5285-5298. DOI: 10.1088/0305-4470/32/28/308, https://doi.org/10.1088/0305-4470/32/28/308
  7. V. V. Mangazeev, 'On the Yang-Baxter equation for the six-vertex model', Nucl. Phys. B 882 (2014) 70-96. DOI: 10.1016/j.nuclphysb.2014.02.019, arXiv:1401.6494, https://arxiv.org/abs/1401.6494
  8. V. V. Mangazeev, 'Q-operators in the six-vertex model', Nucl. Phys. B 886 (2014) 166-184. DOI: 10.1016/j.nuclphysb.2014.06.024, arXiv:1406.0662, https://arxiv.org/abs/1406.0662
  9. G. Bosnjak and V. V. Mangazeev, 'Construction of R-matrices for symmetric tensor representations related to Uq(sl^n)U_q(\widehat{sl}_n)', J. Phys. A: Math. Theor. 49 (2016) 495204. DOI: 10.1088/1751-8113/49/49/495204, arXiv:1607.07968, https://arxiv.org/abs/1607.07968
  10. S. Khachatryan, A. Ferraz, A. Klumper and A. Sedrakyan, 'Integrability in three dimensions: Algebraic Bethe ansatz for anyonic models', Nucl. Phys. B 899 (2015) 444-450. DOI: 10.1016/j.nuclphysb.2015.08.007, https://doi.org/10.1016/j.nuclphysb.2015.08.007
  11. A. Kuniba, S. Maruyama and M. Okado, 'Multispecies TASEP and the tetrahedron equation', J. Phys. A: Math. Theor. 49 (2016) 114001. DOI: 10.1088/1751-8113/49/11/114001, https://doi.org/10.1088/1751-8113/49/11/114001
  12. V. V. Bazhanov and S. M. Sergeev, 'An Ising-type formulation of the six-vertex model', Nucl. Phys. B 986 (2023) 116055. DOI: 10.1016/j.nuclphysb.2022.116055, arXiv:2205.10708, https://arxiv.org/abs/2205.10708
  13. S. Iwao, K. Motegi and R. Ohkawa, 'Tetrahedral L-operators, tensor Schur polynomials and q-deformed loop elementary symmetric functions', arXiv:2604.22141 (2026), https://arxiv.org/abs/2604.22141
  14. V. V. Bazhanov, R. M. Kashaev, V. V. Mangazeev and S. M. Sergeev, 'Quantum Dilogarithms and New Integrable Lattice Models in Three Dimensions', arXiv:2512.23338 (2025), https://arxiv.org/abs/2512.23338
  15. A. Kuniba, 'Quantum Groups in Three-Dimensional Integrability', Springer, Singapore (2022). DOI: 10.1007/978-981-19-3262-5, https://doi.org/10.1007/978-981-19-3262-5