宏观热力学体系的相变临界行为一直是统计力学研究的核心问题之一。临界行为之所以重要和引人入胜,是因为体系在临界点处表现出较非临界点处非常奇特的性质:内部关联长度发散导致体系相较于其它非临界点处出现巨大的热涨落,进而导致刻画体系特征的各种响应函数(对应二阶矩,与涨落有关)在该点呈现出奇异性(发散或间断点等非解析性)。
本文将以气液连续相变,平均场近似下Ising模型的连续相变,以及更为一般的Landau关于连续相变的唯象理论——二级相变理论为例,具体计算在每种情形下各种热力学量和响应函数的临界行为,这个临界行为由一组特征的临界指数标定。然后通过比较这些临界指数我们可以发现这些看似完全不同的体系和处理手法在相变临界点附近呈现出完全相同的物理行为,并由此揭示出相变临界行为的普适性。
最后探讨这种普适性与体系内部对称性间的深刻联系。
在考虑了分子间的相互作用势能后(即近端排斥远端吸引),我们可以对原始的无分子间相互作用的理想气体状态方程中的压强和体积做修正。修正后最简单的物态方程即为范德瓦尔斯方程。现考虑1mol非理想气体所满足的如下的范德瓦尔斯方程。其中参数的值依赖于具体物质(比如水分子有水分子的和,氢气分子有氢气分子的和)。我们将利用该方程来研究在临界点处气液连续相变的问题。
因为相变临界点对应于图上的拐点,所以此时压强对体积的一阶和二阶偏导数均为0。容易解出临界点处的体积,温度和压强分别是:以临界点的状态值作为单位将原始的范德瓦尔斯方程改写成与具体物质无关的普适方程:为直观起见,图1给出了上述普适方程在不同约化温度下的函数图像:图1函数图像。其中纵坐标是约化压强,横坐标是约化体积。
蓝线代表约化温度,对应相变临界点情形,蓝线上面的绿线和红线代表,下方的紫线和黑线代表在临界点附近对约化压强,体积,和温度作小量展开到一阶:其中是临界点附近的小量。所以普适方程可以写成关于小量的形式:为直观起见,图2给出了上式在不同下的函数图像:图2函数图像。其中纵坐标是约化压强在临界点附近展开对应的小量,横坐标是约化体积在临界点附近展开对应的小量。
黑线代表在临界温度时的关系,红线和紫线代表临界温度以上的关系,蓝线和绿线代表临界温度以下的关系。
将上面关于小量的普适方程在附近泰勒展开:从图2可以看出,在临界温度以下时,即,由于曲线的非单调性(即关于的多值性),所以体系必然存在气液相变。在相变发生时,根据两相平衡条件,相和相(这里分别指代液相和气相)的温度,压强,和化学势(吉布斯自由能)必须相同。
根据两相压强相等的力学平衡条件可以得到:根据等温线上两相化学势相等的化学平衡条件可以得到:将间的小量关系代入上式得到:若假设,则上述方程约化成:所以意味着和同号,与图2()矛盾。所以先前的假设错误。正确的关系只能是,也就是(关于纵轴对称)。注意到这个关系对应到图1就是和关于临界体积呈对称分布。
接下来将化学平衡的结果代入到力学平衡的结果得到:这个结果对应于从临界温度下方趋于临界点时,气相和液相间存在体积差,也就是两相间存在密度差(即两相间不对称,处于有序状态!)。但当我们从临界温度上方趋于临界点时,容易看出此时由于曲线单调,为了同时满足热学平衡,力学平衡和化学平衡条件,只能取,对应于此时的气相和液相间不存在体积差,也就是两相间不存在密度差。
所以严格来说,高于临界温度时其实只有一相(叫气相也行,叫液相也行,因为两相间不存在密度差所以两相完全对称,处于无序状态!无法区分谁是气体谁是液体)。所以密度差是个关键的指标,它表征了体系当前进入到了哪种序。我们把这个关键的指标叫做“序参量”!在我们现在考虑的气液相变的例子里,两相的密度差就可以充当体系的序参量。
我们下面看它的临界行为。当从临界温度下方趋于临界点时:当从临界温度上方趋于临界点时:因为,所以。所以如图3所示,当从临界温度下方趋于临界点时,序参量(密度差)是以的方式趋于0的;当从临界温度上方趋于临界点时,序参量(密度差)恒为0。所以我们得出在气液相变问题里序参量的临界指数是1/2!因为此时序参量(密度差,基本对应吉布斯自由能的一阶导数)在临界温度处连续,所以我们也称这种相变叫“连续相变”。
下下面看等温压缩率的临界行为:所以我们发现不管是从临界温度下方还是上方趋于临界点,等温压缩率都是以的方式趋于无穷大(发散的)!所以我们得出在气液相变问题里等温压缩率的临界指数是1!因为等温压缩率基本上是体积对压强的一阶导数,而体积又基本上是吉布斯自由能对压强的一阶导数,所以等温压缩率基本上是吉布斯自由能对压强的二阶导数。也就是说在临界温度处,吉布斯自由能对压强的二阶导数出现奇异性(发散)!
吉布斯自由能对热力学参数的二阶导数一般对应到体系的响应函数。所以此时响应函数在临界点附近出现奇异性。这种连续相变一般被称作“二级相变”。
下面看等容比热的临界行为(由于此时气液相变系统的所有信息是以关于、、的物态方程而不是以自由能或配分函数的形式给出的,所以这里等容比热的计算相当繁琐,这里直接给出结论):所以我们发现不管是从临界温度下方还是上方趋于临界点,等容比热都是以的方式发生奇异性(间断点)的!所以我们得出在气液相变问题里等容比热的临界指数是0!
根据先前关于小量的物态方程易知:当固定在等温线(即)时,偏离压强和体积临界值的小量间满足三次方的关系:或者用原始的量可以把关系表示成:这里的临界指数3也是标定气液相变临界行为的一个特征指标!
考虑一个任意维度下的Ising模型。在平均场近似下,我们把Ising形式的相互作用等效成了外磁场的一部分,使得体系被近似成一个等效的“无相互作用”体系。
在统计力学的框架下通过对Ising形式正则配分函数的操作,容易得到关于平均磁化的自洽方程。当外场时(也就是考虑自发磁化),临界点附近的平均磁化是个小量。所以将自洽方程关于小量展开得到:其中特征温度。上述自洽方程其中一个解是,其余两个解是否存在取决于以下方程右侧是否大于0:当时,方程右侧小于0,故原方程只有这个解;当时,方程右侧大于0,故原方程除了有的解,还有这两个解。
由于对应的Helmholtz自由能比时的自由能要高,所以是非稳定态必须舍去。所以当体系从临界温度下方逼近临界点时,的临界行为是:当体系从临界温度上方逼近临界点时,。所以如图6所示,当从临界温度下方趋于临界点时,序参量(平均自发磁化)是以的方式趋于0的;当从临界温度上方趋于临界点时,序参量(平均自发磁化)恒为。所以我们得出在平均场近似下Ising模型里序参量的临界指数是!
这个结果和之前气液相变里序参量(密度差)的临界指数完全一样!
下下面看等温磁化率的临界行为。由于此时需要考虑磁场的变动,所以我们需要使用原始带的自洽方程:所以体系的磁化率是:根据之前的结论,当外场趋于0且从临界温度下方趋于临界点时,所以:当从临界温度上方趋于临界点时,所以:所以我们发现不管是从临界温度下方还是上方趋于临界点,等温磁化率都是以的方式(类似宏观顺磁性的居里定律)趋于无穷大(发散的)!
所以我们得出在平均场近似下Ising模型里等温磁化率的临界指数是1!这个结果和之前气液相变里等温压缩率的临界指数完全一样!
下下面考虑无外磁场时比热的临界行为。
思路是从Ising体系的正则配分函数出发求得内能,然后通过内能在临界点附近的行为推断出比热的临界行为:当从临界温度下方趋于临界点时,所以:所以当从临界温度下方趋于临界点时比热是:当从临界温度上方趋于临界点时,所以容易看出比热是:所以我们发现不管是从临界温度下方还是上方趋于临界点,比热都是以的方式发生奇异性(间断点)的!所以我们得出在平均场近似下Ising模型里比热的临界指数是0!
这个结果和之前气液相变里等容比热的临界指数完全一样!
最后看固定在临界温度上,当外磁场很弱的情况下,和(或)的关系:因为在临界点附近且很小,所以也很小。所以下里的宗量很小。所以方程右端可以展开成所以如图8所示,当固定在临界温度时(且外场很小时),磁场和平均磁化或平均磁矩间满足三次方的关系:或。这里的临界指数3也是标定平均场近似下Ising模型临界行为的一个特征指标!
这个结果和之前气液相变里偏离压强和体积临界值的小量间满足的关系所对应的临界指数3,即,完全一样!
从上面的分析中我们发现气液相变,平均场近似下的Ising模型,和更为抽象一般的Landau二级相变理论都给出了完全相同的一组临界指数!
看起来相互间没有任何关系的体系(气液相变,Ising铁磁相变体系,和Landau二级相变理论所描述的更为抽象一般的体系),出发点和处理手法也完全不同(气液相变里是从非理想气体所满足的物态方程出发,Ising铁磁相变体系是从平均磁化满足的自洽方程出发,Landau二级相变理论是从作为序参量的多项式函数的Landau自由能出发),竟然在临界点附近的行为表现出惊人的一致性!!!
这种一致性产生的根本原因在于:气液连续相变,Ising铁磁相变体系,和Landau二级相变所描述的更为抽象的体系具有完全相同的对称性,即所谓的对称性!所以它们所描述的体系处于完全相同的普适类里!如果把普适类看成是一类几何/拓扑对象,那么这些临界指数就可以看成是标定这个几何/拓扑对象特征的“拓扑不变量”!
这些“拓扑不变量”和具体物质(可以类比成对象局域的几何性质)无关,而只与体系的对称性(可以类比成对象本身的拓扑性质,比如亏格和缠绕数)有关,具有普适性。
然而值得注意的是,这个普适性也是有限度的。比如这三套理论对于比热临界指数的估计都是0(这与实验中测出的非零的比热临界指数差别很大),这是因为它们完全忽略掉了体系本身的热涨落效应,所以它们只能给出一种特殊的普适类。对于具有别的对称性的体系(比如二维XY模型),它将属于别的普适类。别的普适类将由别的一个组临界指数标定。
所以自然界中千差万别的体系可以被划分进不同的普适类里,每个普适类内部的所有系统在临界点附近都是完全等价的,所以只要研究清楚普适类里其中一个系统在临界点处的行为,我们就可以直接推断出别的和它在同一个普适类里的其它系统在临界点处的行为!所以对不同的具体系统的物理性质的研究可以转化为对不同普适类的抽象数学性质的研究!而这也是相变和临界现象问题里最引人入胜的地方之一。