Multicellular simulations with shape and volume constraints using optimal transport

技術的概要：最適輸送を用いた形状・体積制約を伴う多細胞シミュレーション

問題提起
細胞内構成要素から巨視的組織、細菌コロニーに至るまで、生体システムは体積排除、形状相互作用、変形性によって駆動される複雑なダイナミクスを示す。これらの制約がどのように巨視的な自己組織化をもたらすかを理解することは、発生生物学および生物物理学における根本的な課題である。既存の計算モデルは、物理的精度、数値的コスト、数学的扱いやすさの間にトレードオフを抱えている。一般的なアプローチには以下が含まれる：

• 点粒子系： 数学的に単純であるが、しばしば現実的な形状表現に欠ける。
• 位相場モデルおよびレベルセットモデル： 任意の形状を表現可能であるが、計算コストが高く、接触相互作用に対して困難を伴う場合がある。
• セルオートマトンおよび格子モデル： 空間をボクセルに離散化するため、幾何学的柔軟性が制限される可能性がある。
• ボロノイおよび頂点モデル： 組織様構造に対して効率的であるが、しばしば体積に対してソフト制約に依存するか、特定の形状変形に対する制御が欠如している。

厳密な体積制約、任意の動的形状、および高粒子数（数十から数万件）を、個々の変形性と集団的相互作用の両方を支援する統合フレームワーク内でモデル化する上で、顕著なギャップが存在する。

手法
著者らは、最適輸送理論、特にラグジュア分割（Laguerre tessellations）に基づいたメソスケールシミュレーションフレームワークを導入する。

1. 静的モデル（ラグジュア分割）：
   各粒子 $i$ が占める空間は、ラグジュアセル $L_i$ としてモデル化され、以下の条件を満たす点 $x$ の集合として定義される：
   $$L_i = \{x \in \Omega, c(x, x_i) - w_i \leq c(x, x_j) - w_j \text{ for all } j\}$$
   ここで、$x_i$ は粒子の位置、$v_i$ はその目標体積、$c(x, x_i)$ は粒子の形状および変形性の性質を符号化するコスト関数である。重み $w_i$（カンタロビッチポテンシャル）は、厳密な体積制約 $|L_i| = v_i$ を満たすように一意に決定される。

   • この定式化は、$w_i=0$ かつ $c$ がユークリッド距離であるボロノイ図を一般化し、任意のコスト関数（例えば、パワーカー、異方性コスト）を選択することで、曲がった境界および非球形形状を可能にする。
   • 空の空間は、ゼロコストを持つ追加の粒子としてモデル化される。

2. 動的モデル：
   粒子位置 $x_i$ の進化は、一次の勾配降下方程式に従う：
   $$\dot{x}_i = -\tau_i \nabla_{x_i} T_c + F_i^{\text{point}} + F_i^{\text{surf}} + \dot{\xi}_i$$

   • 非圧縮性力： 項 $-\nabla_{x_i} T_c$ は、総輸送コスト $T_c$ から導かれる反発力として機能する。これは粒子が重なり合わないことを保証し、体積制約を自然に強制する。対力とは異なり、この力はすべての接触面のグローバルな幾何学に依存する。
   • 外部・相互作用力（$F^{\text{point}}$）： 点粒子モデルと同様に、自己推進、重力、または引力を含む。
   • 表面力（$F^{\text{surf}}$）： 界面に沿った圧力および表面張力効果をモデル化し、細胞選別や接着のような現象を可能にする。
   • 形状ダイナミクス： 形状は、パラメータ化されたコスト（例えば、共分散行列で定義される楕円体）や、外部場（例えば、走化性）に応答するレベルセットポテンシャル関数を介して進化し得る。

3. 数値実装：

   • 問題は、離散グリッド（ボクセル）上の半離散最適輸送問題として解かれる。
   • 著者らは、カンタロビッチポテンシャル $w_i$ を効率的に計算するために、分割統治マルチスケール戦略および**準ニュートンソルバ（L-BFGS-B）**を採用する。
   • 実装は KEOPS および GEOMLOSS ライブラリを用いて GPU で大規模並列化されており、最先端のパフォーマンスで 3 次元における最大 50,000 粒子のシミュレーションを可能にする。

主要な貢献

• 統合フレームワーク： この手法は、ボロノイ分割、レベルセット法、頂点モデルの概念を最適輸送の傘下に統合する。単一のフレームワーク内で任意の形状、厳密な体積制約、および複雑な相互作用メカニズムを支援する。
• 厳密な体積強制： ソフト制約または緩和力を使用する多くの既存モデルとは異なり、このアプローチはカンタロビッチポテンシャルの最適化を通じて、常に数学的に体積制約を保証する。
• スケーラビリティ： GPU 加速ソルバは、多くの従来の有限要素法や位相場法では計算的に不可能な、3 次元における大規模システム（最大 50,000 粒子）のシミュレーションを可能にする。
• モデリングの汎用性： このフレームワークは、以下の多様な現象を成功裡にモデル化する：
  • 変形性によって駆動される活性ブラウン粒子における流体 - 固体相転移。
  • 体積排除のみによる群れをなす棒状粒子における創発的な配向秩序。
  • 形状変形によって駆動される走化性移動。
  • 異なる表面張力と軟らかさを持つ 2 種類の細胞の 3 次元不均一凝集体における細胞選別および取り込み現象。これは、界面張力差仮説（DITH）と整合する結果を再現する。

結果
著者らは、いくつかのベンチマークを通じてフレームワークを検証する：

• 成長と分裂： 3 次元凝集体は、単一の細胞から分裂を介して 50,000 細胞へと成長し、約 1 日間の計算でドメインを充填する。
• 相転移： 活性ブラウン粒子の系において、変形性パラメータ $\alpha$ は、流体相（拡散的挙動）から固体相（ジャムした六方最密充填）への転移を駆動する。このモデルは、固体相における実効拡散の消失および平均二乗変位の飽和を捉える。
• 棒状配列： 自己推進する棒状粒子のシミュレーションは、衝突誘発変形に起因する長距離配向秩序の自発的創発を示し、高分子物理学の理論と一致する。
• 3 次元細胞選別： このモデルは、異なる表面張力と軟らかさを持つ 2 種類の細胞の 3 次元混合物において、各種の選別パターン（取り込み、分離、チェッカーボード）を再現し、理論的予測および実験的観察と一致する。
• 群衆ダイナミクス： このモデルは群衆の避難をシミュレートし、変形性粒子（流体様）は安全に避難するが、硬球粒子（固体様）は安定したアーチを形成し出口でジャムすることを示す。

意義と主張
本論文は、詳細な物理記述と計算効率の間のギャップを埋める多細胞システムをモデル化するための新規アプローチを導入すると主張する。

• 検証： 主要な貢献は、古典的計算モデルに対するこの新フレームワークの検証であり、体積制約および 3 次元複雑性の処理において優位性を提供しつつ、確立された結果（細胞選別、相転移など）を再現できることを示す。
• 生物学的関連性： 著者らは、個々の細胞変形と高密度凝集体の挙動の両方が重要な、胚盤胞の発生などの複雑な発生過程の in silico モデリングにおけるこの手法の可能性を強調する。
• 理論的拡張： この研究は、既存の手法（点粒子、レベルセット、頂点）に対する新たな最適輸送の視点を提供し、体積相互作用を有するシステムに対する粗視化連続モデル（平均場近似）の導出への道筋を示す。これは数学的生物学において未解決の課題である。
• アクセシビリティ： 著者らは自由に利用可能な Python 実装を提供しており、計算物理学および生物学におけるさらなる研究を可能にする。

著者らは謙虚であり、フレームワークは新しいものであるが、即座の焦点は既存モデルに対する検証にあったと指摘する。彼らは、アクティブ頂点理論をラグジュア分割に拡張すること、連続極限の導出を精緻化すること、リアルタイムシミュレーションのための複雑な動的コストおよび表面相互作用の処理を改善することなどの将来の方向性を特定している。