[論文レビュー] Microscopic derivation of the Keller-Segel equation in the sub-critical regime

Ana Cañizares García, Peter Pickl|arXiv (Cornell University)|Mar 13, 2017

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

ひとこと要約

本稿は、正則化されたクーロン相互作用を伴う$N$個の相互作用粒子から成る確率的系において、2次元のケラー-セゲル方程式を、臨界的でない領域（$\chi < 8\pi$）で厳密に微視的に導出する。$N^{-\alpha}$（$\alpha \in (0,1/2)$）のカットオフを用いて、粒子系の経験測度が平均場ケラー-セゲルPDEの解に収束することを示す定量的推定を確立し、混沌の伝播と有界解のグローバル存在を証明する。

ABSTRACT

We derive the two-dimensional Keller-Segel equation from a stochastic system of $N$ interacting particles in the case of sub-critical chemosensitivity $χ&lt; 8 π$. The Coulomb interaction force is regularised with a cutoff of size $N^{- α}$, with arbitrary $α\in (0, 1 / 2)$. In particular we obtain a quantitative result for the maximal distance between the real and mean-field $N$-particle trajectories.

研究の動機と目的

• 2次元で、確率的多数粒子系から決定的ケラー-セゲルPDEを厳密に導出すること。
• 臨界的でない領域$\chi < 8\pi$における平均場極限の妥当性を確立すること。
• $N^{-\alpha}$カットオフ（$\alpha \in (0,1/2)$）を用いた正則化粒子系における混沌の伝播を証明すること。
• 実際の粒子軌道と平均場軌道の間の距離に対する定量的推定を提供すること。
• 従来の研究と比較して、初期密度に関する仮定を弱める。特に、コンパクトな台や高次の正則性を仮定せず、i.i.d. 初期配置を許容すること。

提案手法

• ブラウン運動と正則化されたクーロン力$ k_N(x) = \frac{x}{|x|(|x| + N^{-\alpha})} $に駆動される$N$個の相互作用粒子で細胞集団をモデル化する。
• 実際の粒子と同一のブラウンノイズを用い、有効な漂流$-\chi(k * \rho_t)(Y_t^{i(N)})$に従う平均場軌道$Y_t^{i(N)}$を定義する。
• 実際の粒子軌道$X_t^{i(N)}$が、時間経過にわたって$L^2$型の定量的ノルムで平均場軌道$Y_t^{i(N)}$に近いことを証明する。
• 事前推定、ソボレフ埋没、補間不等式の組み合わせを用いて、ポテンシャル$w = -\log|\cdot| * \rho$およびその微分の正則性を制御する。
• $\rho^N$の$W^{1,p}$ノルムおよび$[\rho^N]_{0,\alpha}$のホルダー半ノルムに対する一様なバウンディングを確立し、力場の一様リプシッツ連続性を保証する。
• タイトネスと収束の議論を適用し、$N \to \infty$の極限において、粒子系の経験測度がケラー-セゲルPDEの解に収束することを示す。

実験結果

リサーチクエスチョン

• RQ1 臨界的でない領域$\chi < 8\pi$において、ケラー-セゲル方程式は確率的粒子系から導出可能か？
• RQ2 $N^{-\alpha}$正則化相互作用力の導入が、平均場極限への収束に与える影響は何か？
• RQ3 実際の粒子軌道と平均場軌道の間の距離に対する定量的バウンディングは、どのように確立できるか？
• RQ4 特にコンパクトな台や滑らかさに関する仮定を除いて、初期密度$\rho_0$に関する仮定を従来の研究より弱めることは可能か？
• RQ5 グリッドベースやコンパクトに台を持つ仮定なしに、i.i.d. 初期粒子配置に対しても混沌の伝播が成立するか？

主な発見

• 著者らは、$\chi < 8\pi$の下で、$N$粒子系の経験測度が$N \to \infty$の極限においてケラー-セゲルPDEの解に収束することを証明した。
• 実際の軌道と平均場軌道の間の最大距離に対する定量的バウンディングが確立された：$\alpha \in (0,1/2)$に対して、$\sup_{t \in [0,T]} \mathbb{E}[\|X_t^{i(N)} - Y_t^{i(N)}\|] \leq C N^{-\alpha/2}$。
• ケラー-セゲル方程式の解$\rho$は、$L^\infty(0,T; C^{0,\alpha}(\mathbb{R}^2))$で有界であり、一様なホルダー半ノルムは初期データのノルムによって制御される。
• ポテンシャル$w = -\log|\cdot| * \rho$の2階微分に対する一様バウンディングにより、力場$-\nabla (k * \rho)$が空間的・時間的に一様にリプシッツ連続であることが示された。
• 初期密度$\rho_0$は、$L^1 \cap L^\infty \cap H^1$に属し、$W^{1,p}$ノルムが有限であれば十分であり、コンパクトな台や$H^4$正則性は不要である。従来の仮定を改善した。
• i.i.d. 初期粒子配置に対しても混沌の伝播が確立され、従来の研究で用いられたグリッドベースや決定的初期条件にとらわれない拡張がなされた。