Quantum integrable matrix models of spinor Bose gases in one spatial dimension

Abstract: Degenerate spinor Bose gases with repulsive density-density interaction and anti-ferromagnetic spin-spin coupling in one spatial dimension are shown to be described by a quantum integrable matrix extension of the nonlinear Schrödinger model, whose fundamental fields are described by an $m\,\times\,n$ matrix of bosonic field operators. The eigenstates of this model are constructed for arbitrarily sized matrix field operators by means of algebraic Bethe-ansatz techniques, and the corresponding Bethe equations governing the spectra of conserved quantities are derived. The approach thus generalizes previously chosen techniques to account for arbitrary spin multiplets and their spin-spin interaction. Focusing on the specific case of the $2\times2$ model, which is shown to correspond to a spin-$1$ Bose gas, a set of integral equations is derived, which describe its equilibrium thermodynamic properties. From these, the ground state phase diagram is computed both, numerically and analytically in the parameter plane spanned by the chemical potential and an external magnetic field. Furthermore, the existence of paired bound states is shown to modify the Pauli exclusion principle for interacting bosons in one dimension. In particular, it is found that no two quasiparticle rapidities can coincide, provided that the Lieb parameter satisfies $γ>4/3$.