Matlab蒙特卡洛模拟二维伊辛模型相变过程

icfg66

最新推荐文章于 2026-10-10 16:53:42 发布

关注

VIPMatlab蒙特卡洛模拟二维伊辛模型相变过程

一、什么是伊辛模型伊辛（Ising）模型是描述磁系统相变最简单的模型，但模型里自旋之间的相互作用赋予了它奇妙的特性，最有趣的就是对称性破缺。这一模型可以被推广用于研究连续的量子相变、基本粒子超弦理论、动力学临界行为等，甚至被认为可以描述深林火灾、交通拥堵、舆论传播等社会经济现象。如图，每个格点的方向只有向上或向下两者状态，但临近的自旋之间有相互作用，而且点阵可以是一维、二维、三维、甚至更高维...

开通可读全文 ·2.2w 人阅读

最低0.47元/天解锁文章

一、什么是伊辛模型

伊辛（Ising）模型是描述磁系统相变最简单的模型，但模型里自旋之间的相互作用赋予了它奇妙的特性，最有趣的就是对称性破缺。这一模型可以被推广用于研究连续的量子相变、基本粒子超弦理论、动力学临界行为等，甚至被认为可以描述深林火灾、交通拥堵、舆论传播等社会经济现象。
伊辛模型介绍
如图，每个格点的方向只有向上或向下两者状态，但临近的自旋之间有相互作用，而且点阵可以是一维、二维、三维、甚至更高维，这两个特点让伊辛模型的严格求解成为了世纪难题。为了定量描述这个系统的能量，我们假设第 $i$ 个格点的自旋为 $s_i$ ， $s_i$ 只能取+1或-1，如果相邻两个格点同方向，则它们相互作用的能量更小，设为 $- J$ ，如果为反方向，则为 $J$ ， $J$ 称为耦合系数，通常为正值，代表铁磁系统，如果 $J$ 为负值，则代表反铁磁系统。如果外磁场的强度为B，格点的自旋磁矩为 $\mu$ ，那么可以写出这个体系的哈密顿量：
$H=-J\sum_{<ij>}s_i s_j-\mu B \sum_{i=1}^N s_i$

自发对称性破缺

我们先定性地了解一下这个系统的性质，令外磁场零，当温度 $T\rightarrow 0$ 时，体系为了保持能量最低，所用的格点趋向于同方向，系统整体要么向下，要么向上，呈现强磁性。当温度 $T\rightarrow \infty$ 时，系统热运动占主导地位，格点方向呈现随机性，系统整体不带磁性，从上或从下观察体系，呈现出对称性，或者说无法通过系统磁性区分上或下。现在再考虑，当温度T从 $\infty$ 逐渐降温，那么系统必定存在某个温度 $T_c$ ，高于此温度时系统无磁性，低于此温度时，系统磁性逐渐加强。这个温度就是临界温度，也是相变点，系统从对称磁体转变为非对称磁体，而这就是对称性破缺，因为这种破缺不是外界扰动（如外加磁场）引起的，而是由内部的关联作用力造成的，所以称之为自发的对称性破缺。

上图就是对自发对称性破缺的定性描述，当逐渐降温到 $T_c$ 时，系统磁性开始出现分化，要么向下要么向上，最终平均的磁化强度 $\overline{s}$ 趋向于+1或-1。

我们感兴趣的问题主要有两个：第一，不同维度、不同分布的格点，其临界温度 $T_c$ 是多少；第二， $T_c$ 附近， $\overline{s}$ 随温度 $T$ 以什么样的幂指数 $\alpha$ 趋近于 $T_c$ ：
$\overline{s}\sim(T_c-T)^{\alpha}$

统计物理的思路

我们先考虑简单的9个格点的例子，实际格点数的量级为 $10^{29}$ 。
简单例子
假设格点耦合强度 $J = 1$ ，那么这个9格点体系的能量为：
$E = [(- 1 + 1) + (- 2 + 1) + (2) + (2 - 1) + (- 3 + 1) + (- 1 + 2) + (- 1 + 1) + (- 1 + 2) + 2] / 2 = 2$
第一个格点有两个相邻格点，右边的与其同向，耦合能为-1，下面的与其反向，耦合能为1；第二个格点有三个相邻格点，左和下与其同向，耦合能为-2，右边与其反向，耦合能为1。类似的可以计算其余格点的耦合能。最后除以2是因为每个相互作用重复计算了一次。而这一特定分布以某概率P出现：
$P\varpropto e^{-E/kT}$