We analyze the stability and optimization of the synchronization process between two coupled self-excited systems modeled by the multi-limit cycles van der Pol oscillators through the case of an enzymatic substrate reaction with ferroelectric behavior in brain waves model. The one-way and two-way couplings synchronization are considered. The stability boundaries and expressions of the synchronization time are obtained using the properties of the Hill equation. Numerical simulations validate and complement the results of analytical investigations.