易辛模型-族谱新闻-族谱网

易辛模型

2020-10-16

出处：族谱网

作者：阿族小谱

浏览:248次

转发:0次

评论:0

定义令Λ为所有晶格点的集合，其中每个晶格点都有一个所有和它相邻的晶格点的集合(在数学上称之为图)并使这些晶格点形成一个d维的晶格。对于每个晶格点k∈Λ都有一个离散变数σk，其中σk∈{+1,−1}，代表一个晶格点的自旋。而所有变数的集合σ=(σk)k∈Λ则称作自旋组态。对于两个相邻的晶格点i,j∈Λ，我们可以引入一个交互作用参数Jij，此外，我们可以假设每个自旋j∈Λ都和外加的磁场hj作用。则整个系统的哈密顿量可写成:其中代表晶格点i和晶格点j是相邻的的晶格点。因此哈密顿量的第一项为对每一对相邻晶格点的总和(每一对只算一次)，代表所有自旋之间交互作用的能量，而第二项则是磁场和自旋交互作用的能量。µ是晶格点磁矩的值，值得注意的是，电子的磁矩和他的自旋方向相反，所以哈密顿量的第二项应该要是正号比较合理，但在习惯上，还是会令第二项为负号。该系统的的组态概率P(σ)为在热平衡下某个特定自旋组态σ的...

定义

令Λ为所有晶格点的集合，其中每个晶格点都有一个所有和它相邻的晶格点的集合(在数学上称之为图)并使这些晶格点形成一个d维的晶格。对于每个晶格点 k ∈ Λ 都有一个离散变数σ k ，其中 σ k ∈ {+1, −1}，代表一个晶格点的自旋。而所有变数的集合σ = (σ k ) k ∈Λ 则称作自旋组态。

对于两个相邻的晶格点 i , j ∈ Λ ，我们可以引入一个交互作用参数 J ij ，此外，我们可以假设每个自旋 j ∈ Λ都和外加的磁场h j 作用。则整个系统的哈密顿量可写成:

其中代表晶格点 i 和晶格点 j 是相邻的的晶格点。因此哈密顿量的第一项为对每一对相邻晶格点的总和(每一对只算一次)，代表所有自旋之间交互作用的能量，而第二项则是磁场和自旋交互作用的能量。µ是晶格点磁矩的值，值得注意的是，电子的磁矩和他的自旋方向相反，所以哈密顿量的第二项应该要是正号比较合理，但在习惯上，还是会令第二项为负号。

该系统的的组态概率 P(σ)为在热平衡下某个特定自旋组态 σ 的概率，为波兹曼分布:

P β β --> ( σ σ --> ) = e − − --> β β --> H ( σ σ --> ) Z β β --> , {\displaystyle P_{\beta }(\sigma )={e^{-\beta H(\sigma )} \over Z_{\beta }},}

其中 β = ( k BT ) ，而:

Z β β --> = ∑ ∑ --> σ σ --> e − − --> β β --> H ( σ σ --> ) {\displaystyle Z_{\beta }=\sum _{\sigma }e^{-\beta H(\sigma )}}

是该概率分布的归一化常数，在统计力学中又称做配分函数。对于有为自旋组态函数的物理量f(σ) ，其期望值可表示为:

⟨ ⟨ --> f ⟩ ⟩ --> β β --> = ∑ ∑ --> σ σ --> f ( σ σ --> ) P β β --> ( σ σ --> ) {\displaystyle \langle f\rangle _{\beta }=\sum _{\sigma }f(\sigma )P_{\beta }(\sigma )\,}

参数

H (σ) 中两项前的负号是约定俗成的。因为第一项为负号，因此参数 J ij 的正负号决定了该系统的性质，对于每一对 i 和 j :

除此之外该系统为非铁磁性。

在铁磁性的易辛模型中，相邻自旋同方向时能量较低，因此自旋会倾向于同向排列，反之，在反铁磁性的易辛模型中，的相邻自旋反向的能量较低，因此自旋会顷向于反向排列。

H (σ) 中的第二项为负号，表示自旋顷向于和外加磁场同向，因此 h j 的正负也决定自旋顷向排列的方向。对于所有的j ，如果:

简化

一个常见的简化是假设是假设没有外加的磁场作用在易辛模型上，也就是说，对于所有的 j∈Λ， h j = 0 。利用这项简化，其哈密顿量可以写成:

H ( σ σ --> ) = − − --> ∑ ∑ --> J i j σ σ --> i σ σ --> j . {\displaystyle H(\sigma )=-\sum _{ }J_{ij}\sigma _{i}\sigma _{j}.}

此时易辛模型在反转所有自旋之下是对称的:一个外加的力场会破坏这种对称。

另一个常见的简化是假设所有相邻晶格点的交互作用都是相等的，因此可以设 J ij = J 对于所有相邻的 i , j ∈Λ，而其哈密顿量可以写: H ( σ σ --> ) = − − --> J ∑ ∑ --> σ σ --> i σ σ --> j . {\displaystyle H(\sigma )=-J\sum _{ }\sigma _{i}\sigma _{j}.}

特性

一维平移不变铁磁性在无外加场下的易辛模型的所有组态之概率密度投射在二维平面上。

被研究的最透彻易辛模型是在 d维晶格上，平移不变、铁磁性并且无外加场的模型。也就是：Λ = Z , J ij = 1, h = 0。

易辛在他1924年的博士论文中，解决了在 d=1 时的情况，这个一维的模型可以想像成一排的自旋，而每个自旋都只和它左右两边的自旋交互作用。这个一模型不会产生相变，换句话说，对于所有正值的β，任意两自旋的相关系数 i σ j > 都对 | i − j | 呈指数衰减:

⟨ ⟨ --> σ σ --> i σ σ --> j ⟩ ⟩ --> β β --> ≤ ≤ --> C exp ⁡ ⁡ --> ( − − --> c ( β β --> ) | i − − --> j | ) {\displaystyle \langle \sigma _{i}\sigma _{j}\rangle _{\beta }\leq C\exp(-c(\beta )|i-j|)}

其中c(β)是一个只和β有关的函数。

由此可知。这个系统是无序的。根据一维模型的结论，易辛错误的认为任何维度的易辛模型都不会有相变。但事实上，在二维或二维以上的模型中，该系统可以历经从无序相转变成有序相的相变。基本上在β值小(高温)时，系统处在无序相，而β值大(低温)时，系统处在有序相中。换句话说，当系统在有序相时：

⟨ ⟨ --> σ σ --> i σ σ --> j ⟩ ⟩ --> β β --> ≥ ≥ --> c ′ ′ --> ( β β --> ) > 0. {\displaystyle \langle \sigma _{i}\sigma _{j}\rangle _{\beta }\geq c^{\prime }(\beta )>0.\,}

其中c"(β)也是一个只和β有关的函数。

这个性质是首先被鲁道夫·佩尔斯(Rudolf Peierls)在1936年证明的，他的证法后来被称为佩尔斯论述 (Peierls argument)。

在零磁下的二维方晶格易辛模型的解析解后来在1943年被昂萨格解出，他证明了该模型的相关函数（英语： correlation function ）和自由能可由一个无交互作用的格点费米子场(noninteracting lattice fermion)来界定。昂萨格在1949年发表了一个决定了自发磁化现象的公式，但却没有给出推导过程。后来是杨振宁在1951年发表了第一个正式的推导过程，其中里面用到了包括塞格极限定理（英语： Szegő limit theorems ）和弗雷德霍姆行列式（英语： Fredholm determinant ）等数学工具。

数值模拟

蒙特卡洛方法

易辛模型一般来说很难直接进行数值计算，因为他的自旋组态非常之多。考虑一个拥有 L 个晶格点的模型，每个晶格点 σ j 有 ±1 两种可能，因此所有的自旋组态共有 2 种可能，这个数字会随着 L 的增加而进行指数增长。这也是为什么一般在做易辛模型的数值模拟时，都会采用蒙特卡洛方法（Monte Carlo method）

用蒙特卡洛方法来模拟易辛模型所用的哈密顿量为:

此外，可以假设外加磁场 ( h ) 为零以简化模型，因为大部分的问题都只需要用到零磁下的模型。因此，近一步简化的哈密顿量为

可以从这个模型计算出的物理量包括该系统在定温下的比热和磁化强度等等。

Metropolis 算法

Metropolis–Hastings 算法（英语： Metropolis–Hastings algorithm ）是在数值模拟易辛模型时最常用的一种蒙特卡洛方法。这种方法首先要选定一个选择概率 g (μ, ν)，代表系统在状态μ下，在所有可能状态中选到状态ν的概率。另外还需定义出一个接受概率 A (μ, ν) ，也就是说当系统在状态μ下，接受系统跳到态ν的概率。如此的设计是为了让系统达到细致平衡。如果状态ν被接受了，则整个系统就会跳到状态ν，并且选择和决定下一个要跳到的状态。如果状态ν被没被接受，系统则留在状态μ，一样重新选择下一个要跳到的状态。这样的步骤一直重复直到某些条件达成为止，譬如说整个易辛模型完全被磁化，也就是所有的自旋都只到同一个方向。另外在实行这种算法，有一点需要注意的是需要选到适当 g (μ, ν) 以保证整个过程的遍历性(ergodicity)。

在热平衡时，整个系统的能量只会有小幅度的扰动，这点促成了在演算时采用单一自旋反转法进行计算，也就是说每次系统转换其状态时，只改变其中一个自旋的方向。对自旋数很多的一易辛模型来说，系统在不同的状态之间跳跃时，其能量改变的幅度都很小。事实上，对于每个晶格点都和 c 个晶格点相邻的模型来说，每次能量改变的最大幅度为 2 cJ 。此外，采用这种单一自旋反转法可以保证演算过程的遍历性，因为任意一个状态都可以借由逐次的反转相异的自旋，而变成任意其他状态。

在下面的演算过程我们会采用周期性边界条件使的每个晶格点的相邻数 c 都相等。

演算过程

整个用于数值模拟易辛模型的演算过程，可由下列的方法建立。

因为共有 L 个晶格点，而且单一自旋反转法是唯一系统可以从一个状态跳到另一个状态的方法，所对于任何一个状态 μ，有 L 个新状态 ν 它可以跳去。可以假设 μ 跳到这 L 个新状态的概率是相等的，因此 g (μ, ν) = 1/ L 。为了满足细致平衡（英语： detailed balance ），所以下面的等式要成立:

因此我们希望接受概率满足:

如果 H ν > H μ 则 A (ν, μ) > A (μ, ν)，因此 Metropolis 另 A (μ, ν) 和 A (ν, μ) 中较大的为 1。按照这样的定可以得到 A (μ, ν) 的值为:

A ( μ μ --> , ν ν --> ) = { e − − --> β β --> ( H ν ν --> − − --> H μ μ --> ) , if H ν ν --> − − --> H μ μ --> > 0 1 , otherwise . {\displaystyle A(\mu ,\nu )={\begin{cases}e^{-\beta (H_{\nu }-H_{\mu })},&{\text{if }}H_{\nu }-H_{\mu }>0\\1,&{\text{otherwise}}.\end{cases}}}

最后，整个算法最基本的形式为:

用 g (μ, ν) 选出一个自旋，并且计算所有和其自旋相关的能量贡献。(也就是所有和它相邻自旋交互作用的能量贡献。)

反转其自旋，再计算一次所有和它相关的能量贡献。

如果其能量贡献下降，保持自旋反转。

如果贡献能量上升，则令自旋有 e − − --> β β --> ( H ν ν --> − − --> H μ μ --> ) {\displaystyle e^{-\beta (H_{\nu }-H_{\mu })}} 的概率保持反转。

重复步骤一。

整个系统能量的改变量 H ν − H μ 只和反转的自旋和和它相邻的自旋有关，所以只要相邻的自旋数不要太多，计算的速度算是相当快的。而整个系统会逐渐的趋近于某个平衡。

将易辛模型视为马尔可夫链

将易辛模型比拟为马尔可夫链是一件很容易的事情，因为下一刻状态 ν 的转移概率 P β (ν) 只和目前状态 μ有关。事实上，Metropolis 算法就是马尔可夫蒙特卡洛(Markov Chain Monte Carlo)的一个版本。而且因为演算时采用的是单一自旋反转法，每个状态可视为和另外 L 个状态相连，而每一个状态转移则是将其中一个自旋反转。此外，因为不同状态之间 H σ 变化只和相邻自旋的交互作用 J 有关，所以易辛模型也可视为选民模型（英语： Voter Model ）的一种。

一维易辛模型

在一维易辛模型系统中，假设每个带有自旋的原子分布在一维的圆圈中，且原子仅和邻居发生交互作用，交互作用均为 J {\displaystyle J} ，能量可表示为

其中 h 为外加磁场的强度， J 为相邻原子的耦合强度。该系统的自由能为：

相邻自旋的相关函数（英语： Correlation function (statistical mechanics) ）为：

其中，C(β)和 c (β)是关于温度的函数，当温度 T > 0时取正值。当 T → 0时， c (β)趋于零。

证明

易辛自己的论文中已给出了一维易辛模型配分函数精确解。无外加磁场，即 h = 0时，计算是很简单的。

此时能量简化为：

对于相邻两原子自选朝向的四种情况，之间的能量只有两种状态：同向和反向。应用变量替换，

代入配分函数定义式，得：

因此，系统的自由能为

同样的变量替换法，得间隔N-1个自旋间的相关函数：

T ≠ 0时相关性呈指数衰减。而仅仅在绝对零度时 T = 0下，换句话说，β → ∞ 时，系统保持长程相关性。

如果外磁场 h ≠ 0，配分函数的计算需要引入传递矩阵法。周期性边界条件的近似下，配分函数为

系数 V σ σ --> , σ σ --> ′ {\displaystyle V_{\sigma ,\sigma "}} 可视为2x2矩阵的元素。相邻两自旋各有四种可能状态，其玻尔兹曼因子分别为：

或写作传递矩阵的四个元素：

配分函数的展开式恰好和传递矩阵自乘后对角元素之和（矩阵的迹）一致。而矩阵的迹可通过求解传递矩阵的特征值得出，因而：

其中λ 1 是 V 最大的特征值，λ 2 是该矩阵另一个特征值。因而

由于|λ 2 | < λ 1 ，当 N 很大时对配分函数的贡献可忽略。

结论

系统最低能量为− N ，此时所有自旋朝向相同。而其它构型相较基态的能量增量等于一维序列中自旋转向的次数 k 。例如，能量次低的状态相较基态，能量差为2 k 。由于能量与转向数之间是线性关系，则两个相邻自旋方向相反的概率符合玻尔兹曼分布：

借由统计力学的配分函数可以计算再给定温度下 T {\displaystyle T} ( β β --> = 1 / ( k T ) {\displaystyle \beta =1/(kT)} )的每个原子的磁矩期望值为

所以一维易辛模型并没有居里温度、不会发生相变，即没有自发磁化（英语： spontaneous magnetization ）的现象。

二维易辛模型

设二维方晶格易辛模型横纵两方向的交互作用能分别为 J 1 {\displaystyle J_{1}} and J 2 {\displaystyle J_{2}} 。拉斯·昂萨格求得无外磁矩，即 h = 0 {\displaystyle h=0} 时自由能的解析解：

从自由能的偏导数可得到各种热力学函数。

特别地，二维易辛模型有一个相变点，临界温度 T c {\displaystyle T_{c}} 满足以下方程：

若 J 1 = J 2 {\displaystyle J_{1}=J_{2}} ，则 T c = 2 J k ln ⁡ ⁡ --> ( 1 + 2 ) {\displaystyle T_{c}={\frac {2J}{k\ln(1+{\sqrt {2}})}}} ，或 β β --> = 1 k T c ≈ ≈ --> 0.44 J {\displaystyle \beta ={\frac {1}{kT_{c}}}\approx {\frac {0.44}{J}}} .