Hydrodynamic stability and pattern formation in hexatic epithelial layers

技術要約：六角形エピテリウム層における流体力学的安定性とパターン形成

問題提起
本論文は、変形可能な細胞の密充填に由来する六角形秩序（6回回転対称性）に焦点を当て、エピテリウム層の流体力学的安定性とパターン形成について論じている。極性（1回回転対称）やネマチック（2回回転対称）の系における能動的液晶流体力学は、細胞の集団移動や欠陥ダイナミクスを説明するために広く適用されてきたが、六角形エピテリウムの特有の流体力学的挙動については未だ十分に解明されていない。中心的な課題は、機械的な活性、横方向の接着、および六角形秩序のマルチスケールな性質の相互作用が、静止状態の安定性と集団運動の出現にどのように影響するかを理解することである。不安定性が一般的な能動的ネマチックとは異なり、六角形系では、異なる長さスケールにおけるエネルギー注入と散逸の間のより複雑なバランスが存在する可能性がある。

手法
著者らは、受動的な$p$-atic液晶理論から導出され、能動的な六角形系（$p=6$）へと拡張された連続体流体力学モデルを用いている。

• 支配方程式: モデルは、質量密度、運動量密度、および局所的な六角形配向を表すランク6の秩序パラメータテンソル（$Q_6$）を結合させている。能動応力は、細胞スケールで発生する力を反映して $\sigma^{(a)} = \alpha_6 \nabla^{\otimes 4} \odot Q_6$ とモデル化されている。
• 次元削減: 研究では、以下の2つのレジームに焦点を当てている：
  1. チャネル閉じ込め: 幅 $L$ の2Dチャネルにおいて、系は縦方向に対して不変である。これにより、問題は六角形角度 $\theta(y)$ と縦方向速度 $v_x(y)$ に関する結合方程式へと簡略化される。
  2. 周期境界領域: 完全なベクトル流および配向ダイナミクスを許容する2D周期領域。
• 散逸の正則化: 重要な手法上のステップとして、小スケールにおける標準的な粘性モデルの非物理的な不安定性を対処することが挙げられる。著者らは、エネルギー散逸スケールが小スケールでのエネルギー注入率（$\sim q^6$）と一致するように、ハイパー粘性（高次の速度勾配、具体的には $\eta_6 \nabla^4$）を導入し、発散を防いでいる。これにより、標準的なせん断粘性（$\eta_2$）とハイパー粘性（$\eta_6$）の両方を含む、スケール依存の散逸率が導かれる。
• 解析: 著者らは、静止状態に対する線形安定解析を行い、分散関係を導出した。彼らは、デボラ数（$De$）と、粘性散逸とハイパー粘性散逸のレジーム間の遷移を決定するクロスオーバー長さスケール$\ell_\times = (\eta_6/\eta_2)^{1/4}$ を含む、重要なパラメータを特定している。全非線形方程式の数値シミュレーションを用いて、解析的な予測を検証し、遷移後のダイナミクスを探索している。

主要な貢献と結果

1. スケール依存の安定化メカニズム:
   本論文は、六角形エピテリウムの流体力学的安定性が、能動的ネマチックのようにシステムサイズのみによって決定されるのではなく、長さスケールにわたるエネルギー注入と散逸のバランスによって決定されることを示している。

• 生成的安定性: 横方向の接着相互作用（$\eta_2$ で表される）が強い場合、系はサイズに関わらず生成的に安定である。エネルギーが注入スケールにおいて効率的に散逸されるため、静止状態は微小な摂動に対しても安定に保たれる。
• 不安定性の閾値: 不安定性は、活動性が散逸に対して十分に高い場合にのみ発生し、これはクロスオーバースケール $\ell_\times$ における臨界デボラ数（$De_\times > 16/3$）によって特徴付けられる。

2. パターン形成不安定性の階層構造:
   閉じ込められたチャネル内では、系は**逆方向流レーン（counter-flowing lanes）**の形成に至る一連の不安定性を示す：

• 自発的流動転移: 能動的ネマチクスと同様に、一次的な不安定性が細胞の逆方向の流れを持つ2つのレーンの形成をもたらす。
• レーンの増殖: ネマチクスとは異なり、複数の固有の長さスケールの存在により、さらなる不安定性の階層が存在する。チャネル幅 $L$ が増加するか、あるいは活動性が上昇するにつれて、系は増加する数のレーンへと分割される（$n_0 \approx 2\lfloor q_0/q_{min} \rfloor$）。
• 固定されたレーン幅: 決定的なことに、これらのレーンの幅は、システムサイズに依存せず、クロスオーバー長さ $\ell_\times$ に比例する、スケールに依存しない値に収束する。これは、レーン構造が材料のレオロジーにおける固有の特性であることを示唆している。

3. 二次元ダイナミクスと「能動的乱流」:
   2D周期領域では、1Dのレーン構造は渦の増殖に対して不安定になる。

• 低活動時: 系は、波長 $\sim \ell_\times$ を持つ静止した秩序バンド（ヘリンボーン・パターン）を形成する。
• 高活動時: 系は時空間カオスの状態へと転移する。このレジームは「能動的乱流」と呼ばれるネマチクスで見られる特徴を共有しているが、著者らは明確な相違点についても述べている。このカオス状態は、線形不安定性の「記憶」を保持しており、特定の悪性間葉系腫瘍（例：線維肉腫）で見られるヘリンボーン・パターンに類似した、配向と渦構造のクラスターを示す。

意義と主張
著者らは、本研究が六角形エピテリウムの流体力学的安定性に関する初の包括的な解析を提供すると主張している。主な意義は、六角形秩序がマルチスケールな組織化をもたらし、それが能動的流体の安定性基準を根本的に変えることを確立した点にある。

• ネマチックとの対比: 本論文は、バルクの能動的ネマチクスにおいて、ドメインが十分に大きければいかなる有限の活動性でも系を不安定化させるという事実に対し、鮮明な対比を示している。六角形エピテリウムでは、強い横方向の接着（高い $\eta_2$）によって、システムサイズに関わらず静止状態を無期限に安定化させることができる。
• 生物学的関連性: 本モデルは、エピテリウム組織の機械的安定性が、横方向の接着（例：E-カドヘリンの発現）を通じて調整可能であることを示唆している。固定幅のレーンやヘリンボーン・パターンの出現は、複雑な細胞レベルのルールを導入することなく、移動する細胞クラスターの空間的組織化や、特定の病理で見られる形態学的パターンを理解するための理論的枠組みを提供する。
• 理論的枠組み: 本研究は、$p$-atic秩序とスケール依存の散逸の相互作用が複雑な集団行動の出現を支配するという連続体フレームワークを確立しており、これは極性またはネマチックの典型的な単一秩序パラメータによる記述を超えたものである。