Topological delocalisation of confined 3D active nematics

技術要約：閉じ込められた3次元能動的ネマチックにおけるトポロジカルな非局在化

問題提起
能動的ネマチックは、移動するトポロジカル欠陥の自発的な生成と消滅によって特徴付けられる非平衡系であり、能動的乱流などの現象を駆動する。バルクの3次元能動的ネマチックにおける転位線のダイナミクスは、ポリマーメルトと同様のパーコレーションや絡み合いを示すことが理解されているが、閉じ込められた状態におけるこれらの系の挙動については依然として解明されていない。具体的には、幾何学的制約（曲率）と内部の活性（アクティビティ）との相互作用が、3Dの閉じ込められた設定においてどのように作用するのかが完全には特徴付けられていない。先行研究の多くは、閉じた曲面上の2D膜やバルク系に焦点を当てており、不均一な曲率を持つ閉じた3D幾何学構造内において、欠陥線が境界点とどのように接続し、バルクのダイナミクスを駆動するかという理解には空白が残されている。

手法
著者らは、能動的ネマチコ流体力学の連続方程式を解くために、ハイブリッド格子ボルツマン法を用いた数値シミュレーションを用いている。

• システム構成: シミュレーションでは、平面的な退化アンカリング（ディレクター場が表面に対して接線方向に整列する）を持つ閉じたシリンダー（およびスフェロシリンダー）内に閉じ込められた3次元能動的ネマチックをモデル化している。
• 物理モデル: 系は、シリンダー領域の位相場 $\phi$ と、液晶のテンソル秩序パラメータ $Q$ を含む粗視化された自由エネルギー密度によって記述される。ダイナミクスは、位相場のためのCahn-Hilliard方程式と、$Q$ テンソルのためのBeris-Edwards方程式によって制御され、これらは非圧縮ナビエ・ストークス方程式によって流体速度と結合している。
• パラメータ: 研究では、シリンダー半径 ($R$) と活性強度 ($\zeta$) を変化させている。無次元の活性数 $A = R\sqrt{\zeta/K}$ が、閉じ込めサイズと能動的長さスケール $\ell_\zeta = \sqrt{K/\zeta}$ の比として定義されている。
• 解析: 転位密度テンソルを用いて欠陥線が特定される。主要な指標には、欠陥占有率（局在状態と非局在状態を区別するため）、欠陥輪郭長、巻き付けパターン（ねじれ角 $\beta$ によって特徴付けられる）、および欠陥数と長さのスケーリング則が含まれる。

主要な貢献と結果

1. 弾性と活性の競合:
   本研究は、受動的な弾性と能動性の間の根本的な競合を明らかにしている。

   • 局在状態: 低活性（または平衡状態）では、受動的な弾性が支配的である。弾性変形を最小限に抑えるため、欠陥は曲率の高い領域（シリンダーのエッジやエンドキャップ）にピン留めされる。平衡状態では、2本の短い転位線がエッジ付近の $+1/2$ 端点を結んでいる。
   • 非局在状態: 活性が増加するにつれ、能動的な応力が弾性的な局在障壁を克服する。欠陥は運動性を獲得し、欠陥がシリンダー体積全体をカオス的に動き回る「トポロジカル非局在化」の領域へと至り、これは能動的乱流に類似している。

2. 相転移のシグネチャ:
   局在状態から非局在状態への転移は、臨界活性数 $A_c$ によって特徴付けられる。

   • ゆらぎ: この転移は、欠陥占有率の顕著なゆらぎを伴う。
   • メカニズム: 非局在化のメカニズムはシリンダー半径に依存する。小さな半径 ($R \le 20$) では、非局在化はシリンダー表面に沿った新しい欠陥対の核生成（不安定性駆動）によって起こる。より大きな半径では、欠陥は必ずしも新しい欠陥を即座に生成することなくエンドキャップから脱出し、シリンダーを横断して、最終的に完全に非局在化した状態へと至る。

3. 欠陥の形態と巻き付け:

   • 巻き付けパターン: 局在状態において、欠陥は純粋な $+1/2$ 転位である。非局在状態では、活性によって、ねじれ（$\cos \beta \approx 0$）や中間状態を含む幅広い巻き付けパターンが導入される。
   • 電荷分離: 非局在状態においても、「電荷分離」効果がエンドキャップ付近に持続する。局所的な巻き付け特性 ($\cos \beta$) はエンドキャップからシリンダー中心に向かって指数関数的に減衰し、その振幅は $A^{-2}$ に比例するが、特徴的な減衰長は活性に依存しない。

4. 閉じ込め下でのバルクとのスケーリング則:
   本論文は、バルクと比較した閉じ込め系における明確なスケーリング挙動を確立している。

   • 欠陥数: 平均的な欠陥線の数 $\langle N_d \rangle$ は $\ell_\zeta^{-2}$ に比例する。これは、線がシリンダー表面に拘束されているため、2次元系と同様の挙動を示す。これは、バルク（または弱い閉じ込め）における欠陥ループが $\ell_\zeta^{-3}$ にスケールすることとは対照的である。
   • 欠陥長: バルク系では欠陥長が能動的長さスケールに線形に比例する ($L_d \sim \ell_\zeta$) のに対し、閉じ込めは活性に依存しない特徴的な長さを導入する。平均的な欠陥長はシリンダー半径に従って $L_d \sim R^{5/3}$ とスケールする。
   • ポリマーとの類似性: 閉じ込め下での欠陥長の確率分布は指数関数的に減衰する（バルクにおけるべき乗分布とは異なる）。$R^{5/3}$ というスケーリング指数は、自己回避ランダムウォーク（Flory指数 $\nu = 3/5$）に対応しており、これは閉じ込められた「ブロブ」内の「リビングポリマー」と類似した挙動を示すことを示唆している。

意義
著者らは、不均一な曲率の閉じ込めが、能動的トポロジカルダイナミクスを制御するための多用途なメカニズムを確立したと主張している。活性数と境界幾何学を調整することで、欠陥のピン留め部位、渦形成、および非局在的乱流への転移を制御できる。トポロジカル絶縁体や超伝導体（磁束線がピン留めされるか脱出するか）に似た相転移の特定、および欠陥線に関するポリマー的なスケーリング則の提示は、3次元能動物質を理解するための新しい理論的枠組みを提供する。これらの結果は、幾何学的制約が特定の欠陥構成を安定化させたり、自発的な流れを誘起したりするために利用できることを示唆しており、3次元能動的ネマチックが存在する上皮の形態形成や癌組織などの生物学的システムへの適用可能性を示している。