= Solution
Under the <bipartite spin rotation>, the <staggered magnetization> becomes the uniform transformed $S^z$. In the <Bogoliubov transformation> vacuum $\alpha_k|0_\alpha\rangle=0$,
$$
\langle a_k^\dagger a_k\rangle=\sinh^2\theta_k=\frac12\left(\frac A{\omega_k}-1\right)=\frac12\left(\frac1{|\sin k|}-1\right).
$$
The <quantum depletion of Néel order> therefore gives, with the zero modes regulated,
$$
\boxed{M_s=S-\frac1N\sum_k\sinh^2\theta_k\ \longrightarrow\ S-\int_{-\pi}^{\pi}\frac{dk}{2\pi}\frac12\left(\frac1{|\sin k|}-1\right)}.
$$
Near each zero of $\sin k$, the integrand behaves as $1/(2|k-k_0|)$, producing a <logarithmic divergence>. With a finite-size cutoff of order $N^{-1}$, the depletion grows as $\pi^{-1}\log N+O(1)$. Thus the large-$S$ expansion about a state with finite <Néel order> is not self-consistent in the infinite one-dimensional chain. The divergent expression is not a negative physical magnetization; it signals breakdown of that ordered approximation. It does not determine whether the exact excitation spectrum is gapped.
Back to article page
