[論文レビュー] Karlin Theory On Growth and Mixing Extended to Linear Differential Equations

Lee Altenberg|arXiv (Cornell University)|Jun 16, 2010

Mathematical Biology Tumor Growth参考文献 15被引用数 5

ひとこと要約

この論文は、離散時間系における成長と混合に関するKarlinの1982年証明を連続時間線形微分方程式へと拡張し、行列 𝐷 + 𝑚𝐴 の固有値の実部（スペクトル絶対値）が混合強度 𝑚 の増加とともに減少することを証明する。この結果は、サイト固有の成長/減衰率を表す非一様な対角行列 𝐷 と、混合を表す非負で非可約な非対角行列 𝐴 に対して成り立つ。減少のメカニズムは、行列 𝐴 のスペクトル絶対値に関連するペロンベクトルの不等式によって駆動される。

ABSTRACT

Karlin's (1982) Theorem 5.2 shows that linear systems alternating between growth and mixing phases have lower asymptotic growth with greater mixing. Here this result is extended to linear differential equations that combine site-specific growth or decay rates, and mixing between sites, showing that the spectral abscissa of a matrix D + m A decreases with m, where D does-not-equal c I is a real diagonal matrix, A is an irreducible matrix with non-negative off-diagonal elements (an ML- or essentially non-negative matrix), and m >= 0. The result is based on the inequality: u' A v < r(A), where u and v are the left and right Perron vectors of the matrix D + A, and r(A) is the spectral abscissa and Perron root of A. The result gives an analytic solution to prior work that relied on two-site or numerical simulation of models of growth and mixing, such as source and sink ecological models, or multiple tissue compartment models of microbe growth. The result has applications to the Lyapunov stability of perturbations in nonlinear systems.

研究の動機と目的

• 離散時間系におけるKarlinの1982年結果（成長と混合）を連続時間線形微分方程式へと拡張すること。
• 非一様な成長率に対して、行列 𝐷 + 𝑚𝐴 のスペクトル絶対値が混合パラメータ 𝑚 の増加とともに減少することを確立すること。
• 生態学的・ウイルス学的・コンpartimentalダイナミクスのモデルにおいて、従来は数値的または2サイトシミュレーションに依存していた問題に対する解析的解法を提供すること。
• サイト固有の成長とサイト間混合を伴うシステムへの結果の一般化可能性を示すこと。特に非線形系の安定性解析への応用を含む。

提案手法

• 離散時間写像から線形1階微分方程式へのKarlinの定理5.2の形式的拡張。
• 連続時間系における長期的成長率の主要指標としてスペクトル絶対値の使用。
• 非負の非対角要素を持つ非可約行列（ML行列）に対するペロン＝フロベニウス理論の適用。
• 𝐮ᵀ𝐴𝐯 < 𝑟(𝐴) の不等式の導出。ここで 𝐮 と 𝐯 は行列 𝐷 + 𝐴 の左および右ペロンベクトルであり、𝑟(𝐴) は行列 𝐴 のスペクトル絶対値である。
• 𝐷 ≠ 𝑐𝐼 の場合に、行列 𝐷 + 𝑚𝐴 のスペクトル絶対値が 𝑚 に関して厳密に減少することの証明。凸性および厳密不等式条件に基づく。
• ウイルスダイナミクス、生態学的ソース・サングシステム、および非均一な成長と混合を有するその他のコンpartimentalシステムへの結果の応用。

実験結果

リサーチクエスチョン

• RQ1 非一様な成長率を持つ連続時間系における線形微分方程式において、混合の増加に伴い成長が減少するKarlinの結果は成り立つか？
• RQ2 サイト固有の成長（𝐷）と混合（𝑚𝐴）を組み合わせた行列のスペクトル絶対値が、混合強度 𝑚 の増加に伴い減少することを示せるか？
• RQ3 成長率がサイトごとに非一様である場合、スペクトル絶対値の厳密減少を保証する条件は何か？
• RQ4 この結果は、ウイルス学・生態学・集団動態学における離散時間モデルを超えたシステムへどのように一般化できるか？
• RQ5 数値シミュレーションや2サイト近似に依存せずに、混合に伴う成長率の解析的減少を確立できるか？

主な発見

• 行列 𝐷 + 𝑚𝐴 のスペクトル絶対値は、𝐷 が単位行列のスカラー倍でない限り、混合強度 𝑚 の増加に伴い厳密に減少する。
• 減少は不等式 𝐮ᵀ𝐴𝐯 < 𝑟(𝐴) によって駆動され、ここで 𝐮 と 𝐯 は行列 𝐷 + 𝐴 の左および右ペロンベクトルであり、𝑟(𝐴) は行列 𝐴 のスペクトル絶対値である。
• 従来は数値的シミュレーションや2コンpartiment近似に依存していたウイルス学および生態学分野の問題に対して、解析的解法を提供する。
• 混合レベルにかかわらず、システムの長期的成長率は、最大のコンpartiment固有成長率によって上から抑えられる。
• 混合強度が無限大に近づくにつれ、システムの成長率はコンpartiment成長率の重み付き平均に収束する。
• 成長率が非一様である場合（すなわち 𝐷 ≠ 𝑐𝐼 の場合）、混合に伴う成長率の低下は 𝑚 に関して厳密に凸である。