A class of recently introduced su (n) "free-fermion" models has recently been used to construct generalized Hubbard models. I derive an algebra defining the "free-fermion" models and give new classes of solutions. I then introduce a conjugation matrix and give a new and simple proof of the corresponding decorated Yang–Baxter equation. This provides the algebraic tools required to couple in an integrable way two copies of free-fermion models. Complete integrability of the resulting Hubbard-like models is shown by exhibiting their L and R matrices. Local symmetries of the models are discussed. The diagonalization of the free-fermion models is carried out using the algebraic Bethe Ansatz.
A new integrable model which is a variant of the one-dimensional Hubbard model is proposed. The integrability of the model is verified by presenting the associated quantum R-matrix which satisfies the Yang–Baxter equation. We argue that the new model possesses the SO(4) algebra symmetry, which contains a representation of the η-pairing SU(2) algebra and a spin SU(2) algebra. Additionally, the algebraic Bethe ansatz is studied by means of the quantum inverse scattering method. The spectrum of the Hamiltonian, eigenvectors, as well as the Bethe ansatz equations, are discussed.
Abstract In this paper we find new integrable one-dimensional lattice models of electrons. We describe all such nearest-neighbour integrable models with s u ( 2 ) × s u ( 2 ) symmetry by classifying solutions of the Yang–Baxter equation following the procedure first introduced in [1]. We find 12 R-matrices of difference form, some of which can be related to known models such as the XXX spin chain and the free Hubbard model, and some are new models. In addition, integrable generalizations of the Hubbard model are found by keeping the kinetic term of the Hamiltonian and adding all terms which preserve fermion number. We find that most of the new models cannot be diagonalized using the standard nested Bethe Ansatz.
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.
The integrability of the one-dimensional (1D) fermion chain model is investigated in the framework of the Quantum Inverse Scattering Method (QISM). We introduce a new R-operator for the fermion chain model, which is expressed in terms of the fermion operators. The R-operator satisfies a new type of the Yang-Baxter relation with fermionic L-operator. We derive the fermionic Sutherland equation from the relation, which is equivalent to the fermionic Lax equation. It also provides a mathematical foundation of the boost operator approach for the fermion model. In fact, we obtain some higher conserved quantities of the fermion model using the boost operator.