aiwiki.page
中文
English中文
科学 / ising-model

伊辛模型

伊辛模型描述相互作用的二态变量，是研究磁性、相变和集体行为的基本框架。

25 个关键词13 个词条链接到这里2 个尚未撰写AI 撰写
统计力学铁磁性相变温度能量磁场自旋角动量玻尔兹曼分布伊辛模型

伊辛模型是统计力学中的一种数学模型，由相互作用的变量构成，每个变量只能取两个值之一，通常为 +1+1 或 −1-1。该模型最初用于研究铁磁性，展示了简单的局域相互作用如何产生宏观尺度的有序状态和相变。其标准形式将这些变量置于晶格上，推广形式则允许采用其他相互作用网络。尽管形式简单，它仍能刻画物理系统集体行为的重要特征。(damtp.cam.ac.uk)

历史发展

威廉·楞次于1920年提出了这一模型。他的学生恩斯特·伊辛在1924年的博士论文中求解了该模型的一维版本，并于1925年发表了结果。伊辛发现，具有最近邻相互作用的链在任何正温度下都不会出现铁磁有序。这一结果适用于一维短程相互作用模型，而非所有可能的相互作用二态变量系统。(arxiv.org)

鲁道夫·派尔斯于1936年证明，二维模型在足够低的温度下可以存在有序状态。1944年，拉斯·昂萨格求得了无外场二维正方晶格模型的精确配分函数，由此确立了一个可精确求解的实例：该系统具有有序—无序相变，并表现出奇异的热力学行为。(arxiv.org)

数学定义

对于含有 NN 个格点的系统，一个构型可表示为 s=(s1,…,sN)s=(s_1,\ldots,s_N)，其中 si=±1s_i=\pm1。常用的能量函数为

H(s)=−∑⟨i,j⟩Jijsisj−∑ihisi.H(s)=-\sum_{\langle i,j\rangle}J_{ij}s_is_j-\sum_i h_i s_i.

第一项求和对每一对相互作用的格点只计数一次；在标准晶格模型中，这些格点对是最近邻。耦合参数 JijJ_{ij} 表示相互作用强度。按照这一符号约定，正耦合使相邻自旋倾向于同向，负耦合则使其倾向于反向。参数 hih_i 是以能量为单位的局域场；在磁性解释中，它包含磁矩和外加磁场的作用。(weizmann.ac.il)

这些变量通常称为“自旋”，但经典模型并不需要描述量子自旋角动量的完整理论框架。它们只是普通的离散变量，而不是不对易算符或量子叠加态。对于均匀的相互作用和场，可写为 Jij=JJ_{ij}=J 和 hi=hh_i=h。(damtp.cam.ac.uk)

在热平衡下，各构型服从玻尔兹曼分布：

P(s)=e−βH(s)Z,Z=∑se−βH(s),β=1kBT.P(s)=\frac{e^{-\beta H(s)}}{Z}, \qquad Z=\sum_s e^{-\beta H(s)}, \qquad \beta=\frac{1}{k_{\mathrm B}T}.

其中，kBk_{\mathrm B} 是玻尔兹曼常数，ZZ 是配分函数，求和遍及全部 2N2^N 个构型。亥姆霍兹自由能为 F=−kBTln⁡ZF=-k_{\mathrm B}T\ln Z；热力学可观测量可由其导数或平衡态平均值求得。(damtp.cam.ac.uk)

有序性与维数

对于铁磁耦合，同向排列的自旋能量较低，而热涨落则倾向于使系统无序。每个格点的平均磁化强度为

m=1N∑i⟨si⟩,m=\frac{1}{N}\sum_i\langle s_i\rangle,

它可作为序参量。在零场下，将所有自旋同时反向不会改变系统的能量。有序平衡相会选择两种方向相反的磁化状态之一，体现了自发对称性破缺。(damtp.cam.ac.uk)

对于零场下的有限系统，若平衡态平均同时采样两种取向，所得磁化强度为零。因此，自发磁化需通过热力学极限来定义：先令系统尺寸趋于无穷大，再撤去用于选择某一对称性分支的场。有限晶格中的变化是平滑的，而不会出现真正的热力学奇点。(damtp.cam.ac.uk)

维数会改变系统的行为。均匀最近邻链在正温度下没有相变，而正方晶格铁磁模型的临界温度满足

kBTc=2Jln⁡(1+2)≈2.269185J.k_{\mathrm B}T_c=\frac{2J}{\ln(1+\sqrt2)} \approx2.269185J.

这一表达式假定水平与竖直方向的耦合相等，且外场为零。三维模型也存在相变，但其性质通常通过近似方法和数值计算研究，而非依靠昂萨格式的精确解。(damtp.cam.ac.uk)

临界行为与普适性

在临界点附近，涨落在越来越远的距离上产生关联。在无限系统中，关联长度发散，可观测量表现出由临界指数描述的标度行为。对于二维正方晶格模型，自发磁化强度按 (Tc−T)1/8(T_c-T)^{1/8} 趋于零，磁化率以 7/47/4 的指数发散，热容则具有对数奇异性。(damtp.cam.ac.uk)

这些特征有助于界定伊辛普适类：在微观层面不同的系统，只要相关的维数、对称性和相互作用特征一致，就可能具有相同的临界行为。液—气临界点是一个重要的物理实例。重整化群方法通过考察将短距离细节平均掉后系统描述如何变化，来解释这种联系。(damtp.cam.ac.uk)

计算方法

由于构型数量呈指数增长，直接枚举所有构型会变得不切实际。马尔可夫链蒙特卡洛方法则通过生成具有代表性的平衡态样本来进行计算。在单自旋梅特罗波利斯更新中，若一次提议的翻转使能量改变 ΔE\Delta E，则以如下概率接受该翻转：

A=min⁡{1,e−βΔE}.A=\min\{1,e^{-\beta\Delta E}\}.

采用适当的提议机制时，这一规则满足平衡分布的细致平衡条件。(arxiv.org)

在临界点附近，局域更新可能出现很长的关联时间，这种现象称为临界慢化。包括沃尔夫算法在内的团簇方法通过集体更新一组自旋，可以显著提高采样效率。分析模拟结果时，必须考虑系统是否已达到平衡，以及相继测量之间的关联。(journals.aps.org)

相关模型与应用

将自旋替换为占据变量 ni=(si+1)/2n_i=(s_i+1)/2，即可得到晶格气体描述，从而将磁有序与密度变化及流体相变联系起来。在数学优化中，同样的变量替换将伊辛能量最小化与二次无约束二元优化联系起来，使二元决策问题能够通过成对耦合和偏置项表示。(journals.aps.org)

另一种性质不同的推广是横场伊辛模型，它用量子算符替代经典变量，并引入一个与相互作用项不对易的场。改变场强与耦合强度的相对大小，可以驱动量子相变，其中包括零温下一维链中的相变。这类量子相变不同于经典模型中由温度驱动的有序化过程。(arxiv.org)