On the Bethe states of the one-dimensional supersymmetric t-J model with generic open boundaries

By combining the algebraic Bethe ansatz and the off-diagonal Bethe ansatz, we investigate the supersymmetric t-J model with generic open boundaries. The eigenvalues of the transfer matrix are given in terms of an inhomogeneous T-Q relation, and the corresponding eigenstates are expressed in terms of nested Bethe states which have well-defined homogeneous limit. This exact solution provides basis for further analyzing the thermodynamic properties and correlation functions of the model.


Introduction
The t − J model is one of the cornerstones in the study of high-T c superconductivity [1], which is a large-U limit of the single-band Hubbard model [2][3][4][5]. The Hamiltonian of the model have played essential roles in theoretical study of strongly correlated copperoxide based materials [6]. In general, the Hamiltonian of the supersymmetric t − J model with the general boundary interaction terms is given by P c + j,α c j+1,α + c + j+1,α c j,α P + J L−1 k=1 S k · S k+1 − 1 4 n k n k+1 + L−1 l=1 n l + n l+1 −µN + ξ 1 n 1 + 2h z 1 S z 1 + 2h − 1 S − 1 + 2h + 1 S + 1 + ξ L n L + 2h z L S z L + 2h − L S − L + 2h + L S + L , (1.1) where t is the nearest neighbor hopping of electrons and J is the antiferomagetic exchange; L is the total number of lattice sites; The operators c j,σ and c + j,σ are the annihilation and creation operators of the electron with spin σ = ±1 on the lattice site j, which satisfies anticommutation relations, i.e., {c + i,σ , c j,τ } = δ i,j δ σ,τ . There are only three possible states at the lattice site i due to the factor P = (1 − n j,−σ ) ruled out double occupancies; The operator n j = σ=± n j,σ means the total number operator on site j and n j,σ = c + j,σ c j,σ ; µ is the chemical potential andN = L j=1 n j ; ξ 1,L are the boundary chemical potentials; h z 1,L and h ± 1,L are the boundary fields; The spin operators S − = L j=1 S − j , S + = L j=1 S + j and S z = L j=1 S z j form the su(2) algebra and can be expressed by S − j = c + j,1 c j,−1 , S + j = c + j,−1 c j , S z j = 1 2 (n j,1 − n j,−1 ). (
In this paper, we study the supersymmetric t − J model with generic integrable boundary conditions in grading: bosonic, fermionic and fermionic (BFF). By combining the graded nested algebraic Bethe ansatz and off-diagonal Bethe ansatz, we obtain the Bethe states which have well-defined homogeneous limit and the corresponding eigenvalues of the transfer matrix of the model. Numerical results for the small size systems suggest that the spectrum obtained by the nested Bethe ansatz equations (BAEs) is complete.
The paper is organized as follows. In section 2, the associated graded R-matrix and corresponding generic integral non-diagonal boundary reflection matrices are introduced. In section 3, by using the graded algebraic Bethe ansatz, we derive the eigenvalues of the transfer matrix of the system which related with the eigenvalues of the nested transfer matrix. In section 4, the eigenvalues of the nested transfer matrix are derived by offdiagonal Bethe ansatz, and the Bethe states are also be given. In section 5, we construct the nested inhomogeneous T − Q relation and the nested Bethe ansatz equations of the supersymmetric t − J model. Section 6 contains our results and give some discussions.

Integrability of the model
In this paper we consider J = 2t = 2 which corresponds to the supersymmetric and integrable point [29]. The integrability of the model is associated with the rational R-matrix R(u) given by The R-matrix R(u) possesses the following properties Initial condition: R 12 (0) = ηP 12 , (2.2) Unitarity relation: Crossing Unitarity relation: Here P 12 is the graded permutation operator with the definition

5)
p(α i ) is the Grassmann parities which is one for fermions and zero for bosons. Here, we choose BFF grading which means p(1) = 0, p(2) = p(3) = 1 and R 21 (u) = P 12 R 12 (u)P 12 , st i denotes the super transposition in the i-th space (A st ) ij = A ji (−1) p(i)[p(i)+p(j)] and ist i denotes the inverse super transposition. The functions ρ 1 (u) and ρ 2 (u) are given by Here and below we adopt the standard notations: For any matrix A ∈ End(V), A j is an super embedding operator in the Z 2 graded tensor space V ⊗ V ⊗ · · · , which acts as A on the j-th space and as identity on the other factor spaces. For R ∈ End(V ⊗ V), R ij is an super embedding operator of R in the Z 2 graded tensor space, which acts as identity on the factor spaces except for the i-th and j-th ones. The super tensor product of two operators are defined through (For further details we refer the reader to [30]).
The R-matrix is an even operator (i.e., the parities of the non-zero matrix elements R ac bd of the R-matrix satisfies p(a) + p(b) + p(c) + p(d) = 0) and satisfies the graded quantum Yang-Baxter equation (QYBE) In terms of the matrix entries, it reads Let us now introduce the reflection matrix K − (u) and its dual one K + (u). The former satisfies the graded reflection equation (RE) [31] 9) and the latter satisfies the dual RE which take the form [32] For our case, the dual reflection equation (2.10) reduces to In this paper we consider the generic non-diagonal K-matrices K − (u) (2.14) Here the four boundary parameters c, c 1 , c 2 and ζ are not independent with each other, and satisfy a constraint The dual non-diagonal reflection matrix K + (u) is given by with the constraint In order to show the integrability of the system, we first introduce the "row-to-row" monodromy matrices T 0 (u) andT 0 (u) where {θ j , j = 1 · · · L} are the inhomogeneous parameters and L is the number of sites. The one-row monodromy matrices are the 3 × 3 matrices in the auxillary space 0 and their elements act on the quantum space V ⊗L . The tensor product is in the graded space, so we can write For the system with open boundaries, we need to define the double-row monodromy matrix which satisfies the similar relation as (2.9), in terms of matrix entries, they are Then the transfer matrix of the system is constructed as By using the (2.8), (2.9) and (2.10), we can prove the commutativity of t(u). (For further details about the commuting transfer matrix with boundaries for graded case, we refer the reader to [14,32,33]). The Hamiltonian (1.1) can be constructed by taking the derivative of the logarithm of the transfer matrix t(u) of the system with the parameters chosen as follows:

Nested algebraic Bethe ansatz
The block-diagonal structure of the K-matrix (2.14) permits us to use the nested algebraic Bethe ansatz to construct the associated Bethe state and obtain the eigenvalue as follows.
We first represent the double-row monodromy matrix Then the transfer matrix can be expressed by where k ± ij is the K ± matrix element in the ith row and jth column. Now we use the graded version of the nested algebraic Bethe ansatz method to obtain the eigenvalues of the transfer matrix (3.2). For this purpose, we first define the reference state |Ψ 0 as From the relations (2.21), (3.1) and (3.3), the elements of matrix T 0 (u) acting on the reference state |Ψ 0 give rise to The operators B 1 (u) and B 2 (u) acting on the reference state give nonzero values, and can be regarded as the creation operators of the eigenstates of the system. Following the procedure of the nested algebraic Bethe ansatz, the eigenstates of the transfer matrix can be constructed as where we have used the convention that the repeated indices indict the sum over the values 1,2, and F a 1 ...an is a function of the spectral parameters u j . Moreover, the coefficients where r ij = u + ηP ij , P α 1 α 2 β 1 β 2 = (−1) p(α 1 )p(β 2 ) δ α 1 β 2 δ β 1 α 2 with the grading p (1) = p (2) = 1, and where the corresponding eigenvalue Λ(u) is andΛ (u, {u j }) is the eigenvalue of the nested transfer matrixt (u, {u j }) given bŷ , (3.14) The vector components {F a 1 ...an } allow us to reconstruct the associated Bethe state (3.6), while the eigenvalueΛ (u, {u j }) gives rise to the associated eigenvalue (3.12) of the transfer matrix t(u) of the model. We shall determine the eigenvalueΛ (u, {u j }) and the corresponding eigenstate |F in the next section. The condition that the unwanted terms should be zero gives rise to that the M Bethe roots must satisfy the associated Bethe ansatz equations (BAEs) where Some remarks are in order. It is easy to check that the nested Bethe state |u 1 , . . . , u M ; F given by (3.6) and the eigenvalue Λ(u) given by (3.12) both have well-defined homogeneous limit (i.e., θ j → 0). This implies that in the homogeneous limit, the resulting Bethe states and the eigenvalue give rise to the eigenstate and the corresponding eigenvalue of the super t − J model described by the Hamiltonian (1.1).
Then it was shown in [36,37] that the eigenstate |F in (3.16) can be expressed as where the reference state |0 is provided that the parameters {w j |j = 1, . . . , M } satisfy the BAEs (4.24). The corresponding vector components {F a 1 a 2 ...a M } allow us to reconstruct the eigenstates |u 1 , . . . , u M ; F given by (3.6) of the original system 1 .

Concluding remarks
In this paper, we have studied the one-dimensional supersymmetric t − J model with the most generic integrable boundary condition, which is described by the Hamiltonian (1.1) and the corresponding integrable boundary terms are associated with the most generic nondiagonal K-matrices given by (2.14)-(2.15). By combining the algebraic Bethe ansatz and the off-diagonal Bethe ansatz, we construct the eigenstates of the transfer matrix in terms of the nested Bethe states given by (3.6) and (4.27), which have well-defined homogeneous limit. The corresponding eigenvalues are given in terms of the inhomogeneous T −Q relation (5.1) and the associated BAEs (5.2)-(5.3). The exact solution of this paper provides basis for further analyzing the thermodynamic properties and correlation functions of the model. These are under investigation and results will be reported elsewhere.