Estimating the boundaries of a limit cycle in a 2D dynamical system using renormalization group

While the plausibility of formation of limit cycle has been a well studied topic in context of the Poincare-Bendixson theorem, studies on estimates in regard to the possible size and shape of the limit cycle seem to be scanty in the literature. In this paper we present a pedagogical study of some aspects of the size of this limit cycle using perturbative renormalization group by doing detailed and explicit calculations upto second order for the Selkov model for glycolytic oscillations. This famous model is well known to lead to a limit cycle for certain ranges of values of the parameters involved in the problem. Within the tenets of the approximations made, reasonable agreement with the numerical plots can be achieved.