Correlated and anti-correlated density dependent motility

技術要約：相関および反相関密度依存性運動性

問題提起
本研究は、生物学的なクオラムセンシングに類似したメカニズムである、密度依存性運動性（Density Dependent Motility: DDM）に支配される粒子系の創発的な定常状態と相転移を調査するものである。これまでの研究では、硬い球体モデルを用いて、特に「反相関」領域（高密度領域で運動が抑制される領域）におけるDDMを探索してきたが、「相関」された運動性（高密度領域で運動が促進される領域）を持つ系や、硬いコアではなく軟らかな斥力ポテンシャルを用いて相互作用する系に関する理解には依然として空白が存在する。本論文は、軟らかな斥力相互作用を持つ二次元粒子系において、相関および反相関の両方のDDMルール下でどのような相が出現するかという新しい問いに取り組む。

手法
著者は、Fortranによるカスタムコードを用いたランジュバン動力学シミュレーションを用い、$N=2025$ 個の粒子を、周期境界条件を持つ二次元正方形ボックス（$L=82.16$、全密度 $\rho=0.3$）内でモデル化している。

• 相互作用: 粒子は、短距離の軟らかな斥力ポテンシャル $V(r) = \epsilon r^{-6}$（カットオフ距離 $r_c=3$）を介して相互作用する。この選択は、観察されるクラスター形成が引力ではなく、運動性のルールによって駆動されることを保証するためである。
• DDMルール: 粒子は、局所密度 $\rho_i$（感知半径 $2a$ 内で計算、ここで $a$ は平均粒子間距離）に基づき、活性（active）状態と受動（passive）状態の間を交互に切り替わる。
  • 相関（Correlated）ケース: 粒子は希薄な領域では受動的であり、高密度な領域では活性的である（$\rho_i > \rho_c$）。
  • 反相関（Anti-correlated）ケース: 粒子は希薄な領域では活性的であり、高密度な領域では受動的である。
• 動力学: 粒子の運動は、自己推進速度 $v_p$、回転拡散、および熱雑音を組み込んだランジュバン方程式に従う。系は、自己推進と拡散の関係を示すペクレ数（$Pe$）によって特徴付けられる。
• パラメータ: シミュレーションは、臨界密度 $\rho_c \in [0.15, 0.45]$ と活性パラメータ $b$（活性粒子において $Pe=b$）の位相空間をスキャンする。ここで$b$は$[5, 50]$の範囲にある。系は定常状態に達するまで最大$10^6$ ステップタイムステップ進化させる。

主要な貢献と結果

1. 相関（Correlated）ケース
この領域では、高密度領域が活性化し、低密度領域は受動的なままとなる。

• 相挙動: 中間的な臨界密度において、系は流体-流体相分離を示し、これは活性粒子からなる高密度クラスターと、希薄な受動的流体の共存によって特徴付けられる。
• 転移: 2つの主要な転移が特定されている：
  • $t_1$: 一様な「活性流体」から「流体-流体」状態への転移。これは $\rho_c \approx 0.9$ の活性分率で発生する。
  • $t_2$: $\rho_c$ がさらに増加するにつれて、「流体-流体」から一様な「受動的流体」への転移。
• 活性の影響: 活性（$b$）を増加させると、両方の転移がより低い臨界密度へとシフトする。第一の領域では、活性は相分離を促進する（運動誘起相分離（MIPS）と同様）。第二の領域では、高い活性は粒子を受動的になる希薄な領域へと駆動する。
• クラスターの形態: 活性クラスターは準安定であり、時間とともに構成メンバーや形状を変化させる。活性が増加するにつれ、これらの活性クラス的クラスターのサイズは減少する。全密度 $\rho=0.3$ または $\rho=0.5$ において、六角形秩序を持つ活性流体相（hexatic active fluid phase）は観察されなかった。

2. 反相関（Anti-correlated）ケース
この領域では、希薄な領域が活性的であり、高密度な領域は受動的である。これは、人混みの中で運動が抑制される社会性昆虫のようなシステムを模倣している。

• 相挙動: 系は通常、コンパクトな受動相と活性ガスとの共存を示す。
• 転移: $\rho_c$ の増加に伴い、3つの主要な転移が特定される：
  • $t_1$: 「受動的流体」 $\to$ 「受動的ガラス + 活性ガス」（$\rho_c \approx 0.2$ 付近）。
  • $t_2$: 「受動的ガラス + 活性ガス」 $\to$ 「受動的結晶（六角形秩序） + 活性ガス」（$\rho_c \approx 0.35$ 付近）。
  • $t_3$: 「受動的結晶 + 活性ガス」 $\to$ 「活性流体」（活性分率 $\phi_a > 0.9$ のとき）。
• 構造の詳細: 受動的な固体相は、低い活性では無定形（ガラス状）であり、高い活性では六角形秩序を持つ。高い活性は受動相に圧力を加え、六角形秩序を誘導する。
• 活性構造: 特定の条件下では、活性ガスが円形のクラスター（「活性円」）を形成し、それらは時間とともに単一の大きな円へと合体する。

意義と主張
本論文は、相関および反相関のDDM間の根本的な違いを実証している。相関ケースは流体-流体分離と準安定な活性クラスターを好むのに対し、反相関ケースは、コンパクトな受動固体（六角形秩序相へと秩序化し得る）と活性ガスの共存を好む。

著者は、これらの知見が、生物学的システム（例：細胞質流動を持つ細胞質（Dictyostelium）対 社会性昆虫）における密度感受的な社会的行動を理解するための一般的な枠組みを提供し、分散型ロボットシステムへの応用可能性を示唆していると主張している。本研究は、変形可能な生物学的細胞をより良く表す軟らかなポテンシャルが、硬い球体モデルとは異なる相挙動をもたらすことを強調している。

限界と今後の方向性
本研究は、システムサイズ（単一サイズを使用）および運動性スイッチのバイナリ（二値）性に関する限界を認めている。著者は、今後の調査として以下を提案している：

• 複数の臨界密度や多段階の運動性変化を持つモデル。
• 連続的な運動性関数。
• 「活性円」相の次元効果を評価するための三次元実装。
• 生物学的なロッド（棒状粒子）で見られる形状効果（アスターやストライプの形成など）を捉えるための非球形粒子。