咋用函数计算2D Ising model的总能量

呵呵,好多年前好像也碰到过这样的题,可惜当时懒,没有写。如果是简单情形的2D Ising model 比如正方格子情况以及只考虑最近邻的相互作用,那就直接根据各个格点的自旋方向和最近邻的交换能加和就行了。 初始的构型如果没给直接用基态好了。然后考虑温度会改变spin方向, 就把温度的激发造成的自旋的转向一个个在各个格点更新一下状态,数目,大小和位置就通过统计抽样给出。你这个情况是温度变不变? 如果没搞错的话用正则系综。如果是三角格子,由于有spin frustration(反磁性情况下在相邻三个格点无法满足)就麻烦一点,我也不知道是怎么处理得到平衡态的。有空动手试一下看看怎么样的,当然也可以参考一下文献,搞spin glass啥的人天天就玩这个。


    推荐阅读