The variational method is used to determine rotational-vibrational wave functions for HNS, NSH, and their ions from potential energy functions generated by high-level ab initio calculations, These wave functions describing the nuclear motion are then used in conjunction with ab initio dipole moment functions to predict radiative transition probabilities. Some spectroscopic constants are presented.