This is the Blume–Capel model. In the mean-field approximation, write and replace
Because every site has coordination number , the resulting energy is
The one-site partition function is therefore
With , the mean-field partition function is . Taking gives
Write the spin-one-half Heisenberg antiferromagnet as
on a bipartite lattice of coordination number . A one-spin density operator is , where is its Bloch vector and . In a two-sublattice product state,
The minimum is , attained by pure antiparallel Bloch vectors, so the minimizing product state is a Néel state. In the limit this mean-field approximation becomes exact and the ground-state energy per bond is therefore
The phrase “energy density” requires a coupling convention. For the unscaled Hamiltonian above, every site belongs to bonds and
which diverges as . With the standard Kac normalization
the finite energy density is
If the convention divides by the spatial dimension instead, the answer is per site. A Hamiltonian written with rather than multiplies all these energies by four.