材料怎样改变状态:相变、临界现象与标度

冰融化时吸热,某些共混物接近分相温度时变得浑浊,磁体升温后会失去自发磁化。这些变化怎样从许多分子的运动中产生?从相平衡和序参量出发,逐步建立自由能、涨落、关联长度与标度的联系,并用可操作的自由能模型比较连续转变和不连续转变。

看见变化,也要说清改变的是什么

冬天把冰块放进温暖房间,它先升温,随后一部分变成液态水。在冰和水共存的阶段,继续输入的能量可以主要用于改变状态,温度并不必然持续上升。另一类变化更难凭外形辨认:一杯原本透明的混合液接近某个温度时变得乳白,或一块磁体升温后不再保持原来的自发磁化。相变研究把这些现象放进同一个问题:大量粒子相互作用时,为什么宏观性质能够在特定条件附近出现突变,或者出现异常强烈而连续的变化?

相是材料中具有相同宏观状态的一种均匀组织。冰和水的化学组成可以相同,区别在于分子的排列与运动方式;两个聚合物富集相的区别可以主要在组成,外观上都仍然是液体。一块多晶材料含有许多取向不同的晶粒,它们也可能属于同一种晶相。因此,相的判断要结合组成、对称性、结构和热力学性质,不能简单按“看见几块”来数。界面常使两相更容易辨认,但在临界点附近,两相的差别和界面都可以逐渐消失。

相平衡意味着宏观上没有持续的净热流、净体积改变或净组分迁移。两相温度相等,避免热量持续单向传递;机械平衡需要相应的应力平衡,平直界面且忽略外场时常简化为压力相等;允许交换的每一种组分还要在两相中具有相同化学势。化学势表示增加少量该组分所改变的自由能,像一种决定迁移方向的“每份物质的代价”。这个比喻对应自由能差,而不是分子真的在比较价格。弯曲界面可以存在毛细压力差,因此压力条件也要随几何说明。

平衡并不要求两个相密度相同,或所含物质量相同。一杯水中漂着一小块冰,两相在合适条件下可以保持平衡,尽管质量、体积和密度都不同。温度、压力是不会随样品分成两份而简单减半的强度量;总能量、总熵和物质量则是会随规模变化的广延量。相平衡约束的是合适的强度量,样品总量和组成守恒决定两相各有多少。分清这两类量,能避免把“平衡”误解成一切数值都要相等。

给出温度和压力也不总能唯一决定任意复杂样品的全部状态。多组分体系还需要组成信息,带有外场、应变或固定网络约束的材料还要指明这些条件;在两相共存线上,需要额外的总组成或总量约束来确定相的比例。对于高分子,玻璃化和慢松弛还会使历史留下影响。后面讨论的平衡相变先假定材料有足够时间探索可达状态,再把实际测量的时间条件带回来。

这样建立的语言能让不同现象各归其位。熔化主要比较晶体与液体的自由能,分相比较不同组成分布的自由能,磁有序比较不同平均取向的自由能。真正共同的部分是:找一个能够区分状态的量,计算或测量它怎样随控制条件变化,再解释它背后的集体组织。相变理论从这里开始,而不是从给每条实验曲线贴上“一级”或“二级”的名称开始。

自由能的斜率怎样记录相变

在温度和压力固定、物质数量给定时,平衡状态具有最低的吉布斯自由能。想象分别测得冰和液态水的自由能随温度变化,两条曲线在平衡熔点相交:交点两侧,较低的那一条分别由不同相提供。真实平衡自由能沿着两条曲线的较低部分走,因此本身可以连续,却在交点改变斜率。这里“斜率”就是温度改变一点时自由能改变得有多快。把这个几何特征连接到可测量物理量,便能理解潜热为何出现。

自由能斜率对应熵与体积式 1
符号含义单位
S样品总熵J K⁻¹
G样品吉布斯自由能J
T绝对温度K
P压力Pa
V样品体积m³
n保持固定的各组分物质的量mol
∂G/∂T, ∂G/∂P对温度和压力的偏导数;下标指定保持不变的量分别为J K⁻¹和m³

怎样理解恒定组成的基本关系中,温度微变对应−S乘温度变化,压力微变对应V乘压力变化;按局部变化率读取系数得到这两个导数。

算一个例子若固定压力时,G在很小温区内每升高1 K减少20 J,则该状态S约为20 J K⁻¹。若定温加压1000 Pa使G增加0.001 J,则局部体积估算为10⁻⁶ m³。

使用条件平衡热力学、组成固定,仅考虑通常的压力—体积机械功;附加电场、应力等控制量需保持固定或加入对应项。

对温度求导得到熵的负值,对压力求导得到体积。求导不是额外的物理假设,而是把一小段变化除以引起它的控制量变化,再取越来越小的区间,得到当地变化率。下标提醒读者哪些量保持不变:测温度斜率时固定压力,测压力斜率时固定温度。自由能曲线的温度斜率突然改变,意味着熵发生跳变;压力斜率改变则意味着体积跳变。一级相变至少有合适的一阶导数不连续,并不要求熵与体积在每一个特殊情形都同时发生非零跳变。

熔融潜热把熵跳变变成直观的能量。等温、可逆地把一定量物质从一相变成另一相,吸收的热等于绝对温度乘熵变化。举例说,假想一种物质在300 K发生转变,每摩尔转变熵为10焦耳每摩尔每开尔文,所需热就是3000焦耳每摩尔。摩尔是计量粒子数的单位,一摩尔含有约6.022×10²³个指定实体;它让我们不必逐个数分子,也能比较相同数量物质的转变。这个算例只使用假想参数,说明温度、熵和热量怎样配合单位。

等温相变的摩尔潜热式 2
符号含义单位
Lm指定转变方向的摩尔潜热,吸热为正J mol⁻¹
Ttr平衡转变绝对温度K
Δsm末相减初相的摩尔熵变化J mol⁻¹ K⁻¹

怎样理解在平衡转变温度可逆地交换热时,熵变化等于吸热量除以温度,因此热量等于温度乘熵变化。

算一个例子Ttr=300 K、Δsm=10 J mol⁻¹ K⁻¹时,Lm=3000 J mol⁻¹。反向转变的Δsm与Lm都改变符号。

使用条件给定压力下的平衡、可逆、等温转变;区分每摩尔、每千克与样品总量,不能混用基准。

连续相变的序参量可以连续地趋向零,熵和体积等相关一阶导数在转变处保持连续,而响应函数表现出异常。例如热容衡量温度提高一点需要补充多少热;磁化率衡量施加一点场会改变多少磁化。它们对应自由能更高阶的导数,可能出现有限跳变、尖峰、幂律发散或对数异常。把所有连续转变都限定为“二阶导数有限跳变”,会错过真实临界行为。日常常用的“二级相变”可以指这类连续有序转变,但具体异常仍要逐个说明。

实验中看见一个宽峰,也不必马上认定转变是连续的。有限样品、温度不均、组分分布和测量速率可以把本来尖锐的转变展宽;反过来,连续转变附近也可能出现很尖的热容峰。一级转变在过热或过冷条件下还可能显示滞后,因为新相的形成要创造界面。滞后说明实验经过的路径与成核和松弛有关,本身不等于平衡转变温度。把加热和冷却曲线、速率变化和结构信号结合,能让分类建立在过程的物理原因上。

高分子玻璃化尤其需要这个区分。冷却时,链段运动变得很慢,结构来不及在实验时间内充分调整,材料表现得像固体;玻璃化温度会随时间尺度和热历史改变。它与本章讨论的平衡连续临界点有不同的测量条件。聚合物熔融、共混物分相和玻璃化可以在一段温度扫描中先后出现,但它们分别涉及有序结构、组成组织和动力学冻结。先认清改变的是哪一种状态信息,再使用自由能与响应函数解释曲线。

序参量把集体变化变成一个可读的量

要比较有序与无序,需要选一个能够反映两者差别的量。磁体中的局部磁矩可以偏向两个相反方向,把它们平均后,得到平均磁化;在两个等价子晶格上有序分布的合金中,可以比较两类原子占据不同位置的偏好;在对称二元混合物中,可以使用局部组成偏离一半的程度。这些量称为序参量。选择的依据是它能否区分所关心的相,而不是所有材料都必须使用同一个字母或同一种平均方式。

以允许正负两种等价有序状态的标量序参量m为例,m等于零表示平均对称,正负非零值表示两种互为对称的有序状态。在无外场时,把所有方向翻转不改变材料能量,因此自由能对m与−m应该相同。这个对称性会排除只含奇数次幂的项:把m换成−m时,一次和三次项会改变符号,二次和四次项保持不变。这里“标量”指用一个数字即可描述局部有序程度;具有多个方向分量或张量性质的序参量要采用相应的对称性分析。

Landau理论的思路,是在序参量足够小的范围,把自由能写成其逐次幂的和,保留能够描述局部形状的低阶项。二次项决定原点附近向上还是向下弯,正的四次项则防止自由能在较大序参量处无止境下降。温度可以改变二次项的系数,于是同一个多项式能够从一个最低点变成两个最低点。这个表达不需要先知道每个原子的精确位置,却把对称性与稳定性直接写进了可计算的图形。

正负对称的连续相变模型式 3
符号含义单位
f用固定正能量尺度归一化的自由能函数,省略与m无关的背景1
m无量纲标量序参量1
a正的二次系数比例因子1
b正的四次项系数1
t约化温度1
T绝对温度K
Tc此平均场模型的临界温度K

怎样理解正负对称允许偶次幂。at改变原点曲率,正b使较大幅度的自由能上升。除以Tc把温差变成可比较的无量纲距离。

算一个例子a=b=1、Tc=300 K、T=297 K时,t=−0.01,f=−0.005m²+0.25m⁴。m=0.1时f=−0.000025,比m=0处低。

使用条件无外场、正负对称标量序参量的均匀平均场展开,a和b在临界附近近似常数;空间涨落另行处理。

最低点处的斜率为零,因为往左或往右作极小移动都不能再降低自由能。对这个多项式求导,会出现一个m因子与另一括号的乘积;乘积为零意味着至少一个因子为零,于是得到原点和可能的两个非零候选。再检查曲率为正或比较自由能,就能确定哪些候选是真正的稳定最低点。当约化温度从正变负,原来的最低点变成最高点,两侧最低点连续地从零长出来。

最低点连续离开零的位置式 4
符号含义单位
meq平衡序参量;正负号代表两个等价状态1
t约化温度(T−Tc)/Tc1
a, b前式正的二次比例系数与四次系数1
√平方根,取非负根后再加正负号1
T, Tc绝对温度与临界温度K

怎样理解自由能的一阶导数为atm+bm³。令其为零并检查最低点,得到温度高于或等于Tc时的零解,以及低于Tc时的两个非零解。

算一个例子a=b=1时,t=−0.01给出meq=±0.1;t=−0.0025给出±0.05。到临界点的距离缩小4倍,序参量大小缩小2倍。

使用条件沿用前式的均匀平均场模型。等号t=0处二阶导数为零,但正四次项仍使零点成为最低点。

这种连续出现给出了平均场序参量随温度变化的平方根规律。如果把临界温度以下的距离缩小四倍,序参量就缩小两倍,而不是四倍。平方根说明集体有序的消失具有非线性。可是这个推导把空间中不断变化的局部状态浓缩为同一个平均量,忽略了长距离涨落。它适合建立机制和某些平均场范围内的定量描述;越来越靠近一些真实临界点时,需要让不同区域之间的关联重新进入理论。

对称性也要按照实际对象使用。向列液晶的分子头尾等价,指向矢翻转不改变物理状态;这并不意味着其标量有序度正负号对应完全相同的状态。完整的取向序参量是张量,其允许的三次不变量能使各向同性—向列转变出现不连续跳变。把“所有相变自由能都只能含偶数次项”当作普遍规则,就会丢掉这个机制。Landau方法的第一步始终是确定对象的对称性,再决定允许哪些项。

一个最低点如何让位于另一个

连续转变中,有序最低点从原点平滑出现;一级转变还可以通过另一种方式发生:原来的最低点仍在,远处却先出现新的最低点。随着控制量改变,两者相对深浅互换,平衡状态便从一个位置跳到另一个位置。自由能在两者等深时保持连续,序参量却可以跳变。把图形画出来,比只记住“一级有潜热”更能说明亚稳态、过冷和成核来自哪里。

仍然使用正负对称的标量,可以把四次项系数改为负,并加入正的六次项保证大幅度处的稳定。负四次项允许远离原点出现低谷,六次项则让曲线最终向上。这个模型中的转变不必等到原点失稳才发生;只要远处低谷首先变得更低,平衡态就已改变。原点如果仍向上弯,便可以暂时停留为亚稳态,直到出现足够的扰动或成核事件。

比较连续与一级转变的可调模型式 5
符号含义单位
f省略恒定背景后的无量纲自由能1
m无量纲序参量,可连续取值1
r可调二次项系数,在交互中独立控制1
u四次项系数,交互中选+1或−11

怎样理解正六次项稳定大幅度。u=+1时最低点在r穿过0时连续出现;u=−1时远处最低点可以先与原点等深,再取代原点。

算一个例子u=−1、r=3/16、m²=3/4时,三个贡献分别为9/128、−9/64和9/128,合计0,与f(0)相同;|m|=√(3/4)≈0.866。

使用条件无量纲均匀标量模型;r不是直接以K计的温度,需材料标定才能映射温度。交互比较平衡最低点,不计算空间成核或时间。

对这里选定的无量纲系数,四次项为负时,两类最低点在二次系数等于3/16处等深,非零最低点的序参量大小为√(3/4),约0.866。原点要到二次系数降到零才失去局部稳定性。因此,平衡转变和原点失稳之间存在一段明确的距离。直接数值代入也能核对:在3/16与精确的√(3/4)对应的位置,正二次项、负四次项和正六次项恰好合计为零,与原点自由能相同;使用0.866这个舍入值时会留下很小的舍入差。

图中的曲线只描述均匀序参量的自由能,不包含一个小新相核的表面。真实成核还要在空间中产生边界,既有新相体积带来的收益,又有界面面积带来的代价。足够小的核可能因表面代价而缩小,足够大的核才持续生长。温度低于平衡转变温度,却尚未观测到新相,便可以由这种势垒与有限等待时间解释;继续冷却会改变体积驱动力与成核速率。

绿色点:全局最低;橙色点:亚稳最低。曲线与系数均为无量纲。

同一个六次自由能模型,用四次项符号改变转变方式。滑块调整二次系数,按钮直接到达等深或连续分岔条件;反馈同时列出最低点位置及其全局/局部地位。

在交互中选择正四次项,缓慢把二次系数从正调到负,观察最低点怎样连续离开原点。再切换成负四次项,寻找三个最低点共存的区间:中间代表无序状态,两侧代表两种等价有序状态。标记为全局最低的点决定平衡选择,其余局部最低点表示可以暂存的亚稳状态。这里改变的是自由能函数的控制参数,并未模拟分子随时间跳跃,因此滑块运动速度不会改变图中平衡判定。

这个对照也说明,局部曲率、全局稳定性和转变级数是相关但不同的问题。正曲率保证一个状态能够抵抗附近的微小扰动,不保证它比所有遥远状态更有利;两个状态等深决定平衡共存,不说明样品需要多久才形成界面。相变的理论描述在这些层次之间逐步推进:先用自由能比较状态,再加入空间结构,最后加入时间演化。每一步回答一个清楚的问题,便不会让一张势能图承担它没有计算的结论。

从两种局部选择走向集体有序

Landau图景从平均有序程度出发,Ising模型则把镜头拉回许多相互作用的局部单元。想象一张规则网格,每个位置放一个只能取加一或减一的变量。这个变量可以代表简化的磁矩方向,也可以在另一种映射中代表两种组分的占据选择。相邻位置取相同符号时能量较低,热运动则允许局部翻转。模型简单到每个位置只有两个选择,却仍保留“邻居影响邻居”的集体机制。

最近邻Ising模型的能量式 6
符号含义单位
ℋ一个给定格点构型的总能量J
J相邻同向偏好的正相互作用能量;字母与焦耳单位同形J
h偏好正符号的外场耦合能,包含场与磁矩的乘积J
si, sj第i、第j个位置的两态变量1
i, j格点编号1
Σ⟨i,j⟩对每个最近邻对求和,每条键只计一次1
Σi对所有格点求和1

怎样理解同号乘积为+1,键能为−J;异号为−1,键能为+J。把每条键和每个格点的场贡献相加得到总能量。

算一个例子两个相邻位置且h=0时,(+1,+1)能量为−J,(+1,−1)为+J,相差2J。若J=10⁻²¹ J,差为2×10⁻²¹ J。

使用条件两态最近邻铁磁耦合模型。几何、边界条件和场条件决定具体结果;映射到组成变量时需相应的守恒约束。

式中的求和符号表示把指定范围内的贡献逐一加起来。最近邻对的求和每条相邻连线只数一次,避免同一相互作用重复。若相邻两个变量同号,乘积为加一,对能量贡献为负J;异号时乘积为负一,贡献为正J,因此一条键从同号改成异号需要增加2J能量。外场项偏好其中一种符号;这里把场与磁矩的耦合合并为具有能量单位的h,避免把能量参数直接误读成以特斯拉计的磁场。

能量低的状态更有统计优势,但温度不为零时,高能状态也可以出现。一个局部翻转可能增加邻居不一致的数目,因而增加能量;可供选择的不一致排列却往往很多,提供熵收益。低温时相邻一致的大片区域有利,高温时多种杂乱排列更容易占据统计权重。相变正是这两种趋势在大量相连位置中的集体结果,并非某一个位置独自跨越了固定温度门槛。

这个模型由Lenz在1920年提出的思路发展而来,Ising在1925年研究了一维最近邻链。对没有外场、相互作用有限且只沿一条无限长链连接的模型,任意非零温度都没有自发磁化。一个直观原因是,翻转一大片连续区域只需要在边界付出有限代价,而放置边界的位置有很多。升到二维后,区域边界的长度会随区域增大,几何竞争改变,有限温度有序便可以出现。维数由此成为物理条件,不能把一维结果直接推广到空间中所有模型。

1944年Onsager给出了零外场二维方格Ising模型的精确自由能结果,揭示真实的临界奇异性与简单平均场预测并不相同。对两个方向最近邻耦合相等的无限方格模型,临界热能约为相互作用强度的2.269倍。这个数字属于指定格子、指定耦合和零外场条件,改变邻居连接方式或相互作用范围,就会改变临界温度。更耐人寻味的是,不同微观模型在临界点附近又可能共享相同的指数,这为后面的普适性问题提供了入口。

二维连接让区域边界具有长度局部一致++++++++++++++++++++++++++++++++翻转内部区域++++++++++−−−−++++−−−−++++++++++边界上的异号键需要能量1944:Onsager求解零场二维方格模型
图 1 · 每个圆点代表一个两态变量,连线表示最近邻相互作用。翻转内部区域会在边界产生异号键,边界代价与二维几何有关;图示用有限网格解释机制,Onsager的热力学临界结果对应无限方格模型。

网格图中一个大块同号区域,并不意味着每个位置永久不动。局部翻转始终存在,区域边界也会移动。有限样品在零外场下甚至可能整体从正磁化切换为负磁化,使很长时间的有符号平均接近零。分析自发有序时,要说明取热力学极限和趋零外场的顺序,或使用绝对磁化、分布形状等适合有限样品的观测量。平均值是测量定义,不能只看到一个零就断定没有有序结构。

涨落怎样把远处的区域联系起来

温度固定并不意味着每个小区域始终拥有相同能量或组成。分子不断移动,一小块区域中的密度、组分比例或局部磁化都会围绕平均值变化,这种变化称为涨落。两个区域如果完全独立,一个偏高不提供另一个偏高的线索;如果相互关联,观察到一处偏高,就会改变对另一处状态的预期。关联函数便把这种“共同偏离平均”的程度按距离记录下来,让肉眼难以辨认的集体范围成为可测量的量。

关联函数与一个示意衰减长度式 7
符号含义单位
C(r)相距r两处的序参量关联1
δm局部m减其平均值的偏离1
m无量纲局部序参量1
𝐱, 𝐫位置矢量与间隔矢量m
r间隔矢量长度m
⟨…⟩对平衡样本或等价位置取统计平均1
C0示意指数衰减的幅度1
ξ关联长度m
e自然指数底数,约2.718281
∼表示所讨论近似的衰减形式,不要求对所有r精确相等1

怎样理解共同偏高或共同偏低的偏离乘积为正,独立涨落的乘积平均趋于零。指数形式示意相隔几个关联长度后联系迅速减弱。

算一个例子若ξ=10 nm,r=10 nm时指数因子e⁻¹≈0.368;r=30 nm时为e⁻³≈0.0498。ξ增为30 nm时,相隔30 nm的因子重新增为0.368。

使用条件平移统计近似均匀的平衡体系。纯指数是长度解释示意;完整临界关联可带幂次前因子,在临界点转为长距离幂律。

关联长度给出偏离能够在多远范围内保持显著联系。在远离临界点的许多体系里,超过这个距离后,关联迅速衰减;靠近临界点,越来越大的区域开始共同波动,关联长度增长。指数衰减中的e表示自然指数的底数,距离等于关联长度时,纯指数因子降到起始值的约0.368。这是理解长度的一种简单方式;临界附近的完整关联函数还可以带有距离的幂次因子,不能把所有关联都强行拟合成单一指数。

临界乳光把这件事变得可见。液体混合物中的组成涨落会引起局部折射率变化,当相应结构尺度和响应变得足够大,光被强烈散射,样品便出现浑浊。散射测到的并不是一层已经凝固的颗粒,而是材料内部随时间变化的密度或组成相关结构。通过不同散射角对应的空间尺度,可以估计关联长度以及长波长涨落强度。靠近临界点时散射增强,是涨落和宏观光学现象之间的一座桥。

接近临界点,关联延伸得更远00.510204060关联 / 起始关联间隔距离 / nmξ = 10 nmξ = 30 nm指数曲线用于解释长度;临界点需幂律
图 2 · 相同起始关联下,较大的ξ让相隔较远的区域仍有明显联系。图中比较纯指数衰减以解释长度含义;正文的临界关联还包含距离幂次,不能用这两条示意曲线替代完整临界函数。

涨落也与响应有关。一个体系若在没有外场时就容易出现大范围的自发偏离,施加很小的偏好通常就能显著改变其平均状态。磁化率增大、压缩响应增大和组成响应增强,都可以从合适的涨落—响应关系理解。这里比较的是指定平衡条件下、与外场共轭的变量;不同实验固定不同控制量,得到的响应函数也可能不同。把所测量变量和边界条件写清楚,才能将热容、压缩性或磁化率放到同一统计框架中。

平均场理论把每个位置的邻居替换为平均环境,在关联长度不大或涨落相对较弱时十分有效。临界附近,一个位置的偏离却会与许多区域共同出现,把邻居视为互不相关的平均背景就丢失了这种信息。问题并非平均值没有用途,而是平均值之外的关联也开始决定主导行为。空间维数、相互作用范围和观测距离都会影响什么时候必须保留这些涨落。

任何真实样品都是有限的。当关联长度接近样品厚度或观察窗口,继续增长便受到几何限制,理想无限系统的发散会变成有限的峰和圆滑的过渡。临界区域也常出现慢化:越来越大的关联区域需要更长时间才能重新调整。薄膜、微粒和短时间测量因而会各自截断临界行为。实验设计需要同时给出温度精度、样品尺度与等待时间,才能判断所观察的是临界趋势中的哪一段。

幂律如何比较不同临界现象

为了比较不同转变温度的材料,可以用相对距离代替直接用温差:温度减临界温度,再除以临界温度,得到无量纲的约化温度。这样,300 K附近偏离3 K,与600 K附近偏离6 K具有相同的相对距离。临界点附近,许多量的主要变化可写成这个距离的幂。例如序参量随负约化温度的某个正幂减小,关联长度随绝对距离的某个负幂增大。指数描述曲线接近临界点的方式,前面的系数则保留材料和单位的信息。

温度方向上的三种临界幂律式 8
符号含义单位
t约化温度;|t|为其绝对值1
T, Tc绝对温度与临界温度K
|m|自发序参量大小1
B序参量幅度系数1
χ此处为对无量纲共轭场的序参量响应,不是混合接触参数1
χ0响应幅度系数1
ξ, ξ0关联长度及其幅度系数m
β, γ, ν序参量、响应和关联长度的临界指数1
∼临界附近主要渐近幂律,幅度按临界两侧分别定义1

怎样理解以相对温差比较不同材料。正指数使序参量在临界点消失,负指数使响应或长度增长。幅度保留单位,指数控制倍率关系。

算一个例子若ξ0=1 nm、ν=1,|t|从0.1降到0.01,ξ从10 nm升到100 nm。若B=1、β=1/2,−t从0.04降到0.01,|m|从0.2降到0.1。

使用条件指定普适类、适当临界区间、近零共轭场与充分平衡;有限尺寸、背景和交叉效应需在数据中辨认。

幂律的直观含义是固定倍数规则。如果关联长度的指数为1,把到临界点的相对距离缩小十倍,关联长度就增大十倍;指数为一半时,只增大约3.16倍。把两条坐标都改为对数刻度,乘法变成加法,幂次成为直线斜率。因此双对数图是检查标度的工具,但有限范围内的一条近似直线还要结合背景项、误差和交叉区间分析。临界温度若估计偏了一点,横轴距离在最靠近临界点时便可能出现很大相对误差。

常用的六个静态临界指数分别对应不同观测:β描述自发序参量消失,γ描述零场响应增强,ν描述关联长度增长,δ描述恰在临界温度时序参量如何响应外场,α描述热容的奇异部分,η描述临界关联函数相对于简单距离幂律的修正。这里γ是指数,χ是响应函数,二者不能因为都与“响应”有关而混成同一个量。临界指数也不同于某一材料的临界温度:温度给出转变发生在哪里,指数给出接近它时怎样变化。

热容、临界等温线与空间关联式 9
符号含义单位
csing采用固定热容基准归一化后的奇异热容部分1
A归一化热容的幅度1
t约化温度(T−Tc)/Tc1
T, Tc绝对温度与临界温度K
m无量纲序参量1
h此式为无量纲共轭场,已用固定能量尺度归一化1
D临界等温线幅度1
C(r)序参量关联函数1
K关联幅度;此处为符号,不是开尔文单位1
r, r0间隔距离与固定参考长度m
d空间维数1
α, δ, η热容、临界等温线和关联衰减指数1
∼主导渐近关系;α=0时还需区分对数奇异性与有限跳变1

怎样理解每个指数对应一种控制条件:热容沿温度变化,临界等温线固定t=0而改变场,空间关联则固定临界条件并改变距离。长度比确保幂的底数无量纲。

算一个例子若δ=3、D=1,h从0.008降为0.001,|m|从0.2降为0.1。二维Ising的d=2、η=1/4,距离增加16倍时,关联幂律因子变为原来的一半。

使用条件静态平衡临界渐近式。热容需扣除规则背景;关联式适用于微观长度以上、有限尺寸以下的临界尺度。

这些指数不是六个可以任意调整的拟合数字。在常见的短程相互作用临界体系中,自由能的标度结构使它们满足若干关系。例如序参量、临界等温线和响应指数相互联系,关联长度与关联函数又约束响应指数;在适用超标度的范围,空间维数也进入关系。二维短程Ising模型提供一个可核对例子:β为1/8、γ为7/4、δ为15、ν为1、η为1/4,α为0但热容具有对数奇异性。α为零并不单独说明热容是常数,还要识别是否存在对数或有限跳变。

临界指数之间的标度关系式 10
符号含义单位
α热容指数1
β自发序参量指数1
γ响应指数1
δ临界等温线指数1
ν关联长度指数1
η临界空间关联指数1
d空间维数1

怎样理解这些关系把不同观测约束到同一个奇异自由能与关联标度结构。最后一式为超标度关系,额外要求相关体积的自由能标度适用。

算一个例子二维Ising:α=0、β=1/8、γ=7/4、δ=15、ν=1、η=1/4。可逐项核对0+2/8+7/4=2,(1/8)×14=7/4,1×(2−1/4)=7/4,2=2×1。

使用条件常见静态平衡标度框架;超标度通常用于上临界维数以下,不能把最后一式直接用于更高维的普通平均场指数。对数修正等需独立处理。

平均场理论的常见标量结果是β为1/2、γ为1、δ为3、ν为1/2、η为0,热容在简单双井模型中出现有限跳变。对于短程标量模型,上临界维数为4;在更低维数,临界涨落通常改变指数,在4维会出现对数修正,更高维数可恢复平均场型指数。这里维数指空间的独立方向数,不是把时间加成第四维。它也与序参量有几个分量不同:二维薄膜上的三分量磁矩,空间维数与序参量维数就是两个不同数字。

实际数据常显示交叉:离临界点较远时平均场近似较好,靠近后才进入涨落主导范围,再靠近又受样品有限尺寸或平衡时间约束。用哪一段区间拟合,因而是物理选择的一部分。先说明所拟合的是序参量、背景扣除后的响应还是关联长度,再给温度区间和几何条件,才能判断指数究竟反映临界规律、交叉行为,还是其他过程产生的有效斜率。

为什么不同材料会共享一种临界规律

水、磁体和聚合物共混物的分子并不相同,却可能在合适条件下具有相似的临界指数。原因可以从不断放大观察尺度理解。靠近临界点,关联长度远大于分子尺寸;我们可以把一个小区域里许多局部变量合成一个集体变量,再用这些集体变量描述更大的区域。许多原子尺度的细节在这个过程中改变了数值系数,却不再决定长距离变化的主要形式。长距离上仍保留的对称性、维数和相互作用性质,便成为分类的关键。

Kadanoff在1966年的块自旋思想把这个过程画得很清楚:把许多格点分成小块,每块用其集体磁化表示,再问这些块如何相互影响。块要比原子尺度大,又比关联长度小,才保留正在研究的长距离信息。Wilson随后把逐层改变观察尺度时参数怎样变化发展为重整化群方法。这里“重整化”不是把图像模糊掉便宣布细节无关,而是同时改变观察单位和有效相互作用,再追踪哪些改变会在更大尺度上放大、哪些会减弱。

1966:把局部变量组成集体变量每块用一个集体有序变量表示同时重新描述块之间的有效作用
图 3 · 虚线框把多个局部变量归成一块,下图用一个集体变量表示每块。粗粒化同时改变尺度和有效相互作用;示意中的整齐块结构便于理解操作,不代表临界构型总是规则棋盘。

如果改变尺度后,描述趋向同一种形式,就说到达或接近一个固定点。不同材料的起始参数可以不同,只要在放大尺度时走向同一固定点,便可以共享临界指数,形成同一普适类。临界温度、密度和幅度系数仍取决于化学组成与微观作用,所以普适性没有抹去材料差异。它指出的是某些指数和合适无量纲比值可以共享,而不是任意两种材料的整条相图都相同。

确定普适类时,空间维数和序参量的分量数十分重要,但还要考察对称性、相互作用是否长程、无序或约束是否改变主导行为。单分量短程标量有序常与Ising类相关,两分量平面转动可联系XY类,三分量旋转对称磁矩对应Heisenberg类。聚合物自避链与O(n)模型的联系涉及分量数趋于零的数学极限,并不是“四分量模型”;它把一类路径计数问题与临界统计联系起来,不能按空间中看见几条链来解释这个分量数。

对高分子共混物,组分涨落往往提供一个接近标量的序参量,短程相互作用和三维体相条件下可以接近三维Ising类的临界行为。长链又会改变涨落显著的温度范围,使某些实验区间看起来更接近平均场。薄膜几何、交联约束、长程电作用或明显的动力学冻结则会改变需要研究的问题。模型选择应从组成变量、几何和相互作用出发,而不是因为材料名称里有“高分子”就预先指定一种指数。

一个完整的实验判断可以按这样的逻辑推进:先用结构和组成信息确定两个状态的区别,再选择序参量;用可逆路径和充分等待寻找平衡范围,区分潜热与展宽响应;用散射或空间关联测量补上关联长度,再检验同一温度区间内的指数关系。如果温度继续逼近时数据偏离原来的幂律,就检查有限尺寸、平衡时间和交叉尺度。这样的分析把每一个偏离变成可以测量的物理问题,而不是急于把所有曲线拉到同一条直线上。

相变最有意思的地方,是宏观规律来自许多局部选择之间的联系。一个单元怎样转向或一个链段更愿意接触谁,提供微观规则;关联把这些规则传到更大的范围;自由能与序参量记录集体结果;标度则比较不同材料接近临界点的方式。把这几层连起来,便能同时理解为什么有些变化突然发生,为什么有些变化连续却异常敏感,以及为什么截然不同的材料能够共享一组长距离规律。