单元系与相变
目录
1 粒子数可变的单元系 3
2 热动平衡条件 4
2.1 热动平衡判据 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.2 热动平衡条件 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.3 稳定平衡条件 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3 单元系的复相平衡 8
3.1 一级相变的克拉珀龙方程 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
3.2 蒸气压方程 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
3.3 气液相变与临界点 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.3.1 范德瓦尔斯气体与麦克斯韦等面积法则 . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.3.2 临界点与临界系数 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
4 表面系统的相变 13
4.1 表面张力与表面能 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
4.2 表面系统的热力学函数 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
4.3 液滴的三相平衡 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
5 超导相变 17
6 顺磁-铁磁相变与朗道连续相变理论 19
6.1 序参量与对称性破缺 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1
6.2 临界点的自由能 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
6.3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
6.4 磁化率 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
1 粒子数可变的单元系
单元系是指只有一种粒子的系统. 在粒子数可变系统中, 引入了一个新的独立变量: 粒子数. 零使用
往的结论, 定义对粒子数归一化后的热力学广延量
𝑢 =
𝑈
𝑁
𝑠 =
𝑆
𝑁
, 𝑣 =
𝑉
𝑁
对于单位粒子数的粒子系统, 自然原先的结论是可用的. 有热力学基本方程
d𝑢 = 𝑇d𝑠 𝑃d𝑣
所以对总的系统而言, 内能的微分就是
d𝑈 = d(𝑁𝑢) = 𝑁d𝑢 + 𝑢d𝑁 = 𝑁 (𝑇d𝑠 𝑃d𝑣) + 𝑢d𝑁
再利用 𝑆, 𝑉 的微分
d𝑆 = d(𝑁𝑠) = 𝑁d𝑠 + 𝑠d𝑁, d𝑉 = d(𝑁𝑣) = 𝑁d𝑣 + 𝑣d𝑁
得到
d𝑈 = 𝑇d𝑆 𝑃d𝑉 + (𝑢 𝑇 𝑠 + 𝑃𝑣)d𝑁
所以定义化学势 𝜇
𝜇 = 𝑢 𝑇 𝑠 + 𝑃𝑣
由于它定义式中的是单位粒子数的热力学量, 所以它是一个强度量. 实际上它是单位粒子数 Gibbs
由能
𝜇
=
𝐺
𝑁
那么热力学基本方程就是
d𝑈 = 𝑇d𝑆 𝑃d𝑉 + 𝜇d𝑁
那么由各个自由能的表达式
𝐻 = 𝑈 + 𝑃𝑉
𝐹 = 𝑈 𝑇 𝑆
𝐺 = 𝑈 𝑇 𝑆 + 𝑃𝑉
可以给出它们的微分形式
d𝐻 = 𝑇d𝑆 + 𝑉d𝑃 + 𝜇d𝑁
d𝐹 = 𝑆d𝑇 𝑃d𝑉 + 𝜇d𝑁
d𝐺 = 𝑆d𝑇 + 𝑉d𝑃 + 𝜇d𝑁
所以化学势可以用偏导数表示为
𝜇 =
𝜕𝐺
𝜕𝑁
𝑇,𝑃
=
𝜕𝐹
𝜕𝑁
𝑇,𝑉
=
𝜕𝐻
𝜕𝑁
𝑆, 𝑃
如果利用勒让德变换将热力学基本方程中的 (𝑆, 𝑉, 𝑁) 换成 (𝑇, 𝑉, 𝜇) 则得到巨势
Ψ = 𝑈 𝑇𝑆 𝜇𝑁
其微分为
dΨ = 𝑆d𝑇 𝑃d𝑉 𝑁d𝜇
2 热动平衡条件
2.1 热动平衡判据
孤立系熵不减小, 所以孤立系的平衡态熵达到极大, 得到热动平衡的充要条件
𝑆 = 𝑆
max
孤立系处于平衡态
如果系统有 𝑛 个自由度, 那么熵由 𝑛 个独立变量给出
𝑆 𝑆(𝑥
1
, 𝑥
2
, · · · , 𝑥
𝑛
)
那么得到熵判据: 系统的任何变化都会导致熵减小
𝛿𝑆 = 0
𝛿
2
𝑆 < 0
𝛿𝑈 = 0, 𝛿𝑉 = 0, 𝛿𝑁 = 0
其中
𝛿𝑆 =
𝑛
Õ
𝑖=1
𝜕𝑆
𝜕𝑥
𝑖
𝑥
0
1
,𝑥
0
2
,··· , 𝑥
0
𝑛
𝛿𝑥
𝑖
, 𝛿
2
𝑆 =
𝑛
Õ
𝑖=1
𝑛
Õ
𝑗=1
𝜕
2
𝑆
𝜕𝑥
𝑖
𝜕𝑥
𝑗
𝑥
0
1
,𝑥
0
2
,··· , 𝑥
0
𝑛
𝛿𝑥
𝑖
𝛿𝑥
𝑗
第一条表明熵达到极值, 第二条表明熵达到极大值. 方程可能由多个, 其中熵最大的对应, 其他
的解对应于亚稳态
如果 𝛿
2
𝑆 = 0, 么稳态要求 𝛿
3
𝑆 = 0, 为三次项是奇的, 若其非零则正向和负向的变化必然会导致一侧
熵减小. 此时的稳态要求四次项为负
𝛿
4
𝑆 < 0
也存在其他的判据
自由能判据: 等温等容系统处于平衡态时, 自由能达到极小
吉布斯自由能判据: 等温等压系统处于平衡态时, 吉布斯自由能达到极小
内能判据: 等容等熵等粒子数系统处于平衡态时, 内能达到极小
2.2 热动平衡条件
热动平衡条件是指
热平衡: 物体内部不发生热量交换
力平衡: 物体内部不发生宏观位移
相平衡: 物体内部不发生物质转移 (即相变)
化学平衡
:
不发生化学反应
设系统由两个均匀部分组成
,
对整个系统来说
,
各个广延量是两个子系统广延量之和
𝑆 = 𝑆
1
+ 𝑆
2
, 𝑈 = 𝑈
1
+ 𝑈
2
, 𝑉 = 𝑉
1
+ 𝑉
2
, 𝑁 = 𝑁
1
+ 𝑁
2
取各个广延量 𝑈, 𝑉 , 𝑁 为自变量, 那么考察某个子系统的熵 𝑆
𝛼
的变分 (𝛼 = 1, 2), 对于每个子系统, 有热
力学基本方程
d𝑈
𝛼
= 𝑇
𝛼
d𝑆
𝛼
𝑃
𝛼
d𝑉
𝛼
+ 𝜇
𝛼
d𝑁
𝛼
d𝑆
𝛼
=
1
𝑇
𝛼
d𝑈
𝛼
+
𝑃
𝛼
𝑇
𝛼
d𝑉
𝛼
𝜇
𝛼
𝑇
𝛼
d𝑁
𝛼
对于孤立系统而言, 总内能, 总体积, 总粒子数都是常数, 即有约束条件
𝛿𝑈 = 𝛿𝑉 = 𝛿𝑁 = 0
𝛿𝑈
1
= 𝛿𝑈
2
, 𝛿𝑉
1
= 𝛿𝑉
2
, 𝛿𝑁
1
= 𝛿𝑁
2
由于 𝑆 = 𝑆
1
+ 𝑆
2
, 所以总熵的变分是两个子系统变分的和
𝛿𝑆 = 𝛿𝑆
1
+ 𝛿𝑆
2
代入即得
𝛿𝑆 =
1
𝑇
1
1
𝑇
2
𝛿𝑈
1
+
𝑃
1
𝑇
1
𝑃
2
𝑇
2
𝛿𝑉
1
𝜇
1
𝑇
1
𝜇
2
𝑇
2
𝛿𝑁
1
式中的 𝛿𝑈
1
, 𝛿𝑉
1
, 𝛿𝑁
1
是独立的, 为了满足熵平衡条件
𝛿𝑆 = 0
就必须有
𝑇
1
= 𝑇
2
, 𝑃
1
= 𝑃
2
, 𝜇
1
= 𝜇
2
其中的 𝑇
1
= 𝑇
2
是热平衡条件. 如果令 𝛿𝑉
1
= 𝛿𝑁
1
= 0, 那么
d𝑆 =
1
𝑇
1
1
𝑇
2
d𝑈
1
如果此时 𝑇
1
> 𝑇
2
, 为了使得 d𝑆 > 0, 就必须有 d𝑈
1
< 0, 即热量从高温物体传递到低温物体
𝑃
1
= 𝑃
2
是力平衡条件. 如果令 𝛿𝑈
1
= 𝛿𝑁
1
= 0, 𝑇
1
= 𝑇2, 那么
d𝑆 =
𝑃
1
𝑇
1
𝑃
2
𝑇
2
d𝑉
1
如果此时 𝑃
1
> 𝑃
2
, 为了使得 d𝑆 > 0, 就必须有 d𝑉
1
> 0, 即压强高的子系统膨胀
𝜇
1
= 𝜇
2
是化学平衡条件. 如果令 𝛿𝑈
1
= 𝛿𝑉
1
= 0, 𝑇
1
= 𝑇
2
, 那么
d𝑆 =
𝜇
1
𝑇
1
𝜇
2
𝑇
2
d𝑁
1
如果此时 𝜇
1
> 𝜇
2
, 为了使得 d𝑆 > 0, 就必须有 d𝑁
1
> 0, 即化学势高的子系统粒子减少
使用自由能判据也能得到相同的平衡条件
2.3 稳定平衡条件
平衡除了要求熵达到极值外, 还要求达到极大值, 即要求二阶导数为负
𝛿
2
𝑆 < 0
不够由于热力学基本方程是关于 𝑈 , 因而使用内能条件
𝛿
2
𝑈 > 0
对两个子系统分布使用热力学基本方程
𝛿𝑈
𝛼
= 𝑇
𝛼
𝛿𝑆
𝛼
𝑃
𝛼
𝛿𝑉
𝛼
+ 𝜇
𝛼
𝛿𝑁
𝛼
将各个广延量对粒子数归一化
𝑆
𝛼
= 𝑁
𝛼
𝑠
𝛼
, 𝑉
𝛼
= 𝑁
𝛼
𝑣
𝛼
它们的变分为
𝛿𝑆
𝛼
= 𝑠
𝛼
𝛿𝑁
𝛼
+ 𝑁
𝛼
𝛿𝑠
𝛼
, 𝛿𝑉
𝛼
= 𝑣
𝛼
𝛿𝑁
𝛼
+ 𝑁
𝛼
𝛿𝑣
𝛼
代入即得
𝛿
2
𝑈
𝛼
= 𝑁
𝛼
(
𝛿𝑇
𝛼
𝛿𝑠
𝛼
𝛿𝑃
𝛼
𝛿𝑣
𝛼
)
+
(
𝑠
𝛼
𝛿𝑇
𝛼
𝑣
𝛼
𝛿𝑃
𝛼
+ 𝛿𝜇
𝛼
)
𝛿𝑁
𝛼
利用化学势的微分
d𝜇 = 𝑠d𝑇 + 𝑣d𝑃
𝛿𝜇
𝛼
= 𝑠
𝛼
𝛿𝑇
𝛼
+ 𝑣
𝛼
𝛿𝑃
𝛼
代入即得
𝛿
2
𝑈
𝛼
= 𝑁
𝛼
(
𝛿𝑇
𝛼
𝛿𝑠
𝛼
𝛿𝑃
𝛼
𝛿𝑣
𝛼
)
由于 𝑈 = 𝑈
1
+ 𝑈
2
, 所以
𝛿
2
𝑈 = 𝛿
2
𝑈
1
+ 𝛿
2
𝑈
2
平衡条件即
Õ
𝛼
𝑁
𝛼
(
𝛿𝑇
𝛼
𝛿𝑠
𝛼
𝛿𝑃
𝛼
𝛿𝑣
𝛼
)
> 0
该式不包 𝛿𝑁
𝛼
, 其余的变分也都是强度量, 这说明平衡条件与粒子数的变化无关, 所以对于任意的粒
, 只要强度量满足平衡条件, 𝑁
𝛼
取任何值都是平衡态. 这要求两相分别满足
𝛿𝑇𝛿𝑠 𝛿𝑃𝛿𝑣 > 0
至此粒子数这一变量被消去了, 系统回到了熟悉的二自由度情形. (𝑇, 𝑣) 为自变量, 那么有变分
𝛿𝑠 =
𝜕𝑠
𝜕𝑇
𝑣
𝛿𝑇 +
𝜕𝑠
𝜕𝑣
𝑇
𝛿𝑣, 𝛿𝑃 =
𝜕𝑃
𝜕𝑇
𝑣
𝛿𝑇 +
𝜕𝑃
𝜕𝑣
𝑇
𝛿𝑣
利用 Maxwell 关系
𝜕𝑠
𝜕𝑣
𝑇
=
𝜕𝑃
𝜕𝑇
𝑣
代入即得
𝜕𝑠
𝜕𝑇
𝑣
(𝛿𝑇)
2
𝜕𝑃
𝜕𝑣
𝑇
(𝛿𝑣)
2
> 0
利用偏导数关系
𝜕𝑠
𝜕𝑇
𝑣
=
𝜕𝑢
𝜕𝑇
𝑣
·
𝜕𝑠
𝜕𝑢
𝑣
=
𝐶
𝑣
𝑇
那么平衡条件就是
𝑐
𝑣
𝑇
(𝛿𝑇)
2
𝜕𝑃
𝜕𝑣
𝑇
(𝛿𝑣)
2
> 0
鉴于后面两个变分都是平方项, 本身就是正的, 所以平衡条件要求它们的系数是正的
𝑐
𝑣
𝑐
𝑣
> 0
𝜕𝑃
𝜕𝑣
𝑇
< 0
即要求恒压热容和等温压缩系数为正
𝑐
𝑣
> 0, 𝜅
𝑇
=
1
𝑣
𝜕𝑣
𝜕𝑃
𝑇
> 0
3 单元系的复相平衡
如果单元系有两种相 𝛼, 𝛽, 那么两相平衡时应当满足平衡条件
𝑇
𝛼
= 𝑇
𝛽
𝑃
𝛼
= 𝑃
𝛽
𝜇
𝛼
= 𝜇
𝛽
𝑇 两相的共同温度, 𝑃 两相的共同压强. 由于相平衡与粒子数无关, 么平衡的系统自由度 2,
选取 (𝑇 , 𝑃) 为自变量, 那么相平衡条件就是
𝜇
𝛼
(𝑇, 𝑃) = 𝜇
𝛽
(𝑇, 𝑃)
𝑇 𝑃 平面上有两个自由变量和一个约束, 将得到一条曲线, 相平衡曲线. 若是有三相平衡,
𝜇
𝛼
(𝑇, 𝑃) = 𝜇
𝛽
(𝑇, 𝑃) = 𝜇
𝛾
(𝑇, 𝑃)
𝑇 𝑃 平面上有三个自由变量和两个约束, 将得到一个点, 称为三相点
3.1 一级相变的克拉珀龙方程
在相平衡曲线上的一个点应有化学势相等
𝜇
𝛼
(𝑇, 𝑃) = 𝜇
𝛽
(𝑇, 𝑃)
假定沿着相平衡曲线移动了一个微元, 那么
𝜇
𝛼
(𝑇 + d𝑇, 𝑃 + d𝑃) = 𝜇
𝛽
(𝑇 + d𝑇, 𝑃 + d𝑃)
展开即得
d𝜇
𝛼
= d𝜇
𝛽
利用化学势的微分
d𝜇 = 𝑠d𝑇 + 𝑣d𝑃
由于平衡时两相具有相同的温度和压强, 所以
𝑠
𝛼
d𝑇 + 𝑣
𝛼
d𝑃 = 𝑠
𝛽
d𝑇 + 𝑣
𝛽
d𝑃
d𝑃
d𝑇
=
𝑠
𝛼
𝑠
𝛽
𝑣
𝛼
𝑣
𝛽
沿着平衡曲线只有一个自由度, 所以不妨设
𝑃 𝑃(𝑇)
上式即给出了 𝑃 𝑇 曲线的斜率. 定义相变潜热 𝜆
𝛼𝛽
为保持温度和压强不变时单位粒子数由相 𝛽 变为相
𝛼 时吸收的热量
𝜆
𝛼𝛽
= d¯𝑄 : 𝛽 𝛼
这应当是一个准静态过程, 所以
𝜆
𝛼𝛽
= 𝑇Δ𝑠 = 𝑇 (𝑠
𝛼
𝑠
𝛽
)
因而上式可以写为
d𝑃
d𝑇
=
𝜆
𝛼𝛽
𝑇 (𝑣
𝛼
𝑣
𝛽
)
克拉珀龙方程. 为了使得分母有意义, 需要相变时体积发生改变, 即发生一级相变. 在平衡时两相具
相同的化学势
𝜇
𝛼
(𝑇, 𝑃) = 𝜇
𝛽
(𝑇, 𝑃)
展开来就是
𝑢
𝛼
𝑇 𝑠
𝛼
+ 𝑃𝑣
𝛼
= 𝑢
𝛽
𝑇 𝑠
𝛽
+ 𝑃𝑣
𝛽
若定义对粒子数归一化的焓
=
𝐻
𝑁
= 𝑢 + 𝑃𝑣
那么相变潜热实际上就是两相的焓差
𝜆
𝛼𝛽
=
𝛼
𝛽
3.2 蒸气压方程
与凝聚相 (液相和固相) 达到平衡的蒸气称为饱和蒸气. 饱和蒸气的压强与温度的关系称为蒸气压方程
𝜆 为由凝聚相 (液相或固相) 变为气相的相变潜热, 那么有克拉珀龙方程
d𝑃
d𝑇
=
𝜆
𝑇 (𝑣 𝑣
)
其中 𝑣 为气相的摩尔体积, 𝑣
为凝聚相的摩尔体积. 用焓表示相变潜热
𝜆 =
沿着相平衡曲线只有一个自由度, 设为 𝑇, 那么希望求 𝜆( 𝑇 ). 考察它的微分
d𝜆
d𝑇
=
d
d𝑇
d
d𝑇
取焓函数的自变量为 (𝑃, 𝑇), 再施加约束 𝑃 = 𝑃(𝑇), 那么 (𝑇, 𝑃(𝑇)), 所以
d
d𝑇
=
𝜕
𝜕𝑇
𝑃
+
𝜕
𝜕𝑃
𝑇
d𝑃
d𝑇
前一项就是恒压热容 𝑐
𝑝
. 为了考察后一项, 利用焓形式的热力学基本方程
d = 𝑇d 𝑠 + 𝑣d𝑃
将自变量换为 (𝑇 , 𝑃), 那么 𝑠 = 𝑠(𝑇 , 𝑃), 𝑠 的全微分为
d𝑠 =
𝜕𝑠
𝜕𝑇
𝑃
d𝑇 +
𝜕𝑠
𝜕𝑃
𝑇
d𝑃
Maxwell 关系
𝜕𝑠
𝜕𝑃
𝑇
=
𝜕𝑣
𝜕𝑇
𝑃
那么
d𝑠 =
𝜕𝑠
𝜕𝑇
𝑃
d𝑇 𝑣
𝜕𝑣
𝜕𝑇
𝑃
d𝑃
代入即得
d = 𝑇
𝜕𝑠
𝜕𝑇
𝑃
d𝑇 +
𝑣 𝑇
𝜕𝑣
𝜕𝑇
𝑃
d𝑃
所以
𝜕
𝜕𝑃
𝑇
= 𝑣 𝑇
𝜕𝑣
𝜕𝑇
𝑃
所以得到焓对温度的全微分
d
d𝑇
= 𝑐
𝑝
+ 𝑣
d𝑃
d𝑇
𝑇
𝜕𝑣
𝜕𝑇
𝑃
d𝑃
d𝑇
所以
d𝜆
d𝑇
= 𝑐
𝑝
𝑐
𝑝
+ 𝑣
d𝑃
d𝑇
𝑣
d𝑃
d𝑇
𝑇
𝜕𝑣
𝜕𝑇
𝑃
d𝑃
d𝑇
+ 𝑇
𝜕𝑣
𝜕𝑇
𝑃
d𝑃
d𝑇
代入克拉珀龙方程
d𝑃
d𝑇
=
𝜆
𝑇 (𝑣 𝑣
)
得到
d𝜆
d𝑇
= 𝑐
𝑝
𝑐
𝑝
+
𝜆
𝑇
𝜕𝑣
𝜕𝑇
𝜕𝑣
𝜕𝑇
𝜆
𝑣 𝑣
由于凝聚相的摩尔体积小于气相的摩尔体积, 𝑣
近似为 0. 固体的体积变化显然也远小于气体, 所以 𝑣
偏导数也可以忽略; 再将蒸气看作理想气体,
𝑣 =
𝑅𝑇
𝑃
克拉珀龙方程就变为
d𝑃
d𝑇
=
𝜆
𝑅𝑇
2
𝑃
替换相变潜热中的导数得到
d𝜆
d𝑇
= 𝑐
𝑝
𝑐
𝑝
再对其积分
𝜆 =
𝑇
𝑇
0
(𝑐
𝑝
𝑐
𝑝
)d𝑇 + 𝜆
0
其中 𝑇
0
是参考温度, 𝜆
0
是在 𝑇
0
处的相变潜热. 那么克拉珀龙方程就变为
d𝑃
d𝑇
=
𝑃
𝑅𝑇
2
𝑇
𝑇
0
(𝑐
𝑝
𝑐
𝑝
)d𝑇 + 𝜆
0
如果蒸气压变化不大, 对凝聚态来说近似为恒压, 从而 𝑐
𝑝
𝑃 无关, 那么上面的方程就有解
ln 𝑃 =
𝜆
0
𝑅𝑇
+
𝑇
𝑇
0
d𝑇
𝑅𝑇
2
𝑇
𝑇
0
(𝑐
𝑝
𝑐
𝑝
)d𝑇 + 𝐴
0
其中 𝐴
0
是积分常数. 果进一步认为温度也变化不大, 那么比热进一步与温度无关成为常数, 上式就进
一步简化为
ln 𝑃 = 𝐴
𝐵
𝑇
+ 𝐶 ln 𝑇
3.3 气液相变与临界点
3.3.1 范德瓦尔斯气体与麦克斯韦等面积法则
利用范德瓦尔斯气体方程
𝑃
+
𝑎
𝑣
2
(𝑣 𝑏) = 𝑅𝑇
如果画出它的等温线
它看起来与实际的等温线非常接近, 不过在温度低 𝑇
𝐶
时出现了一个 𝑃 对应三个 𝑉 的情况. 并且稳定
平衡条件要求
𝜕𝑃
𝜕𝑉
𝑇
< 0
显然在 𝑉
𝑚1
𝑉
𝑚2
区间内违背了这一条件. 为了找到符合物理的等温线, 考察化学势的微分
d𝜇 = 𝑠d𝑇 + 𝑣d𝑃
由于等温线上 d𝑇 = 0, 所以
d𝜇 = 𝑣d𝑃
也就是说沿着等温线的化学势变化可以由积分得到
𝜇 𝜇
0
=
𝑃
𝑃
0
𝑣d𝑃
这实际上是等温线和 𝑃 轴围成的面积. 画出 𝜇 𝑃 图即
鉴于化学势是对粒子数归一化的吉布斯自由能, 两相平衡条件要求吉布斯自由能最小, 因而合法的变化曲
线应当是
𝑄 𝐾 𝐴 𝐵 𝑀 𝑅
𝐴 𝐵 上面那段认为是不符合物理规律的, 舍去. 𝐴 𝐵 的过程即发生了相变, 𝐴 点是完全的液体,
𝐵 点是完全的气体. 在相变过程中压强和化学势保持不变, 这意味着需要在等温线中找到一根平行于 𝑣
的直线段 𝐴𝐵 代替原来的 𝐴𝐵 , 这条直线段的选取应当使得 𝐴, 𝐵 两点具有相同的化学势,
𝜇
𝐵
𝜇
𝐴
=
𝐵
𝐴
𝑣d𝑃 = 0
这也就是说两块面积 𝐴𝐽𝐷 𝐷 𝑁 𝐵 相等, 称为 Maxwell 等面积法则
3.3.2 临界点与临界系数
虽然范德瓦尔斯等温线上 𝐽 𝑁 的过程不满足稳定平衡条件, 𝐴 𝐽 𝑁 𝐵 的过程是满足的,
分别对应过热液体和过饱和蒸气. 当温度升高到一定程度时极小值点 𝐽 和极大值点 𝑁 重合为点 𝐶, 这个
点称为临界点, 作为拐点它满足方程
𝜕𝑃
𝜕𝑣
𝑇
= 0,
𝜕
2
𝑃
𝜕𝑣
2
𝑇
= 0
代入范德瓦尔斯方程即得到临界点
𝑇
𝐶
=
8𝑎
27
𝑏𝑅
, 𝑃
𝐶
=
𝑎
27
𝑏
2
, 𝑣
𝐶
= 3𝑏
由此得到三者的关系
𝑅𝑇
𝐶
𝑃
𝐶
𝑣
𝐶
=
8
3
称其比值为临界系数. 如果将各个变量对临界点做归一化
𝑣
=
𝑣
𝑣
𝐶
, 𝑃
=
𝑃
𝑃
𝐶
, 𝑇
=
𝑇
𝑇
𝐶
称它们为对比温度, 对比压强, 对比体积. 那么范德瓦尔斯方程就变为范氏对比方程
𝑃
+
3
𝑣
2
(𝑣
1
3
) =
8
3
𝑇
它不含物质性质有关的常量, 具有普适性
4 表面系统的相变
4.1 表面张力与表面能
对于表面张力有表面张力系数 𝜎, 它定义为
d 𝑓 = 𝜎d𝑙
其中 d 𝑓 d𝑙 长度上的表面张力, 力的方向垂直于线元朝向一侧的液面. 考察一个球形的液面, 对于液面
上的一个小圆形面积元 𝐴
𝐴 = 𝜋𝑟
2
由于 𝑟 很小, 圆心与球心连线和球心到边缘连线的夹角就是
𝜃 =
𝑟
𝑅
其中 𝑅 球的半. 于边缘上的一小段弧 d𝑙, 它受到其他液面沿着切线方向的表面张力 d 𝑓 = 𝜎d𝑙,
其中垂直于切线方向的分量为
d 𝑓
= 𝜎d𝑙 sin 𝜃 =
𝜎𝑟
𝑅
d𝑙
考察压强只关心垂直分量, 那么总的垂直分量大小由积分给出
𝑓
=
2 𝜋𝑟
0
𝜎𝑟
𝑅
d𝑙 =
2𝜋𝜎𝑟
2
𝑅
所以小面元上的内部压强为
𝑃 =
𝑓
𝐴
=
2𝜎
𝑅
这是单层液面的压强, 由于肥皂泡是双层液面, 所以内部压强还要乘以 2
4.2 表面系统的热力学函数
当表面由于外力在准静态过程中扩张了 d 𝐴 , 表面张力做负功
d¯𝑊 = 𝜎d𝐴
热力学第一定律
d𝑈 = d¯𝑄 + d¯𝑊
由于可逆热机的效率只与温度有关, 那么熵的定义依然成立. 也就是说
d𝑈 = 𝑇d𝑆 + 𝜎d𝐴
做勒让德变换得到 Helmholtz 自由能 𝐹 = 𝑈 𝑇 𝑆, 那么
d𝐹 = 𝑆d𝑇 + 𝜎d 𝐴
所以
𝑆 =
𝜕𝐹
𝜕𝑇
𝐴
, 𝜎 =
𝜕𝐹
𝜕 𝐴
𝑇
注意到第二式, 由于 𝜎 然与面积无关, 所以积分得到 𝐹 = 𝜎𝐴 + 𝐶. 由于 𝐴 = 0 系统不存, 此时 𝐹
自然为 0, 所以 𝐶 = 0,
𝐹 = 𝜎 𝐴
代入第一式即得
𝑆 = 𝐴
d𝜎
d𝑇
进一步由 𝑈 = 𝐹 + 𝑇 𝑆 得到系统的内能
𝑈 = 𝐴
𝜎 𝑇
d𝜎
d𝑇
4.3 液滴的三相平衡
为了考察一个液滴, 设液体相为 𝛼, 气体相为 𝛽, 表面相为 𝛾, 那么有三相的热力学基本方程
d𝐹
𝛼
= 𝑆
𝛼
d𝑇 𝑃
𝛼
d𝑉
𝛼
+ 𝜇
𝛼
d𝑁
𝛼
d𝐹
𝛽
= 𝑆
𝛽
d𝑇 𝑃
𝛽
d𝑉
𝛽
+ 𝜇
𝛽
d𝑁
𝛽
d𝐹
𝛾
= 𝑆
𝛾
d𝑇 + 𝜎d𝐴
视表面为一个几何相, 自然没有粒子数的概念. 在一定温度和总体积下, 上式写为
d𝐹
𝛼
= 𝑃
𝛼
d𝑉
𝛼
+ 𝜇
𝛼
d𝑁
𝛼
d𝐹
𝛽
= 𝑃
𝛽
d𝑉
𝛽
+ 𝜇
𝛽
d𝑁
𝛽
d𝐹
𝛾
= 𝜎d𝐴
在此时对系统微扰, 扰动应满足总体积守恒与粒子数守恒
𝛿𝑛
𝛼
+ 𝛿𝑛
𝛽
= 0, 𝛿𝑉
𝛼
+ 𝛿𝑉
𝛽
= 0
三相自由能的变化为
𝛿𝐹
𝛼
= 𝑃
𝛼
𝛿𝑉
𝛼
+ 𝜇
𝛼
𝛿𝑁
𝛼
𝛿𝐹
𝛽
= 𝑃
𝛽
𝛿𝑉
𝛽
+ 𝜇
𝛽
𝛿𝑁
𝛽
𝛿𝐹
𝛾
= 𝜎𝛿𝐴
结合约束得到总的自由能变化
𝛿𝐹 =
𝑃
𝛼
𝑃
𝛽
𝛿𝑉
𝛼
+
𝜇
𝛼
𝜇
𝛽
𝛿𝑁
𝛼
+ 𝜎𝛿 𝐴
假定液滴是半径为 𝑟 的球体
𝑉
𝛼
=
4𝜋
3
𝑟
3
, 𝐴
𝛼
= 4𝜋𝑟
2
𝛿𝑉
𝛼
= 4𝜋𝑟
2
d𝑟, 𝛿𝐴
𝛼
= 8𝜋𝑟d𝑟
所以自由能变化为
𝛿𝐹 =
𝑃
𝛼
𝑃
𝛽
2𝜎
𝑟
𝛿𝑉
𝛼
+
𝜇
𝛼
𝜇
𝛽
𝛿𝑁
𝛼
鉴于有液面压强公式
𝑃
𝑠𝑢𝑟 𝑓
=
2𝜎
𝑟
由表面力学平衡得到
𝑃
𝛼
= 𝑃
𝛽
+
2𝜎
𝑟
注意此时两相的压强并不相等, 那么两相平衡条件 𝛿𝐹 = 0 就变为
𝜇
𝛼
𝑃
+
2𝜎
𝑟
, 𝑇
= 𝜇
𝛽
(
𝑃
, 𝑇
)
两相具有相同的温度是热平衡的要求. 希望对平面界面的饱和蒸汽压做修正,
𝜇
𝛼
(𝑃, 𝑇) = 𝜇
𝛽
(𝑃, 𝑇)
压强变化不大时认为液体性质不变, 将化学势按压强展开至一阶项
𝜇
𝛼
𝑃
+
2𝜎
𝑟
, 𝑇
= 𝜇
𝛼
(𝑃
, 𝑇 ) +
𝑃
𝑃 +
2𝜎
𝑟
𝜕𝜇
𝛼
𝜕𝑃
𝑇
鉴于
𝜕𝜇
𝛼
𝜕𝑃
𝑇
= 𝑣
𝛼
代入即得
𝜇
𝛼
𝑃
+
2𝜎
𝑟
, 𝑇
= 𝜇
𝛼
(𝑃, 𝑇) +
𝑃
𝑃 +
2𝜎
𝑟
𝑣
𝛼
𝜇
𝛽
(
𝑃
, 𝑇
)
= 𝜇
𝛽
(𝑃, 𝑇) +
𝑃
𝑃 +
2𝜎
𝑟
𝑣
𝛽
视蒸气为理想气体, 鉴于理想气体有化学势
𝜇
𝛽
= 𝑅𝑇
(
𝜑 + ln 𝑃
)
其中 𝜑 𝜑(𝑇) 是与温度有关的函数, 那么
𝜇
𝛽
(
𝑃
, 𝑇
)
= 𝜇
𝛽
(
𝑃, 𝑇
)
+ 𝑅𝑇 ln
𝑃
𝑃
所以
𝑃
𝑃 +
2𝜎
𝑟
𝑣
𝛽
= 𝑅𝑇 ln
𝑃
𝑃
认为 𝑃
𝑃 差不多, 那么认为 𝑃
𝑃 可以忽略, 上式简化为
2𝜎
𝑟
𝑣
𝛽
= 𝑅𝑇 ln
𝑃
𝑃
也就是
ln
𝑃
𝑃
=
2𝜎𝑣
𝛽
𝑟 𝑅𝑇
从而得到平衡液滴半径 𝑟
𝑐
, 中肯半径. 鉴于
𝜇
𝛼
𝑃
+
2𝜎
𝑟
, 𝑇
= 𝜇
𝛼
(𝑃
, 𝑇 ) +
𝑃
𝑃 +
2𝜎
𝑟
𝑣
𝛼
如果在 𝑟
𝑐
的基础上稍微增大液滴半径, 则使得 𝜇
𝛼
减小, 那么气相的粒子将转移到液相, 使得液滴半径继
续增大; 若是稍微减小液滴半径, 则使得 𝜇
𝛼
增大, 那么气相的粒子将转移到液相, 使得液滴半径继续减小
逐渐消失
所以如果蒸气非常干净, 则蒸气压可以超过饱和蒸汽压, 成为过饱和蒸气
气泡, 只需, ,
𝑟 𝑟, 代入上式即得
ln
𝑃
𝑃
=
2𝜎𝑣
𝛽
𝑟 𝑅𝑇
> 0
即要求气泡中蒸气压小于于时的饱和蒸气, 同时气泡的学平衡又要求泡的气体压强于液
压强. 鉴于沸点的定义
沸点是液体的饱和蒸气压等于外界压强时的温度
在液体不多时, 液体压强等于外界压强, 此时饱和蒸气压等于液体压强. 那么
温度小于沸点: 饱和蒸气压小于液体压强
温度等于沸点: 饱和蒸气压等于液体压强
如果液体中溶有其他气体, 由它们提供缺失的压强, 则气泡会在低于沸点的温度出现. 则这部分压强要
由蒸气提供, 这时的温度就要大于沸点, 出现过热液体
5 超导相变
在电动力学中, 当磁介质的磁感应强度增加 dB , 外界需要对介质做功
d¯𝑊 =
HdBd𝑉
如果磁场是均匀的, 那么上式写为
d¯𝑊 = 𝑉H · dB
鉴于 B = 𝜇
0
(H + M ), 那么
d¯𝑊 = 𝑉Hd
1
2
𝜇
0
𝐻
2
+ 𝜇
0
𝑉H · dM
当金属变为超导态后, 其内部磁通完全排出, 在超导体内磁感应强度为零
B = 𝜇
0
(H + M ) = 0
所以此时有
M = H
实验表明, 磁场超过界磁场 𝐻
𝑐
, 超导体将会回到正常. 𝐻
𝑐
是温度的函数, . 超导相变中体积变
化很小, 所以忽略膨胀功, 并取单位体积. 此时的热力学基本方程为
d𝑈 = 𝑇d𝑠 + 𝜇
0
H · dB
在各向同性的超导体中, 磁场与磁感应强度是平行的, 所以
d𝑈 = 𝑇d𝑠 + 𝜇
0
𝐻d𝐵
通过勒让德变换给出吉布斯自由能
𝐺 𝑈 𝑇 𝑠 𝐻𝐵, d𝐺 = 𝑠d𝑇 𝐵d𝐻
那么
𝜕𝐺
𝜕𝐻
𝑇
= 𝐵
积分即得
𝐺 (𝑇, 𝐻) = 𝐺(𝑇, 0)
𝐻
0
𝐵(𝑇, 𝐻
)d𝐻
正常相的金属磁化很弱, 认为 𝑀 = 0, 所以
𝐵 = 𝜇
0
𝐻
所以代入即得正常相的单位体积吉布斯自由能
𝐺
𝑁
(𝑇, 𝐻) 𝐺
𝑁
(𝑇, 0) =
𝜇
0
2
𝐻
2
当金属变为超导态后, 其内部磁通完全排出, 在超导体内磁感应强度为零
B = 0
所以超导相的单位体积吉布斯自由能为
𝐺
𝑆
(𝑇, 𝐻) 𝐺
𝑆
(𝑇, 0) = 0
因而
𝐺
𝑆
(𝑇, 𝐻) 𝐺
𝑁
(𝑇, 𝐻) = 𝐺
𝑆
(𝑇, 0) 𝐺
𝑁
(𝑇, 0) +
𝜇
0
2
𝐻
2
(𝑇, 𝐻
𝑐
) 两相平, 平衡条件要求化学势相等. 鉴于化学势是对粒子数归一化的吉布斯自由, 由于
相变体积变化可以忽略, 所以等价于两相的单位体积吉布斯自由能相等
𝐺
𝑆
(𝑇, 𝐻
𝑐
) = 𝐺
𝑁
(𝑇, 𝐻
𝑐
)
代入上式即得
𝐺
𝑆
(𝑇, 𝐻) 𝐺
𝑁
(𝑇, 𝐻) =
𝜇
0
2
𝐻
2
𝐻
2
𝑐
需要注意 𝐻
𝑐
是温度依赖的. 𝐻 > 𝐻
𝑐
, 正常相的吉布斯自由能较小, 较稳定; 反之则超导相稳定.
𝐻 = 𝐻
𝑐
时两相共存. 在平衡曲线上应有两相吉布斯自由能相等, 沿着平衡曲线则有
d𝐺
𝑆
= d𝐺
𝑁
鉴于有吉布斯自由能的全微分
d𝐺 = 𝑠d𝑇 𝐵d𝐻
所以
d𝐻
𝑐
d𝑇
=
𝑠
𝑁
𝑠
𝑆
𝐵
𝑁
𝐵
𝑆
利用 𝐻 = 𝐻
𝑐
𝐵
𝑁
= 𝜇
0
𝐻
𝑐
, 𝐵
𝑆
= 0
代入上式即得
𝑠
𝑁
𝑠
𝑆
= 𝜇
0
𝐻
𝑐
d𝐻
𝑐
d𝑇
实验上测得
𝐻
𝑐
(𝑇) = 𝐻
𝑐
(0)
"
1
𝑇
𝑇
𝑐
2
#
, 𝑇 𝑇
𝑐
其中 𝑇
𝑐
是没有外磁场时超导相变的温度. 这表明
d𝐻
𝑐
d𝑇
< 0
所以超导相的熵比正常相的熵小, 从正常相转变到超导相将放热, 相变潜热为
𝐿
𝑆 𝑁
= 𝑇 (𝑠
𝑆
𝑠
𝑁
) = 𝜇
0
𝐻
𝑐
d𝐻
𝑐
d𝑇
𝐻
𝑐
0 时相变潜热不为零, 是一级相变; 𝑇 = 𝑇
𝑐
𝐻
𝑐
= 0, 相变潜热为零, 是二级相变
6 顺磁-铁磁相变与朗道连续相变理论
6.1 序参量与对称性破缺
微观上铁磁体是由大量自旋磁矩组成的, 相邻自旋方向相同时具有最小的能量, 因而自旋倾向于排列在相
同的方向上, 使得无外场时具有非零的磁化强度, 成为自发磁化
当温度升高, 热运动倾向于使得自旋方向混. 热运动占主导地位时, 自旋的排列方向随机, 使得
后强度为零. 存在一个临界温度 𝑇
𝑐
𝑇 < 𝑇
𝑐
: 自发磁化, 铁磁相, 𝑀 0
𝑇 > 𝑇
𝑐
: 无自发磁化, 顺磁相, 𝑀 = 0
随着温度的升高, 磁化强度连续低降低为. 称高温顺磁相称为无序相, 低温铁磁相为序相. 序相在
空间的任何位置都没有区别, 而有序相则在某些方向上有特殊性, 称为对称性破缺. 使得对称性破缺的外
场称为对称性破缺场, 不是由外场引起的对称性破缺称为自发对称性破缺
磁化强度反映了系统的对称性, 选取它作为序参量
无序相具有高对称性, 序参量为零
有序相具有低对称性, 序参量不为零
当系统从有序相趋近于临界点时, 序参量连续地变为零
6.2 临界点的自由能
在临界点附近, 𝑀 是小量. 由于二级相变中体积不变, 略去体积. 朗道假定自由能可以展开
𝐹 ( 𝑇 , 𝑀) = 𝑎
0
(𝑇) + 𝑎
1
(𝑇)𝑀 + 𝑎
2
(𝑇)𝑀
2
+ 𝑎
3
(𝑇)𝑀
3
+ · · ·
在没有外场 𝐻 = 0 , 磁化强度的方向不应影响自由能, 所以 𝐹(𝑇 , 𝑀) 𝑀 的偶函, 这要求奇磁幂
系数为零,
𝐹 ( 𝑇 , 𝑀) = 𝑎
0
(𝑇) + 𝑎
2
(𝑇)𝑀
2
+ 𝑎
4
(𝑇)𝑀
4
+ · · ·
𝑎
0
(𝑇) 实际上就是 𝑀 = 0 时的自由能. 两相稳定平衡时要求自由能极小,
𝜕𝐹
𝜕𝑀
𝑇
= 0,
𝜕
2
𝐹
𝜕𝑀
2
> 0
略去四次以上的高次项, 代入得第一式
𝜕𝐹
𝜕𝑀
𝑇
= 𝑀
2𝑎
2
+ 4𝑎
4
𝑀
2
= 0
其解为
𝑀 = 0, 𝑀 = ±
r
𝑎
2
2
𝑎
4
为零考察其稳定性, 代入第二式得到
𝜕
2
𝐹
𝜕𝑀
2
= 2𝑎
2
+ 12𝑎
4
𝑀
2
> 0
代入解得到
𝜕
2
𝐹
𝜕𝑀
2
= 2𝑎
2
, 𝑀 = 0
𝜕
2
𝐹
𝜕𝑀
2
= 4𝑎
2
, 𝑀 = ±
s
𝑎
2
(𝑇 𝑇
𝑐
)
2𝑎
4
在前面的分析中, 𝑇 > 𝑇
𝑐
𝑀 = 0, 𝑇 < 𝑇
𝑐
𝑀 0, 因此稳定性要求
𝑎
2
> 0, 𝑇 > 𝑇
𝑐
𝑎
2
< 0, 𝑇 < 𝑇
𝑐
又有稳定性条件要求四阶导数大于零, 所以在一阶近似下, 不妨选取
𝑎
1
(𝑇) = 𝑎
20
(𝑇 𝑇
𝑐
), 𝑎
4
(𝑇) = 𝑎
4
, 𝑎
20
> 0, 𝑎
4
> 0
那么得到稳定解
𝑀 (𝑇) =
0, 𝑇 > 𝑇
𝑐
±
r
𝑎
20
(𝑇 𝑇
𝑐
)
2𝑎
4
, 𝑇 < 𝑇
𝑐
6.3
鉴于磁场变化时外界应做功
d¯𝑊 = 𝑉d
1
2
𝜇
0
𝐻
2
+ 𝜇
0
𝑉H · dM
前半部分是真空中的磁场能量变化, 后半部分是介质被磁化时增加的能量. 磁场排除在系统外, 有外
界对系统做功
d¯𝑊 = 𝜇
0
𝑉H · dM
进一步忽略系统体积的变化, 在各向同性的介质中 𝐻 𝑀 平行, 则热力学基本方程就写为
d𝑈 = 𝑇d𝑆 + 𝜇
0
𝐻d𝑀
进一步做勒让德变换得到 Helmholtz 自由能
𝐹
=
𝑈
𝑇 𝑆,
d
𝐹
=
𝑆
d
𝑇
+
𝜇
0
𝐻
d
𝑀
所以
𝑆 =
𝜕𝐹
𝜕𝑇
𝑀
代入 𝐹 的表达式
𝐹 ( 𝑇 , 𝑀) = 𝐹
0
(𝑇) + 𝑎
20
(𝑇 𝑇
𝑐
)𝑀
2
+ 𝑎
4
𝑀
4
得到
𝑆 =
𝜕𝐹
0
𝜕𝑇
+ 𝑎
20
𝑀
2
𝑆
0
(𝑇) + 𝑎
20
𝑀
2
进一步代入 𝑀 (𝑇) 的表达式, 得到
𝑆 =
𝑆
0
(
𝑇
)
,
𝑇 > 𝑇
𝑐
𝑆
0
(𝑇)
𝑎
20
(𝑇 𝑇
𝑐
)
2𝑎
4
, 𝑇 < 𝑇
𝑐
𝑇 < 𝑇
𝑐
时的熵比 𝑇 > 𝑇
𝑐
时的熵小, 也就是有序相的熵比无序相小
6.4 磁化率
考察相变平衡时的磁化率. 𝑇 = 𝑇
𝑐
, 𝐻 0 时磁化率依赖于磁场
𝑀 𝑀 (𝑇, 𝐻)
由于
d𝐹 = 𝑆d𝑇 + 𝜇
0
𝐻d𝑀
所以
𝜇
0
𝐻 =
𝜕𝐹
𝜕𝑀
𝑇
(保留至四阶, 并对二阶项做一阶近似)
𝐹 ( 𝑇 , 𝑀) = 𝐹
0
(𝑇) + 𝑎
20
(𝑇 𝑇
𝑐
)𝑀
2
+ 𝑎
4
𝑀
4
所以
𝜇
0
𝐻 = 2𝑎
20
(𝑇 𝑇
𝑐
)𝑀 + 4𝑎
4
𝑀
3
𝑇 = 𝑇
𝑐
, 𝑎
2
为零, 所以
𝑀 =
3
r
𝜇
0
𝐻
4𝑎
4
磁化率定义为
𝜒 =
𝜕𝑀
𝜕𝐻
𝑇
此时只是保持 𝑇 不变, 并不是说 𝑇 = 𝑇
𝑐
, 所以不能使用 𝑇 = 𝑇
𝑐
时的 𝑀. 直接对方程两边对 𝐻 微分得到
𝜇
0
= 2𝑎
20
(𝑇 𝑇
𝑐
) 𝜒 + 4𝑎
4
𝑀
2
𝜒
所以
𝜒 =
𝜇
0
2𝑎
20
(𝑇 𝑇
𝑐
) + 4𝑎
4
𝑀
2
再代入 𝐻 = 0 时的 𝑀
2
, 得到
𝜒 =
𝜇
0
2𝑎
20
(𝑇 𝑇
𝑐
)
, 𝑇 > 𝑇
𝑐
𝜇
0
4𝑎
20
(𝑇 𝑇
𝑐
)
, 𝑇 < 𝑇
𝑐