Integrable construction of a two-dimensional lattice model with anisotropic Hubbard couplings

By defining a graded global R-operator $\mathbb{R}_{ab}^{(2D,2S)}$ that couples free-fermion structures and incorporates anisotropic Hubbard interactions while satisfying the Yang--Baxter equation, we construct a strictly solvable two-dimensional lattice model. We then build the layer-to-layer transfer matrix through a bidirectional-monodromy construction and prove the model's integrability via the associated global RTT relations. Using the nested algebraic Bethe ansatz, we obtain the exact eigenvalues of the transfer matrix and derive the corresponding first- and second-level Bethe equations. Finally, by taking the logarithmic derivative of the transfer matrix at the regular point, we recover explicitly a local Hamiltonian that features anisotropic hopping, an on-site Hubbard interaction, and orbital-coupling contributions.