Robustness of Turing models and gene regulatory networks with a sweet spot

Roozbeh H. Pazuki,Robert G. Endres
DOI: https://doi.org/10.1103/physreve.109.064305
IF: 2.707
2024-06-12
Physical Review E
Abstract:Traditional linear stability analysis based on matrix diagonalization is a computationally intensive process for high-dimensional systems of differential equations, posing substantial limitations for the exploration of Turing systems of pattern formation where an additional wave-number parameter needs to be investigated. In this paper, we introduce an efficient and intuitive technique that leverages Gershgorin's theorem to determine upper limits on regions of parameter space and the wave number beyond which Turing instabilities cannot occur. This method offers a streamlined avenue for exploring the phase diagrams of other complex multi-parametric models, such as those found in gene regulatory networks in systems biology. Due to its suitability for the asymptotic limit of infinitely large systems, it predicts the existence of a sweet spot in network size for maximal Jacobian stability. https://doi.org/10.1103/PhysRevE.109.064305 Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI. Published by the American Physical Society
physics, fluids & plasmas, mathematical
What problem does this paper attempt to address?