 Solution

Let $[a,a^{†}]=1$, $n=a^{†}a$, and take $S^{+}=(S^{-}{)}^{†}=\sqrt{2S-n}a$. The physical Holstein–Primakoff occupation constraint limits the Hilbert space, which contains only $n=0,1,\dots ,2S$. Acting in this number basis givesThe upper and lower endpoints are annihilated as required. ConsequentlyThe change in $n$ also gives $[S^z,S^{\pm }]=\pm S^{\pm }$. These matrix elements prove the spin algebra on the full physical spin-$S$ space. Extending the oscillator space to arbitrary occupation would not give a physical finite spin: the occupation constraint and operator ordering in the square root matter. bosons on distinct sites have $[a_m,a_n^{†}]={\delta }_{mn}$ and all other site commutators zero.