相变和临界现象

L σ iσ j[cosh N βJ + cosh N −1 βJ sinh βJ (σ1σ 2 + σ 2σ 3 +L+ σ Nσ1)
σ1 σ2 σ N
+ cosh N −2 βJ sinh2 βJ (σ1σ 2σ 2σ 3 + σ 2σ 3σ 3σ 4 +L
+ σ1σ 2σ 3σ 4 + σ1σ 2σ 4σ 5 +L)
这样一来,
Z = ∑ ∑L∑(cosh βJ + σ1σ 2 sinh βJ )(cosh βJ + σ 2σ3 sinh βJ ) σ1 σ2 σ N L(cosh βJ + σ Nσ1 sinh βJ )
展开来有
∑∑ ∑ Z =
L [cosh N βJ + cosh N −1 βJ sinh βJ (σ1σ 2 + σ 2σ 3 +L+ σ Nσ1)
研究一下其中的各项 会发现, 只有两项目对自旋变量求和后不为零, 即
Z = 2N (coshN βJ + sinh N βJ );
lim Z = 2N cosh N βJ (1+ tanh N βJ )
N →∞
我们可以把其中的各项用图形表示出来
= 2N cosh N βJ
σ1σ 2 : 1 2 3 4 5
相互作用 J 起源于原子之间的交叠积分. 考虑两个原子, 位于晶格上邻近的两个
格点上. 设在格点上的单原子波函数为 ψ i 和 ψ j , 用下面的示意图表示.
ψi ψ j
i
j
为简单起见, 考虑原子最外层只有一个电子, 如果这两个原子中的最外层的电子
的自旋构成一个三重态: |↑>i|↑> j ; |↓>i|↓> j ; (|↑>i|↓> j + |↓>i|↑> j ) / 2
容易得到 Z = 2(eJ + e−J )
关联函数定义为 显式地写出来为
∑ σ1σ 2 = Z −1Tr(σ1σ 2e−βH ) = σ1σ 2e−βH {σ i }
∑ ∑ σ1σ 2 = Z −1
σ1σ 2eβJσ1σ 2
σ1 σ2
= Z −1[1×1× eβJ +1× (−1) × e−βJ + (−1) ×1× e−βJ + (−1) × (−1) × eβJ ]
伊辛模型很简单(至少看起来很简单), 一,二维情况有精确解; 对于这个模型的各方 面的性质的研究, 很多时候, 可以给出直观的图象. 现代的 Monte Carlo模拟可以 很容易地研究这样的系统. 它是我们理解相变和临界现象的一个最好的入门模型.
对于学统计物理的人来说, 伊辛模型的地位相当于我们学英语时的字母歌.
∫ 如果观察时间T远小于渡越时间, 则有
T 0
σ1
≠
0
≠
σ
,
这也可以称为各态历经的破缺.
(3)如果系统是二维, 三维的很大的系统 , TC ≠ 0. 平衡态时, 自由能取极小. 而 自由能为 F = U − TS . 能量和熵是两个竞争的因素. 有序的态能量小, 但 − TS 大, 而无序意味着能量大, 但 − TS 小. T < TC 时, 能量起主导作用, 所以系统是 有序的, T > TC , 熵起主导作用, 系统是无序的.
则空间波函数是反对称的, 为 ψ triplet = [ψ i (r1)ψ j (r2 ) −ψ i (r2 )ψ j (r1)]/ 2
两电子自旋的点积大小为 < sˆi ⋅ sˆ j >triplet = h2 / 4
系统的能量为
∫ Etriplet =
d
3r1d
3r2ψ
* triplet
Hˆψ
triplet
相变和临界现象
¾相变和临界现象的历史回顾 ¾自旋模型 ¾一维伊辛模型的解
p
一个多世纪以前，英国学者安德鲁斯便开
始了临界现象的研究。他考察了不同温度
临界点
条件下，二氧化碳的压力如何随体积而改 变，他发现，P-V图上的等温线的平直部
分随温度升高而变短，当达到某一温度
时，平直部分的长度变为零，这就是所谓
| 4 >= (|↑>1|↓>2 − |↓>1|↑>2 ) / 2; s = 0, s z = 0, E = −3J / 4
它的配分函数为 作业: 求关联函数
∑ Z = Tr(e−βHˆ ) = 4 e−βEi = e−βJ / 4 + 3e3βJ / 4 i=1
∑ sˆ1xsˆ2x = Z −1Tr(sˆ1xsˆ2xe−βHˆ ) = Z −1 4 < i | sˆ1xsˆ2x | i > e−βEi i=1
<ij >
<ij >
前面的常数项对热力学 没有贡献.
海森堡模型是个量子模型. 它的配分函数可以写为
∑ Hˆ = J sˆi ⋅ sˆ j
<ij >
Z = Tr(e−βHˆ )
为简单起见, 我们考虑两个格点的海森堡模型
Hˆ = Jsˆ1 ⋅ sˆ2
它的希尔伯特空间有四个基: |↑>1|↑>2; |↑>1|↓>2; |↓>1|↑>2; |↓>1|↓>2
由于 sˆ1 ⋅ sˆ2 = [(sˆ1 + sˆ2 )2 − sˆ12 − sˆ22 ]/ 2 = [sˆ2 − 3 / 2]/ 2 其中, sˆ = sˆ1 + sˆ2
哈密顿量的本征态和本征值为 |1 >=|↑>1|↑>2; s = 1, s z = 1, E = J / 4
| 2 >= (|↑>1|↓>2 + |↓>1|↑>2 ) / 2; s = 1, s z = 0, E = J / 4 | 3 >=|↓>1|↓>2; s = 1, s z = −1, E = J / 4
* j
(r1)
Hˆψ
i
(r1)ψ
j (r2 )
可以把三重态的能量和单重态的能量用它们对应的两电子自旋的点积表示出来, Etriplet = E0 − J < sˆi ⋅ sˆ j >triplet = E0 − Jh2 / 4; Esinglet = E0 − J < sˆi ⋅ sˆ j >singlet = E0 − J (−3h2 / 4);
铁磁相变与气液相变有一个重要区别: 铁磁相变没有两相共存态. 这里涉及一个重要的概念, 对称破缺.
自旋模型
∑ 海森堡模型 Hˆ = −J sˆi ⋅ sˆ j <ij >
其中 sˆi 是第i个格点上的自旋的自旋算符,
Heisenburg W. (1928) Z. Physik 49,619
< ij >表最近邻求和.
+ cosh N −3 βJ sinh3 βJ (σ1σ 2σ 2σ 3σ 3σ 4 +L) +L
+ sinh N βJ (σ1σ 2σ 2σ 3σ 3σ 4 Lσ N −1σ Nσ Nσ1
的临界点，当温度高于此温度时，气体和
液体便不可区分了，正如安德鲁斯所说：
“如果这时有人问你，系统究竟是气体还
是液体，我相信这一问题是不容许有答案
的。”一般地，人们把一类相变的终点称
为临界点，如气液临界点和铁磁相变临界
点。与临界点有关的现象称为临界现象 .
v
也称作连续相变.
从气体到液体可以不连续地变化过去, 比如 沿图中橙色线所示, 发生气液相变时, 比容有 跳跃. 这种情况下, 两相有确且的含义. 气体到液体可以不连续地变化过去, 比如 沿图中蓝色线所示,这种情况下, 两相有确且 的含义.
一维伊辛模型的一个简单解法
考虑一个一维的闭合自旋链, 哈密顿为
N
∑ H = −J σ iσ i+1 = −J (σ1σ 2 + σ 2σ 3 + σ 3σ 4 +L+ σ N −1σ N + σ Nσ1) i=1
它的配分函数为
N
∑ ∑∑ ∑ Z =
βJ ∑σ iσ i+1
e i=1
=
L eβJσ1σ 2 eβJσ 2σ3 LeβJσ Nσ1
系统的能量为
∫ Esin glet =
d
3r1d
3r2ψ
* sin
glet
Hˆψ
sin
glet
简单地可证明 Etriplet = K0 − K1;
Esin glet = K0 + K1
其中
∫ K0 =
d
3r1d
3r2ψ
* i
(r1)ψ
* j
(r2
)Hˆψ
i
(r1)ψ
j
(r2 )
和
∫ K1 =
d 3r1d 3r2ψ i*(r2 )ψ
sˆ
y j
+
J3sˆiz
sˆ
z j
)
中
J1 = J2
= 0 时, 推广的海森堡模型
<ij >
变为了伊辛模型.
若J > 0 , 称为铁磁型的; 若 J < 0 , 称为反铁磁型的.
J > 0 时系统的基态
J < 0 时系统的基态
伊辛模型, (1)可以给出许多磁性系统相变的定性特征: 高温时无序, 没有宏观磁矩; 低温时有序, 有宏观磁矩. (2)可以描述二元金属的结构相变. (3)它和气液相变同属 于一个普适类.
=