On Sat, 24 Jan 2009 12:00:02 +0100, comp-phys-alps-users-request wrote
- Re: uniform susceptibility (Matthias Troyer)
Message: 4 Date: Fri, 23 Jan 2009 09:21:06 -0800 From: Matthias Troyer troyer@phys.ethz.ch Subject: Re: [ALPS-users] uniform susceptibility To: comp-phys-alps-users@phys.ethz.ch Message-ID: 71F122F5-D527-4E3A-B308-8A86A23FD70C@phys.ethz.ch Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes
On 23 Jan 2009, at 08:58, Tadeusz Wasiutynski wrote:
Thank you for prompt answer. There is still something unclear for me. When running spinmc Ising model with say L=8 and T=1, square lattice, update cluster, sweeps 100000 I get: Magnetization^2 0.99862084 |Magnetization| 0.999285325 Susceptibility 63.911734 So it looks like Susceptibility=N * Magnetization^2 On the other hand given susceptibility at Tc scales very well with size. Where do I make mistake?
The "susceptibility" is measuredas N*Magnetizatio^2/T. Since T=1 here, you get what you observe
Matthias
In the other words it is not (<Magnetization^2>-<|Magnetization|>^2)/T as one could think. tadeusz