Shim, Sun, Wang and Yagi (arXiv:2602.22060, JHEP 06 (2026) 141) constructed three-dimensional bipartite lattice models whose R-matrices solve the bicolored tetrahedron equations (BTEs) via line defects in state integral models on shaped triangulations, with an explicit solution produced by Teichmüller TQFT. Integrability of these models, however, is unestablished: the R-matrix parameters originate in shape gauge transformations and cancel from the layer transfer matrix under simple periodic boundary conditions, and the invertibility, nondegeneracy and diagonalizability conditions needed to deduce commuting transfer matrices from the BTEs are unverified. The identification of the quantum-group-like algebraic structure conjectured to underlie BTE solutions, and the suggested three-dimensional quantum gravity interpretation of the models, are likewise open. This record poses the completion of that program: introducing genuine spectral parameters (for example by symmetry-twisted boundary conditions), establishing or refuting the integrability criteria, characterizing the underlying algebra, and grounding the gravity interpretation. A 2026 follow-up by one of the authors solves the spectral-parameter problem for a modified edge-formulation construction, but the questions for the original BTE models remain open.

Background

Two-dimensional integrable lattice models are governed by the Yang–Baxter equation, a consistency condition on the local Boltzmann weights (encoded in an R-matrix) under which the row-to-row transfer matrices commute for different values of a spectral parameter — the property from which exact solvability follows. Zamolodchikov's tetrahedron equation (1981) is the three-dimensional analogue: for interaction-round-a-cube (IRC) models, i.e. classical spin models on a cubic lattice whose eight corner spins interact around each cube, it expresses the equality of two decompositions of a rhombic dodecahedron into four cubes, and its solutions are correspondingly rare. In two dimensions the algebraic origin of most Yang–Baxter solutions is understood through quantum groups (quantized universal enveloping algebras), whose representations produce R-matrices; no comparable structure is known for the tetrahedron equation. An early idea of Maillet related the tetrahedron equation to pentagon identities for an operator represented by a tetrahedron.

Shim, Sun, Wang and Yagi proposed a geometric route to new solutions. State integral models — of which Teichmüller TQFT, a topological quantum field theory constructed by Andersen and Kashaev from quantum Teichmüller theory, is the central example — are statistical-mechanical models defined on triangulated 3-manifolds whose tetrahedra carry the shape of an ideal hyperbolic tetrahedron, i.e. dihedral angles assigned to opposite edge pairs with sum π. Their Boltzmann weights are built from the quantum dilogarithm (a q-deformation of the classical dilogarithm, characterized by its inversion and pentagon identities) and satisfy the shaped pentagon identity, which makes the partition function invariant under shaped 2–3 moves (recompositions of two tetrahedra into three that preserve total angles around edges). By introducing line defects — lines around which the total dihedral angle differs from 2π — into state integral models on shaped triangulations of the rhombic dodecahedron, the authors obtained Boltzmann weights satisfying a pair of equations they call the bicolored tetrahedron equations (BTEs), a variant of the tetrahedron equation adapted to a bipartite (two-color) cubic lattice, earlier literature's 'modified tetrahedron equations'. Teichmüller TQFT provides an explicit solution, and these weights define three-dimensional bipartite lattice models.

The BTEs imply integrability only under additional hypotheses, spelled out in the source paper. Slicing the periodic lattice into layers, the layer transfer matrix (a partial trace of a product of R-matrices, acting on tensor products of infinite-dimensional spaces of continuous edge variables) satisfies: if certain reduced R-matrices obtained by trace reduction are invertible, the layer transfer matrices of two such models commute (Proposition 2.2); and if one combined transfer matrix is diagonalizable with nondegenerate eigenvalues, then all transfer matrices commute (Lemma 2.3), yielding integrability in the standard sense of commuting conserved quantities. Two obstacles block the application of these criteria. First, the parameters of the constructed R-matrices are introduced by shape gauge transformations — shifts of the dihedral angles that leave the total angle around every internal edge invariant — so with simple periodic boundary conditions these parameters cancel from the layer transfer matrix, which therefore has no spectral parameters. Twisting the boundary conditions by symmetries of the R-matrices can introduce parameter dependence, as shown for the quantized six-vertex model on a torus by Inoue, Kuniba, Terashima and Yagi, where layer transfer matrices depending on two spectral parameters associated with the homology cycles of the torus were proven to commute using tetrahedron equations and inversion relations; but that argument is tied to the specific model. Second, the required invertibility, nondegeneracy and diagonalizability properties have not been verified for the Teichmüller-TQFT-derived R-matrices. Beyond integrability, the authors conjecture that the BTE solutions are governed by algebraic structures generalizing the quantum groups of the Yang–Baxter equation, and that the resulting lattice models inherit a three-dimensional quantum gravity interpretation from Teichmüller TQFT; both directions are unexplored.

Subsequent work has clarified the landscape without answering these questions. Bazhanov, Kashaev, Mangazeev and Sergeev constructed, by a different route, a new class of integrable 3D lattice models with continuous families of commuting layer-to-layer transfer matrices, from quantum dilogarithms satisfying inversion and pentagon identities. Yagi (2026) showed that the solutions of the source paper satisfy only the weaker (bicolored) form of the tetrahedron equation with gauge-removable parameter dependence, and gave a modified construction — mapping the edge state variables of a state integral model to the cube vertices of an IRC model — in which the tetrahedral weight solves the strict tetrahedron equation with genuine spectral parameters (invariant under shape gauge transformations up to overall shifts), producing integrable IRC models from state integral models, including the edge formulation of Teichmüller TQFT. The integrability of the original BTE models themselves, the algebra behind their solutions, and their gravity interpretation remain open.

Problem Statement

For the three-dimensional bipartite lattice models whose R-matrices are the bicolored tetrahedron equation (BTE) solutions constructed by Shim–Sun–Wang–Yagi from line defects in state integral models on shaped triangulations — concretely, the solution produced by Teichmüller TQFT — complete the integrability and characterization program that the construction leaves open:

(1) Spectral parameters. Give a concrete prescription within the state-integral/geometric framework — for example, boundary conditions twisted by suitable symmetries of the R-matrices — under which the layer transfer matrix of these models acquires genuine spectral-parameter dependence that cannot be removed by shape gauge transformations; the prescription should be model-independent, or at least explicitly implemented for the Teichmüller-TQFT solution.

(2) Integrability criteria. Verify or refute, for these R-matrices and transfer matrices, the algebraic conditions under which the BTEs imply integrability: invertibility of the relevant trace-reduced R-matrices, and the existence of a combined layer transfer matrix that is diagonalizable with nondegenerate eigenvalues, from which commutativity of the transfer matrices follows.

(3) Algebraic structure. Identify and characterize the algebraic object playing, for the (bicolored) tetrahedron equation, the role that quantum groups play for the Yang–Baxter equation: define it precisely and characterize its relation to the state-integral BTE solutions.

(4) Gravity interpretation. Make precise the proposed interpretation of the resulting lattice models via the relation of Teichmüller TQFT to three-dimensional quantum gravity.

These components are linked by a single objective: turning the state-integral solutions of the bicolored tetrahedron equations into genuinely integrable — or provably non-integrable — three-dimensional lattice models whose algebraic origin and physical meaning are understood. A complete resolution establishes or rules out integrability through (1) and (2) — where ruling it out requires an independent obstruction (a proof that no genuine spectral parameters can be introduced under any admissible boundary conditions, or that no commuting family of transfer matrices can exist), since the criteria in (1) and (2) are sufficient conditions for integrability and their failure alone does not rule it out — and characterizes the construction through (3) and (4); partial answers are judged on the components they address, with the common objective in view.

References

  1. Myungbo Shim, Xiaoyue Sun, Hao Ellery Wang, Junya Yagi, 'Solving the tetrahedron equation by Teichmüller TQFT', JHEP 06 (2026) 141, doi:10.1007/JHEP06(2026)141, arXiv:2602.22060, https://arxiv.org/abs/2602.22060
  2. Junya Yagi, 'State integral models and the tetrahedron equation', J. Phys. A: Math. Theor. 59 (2026) 285203, doi:10.1088/1751-8121/ae848f, arXiv:2603.09230, https://arxiv.org/abs/2603.09230
  3. Vladimir V. Bazhanov, Rinat M. Kashaev, Vladimir V. Mangazeev, Sergey M. Sergeev, 'Quantum Dilogarithms and New Integrable Lattice Models in Three Dimensions', arXiv:2512.23338, https://arxiv.org/abs/2512.23338
  4. Rei Inoue, Atsuo Kuniba, Yuji Terashima, Junya Yagi, 'Quantized six-vertex model on a torus', arXiv:2505.08924, https://arxiv.org/abs/2505.08924
  5. Shinsuke Iwao, Kohei Motegi, Ryo Ohkawa, 'Tetrahedral L-operators, tensor Schur polynomials and q-deformed loop elementary symmetric functions', arXiv:2604.22141, https://arxiv.org/abs/2604.22141
  6. A. B. Zamolodchikov, 'Tetrahedron equations and the relativistic S-matrix of straight-strings in 2+1-dimensions', Communications in Mathematical Physics 79 (1981) 489-505, doi:10.1007/BF01209309, https://doi.org/10.1007/BF01209309
  7. Jørgen Ellegaard Andersen, Rinat Kashaev, 'A TQFT from quantum Teichmüller theory', Communications in Mathematical Physics 330 (2014) 887-934, doi:10.1007/s00220-014-2073-2, arXiv:1109.6295, https://arxiv.org/abs/1109.6295
  8. J. M. Maillet, 'On Pentagon And Tetrahedron Equations', St. Petersburg Math. J. 6 (1995) 375-383, arXiv:hep-th/9312037, https://arxiv.org/abs/hep-th/9312037