あとなぜ無視したのか分からんが副格子仮定したうえでiPEPSなり使えば熱力学極限の計算ができる