The small polaron, an one-dimensional lattice model of interacting spinless fermions, with generic non-diagonal boundary terms is studied by the off-diagonal Bethe ansatz method. The presence of the Grassmann valued non-diagonal boundary fields gives
rise to a typical $U(1)$-symmetry-broken fermionic model. The exact spectra of the Hamiltonian and the associated Bethe ansatz equations are derived by constructing an inhomogeneous $T-Q$ relation.
The Izergin-Korepin model with general non-diagonal boundary terms, a typical integrable model beyond A-type and without U(1)-symmetry, is studied via the off-diagonal Bethe ansatz method. Based on some intrinsic properties of the R-matrix and the K-
matrices, certain operator product identities of the transfer matrix are obtained at some special points of the spectral parameter. These identities and the asymptotic behaviors of the transfer matrix together allow us to construct the inhomogeneous T-Q relation and the associated Bethe ansatz equations. In the diagonal boundary limit, the reduced results coincide exactly with those obtained via other methods.
