“…However, the case of nondiagonal boundary terms [3] has resisted solution for many years (see, e.g., [4]). A Bethe Ansatz solution for the latter case has recently been proposed in [5,6] (see also [7]). In terms of the parameters introduced there, the Hamiltonian is given by where σ x , σ y , σ z are the usual Pauli matrices, η is the bulk anisotropy parameter, α ± , β ± , θ ± are boundary parameters, and N is the number of spins.…”