After a bipartite spin rotation, a leading Holstein–Primakoff transformation expansion gives . The bosonic Bogoliubov diagonalization gives the displayed spectrum and vacuum energy . The anisotropy gap is . At the finite-angle transformation is singular at its zero modes, and quantum depletion of Néel order is infrared divergent in one dimension.
With unit lattice spacing and , a nearest-neighbour Heisenberg model has quadratic low-momentum ferromagnetic magnons and linear antiferromagnetic modes in the linear spin-wave approximation. The ferromagnetic polarized vacuum is an exact ground state, whereas the assumed antiferromagnetic ordered vacuum has quantum depletion of Néel order. In one dimension its logarithmic infrared divergence invalidates finite ordered magnetization and prevents the linear approximation from deciding the exact gap.
Hypercubic antiferromagnetic spin-wave dispersion Created 2026-10-07 Updated 2026-10-20
For a nearest-neighbour Heisenberg antiferromagnet on a hypercubic lattice, the linear spin-wave approximation gives and the displayed frequency. Near a gapless point its spin-wave velocity is . Its ground-state boson occupation is , producing quantum depletion of Néel order. Exact zero modes need an infrared regulator or separate global-rotation treatment.
Linear spin-wave approximation Created 2026-10-06 Updated 2026-10-18
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- excitation energies above the order- 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 vacuum of the diagonal quasiparticles contains original spin bosons: . Thus the quantum depletion of Néel order gives
In one dimension . Removing the global zero modes with an infrared cutoff, the continuum integral contains near both gapless points and diverges logarithmically as the cutoff is removed. This is a breakdown of the assumed Néel order reference, not a physical infinitely negative magnetization. It indicates that the one-dimensional ordered spin-wave expansion cannot maintain a finite staggered moment; it does not determine the exact spectral gap or all properties of the spin chain for arbitrary .
At zero temperature in dimensions the singular contribution scales as . It is infrared finite for , permitting a finite quantum reduction and a self-consistent ordered spin-wave description in an appropriate regime. This differs from the positive-temperature contribution, where yields . Its divergence for agrees with the Mermin-Wagner theorem for short-range continuous-symmetry models. Finite-temperature nonordering and the zero-temperature one-dimensional depletion argument are separate statements.
Under the bipartite spin rotation, the staggered magnetization becomes the uniform transformed . In the Bogoliubov transformation vacuum ,
The quantum depletion of Néel order therefore gives, with the zero modes regulated,
Near each zero of , the integrand behaves as , producing a logarithmic divergence. With a finite-size cutoff of order , the depletion grows as . Thus the large- 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.
Use dimensionless spin operators and nearest-neighbour exchange on a regular bipartite graph with N sites and coordination number of a lattice z. Count each bond once and define , real for the inversion-symmetric examples considered here.
For the ferromagnetic Heisenberg model, . Expanding the Holstein–Primakoff transformation about all spins along +z gives
up to interactions of higher order in boson density, with . The isotropic spin-wave vacuum energy comparison starts with this normal ordering. The linear spin-wave approximation then gives
Each ferromagnetic magnon lowers total by one. The fully aligned state is an exact eigenstate and ground state; there is no harmonic vacuum depletion or added zero-point energy in this normal-ordered Hamiltonian. The global spin-rotation multiplet remains degenerate. On a hypercubic lattice of spacing a, and , so . The exact one-magnon sector agrees with this result even at finite S; interactions enter for multiple magnons.
For the Heisenberg antiferromagnet, take and expand around a Néel state. On A, and ; on B, and . A bond becomes
With N/2 momenta in the magnetic Brillouin zone,
The bosonic Bogoliubov transformation has and . For example and . It gives two branches in this reduced zone,
The reduced-zone counting of bipartite spin waves gives N oscillators in total, not 2N. The negative vacuum correction and quantum depletion of Néel order, , distinguish the antiferromagnet from the ferromagnet. At nonzero temperature replace by . Hypercubic small-momentum modes have , hence linear physical frequency .
Self-consistency requires depletion much smaller than S and respects the Holstein–Primakoff occupation constraint. In one dimension the antiferromagnetic vacuum depletion diverges logarithmically, invalidating ordered spin-wave theory; it does not alone prove an exact gap or its absence. At positive temperature, isotropic short-range models in d≤2 cannot sustain such magnetic order by the Mermin-Wagner theorem. In three dimensions the quadratic ferromagnetic modes give the low-temperature magnetization reduction proportional to , whereas the antiferromagnetic thermal reduction scales as . Anisotropy, fields, frustration and magnon interactions alter these conclusions. [ref]
The Holstein–Primakoff transformation converts a spin algebra into bosonic operators while retaining the finite physical spin Hilbert space. With , write
The Holstein–Primakoff occupation constraint is . Acting on , the raising matrix element is and the lowering element is . Their squared difference is , reproducing . Operator ordering inside the square roots is essential. For low boson occupation relative to , expand them to obtain and . Keeping quadratic terms gives the linear spin-wave approximation; higher terms describe magnon interactions. Large helps control the expansion, but a low-dimensional infrared divergence can invalidate it even at large finite .
For the ferromagnet choose , , with unit lattice spacing and . Expand around the all-up state. A longitudinal bond contributes , and the transverse exchange contributes . After a Fourier transform,
These are ferromagnetic magnons. The polarized vacuum is an exact ground state of the full ferromagnetic Hamiltonian: every bond is in its maximum total-spin sector, minimizing that bond energy. It therefore has no quantum zero-point reduction of its magnetization. The one-magnon dispersion is also exact, although interacting many-magnon states need not be free bosons. Rotating the polarization generates the maximum-total-spin ground-state multiplet. Selecting a definite direction in a thermodynamic symmetry-breaking description breaks to ; the two transverse spin components form a canonical pair, producing a single quadratic type-B Goldstone boson. A small aligning field gaps it by the Zeeman energy.
At nonzero temperature the one-dimensional ferromagnet is different: its thermal magnon number contains , with infrared behaviour . The divergent magnetization reduction signals the absence of finite-temperature long-range order in an infinite short-range isotropic chain, consistent with the Mermin-Wagner theorem. It does not disallow the polarized zero-temperature ground states just constructed.[ref]
For the antiferromagnet use , still , on an even bipartite chain. The classical reference state alternates up and down. Rotate one sublattice through so that both local quantization axes point along their reference spins. In unrotated physical components, to leading order , , and . Choosing Fourier phases using the actual site positions gives
The reduced Brillouin zone has momenta and length . The pair creation and annihilation terms express the fact that an antiparallel reference product state is not an eigenstate. A bosonic Bogoliubov transformation with , , diagonalizes each two-mode block. With the sign convention and , the result is
Thus the harmonic antiferromagnetic modes have linear, rather than quadratic, long-wave dispersion. In an unfolded Brillouin-zone description the same result is represented by zeros at both zero and the antiferromagnetic wavevector; these are not extra independent degrees of freedom beyond the two-sublattice accounting.
The Bogoliubov vacuum has no or excitations but contains original spin bosons, with . Its staggered magnetization is . This quantum depletion of Néel order diverges logarithmically at small in one dimension. Linear spin-wave theory therefore does not establish a stable ordered antiferromagnetic ground state for this chain. Its nominal gapless dispersion is a result of expanding around an invalid ordered reference, not a general proof of the exact excitation spectrum. For the nearest-neighbour isotropic chain, integer-spin systems have a Haldane gap, whereas half-integer-spin systems are gapless; the spin-one-half chain in particular has fractional spin excitations rather than a simple dilute gas of ordered-state magnons. This distinction requires physics beyond the harmonic expansion, including the topological term in the continuum theory.[ref]
Figure 1.
Original harmonic chain dispersions; the antiferromagnetic curve is a spin-wave approximation, not the exact spectrum for every spin
. Details.
Two sublattices give two bosonic modes for each of N/2 wavevectors in a magnetic Brillouin zone. For a symmetric quadratic pairing block with number coefficient A and anomalous coefficient B, a bosonic Bogoliubov transformation gives two equal positive energies and the displayed vacuum shift. There are N modes overall. The quantum depletion of Néel order per site is , where . Equivalently one can use one branch over a full zone with the corresponding factors of one half. Mixing these conventions doubles the vacuum energy or depletion.