[論文レビュー] On spatially uniform behavior in reaction-diffusion PDE and coupled ODE systems

Murat Arcak|ArXiv.org|Aug 18, 2009

Gene Regulatory Network Analysis参考文献 37被引用数 3

ひとこと要約

本稿では、反応拡散PDEおよび結合されたODE系における空間的一致性のための新しいリャプノフ基準を提示する。従来のグローバルリプシッツ仮定を緩め、反応項のヤコビアンとラプラシアンの2番目のノイマン固有値を含むより緩い不等式に置き換える。主な貢献は、線形計画法（LMIs）を用いた検証可能な条件であり、これは極限円周運動を示す系に対しても、漸近的な空間的一致性を保証する。また、拡散結合されたODEネットワークにおける同期化にも拡張可能である。

ABSTRACT

We present a condition which guarantees spatial uniformity for the asymptotic behavior of the solutions of a reaction-diffusion PDE with Neumann boundary conditions. This condition makes use of the Jacobian matrix of the reaction terms and the second Neumann eigenvalue of the Laplacian operator on the given spatial domain, and replaces the global Lipschitz assumptions commonly used in the literature with a less restrictive Lyapunov inequality. We then present numerical procedures for the verification of this Lyapunov inequality and illustrate them on models of several biochemical reaction networks. Finally, we derive an analog of this PDE result for the synchronization of a network of identical ODE models coupled by diffusion terms.

研究の動機と目的

• ノイマン境界条件を満たす反応拡散PDEの解が空間的に一様な漸近的挙動を示す条件を同定すること。
• 先行研究におけるグローバルリプシッツ仮定を、より緩いリャプノフ不等式条件に置き換えること。
• リャプノフ不等式の数値的検証手順を、線形計画法（LMIs）を用いて開発すること。
• PDEの結果を、拡散項に類する形で結合された同一のODEネットワークにおける同期化に拡張すること。
• 生化学的モデル、特にグッドウィン型およびゴールドベーター型の周期的系を含む、実例への適用を示すこと。

提案手法

• 空間領域上のラプラシアン演算子の2番目のノイマン固有値 $ \lambda_2 $ とヤコビアン行列 $ J(x) $ を含むリャプノフ不等式条件を提案する。
• グローバルリプシッツ境界を避けるために、行列 $ J(x) - \lambda_2 D $ がリャプノフ不等式を満たすことが、空間的一致性の十分条件であることを導出する。
• 2つのLMIsに基づく検証手順を開発する：1つ目は定数ヤコビアン行列の凸包（または凸包）を用いる（定理2）、2つ目は特別な凸集合を用いてLMIsの次元を低減する（定理3）。
• S手続きを用いてリャプノフ不等式を解けるLMIs形式に変換し、数値的検証を可能にする。
• PDEの結果を結合されたODE系に拡張する際、ラプラシアン演算子の固有値を結合グラフのラプラシアン行列の2番目に小さい固有値に置き換える（定理4）。
• ラプラシアン行列の性質がラプラシアン演算子と類似していることを利用し、PDEの場合と同一の条件下で同期化が成立することを証明する。

実験結果

リサーチクエスチョン

• RQ1 ノイマン境界条件を満たす反応拡散PDEが、どのような条件下で空間的に一様な漸近的挙動を示すか？
• RQ2 従来の均一性証明で一般的に使われるグローバルリプシッツ仮定を、空間的一致性を保証する上で緩めることは可能か？
• RQ3 提案されたリャプノフ不等式条件は、現実的な生化学的モデルに対してどのように数値的に検証可能か？
• RQ4 同じ条件は、拡散項に類する形で結合された同一のODEネットワークにおける同期化を保証するか？
• RQ5 固定点吸引子でない系、例えば極限円周運動を示す系に対しても、この手法は適用可能か？

主な発見

• ヤコビアン行列 $ J(x) - \lambda_2 D $ を含むリャプノフ不等式に基づく提案条件は、反応項のグローバルリプシッツ連続性を仮定せず、反応拡散PDEの解の空間的一致性を保証する。
• グローバルリプシッツに基づくアプローチに比べ、本手法は著しく優れた推定値を達成し、例2では減衰率推定値が桁違いに改善されている。
• ヤコビアン行列が定数行列の凸包に属する系では、検証条件が線形計画法（LMI）に簡略化され、数値的に解ける。
• ヤコビアンが特別な凸集合（例：ランク1構造）に属する場合、LMIsの次元が低減され、例2のグッドウィン型オシレーターで示されるように、解析的可解性テストが可能になる。
• 同一のリャプノフ条件は、拡散結合された同一のODEネットワークにおける同期化を保証し、$ \lambda_2 $ が結合グラフのラプラシアン行列の2番目に小さい固有値に置き換えられる。
• 数値的検証により、生物時計のゴールドベーター・モデルに対してもLMIs条件の実行可能性が確認され、複雑な生物学的オシレーターへの適用可能性が示された。