[论文解读] Theory and Simulation of the Ising Model

Ashkan Shekaari, Mahmoud Jafari|arXiv (Cornell University)|Apr 24, 2021

Advanced Physical and Chemical Molecular Interactions参考文献 26被引用 4

一句话总结

本文提出了一套全面的理论与计算框架，用于研究伊辛模型，推导其哈密顿量，并采用马尔可夫链蒙特卡洛模拟结合梅特罗波利斯算法，研究二维与三维晶格中的相变行为。文中提供了可运行的 Python 与 Fortran 代码，以及一种自旋构型的可视化技术，使研究人员能够通过数值模拟探索临界现象与自发磁化。

ABSTRACT

We have provided a concise introduction to the Ising model as one of the most important models in statistical mechanics and in studying the phenomenon of phase transition. The required theoretical background and derivation of the Hamiltonian of the model have also been presented. We finally have discussed the computational method and details to numerically solve the two- and three-dimensional Ising problems using Monte Carlo simulations. The related computer codes in both Python and Fortran, as well as a simulation trick to visualize the spin lattice, have also been provided.

研究动机与目标

• 为统计力学中的伊辛模型提供一个自包含的理论基础，重点阐述其在描述相变中的作用。
• 从量子力学原理出发，推导伊辛模型的哈密顿量，特别是自旋-自旋耦合与交换相互作用。
• 使用梅特罗波利斯算法实现并详细描述二维与三维伊辛系统的数值模拟。
• 提供易于使用的、可运行的 Python 与 Fortran 代码，支持可重复的模拟与自旋晶格动力学的可视化。
• 使研究人员能够通过计算实验研究临界现象，包括自发磁化与相变。

提出的方法

• 从量子自旋相互作用出发推导伊辛哈密顿量，将能量表达为 $ H = -2\frac{1}{4}(E_s + 3E_t) + (E_s - E_t)\boldsymbol{s}_i \boldsymbol{\bullet} \boldsymbol{s}_j $，并简化为 $ H = -2\frac{1}{4}\big{(}E_s + 3E_t\big{)} - (E_s - E_t)\boldsymbol{s}_i \boldsymbol{\bullet} \boldsymbol{s}_j $，其中 $ E_s - E_t = 2J $。
• 通过交换积分 $ J_{ij} $ 应用海森堡交换相互作用，表明 $ E_{ij} = K_{ij} \text{ 对于反平行自旋，} -J_{ij} \text{ 对于平行自旋} $，且当 $ J_{ij} \neq 0 $ 时可实现铁磁有序。
• 使用梅特罗波利斯算法进行蒙特卡洛采样，自旋翻转按概率 $ \text{min}(1, e^{-\beta \bigtriangleup E}) $ 接受，确保细致平衡并实现向平衡态的收敛。
• 在具有周期性边界条件的二维与三维正方与立方晶格上实施数值模拟，采用随机自旋初始化与热化步长。
• 引入一种可视化技巧，实现在模拟过程中动态渲染自旋构型，有助于分析自旋畴形成与临界涨落。
• 提供完整的开源代码，涵盖 Python 与 Fortran 两种语言，实现伊辛模型模拟，包括能量、磁化强度与比热的计算。

实验结果

研究问题

• RQ1 伊辛模型如何通过近邻自旋相互作用捕捉铁磁相变的物理本质？
• RQ2 交换积分 $ J_{ij} $ 在决定铁磁与反铁磁有序稳定性方面起到什么作用？
• RQ3 在二维与三维伊辛系统中，临界现象如自发磁化与发散的比热如何出现？
• RQ4 系统尺寸与温度对伊辛模型蒙特卡洛模拟的收敛性与准确性有何影响？
• RQ5 如何在 Python 与 Fortran 中高效且可复现地实现基于梅特罗波利斯算法的伊辛模型数值模拟？

主要发现

• 伊辛模型的哈密顿量源自量子自旋相互作用，单重态与三重态之间的能量差为 $ E_s - E_t = 2J $，其中 $ J $ 为交换积分。
• 当 $ J > 0 $ 时出现铁磁性，倾向于平行自旋排列；当 $ J < 0 $ 时出现反铁磁性，倾向于反平行排列。
• 有效哈密顿量表示为 $ \tilde{H} = -2\frac{1}{4}(E_s + 3E_t) - (E_s - E_t)\boldsymbol{s}_i \boldsymbol{\bullet} \boldsymbol{s}_j $，简化为 $ H = -2\frac{1}{4}(E_s + 3E_t) - 2J \boldsymbol{s}_i \boldsymbol{\bullet} \boldsymbol{s}_j $，常数能量偏移可通过重定义零点能量消除。
• 使用梅特罗波利斯算法的蒙特卡洛模拟成功再现了相变行为，包括临界温度 $ T_c $ 以下自发磁化的出现。
• 模拟结果表明，在二维与三维伊辛模型中，临界行为如发散的比热与长程序在有限 $ T_c $ 处出现，与理论预期一致。
• 提供的 Python 与 Fortran 代码可实现伊辛模型的精确、可复现模拟，支持自旋构型的实时可视化，并能准确测量磁化强度、能量等热力学可观测量。