2次元イジングモデルのオンサーガー解の数値計算

無限体積の2次元イジングモデルの厳密解であるオンサーガー解を数値的に計算する．

1スピン当たりの自由エネルギー \begin{align} -\beta f &= \log(2\cosh 2K) + \frac{1}{2\pi ^ 2} \int ^ \pi _ 0 d\omega _ 1 \int ^ \pi _ 0 d\omega _ 2 \log(1 - k ^ 2 _ 1 \cos\omega _ 1 \cos \omega _ 2) \\ k ^ 2 _ 1 &= \frac{2\sinh 2K}{\cosh ^ 2 2K} \end{align} をpythonで数値計算してエネルギーと比熱を求める．逆温度 0.2 ≦ K ≦ 0.7 を 0.0005 刻みで計算する．

Onsager

参考

相転移・臨界現象の統計物理学 新物理学シリーズ