Ordering and Defect Dynamics in Passive and Active Nematopolars

技術要約：パッシブおよびアクティブなネマトポーラーにおける秩序形成と欠陥ダイナミクス

問題提起
本論文は、極性（head-tail非対称）とネマティック（head-tail対称）の相互作用が共存する「ネマトポーラー（nematopolar）」系における、非平衡集団行動の理論的理解に取り組んでいる。このような系は、強誘電液晶、リビング・リキッドクリスタル、細菌コロニー、マイクロチューブル・モーター混合物など、多様な生物学的および合成的コンテキストで観察されている。従来の理論的アプローチは、二成分モデル（極性ベクトルとネマティックテンソルの結合）を利用するか、あるいはアクティブ・タービュランスに焦点を当ててきたが、乾燥（dry）ネマトポーラー系における最小単一場モデルを用いた、相秩序化および粗大化（coarsening）ダイナミクスの系統的な特性評価については、未解明のままである。具体的には、競合する対称性の相互作用と、自己駆動（アクティビティ）が欠陥ダイナミクスおよびドメイン成長則に果たす役割が探索されていない。

手法
著者らは、ネマトポーラー系のための最小限の単一場「乾燥」モデルを提案し、数値的に調査している。

• モデル: 系は、パッシブ極限において純粋な緩和ダイナミクス（モデルA様）に従って進化する単一の極性ベクトル場 $\mathbf{p}$ によって記述される。自由エネルギー汎関数には、バルク偏極のためのダブルウェル・ポテンシャル、極性弾性項（$\propto |\nabla \mathbf{p}|^2$）、および $\mathbf{p}$ から構成されるテンソル（$\hat{P} \propto \mathbf{p}\mathbf{p}^\top - |\mathbf{p}|^2 I/2$）に由来するネマティック的な項が含まれる。
• アクティビティ: 非平衡挙動は、運動方程式における自己随伴項 $\Lambda(\mathbf{p} \cdot \nabla)\mathbf{p}$ を通じて導入される。この項により、偏極は秩序パラメータであると同時に局所的な速度場としても機能することが可能となり、これは乾燥アクティブ物質のトナー・テュー（Toner-Tu）理論と整合している。
• シミュレーション: 著者らは、周期境界条件を持つ2次元正方格子（$N=512, 1024$）上で有限差分法を用いている。
• 観測量: 研究では、トポロジカル欠陥（整数および半整数電荷）の位置、数、および型、ならびに極性およびネマティックの相関関数（$C_p, C_n$）、および特徴的な相関長さ（$L_p, L_n$）を追跡している。
• アプローチ: 調査は2つの段階で行われる。(1) 単独の欠陥ペアとループの制御されたシミュレーションによる、素的な緩和メカニズムの解析、および (2) 無秩序な初期状態から開始される大規模シミュレーションによる、グローバルな粗大化ダイナミクスの研究である。

主要な貢献と結果

1. 欠陥形態とストリングを介した相互作用（パッシブの場合）:

   • 極性とネマティックの結合が均衡しているパッシブ領域（$\Lambda=0$）では、系は自発的に脱偏極ストリング（偏極が消失する線）を形成し、これが半整数欠陥（$\pm 1/2$）を接続したり、閉じた脱偏極ループを形成したりする。
   • 欠陥ペア: 反対の電荷を持つ半整数欠陥がストリングによって接続されている場合、相対的な極性/ネマティック強度比（$\bar{k} = k_p/k_n$）の関数として、非単調な消滅速度を示す。これは、有効弾性定数の減少と、ストリング自体のエネルギーコストとの間の競争に起因する。
   • 同符号の欠陥: 同符号の半整数欠陥が、純粋なネマティック系では無限に反発し続けるのに対し、本モデルでは、クーロン的な斥力と接続するストリングの張力のバランスによって決定される有限の平衡間隔に達する。
   • ループのダイナミクス: 脱偏極ループは、パラメータに応じて2つの異なるメカニズムで緩和される：(i) 収縮（曲率駆動であり、ネマティック優位の領域で支配的）、および (ii) 回転緩和（ループ内部の偏極場が外部と一致するように再配向するもので、極性優位の領域で支配的）。極端な極性領域では、ループは整数電荷のストリングへと破断することがある。

2. 粗大化ダイナミクス（パッシブの場合）:

   • 無秩序な状態から開始すると、系は粗大化を起こす。特徴的な長さスケール $L(t)$ は $L(t) \sim (t/\ln t)^{1/2}$ として成長する。
   • このスケーリングは、連続的な回転対称性と点状の欠陥を持つ系（2D XYモデルやO(2)/Z2 ネマティック系に類似）と一致しており、ストリングやループの存在が基本的なダイナミック・スケーリング則を変えることはないが、純粋なネマティック系と比較して秩序化の絶対的な速度を遅らせることを裏付けている。
   • 全欠陥数は $n_{tot} \sim \ln t / t$ と減少する。

3. アクティブ・ダイナミクスと粗大化の停止:

   • 自己随伴項（$\Lambda \neq 0$）を導入すると、系は運動誘起電荷対称性の破れ（MICS: Motility-Induced Charge Symmetry Breaking）を示す。正負の欠陥の間の対称性が破れ、正の整数欠陥（アスター型）と負の半整数欠明（トレフォイル型）の共存が生じる。
   • 負の整数欠陥と正の半整数欠陥は、安定な種へと分解され、動的に抑制される。
   • 粗大化の停止: 十分に強いアクティビティは、特徴的な長さスケールが飽和（プラトーに到達）し、欠陥密度が有限に保たれる定常状態をもたらす。系は「粗大化の停止（arrested coarsening）」領域に入り、欠陥の進化が強く抑制されることで、一様な秩序状態への到達が妨げられる。

意義
本論文は、ネマトポーラー系における秩序形成と欠陥ダイナミクスを理解するための、単純かつ統一的な枠組みを提供すると主張している。最小限の単一場モデルを用いることで、著者らは以下を実証している：

• 複数の秩序パラメータを必要とせず、極性とネマティックの自由エネルギー項の競合から、複雑な欠陥構造（ストリング、ループ、および混合電荷の共存）が創発し得ること。
• 単純な自己随伴項として実装されたアクティビティが、欠陥を介した秩序形成を再形成し、非平衡状態（粗大化の停止）を安定化させ、パッシブ系には存在しない対称性の破れを誘起する根本的なメカニズムとなること。
• これらの結果は、強誘電液晶、リビング・リキッドクリスタル、マイクロチューブル・モーター混合物など、両方の種類の秩序が同一の自由度から創発する様々な生物学的および合成的システムの振る舞いを解釈するための理論的基礎を提供する。

著者らは、本モデルは最小限のものであるものの、これらの系の本質的な現象論を捉えており、競合する対称性とアクティビティの役割を明確に分離していると明記している。彼らは、このフレームワークが、競合する配向秩序を持つより広範なクラスのアクティブ物質を調査するための出発点として機能することを提案している。