Aspect-ratio-dependent void formation in active rhomboidal and elliptical particle systems

技術要約：能動的な菱形および楕円形粒子系におけるアスペクト比依存的な空隙形成

問題提起
自己駆動粒子が化学エネルギーを運動エネルギーに変換する能動的ネマティック系は、複雑な時空間的秩序化と密度不均一性を示す。先行研究では、トポロジカル欠陥がこれらのダイナミクスを駆動すること、および粒子の形状（アスペクト比）が流体モデルにおける集団行動に影響を与えることが確立されているが、「ドライ」な能動的ネマティック（相互作用が直接的な衝突や排除体積効果によって支配される系）における粒子形状の具体的な役割については、まだ十分に理解されていない。特に、中密度系における「空隙領域（void regions）」（局所的な粒子密度がほぼゼロになる領域）の出現に関して、粒子の形状（単なるアスペクト比を超えた点）や直接的な接触力学が、空隙のサイズと安定性をどのように決定するかについては、完全には解明されていない。

手法
著者らは、周期境界条件を持つ正方形の系を用いた、二次元能動的ネマティックモデルの数値シミュレーションを行った。系は、$N=4140$個の自己駆動粒子で構成されている。

• 粒子の形状： 菱形と楕円形の二つの異なる幾何学的形状を比較した。
• パラメータ： 粒子の面積（$s = \pi/23$）は一定に保たれ、アスペクト比（$a$）が変化させた。これに伴い、各形状の長軸方向の長さも調整された。
• ダイナミクス： 粒子の運動は、自己駆動、ランダムな回転、およびランダムな変位を含む試行ステップを用いてモデル化された。重なりは拒絶され、排除体積の制約が強制された。メモリ変数（$\Lambda$）を用いて重なりの履歴を追跡し、その後の運動に影響を与えた。
• 解析： 本研究では、局所数密度、局所ネマティック秩序、トポロジカル欠陥の巻き数、および空隙領域の特性を解析した。空隙領域は、局所密度が系の平均の30%を下回るグリッドポイントの連結クラスターとして定義された。主要な指標には、平均空隙面積、空隙サイズの分布、および全域的および空隙周辺の粒子に特化した配向相関長が含まれる。

主な結果

1. 空隙の出現： 両方の系において、アスペクト比が臨界閾値（菱形では $a \approx 3.0$、楕円形では $a \approx 4.0$）を超えた場合にのみ、空隙領域が出現した。局所密度の分散と空隙領域の総面積はアスペクト比とともに増加し、最終的に飽和した。
2. 形状依存的な空隙特性：
   • 菱形系： 特徴的な空隙サイズを示し、それはアスペクト比とともに単調に増加した。空隙サイズの分布は狭く、総空隙面積は $a=5.0$ 付近で飽和した。
   • 楕円形系： 明確な特徴的スケスを伴わない、広範な空隙サイズの分布を示した。空隙の数密度は菱形系よりも低く、総空隙面積はアスペクト比に対して菱形系のような明確な飽和挙動を示さなかった。
3. トポロジカル欠陥： 両方の系は $\pm 1/2$ の巻き数を持つ欠陥によって支配されていた。しかし、菱形系はより高い密度の欠陥を含んでいた。空隙領域は両方の系において主に巻き数 $+1$ と関連しており、多くの場合、空隙の周囲に対して接線方向に整列した粒子に囲まれていた。
4. 配向相関： 配向相関長（$\lambda$）は、両方の系においてアスペクト比とともに増加した。極めて重要なことに、特に空隙の周囲に位置する粒子の相関長（$\lambda^{(void)}$）は、楕円形系よりも菱形系において有意に高かった。空隙の周辺における粒子の整列（$\langle \bar{\beta} \rangle$ で定量化）は、菱形系においてより強かった。

意義および主張
本論文は、ドライな能動的ネマティック系において、粒子の具体的な幾何学的形状（菱形か楕円形か）が、アスペクト比のみの効果とは異なる、マクロな集団パターンを決定する上で不可欠な役割を果たすと主張している。

• 空隙サイズのメカニズム： 著者らは、菱形系における空隙領域の特徴的なサイズは、空隙を取り囲む粒子場の配向相関長によって決定されると提案している。これは、観測された平均空隙面積と、配向相関長から推定された面積の比が、菱形においてはほぼ一定に保たれているという知見によって支持されている。
• 楕円形系との対比： 対照的に、楕円形粒子の場合は空隙サイズの分布が広いため、空隙サイズは局所的な配向相関長ではなく、系全体のサイズによって決定される可能性が示唆されている。
• モデルへの示唆： これらの結果は、アスペクト比のみに依存する（流体モデルで一般的である）能動的ネマティックのモデルは、直接的な衝突が支配的な系においては不十分である可能性を示唆している。粒子が衝突する際のパッキングや整列の仕方に影響を与える、粒子の詳細な幾何学的形状が、空隙形成を含む時空間パターンの出現における決定的な要因となる。

本研究は、アスペクト比と粒子形状の両方が集団運動に影響を与える一方で、空隙サイズを支配する具体的なメカニズムは、配向相関とパッキングの制約の違いにより、菱形粒子と楕円形粒子の間で根本的に異なるという結論を下している。