The Mermin-Wagner theorem has pointed out that there should be no spontaneous magnetization at finite temperature in 2D-isotropic Heisenberg model, or there will be no long-range ferromagnetic order. However, under classical Monte Carlo method with Metropolis configuration-updating strategy, I have found the equilibrium state at a low temperature exhibits the configuration shown in the attachment, in which a vortex structure similar to that in the XY model appears. Is this a normal result? Or is there something wrong with my program?
The configuation has been shown in the attachment. My code file has been added as well.