Reaching the ground state of a large 2D classical spin model (e.g. classical Heisenberg model) might be a relatively difficult task while using conventional "flip/reject" Monte Carlo method. The system may be easily trapped in a local minima instead of the global minima, for example, the spin domains might appear while you are trying to get a pure ground state.
As far as know there are some techniques can help solve this issue includes "simulated annealing" and "parallel temperature", I would like to know if there are some other techniques might be helpful?
Secondly, many people use Langevin dynamics as an alternative to Monte Carlo method, but I have not tried it yet, I would like to know if the dynamics method has the same "local minima" issue as well?
Any comments and suggestions will be very welcome.