Recent content by Berti

Thank you very much for all of this clarification and advice.

Sorry. The part of the code for the magnetization M is : sigmaxop = [] sites = [] for i in range (L): # \ sum_i sigma ^ z (i) sigmaxop.append ((sigmaz) .ToList ()) sites.append () And now I'd like to get the fluctuations <M^2> . So I tried to modify this part of the code to get...

I wanted to simulate the Ising model and it was okay until I wanted to get the fluctuations <M^2>. In fact, first, I wanted to obtain the magnetization M : $$ M = \sum_i \sigma ^ z_i $$ and it worked, I got indeed the magnetization M by writing these lines : sigmaxop = [] sites = []...

I am currently a student doing a Master in physics. I like the idea of having a forum where we can help each other, feel less alone and sometimes can "accelerate our knowledge" ... It also helps not to get stuck with a problem that seems a priori impossible while sometimes, an advice is enough...