[論文レビュー] Complex-balanced equilibria of generalized mass-action systems: Necessary conditions for linear stability

Balázs Boros, Stefan C. Müller+1 more|arXiv (Cornell University)|Jun 28, 2019

Nonlinear Dynamics and Pattern Formation参考文献 24被引用数 18

ひとこと要約

本稿では、一般化質量作用速度系における複雑平衡の線形安定性の必要条件を確立し、線形安定性が一意性を意味すること、およびヤコビ行列のD安定性とD半安定性によって特徴づけられることを示している。サイクル構造のネットワークでは、線形安定性はD安定性と同値であり、弱可逆ネットワークでは、すべてのサイクルのヤコビ行列に対するD半安定性が要求される。古典的質量作用速度系は対角安定であることが示され、これにより線形安定性がロバストに保たれることが保証される。

ABSTRACT

It is well known that, for mass-action systems, complex-balanced equilibria are asymptotically stable. For generalized mass-action systems, even if there exists a unique complex-balanced equilibrium (in every stoichiometric class and for all rate constants), it need not be stable. We first discuss several notions of matrix stability (on a linear subspace) such as D-stability and diagonal stability, and then we apply our abstract results to complex-balanced equilibria of generalized mass-action systems. In particular, we show that linear stability (on the stoichiometric subspace and for all rate constants) implies uniqueness. For cyclic networks, we characterize linear stability (in terms of D-stability of the Jacobian matrix); and for weakly reversible networks, we give necessary conditions for linear stability (in terms of D-semistability of the Jacobian matrices of all cycles in the network). Moreover, we show that, for classical mass-action systems, complex-balanced equilibria are not just asymptotically stable, but even diagonally stable (and hence linearly stable). Finally, we recall and extend characterizations of D-stability and diagonal stability for matrices of dimension up to three, and we illustrate our results by examples of irreversible cycles (of dimension up to three) and of reversible chains and S-systems (of arbitrary dimension).

研究の動機と目的

• 一般化質量作用速度系における複雑平衡の線形安定性の必要条件を同定すること。
• 一般化動力学における漸近的安定性と線形安定性の違いを明確にすること。
• ストォイキオメトリック部分空間上で、D安定性や対角安定性といった行列安定性の概念を用いて線形安定性を特徴づけること。
• 線形安定性が複雑平衡の一意性を意味することを確立すること。
• 古典的質量作用速度系に関する既知の結果を拡張し、複雑平衡が対角安定であることを証明すること。

提案手法

• ストォイキオメトリック行列と動力学順序行列を用いた常微分方程式による一般化質量作用速度系の分析。
• 複雑平衡を反応ネットワークのラプラシアン行列Akを用いて定義し、Ak x^Y = 0 の解として表す。
• 線形化を用いて安定性を検討し、ヤコビ行列 J(x) = Y Ak diag(x^Y) eY^T diag(x^{-1}) に注目する。
• ストォイキオメトリック部分空間S上で、D安定性、対角安定性、D半安定性といった行列安定性の概念を適用する。
• 3次元までの行列についてのD安定性および対角安定性の特徴付け（トレースや行列式の条件を含む）を用いる。
• 具体的なネットワークタイプ（不可逆サイクル、可逆鎖、Sシステム）に応用し、行列不等式から明示的な安定性基準を導出する。

実験結果

リサーチクエスチョン

• RQ1 すべての速度定数に対して、一般化質量作用速度系における複雑平衡が線形安定となる条件は何か？
• RQ2 一般化質量作用速度系における線形安定性は、複雑平衡の一意性を意味するか？
• RQ3 サイクル構造のネットワークにおいて、複雑平衡の線形安定性を特徴づける行列的性質（例：D安定性）は何か？
• RQ4 弱可逆ネットワークにおいて、線形安定性を保証するサイクルヤコビ行列に必要な条件は何か？
• RQ5 複雑平衡を有する古典的質量作用速度系は対角安定であるか？その結果、ロバスト性にどのような意味があるか？

主な発見

• 一般化質量作用速度系における複雑平衡の線形安定性は、すべてのストォイキオメトリッククラスにおいてその平衡の一意性を意味する。
• サイクル構造のネットワークでは、線形安定性はストォイキオメトリック部分空間上でのヤコビ行列のD安定性と同値である。
• 弱可逆ネットワークでは、線形安定性はネットワーク内すべてのサイクルについて、そのヤコビ行列がD半安定であることと同値である。
• 古典的質量作用速度系では、複雑平衡が対角安定であることが示され、これにより線形安定性および動力学順序の微小な摂動に対するロバスト性が保証される。
• 3次元の不可逆サイクルでは、線形安定性が成り立つための必要十分条件は、ヤコビ行列の行列式が正であること、かつすべての対角成分が非正であり、かつ少なくとも1つが厳密に負であることである。
• Sシステムおよび可逆鎖では、任意の可逆反応の動力学順序が対応する物質の化学计量係数を超える場合、線形不安定性を検出できる。