[論文レビュー] Variational formulation of active nematics: theory and simulation

Waleed Mirza, Alejandro Torres-Sànchez+2 more|arXiv (Cornell University)|Jun 2, 2023

Micro and Nano Robotics被引用数 5

ひとこと要約

本稿は、Onsagerの形式に基づく可積み2次元アクティブネマチックスの変分フレームワークを導入し、自由エネルギー、散逸、活性を統合して支配方程式と適合する有限要素法を導出する。密度依存のアクティブネマチックダイナミクスの正確なシミュレーションを可能にし、治癒過程や閉じた細胞集団における欠陥を介した自己組織化を明らかにする。

ABSTRACT

The structure and dynamics of important biological quasi-two-dimensional systems, ranging from cytoskeletal gels to tissues, are controlled by nematic order, flow, defects and activity. Continuum hydrodynamic descriptions combined with numerical simulations have been used to understand such complex systems. The development of thermodynamically consistent theories and numerical methods to model active nemato-hydrodynamics is eased by mathematical formalisms enabling systematic derivations and structured-preserving algorithms. Alternative to classical nonequilibrium thermodynamics and bracket formalisms, here we develop a theoretical and computational framework for active nematics based on Onsager's variational formalism to irreversible thermodynamics, according to which the dynamics result from the minimization of a Rayleighian functional capturing the competition between free-energy release, dissipation and activity. We show that two standard incompressible models of active nemato-hydrodynamics can be framed in the variational formalism, and develop a new compressible model for density-dependent active nemato-hydrodynamics relevant to model actomyosin gels. We show that the variational principle enables a direct and transparent derivation not only of the governing equations, but also of the finite element numerical scheme. We exercise this model in two representative examples of active nemato-hydrodynamics relevant to the actin cytoskeleton during wound healing and to the dynamics of confined colonies of elongated cells.

研究の動機と目的

• 密度の変動と分解・再合成を伴う可積み2次元アクティブネマチック系の熱力学的に整合したフレームワークの構築。
• 非平衡熱力学に根ざした変分原理に基づき、アクティブネマチックモデルの物理的解釈を明確化すること。
• 支配方程式および適合する有限要素数値スキームの体系的導出。
• アクチン細胞骨格の治癒過程におけるダイナミクスや閉じた細長い細胞集団における欠陥の進化といった生物学的に関連する状況のシミュレーション。
• メッシュサイズおよび時間ステップにおける空間的・時間的離散化の収束解析を通じた数値手法の妥当性検証。

提案手法

• 時間発展が自由エネルギーの変化、散逸、活性の組み合わせを最小化する汎関数を含む、Onsagerの変分原理に基づくダイナミクスの定式化。
• 密度 ρ、ネマチック秩序パラメータ S、配向場 n を状態変数として導入し、それらに関連するフラックスおよび時間変化率を定義。
• 自由エネルギー、散逸、活性仕事の項を含むラグランジアンに変分原理を適用し、支配方程式を導出。
• 変分構造を保つ有限要素法（FEM）を構築し、エネルギー安定性と一貫性のある空間離散化を確保。
• 時間に陽的でないガラーキンFEMスキームを用い、適応的時間ステップを採用。メッシュサイズおよび時間ステップにおける収束解析により妥当性を検証。
• 2つの生物学的系にフレームワークを適用：大型卵細胞における治癒過程と、閉じた細胞集団における欠陥ダイナミクス。

実験結果

リサーチクエスチョン

• RQ1 密度依存効果を伴う可積み2次元アクティブネマチックスの統一的で熱力学的に整合したフレームワークは、どのように構築可能か？
• RQ2 Onsagerの変分形式を用いることで、アクティブネマチックスにおける運動方程式と数値スキームの両方の導出にどのような影響を与えるか？
• RQ3 可積み性と分解・再合成を組み込むことで、アクティブ系におけるネマチック構造の自己組織化にどのような影響が生じるか？
• RQ4 トポロジカル欠陥とアクティブ流れは、大型卵細胞の治癒ダイナミクスにおいて果たす役割は何か？
• RQ5 活性強度と異方性は、閉じたアクティブネマチック系における欠陥ダイナミクスと流れのパターンにどのように影響を与えるか？

主な発見

• 変分フレームワークは、標準的な非可積みアクティブネマチックスモデルを正確に回復するとともに、可積み性と密度依存性を有する系へと拡張している。
• 有限要素スキームは最適な空間収束（収束率 ∼h²）と時間収束（収束率 ∼Δt）を示し、数値的安定性を確認している。
• 治癒モデルでは、アクティブネマチック流れが欠陥の移動と再組織化を駆動し、大型細胞における細胞骨格の再組織化を促進している。
• 閉じた集団では、活性度の増加に伴いアクティブ乱流が発現し、欠陥ダイナミクスは ℓₐ/ℓ₀ の比に支配されている。ここで ℓₐ はアクティブネマチック長尺度である。
• 異方性活性（κ = λ_aniso/λ で定量化）は欠陥の運動および流れのパターンを顕著に変化させ、κ の値が高くなると指向性の流れと欠陥ペアリングが促進される。
• ネマチック秩序、活性、流体力学的流れの相互作用によって、アスターレ、タクトイド、ネマチックバンドの自己組織化が自然に出現するのをモデルが捉えている。