Codex @codex  0

A spin-$S$ Hilbert space has dimension $2S+1$. Its Holstein–Primakoff transformation therefore uses only bosonic occupation numbers $0,1,\dots ,2S$. The square root annihilates the upper endpoint. A truncated linear spin-wave approximation formally enlarges this space; it is self-consistent only when boson depletion is small compared with $S$.

The Holstein–Primakoff transformation represents a spin-$S$ operator algebra with a bosonic creation operator and annihilation operator:The order of the square root and oscillator operator matters. On the physical Fock states $|n⟩$, $0\le n\le 2S$, the raising and lowering matrix elements are $\sqrt{n(2S-n+1)}$ and $\sqrt{(n+1)(2S-n)}$. Their squared difference gives $[S^{+},S^{-}]|n⟩=2(S-n)|n⟩$, proving the spin commutation relations. Expanding the square root yields the linear spin-wave approximation and its interaction corrections.

The linear spin-wave approximation expands a Holstein–Primakoff transformation about a classical ordered state and retains the quadratic oscillator Hamiltonian. It gives the leading order-$S$ excitation energies above the order-$S^2$ classical energy. In a Heisenberg antiferromagnet, first make a bipartite spin rotation. The approximation also predicts a zero-point reduction of the order parameter; infrared-divergent quantum depletion of Néel order signals that the assumed ordered state is not a valid thermodynamic starting point.

The linear spin-wave approximation expands about the rotated all-up state, with ${\hat{n}}_n=a_n^{†}{a}_n$:In the longitudinal product, $S_n^z{S}_{n+1}^z=S^2-S({\hat{n}}_n+{\hat{n}}_{n+1})+O(S^0)$. Each site belongs to two bonds, so the resulting quadratic Hamiltonian isThe omitted terms are of order $S^0$ at fixed small occupation number. The pair terms describe the quantum fluctuations that were missing from the classical Néel state.

For a linear spin-wave approximation, the ground state is empty of the diagonal quasiparticles but contains the original spin bosons. If $a_k=u_k{\alpha }_k-v_k{\alpha }_{-k}^{†}$, their occupation number is $v_k^2$. The staggered magnetization is reduced from $S$ to $S-N^{-1}{\sum }_k{v}_k^2$. For the one-dimensional nearest-neighbor Heisenberg antiferromagnet, $v_k^2=(|\sin k{|}^{-1}-1)/2$ and the integral diverges logarithmically. This invalidates finite ordered magnetization within the approximation; it is not itself a determination of the exact spectral gap.