This work is concerned with the quasi-classical limit of the boundary quantum inverse scattering method for the twisted sl (2|1) (2) vertex model with diagonal K -matrices. In this limit Gaudin’s Hamiltonians with diagonal boundary terms are presented and diagonalized.
Richardson–Gaudin models are a class of quantum integrable models connected to many physical systems, including pairing Hamiltonians from the theory of superconductivity. They can be obtained in the quasi-classical limit of the Quantum Inverse Scattering Method, which is based on an R-matrix and a Lax operator satisfying the Yang-Baxter equation. They can also be obtained from the Boundary Quantum Inverse Scattering Method, which relies on solutions of the reflection equations known as K-matrices. In this thesis we study these latter models systematically, explore the connections betweenthem and investigate the interpretation of the “boundary”.First of all, we consider Richardson–Gaudin models obtained from the spin-1/2 su(2) Boundary Quantum Inverse Scattering Method with diagonal K-matrices. We prove that the trigonometric boundary construction is equivalent to its rational limit, through a change of variables, rescaling, and a basis transformation. Moreover, we prove that the twisted-periodic and boundary constructions are equivalent in the trigonometric case, but not in the rational limit. Thus, including the “boundary” does not lead to a new model in this case.Next, we investigate Richardson–Gaudin models obtained from the spin-1/2 su(2) Boundary Quantum Inverse Scattering Method with non-diagonal K-matrices. Here the situation is different. The conserved operators in the boundary construction are no longer equivalent to the ones in the twisted-periodic construction. Also, the rational and the trigonometric boundary constructions are not equivalent. In the rational case this allows us to construct a generalisation of the p+ip pairing Hamiltonian with external interaction terms. In the trigonometric case the expressions for the conserved operators involve several free parameters, which can be adjusted to construct a variety of Hamiltonians. This result offers opportunities for future investigations.Finally, we study the case of the q-deformed bosonic Lax operator. This case is much more challenging than the case of the spin Lax operator. It is not straightforward to define the rational and quasi-classical limits of the bosonic Lax operator. Even after making modifications to the Lax operator for these limits to be well defined, it turns out that the limits do not commute. We state some open questions for future work.
We discuss a generalised version of Sklyanin's Boundary Quantum Inverse Scattering Method applied to the spin-1/2, trigonometric sl ( 2 ) case, for which both the twisted-periodic and boundary constructions are obtained as limiting cases. We then investigate the quasi-classical limit of this approach leading to a set of mutually commuting conserved operators which we refer to as the trigonometric, spin-1/2 Richardson–Gaudin system. We prove that the rational limit of the set of conserved operators for the trigonometric system is equivalent, through a change of variables, rescaling, and a basis transformation, to the original set of trigonometric conserved operators. Moreover, we prove that the twisted-periodic and boundary constructions are equivalent in the trigonometric case, but not in the rational limit.
Following Sklyanin's proposal in the periodic case, we derive the generating function of the Gaudin Hamiltonians with boundary terms. Our derivation is based on the quasi-classical expansion of the linear combination of the transfer matrix of the XXX Heisenberg spin chain and the central element, the so-called Sklyanin determinant. The corresponding Gaudin Hamiltonians with boundary terms are obtained as the residues of the generating function. By defining the appropriate Bethe vectors which yield strikingly simple off shell action of the generating function, we fully implement the algebraic Bethe ansatz, obtaining the spectrum of the generating function and the corresponding Bethe equations.