Ising Model: Wärmekapazität und Suszeptibilität

mzh · Anmeldungsdatum: 31.05.2010 Beiträge: 44

Hallo zusammen
Wir sollen da ein Ising Model programmieren. Läuft eigentlich auch recht gut, jetzt sollen wir zeigen, wie die Wärmekapazität und die Suszeptibilität am kritischen Punkt (beim Phasenübergang) divergiert.
Für die Wärmekapazität ist folgende Formel gegeben:

Hier ist N die Seitenlänge des Spingitters (= 20), J die Kopplungskonstante (=-1.0), E die Energie und kB die Boltzmannkonstante (=1.0). Dafür kriege ich auch den Phasenübergang (siehe unteres Bild).

Nun sollen wir noch die Suszeptibilität über die Temperatur auftragen.
Die Formel, die uns gegeben wurde ist

Hier ist M die Magnetisierung des Spingitters

, wobei S der Betrag des Spins ist. N ist wieder die Seitenlänge des Gitters.
Wenn ich das nun ebenfalls gegen die Temperatur auftrage erhalte ich das obere Bild (sorry wegen der Reihenfolge, kommt vom hochladen).
Ich hab mir das mal ein bisschen an anderen Stellen angeschaut, die Wärmekapazität scheint qualitativ zu stimmen, aber die Suszeptibilität sollte glaub positiv divergieren. Kann das jemand bestätigen? Haben wir eine falsche Formel? Hab leider keine Quelle für die Herkunft der Formel.

Besten Dank für Hinweise.