[論文レビュー] 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型の配向を伴う自己推進が、イベント駆動分子動力学を用いた密度の高い硬円系の相挙動にどう影響するかを研究し、極性転移と方向性転移の変位を明らかにし、それらを局所構造と遊離体積幾何に結びつける。

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型整列が排除体積硬円物理とどのように相互作用するかを調べる。
• 極性秩序の秩序喪失と六方 Orientational（Alder）転移という競合する秩序パラメータを特徴づける。
• 高充填率における剛性不変性と自己推進が固- fluid転移へ与える影響を理解する。
• 局所構造と遊離体積幾何を通じて異常な転移点シフトの微視的起源を分析する。

提案手法

• 硬円とVicsek型整列を持つ活性物質モデルを、Δt*で固定間隔更新するイベント駆動分子動力学（EDMD）枠組みで用いる。
• 排除体積は硬円衝突によって取り込み；速度方向は近傍速度のベクトル和とノイズで更新され（式1）、速度スカラーはそれに応じて更新される（式2）。
• 二つの秩序パラメータを探索する：極性秩序 ψ = |sum v_i| / sum|v_i| と六方 Orientational秩序 Φ6^G および Φ6^L（式3–6）。
• 局所構造は遊離体積 v_f^i、遊離表面 s_f^i、局所円性 c_i^*、全体的な円形性 C^*（式7–8）で分析する。
• ν（充填率）とノイズηの下でΔ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型整列は、2Dの硬円排除体積効果とどのように相図を修正するか？
• RQ2 dense硬円系における自己推進頻度（Δt*）が極性・方向性秩序に与える影響は？
• RQ3 局所構造特徴（遊離体積幾何・円形性）は転移点の異常シフトとどう関係するか？
• RQ4 dense硬円の不可Compressibilityは motility-induced相分離を抑制し、凍結/融解（Alder）転移へどう影響するか？
• RQ5 この活性硬円系における局所構造と全体的な相挙動を結ぶ微視的機構は何か？

主な発見

• Vicsek相互作用を伴う密度の高い硬円系では η ≈ 1.2π の周辺で極性秩序の喪失・発生転移が生じる。
• Vicsek相互作用の頻度を高める（Δt*が小さくなる）と、弾性的衝突との競合により転移点が移動する。
• 全体的な六方向き秩序 Φ6^G は η ≈ 1.2π の接近で鋭い谷（カスプ）を示し、ηが小さくなると ν が高くなる方向へ転移点が移動することを示し、ノイズ・Vicsek更新・結晶化の結合を示唆する。
• 局所的な向き秩序 Φ6^L は特定のηでΦ6^Gほど影響を受けにくく、結晶粒境界を形成する多結晶クラスターを伴う。
• 遊離体積幾何は動力学に大きく影響する。二つの異なる円形性範囲（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