Approximating Parameter Regions for Multistationarity and Multistability in Mass-Action Systems
Abstract
Identifying kinetic parameters and conserved quantities that support multistationarity and multistability remains challenging in chemical reaction network analysis. Structural theory provides existence and exclusion results, but does not always yield explicit, interpretable descriptions of the corresponding parameter regimes. We present a computational and machine-learning framework for discovering and characterizing parameter regimes with numerical evidence of multistationarity and multistability in mass-action systems. We refer to these numerically detected behaviors as approximate multistationarity and approximate multistability. The framework combines Latin hypercube sampling, adaptive sampling, and interpretable surrogate modeling. Logistic-regression and neural-network surrogates are trained separately on accumulated detection labels and independently guide later sampling toward regions with high scores for approximate multistationarity or multistability. Boundary refinement targets transitions between detected and undetected numerical classes. Buckingham--Pi groups are then used with CN2 rule induction to obtain interpretable empirical descriptions of detected regimes. Illustrative networks yield rules consistent with known multistationarity conditions, while the Schlögl model demonstrates adaptive enrichment of multistationarity detections. In the hybrid histidine kinase model and the sequential dual phosphorylation cycle, adaptive sampling substantially increased the number of parameter - total pairs with detected multistability relative to global sampling under equal sampling budgets. The framework complements analytical chemical reaction network theory by providing efficient, parameter-resolved discovery and characterization of multistationary and multistable regimes.
Related
- Double-layer curtain canopy systems modulate volatile compounds and sensory attributes of ‘Syrah’ grapes and wines in hot-arid regions
- Parameter optimization of IPT systems based on genetic algorithms
- Diffusion learning reveals viable parameter manifolds and compensation geometry in biological dynamical systems
- An optimized machine learning approach for reliable agronomic parameter prediction in precision farming systems
- Neuromodulation-inspired gated associative memory networks: Extended memory retrieval and emergent multistability
BibTeX
@misc{talabis2026approximating,
title = {{Approximating Parameter Regions for Multistationarity and Multistability in Mass-Action Systems}},
author = {Dylan Antonio S. J. Talabis and Victoria May P. Mendoza and Bryan S. Hernandez},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.05242},
eprint = {2610.05242},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.05242},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Approximating Parameter Regions for Multistationarity and Multistability in Mass-Action Systems
AU - Dylan Antonio S. J. Talabis
AU - Victoria May P. Mendoza
AU - Bryan S. Hernandez
PY - 2026
DA - 2026/10/04
PB - arXiv
DO - 10.48550/arxiv.2610.05242
UR - https://arxiv.org/abs/2610.05242
AB - Identifying kinetic parameters and conserved quantities that support multistationarity and multistability remains challenging in chemical reaction network analysis. Structural theory provides existence and exclusion results, but does not always yield explicit, interpretable descriptions of the corresponding parameter regimes. We present a computational and machine-learning framework for discovering and characterizing parameter regimes with numerical evidence of multistationarity and multistability in mass-action systems. We refer to these numerically detected behaviors as approximate multistationarity and approximate multistability. The framework combines Latin hypercube sampling, adaptive sampling, and interpretable surrogate modeling. Logistic-regression and neural-network surrogates are trained separately on accumulated detection labels and independently guide later sampling toward regions with high scores for approximate multistationarity or multistability. Boundary refinement targets transitions between detected and undetected numerical classes. Buckingham--Pi groups are then used with CN2 rule induction to obtain interpretable empirical descriptions of detected regimes. Illustrative networks yield rules consistent with known multistationarity conditions, while the Schlögl model demonstrates adaptive enrichment of multistationarity detections. In the hybrid histidine kinase model and the sequential dual phosphorylation cycle, adaptive sampling substantially increased the number of parameter - total pairs with detected multistability relative to global sampling under equal sampling budgets. The framework complements analytical chemical reaction network theory by providing efficient, parameter-resolved discovery and characterization of multistationary and multistable regimes.
ER -