Hi whatever-your-name-is, I cannot run your Mathematica notebook at the moment but here are two questions: - does your Mathematica notebook treat the ferromagnet? - does your Mathematica notebook use periodic boundary conditions like your ALPS I put or not? If the answer is yes to both questions then first compared the ground state energies of your code and ALPS Best regards Matthias Troyer

On Nov 21, 2014, at 13:04, <konefvitbka53@gmail.com> wrote:

First , I was obtained density of magnetization and density of energy curves of isotropic Heisenberg model on square lattice (2x2) using worm and dirloop_sse. After I was obtained these curves with wolfram mathematica. Curves for magnetization density(density of energy) obtained with mathematica and worm(sse) is strongly differ one from another.They must match. I was attached the files of mathematica and Heis2x2.py. <Model.zip>