[论文解读] Liquid, Glass and Crystal in Two-dimensional Hard disks

Ludger Santen, Werner Krauth|arXiv (Cornell University)|Jul 23, 2001

Material Dynamics and Properties参考文献 1被引用 9

一句话总结

本研究利用先进的蒙特卡罗方法，研究了二维多分散硬盘的热力学与动力学相行为。结果表明，动态液-玻璃转变发生在密度 ρᴳ 处，该密度与多分散度无关，且无伴随的热力学转变；而由多分散度驱动的晶体-无序相变是连续的（Kosterlitz-Thouless 类型），挑战了以往认为其为一级相变的假设。

ABSTRACT

We study the thermodynamic and dynamic phase transitions in two-dimensional polydisperse hard disks using Monte Carlo methods. A conventional local Monte Carlo algorithm allows us to observe a dynamic liquid-glass transition at a density, which depends very little on the degree of polydispersity. We furthermore apply Monte Carlo methods which sample the Boltzmann equilibrium distribution at any value of the density and polydispersity, and remain ergodic even far within the glass. We find that the dynamical transition is not accompanied by a thermodynamic transition in this two-dimensional system so that the glass is thermodynamically identical to the liquid. Moreover, we scrutinize the polydispersity-driven transition from the crystal into the disordered phase (liquid or glass). Our results indicate the presence of a continuous (Kosterlitz-Thouless type) transition upon increase of the polydispersity.

研究动机与目标

• 解决长期以来关于二维硬盘动态玻璃转变是否伴随热力学相变的争议。
• 研究在二维硬盘中随多分散度增加而引起的晶体-无序相变的本质。
• 开发并应用非局部蒙特卡罗算法，即使在玻璃相深处也能保持遍历性，从而获取平衡热力学性质。
• 精确确定不同多分散度下动态玻璃转变的位置与特征。
• 研究终端多分散度在抑制结晶并实现平衡模拟中玻璃相访问中的作用。

提出的方法

• 采用非局部蒙特卡罗算法，即使在玻璃态区域也能实现玻尔兹曼分布的遍历采样，从而精确计算平衡热力学可观测量。
• 使用固定概率增量模型根据均匀累积分布分配粒子半径，确保系统中具有明确定义的多分散度 ε。
• 在 [NV] 系综中结合局部蒙特卡罗动力学，计算有效扩散系数 Dᵢ 作为粒子尺寸和密度的函数。
• 利用标度形式 Dᵢ(ρ) ∼ (ρᴳᵢ − ρ)^α 将扩散系数 Dᵢ(ρ) 外推至零，以定义玻璃转变密度 ρᴳ。
• 测量键取向序参量 ⟨Ψ₆⟩ 并分析其概率分布，以检测一级相变或连续相变的特征。
• 在一系列密度 ρ 和多分散度 ε 下进行系统性模拟，包括终端多分散度极限 ε = 0.58，以绘制完整的相图。

实验结果

研究问题

• RQ1 二维多分散硬盘的动态液-玻璃转变是否伴随热力学相变？
• RQ2 玻璃转变密度 ρᴳ 如何依赖于多分散度 ε？
• RQ3 随着多分散度增加，从晶体相到无序相（液相或玻璃相）的相变本质是什么？
• RQ4 系统是否在序参量中表现出滞后现象或双峰分布，表明其为一级相变？
• RQ5 能否利用非局部蒙特卡罗方法在玻璃相深处可靠计算平衡热力学性质，包括状态方程和压缩因子？

主要发现

• 动态液-玻璃转变发生在密度 ρᴳ 处，该密度几乎与多分散度 ε 无关，且 ρᴳ 显著高于单分散系统结晶密度。
• 动态玻璃转变不伴随热力学转变；玻璃相在热力学上与液相完全相同，表明无自由能不连续性，仅存在遍历性破缺。
• 由多分散度增加驱动的晶体-无序相变是连续的，体积及键取向序参量 ⟨Ψ₆⟩ 的概率分布均无滞后现象。
• 平均键取向序参量 ⟨Ψ₆⟩ 随 ε 增加而平滑减小，其分布保持单峰，支持连续相变而非一级相变。
• 对于小 ε，外推得到的玻璃转变密度 ρᴳ 在不同粒子尺寸下具有鲁棒性，但当 ε 极大时，由于有限尺寸效应及粒子通过空隙逃逸，出现轻微偏差。
• 确认终端多分散度极限 ε = 0.58 为固定概率增量模型中的最大可实现值，超过此值系统将无法维持晶体有序。