I have read a few discussions about the difference between the Heisenberg model (using quantum spin operators) and Ising model (with spins $\pm 1$), notably this one or this Quora post.
All the responses seem to be of theoretical nature, highlighting, e.g., the difference in symmetry or the Mermin-Wagner theorem stating the Heisenberg model does have a phase transition at $d\leq 2$.
I am interested in practical aspects. Imagine I want to simulate some real material (let us consider a $3d$ lattice only). Even an undergraduate can write and run a Monte Carlo code simulating the Ising model a few hours. On the other hand, the Heisenberg model involves a true effort. Writing it requires immense knowledge, intuition and skill (constructing the Hamiltonian and total spin operators and computing the expectation values), and running it demands powerful computers even for a fairly low number of atoms and approximation schemes (see this review of methods).
Given the amount of effort, what is the justification for using the quantum Heisenberg model? What new vital features observed in experiment does it provide? If any, can they not be reproduced by the Ising model with vectorised but classical 3D spins, $-J \sum_{\langle ij\rangle} \mathbf S_i \cdot \mathbf S_j$?
Edit: To clarify that what I am trying to figure out is different from this question, I am not interested in theoretical features, but in practical/experimental results in three dimensions. That is, to what extent do Heisenberg/Ising models correctly predict:
- critical exponents of magnetisation, heat capacity and susceptibility,
- collective excitations (can these be captured by 3d Ising model?),
for a range of ferromagnetic materials, and where do they fail?