Efficient Methods to Compute Hopf Bifurcations in Chemical Reaction Networks Using Reaction Coordinates
Résumé
We build on our previous work to compute Hopf bifurcation fixed points for chemical reaction systems on the basis of reaction coordinates. For determining the existence of Hopf bifurcations the main algorithmic problem is to determine whether a single multivariate polynomial has a zero for positive coordinates. For this purpose we provide heuristics on the basis of the Newton polytope that ensure the existence of positive and negative values of the polynomial for positive coordinates. We apply our method to the example of the Methylene Blue Oscillator (MBO).