标签：---- Metropolis -++- +-++ 组态 ++++ Ising
Ising模型哈密顿量为
\[\hat{H} = - J \sum_{i,j相邻} s_i s_j - B \sum_i s_i \]若取 \(J=1, B=1\)，做 4 x 4 一共 16 个 spin 的 Ising 模型模拟，可以得到很平缓的相变曲线，如下图
曲线最左侧有些点低于低于1.0，增大 Markov 链的根数，增大链的长度，都没有改善这一点。所以仔细输出中间结果仔细观察了一下，发现温度特别低的时候存在一些有趣的（不好的）情况。
1. 组态阻塞
做统计模拟的人大概喜欢叫做微观状态，我习惯了组态这个词，就这么用了。如果 Markov 链产生如下组态
----
++++
++++
----
那么，翻转任何一个spin，都会导致体系上升。所以这实际上是一个模型的亚稳态，在温度较低的时候，e.g. kT=0.1，接受概率是 \(e^{-2/0.1} = 2e-9\)，所以翻不过去。这样的组态比基态温度高很多，按照Metropolis模拟的思想，这样的组态出现概率理应极低。但是，因为低温下的阻塞，Metropolis模拟的精神得不到体现。
2. 组态振荡
最低能组态 E = -32，下面这些 E = -8 的组态在低能下会不断循环。
组态1：
-++-
++++
++-+
----
组态2：
++-+
++-+
+---
++-+
组态3:
++++
+-++
---+
+-++
组态4:
+-++
----
+-++
++++
组态5:
-++-
--++
-+++
-++-
组态6:
----
-++-
+++-
++--
总而言之，4 x 4 的格子做低温 kT/J = 0.1, ..., 0.5 的时候，会有阻塞和循环现象，导致结果中有些非物理的内容。
自然有如下问题：
- 更大的格子上有没有类似的现象？我猜是有。
- 一个重要的问题是，不同能量的组态数时怎么样的？或者说，态密度是怎么样的？这个问题的答案对于 Ising Model 的结果无疑是很重要的。
标签：----,Metropolis,-++-,+-++,组态,++++,Ising 来源： https://www.cnblogs.com/luyi07/p/15519317.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。