the interaction among oscillators in a system not only leads them into a cooperative dynamics but also often quenches their oscillations .
there are two dynamically different oscillation quenching phenomena which are termed as amplitude death ( ad ) and oscillation death ( od ) @xcite . in ad , the amplitude of oscillation quenches to zero , whereas od is caused by quenching in the frequency of oscillation @xcite .
it is also defined that ad occurs via the stabilization of a homogeneous steady state ( hss ) while od occurs via the stabilization of an inhomogeneous steady state @xcite .
the mechanisms underlying these two quenching phenomena have been identified recently and the results show that the parametric mismatch @xcite , dynamic @xcite , time delay @xcite and nonlinear couplings @xcite underly ad , while od occurs mainly because of the symmetry breaking coupling in the system @xcite . recently ,
diverse routes of transition from the ad to od have also been reported @xcite .
experimental and theoretical studies show definitive evidence of oscillation quenching in realistic systems ranging from biological @xcite , chemical @xcite , electronic @xcite and laser @xcite systems to climate @xcite systems .
such oscillation quenching in many cases leads to undesirable situations . in the interaction between neuronal dynamics and brain metabolism , the decrease in cerebral metabolic rate ,
coupled with the stabilizing properties of atp - gated potassium channels , leads to a burst suppression in the eeg pattern which symbolizes inactivated brain @xcite .
this type of suppression results in hypothermia , coma and ohtahara syndrome , a type of early infantile encephalopathy and is also observed during deep levels of anesthesia .
the suppression of normal sinus rhythm of pacemaker cells causes cardiac arrest @xcite .
owing to the fatal consequences due to oscillation suppression , interesting efforts have been undertaken to retrieve / resurrect oscillations of the system @xcite . in @xcite , the oscillation death in diffusively
coupled oscillators has been found to be eliminated through a spatial disorder in the form of parametric mismatch and in ref .
@xcite processing delay is used to revoke oscillations successfully in delay coupled systems .
regarding the above mentioned issues , in this article , we demonstrate that the problem can also be well resolved by providing a suitable feedback in the system .
the latter can be found to be present in most of the natural systems , including neural networks @xcite , genetic networks @xcite , vision systems @xcite , etc .
the vital role of feedback in controlling the dynamics of the given system and the control over synchronization @xcite are already known , which can be seen in a variety of fields ranging from electronics @xcite , biology to quantum information @xcite .
for example , the feedback control of deep brain simulation has been found to be the most effective treatment for chronic neural diseases like essential tremor , dystonia and parkinson s disease @xcite .
also , the feedback generated by the voltage - gated ion channels in neural cells is found to be crucial in generating neural signals @xcite .    in this article
, we show the applicability of the feedback technique in resurrecting oscillations in a wide range of systems .
we show both numerically and analytically that the addition of feedback destabilizes the stable attractors which results in a wiping out of the oscillation quenching and inducing a resurrection of oscillations .
in addition , by considering a rather general system , we prove analytically the above destabilizing nature over the ad state .
further , it will be more important to develop an adaptable mechanism thereby improving the ones available in the literature at present .
this is because , for example to use the available parametric mismatch method , one needs to tune the internal parameters of the system , while the processing delay also depends on the underlying process of the system where that process may be unknown in many situations so as to hinder the efficiency .
in contrast , the feedback method suggested here can be given more easily which is already in practice under different contexts such as in deep brain simulation @xcite .
in addition to the above adaptable nature of feedback , with the aid of numerical and analytical studies , we show the important fact that this method does not impose a restriction that the output of all the oscillators need to be fedback . from the output of only a few of the system oscillators ,
we show that in typical systems the resurrection of oscillations can be achieved easily .
the structure of the paper is as follows . in section [ genr ]
we present the general form of the system that we consider . in section [ diffu ]
, we illustrate the role of feedback on two important oscillation quenching scenarios , namely the symmetry breaking coupling and the parametric mismatch in a system of diffusively coupled stuart - landau oscillators , through numerical analysis . in section
[ anal ] , we present suitable analytical support of the numerical results based on an appropriate linear stability analysis . in section [ dynam ] , we illustrate the role of the considered feedback in indirectly coupled or dynamically coupled stuart - landau oscillators . a realistic chemical oscillator model , namely the brusselator model , is considered in section [ brusec ] .
we have also proved the applicability of the technique in more general situations over amplitude death state in appendix [ gen_an ] .
in addition , in the appendix [ more ] we illustrate our method with different coupling schemes and with different models such as van der pol oscillator , rsseler system and so on .
appendix [ app2_an ] includes the details in obtaining the boundary curves of the ad region which were given in sec .
a summary of our results and conclusions are presented in section [ conclu ] .
consider a system of coupled dynamical systems , @xmath0 where @xmath1 characterizes the dynamics of the isolated @xmath2-th system , @xmath3 is a @xmath4 dimensional state vector of the system @xmath2 , @xmath5 is the @xmath6 coupling matrix of the network , @xmath7 are respectively the uniform coupling and feedback strengths , @xmath8 is a coupling function and @xmath9 is the feedback term which can be written as @xmath10 . here @xmath11 is simply a @xmath12 constant matrix and @xmath13 characterizes the feedback and it depends on the state vectors of the system .
such a dependence of @xmath14 on the state vectors of the system may be linear ( example : @xmath15 ) or nonlinear ( example : @xmath16 ) , where @xmath17 s represent weight factors which can take values from @xmath18 to @xmath19 and @xmath20 is a suitable number . in our following study
, we consider the form of @xmath14 as @xmath15 .    in the appendix [ gen_an ] , we have considered rather general forms for @xmath21 and @xmath22 and shown analytically that the trivial ad state which appears in the system could be wiped out through the strengthening of @xmath23 so as to resurrect oscillations .    ) with respect to the introduced two different forms of feedback ( i ) @xmath24 and ( ii ) @xmath25 , @xmath26 , for the choice @xmath27 and @xmath28 .
( b)for two choices of @xmath29 , @xmath30 and @xmath31 , the change in the hb point by the increment of @xmath32 ( which appears in the expression of @xmath17 ) with @xmath33 in ( [ sys1 ] ) .
inset in ( b ) projects the od and oscillatory ( os ) regions in the ( @xmath34 ) space for @xmath35 .
( c),(d ) and ( e ) : temporal behaviors of the system in the os , od and revoked os state respectively for the values of ( @xmath34 ) @xmath36 ( @xmath37 ) , ( @xmath38 ) , and ( @xmath39 ) . the color bar in figs .
( c ) , ( d ) and ( e ) gives information on the value of @xmath40 corresponding to different colors.,width=362 ]
to start with , we use the paradigmatic model known as the coupled stuart - landau oscillators for the purpose of illustration for the advocated feedback method .
it is well known that the dynamical equation defining the stuart - landau oscillator can be obtained from a general ordinary differential equation near a hopf bifurcation point @xcite . as
the hopf bifurcation arises widely in the literature , the stuart - landau oscillator helps to model a variety of systems in different areas ranging from biology @xcite to lasers @xcite and is also used in the reaction - diffusion process @xcite .
this model is often used in neural networks to model spiking neurons @xcite .
the first reason for using the model in neural networks is that the periodically spiking neurons have an exponentially stable limit cycle attractor and secondly the real part of the complex amplitude of the stuart - landau oscillator can describe the membrane voltage in the neurons and the imaginary part can be related to the recovery variable embedding the effects of the other variables of physiological neuron models @xcite .
thus , in the literature we can find the use of this model in studying the effects of synchronization and desynchronization in neural networks @xcite and also various collective dynamical states , including chimeras @xcite . in our study
, we also include other useful models such as the van der pol oscillator , rssler system and brusselator model .
the corresponding results are briefly indicated in the appendix b and sec [ brusec ] .
now , we first consider a system of coupled stuart - landau oscillators which is characterized by    @xmath41    where @xmath42 , @xmath43^{t}$ ] , @xmath44 is the kronecker delta ( @xmath44 @xmath45 , if @xmath46 and @xmath47 , if @xmath48 ) and * i * represents the identity matrix , in eq .
( [ sys1 ] ) . here
the diffusive coupling acts only on the first half of the evolution equations ( @xmath49- variable alone ) which breaks the rotational symmetry @xcite and consequently induces oscillation death in the system .    to elucidate clearly the role of feedback in ( [ sys1 ] ) ,
we first consider the case @xmath27 , with @xmath50 .
in the case @xmath51 , @xmath52 , we note that the individual systems in ( [ sys1 ] ) show limit cycle oscillations with @xmath53 .
the introduction and strengthening of diffusive coupling ( @xmath54 ) stabilizes the symmetric pair of nontrivial equilibrium points @xmath55 @xmath56 , where @xmath57 , @xmath58 , @xmath59 and @xmath60 , through a sub - critical hopf bifurcation @xcite which gives rise to od . with the introduction of feedback @xmath61 ,
these nontrivial equilibrium points soon lose their stability via a super - critical hopf bifurcation . in our study
, we introduced such a feedback in two ways : ( i ) @xmath62 and ( ii ) @xmath63 , @xmath26 ( again here @xmath64 denotes the kronecker delta ) .
first by setting @xmath62 , for different values of @xmath29 we traced the hopf bifurcation points and these points are collectively shown as the hb curve@xmath65 in fig .
[ am1](a ) .
the region lying under this curve is an od region ( denoted by od@xmath65 ) and the parametric region above this curve is free from od and corresponds to oscillatory states ( os ) .
similarly for @xmath63 , @xmath26 , the od region ( od-@xmath66 which includes od-@xmath19 also ) and the curve of hopf bifurcation points ( hb curve@xmath67 ) have been shown in fig.[am1](a ) , which show that the uniform distribution of @xmath17 as in ( i ) helps to redeem from the od state sooner than in the case ( ii ) .    by extending the constituents of the network to @xmath33
, we have verified that this technique can work as well with larger @xmath68 .
further , we have checked whether the feedback needs contributions from all the constituents of the network .
this is vital as in a practical situation we can not assure or impose all the constituents to contribute to the feedback .
thus , we have distributed @xmath17 s as @xmath69 , @xmath70 , where @xmath32 determines the number of contributing components of the network . by varying @xmath32
, we have drawn the hb curves separating the od state with the os state for two different values of @xmath29 , @xmath30 and @xmath31 in fig .
[ am1](b ) .
interestingly , these curves demonstrate clearly that the contribution from even a single oscillator is sufficient to revoke oscillations in the network .
secondly , the critical value of @xmath23 ( hb point ) above which the os state arises gets decreased sharply with that of @xmath32 .
these facts prove that this technique can work well regardless of the number of oscillators present in the network and the number of them which contributes towards the feedback . for simplicity , we chose @xmath35 in the following studies .     in the system ( [ sys1 ] ) endowed with a parametric mismatch of @xmath71 , @xmath72 and @xmath27 .
( c ) the reduction of ad regions in the ( @xmath73 ) space with the increase of @xmath23 and by keeping @xmath74 and @xmath27 .
( d ) and ( e ) the emergence of os state from ad and od states in @xmath33 case for @xmath71 and @xmath72 which is obtained for two different set of initial conditions .
ad region in fig .
( e ) is also in conformity with the analytical results presented in the sec .
[ anal].,width=362 ]    the inset of fig .
[ am1](b ) depicts the information about od and os regions with @xmath33 oscillators . correspondingly , in figs .
[ am1](c ) , [ am1](d ) and [ am1](e ) , we have captured the temporal behaviors of the system ( in the @xmath75-variables ) for different sets of ( @xmath76 ) , for a finite time interval after leaving out sufficiently large transients .
the first one ( fig .
[ am1](c ) ) shows the temporal behavior at an os state of the system when @xmath52 and @xmath77 , where the value of @xmath29 is not sufficient to induce od .
now increasing @xmath29 to @xmath78 ( while keeping @xmath52 ) , the subsequent figure ( fig .
[ am1](d ) ) shows the quenching of this oscillation . now switching @xmath23 on , the temporal behavior in fig .
[ am1](e ) shows the resurrected oscillations for @xmath79 .
rubchinsky and sushchik @xcite have introduced a disorder in the form of parametric mismatch which revokes oscillations in ( [ sys1 ] ) . on the other hand , just like symmetry breaking which is predominant in inducing od , the parametric mismatch is also a key candidate that induces ad in the system .
recently , koseska _ et al .
_ in @xcite have shown that an increase in this inhomogeneity not only induces ad but also od , whereas the feedback mechanism that we consider here does not show such a behavior in the system .
this feature provides a definitive advantage over parametric mismatch .
now , we augment the system with a parametric mismatch . for the mismatch in the parameter @xmath80 , @xmath81 ,
in the above coupled stuart - landau oscillators ( [ gen])-([sys1 ] ) with @xmath27 , we find the existence of both ad ( @xmath82 ) and od ( @xmath30 ) while @xmath52 . now ,
from these death states the transitions towards os state by @xmath23 are demonstrated in figs .
[ bif1](a ) and [ bif1](b ) , which show the destabilization of both the ad and od states via supercritical hopf bifurcations .
further , the role of @xmath23 over @xmath83 and @xmath29 is more clearly demonstrated in fig .
[ bif1](c ) , where the colored islands denote the ad regions for the values of @xmath23 @xmath84 and @xmath85 .
one can check that similar phenomenon occurs to od regions also which we do not depict here explicitly .    now considering the case of @xmath33 globally coupled oscillators with @xmath86 for @xmath87 and @xmath88 for @xmath89 in ( [ sys1 ] ) , we have demonstrated the reduction in the ad and od regions in figs .
[ bif1](d ) and [ bif1](e ) .
the above two figures are plotted for two different sets of initial conditions . among them , the ad region in fig .
[ bif1](e ) is the analytically relevant region which has been obtained in the next section .    in the next section , the above
obtained numerical results on the system ( [ sys1 ] ) are verified through analytical results wherever possible .
also , we illustrate the applicability of the technique to several situations in appendix [ more ] , where we considered ( i ) repulsive link : stuart - landau ( sl ) oscillator , ( ii ) conjugate coupling : sl oscillator , ( iii ) repulsive link : van der pol oscillator , ( iv ) directly and indirectly coupled rssler system and ( v ) other chaotic oscillators ( sprott and lorenz systems ) .
) space for @xmath27 oscillators in system ( [ sys1 ] ) with @xmath90 . with parametric mismatch case : fig .
( b ) ad regions of system ( [ sys1 ] ) for @xmath91 , @xmath74 and @xmath92 bounded by the curves defined in ( [ condi ] ) . fig .
( c ) analytically obtained boundaries of ad regions in ( @xmath83 , @xmath29 ) for different values of @xmath23 ( = @xmath93 ) , where @xmath92 .
( d ) hb curves for the three different cases ( i ) @xmath92 , ( ii ) @xmath94 and ( iii ) @xmath95 for @xmath91 and @xmath72.,width=366 ]    in this section , we present relevant analytical confirmations of the numerical results corresponding to the system ( [ sys1 ] ) for the cases with and without parametric mismatch discussed earlier in sec .
[ diffu ] .
in addition through the obtained analytical results , we show the effectiveness of the technique with the feedback contribution coming from a few number of oscillators in the network .      in the absence of any parametric mismatch ,
as pointed out earlier the system in ( [ sys1 ] ) has a trivial equilibrium point @xmath96 : @xmath97 and two pairs of non - trivial equilibrium points as @xmath98 where @xmath60 and @xmath99 .
the linear stability of these fixed points is determined by the eigenvalues of the jacobian matrix    @xmath100    where @xmath101 , @xmath102 , @xmath103 . while @xmath52 , the eigenvalues of @xmath104 corresponding to the trivial equilibrium point @xmath96 are @xmath105 the eigenvalues corresponding to the equilibrium points @xmath106 and @xmath107 are given by    @xmath108    similarly , @xmath109 and @xmath110 have the set of eigenvalues    @xmath111    from eqs .
( [ 0eig]-[2eig ] ) we can note that among the five equilibrium points , @xmath106 and @xmath107 are found to have all their eigenvalues satisfying the condition @xmath112<0 $ ] for the parametric range @xmath113 while @xmath114 , @xmath115 while @xmath116 and @xmath117 while @xmath118 and thus they are stable in this range . on the other hand the other equilibrium points @xmath96 , @xmath119 and @xmath120 can never become stable for any choice of parametric values .
thus the stabilization of @xmath121 and @xmath122 essentially gives rise to oscillation death in the system .
now we introduce the feedback in such a way that @xmath123 and @xmath124 of @xmath14 in ( [ sys1 ] ) take the values @xmath125 and @xmath126 . with such a choice , we find that the stability determining eigenvalues corresponding to the equilibrium points change as @xmath127 the above equation shows that an increase in @xmath23 can destabilize the equilibrium points @xmath121 and @xmath122 through a hopf bifurcation . in the cases where @xmath128 , for all values of @xmath29 the hopf bifurcation occurs at @xmath129 whereas in the case of @xmath130 , if @xmath131 , the hopf bifurcation occurs at @xmath132 where @xmath133 is given in ( [ c11 ] ) . if @xmath134 , the hopf bifurcation occurs as in ( [ etac1 ] ) .
the other equilibrium points @xmath96 , @xmath119 and @xmath120 are found to remain unstable .
for the case of @xmath28 , the curve of hopf bifurcation points ( @xmath135 ) separating the death regimes with the oscillatory regimes is shown in fig .
[ anal1](a ) which matches exactly with the one obtained numerically ( hb curve - @xmath19 in fig .
[ am1](a ) ) .
the analytical treatment of the other case corresponding to the hb curve-2 can also be investigated in a similar manner , though the results can not be written down in such a transparent manner .
so we do not present the details here .
next with the introduction of a parametric mismatch in the system ( [ sys1 ] ) , we prove the validity of the technique for the trivial ad state analytically for the @xmath27 case . as mentioned earlier in sec .
[ parame ] , an increase in the value of the parametric mismatch parameter @xmath83 causes the stabilization of the trivial equilibrium point @xmath97 and thus introduces ad in the system . to destabilize the latter ,
we first introduce the feedback in such a way that @xmath136 , where @xmath137 .    as before , through a linear stability analysis of the system ( [ sys1 ] ) with parametric mismatch included , we look for the stable regions of the equilibrium point @xmath97 , which can be studied from the characteristic equation of the linear eigenvalue problem of the system as @xmath138    ) with @xmath33 oscillators in the os , od and revoked os states respectively , for the values of ( @xmath34 ) @xmath36 ( @xmath139 ) , ( @xmath140 ) , and ( @xmath141 ) . here , we have taken @xmath142 , @xmath143 , @xmath144 , @xmath145 , @xmath72 and @xmath146 . the color bars in figs .
( a ) , ( b ) and ( c ) represent the values of the variables @xmath40.,width=362 ]    the coefficients @xmath147 , @xmath148 , @xmath149 and @xmath150 in the above equation are given by    @xmath151 , \nonumber \\
a_{02}&=&\frac{1}{4}\left[(4(1-\epsilon)+\epsilon \omega_1 \omega_2 ) \tilde{a}-(\epsilon-4)(\omega_1 ^ 2 a_2+\omega_2 ^ 2 a_1 ) \right ] , \nonumber \\ a_{03}&=&(1+\omega_1 ^
2)(1+\omega_2 ^ 2)-\frac{\epsilon}{2 } ( 2+\omega_1 ^ 2+\omega_2 ^ 2 ) .
\label{coeff}\end{aligned}\ ] ]    here @xmath152 .
since the characteristic equation for the eigenvalues is quartic in nature , we use the well known routh - hurwitz criteria @xcite to obtain the stable ad regions of the system . by doing so
, we find that the ad regions are bounded by the curves @xmath153 the details of obtaining the boundary curves from the r - h criteria are presented in the appendix [ app2_an ] , and the ad region bounded by the curves ( [ condi ] ) has been shown in fig .
[ anal1](b ) for @xmath125 and @xmath126 .
[ anal1](c ) which portrays the boundaries of the ad regions in the ( @xmath73 ) space clearly shows that the analytical results match nicely with that of the numerical results given in fig .
[ bif1](c ) .
next by varying the nature of oscillators contributing towards feedback we have plotted fig .
[ anal1](d ) , where we considered three cases ( i ) @xmath125 , @xmath126 ( both the oscillators contributing ) ( ii ) @xmath125 , @xmath154 ( only the high frequency oscillator contributing ) and ( iii ) @xmath155 , @xmath126 ( only the low frequency oscillator contributing ) .
the analytically obtained hopf bifurcation curves for all the above three cases have been presented in fig .
[ anal1](d ) . from the figure
, we can note that for the first two cases ( i ) and ( ii ) the os state gets revoked from the ad state for even small values of @xmath23 and the hopf bifurcation curves of these two cases are closer to each other .
but , in the case where the low frequency oscillator alone is contributing , we find that comparatively higher values of @xmath23 are needed to revoke oscillations .
this shows that for a quicker resurrection of oscillations , the feedback from the high frequency oscillator is preferable .
however , the resurrection of oscillations is possible even if one of the oscillators is contributing towards the feedback .
next , we extend our studies on the revocation of oscillations from the ad state for the case of @xmath33 oscillators , where the parametric mismatch in the system is introduced in such a way that the system has two groups of oscillators .
the first group contains @xmath156 oscillators with @xmath157 , @xmath158 and the second group has @xmath159 oscillators with @xmath160 , @xmath161 with @xmath162 and @xmath163 .
also , we consider that among the @xmath156 oscillators in the first group only the output of a sub - group of @xmath164 oscillators is fed back , in other words , the @xmath17 s of @xmath14 in ( [ sys1 ] ) take the values as @xmath165 similarly in the second group of @xmath159 oscillators only the output of @xmath166 oscillators is fedback or @xmath167 the total number of oscillators contributing towards feedback is @xmath168 .
the results corresponding to @xmath169 oscillators presented earlier in sec .
[ parame ] clearly demonstrate the appearance of ad state in the system . the temporal behavior of the system in the original os state , ad state and revoked os state are shown in fig .
[ sync ] , which shows the coherent nature among the oscillators in the first and second group .
next to study the case of the @xmath33 oscillators analytically , we first try to reduce the problem to a simpler level .    )
space for the values of ( @xmath170)@xmath36 @xmath171 , @xmath172 and @xmath173 .
( b ) shows the curves of @xmath174 for various combinations of ( @xmath164 , @xmath166 ) with fixed @xmath175 and @xmath176 .
( c ) fixing @xmath177 or @xmath178 , the values of @xmath174 for different values of @xmath166 or @xmath164 is plotted with @xmath176 .
( d ) shows the value of @xmath174 for different values of @xmath164 and @xmath166 for @xmath179 , where the values of @xmath174 are represented by the color function . in the above figures , we make the choice @xmath169 , @xmath146 and @xmath72,width=359 ]    due to the existence of coherence among the oscillators in the first and second group
, we represent the state of the oscillators in the first group by @xmath180 and the state of the oscillators in the second group by @xmath181    @xmath182    where @xmath183^{t}$ ] and @xmath184 , @xmath57 .
the coefficients in the feedback take the form @xmath185 and @xmath186 , where @xmath187 , which is the ratio of oscillators present in the first group , and @xmath188 is the ratio of oscillators present in the second group . @xmath189 and @xmath190 are the ratio of the oscillators contributing towards feedback from the first and second groups , respectively . thus @xmath191 and @xmath192 , where @xmath164 and @xmath166 are respectively the number of oscillators in the first and second groups that are contributing towards feedback .
the stability of the system ( [ sys_app ] ) around the equilibrium point @xmath97 can be studied similar to the previous case and one can obtain the characteristic equation for the eigenvalues of the above system ( [ sys_app ] ) as @xmath193 the coefficients in the above equation are given by @xmath194 where @xmath195 .
as in the previous case , using the r - h criteria we determine the ad regions of the system .
consequently , we find that the ad regions are bounded by the curves defined by @xmath196 next , using the above relations we find the ad regions in the different cases of the system .
we consider here only two cases , ( i ) @xmath197 , ( ii ) @xmath198 .    _ case-1 : @xmath169 .
_  in this case , the population of the high frequency oscillators ( @xmath156 ) and the population of the low frequency oscillators ( @xmath159 ) are equal ( note that @xmath162 ) .
first , we check the consistency of the obtained analytical results with the numerical results . for @xmath169 ,
we have plotted the boundaries of the ad regions obtained from analytical and numerical studies for different values of @xmath164 and @xmath166 in fig .
[ an_ne](a ) .
( note that for @xmath199 and @xmath200 , the numerical results have been given in fig .
[ bif1](e ) of sec .
[ parame ] ) .
the latter shows the consistency between the numerical and analytical results .
now , we look at the preferential feedback configuration for quicker resurrection of oscillations . for the purpose
, we fix the total number of oscillators contributing towards feedback as @xmath201 and vary the number of oscillators contributing from the first group ( @xmath164 ) and from the second group ( @xmath166 ) . fig .
[ an_ne](b ) shows the boundaries of the ad regions for different values of @xmath164 ( also @xmath202 ) , where we can observe that on increasing the contributions from the high frequency oscillators ( or @xmath164 ) , the resurrection of oscillations occur for lower values of @xmath23 . again in fig .
[ an_ne](c ) , we first fixed @xmath177 and plotted the critical value of @xmath23 needed for resurrection of oscillations ( @xmath174 ) for different values of @xmath166 .
similarly , we fixed @xmath178 and plotted @xmath174 for different values of @xmath164 in the same fig .
[ an_ne](c ) . from the figure
, we find that a considerable decrease in the value of @xmath174 occurs only when @xmath164 is varied .
the above results show that a feedback from the high frequency oscillators is more preferable for a quicker resurrection of oscillations .
they are also evident from the fig .
[ an_ne](d ) , where @xmath174 is represented by a color function . in the figure , one can find that the rate of change in the value of @xmath174 is larger along the @xmath164 direction .
_ case-2 : @xmath198 .
_  just as the technique shows preference over the feedback of high frequency oscillators , it shows dependence over the population in the two groups of oscillators . to illustrate the above
, we considered two situations ( i ) @xmath156 @xmath203 @xmath159 , ( ii ) @xmath156 @xmath204 @xmath159 .
[ an_ne2](a ) has been plotted for the case ( i ) where @xmath205 and @xmath206 .
similarly fig .
[ an_ne2](b ) corresponds to the case ( ii ) where @xmath207 and @xmath208 .
the values of @xmath174 for different values of @xmath164 and @xmath166 are presented in fig .
[ an_ne2 ] . comparing fig .
[ an_ne2](a ) with fig .
[ an_ne2](b ) , we find that the values of @xmath174 are larger in the case of @xmath209 than that in the case @xmath210 .
thus in the case where lower frequency oscillators are highly populated than the high frequency oscillators , we require stronger feedback to revoke oscillations .     for different values of @xmath211 for the case ( a ) @xmath210 ( @xmath205 , @xmath206 ) , ( b ) @xmath212 ( @xmath207 , @xmath208 ) , with @xmath146 and @xmath72 .
, width=366 ]
following the studies on directly coupled systems , we turn to check the validity of the proposed scheme to indirectly coupled oscillatory systems . for this purpose , we consider a collective system coupled to a dynamic environment , where the combined set of dynamical equations is characterized by    ) space corresponding to the system ( [ sys2 ] ) with @xmath27 and @xmath213 , ( b ) the reduction in the ad and od regions in the system with @xmath33 coupled oscillators in the ( @xmath76 ) space for @xmath213 .
( c ) , ( d ) , ( e ) and ( f ) the corresponding temporal behavior of the system ( [ sys2 ] ) of @xmath33 oscillators in the ad , od and the revoked oscillatory state.,width=362 ]    @xmath214    here , the state of the system along with the environment is defined by the state vector @xmath215^t$ ] , where the variables @xmath216 and @xmath40 correspond to the system and @xmath217 represent the environment . while @xmath52 , the increase in @xmath29 causes a stabilization of the trivial equilibrium point ( @xmath218)=(@xmath219 ) and a further increase in @xmath29 stabilizes a pair of nontrivial equilibrium points defined by ( @xmath218)=(@xmath220 ) , where @xmath221= @xmath222 , @xmath223 , @xmath224 and @xmath225 via pitchfork bifurcation @xcite .
the introduction of the feedback , @xmath61 , destabilizes both the ad and od states via hopf bifurcation .
such a reduction in the territories of ad and od states is illustrated in the ( @xmath34 ) space in fig.[bif2](a ) .
the curve made up of the hopf bifurcation points separates out the ad and od regions with the os region .    by extending @xmath68 to @xmath226 , fig .
[ bif2](b ) depicts the ad and od regions of the system . when @xmath52 , the ad state which arises by increasing @xmath29 is found to disappear with an increase of @xmath29 .
but for larger @xmath29 , we find the appearance of the od state . now by switching @xmath23 on , both the ad and od states are shown to be wiped out simultaneously .
the temporal behaviors of the system in the ad state and the resurrected oscillatory state which arise through the enhancement of @xmath23 are shown in figs .
[ bif2](c ) and [ bif2](d ) .
similarly , figs .
[ bif2](e ) and [ bif2](f ) show the behavior of the system at the od state and the resurrection of oscillations by an increase in @xmath23 .
) for @xmath227 and @xmath228 .
( b ) the occurrence of od through the coupling of the system for @xmath229 and @xmath52 , and ( c ) the resurrection of oscillations by introducing feedback @xmath230.,width=362 ]    in this section , we consider an interesting coupled chemical oscillator modeled by the brusselator model @xcite .
we consider the case where two identical cells are coupled .
in such a case , the functions in eq .
( [ gen ] ) describing the system are given by    @xmath231    in the absence of the coupling , the system shows stable limit cycle oscillations which has been illustrated in fig . [ br](a ) . by introducing the coupling
, the system tends to an inhomogeneous steady state @xcite .
for example , for @xmath227 , @xmath228 and @xmath229 , the system tends to an inhomogeneous steady state which has been illustrated in fig .
[ br](b ) .
in such a realistic example , by introducing the feedback , we found that the oscillations are revoked by increasing @xmath23 and it has been illustrated in fig .
[ br](c ) .
from a knowledge of the role of feedback in controlling the dynamics and the coherent activities of the system such as synchronization @xcite , we have here analyzed whether it can control oscillation quenching tendencies .
for this purpose , we have demonstrated the effect of feedback over quenching induced by parametric mismatch and symmetry breaking ( the key candidates for inducing ad and od ) , indirect coupling and for some more cases ( given in appendix [ more ] ) , through numerical as well through analytical studies wherever possible .
further , through analytical studies on ad state of a more general system ( appendix [ gen_an ] ) , we found the general applicability of the mechanism where the feedback resurrects oscillations from the ad state . in the case of od , a proper / suitable form of linear feedback would be helpful to resurrect oscillations or one can also explore the role of nonlinear feedback in such nontrivial od states . from the results obtained for different cases , we find that the trivial ad state is found to be destabilized through hopf bifurcation whereas the nontrivial od states are found to be destabilized even through saddle node type bifurcation ( see appendix [ more ] ) .
in addition to the adaptability of the technique in practical situations , we have illustrated here one more important feature of the technique namely that it does not put any restriction over the number of oscillators contributing towards feedback . even with the feedback from a few number of oscillators we can break the death state of the system and thus provide an attractive methodology in practical situations . considering a two population network , the contribution from the high frequency oscillators
are found to be more preferable compared to the feedback from the low frequency oscillators .
the work of vkc forms part of a research project sponsored by insa young scientist project .
the work forms part of an irhpa project of ml , sponsored by the department of science technology ( dst ) , government of india , who is also supported by a dae raja ramanna fellowship .
sk thanks the department of science and technology ( dst ) , government of india , for providing a inspire fellowship .
in this appendix , we show the applicability of the feedback technique over the ad state of a general two coupled system . for this purpose , we assume @xmath27 and the general forms for @xmath21 , @xmath232 and @xmath22 in eq .
( [ gen ] ) , where @xmath21 can be chosen as a polynomial in @xmath233 ,    @xmath234    in the above , @xmath235 and @xmath236 are nonlinear functions in @xmath216 and @xmath40 and the constants @xmath237 , @xmath238 , @xmath239 and @xmath240 are system parameters .
the function @xmath241 can be written as @xmath242 where @xmath243 is a nonlinear function in @xmath244 that introduces nonlinear coupling in the system .
the coupling matrix @xmath22 can be taken as    @xmath245    where @xmath246 , @xmath247 and @xmath248 .
the systems are coupled through both direct and conjugate variables , where @xmath249 introduces direct coupling , @xmath250 and @xmath251 introduce conjugate coupling .
@xmath252 , @xmath189 and @xmath190 are coupling strengths .
the system has a trivial equilibrium point at @xmath97 , which may become stable due to the parametric mismatch in the system or due to the coupling in the system .
for example , in @xcite the conjugate coupling in the system induces ad even when the oscillators are identical .
the feedback can be given as @xmath253 , where one can also add nonlinear terms in the feedback , if needed :    @xmath254    first , considering the case of coupled identical oscillators @xmath255 , @xmath256 , @xmath257 and @xmath258 ( @xmath57 ) , the linearization around the trivial equilibrium point @xmath97 can be done .
the eigenvalues corresponding to the case can be obtained easily ( note that the nonlinear terms in ( [ fgen ] ) and ( [ hgen ] ) do not play any role in the linearized equation for the trivial equilibrium point ) .
the obtained eigenvalues corresponding to the case are of the form @xmath259 in the above @xmath260 and @xmath261 are eigenvalues corresponding to the trivial equilibrium point of the system when the feedback is absent ( @xmath52 ) .
they are given by    @xmath262    depending on the values of the system and coupling parameters , the real part of the eigenvalues @xmath263 and @xmath264 in ( [ 12 ] ) and ( [ lamgen ] ) are positive or negative while @xmath52 .
when all the eigenvalues in ( [ 12 ] ) and ( [ lamgen ] ) have negative real parts , the stabilization of the equilibrium point gives rise to ad in the system . from eq .
( [ 12 ] ) , we notice that the increase in @xmath23 causes the eigenvalues @xmath265 to be more positive .
thus a destabilization of the equilibrium point @xmath97 occurs or it wipes off ad . if the equilibrium point was unstable while @xmath52 , the increase in @xmath23 never stabilizes the equilibrium point .
thus the above analysis makes clear the role of @xmath23 in destabilizing the attractor at @xmath97 .
we can observe a similar effect even in the case where parametric mismatch is present ( the case where @xmath266 , @xmath267 , @xmath268 and @xmath269 ) in the system .
the eigenvalues of the system can be obtained by solving the equation @xmath270 where @xmath271 . as @xmath148 , @xmath149 , @xmath150 are not simple in their form , we do not present them here .
although the eigenvalues @xmath272 obtained from ( [ eigll ] ) are not of the simple form , the stability of the equilibrium point @xmath97 in the different parametric regions can be found through the routh - hurwitz ( r - h ) criteria . from the r - h criteria ,
an equilibrium point is said to be stable only when all the conditions given below are satisfied by the coefficients in the eigenvalue equation ( [ eigll ] ) .
the r - h criteria are given as @xmath273 if the coefficients in the characteristic eigenvalue equation ( [ eigll ] ) fail to satisfy any one of the condition given above , the equilibrium point becomes unstable . in this aspect ,
we consider one of the simpler condition in ( [ acon ] ) , namely @xmath274 .
the condition @xmath274 is broken when @xmath275 , thus this clearly shows that an increase in @xmath23 destabilizes the equilibrium point @xmath97 .
further , more clear analytical illustration on the role of @xmath23 in the parameter mismatched system is given in sec .
[ par_ann2 ] and [ par_ann100 ] with stuart - landau model as an example .
as the above type of proof for the non - trivial od state is too cumbersome , we have illustrated the role of feedback over the state with more examples in the body of the paper as well in appendix [ more ] both numerically and analytically ( in some cases ) . from the above illustrations
, one can also notice that in the case of ad , the nonlinear feedback terms can not play any role ( as they lose their significance in the linearized limit ) and do not provide any control over it , whereas in the case of od , the nonlinear feedback also can provide a control over it .
we consider the case of two stuart - landau oscillators coupled diffusively with a repulsive link ( @xmath27 ) , as studied in @xcite .
the functions characterizing this equation have the forms ,    ) with @xmath23 for @xmath213 .
( b ) and ( c ) show the corresponding transition routes from ad and od states to os state .
( d ) the emergence of os regions in the system ( [ sys4 ] ) with the increase of @xmath23 for @xmath276 , @xmath213 and @xmath277 in ( [ sys4 ] ) . , width=355 ]    @xmath278    the eigenvalues corresponding to the trivial equilibrium point @xmath97 of the system are    @xmath279    we demonstrate the destabilization of this ad state as well as the od state corresponding to the system ( [ sys3 ] ) with respect to the feedback in fig .
[ re](a ) , where we can find that an increase in @xmath23 causes the reduction in ad and od regions of the system .
the transition route followed by the system , as it transits to oscillatory state , is shown in figs .
[ re](b ) and [ re](c ) .
these figures show that the oscillations are resurrecting from the ad state via hopf bifurcation as in the previous cases , whereas the resurrection of oscillations from od state occurs through saddle node bifurcation .    ) in the ( a ) ad state for ( @xmath34 ) @xmath36 ( @xmath280 ) , ( b ) the revoked os state for ( @xmath34 ) @xmath36 ( @xmath281 ) , ( c ) od state for ( @xmath34 ) @xmath36 ( @xmath282 ) , and ( d ) the revoked os state for ( @xmath34 ) @xmath36 ( @xmath283 ) .
, width=359 ]    recently , ad and od in the same system with large number of oscillators ( @xmath68 ) coupled globally has been seen in @xcite , whose equation is defined by    @xmath284    even with @xmath276 oscillators in system ( [ sys4 ] ) , we have shown the reduction in ad and od regions of the system in fig .
[ re](d ) .
the temporal behavior of the system for different values of @xmath29 and @xmath23 are shown in figs .
[ re2](a)-(d ) .
next , we consider the case of @xmath27 oscillators coupled through a conjugate coupling described by    @xmath285    this system has a trivial equilibrium point at @xmath286 : @xmath97 and has pairs of non - trivial equilibrium points for @xmath287 , @xmath288 : @xmath289 ) , where @xmath290 , @xmath291 and @xmath292 : @xmath293 ) , where @xmath294 , @xmath295 in which @xmath296 .
the trivial equilibrium point @xmath286 has the eigenvalues    @xmath297    for @xmath52 , the equilibrium point @xmath286 is unstable for all values of @xmath29 and for @xmath298 a saddle node type bifurcation occurs which stabilizes @xmath299 and @xmath300 as shown in fig .
[ con](a ) . for @xmath61 ,
the eigenvalues corresponding to the equilibrium point @xmath286 are still unstable , and thus the system is still free of ad .
then , the equilibrium points @xmath299 and @xmath300 are also destabilized through hopf bifurcations which is demonstrated in fig .
[ con](b ) .
then , the od regions of the system in the ( @xmath34 ) space is shown in fig .
[ con](c ) , which clearly demonstrates the destabilization of od with the introduction of feedback .     in the case of conjugate coupling .
( b ) the destabilization of the od state with the introduction of feedback via hopf bifurcation ( @xmath301 ) .
( c ) the reduction in the od regions of the system in the ( @xmath34 ) space for @xmath302.,width=359 ]      next , we illustrate the role of feedback in the case two van der pol oscillators coupled diffusively through a repulsive link .
the corresponding dynamical equations are defined through @xcite    @xmath303    the ad and od regions of the system in the ( @xmath34 ) space are given in fig .
[ vd ] which show the reduction in the ad and od regions with respect to @xmath23 .
the transition from ad to oscillatory state occurs through a hopf bifurcation . on the other hand considering the transition from od to oscillatory state , the od state
is transformed to ad through inverse pitchfork bifurcation and the oscillatory state arises from the ad state through hopf bifurcation .
this shows that the role of the feedback in setting oscillations back in the system is not restricted to any particular oscillator . in the next example , we show that the feedback can destabilize oscillation quenching scenario even in chaotic oscillators .    ) with respect to @xmath29 and @xmath23 for @xmath304 .
( b ) the transition route from od state to os state via ad state.,width=355 ]      we consider the case of the @xmath27 coupled chaotic rssler system @xcite defined by    @xmath305    here , the oscillators are coupled directly by diffusive type coupling and are also coupled indirectly to an environment defined by the variable @xmath306 .
this system ( [ ros ] ) has only non - zero equilibrium points which include @xmath288 : @xmath307 @xmath36 @xmath308 , where , @xmath309 .
for the parametric choice @xmath310 , @xmath311 , @xmath312 and @xmath313 , the equilibrium point @xmath300 is stable .
the od regions corresponding to the above equilibrium points are given in figs .
[ rtt](a ) and [ rtt](b ) with respect to the direct coupling strength ( @xmath29 ) and with respect to the indirect coupling strength ( @xmath272 ) , respectively .
the above figures clearly show that the feedback is applicable even for the case of chaotic oscillators .    )
( a ) in the ( @xmath314 ) space for @xmath315 , ( b ) in the ( @xmath76 ) space for @xmath316 . in both the figures ,
we fixed @xmath310 , @xmath311 , @xmath317 and @xmath318.,width=321 ]      to illustrate further the role of feedback in chaotic oscillators , we consider the @xmath27 coupled sprott and lorentz oscillators which are defined respectively by + ( i)two coupled sprott systems with repulsive link @xcite ,    @xmath319     in the case of ( a ) the sprott system ( [ spr ] ) for @xmath320 and ( b ) the lorenz system ( [ lor ] ) for @xmath321 , @xmath322 and @xmath323.,width=355 ]    and ( ii)diffusive coupling among two lorenz oscillators @xcite ,    @xmath324    the ad and od regions corresponding to the coupled sprott and lorenz systems ( ( [ spr ] ) and ( [ lor ] ) ) are given in fig .
[ s ] , which again confirm that the feedback wipes out ad and od in the system .
the coefficients in the characteristic eigenvalue equation ( [ eveq ] ) are given by    @xmath151 , \nonumber \\
a_{02}&=&\frac{1}{4}\left[(4(1-\epsilon)+\epsilon \omega_1 \omega_2 ) \tilde{a}-(\epsilon-4)(\omega_1 ^ 2 a_2+\omega_2 ^ 2 a_1 ) \right ] , \nonumber \\ a_{03}&=&(1+\omega_1 ^
2)(1+\omega_2 ^ 2)-\frac{\epsilon}{2 } ( 2+\omega_1 ^ 2+\omega_2 ^ 2 ) .
\label{coeff_ap}\end{aligned}\ ] ]    here @xmath152 . the eigenvalue equation in ( [ eveq ] ) can be solved directly to get the stable regions of the steady state @xmath97 . on the other hand
, we can use the r - h criteria to have a closer look at the stable regions of the system ( [ sys1 ] ) . according the r - h criteria
, the stable region corresponds to the region in which @xmath325    ) are satisfied is depicted in the above figure .
the light - blue shaded region ( i ) which lies under the light - blue curve ( curve -@xmath19 ) corresponds to @xmath274 , the region under the pink curve ( curve -@xmath66 ) which is denoted by region ( ii ) corresponds to @xmath326 , green shaded region ( iii ) corresponds to @xmath327 .
the regions satisfying the conditions ( ii ) and ( iii ) in ( [ rhcri ] ) are the gray shaded region ( iv ) ( the region under curve -@xmath328 ) and red shaded region ( v ) ( the region under curve -@xmath329 ) , separately .
the other condition @xmath330 is satisfied everywhere in the considered region of ( @xmath76 ) .
the region in which all the conditions in ( [ rhcri ] ) are satisfied is denoted by blue shaded region ( vi ) ( the region under curve -@xmath331).,width=359 ]    now we consider the above criteria one by one and obtain the required stable region of the steady state .
+ ( i)a .
@xmath274 :  this condition will be satisfied when @xmath332 this simple condition promises that above the value of @xmath333 , the trivial equilibrium point can never be stable and thus oscillation can be revoked by an increase in the value of @xmath23 .
the region in the @xmath334 space in which the above condition is satisfied is denoted by i in fig .
[ ap_fig ] .
+ ( i)b .
@xmath330 :   when @xmath335 , if the condition ( @xmath330 ) is satisfied for @xmath52 then it will be satisfied for all values of @xmath23 . if the condition @xmath330 is not satisfied while @xmath52 then the condition will not be satisfied for any value of @xmath23 .
when @xmath336 , the condition @xmath330 will be satisfied only for the values of @xmath23 given by @xmath337 , \nonumber \\ \eta&<&\frac{2}{\tilde{a}}\left[\epsilon-3 - \sqrt{\epsilon^2 - 3 \epsilon+3-\omega_1 ^ 2-\omega_2 ^ 2}\right ] .
\label{ic}\end{aligned}\ ] ] this condition @xmath330 is satisfied in the whole region considered in fig .
[ ap_fig ] .
+ ( i)c .
@xmath326 :  when @xmath338 , if the condition is satisfied for @xmath52 then it will continue to be satisfied for all values of @xmath23 , if the condition is not satisfied for @xmath52 , by the variation of @xmath23 also the condition remains to be unsatisfied .
when @xmath339 , it will be satisfied only for the values of @xmath23 given by @xmath340 the region of @xmath334 in which the above condition is satisfied is denoted by region ( ii ) which lies under pink curve ( or curve -@xmath66 ) in fig .
[ ap_fig ] .
+ ( i)d .
@xmath327 :  similar to the previous conditions , if @xmath341 , the condition will be satisfied when @xmath342 otherwise it will be satisfied for all values of @xmath23 only if the condition is satisfied for @xmath52 .
the region in which the above condition is satisfied is shown in fig .
[ ap_fig ] as region iii .
+ ( ii )  @xmath343 :  we can find that @xmath344 is a cubic polynomial in @xmath23 , and the region in which the above condition is satisfied is given by gray shaded region ( iv ) .
the real roots of @xmath23 satisfying the equation @xmath345 forms the boundary of the region .
+ ( iii )  @xmath346 : the region satisfying this condition is shown by red shaded region ( v ) , whose boundary is the solution of the quintic equation @xmath347 .    from fig .
[ ap_fig ] , we can find that the region satisfying all the criteria given in ( [ rhcri ] ) is the region bounded between the curves @xmath348 the region is denoted by blue shaded region vi in fig .
[ ap_fig ] . in this region ,
the equilibrium point @xmath97 is found to be stable or ad occurs in the region .
99 g. saxena , a. prasad , and r. ramaswamy , phys . rep . * 521 * 205 ( 2012 ) .
a. koseska , e. volkov , and j. kurths , phys . rep . * 531 * 173 ( 2013 ) .
a. koseska , e. volkov , and j. kurths , phys .
lett . * 111 * 024103 ( 2013 ) .
t. banerjee and d. ghosh , phys .
e * 89 * 062902 ( 2014 ) .
k. konishi , phys .
e * 68 * 067202 ( 2003 ) .
w. zou , d. v. senthilkumar , j. duan , and j. kurths , phys .
e * 90 * 032906 ( 2014 ) .
f. m. atay , phys .
91 , 094101 ( 2003 ) .
w. zou , d. v. senthilkumar , y. tang , y. wu , j. lu , and j. kurths , phys .
e * 88 * 032916 ( 2013 ) .
m. lakshmanan and d.v .
senthilkumar , _ dynamics of nonlinear time - delay systems _ , ( springer , berlin , 2010 ) .
a. prasad , m. dhamala , b.m .
adhikari , and r. ramaswamy , phys .
e * 81 * 027201 ( 2010 ) .
w. zou , d. v. senthilkumar , a. koseska , and j. kurths , phys .
e * 88 * 050901(r ) ( 2013 ) . c. r. hens , p. pal , s. k. bhowmick , p. k. roy , a. sen and s. k. dana , phys .
e * 89 * 032901 ( 2014 ) .
g. b. ermentrout and n. kopell , siam j. appl .
math * 50 * 125 ( 1990 ) .
r. curtu , _ physica d _ * 239 * 504 ( 2010 ) .
m. dolnik and m. marek , j. phys .
chem * 92 * 2452 ( 1988 ) .
m. toiya , v. k. vanag , and i. r. epstein , angew .
, int . ed . * 47 * , 7753 ( 2008 ) .
m. heinrich , t. dahms , v. flunkert , s. w. teitsworth , e. schll , new j. phys * 12 * 113030 ( 2010 )
. d. ruwisch , m. bode , d. volkov , and e. volkov int .
j. bifurcation chaos appl .
* 9 * 1969 ( 1999 ) .
b. gallego and p. cessi , j. clim .
* 14 * 2815 ( 2001 ) .
s. ching , p. l. purdon , s. vijayan , n. j. kopell , and e. n. brown , pnas , * 109 * 3095 ( 2012 ) .
j. jalife , r. a. gray , g. e. morley , and j. m. davidenko , chaos * 8 * 79 ( 1998 ) .
l. rubchinsky and m. sushchik , phys .
e * 62 * 6440 ( 2000 ) .
w. zou , d. v. senthilkumar , m. zhan , and j. kurths , phys .
lett * 111 * 014101 ( 2013 ) .
s. kak , circuits , systems , and signal processing * 12 * 263 ( 1993 ) .
a. becskei , b. sraphin , and l. serrano , embo .
j. * 20 * 2528 ( 2001 ) .
s. draghici , int . j. neural syst .
* 8 * 113 ( 1997 ) .
v. k. chandrasekar , j. h. sheeba , and m. lakshmanan , chaos * 20 * 045106 ( 2010 ) m. g. rosenblum and a. s. pikovsky , phys .
* 92 * 114102 ( 2004 ) .
o. v. popovych , c. hauptmann , and p. a. tass , phys .
lett . * 94 * 164102 ( 2005 ) .
g. yuan , x. zhang , and z. wang , optik - int .
j. light electron opt . * 125 * 1950 ( 2014 ) .
g. f. franklin , j. d. powell , and a. emami - naeini , _ feedback control of dynamical systems _ ( addison - wesley , 3rd ed . , uk , 1994 ) . j. e. gough , r. gohm , and m. yanagisawa , phys . rev . a * 78 * 062104 ( 2008 ) .
s. lloyd , phys .
a * 62 * 022108 ( 2000 ) .
s. little and p. brown , ann . n. y. acad .
sci * 1265 * 9 ( 2012 ) .
s. santaniello , g. fiengo , and l. glielmo , _ in proceedings of the ieee international conference on control application ( cca 2008 ) _
( ieee , new york , 2008 ) pp .
666 - 671 .
j. a. anderson , _ an introduction to neural networks _ ,
( a bradford book , 3rd ed . , ( 1997 ) ) . y. kuromato , _ chemical oscillations , waves and turbulence _ , ( springer , berlin , 1984 ) . c. u. choe , t. dahms , p. hvel , and e. schll , phys
e * 81 * 025205(r ) ( 2010 ) .
t. aoyagi , phys .
. lett . * 74 * 4075 ( 1995 ) .
s. uchiyama , physica a * 391 * 2807 ( 2012 ) . n. tukhlina and m. rosenblum , j. biol .
* 34 * 301 ( 2008 ) .
o. dhuys , i. fischer , j. danckaert , and r. vicente , phys .
e * 83 * 046223 ( 2011 ) .
k. ldge , _ nonlinear laser dynamics _ , ( wiley - vch , germany , 2011 ) p. r. bandyopadhyay and a. m. hellum , sci .
* 4 * 6650 ( 2014 ) .
e. m. izhikevich , _ dynamical systems in neuroscience : the geometry of excitability and bursting _ , ( mit press , cambridge , 2007 ) . n. tukhlina , m. rosenblum , a. pikovsky , and j. kurths , phys . rev .
e * 75 * 011918 ( 2007 ) .
o. v. popovych , c. hauptmann , and p. a. tass , biol
. cybern . * 95 * 69 ( 2006 ) .
k. pyragas , o. v. popovych , and p. a. tass , euro .
. lett . * 80 * 40002 ( 2007 ) .
j. h. sheeba , v. k. chandrasekar , and m lakshmanan , phys .
lett * 103 * 074101 ( 2009 ) .
g. c. sethia and a. sen , phys .
lett * 112 * 144101 ( 2014 ) .
l. schmidt and k. krischer , phys .
lett * 114 * 034101 ( 2015 ) .
a. zakharova , m. kapeller , and e. schll , phys .
lett . * 112 * 154101 ( 2014 ) .
a. koseska , e. volkov , and j. kurths , europhys . lett . *
85 * 28002 ( 2009 ) .
a. koseska , e. volkov , and j. kurths , chaos * 20 * 023132 ( 2010 ) .
m. lakshmanan , s. rajasekar , _ nonlinear dynamics : integrability , chaos and patterns _ , ( springer - verlag , berlin , 2003 ) k. bar - eli , physica d * 14 * 242 ( 1985 ) . c. r. hens , o. i. olusola , p. pal , and s. k. dana , phys .
e * 88 * 034902 ( 2013 ) .
m. nandan , c. r. hens , p. pal , and s. k. dana , chaos * 24 * 043103 ( 2014 ) .
v. resmi , g. ambika , and r. e. amritkar , phys .
e * 84 * 046212 ( 2011 ) .