[论文解读] Dynamic transition in an atomic glass former: a molecular dynamics evidence

Estelle Pitard, Vivien Lecomte+1 more|Springer Link (Chiba Institute of Technology)|May 12, 2011

Glass properties and applications被引用 7

一句话总结

本研究通过大偏差统计方法分析轨迹活性，利用分子动力学模拟为Lennard-Jones原子玻璃形成体（Kob-Andersen混合物）提供了首次一阶动态相变的证据。结果识别出临界点 s_c(∞) ≈ 0，表明该动态相变在实验中可被探测，并证实了动力学受限模型理论在原子系统中的适用性。

ABSTRACT

We find that a Lennard-Jones mixture displays a dynamic phase transition between an active regime and an inactive one. By means of molecular dynamics simulations and of a finite-size study, we show that the space time dynamics in the supercooled regime coincides with a dynamic first order transition point.

研究动机与目标

• 通过有限尺寸标度方法，确立原子玻璃形成体中动态相变的存在性。
• 精确确定热力学极限下动态相变点的位置，解决其是否处于实验可探测范围的问题。
• 基于有效力定义一个物理解释清晰的活性可观测量，避免经验性假设。
• 确认动力学受限模型理论中的动态相变情景在真实原子系统中的适用性。

提出的方法

• 采用克隆算法的分子动力学版本，采样活性出现大偏差的罕见轨迹。
• 将活性定义为时间平均的有效力，将其与粒子间的物理作用力关联。
• 使用累积量生成函数 ψ(s) 分析不同系统尺寸下轨迹活性的统计特性。
• 通过样条插值与二阶导数最大化的有限尺寸标度方法，定位尺寸为 N 时的相变点 s_c(N)。
• 分析自由能 ψ(s) 的标度行为：当 s<0 时 ∼N，当 s>0 时 ∼N^{1+α}（α≈0.43），表明存在不同的动力学相。
• 通过重标度与有限尺寸外推，推断无限尺寸极限下的相变点 s_c(∞)。

Figure 2: The dynamical free energy per particle $\frac{\psi(s,N)}{N}$ is shown as a function of $s$ for increasing sizes $N=45$ , $82$ , $155$ , $250$ and $393$ . While the curves collapse onto a single master curve at $s<0$ , a strong size-dependent behavior is observed at $s>0$ . On the $s>0$ sid

实验结果

研究问题

• RQ1 在如Kob-Andersen Lennard-Jones混合物这样的原子玻璃形成体中，是否存在动态相变？
• RQ2 该动态相变的临界点在热力学极限下的位置在哪里？是否可在实验中探测？
• RQ3 能否基于物理作用力而非经验性构造，定义动力学活性的概念？
• RQ4 自由能随系统尺寸的标度行为如何区分活跃与非活跃动力学相？
• RQ5 在大偏差系综中，非活跃构型的几何与拓扑特性为何？

主要发现

• 系统在活跃与非活跃态之间表现出一阶动态相变，由活性分布的双模态特征及自由能 ψ(s) 的拐点所证实。
• 相变点 s_c(N) 随系统尺寸 N 增大而减小，且在无限尺寸极限下趋于接近零的值。
• 自由能 ψ(s) 在 s<0 时标度为 ∼N（活跃相），在 s>0 时标度为 ∼N^{1+α}（α≈0.43，非活跃相），表明非活跃构型具有非平凡几何结构。
• 非整数指数 α≈0.43 暗示非活跃相中存在长程有效相互作用及复杂的空间结构。
• 相变点位于 s_c(∞)≈0，使其在探测轨迹统计的实验中具有潜在可观测性。
• 通过有效力定义的活性可观测量，提供了物理解释清晰且可测量的动力学异质性表征。

Figure 3: The average activity per particle $\langle V_{\text{eff}}\rangle/N$ is measured in an ensemble of trajectories biased by a weight $\text{e}^{-sK}$ and is plotted as a function of $s$ . While the $s<0$ regime shows convergence to a size independent limit of an order comparable to its equili