[论文解读] Phase Transition of Hard Disk Systems with Vicsek-type Interactions

Nobuaki Murase, Masaharu Isobe|arXiv (Cornell University)|Mar 18, 2026

Material Dynamics and Properties被引用 0

一句话总结

本文研究自推进在 Vicsek-type 对齐下如何影响密集硬盘系统的相行为，利用事件驱动分子动力学揭示极化与定向转变的位移，并将其与局部结构及自由体积几何联系起来。

ABSTRACT

The phase diagram of self-propelled hard disk systems with Vicsek-type alignment interactions was investigated by event-driven molecular dynamics simulations. The model incorporates two competing order parameters: the polar order-disorder transition associated with collective velocity alignment (Vicsek model) and the orientational order arising from solid-fluid transitions (Alder transition) induced by excluded volume effects. The incompressibility of hard disks suppresses motility-induced phase separation at high packing fractions. Distinctive fluctuations were observed near the transition point, accompanied by anomalous shifts in the transition point as functions of noise intensity and packing fraction. Analysis of local configurational parameters -- specifically, orientational order and circularity of free volume -- provides insight into the microscopic origins of these anomalous phase transition shifts.

研究动机与目标

• 研究 Vicsek-type 对齐如何与两维排斔体体积硬盘物理相互作用
• 表征竞争的序参数：极化序/无序与六重定向（Alder）转变
• 理解高装填分数下不可压缩性与自推进如何影响固-流相变
• 通过局部结构与自由体积几何分析异常转变位移的微观起源

提出的方法

• 在一个具有硬盘和 Vicsek-type 对齐的活性物质模型中，在固定时间步 Δt* 更新，采用事件驱动分子动力学（EDMD）框架。
• 通过硬盘碰撞引入排斔体积；速度方向由邻近速度的向量和（带噪声）更新（式1），并相应更新速度（式2）。
• 探索两个序参数：极化序ψ = |sum v_i| / sum|v_i|与六方定向序Φ6^G与Φ6^L（式3–6）。
• 通过自由体积 v_f^i、自由表面 s_f^i、局部圆度 c_i^* 与全局圆度 C^*（式7–8）分析局部结构。
• 在不同 packing fraction ν 与噪声 η 下进行模拟，改变 Δt* 以研究碰撞与 Vicsek 更新之间的竞争。
• 通过追踪 ψ、Φ6^G、Φ6^L 与局部结构量来识别相行为及转变点。

Figure 1: (Color online) The polar order parameter $\psi$ as a function of noise $\eta$ in densely packed hard disk systems with $\nu=0.72$ , for various values of the Vicsek interaction time interval $\Delta t^{*}$ . The inset shows the packing fraction $\nu$ dependence of $\psi$ at a fixed $\Delta

实验结果

研究问题

• RQ1 Vicsek-type 对齐如何与硬盘排斔效应在二维中修改相图？
• RQ2 自推进频率（Δt*）对密集硬盘系统中的极化与定向有序有什么影响？
• RQ3 局部结构特征（如自由体积几何与圆度）如何与转变点的异常位移相关？
• RQ4 密集硬盘的不可压缩性是否抑制运动诱导的相分离，这又如何影响凝固/熔化（Alder）转变？
• RQ5 在该活性硬盘系统中，局部结构与全局相行为之间的微观机制连接是什么？

主要发现

• 在与 Vicsek 相互作用耦合的密集硬盘系统中，极化有序-无序转变在 η ≈ 1.2π 附近发生。
• 增加 Vicsek 作用频率（减小 Δt*）由于与弹性碰撞的竞争而将转变点向不同方向移动。
• 全局六方定向序Φ6^G 在 η ≈ 1.2π 附近呈现峰值，并随着 η 的减小向更高 ν 移动，表明噪声、Vicsek 更新与晶化之间存在耦合。
• 局部定向序Φ6^L 在某些 η 下相对不那么受影响，出现多晶簇聚形成晶粒边界。
• 自由体积几何显著影响动力学：存在两段明显的圆度区间（c_i^* ≈ 1.4 与 ≈ 2.3）与矩形自由体积形状相关，促使大位移跳跃与液化。
• 自由体积的形状不仅仅是大小，决定了运动性与液化，从几何角度揭示了转变位移的机制。

Figure 2: (Color online) Contour map of the global orientational order parameter $\Phi_{6}^{\rm G}$ in the $\nu-\eta$ plane at $(N,\Delta t^{*})=(4096,10)$ . The snapshots color-coded by orientation direction $arg(\phi_{6}^{i})$ , at $\nu=0.72$ for $\eta=0.1,1.2\pi$ , and $1.9\pi$ (from left to righ