Quantum simulations of many-body systems using 2-qubit Yang-Baxter gates offer a benchmark for quantum hardware. This can be extended to the higher dimensional case with
n-qubit generalisations of Yang-Baxter gates called
n-simplex operators. Such multi-qubit gates potentially lead to shallower and more efficient quantum circuits as well. Finding them amounts to identifying unitary solutions of the
n-simplex equations, the building blocks of higher dimensional integrable systems. These are a set of highly non-linear and over determined system of equations making it notoriously hard to solve even when the local Hilbert spaces are spanned by qubits. We systematically overcome this for higher simplex operators constructed using two methods: from Clifford algebras and by lifting Yang-Baxter operators. The
n=3 or the tetrahedron case is analyzed in detail. For the qubit case our methods produce 13 inequivalent families of unitary tetrahedron operators. 12 of these families are obtained by appending the 5 unitary families of 4 by 4 constant Yang-Baxter operators of Dye-Hietarinta, with a single qubit operator. As applications, universal sets of single, two and three qubit gates are realized using such unitary tetrahedron operators. The ideas presented in this work can be naturally extended to the higher simplex cases.
The circuit model of quantum computation can be interpreted as a scattering process. In particular, factorised scattering operators result in integrable quantum circuits that provide universal quantum computation and are potentially less noisy. These are realized through Yang-Baxter or 2-simplex operators. A natural question is to extend this construction to higher qubit gates, like the Toffoli gates, which also lead to universal quantum computation but with shallower circuits. We show that unitary families of such operators are constructed by the 3-dimensional generalizations of the Yang-Baxter operators known as tetrahedron or 3-simplex operators. The latter satisfy a spectral parameter-dependent tetrahedron equation. This construction goes through for
n-Toffoli gates realized using
n-simplex operators.
A new realization of doubling degeneracy based on emergent Majorana operator Γ presented by Lee-Wilczek has been made. The Hamiltonian can be obtained through the new type of solution of Yang-Baxter equation, i.e. -matrix. For 2-body interaction, R(θ) gives the "superconducting" chain that is the same as 1D Kitaev chain model. The 3-body Hamiltonian commuting with Γ is derived by 3-body R₁₂₃-matrix, we thus show that the essence of the doubling degeneracy is due to [R(θ), Γ=0]. We also show that the extended Γ'-operator is an invariant of braid group BN for odd N. Moreover, with the extended Γ'-operator, we construct the high dimensional matrix representation of solution to Yang-Baxter equation and find its application in constructing 2N-qubit Greenberger-Horne-Zeilinger state for odd N.
Integrating infinitesimal generators of local dynamics into finite-time evolution is generally harder than extracting a local Hamiltonian from a quantum gate. For integrable one-dimensional quantum many-body systems, a recent bootstrap constructs one-parameter
R-operators directly from Hamiltonians satisfying the Reshetikhin condition. However, this condition excludes important models such as the Hubbard model, whose two-parameter
R-operator was obtained by Shastry from a decorated Yang--Baxter equation obeyed by its free-fermion sector. Here, I extend the bootstrap to
R-operators constrained by both the ordinary and decorated Yang--Baxter equations. Their joint expansion produces additional Hamiltonian-level necessary conditions and recursively determines candidate higher-order terms, subject to two-site locality at each order. I also derive compatibility and factorization conditions under which a Shastry-type superposition solves the
RLL relation after deformation by an onsite involution. The construction recovers the known two-parameter
R-operators of the Hubbard and XYh models, whose full Yang--Baxter equations follow from established tetrahedral Zamolodchikov algebra relations. For a superconducting Hubbard deformation, the scalar equations are generically incompatible, placing it outside this framework. These results provide diagnostics and a systematic route to fermionic
R-operators beyond the Reshetikhin condition. General sufficiency, convergence of the formal bootstrap, and sufficiency of the deformation conditions remain open.