 Solution

Inside the infinite square well, the orbital and spin operators act on different factors of the tensor product. Write $\mathbf{S}={\mathbf{S}}_1+{\mathbf{S}}_2$. The identityshows that the exchange interaction commutes with $S^2$ and $S_z$. The sum of the one-particle kinetic energies also commutes with them and with particle exchange. Thus one may use orbital quantum numbers together with the total-spin sector $S=0,1$ and projection $M=-S,\dots ,S$, with physical $S_z=\hslash M$.

For identical spin one-half particles the Pauli exclusion principle further restricts the states. The spin-one-half singlet state is antisymmetric in spin and must multiply a symmetric spatial wavefunction; a spin-one-half triplet state must multiply an antisymmetric spatial wavefunction. Thus $n_1,n_2$ label unordered orbital occupations in symmetry-adapted states, not arbitrary ordered product states. This gives the labels $|n_1,n_2,S,M⟩$, with the paper's last label expressed in units of $\hslash$.