液体没有结晶，为什么也会变得几乎无法流动？

看起来仍然无序，为什么已经很难流动

如果把液体缓慢冷却，分子有时会排列成晶体。但如果材料和冷却条件使结晶来不及发生，另一条路径就可能出现：排列仍无长程周期，结构改变却越来越慢，最后在观察时间里表现得像固体。

这时不是所有分子都停止运动。分子仍会振动，也可能在原位置附近晃动；困难在于逃出原有邻域，完成足以松弛整体结构的重排。几次小幅运动，与成功更换周围邻居，是不同的过程。1

因此，“没有结晶”只告诉我们它没有形成那种周期秩序，并没有告诉我们改变排列需要付出多大代价。想象一个人挤在拥挤人群中：他的脚可以动，但要换到另一处位置，往往需要旁边的人也让出空间。玻璃的问题还多一层：周围环境会储存和释放应力，一处变化可能帮助别处变化，也可能受到远处边界限制。

这只是理解问题的入口。分子没有意图，液体也不是一群预先固定的人；真正的研究要把这些比喻换成可测量的松弛时间、结构相关、重排事件与响应曲线。

本篇的主问题是：局部结构、允许的变化路径和周围材料的响应，分别怎样参与慢化；改变条件后，哪些证据能区分它们？

先分清三种容易混在一起的过程

过冷液体随降温变慢、已经形成的玻璃继续老化、稳定玻璃加热后转向液态，关系密切，却不是同一个实验。

玻璃转变温度$T_g$与观察方式和时间尺度有关。若某材料在几秒内来不及明显改变，在更长窗口中可能已经发生可测松弛；但不能由此承诺任意真实玻璃在一个可实现的等待时间后都能完全流平。1

最简单的Maxwell黏弹性模型可以说明时间尺度的作用。给材料一个小的固定应变后，模型中的应力按

衰减。这里$\eta$是黏度，$G$是模型中的弹性模量。观察时间远小于$\tau$时，应力还保留着；远大于$\tau$时，大部分应力已经释放。这个单时间模型没有描述真实玻璃的全部宽谱松弛，但它说明，“能承受一阵子剪切”与“结构永远不能改变”不是同一判断。

接下来讨论的几种解释，都需要回答为什么有关的$\tau$会变得很大。

不同材料的温度依赖也不一样。一些玻璃形成液体较接近固定斜率的Arrhenius关系，另一些在降温时明显弯曲。文献中的“强”与“脆”常用来描述这种温度敏感性，不能按字面理解为玻璃是否结实、是否容易摔碎。1

一次重排，要跨过怎样的困难

一类解释从结构和可选择的配置出发。如果局部排列被周围环境强烈约束，要换成另一种排列，就需要更大区域共同调整，或者付出与周围不匹配的代价。RFOT等理论把这种关系发展成关于无定形结构、配置数和长度尺度的描述，但在真实有限维材料中，具体对应仍是研究问题。16

另一类解释从动力学规则出发。一处可以重排的区域，能够创造邻近区域继续重排的条件。随着活动变稀少，运动需要等待、接续或者绕行，整体就慢下来。这通常叫“动力学促发”：某处的运动使别处更容易运动。

还有一类解释强调弹性代价。分子移动时，周围致密材料不能完全不受影响；局部变形产生的应力，需要通过周围材料调整或释放。如果这个环境越来越难以响应，重排所面对的代价就可能增加。

这些描述并不天然互斥。 应力变化可以成为促发的物理媒介，结构也可以决定哪里容易产生重排。问题在于：各自使用的量能否测到，能否预测改变条件后的结果，以及是否还需要未测量的假设。

例如，“这个区域看起来比较松，因此容易动”，只是一种相关解释。若控制温度和初始统计结构，只改变邻近区域的运动能力，目标区域的松弛随之改变，就提供了更直接的动力学信息。

一个实际计算的例子：相同分布，可以有不同速度

为把“状态”和“改变状态的路径”分开，本篇计算了一条只有8个可变位置的二态链。每个位置可以处于0或1，1代表一种允许邻近位置较容易改变的活动状态。它不是一个真实分子，也不是实测缺陷计数。

设每个位置在平衡时处于1的概率为$c$。比较两个模型：

• 独立模型：每个位置都可以自行变化，$0\to1$的速率为$c$，$1\to0$的速率为$1-c$。
• 相邻约束模型：保留上述速率，但只有左邻处于1时，这个位置才允许变化。最左侧另设一个始终活动的边界，使有限链不会卡在完全无法变化的状态。

后一个是有限East模型的一个版本。两种规则都满足详细平衡，并具有相同的平衡分布：

这里相邻位置的平衡统计没有额外相关性。然而，当一串位置都处于0时，受约束模型必须先通过一系列允许的变化，才有机会使较远位置改变。观察到同一种状态的概率，与到达或离开它的速度，可以分开变化。

图比较总活动比例的归一化自相关：初始活动的偏差在之后还能保留多少。自相关曲线下面积，是本例采用的积分松弛时间。

对256个状态的生成矩阵实际求解后，在$c=0.1$时，独立模型的积分松弛时间为1，受约束模型约为898；在$c=0.2$时分别为1与93；在$c=0.35$时分别为1与15。时间使用模型单位，两个模型在每个$c$下的平衡分布完全相同。

这里不是通过人为给整个受约束系统乘一个统一慢速系数得到结果。两个模型使用相同的基本变化速率，差别在于哪些变化当时允许发生。因此结果具体说明：即使完整静态分布相同，改变路径仍能使动力学明显不同。

这也不是对真实玻璃理论的最终裁决。小链没有分子位置、应力场、真实配置熵和材料黏度，数值还依赖链长与边界。它只检验一个推理：仅凭静态图像相似，不能断定系统应该同样快；也不能从这个反例直接否定结构和热力学在真实材料中的作用。

在分子模拟里，怎样把邻近活动的作用单独拿出来

Herrero和Berthier的2024年研究采用了更接近粒子体系的比较：先让整个模型液体在同一温度达到平衡，再令一个空间区域使用能够加快平衡采样的交换算法，另一部分保持通常的分子动力学。两边没有被设成一热一冷。5

随后观察仍使用通常动力学的区域：靠近快速区域的地方，邻居更换和结构松弛发生得更快；距离较远时，逐渐接近该模型本来的体相行为。低温下，这种活动向外传播变得更慢。

这比看到“快区域挨着快区域”更进一步，因为研究者主动改变了邻近区域的运动规则。补充材料也检查了几何、交换算法设置和结构分布等因素。

但实验对象仍是指定相互作用与算法下的模型粒子。特别低温时，远处体相的完整松弛太慢，研究不得不使用外推来估计部分动态长度；幂律传播和对数形式在可观察范围中也都能描述部分数据。不能把有限窗口中的拟合直接说成唯一的长时定律。

还有一个值得保留的结果：补充材料中的短时笼内振动指标，其空间影响范围没有随着降温像长时松弛范围那样增长。这说明，某个短时“软硬”指标并不自动等于决定长时重排的全部弹性机制。它没有直接否定所有弹性解释，却提高了这些解释与观测对应的要求。5

2026年的软硬边界实验，改变了什么

2026年Luo等人的研究把问题推进到真实的分子玻璃薄膜。他们使用同一种有机分子TPD，在基底上先放一层由液体冷却形成的玻璃，再在其上形成一层更稳定的玻璃。这样，要观察的转化界面与基底之间，隔着约50 nm的同种分子材料。2

比较的基底包括刚性硅，以及覆盖薄橡胶层的较软边界。研究者用椭偏测量和光学分层模型，跟踪埋在材料内部的界面位置。

这个设计的价值是拉开了“直接接触”和“远处影响”。如果只有紧邻基底的分子受到作用，一个隔着缓冲层的界面不应轻易出现同样明显的变化。论文报告的速率差异能在更远距离上持续，给出了非局部影响的证据。

同时，作者做了几项有用的控制：底层玻璃在两种边界上的$T_g$接近；处理过的基底与未额外清洗的基底得到相似结果；光学拟合比较了表面和内部两处转化前沿，以及是否允许底层折射率随老化改变。相关控制与结果在方法及扩展数据中可以回查。2

它仍不是一个所有结构量都完全相同、只改变弹性参数的实验。软边界上形成的稳定玻璃，虚拟温度和分子取向也与刚性边界上的样品不同。作者指出，更高稳定性通常会使转化更慢，而这里软边界的转化反而更快，这使“仅由稳定性差异解释”受到约束；但把所有局部重排能垒认作完全相同，仍包含推论。

这里的“虚拟温度”用于描述玻璃与某个过冷液体状态的热力学对应，不是样品中另有一块区域实际处在那个温度。温度、制备历史和当前结构，不能用一个数互相替代。

界面向上走，未必全是前沿推进

这项实验有一个容易漏掉的读出问题：底层玻璃也会老化、变得更致密，从而收缩。埋藏界面的位置变化因此包含两个方向相反的过程：

前沿推进增加已转化区域的厚度，底层老化收缩则使界面位置往下移。所以，界面几乎不动，也可能是两个非零速率互相抵消；向下移动也不等于发生“反向液化”。

论文在低温、较硬边界的样品上确实观察到净位移很小或向下的情况。不能把这些条件在公开速率表中的空缺补成零，更不能把它们当作已测到一个精确的负转化速度。

本篇从原始扩展数据表提取了7个温度行，复算可直接配对的前沿速率。软边界相对于刚性边界的校正后速率比，在343 K约为2.27，在333 K为3.00，在328 K约为9.75。低温时差异扩大，但更低温度下刚性边界缺少可直接比较的正前沿速度，不能继续按同样方式给出实测比值。

对软边界的5个温度点，将对数速率对逆温度做无权重线性回归，得到表观活化能约407.6 kJ/mol，回归斜率标准误约15.0 kJ/mol，与原文的408 ± 15相符。若直接拿净界面位移率去拟合，结果则约为439.5 ± 11.5 kJ/mol，与扩展数据中的440 ± 12相符。2

差别集中体现在最低温点：318 K时，表中净位移率为$2.3\times10^{-7}$ nm/s，估计的老化收缩率为$4.8\times10^{-7}$ nm/s，相加得到的前沿速率为$7.1\times10^{-7}$ nm/s，约是净位移率的3.09倍。

因此，“多数温度点老化修正较小”可以成立，却不能把这项修正从整个低温拟合中直接删去。

这次复算使用的是论文公开表中已经汇总、舍入的数字，没有重新拟合原始椭偏信号。上面的±也只是这组回归的标准误，并未包含光学模型、老化率估计和全部实验误差的传播。

“测到能垒”究竟测到了哪一层

把温度相关的速率画成图，可以从斜率定义表观活化能。若用松弛时间表示，一种常用定义是：

它回答的是：稍微改变逆温度，松弛时间的对数改变得有多快。它不是直接拍到单个分子沿反应路径爬过了多高的能垒。

特别地，若写成$\tau=\tau_0\exp[E(T)/(RT)]$并假设前因子$\tau_0$不随温度改变，那么：

当$E$本身随温度改变时，斜率得到的量与指数里的$E(T)$就不是同一个量。前因子若也变化，还会增加一项。

Luo等人用软边界的近似Arrhenius斜率估计局部贡献，再把刚性边界或体相的表观能垒与它相减，估计非局部弹性贡献。这个分解需要认为：软边界充分释放了非局部代价，而且不同条件下的局部部分足够可比。2

本篇也按论文印出的体相VFT拟合复算了斜率：

这里采用原文的$B=3493$ K、$T_0=258.9$ K和$\log_{10}(\tau_0/\mathrm{s})=-19.2$，并将其与软边界回归斜率比较。

按这组假设，两部分在约319.1 K处相等，能够复现论文约319 K的交点。然而，复现一个曲线交点，并没有单独证明该交点就是材料产生玻璃化的原因。

还有两点限制。第一，刚性边界前沿速度与体相松弛速率的对应，并非在所有条件下独立确定的普适关系；论文扩展数据中的一条比较线就固定了斜率。第二，向更低温度推算出很大的速度差，是外推结果，不能与实验窗口中的直接比较混写。

因此，这项工作较有力地支持：远处边界能够改变这里测到的转化动力学，弹性释放是一个有具体对照支持的解释。 至于从两条温度曲线分解出唯一的局部与弹性能垒，以及据此解释所有玻璃形成体，证据还需要继续连接。

看到液态区域长大，就一定是普通成核吗

2023年的显微研究在高稳定TPD玻璃上观察到了局部软化区域和向外扩展的形貌变化。其方法很有意思：玻璃被更难软化的材料夹住，局部转向更易变形的状态后，约束与应力使表面出现可追踪的褶皱。3

研究结合原子力显微镜、光学观察与力学模拟；补充材料还用衍射检查了样品保持无定形，并比较了去掉一侧包覆层的情况。因此，这些图像有明确的物理依据，但依然是通过形貌和力学响应推断内部转化，不是逐个分子运动的实时录像。区域半径的误差、温度下的几何校准和有限成像时间，也限定了能分辨的细节。

看到一个区域出现、长大，并形成转化曲线，很容易联想到晶体的成核与生长。但这还不是唯一解释：已经活动的区域也可能不断使周边更容易活动，从而形成相似的前沿。

Chacko等人的2024年研究提出了更有区分力的比较：对同一个模型玻璃，不只研究从低温初态加热，也研究反方向的冷却，并与真正发生结晶的模型对照。它们在玻璃加热时看到区域长大，在冷却时却没有看到对应的区域生长；一些只有动力学约束的模型也产生类似不对称。4

这对“同一对状态之间就是传统界面自由能控制的成核过程”构成了限制。但论文讨论也明确指出，这不排除有限的无定形结构相关长度或RFOT式结构图景。前沿粗糙度在几个模型中的表现还不同，不能把某一种粗糙度增长律当成所有促发机制的固定签名。

这组比较的价值在于：从一条漂亮的转化曲线，前进到正反方向与替代模型的共同检验。

目前能怎样把这些机制放在一起

本篇没有得到一把能从材料名字直接算出全部松弛时间的钥匙。但可以把问题拆得更清楚。

结构决定哪些排列和局部环境能够存在；变化规则决定它们之间怎样转换；周围材料和边界又会改变转换的代价与传播。这些层面彼此联系，任何一个被拿来解释具体材料时，都要通过相应读数和条件变化检验。

有限East链说明静态分布并不足以决定动力学；它不具有真实玻璃的全部热力学性质，不能据此宣布结构解释已经失效。RFOT讨论也指出，在有限维体系里，某一组状态能否被视为长期亚稳，本身就可能依赖允许的动力学和观察时间。6

粒子模拟进一步支持邻近活动影响松弛，但低温长度和传播律仍有测量窗口与外推问题。分子薄膜实验则让远处边界的作用变得可测，同时提醒我们，制备历史、界面运动与局部能垒的可比性需要一起核查。

接下来最有区分力的工作，是在更多材料与边界条件中检查：影响能传播多远，哪一个弹性或动力学读出能预测这种距离，局部结构匹配后速率差是否仍存在，以及同一机制能否同时解释平衡松弛、老化和正反方向的转化。

玻璃难以改变，不要求它先长成晶体。 真正需要解释的是一次次重排怎样被环境限制、彼此接续，又如何在观察时间里累积成流动或保持。当前证据已经能够让这些问题接受具体比较，统一解释仍有尚未连接的环节。

可回查：原表数据、VFT斜率计算、有限链自相关、复算结果与边界。