[论文解读] 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 x^Y = 0 的解，其中 Ak 是反应网络的拉普拉斯矩阵。
• 通过线性化研究稳定性，重点关注雅可比矩阵 J(x) = Y Ak diag(x^Y) eY^T diag(x^{-1})。
• 在化学计量子空间 S 上应用矩阵稳定性概念——D-稳定性、对角稳定性与D-半稳定性。
• 利用至多三维矩阵的D-稳定性和对角稳定性刻画，包括迹与行列式条件。
• 将结果应用于特定网络类型：不可逆环、可逆链与S-系统，从矩阵不等式中推导出显式稳定性准则。

实验结果

研究问题

• RQ1 在广义质量作用系统中，对于所有速率常数，复杂平衡态在何种条件下是线性稳定的？
• RQ2 广义质量作用系统中的线性稳定性是否蕴含复杂平衡态的唯一性？
• RQ3 对于循环网络，何种矩阵性质（如D-稳定性）刻画了复杂平衡态的线性稳定性？
• RQ4 对于弱可逆网络，何种关于环路雅可比矩阵的必要条件可确保线性稳定性？
• RQ5 经典质量作用系统若具有复杂平衡态，是否为对角稳定？这对鲁棒性有何含义？

主要发现

• 广义质量作用系统中复杂平衡态的线性稳定性，意味着在每个化学计量类中平衡态是唯一的。
• 对于循环网络，线性稳定性等价于雅可比矩阵在化学计量子空间上的D-稳定性。
• 对于弱可逆网络，线性稳定性蕴含网络中每个环路的雅可比矩阵均为D-半稳定。
• 在经典质量作用系统中，复杂平衡态是对角稳定的，这蕴含线性稳定性和对质量作用阶微小扰动的鲁棒性。
• 对于三维不可逆环，线性稳定性成立当且仅当雅可比矩阵的行列式为正，且所有对角元非正，且至少有一个严格为负。
• 在S-系统和可逆链中，若任意可逆反应的速率阶超过对应物种的化学计量系数，则可检测到线性不稳定。