Interface crossing behavior of prolate microswimmers: thermo and hydrodynamics

技術要約：長球形マイクロスイマーの界面横断挙動

問題提起
水圏環境における微生物の運動性は、周囲の環境（流体-流体界面を含む）からのフィードバックメカニズムに大きく影響される。界面とマイクロスイマー間の流体力学的相互作用については十分に文書化されているが、清浄な液-液界面を横断する能動的粒子のダイナミクスにおける熱力学的トラッピング力（特にピッカリング効果）の役割については、まだ十分に解明されていない。近年のCheonら（2024年）による実験では、桿菌状のBacillus subtilis（枯草菌）が、その速度と初期投影角に応じて、等方性-ネマチック液-液界面を横断するか、あるいは捕捉されるかのいずれかになることが観察された。本研究は、清浄な液-液界面における長球形（棒状）マイクロスイマーのトラッピング・ダイナミクスを調査し、界面の熱力学的力と能動的な流体力学的力の間の競争を合理的に説明することを目的とする。

手法
著者らは、清浄な液-液界面と長球形マイクロスイマーの相互作用をモデル化するために、大規模な3D格子ボルツマン法を用いたシミュレーションを採用した。

• スイマーモデル： マイクロスイマーは、ストークス流の下での長球形スカーマーモデルを用い、アスペクト比 $AR = 3$ の固体の長球形粒子としてモデル化された。スイマーは、Bacillus subtilisを模倣するため、「プッシャー」型（$\beta = -1$）として定義され、その境界には規定された滑り速度が設定された。
• 流体界面： 二成分混合流体は、相進化のためのCahn-Hilliard方程式および流体速度のためのNavier-Stokes方程式と結合した、Ginzburg-Landau自由エネルギー汎関数を用いてモデル化された。界面は中立的な濡れ性（接触角 $\approx 90^\circ$）として扱われた。
• シミュレーションパラメータ： シミュレーションは低レイノルズ数（$Re \sim O(10^{-2})$）で実施された。毛管数（$Ca = \eta u_0 / \sigma_0$）は、界面張力（$\sigma_0$）または遊泳速度（$u_0$）のいずれかを変化させることで変動させた。界面に対するスイマーの初期投影角（$\phi_0$）も系統的に変化させた。
• 解析： 粒子軌跡、配向ダイナミクス、および速度ゆらぎを解析した。熱力学的なトラッピングエネルギー（界面面積の減少に由来）と、界面を横断するために必要な能動的エネルギーとの間のエネルギーバランスに基づいた理論的枠組みを構築した。

主要な貢献と結果

1. 横断 vs トラッピングの相図： シミュレーションにより、毛管数（$Ca$）と初期配向角（$\phi_0$）の空間における相図を確立した。結果は、中立的な濡れ性を持つスイマーが、能動的な遊泳力と界面張力の競争に応じて、界面を横断するか、あるいは捕捉されるかのいずれかになることを示している。

   • トラッピング： 界面力が支配的な場合（通常は低 $Ca$ または浅い初期角度の場合）に発生する。スイマーは界面に対して平行に配向する。
   • 横断： 能動的な力が支配的な場合（高 $Ca$ の場合）に発生し、スイマーが界面を通過することを可能にする。
   • 検証： 数値的な相境界は、異なる物理的メカニズムが提案されているにもかかわらず、等方性-ネマチック界面を横断するBacillus subtilisの実験的観察結果（Cheon et al., 2024）と極めて良く一致している。

2. 熱力学的トラッピングメカニズム： 本研究は、トラッピングがピッカリング効果によって駆動されることを提案している。長球形のスイマーにとって、界面から粒子を取り除くためのエネルギーコストは接触面積に依存し、これは配向角 $\phi$ によって変化する。著者らは、能動的な仕事量に対する熱力学的トラッピングエネルギー（$\Delta E_{th} \propto \sigma_0 \cdot \text{Area}(\phi)$）のバランスに基づき、臨界毛管数（$Ca^*$）を導出した。この解析的な予測（$Ca^*$）は、シミュレーションで観察されたトラップ状態と横断状態の境界を正確に定義している。

3. 再配向ダイナミクス： 界面に対して平行になるスイマーの再配向は、流体力学的成分と熱力学的成分の両方を持つ複合的なトルクによって駆動されることが明らかになった。

   • 流体力学的トルク： プッシャー型スイマーの力双極子（$B_2$）から生じ、毛管数に比例する。
   • 熱力学的トルク： 界面自由エネルギーの勾配（$\tau_{thermo} \sim -\partial \Delta E_p / \partial \phi$）から生じ、配向角に依存する。
   • 観察： 浅い角度で初期配向したスイマーは、界面に平行にトラップされ再配向される前に、界面法線方向へ向かう初期の再配向を示すことがある（遠視野流体力学による）。再配向の速度は $Ca$ に対して線形に増加する。

4. ポテンシャルエネルギー地形： 速度ゆらぎを積分し、自由エネルギー差を計算することにより、著者らはトラッピングポテンシャルをマッピングした。ポテンシャルは界面付近で放物線状であり、その深さは固定された界面張力に対して比較的一定である。これは、低 $Ca$ における熱力学的トラッピングの支配を裏付けている。

意義と主張
本論文は、棒状マイクロスイマーの界面横断挙動に関する統一的な理論的および数値的な説明を提供すると主張している。その主な意義は以下の通りである：

• メカニズムの特定： 中立的な濡れ性を持つ長球形スイマーのトラッピングがピッカリング効果（界面エネルギーの熱力学的減少）によって支配されていることを示し、界面変形によるバネのような力をトラッピングの原因とした従来の実験的解釈と対比させた。
• 予測能力： エネルギーバランス方程式に基づく導出された臨界毛管数は、トラップ状態と横断状態の間の遷移を正確に予測しており、著者らのシミュレーションおよび独立した細菌実験の両方と優れた一致を示している。
• 二重のトルク起源： マイクロスイマーの界面における再配向は、単なる流体力学的な現象（球形スカーマーでしばしば仮定されるもの）ではなく、非球形粒子においては重要な熱力学的成分を伴うことを確立した。

著者らは、彼らのモデルがトラッピング力の具体的な物理的起源においてCheonらによる実験とは異なるものの、液-液界面におけるマイクロスイマーのマクロなダイナミクスと状態遷移を捉えることに成功していると結論付けている。