 Hopf criticality for an asymmetric Lienard center

For $\ddot{q}=-{\omega }^{2}q+g_2{q}^2+g_3{q}^3+(\tau +f_{1}q+f_2{q}^2)\dot{q}+\cdots$, the Hopf bifurcation amplitude equation has cubic coefficient $(f_2+f_1{g}_2/{\omega }^2)/8$ in the leading position amplitude. The quadratic restoring term produces mean shift $g_2{r}^2/(2{\omega }^2)$ and second harmonic $-g_2{r}^2\cos (2\omega t)/(6{\omega }^2)$. Consequently $⟨q{\dot{q}}^2⟩=g_2{r}^4/8$ and $⟨q^2{\dot{q}}^2⟩={\omega }^2{r}^4/8$, giving the coefficient by energy averaging. A negative value gives a supercritical stable limit cycle; a zero value requires higher-order analysis.