Two Nested Limit Cycles in Two-Species Reactions

We search for limit cycles in the dynamical model of two-species chemical reactions that contain seven reaction rate coefficients as parameters and at least one third-order reaction step, that is, the induced kinetic differential equation of the reaction is a planar cubic differential system. Symbolic calculations were carried out using the Mathematica computer algebra system, and it was also used for the numerical verifications to show the following facts: the kinetic differential equations of these reactions each have two limit cycles surrounding the stationary point of focus type in the positive quadrant. In the case of Model 1, the outer limit cycle is stable and the inner one is unstable, which appears in a supercritical Hopf bifurcation. Moreover, the oscillations in a neighborhood of the outer limit cycle are slow-fast oscillations. In the case of Model 2, the outer limit cycle is unstable and the inner one is stable. With another set of parameters, the outer limit cycle can be made stable and the inner one unstable.