Implications of network topology on stability.
Bottom Line:
Visualization of the Hurwitz determinants over the reduced parameters defines the network's stability phase space, delimiting the range of its dynamics (specifically, the possible numbers of unstable roots at each steady state solution).Any further explicit algebraic specification of the network will project onto this stability phase space.Stability analysis via this hierarchical approach is demonstrated on classical networks from multiple fields.
View Article:
PubMed Central - PubMed
Affiliation: Department of Systemic Cell Biology, Max Planck Institute of Molecular Physiology, Dortmund, Germany.
ABSTRACT
In analogy to chemical reaction networks, I demonstrate the utility of expressing the governing equations of an arbitrary dynamical system (interaction network) as sums of real functions (generalized reactions) multiplied by real scalars (generalized stoichiometries) for analysis of its stability. The reaction stoichiometries and first derivatives define the network's "influence topology", a signed directed bipartite graph. Parameter reduction of the influence topology permits simplified expression of the principal minors (sums of products of non-overlapping bipartite cycles) and Hurwitz determinants (sums of products of the principal minors or the bipartite cycles directly) for assessing the network's steady state stability. Visualization of the Hurwitz determinants over the reduced parameters defines the network's stability phase space, delimiting the range of its dynamics (specifically, the possible numbers of unstable roots at each steady state solution). Any further explicit algebraic specification of the network will project onto this stability phase space. Stability analysis via this hierarchical approach is demonstrated on classical networks from multiple fields. No MeSH data available. Related in: MedlinePlus |
Related In:
Results -
Collection
License getmorefigures.php?uid=PMC4380337&req=5
Mentions: Sel’kov in 1968 [30] proposed the following simple model to account for glycolytic oscillations:x˙1=−k1x1+σ1k2x12x2+σ2k3x2(122)x˙2=k0−k2x12x2−k3x2,(123)with σ1 and σ2 equal to 1. The generalized version is:x˙1=−V11+σ1V122+σ2V23(124)x˙2=V0−V122−V23.(125)The influence topology is shown in Fig. 10A, with principal minors:b1=−1+σ1ρ1−ρ2−ρ3(126)b2=ρ2+ρ3−σ1ρ1ρ3−σ2ρ1ρ3,(127)and Hurwitz determinants:Δ1=1+ρ2+ρ3−σ1ρ1(128)Δ2=(1+ρ2+ρ3−σ1ρ1)(ρ2+ρ3−σ1ρ1ρ3−σ2ρ1ρ3).(129)This result can also be obtained directly from the cycle-based definitions (where I have already removed terms that are zero):Δ1=−c1(130)Δ2=−c1·c1c1¯+c1·c2.(131)The corresponding stability phase space is displayed in Figs. 10B-D for σ1 = σ2 = 1 and different values of ρ3. |
View Article: PubMed Central - PubMed
Affiliation: Department of Systemic Cell Biology, Max Planck Institute of Molecular Physiology, Dortmund, Germany.
No MeSH data available.