Dynamics of a generalized Rayleigh system

Consider the first order differential system given by \begin{equation*} \begin{array}{l} \dot{x}= y, \qquad \dot{y}= -x+a(1-y^{2n})y, \end{array} \end{equation*} where $a$ is a real parameter and the dots denote derivatives with respect to the time $t$. Such system is known as the generalized Rayleigh system and it appears, for instance, in the modeling of diabetic chemical processes through a constant area duct, where the effect of adding or rejecting heat is considered. In this paper we characterize the global dynamics of this generalized Rayleigh system. In particular we prove the existence of a unique limit cycle when the parameter $a\ne 0$.