各位写过Ising Model模拟算法的大手们~大家画出来的比热和磁化率曲线和理论曲线相符吗?
在下用C写了一个算法, 用的是Metropolice算法, 但是画出来的比热和磁化率的曲线在温度0~1的时候有一个小峰, 如图所示. 各位请忽略峰位不太对, 这个图不是方格子的, 方格子的图一时没照出来, 但是问题是一样的...不同曲线是加了不同的外磁场.
苦思许久没有找出程序的错误, 窃以为是能量和磁化强度的涨落在低温时都很小, 相近的温度涨落值是很接近的, 从比热的计算公式
Cv~(-^2)/T^2可见, 分子很接近, 分母越小, Cv越大....T-T 磁化率也有相似的情况, 但是在低温时系统不应该有很大的比热(代表涨落)T^T
不知道大家遇到这个情况呢, 如果遇到了大家是怎么解决的? 还是只能是因为我的程序出错了?
谢谢大家了~~
![]() |
回复此楼
+3/873
马pi精
