Khachatryan's combined (restricted D-simplex) local integrability equations extend the Algebraic Bethe Ansatz to three-dimensional classical vertex models, but every solution exhibited in the founding paper fails at least one physical desideratum: spectral-parameter-dependent solutions factorize into 2D Yang-Baxter solutions or describe interactions along a single lattice direction, constant solutions carry no spectral parameter, and the genuinely three-dimensional rational solutions produce non-local Hamiltonians. The source explicitly leaves open the construction of genuinely three-dimensional, irreducible, spectral-parameter-dependent vertex R-matrices whose transfer matrices commute at all lattice sizes and whose logarithmic derivative yields a local, Hermitian two-dimensional quantum lattice Hamiltonian, and identifies as an interesting forward direction a general criterion under which commutativity verified at finite lattice sizes extends to all sizes via explicit invertible intertwiners. The problem is open: no genuinely three-dimensional, irreducible, spectral-parameter-dependent vertex R-matrix with a local Hermitian Hamiltonian is known.
Background
In a two-dimensional (2D) classical vertex model, the Boltzmann weights are collected into an R-matrix acting on the tensor product of two vector spaces and depending on a spectral parameter . The Yang-Baxter equation (YBE), , is the local compatibility condition guaranteeing that the transfer matrices , defined as traces of monodromy matrices (ordered products of -matrices along a lattice row), commute for every lattice size. This commutativity drives the Algebraic Bethe Ansatz (ABA), the standard constructive method for exact spectra of one-dimensional quantum chains: at a regular point where is a permutation operator, the quantum Hamiltonian is recovered as the logarithmic derivative , and regularity makes a sum of local nearest-neighbor densities.
Generalizing this structure to three classical dimensions (equivalently, to -dimensional quantum evolution) has been a persistent difficulty. Zamolodchikov's tetrahedron equation is the three-dimensional analogue of the YBE, introduced for the factorized scattering of string-like objects in dimensions; its vertex-type versions include the interaction-round-a-cube (IRC) equations, and the homogeneous generalizations to dimensions are the D-simplex equations. Essentially all known vertex-type solutions are constant in the spectral-parameter sense; for instance, the q-oscillator (q-deformed harmonic-oscillator algebra) solutions of Bazhanov and Sergeev define integrable 3D lattice models that, viewed in 2D, only enlarge the internal state spaces by the third direction. Simplified systems, such as the cube equations of Khachatryan, Ferraz, Klumper and Sedrakyan (whose one essentially new physical application is Kitaev's 2D quantum model, the exactly solvable honeycomb-lattice spin model), sharpen the same bottleneck: as the source notes, it is not an easy task to obtain solutions that produce local and Hermitian 2D quantum Hamiltonian operators, the obstacle being a huge number of coupled equations for undefined functions of several spectral variables.
Khachatryan's combined equations (arXiv:2210.16753, published as Eur. Phys. J. Plus 138 (2023) 1058) attack this with two devices. First, inhomogeneous local integrability equations, called combined or restricted D-simplex equations, in which the quantum R-matrix is accompanied by intertwiner matrices of smaller dimension than the quantum R-matrix, generalizing the inhomogeneous YBE . Second, a stepwise method: rather than solving the local equations directly, one studies the commutation equations of transfer matrices on lattices of increasing size ; after finitely many steps the arising constraints either repeat (closing the system) or exclude nontrivial solutions, and in the closing case one can attempt to fix intertwiners through the inhomogeneous equations. Within this framework the paper proves a no-go for spectral-parameter-dependent solutions built from the 3D identity operator plus full projection operators; its three-state solutions of 20-vertex type (so named for the twenty nonzero vertex configurations of the underlying three-state model) are constant or factorize into 2D YBE solutions; its spectral-parameter-dependent solutions describe interactions along a single lattice direction, hence effectively one-dimensional quantum models; and its genuinely three-dimensional (non-factorizable) rational solutions are explicitly non-local. The 3D chiral Potts R-matrix (a three-dimensional lift of the chiral, that is direction-asymmetric, Potts spin-model vertex weights) constructed earlier by the same authors is flagged in the source's summary as the natural next test case, with only preliminary indications that its integrability conditions reduce to 2D projections.
Problem Statement
The problem, as left open by Khachatryan's combined-equations framework, has as its core objective a certified route from genuinely three-dimensional classical vertex models to physically meaningful 2D quantum integrable models:
(1) Construct a genuinely three-dimensional irreducible vertex R-matrix depending nontrivially on a spectral parameter that satisfies the combined (restricted D-simplex / inhomogeneous YBE-type) local integrability equations with invertible intertwiners , such that (i) the corresponding transfer matrices commute, , for all lattice sizes , and (ii) the logarithmic derivative at a suitable point is a Hermitian Hamiltonian of a 2D quantum lattice model with local, finite-range interactions. The objective is a construction rather than a classification: the three-state 20-vertex-type and four-state classes considered in the source are natural first settings, but any class meeting the stated desiderata answers the open question. Genuinely three-dimensional means that does not factorize into 2D YBE solutions and that the induced interactions do not collapse onto a single lattice direction; irreducible means that cannot be brought, by a change of basis on its spaces, into a direct sum of smaller vertex blocks; all solutions exhibited in the source fail at least one of the stated desiderata. Alternatively, establish that no such R-matrix exists within the combined-equations framework.
(2) A closely related direction explicitly discussed in the source is a general criterion giving conditions under which commutativity of transfer matrices verified at finite lattice sizes, the positive output of the stepwise method, extends to commutativity for all sizes; equivalently, a systematic procedure producing invertible intertwiners satisfying the appropriate inhomogeneous Yang-Baxter or D-simplex equations uniformly in the lattice size. The source presents the stepwise scheme as a case-by-case test method and expresses interest in realizing this concept where the direct consideration of the local integrability conditions is not effective; part (2) is a related strengthening, not a required component of an answer to part (1).
A complete answer settles the load-bearing question whether the combined-equations route to genuinely 3D spectral-parameter integrability with local Hermitian quantum Hamiltonians succeeds or fails; a general finite-size-to-all-size commutativity criterion, if achieved, would additionally supply the certification tool by which individual constructions can be validated systematically.
References
- S. A. Khachatryan, 'Integrability in dimensions: combined local equations and commutativity of the transfer matrices', Eur. Phys. J. Plus 138 (2023) 1058. DOI: 10.1140/epjp/s13360-023-04711-w; arXiv:2210.16753. https://arxiv.org/abs/2210.16753
- P. Padmanabhan and V. Korepin, 'Majorana fermions solve the tetrahedron equations as well as higher simplex equations', Nucl. Phys. B 1014 (2025) 116865. DOI: 10.1016/j.nuclphysb.2025.116865; arXiv:2410.20328. https://arxiv.org/abs/2410.20328
- Sh. Khachatryan and A. Sedrakyan, 'Explicit R-matrices for inhomogeneous 3D chiral Potts models: Integrability and the action formulation for IM', Ann. Phys. 437 (2022) 168735. DOI: 10.1016/j.aop.2021.168735; arXiv:2112.09892. https://arxiv.org/abs/2112.09892
- A. Kuniba, S. Matsuike and A. Yoneyama, 'New solutions to the tetrahedron equation associated with quantized six-vertex models', Commun. Math. Phys. 401 (2023) 3247-3276. DOI: 10.1007/s00220-023-04711-y; arXiv:2208.10258. https://arxiv.org/abs/2208.10258
- S. A. Khachatryan, 'New series of multi-parametric solutions to GYBE: quantum gates and integrability', Nucl. Phys. B 996 (2023) 116375. DOI: 10.1016/j.nuclphysb.2023.116375; arXiv:2308.03564. https://arxiv.org/abs/2308.03564
- 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
- V. V. Bazhanov and S. M. Sergeev, 'Zamolodchikov's tetrahedron equation and hidden structure of quantum groups', J. Phys. A 39 (2006) 3295-3310. DOI: 10.1088/0305-4470/39/13/009; arXiv:hep-th/0509181. https://arxiv.org/abs/hep-th/0509181
- J. Yagi, 'State integral models and the tetrahedron equation', arXiv:2603.09230 (2026). https://arxiv.org/abs/2603.09230
- M. Shim, X. Sun, H. E. Wang and J. Yagi, 'Solving the tetrahedron equation by Teichmuller TQFT', arXiv:2602.22060 (2026). https://arxiv.org/abs/2602.22060
- Sh. Khachatryan, A. Ferraz, A. Klumper and A. Sedrakyan, 'Integrability of three dimensional models: cubic equations', arXiv:1502.04055 (2015). https://arxiv.org/abs/1502.04055
- A. Sinha, P. Padmanabhan and V. Korepin, 'Toffoli gates solve the tetrahedron equations', arXiv:2405.16477 (2024). https://arxiv.org/abs/2405.16477