Dimension-dependent continuum limits in tissue mechanics

技術要約：組織力学における次元依存的な連続体極限

問題提起
離散的な細胞ベースのモデルは、微視的な力を組織スケールの進化へと効果的に結びつけるが、細胞数が増えると計算量的に困難になることが多く、パラメータ化も容易ではない。連続体モデルは、集団レベルの推論や創発的な輸送特性の解明のための実用的な代替案となる。しかし、上皮組織力学における1次元の連続体極限はよく研究されている一方で、対応する高次元の問題はほとんど未開拓のままである。空間次元が、細胞スケールの力と組織スケールの輸送との間のリンクにどのように影響するか、特に $D(\rho) \propto \rho^\gamma$ 型の非線形拡散係数に関して、決定的なギャップが存在する。既存の流体力学や集団生物学における現象論的な応用は、その多くが $\gamma > 0$ に限定されているが、上皮力学は $\gamma > 0$ と $\gamma < 0$ の両方を動機付ける。著者らは、マクロな輸送指数 $\gamma$ が微視的な力学のみによって決定されるのか、それとも空間次元によって根本的に変化するのかを調査している。

手法
著者らは、離散的な力学モデルとマクロな偏微分方程式（PDE）モデルの間で、定常状態付近における長期的な指数関数的緩和率を一致させる戦略を採用している。

1. 離散モデル:

   • 1D: 組織は、固定された境界を持つ $N$ 個の機械的なリンク（細胞）の連鎖としてモデル化される。節点は、相互作用力則 $f(\ell)$ に従って駆動される過減衰モデルに従って進化する。著者らは、パラメータ化された力則のファミリー $f(\ell) = k(\ell^\nu - a^\nu)/(\nu a^{\nu-1})$ （$\nu=0$ の場合は対数限界）を利用しており、これは静止長 $a$ の近くでフックの法則に線形化される。
   • 2D: 細胞は、節点が隣接する節点間で二粒子的な機械的力を介して相互作用する三角形分割ネットワークとしてモデル化される。動力学は、隣接する節点への力を合計する同様の過減衰方程式によって支配される。
   • 解析: 著者らは、離散系を等間隔の定常状態の周りで線形化する。長期的な緩和は、系のヤコビ行列（または2Dにおけるラプラシアン行列）の主要な固有値によって支配され、指数関数的な減衰率をもたらす。

2. マクロモデル:

   • 組織密度 $\rho(\mathbf{r}, t)$ は、非線形拡散方程式 $\partial_t \rho = \nabla \cdot (D(\rho) \nabla \rho)$ によって記述される。
   • 著者らは、一般的な形式の非線形拡散係数 $D(\rho) = a^\beta \rho^\gamma / \tau$ を仮定している。ここで、$a$ は微視的な長さスケール、$\rho$ は密度、$\tau$ は特性時間スケールである。
   • 定常状態 $\bar{\rho}$ の周りでPDEを線形化することにより、著者らは緩和率が $D(\bar{\rho})$ とドメインの幾何学的形状（具体的にはラプラシアンの固有値）に依存することを導き出す。

3. 粗視化手順:

   • コアとなる手法は、離散シミュレーションから得られた緩和率（微視的パラメータ $k, \eta, a, \nu$ および次元 $d$ に依存）と、マクロなPDEから予測される緩和率を等置することである。
   • 微視的パラメータを一定に保ったまま、定常密度 $\bar{\rho}$ を変化させることによって（$N$ またはドメインサイズ $L$ を変更することで）、$D(\rho)$ の関数形式および指数 $\gamma$ を推定する。

主な結果

• 指数の次元依存性: 本研究は、非線形拡散係数の指数 $\gamma$ が、微視的な力則の指数 $\nu$ と空間次元 $d$ によって共同で決定されることを確立している。
• スケーリング関係: 次元解析とシミュレーションデータは、普遍的なスケーリング関係 $\gamma = \beta - 2/d$ を明らかにしている。ここで、$\beta$ は微視的な力則に関連するスケーリング指数である。
  • $\nu=1$（フックの法則）を持つ力則を用いた1D（$d=1$）では、結果として得られる拡散係数は $D(\rho) \propto \rho^{-2}$ となる。
  • 同じ微視的力則（$\nu=1$）を用いた2D（$d=2$）では、$D(\rho) \propto \rho^{-1.11}$ を与える。
  • 一般に、次元 $d$ が増加するにつれて、非線形性の大きさ $|\gamma|$ は減少する。
• 線形拡散の出現: $d \to \infty$ とすると $\gamma \to 0$ となり、同じ微視的力学に対して、高次元では線形拡散がますます重要になることを示唆している。
• 検証: 1Dおよび2Dにおける数値シミュレーションは、推定された $D(\rho)$ が離散的な緩和ダイナミクスを正確に再現することを裏付けている。これらの結果は、様々な初期条件やドメインサイズにわたって成立しており、マクロな輸送がグローバルなドメインスケールではなく、局所的な長さスケール（$a$ および $\rho$）に依存することを裏付けている。

意義および主張
本論文は、組織力学における創発的なマクロ輸送現象のための新しい物理的メカニズムを明らかにしていると主張している。主な意義は、次元は単なる幾何学的な詳細ではなく、粗視化後に得られるマクロな法則を根本的に変えるものであることを示した点にある。

• 1Dモデルの限界: 1Dモデルから得られた知見は、高次元においてそのまま保持されるとは限らない。具体的には、線形モデル（フックの法則）の使用は、微視的な長さスケールがどのように創発的な組織スケールの輸送へとマップされるかについての重要な詳細を隠蔽してしまう。
• 新しい粗視化手法: 著者らは、緩和率を比較することに基づいて、微視的なダイナミクスを非線形拡散関数へと粗視化する手法を提案している。このアプローチは、伝統的な連続体極限に依存しないため、有限の系に対して特に適している。
• 物理的メカニズム: 本研究は、なぜ多孔質媒体型の非線形拡散係数が組織力学において生じるのか、そしてどのように特定の指数が微視的な力学と空間次元の相互作用によって固定されるのかについて、物理的な説明を提供している。

著者らは、同じ微視的な相互作用則が次元に依存しない連続体極限を生成するのではなく、次元が有効な非線形拡散率の形状を変え、$D(\rho) \propto \rho^\gamma$ における指数を変化させるという結論を下している。この現象は、ランダムウォーク理論や繰り込み群の手法における次元依存的な振る舞いに類似している。