跳转到内容
← 返回知识库
热力学约 20 分钟阅读

相变与临界现象

相变临界现象序参量普适类伊辛模型重整化群超导超流临界乳光

关键词

相变; 一级相变; 连续相变; 序参量; 临界点; 临界指数; 伊辛模型; 普适性; 朗道理论; 重整化群; 超导相变; 玻色-爱因斯坦凝聚; 相图

"相变"不只是水的三态变化

"相变"在中学物理中以水的三态转变(固、液、气)为典型例子介绍。但现代物理学中,相变的概念远比这广泛:超导体中(金属态→超导态)、铁磁体中(顺磁态→铁磁态)、液晶中(各向同性→向列相)、超流氦中(正常流体→超流体)……几乎每一种涉及自发对称性破缺(spontaneous symmetry breaking)的转变都是相变。

关键区别在于序参量(order parameter):描述新相秩序程度的宏观量,在高温相为零,在低温相非零。水的固液相变序参量是晶体密度波的振幅;铁磁相变序参量是自发磁化强度;超导相变序参量是超导凝聚波函数的振幅(复数!);宇宙早期的电弱对称性破缺序参量是希格斯场的真空期望值。希格斯机制(产生基本粒子质量的机制)本质上是宇宙尺度的"相变",与铁磁相变在数学结构上深度类似——安德森(Philip Anderson)在1962年首先指出这一类比,为温伯格(Steven Weinberg)和萨拉姆(Abdus Salam)在1967—1968年建立电弱统一理论提供了关键思路。

日常说的冰融、水沸,多半是一级相变:有潜热、有两相共存,序参量可以跳变。临界点上的液气转变、铁磁居里点、常规超导转变,则是连续相变:没有潜热平台,关联长度发散,涨落自己变成主角。分类写在自由能与化学势的导数上,也写在有没有潜热这条实验判据上。更偏日常的三态与克拉佩龙方程,见相变;本文要追的是临界点附近那一套指数与普适性。

安德鲁斯的CO₂实验:临界乳光让观察者震撼

临界现象的实验发现可以追溯到托马斯·安德鲁斯(Thomas Andrews)在1869年对CO₂液气相变的研究。安德鲁斯发现:在约31°C(临界温度)、约73.8 bar(约72.8 atm,临界压强)以上,CO₂的液态和气态之间不再有明显界面,两者合并为单一的"超临界流体"(supercritical fluid)。临界点附近,CO₂散射光线极强(临界乳光,critical opalescence),使原本透明的气体变为浑浊——这是密度涨落在所有尺度(从分子到厘米量级)同时显著的结果,在当时令观察者震撼。

范德华(Johannes van der Waals)在1873年给出了方程 (p+a/V2)(V−b)=RT(p + a/V^2)(V - b) = RT,能描述从气态到液态的连续变化,预测了临界点的存在(Tc=8a/(27bR)T_c = 8a/(27bR),pc=a/(27b2)p_c = a/(27b^2),Vc=3bV_c = 3b),获1910年诺贝尔物理学奖。但范德华理论给出的临界指数(如密度差 Δρ∝∣T−Tc∣1/2\Delta\rho \propto |T - T_c|^{1/2},指数 β=1/2\beta = 1/2)与实验不符(实验值约0.326)——这一偏差困扰物理学界超过半个世纪,直到重整化群(1971—1972年)才得到解释。

临界乳光不是灰尘把光挡住,而是流体自己的密度涨落把光散射出去。一九〇八年,马里安·斯莫卢霍夫斯基把浑浊归因于临界点附近压缩率发散:同样的热扰动,在几乎无限柔软的流体里可以堆出接近可见光波长的密度块。爱因斯坦一九一〇年把这件事接到瑞利散射的定量公式上:散射强度随等温压缩率一起涨,临界点处若还假定各体积元彼此独立,强度会趋向无穷。

奥恩斯坦与泽尔尼克一九一四年补上空间关联:涨落在距离上并不独立,关联函数有一个随温度变长的尺度。有了关联,散射在精确临界点仍保持有限,角度分布却会变尖。乳光因此是关联长度涨到光波长附近时的直接视觉后果,不是杂质,也不是气液界面的雾。同一套涨落语言也解释了洁净空气为何仍能让天空发蓝,只是远离临界点时幅度小得多,见瑞利散射与天空的颜色。

伊辛模型(Ising model)是相变理论的最重要格点模型:每个格点上的自旋 si=±1s_i = \pm 1,相邻格点相互作用能为 −Jsisj-Js_is_j,哈密顿量 H=−J∑⟨ij⟩sisj−h∑isiH = -J\sum_{\langle ij\rangle}s_is_j - h\sum_i s_i($h$ 为外磁场)。一维伊辛模型无相变(恩斯特·伊辛,Ernst Ising,1925年,他的博士论文结论);二维伊辛模型在 TcT_c 处有相变,由昂萨格(Lars Onsager)在1944年精确求解,是统计力学中最重要的精确解之一,证明了相变在热力学极限下的数学存在性。

临界指数、标度律与普适性:为何不同系统在临界点"殊途同归"

在临界点 TcT_c 附近,各物理量以幂律形式发散或趋零,指数称为临界指数(critical exponents):

  • 序参量:m∝(Tc−T)βm \propto (T_c - T)^\beta(T<TcT < T_c)
  • 比热:C∝∣T−Tc∣−αC \propto |T - T_c|^{-\alpha}
  • 关联长度:ξ∝∣T−Tc∣−ν\xi \propto |T - T_c|^{-\nu}(在临界点 ξ→∞\xi \to \infty)
  • 磁化率:χ∝∣T−Tc∣−γ\chi \propto |T - T_c|^{-\gamma}
  • 关联函数在临界点:G(r)∝r−(d−2+η)G(r) \propto r^{-(d-2+\eta)}($d$ 为空间维数)

这些指数不是独立的,满足标度律(scaling relations):α+2β+γ=2\alpha + 2\beta + \gamma = 2(Rushbrooke不等式,在临界点变等式),γ=ν(2−η)\gamma = \nu(2 - \eta),dν=2−αd\nu = 2 - \alpha(超标度关系,hyperscaling)等。这些关系来自临界点处标度不变性的假设。

普适性(universality):相同空间维数和序参量对称性的系统,无论微观细节,具有相同的临界指数——称为同一普适类(universality class)。例如,三维伊辛普适类(Z2\mathbb{Z}_2对称性,d=3)包括:液气相变(CO₂、Xe等)、单轴铁磁体(Fe等)、铁电体等,它们的临界指数均为 β≈0.326\beta \approx 0.326,ν≈0.630\nu \approx 0.630,γ≈1.237\gamma \approx 1.237。普适性是重整化群的最深刻预言,意味着微观细节在临界点处"无关紧要"(被称为"无关变量"被RG流驱动到零)。

标度律把上面那组指数锁在一起,不是又测出一套新数字。一九六五年,威德姆假定自由能是齐次函数,温度偏离与外场按特定幂次同时缩放,Rushbrooke 等式就从热力学不等式变成临界点上的等式。卡达诺夫一九六六年把格点自旋打成块:块与块的有效耦合随尺度跑,临界点是跑动里的不动点。微观晶格怎么切、次近邻有多强,只要对称性和维度不变,都会流到同一类,指数因此共用,见伊辛模型。

这里沿用已经写过的三维伊辛类数值:β 约 0.326、ν 约 0.630、γ 约 1.237,平均场则给出 β=1/2。不要把它们读成又测出了一批新指数——它们是同一套实验与数值共识。

金兹堡判据说明平均场何时失效:涨落修正相对朗道自由能不可忽略的温度窗口,在三维伊辛类里相当宽,所以范德华和朗道都会把 β 报成 1/2,实验却停在约 0.326。四维及以上涨落变弱,平均场才定量可靠。动力系统里的流,在这里就是耦合常数随尺度的跑动。

临界乳光的强度随温度的变化,测的其实是压缩率,也就是指数 γ。不必另报一个新指数:三维伊辛类里 γ 约 1.237,平均场是 1。乳光变亮的速度,和磁化率在居里点附近发散的速度,应当对上同一条标度。

实验室里真正棘手的是有限尺寸:样品再干净,关联长度也不能超过容器尺度,幂律会被截断成平滑的峰。有限尺寸标度把峰高、峰宽与线性尺寸绑在一起,外推到热力学极限后,才能把测到的有效指数与普适类预言对照。高压流体池、超薄铁磁膜与光镊俘获的胶体晶体,走的都是这条外推路径,而不是「温度扫一遍就读出 β」。

超导、BKT相变与拓扑序:朗道范式之外的世界

超导相变是量子相变的典范。BCS超导(约翰·巴丁、莱昂·库珀、约翰·罗伯特·施里弗,1957年,共获1972年诺贝尔物理学奖):低温下,晶格振动(声子)在费米面附近的电子之间产生微弱的吸引相互作用(库珀对,Cooper pairs),使一个BCS序参量(具有时间反演对称性的配对波函数)在某一临界温度 TcT_c 以下变为非零——发生连续相变,系统进入超导态。

高温超导体(high-temperature superconductors,HTS)自1986年柏诺兹(Georg Bednorz)和缪勒(K. Alex Müller)发现铜氧化物超导体(Tc≈35T_c \approx 35 K,后续材料提升至~134 K)以来,其配对机制至今仍是凝聚态物理中最大的未解问题之一。铜氧化物超导体中的超导配对可能来自电子-电子强关联(不是声子),"赝能隙"(pseudogap)相变、反铁磁序与超导的竞争是当前研究热点。

BKT(Berezinskii-Kosterlitz-Thouless)相变是二维系统中一种奇特的相变,没有通常意义下的局部序参量(Mermin-Wagner定理禁止二维系统在有限温度下有长程序),但通过涡旋-反涡旋对的束缚/解离,存在从"准长程序"到"短程序"的拓扑相变。科斯特利茨(John Kosterlitz)和撒利斯(David Thouless)因此与霍尔丹(Duncan Haldane)共获2016年诺贝尔物理学奖,这是拓扑物理进入主流物理学的标志性时刻。

朗道理论假定自由能可以按序参量做解析展开,临界点只是展开系数改号。金兹堡指出,在三维,关联长度发散会把这个展开冲垮,平均场指数因此只是远离临界点的近似。连续相变的正确分类后来改由对称性、维度和是否有长程相互作用来指定普适类,而不是只看自由能多项式的几次方。

还有一类转变根本写不出局域序参量。二维的 BKT 靠涡旋束缚,某些拓扑相靠整体不变量,朗道的“有序程度从零长出来”在这里用不上。超导本身仍可用复数序参量,高温铜氧化物却可能同时经历赝能隙、电荷序与配对,竞争让单一朗道序不够用。机制细节见超导体,这里只强调:相变语言大于朗道范式,不等于每一个转变都有一个可局域测量的序参量。

实验上区分「朗道式连续相变」与「拓扑解离」的一条硬判据,是关联函数的衰减形态:前者在临界点呈幂律并伴随比热奇异,后者在整个准长程有序相内都是代数衰减,比热可以几乎无奇异。氦膜超流、二维约瑟夫森结阵列与冷原子中的二维玻色气体,都用这条判据把 BKT 从普通连续相变里抠出来。

从宇宙早期的电弱相变到量子计算机:相变语言的无处不在

宇宙学中的相变是宇宙演化的关键节点。大爆炸后约 10−3510^{-35} 秒,宇宙经历暴胀(inflation)——可能与某个标量场的相变相关;约 10−1210^{-12} 秒,温度降至约100 GeV,电弱对称性自发破缺(SU(2)×U(1)→U(1)emSU(2)\times U(1)\to U(1)_{\rm em}),W和Z玻色子获得质量,希格斯场取得非零真空期望值;约 10−610^{-6} 秒,夸克-强子相变(QCD相变),自由夸克和胶子"束缚"成质子、中子等强子;约20分钟,大爆炸核合成;约38万年,质子俘获电子形成中性氢,宇宙变透明(CMB的起点)。这些宇宙学相变在理论上与凝聚态物理中的相变使用相同的数学语言(有效场论、序参量、对称性破缺),这是当代基础物理学"有效场论普遍性"的深刻体现。电弱破缺的序参量是希格斯真空期望值,数学骨架与铁磁相变同类,见希格斯机制。暴胀是否对应一次真正的热相变,模型之间并不统一,不宜写成已经定论的宇宙史年表。

拓扑绝缘体(topological insulators)和拓扑超导体是2000年代以来凝聚态物理最重要的发现之一。这些材料中存在由时间反演对称性或其他离散对称性保护的拓扑相,其相变不能用朗道序参量理论描述(拓扑不变量的改变),体现了超越朗道范式的"拓扑相变"(topological phase transition)。拓扑超导体可能托管马约拉纳费米子(Majorana fermions,自身的反粒子)激发,是拓扑量子计算(topological quantum computing)的候选平台,对噪声具有天然抗干扰能力(因为拓扑保护)——微软于2025年2月发布的"Majorana 1"拓扑量子处理器正追求这一路线(其马约拉纳准粒子的确证仍受到学界质疑)。

活性物质(active matter)是近年来统计物理的新兴领域:细菌群落、鸟群(椋鸟群舞,murmurations)、细胞骨架网络等由大量自驱动粒子组成的系统,展示出类似相变的集体行为(从无序到有序游动的"活性Ising模型")。这些系统在根本上是非平衡的,但仍可以用相变语言(序参量、临界现象、相图)描述,是统计力学扩展到生命科学的前沿领域。

连接节点:统计力学与玻尔兹曼分布 → 相变与临界现象(系综理论到临界现象);相变与临界现象 → 超导体(超导作为典型相变);相变与临界现象 → 宇宙暴胀(宇宙早期相变)

事实卡

  1. 卡1

    昂萨格(1944年)精确求解了二维伊辛模型,证明了热力学极限下相变的数学存在性,是统计力学最重要的精确解。

  2. 卡2

    高温超导体(1986年,柏诺兹和缪勒发现铜氧化物,Tc≈35T_c \approx 35 K)至今配对机制未解,是凝聚态物理最重要的未解问题之一。

  3. 卡3

    BKT拓扑相变(2016年诺贝尔奖)揭示了二维系统中没有局部序参量的拓扑相变,开创了拓扑物理的研究方向。

  4. 卡4

    宇宙学相变(电弱相变、QCD相变等)与凝聚态相变使用相同的有效场论语言,体现了物理学的深层统一性。

引用

"自然界的所有相变——从冰的熔化到磁铁的磁化——都是同一个深刻数学结构(对称性破缺)的不同表现。"
菲利普·安德森,《多则异》(More is Different,1972年,意译)

跨域连接

  • 伊辛模型:昂萨格对二维伊辛模型的精确解,首次在有相互作用的多体系统中严格证明了相变的存在,并给出比热对数发散这一平均场完全错过的结果。推论:精确解是所有近似方法的裁判——平均场在低维的系统性失败,正是被它钉死的。
  • 希格斯机制:电弱对称性破缺本质上是宇宙尺度的相变:希格斯场的真空期望值扮演序参量,与铁磁体的自发磁化共享同一数学结构。推论:粒子质量随"序参量"的出现而产生——若早期宇宙的相变次序不同,基本粒子的质量谱就会不同。
  • 动力系统:重整化群把"改变观测尺度"变成耦合常数空间里的一个流:参数沿流向不动点,无关方向在流中被甩掉,剩下的决定临界指数。推论:同一普适类系统的临界指数必须一致——液气临界点与单轴铁磁的指数相等,是这条数学机制最惊人的实验兑现。
  • 分形:临界点处关联长度发散,涨落结构在所有放大倍数下看起来都一样——自相似正是分形的定义性特征。推论:临界系统的关联函数必须呈幂律而非指数衰减——幂律就是"没有特征尺度"的数学签名,拟合出指数律即可排除临界点。
  • 活性物质与软物质:鸟群、菌群与细胞骨架从无序到有序游动的转变,可借用序参与相图语言描述,尽管它们本质是非平衡系统。推论:若集体运动确属某个普适类,临界涨落的标度行为应与平衡相变定量相同——这是可检验、也常被检验的预言。

参考文献

  1. Onsager, Lars. "Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition." Physical Review 65 (1944): 117—149.(二维伊辛模型的精确解)
  2. Wilson, Kenneth G.; Kogut, John. "The Renormalization Group and the Epsilon Expansion." Physics Reports 12 (1974): 75—200.(重整化群的经典综述)
  3. Goldenfeld, Nigel. Lectures on Phase Transitions and the Renormalization Group. Westview Press, 1992.(相变与RG最清晰的教材)
  4. Kosterlitz, J. M.; Thouless, D. J. "Ordering, Metastability and Phase Transitions in Two-Dimensional Systems." Journal of Physics C 6 (1973): 1181—1203.
  5. Bednorz, J. G.; Müller, K. A. "Possible High Tc Superconductivity in the Ba-La-Cu-O System." Zeitschrift für Physik B 64 (1986): 189—193.
  6. 汪志诚. 《热力学·统计物理》(第五版). 高等教育出版社, 2013.
  7. Smoluchowski, M. "Molekular-kinetische Theorie der Opaleszenz von Gasen im kritischen Zustande, sowie einiger verwandter Erscheinungen." Annalen der Physik 330 (1908): 205—226.
  8. Einstein, A. "Theorie der Opaleszenz von homogenen Flüssigkeiten und Flüssigkeitsgemischen in der Nähe des kritischen Zustandes." Annalen der Physik 338 (1910): 1275—1298.
  9. Ornstein, L. S.; Zernike, F. "Accidental deviations of density and opalescence at the critical point of a single substance." Proceedings of the Royal Netherlands Academy of Arts and Sciences 17 (1914): 793—806.
  10. Widom, B. "Equation of State in the Neighborhood of the Critical Point." Journal of Chemical Physics 43 (1965): 3898—3905.
  11. Kadanoff, Leo P. "Scaling Laws for Ising Models near \(Tc\)." \_Physics Physique Fizika 2 (1966): 263—272.