Collective dynamics of higher-order Vicsek model emerging from local conformity interactions

技術要約：局所的な適合性相互作用から生じる高次ヴィセック・モデルの集団力学

問題提起
ヴィセック・モデルは、自己駆動粒子系（アクティブマター）における集団運動を理解するためのパラダイメティックな枠組みとして機能している。標準的な連続時間定式化およびそのほとんどの拡張において、整列（アライメント）は厳密にペアワイズ（二体）相互作用に基づき、粒子の方向は近傍の角度の線形結合に基づいて更新される。しかし、生態学や神経科学における最近の証拠は、集団行動がしばしば、三体効果や知覚を介したルールのような、高次（非ペアワイズ）相互作用を伴うことを示唆している。「スウォーマレーター（swarmalator）」モデルは高次相互作用を含むように拡張されているが、そのような高次相互作用がどのようにしてアクティブマター系において自然に発生するのか、そのメカニズムは依然として不明である。本論文は、単純かつ物理的根拠に基づいた局所的なルールから、いかにして高次相互作用が創発し得るかという理解のギャップに対処し、それらがアクティブ粒子の集団力学に与える影響を調査するものである。

手法
著者らは、相互作用の重みが一定ではなく、近傍の局所的な適合性に依存する、連続時間重み付きDADAM（Dry, Aligning, Dilute, Active Matter）モデルを提案している。

1. モデルの導出: 核となるメカニズムは「局所的な適合性（local conformity）」ルールである。各粒子 $j$ は、近傍の粒子 $k$ が局所的な多数派のコンセンサス（局所整列場 $\psi_j$）にどの程度一致しているかに基づいて、重み $w_{jk}$ を割り当てる。この重みは、$w_{jk} = a + b r_j \cos(\psi_j - \theta_k)$ と定義される。ここで、$a$ と $b$ は定数、$r_j$ は局所的な極性秩序、$\theta$ は粒子の角度を表す。
2. 高次項の創発: これらの重みを標準的な整列方程式に代入することで、ペアワイズ項（$\alpha$ に比例）と高次（三体）項（$\beta$ に比例）の両方を含む力学系が得られる。結果として得られる角度 $\theta_j$ の運動方程式には、3つの角度の積を含む項 $\sin(2\theta_k - \theta_l - \theta_j)$ が含まれる。
3. 解析手法:
   • 数値シミュレーション: 著者らは、2D周期境界ボックス内の $N=2048$ 個の粒子を用いたオイラー・丸山法によるシミュレーションを行った。決定論的（$\eta=0$）および確率論的（$\eta > 0$）の両方のレジームにおいて、結合強度 $\alpha$（ペアワイズ）と $\beta$（三体）、およびノイズ強度 $\eta$ を変化させて解析を行った。
   • 平均場理論: 解析的な洞察を得るために、著者らはグローバル結合極限（$r_0 \to \infty$）を導出し、システムを高次相互作用を伴う確率的クラマース・モデルへと簡略化した。そして、定常解とその安定性を決定するためにフォッカー・プランク方程式を解析した。

主な貢献

• 高次相互作用の自然な創発: 本論文は、高次（三体）相互作用が単なるアドホックな追加ではなく、粒子が局所的な多数派と一致する近傍粒子を優先的に選択するという、単純な局所的適合性メカニズムから自然に創発することを実証している。
• 新規な双方向秩序相: 著者らは、高次モデルに特有の新しい定常相である「双方向秩序相（bidirectionally ordered phase）」を特定した。この状態では、粒子は互いに逆方向に進む2つの明確なグループへと自己組織化する。この相は高次相互作用によって安定化されており、標準的なペアワイズ・ヴィセック・モデルではこのような状態は不安定であるが、本モデルではランダムな初期条件から到達可能である。
• 遷移メカニズム: 本研究は、秩序・無秩序転移の性質が、ペアワイズ相互作用（$\alpha$）と三体相互作用（$\beta$）の相対的な強さに決定的に依存することを明らかにしている。

結果

• 決定論的力学 ($\eta = 0$):
  • システムは、均質無秩序、ネマティック無秩序、完全秩序、および双方向秩序の4つの異なる相を示す。
  • ペアワイズ項と高次項の競合により、「部分的秩序相」も出現する。これはノイズがない場合でも発生する。
  • 双方向秩序相は、高次相互作用が支配的な場合（$\beta > 0, \alpha \approx 0$）に安定する。
• ノイズの影響 ($\eta > 0$):
  • 連続的 vs 突発的転移: ペアワイズ相互作用が支配的な場合、秩序から無秩序への転移は連続的（スムーズ）であり、古典的なヴィセック・モデルに類似している。しかし、三体相互作用が支配的な場合、転移は突発的（不連続）となり、極性秩序パラメータ $R$ の突然の跳躍を特徴とする。
  • 双安定性: ペアワーズと高次項が競合するレジーム（例：負の $\alpha$ と正の $\beta$）では、システムは双安定性を示す。初期条件に応じて、システムは同じノイズ強度に対して、秩序状態または無秩序状態のいずれかに落ち着くことができる。
  • 相図: 著者らは、多安定領域を示す部分的な相図と、双方向秩序が持続する条件をマッピングした（これは非常に弱いノイズ下でのみ成立する）。
• 理論的検証: グローバル結合（平均場）理論は、これらの相の存在と分岐メカニズム（突発的転移におけるサドル・ノード、連続的転移におけるピッチフォーク）を定性的に捉えているが、空間的クラスタリングを考慮しない限り、双方向秩序のような空間的に不均質な相を完全に記述することはできない。

意義と主張
本論文は、局所的な適合性相互作用の導入が高次のアクティブマターの集団力学を根本的に変えることを主張している。

• 明確な臨界挙動: 著者らは、高次ヴィセック・モデルが、突発的な秩序・無秩序転移と双安定領域の存在によって、標準的なペアワイズDADAMモデルとは異なる臨界挙動を示すことを示唆している。
• 普遍性クラス: 著者らは、熱力学的極限（$N \to \infty$）においては、臨界挙動が既知のモデルと同等である可能性があると指摘しつつも、有限の集団（現実世界のシステムに関連する）においては、高次相互作用の存在が、異なる普遍性クラスを示す可能性のある明確な挙動をもたらすと主張している。これはさらなる調査を要する未解決の課題である。
• 実験との理論的関連: 本研究は、動物の集団（魚の群れやマウスの群れなど）における非ペアワイズ相互作用の観察を、アクティブマター物理学へと結びつける最小限の理論的枠組みを提供している。これは、文脈依存的な局所的相互作用が、生物学的データに見られる三体効果を効果的に生成し得ることを示唆している。

著者らは、高次相互作用を理解することが、群れ行動のより完全な理解に向けた重要なステップであると結論付けているが、完全な記述のためには、完全な流体力学的理論と有限サイズ効果に関するさらなる調査が必要である。