[論文レビュー] A tumor growth model of Hele-Shaw type as a gradient flow

Lénaïc Chizat, Simone Di Marino|arXiv (Cornell University)|Dec 17, 2017

Mathematical Biology Tumor Growth参考文献 20被引用数 6

ひとこと要約

本稿では、非負測度の空間に最適輸送・成長距離を備えた計量を備えた、ヘレーショウ型腫瘍成長モデルを、初期データの正則性仮定なしに解の存在・一意性を保証する勾配流として定式化する。主な貢献は、最小化移動によるモデルの変分的特徴付けであり、エントロピー正則化に基づく数値的に扱いやすいスキームが提案される。

ABSTRACT

In this paper, we characterize a degenerate PDE as the gradient flow in the space of nonnegative measures endowed with an optimal transport-growth metric. The PDE of concern, of Hele-Shaw type, was introduced by Perthame et. al. as a mechanical model for tumor growth and the metric was introduced recently in several articles as the analogue of the Wasserstein metric for nonnegative measures. We show existence of solutions using minimizing movements and show uniqueness of solutions on convex domains by proving the Evolutional Variational Inequality. Our analysis does not require any regularity assumption on the initial condition. We also derive a numerical scheme based on the discretization of the gradient flow and the idea of entropic regularization. We assess the convergence of the scheme on explicit solutions. In doing this analysis, we prove several new properties of the optimal transport-growth metric, which generally have a known counterpart for the Wasserstein metric.

研究の動機と目的

• 新しい最適輸送・成長距離を用いた、ヘレーショウ型腫瘍成長モデルの変分的・勾配流的定式化を提供すること。
• 初期条件の正則性を要件としない、剛性極限系の解の存在を確立すること。
• 滑らかでない領域上での解の一意性を、変動不等式（EVI）を用いて証明すること。
• 勾配流離散化のエントロピー正則化に基づく、安定かつ収束する数値スキームの開発。
• ワッサースタイン幾何学を超えて、非平衡輸送設定への勾配流理論の拡張。

提案手法

• 非負測度の空間において、関数 $ G(\rho) = -\lambda\rho(\Omega) $ が $ \rho \ll \mathcal{L}^d $ かつ $ \frac{d\rho}{d\mathcal{L}^d} \leq 1 $ を満たす場合に限り定義され、それ以外では $ +\infty $ となる関数の勾配流として腫瘍成長モデルを形式化する。
• 質量移動と成長／縮小の両方を考慮する最適輸送・成長距離を導入し、ワッサースタイン距離を一般化する。
• 時間離散化による最小化移動スキーム（時間方向の陰的オイラー法）を適用し、解を構成する：$ \tilde{x}_{k+1} \in \arg\min_x \left\{ G(x) + \frac{1}{2\tau} d^2(x, \tilde{x}_k) \right\} $。
• 非平衡輸送に適応したJKOスキームフレームワークを用いて、安定かつ収束する数値スキームを導出する。
• 離散化された最小化問題を計算的に扱いやすくするためにエントロピー正則化を用いる。
• 安定性および一致性の議論を用いて、明示解上での数値スキームの収束を証明する。

実験結果

リサーチクエスチョン

• RQ1 ヘレーショウ型腫瘍成長モデルは、非負測度の空間における距離計量上での勾配流として解釈可能か？
• RQ2 提案された最適輸送・成長距離は、初期データの正則性仮定なしに解の存在・一意性を支持するか？
• RQ3 勾配流定式化は、安定かつ収束する数値スキームにどのように離散化可能か？
• RQ4 最適輸送・成長距離とワッサースタイン距離との構造的類似点および相違点は何か？
• RQ5 モデルの変分的構造は、腫瘍成長モデリングのための新たな解析的・数値的ツールの導出にどのように活用可能か？

主な発見

• ヘレーショウ型腫瘍成長モデル（1.2式）は、非負測度の空間に最適輸送・成長距離を備えた計量を備えた勾配流として、関数 $ G $ の勾配流として特徴付けられる。
• 初期条件 $ \rho_0 \in L^1_+(\Omega) $ に対して、正則性仮定なしに最小化移動スキームを用いて解の存在が確立される。
• コンパクトな凸領域上では、変動不等式（EVI）の検証により解の一意性が証明される。
• 最適輸送・成長距離の新たな性質が導出され、時間離散近似の連続性および収束性が確認される。
• 勾配流のエントロピー正則化に基づく数値スキームが提案され、明示解上での収束性が検証される。
• 変分的解釈により、機械的直感が得られる：モデルは最大密度制約下で、質量移動と成長率の両方を最小化する最適な腫瘍の成長経路を記述する。