🔬 condensed matter

Glassy dynamics, crossover temperature and density scaling in fragile glass-formers

本研究结合理论与大规模分子动力学模拟，证明了反幂律系统的热力学与玻璃动力学可以通过特定的交叉温度和密度-温度定标进行统一，从而实现通过单个状态点预测广泛密度范围内的弛豫行为。

原作者： Ankit Singh, Vinay Vaibhav, Swarn Lata Singh, Yashwant Singh

发布于 2026-09-07

📖 1 分钟阅读☕ 轻松阅读

原作者： Ankit Singh, Vinay Vaibhav, Swarn Lata Singh, Yashwant Singh

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

当液体冷却时，它通常会冻结成固体晶体，将其原子排列成整齐、重复的网格。但某些液体，例如瓶中蜂蜜里的糖浆或窗户玻璃中的二氧化硅，却拒绝结晶。相反，随着它们变冷，其分子速度减慢，直到它们实际上被固定在原地，却仍然保持着一种无序的、类似液体的排列方式。这种被称为“玻璃化转变”的现象是物理学中最令人费解的现象之一。不同于水变成冰（这发生在特定温度且伴随明显的结构变化），玻璃的形成是一个逐渐减慢的过程，材料在没有明确结构预警信号的情况下变得极其粘稠。科学家们长期以来一直致力于理解为什么会发生这种情况，分子的运动速度如何与温度相关联，以及是否存在一条隐藏的规则，支配着这些材料在接近这种冻结状态时的行为。

在这项研究中，研究人员试图通过一种特定类型的模型液体来揭示这种减速过程背谱后的隐藏机制。他们专注于一种由两类粒子组成的混合物，这些粒子通过一种随着距离缩短而增强的力量进行相互作用，这种关系被称为反幂律势（inverse-power-law potential）。通过运行涉及一万个粒子的大规模计算机模拟，团队追踪了这些粒子在不同温度和密度下的运动及排列方式。他们测量了粒子重新排列位置所需的时间（称为结构弛豫时间），并将这些观测结果与一个旨在预测这种行为的新理论框架进行了对比。其目标是找到一种统一的方法来描述这些玻璃形成液体的能量和运动，观察是否可以用一套统一的规则来解释它们在广泛条件下的行为。

研究人员发现，理解这种减速的关键在于每个粒子周围形成的特定局部有序性。随着液体冷却，粒子并不仅仅是随机减速；它们被周围邻居形成的“笼子”所困。该理论识别出一个微小的、稳定的粒子簇，作者称之为“协同重排簇”（cooperatively rearranging cluster），它构成了运动的基本单元。对于液体而言，要实现流动或弛豫，整个簇必须同时进行重组。研究发现，随着温度下降，这个稳定簇中的粒子数量会增加，且打破这些键所需的能量也会剧烈上升。这种不断增长的能量壁垒正是导致材料发生剧烈减速、从而从流动的液体转变为坚硬玻璃的原因。

这项工作的关键发现是确定了一个特定的交叉温度，它标志着从高温行为向缓慢的玻璃动力学转变的转折点。在此温度之上，粒子移动相对自由，系统表现得具有可预测性。在此温度之下，一种新的动力学占据了主导，系统中能量的波动开始起作用，导致稳定簇的形成以及弛豫时间的激增。研究人员定义了一个参数，用于衡量这些波动如何影响粒子簇的稳定性。他们发现，在高温度下，该参数保持不变，但一旦系统冷却到交叉点以下，该参数就会大幅下降。这种下降标志着玻璃态的开启，并为材料行为发生根本性变化的时刻提供了一个清晰、可测量的标记。

研究还揭示了一个统一这些液体行为的强大定标律（scaling law）。研究人员发现，如果将密度和温度合并为一个单一变量，来自不同条件下的所有数据都会坍缩到一条单一的主曲线之上。这意味着，在高密度和低温度下的复杂液体行为，通过这个组合视角来看，在数学上与其在低密度和高温度下的行为是完全一致的。这一定标律对热力学性质（如内能和压力）以及动力学性质（如材料弛豫所需的时间）均成立。这表明，粒子间的排斥力是这些系统中玻璃化转变的主要驱动力，超越了相互作用的其他细节。

为了验证他们的理论，团队将预测结果直接与计算机模拟结果进行了对比。他们根据稳定簇的大小以及打破它们所需的能量，计算了粒子重新排列所需的时间。理论预测值与他们在各种密度和温度下测试的模拟结果高度吻合，精度极高。他们还验证了理论预测的交叉温度与模拟数据表现出明显行为变化的温度高度一致。这种理论与模拟之间的一致性提供了强有力的证据，证明这些稳定、协同簇的形成确实是驱动脆弱玻璃形成者（fragile glass-formers）剧烈减速的机制。

研究人员进一步探索了这些关键簇的大小。通过分析数据，他们估计在材料实际上变成玻璃的点，稳定簇仅包含大约十四个粒子。这是一个令人惊讶的微小数字，主要局限于中心粒子的直接邻居。然而，这些稳定粒子被嵌入到一个更大、更松散的粒子群中，后者可以更自由地移动。研究表明，虽然这个微小的稳定簇是控制液体整体速度的瓶颈，但较大的粒子群会对这个小核心的重组做出响应。这种“两尺度”图景——即一个小而坚硬的核心被一个较大的移动云团所包围——为玻璃的形成提供了一个清晰的物理图像。

最终，这项工作为玻璃形成液体如何表现提供了一个统一的描述。它将液体的静态结构、运动所需的能量以及弛豫所需的时间连接到了一个单一且连贯的框架之中。通过确定交叉温度和结合了密度与温度的定标变量，研究人员证明了玻璃形成的复杂动力学可以从单一状态点收集的数据中进行预测。这意味着，如果科学家了解一种材料在某一组特定条件下的行为，他们就能准确预测其在许多其他条件下的行为。这项研究不仅描述了发生了什么，还解释了为什么会发生，指出小规模粒子的协同重排是这些系统中玻璃化转变的根本原因。

技术摘要：玻璃态动力学、交叉温度与脆弱玻璃形成物的密度定标

问题陈述
随着系统接近玻璃转变，动力学发生的剧烈减慢仍是凝聚态物理学中的一个核心挑战。与传统的临界现象不同，玻璃转变缺乏一个普遍接受的热力学奇点或明确的发散相关长度。尽管先前的理论框架试图将动力学减慢与热力学量（例如 Adam-Gibbs 理论）或特定的结构基元联系起来，但建立一个能够准确预测不同密度和温度下各种脆弱玻璃形成物弛豫时间的统一描述仍然难以实现。具体而言，理解分子间相互作用如何决定从高温行为到低温行为的交叉，以及这些动力学如何随密度和温度进行定标，需要进一步的研究。

方法论
作者研究了一种通过反幂律（IPL）势能 $u_{\alpha\gamma}(r) = \epsilon_{\alpha\gamma}(\sigma_{\alpha\gamma}/r)^n$ 进行相互作用的 50:50 二元粒子混合物。该研究结合了大规模分子动力学（MD）模拟与作者此前开发的理论框架，采用双重方法进行研究。

模拟： 使用 LAMMPS 软件包对 $N=10,000$ 个粒子进行了 MD 模拟，涵盖三种密度（ $\rho = 0.75, 0.80, 0.85$ ）及一系列温度。系统在 NVT 系综下通过 Nosé–Hoover 热浴进行平衡。测量的关键量包括：
- 静态对相关函数 ( $g_{\alpha\gamma}(r)$ )。
- 通过热力学积分和共轭梯度最小化得到的构型熵 ( $S_{conf}$ ) 和固有结构能量 ( $e_{IS}$ )。
- 结构弛豫时间 ( $\tau_\alpha$ )，通过自相关中间散射函数 $F_s(q,t)$ 和自重叠相关函数的衰减来确定。
- 通过分析 $F_s(q,t)$ 对 $\log t$ 的导数确定的动态交叉温度 ( $T_h$ )。
理论： 作者利用一种识别“协同重组簇”（CRC）的理论。该理论计算了围绕中心粒子形成稳定簇的粒子数量，这些粒子由能够抵抗热涨落的“稳定键”（s-粒子）定义。理论参数 $\psi(T)$ 量化了浴涨落对稳定 CRC 大小的影响。该理论假设，当 CRC 发生不可逆重组时，结构弛豫就会发生，其激活能 $\beta E^{(s)}$ 由打破这些稳定键所需的能量决定。

主要贡献与结果

识别交叉温度： 研究识别并比较了三种不同的交叉温度：
- $T_{IS}$ ：源于固有结构能量偏离高温行为。
- $T_h$ ：源于 $F_s(q,t)$ 的时间导数出现双极小值结构。
- $T_a$ ：一个理论预测的交叉温度，在该温度下，涨落参数 $\psi(T)$ 从 1（高温）转变为递减函数（低温）。
  结果显示， $T_{IS}$ 和 $T_a$ 高度一致，而 $T_h$ 高出约 12%，这表明动力学性质对涨落变化的敏感度高于热力学性质。
验证 CRC 理论： 该理论框架成功预测了结构弛豫时间 $\tau_\alpha$ 。理论识别了一种局部结构有序性（即 CRC），它定义了弛豫的协同性。计算出的激活能 $\beta E^{(s)}$ 遵循关系式 $\beta E^{(s)} \propto [1 + n^{(s)}]^\delta$ ，其中 $n^{(s)}$ 是 CRC 中稳定粒子的数量。理论对 $\tau_\alpha$ 的预测与 MD 模拟结果在所有密度下均表现出极佳的一致性。
密度-温度定标： 该研究的一个重要发现是 IPL 系统的定标行为。研究表明，当绘制为复合变量 $\Gamma$ （针对 IPL 势能的密度和温度函数）或约化温度 $T/T_a$ 时，热力学量（内能 $U$ 、比热 $C_V$ 、约化压强 $Z$ ）和动力学量（ $\tau_\alpha$ 、 $\beta E^{(s)}$ ）都会塌缩到主曲线（master curves）上。
- 热力学量在作为 $1/\Gamma$ 的函数绘制时发生塌缩。
- 动力学量和参数 $\psi(T)$ 在作为 $T/T_a$ 或 $1/\Gamma$ 的函数绘制时发生塌缩。
- 交叉温度 $T_a$ 随密度按 $T_a \propto \rho^4$ 进行定标，这与 IPL 指数 $n=12$ （其中 $T_a \propto \rho^{n/3}$ ）一致。
减慢的本质： 分析表明，剧烈的减慢是由 CRC 的增长驱动的。在玻璃转变处（定义为 $\tau_\alpha \sim 10^{14}$ ），CRC 中的粒子数被发现很少（约 14 个），局限于前两个配位壳层内。然而，这个 CRC 嵌入在一个由数百个松散结合的“亚稳态”（m）粒子组成的更大簇中。小规模 CRC 的重组会触发整个更大簇的重组，从而导致观察到的宏观弛豫。

意义与主张
本文声称为 IPL 系统提供了一个统一的热力学与动力学描述。通过建立静态对相关函数、涨落诱导参数 $\psi(T)$ 与结构弛豫时间之间的直接联系，该理论提供了一个预测框架。其主要意义在于，利用观察到的密度-温度定标特性，能够通过来自单个状态点的数据来预测广泛密度范围内的弛豫行为。作者得出结论，由静态相关性和系统涨落决定的 CRC 是脆弱玻璃形成物发生剧烈减慢的主要驱动力，为观察到的热力学和动力学性质的交叉行为提供了机理性的解释。

类似论文