Gnral mchansm for ampltud dath n coupld systms V. Rsm and G. Ambka Indan Insttut of Scnc Educaton and Rsarch, Pun - 4111, Inda R. E. Amrtkar Physcal Rsarch Laboratory, Ahmdabad - 389, Inda arxv:111.534v [nln.cd] 6 Oct 11 W ntroduc a gnral mchansm for ampltud dath n coupld synchronzabl dynamcal systms. It s known that whn two systms ar coupld drctly, thy can synchronz undr sutabl condtons. Whn an ndrct fdback couplng through an nvronmnt or an xtrnal systm s ntroducd n thm, t s found to nduc a tndncy for ant-synchronzaton. W show that, for suffcnt strngths, ths two comptng ffcts can lad to ampltud dath. W provd a gnral stablty analyss that gvs th thrshold valus for onst of ampltud dath. W study n dtal th natur of th transton to dath n svral spcfc cass and fnd that th transtons can b of two typs - contnuous and dscontnuous. By choosng a varty of dynamcs for xampl, prodc, chaotc, hypr chaotc, and tm-dlay systms, w llustrat that ths mchansm s qut gnral and works for dffrnt typs of drct couplng, such as dffusv, rplacmnt, and synaptc couplngs and for dffrnt dampd dynamcs of th nvronmnt. PACS numbrs: 5.45.Gg; 5,45.Pq; 5.45.Xt 1. INTRODUCTION Th dynamcs of many complx systms can b undrstood as th collctv bhavor of a larg numbr of dynamcal unts coupld va thr mutual ntractons. Th dynamcs of such connctd systms has bn an ntrstng topc ofstudy spcally du to ts rlvanc n undrstandng a larg varty of natural systms. Basd on th th natur of ntractons among th coupld unts, thy can xhbt many mrgnt phnomna such as synchronzaton, hystrss, phas lockng, ampltud dath, and oscllator dath [1 3]. Among ths, th phnomnon of synchronzaton s th most wdly studd and has rlvanc n many contxts, such as nuronal ntworks, communcaton, lasr systms, tc. [1]. So also, th qunchng or supprsson of dynamcs calld ampltud dath s anothr mrgnt phnomnon of qual rlvanc n such systms. Ths can lad to ntrstng slf-adjustabl control mchansms and plays a promnnt rol as an ffcnt rgulator of th dynamcs. Th occurrnc of ampltud dath has bn rportd n many cass, such as chmcal ractons [4 7], bologcal oscllators [8 1], coupld lasr systms [11, 1] and rlatvstc magntrons [13]. W would lk to projct th mportanc of th phnomnon of ampltud dath n coupld systms n two contxts: on, as a dsrabl control mchansm n cass such as coupld lasrs whr t lads to stablzaton [14, 15] and two, as a pathologcal cas of oscllaton supprsson or dsrupton n cass lk nuronal dsordrs such as Alzhmr s dsas, Parknson s dsas, tc. [16 18]. Th mchansms so far rportd to nduc ampltud Elctronc addrss: v.rsm@srpun.ac.n Elctronc addrss: g.ambka@srpun.ac.n Elctronc addrss: amrtkar@prl.rs.n dath n coupld systms ar d-tunng of oscllators undr strong couplng [19 1], couplng through conjugat varabls [, 3], dynamc couplng [4], and dlay n couplng du to fnt propagaton or nformaton procssng spds [5 8]. Dstrbutd dlays rathr than dscrt or constant dlays hav bn proposd as mor ralstc modls n cology and nurobology, whr th varanc of th dlay plays a rlvant rol [9]. So also, ampltud dath has bn studd n th contxt of attractv and rpulsv couplngs n two chaotc Lornz systms [3]. In all ths mchansms, dath occurs dynamcallydu to th targtngofth unts to on ormor ofth qulbrum statsordu to th stablzaton ofon of ths stats. Th qulbrum stats or fxd ponts can b thr that of th uncoupld systm or thos volvd by couplng. Whl ths mchansms can modl th ampltud dath obsrvd n coupld systms of oscllators, w fnd that all ths mthods ar systm spcfc and may not work n a gnral cas. In th cas of dath by dlay couplng, th lmtatons of th mthod hav bn rportd n svral cass such as prodc [4] and chaotc systms [31, 3]. Morovr, thr ar many cass such as nuronal dsordrs whr dprsson of actvty or daths dutoth prsncofanothragncyormdum [16, 17]. For such cass, th mchansm of dath s stll not fully undrstood and non of th abov mntond mchansms so far rportd s applcabl. In ths work, w ntroduc a mchansm for ampltud dath causd by an ndrct fdback couplng through a dynamc nvronmnt, n addton to drct couplng.w ssntally projct th rol of th nvronmnt n controllng th dynamcs of connctd systms. W fnd that, whl t ssntally xplans qunchng of actvty or supprsson nducd by an xtrnal mdum or agnt, ths mthod can also srv as a gnral mchansm to nduc dath n coupld synchronzabl systms. Its gnralty ls n th fact that t sms to work n any coupld sys-
tm that can synchronz. It s ffctv n qunchng dynamcs n a varty of systms such as prodc oscllators, chaotc systms, hyprchaotc systms, and dlay systms. W show that ths mthod also nducs ampltud dath n systms wth dffrnt forms of drct couplng ntractons, lk dffusv, rplacmnt, synaptc couplng, tc. Spcfcally, w dmonstrat that th prsnt mthod, wth a vard modl for th nvronmnt, works n th cas of hyprchaotc systms for whch dlay couplng s not ffctv to nduc ampltud dath. As such, t s an mportant stp n mthodology toward achvng controls or stablzaton to dsrabl prformanc n many practcal cass. Th rlvanc of ths mthod ls n th fact that dath can b ngnrd and can b asly mplmntd n any systm wth coupld synchronzabl unts. In th prsnt work w us th ndrct fdback couplng through th nvronmnt of our arlr work whr w showd that such a couplng can nduc ant-phas- (or ant) synchronzaton n two systms whch ar not drctly connctd[33]. Consdr two systms coupld drctly such that wth adquat strngth of couplng thy can xhbt synchronous bhavor. Thn f w ntroduc an addtonal ndrct fdback couplng through th nvronmnt or anothr xtrnal systm such that t nducs a tndncy for ant-synchronzaton, thn for suffcnt strngths, ths two comptng tndncs can lad to ampltud dath. W fnd that, n th stat of ampltud dath, th subsystms stablz to a fxd pont of th coupld systm. W also show that th mthod ntroducd hr can nduc ampltud dath n coupld systms wth dffrnt typs of dynamcs for th nvronmnt and for dffrnt typs of drct couplng. W dvlop an approxmat stablty analyss whch provds th thrshold or crtcal valus of th couplng strngth for ampltud dath n th gnral contxt. Drct numrcal smulaton gvng th rgons of ampltud dath n th spac of couplng strngths agrs wll wth th transton curvs obtand from th stablty analyss. W also analyz n dtal th natur of th transton to th ampltud dath stat.w fnd that all th spcfc cass studd xhbt thr contnuous or dscontnuous transtons to dath. In th contnuous cas, as llustratd by two coupld Rösslr systms, durng th transton th full rvrs prod-doublng scnaro s obsrvd, and th systm rachs a on-cycl stat bfor ampltud dath occurs. Th transton to dath thn occurs du to a supr crtcal Hopf bfurcaton. In th dscontnuous cas, th transton s suddn du to th dsapparanc of a dstant attractor and stablzaton of a fxd pont. For two coupld Lornz systms, w fnd that th transton to dath s probably va a sub-crtcal Hopf bfurcaton wth long transnts, and pror to ths, th systms go through a stat of frustraton btwn synchronzd and ant-synchronzd bhavor.. AMPLITUDE DEATH VIA DIRECT AND INDIRECT COUPLING W start wth two systms coupld mutually wth two typs of couplng, namly a drct dffusv couplng and an ndrct couplng through an nvronmnt. Th dynamcs can b wrttn as ẋ 1 = f(x 1 )+ε d β(x x 1 )+ γy, ẋ = f(x )+ε d β(x 1 x )+ γy, (1a) (1b) ẏ = κy γt (x 1 +x ). (1c) Hr, x 1 and x rprsnt two m-dmnsonal oscllators whos ntrnsc dynamcs s gvn by f(x 1 ) and f(x ) rspctvly. Th systms ar mutually coupld usng dffusv couplng [ th scond trm n Eqs. (1a] and (1b)). Th nvronmnt s modld by a on-dmnsonal ovrdampd oscllator y wth a dampng paramtr κ. Th nvronmnt s kpt actv by fdback from both systms as gvn by th last trm n Eq. (1c). Both systms also gt fdback from y [ th last trm n Eqs. (1a and 1b)]. β s a matrx (m m) wth lmnts and 1 and dfns th componnts of x 1 and x that tak part n th dffusv couplng. For smplcty, w tak β to b dagonal, β = dag(β 1,β,...,β m ), and n numrcal smulatons only on componnt β 1 s assumd to b nonzro. γ s a column matrx (m 1), wth lmnts zro or 1, and t dcds th componnts of x 1 and x that gts fdback from th nvronmnt. γ T s th transpos of γ and t dcds th componnts of x 1 and x that gv fdback to th nvronmnt. W tak ε d to b th strngth of drct dffusv couplng btwn th systms, and th strngth of fdback couplng btwn th systms and th nvronmnt. Th drct couplng ε d, gvs a synchronzng tndncy btwn th two systms whl th couplng through th nvronmnt, gvs an ant-synchronzng tndncy. Thus, whn both th couplngs ar abov thr crtcal valus, thr s a comptton btwn th two tndncs, and th nt rsult s th ampltud dath. In th ampltud dath stat, th coupld systms ar drvn to a fxd pont. W llustrat th abov schm for two coupld chaotc Rösslr systms rprsntd by th followng quatons (,j = 1,, j): ẋ 1 = x x 3 +ε d (x j1 x 1 )+ y, ẋ = x 1 +ax, ẋ 3 = b+x 3 (x 1 c), ẏ = κy ε x 1. () Th rsultng tm srs for a synchronzd stat wth only drct couplng, an ant-phas-synchronzd stat wth only ndrct couplng, and th ampltud dath stat wth both drct and ndrct couplngs ar shown n Fg. 1. Whn =, and ε d s suffcntly larg, w
3 x 11, x 1 x 11, x 1 x 11, x 1 5-5 5-5 5-5 5 5 Tm FIG. 1: Tm srs of th frst varabls x 1, = 1, of two coupld Rösslr systms. Synchronzaton for (ε d,) = (.,.). Ant-phas synchronzaton (ε d,) = (.,1.).(c) Ampltud dath for (ε d,) = (.,1.). Hr, th Rösslr paramtrs ar a = b =.1,c = 18. Th dampng paramtr of th nvronmnt s takn to b κ = 1. (c) x 11, x 1 x 11, x 1 x 11, x 1 15-15 15-15 15-15 1 Tm FIG. : Tm srs of th frst varabls x 1 of th coupld Lornz systms for th paramtrs (σ = 1,r = 8,b = 8/3) and κ = 1. Synchronzaton for (ε d,) = (5,). Antsynchronzaton for (ε d,) = (,1). (c) Ampltud dath for (ε d,) = (5,1). (c) obsrv synchronzaton[fg. 1]. Whn s ncrasd for ε d =, th systms ar n an ant-phas synchronzd stat [Fg. 1]. Whn both and ε d ar suffcntly larg, th systms stablz to a stat of ampltud dath [Fg. 1(c)]. W apply th sam schm to two coupld chaotc Lornz systms as gvn by th followng quatons: ẋ 1 = σ(x x 1 )+ε d (x j1 x 1 )+ y, ẋ = (r x 3 )x 1 x, ẋ 3 = x 1 x bx 3, ẏ = κy ε x 1. (3) W fnd that ampltud dath occurs n ths cas also. Ths s llustratd n Fg., whr tm srs for a synchronzd stat [Fg. ], ant-synchronzd stat [Fg. ] and an ampltud dath stat [Fg. (c)] ar shown. So far w hav prsntd th mthod for dntcal systms. Howvr, th mthod also works for nondntcal systms. In gnral, for nondntcal systms th drct couplng wll gv a gnralzd synchronzaton btwn th coupld systms. Smlarly, th ant-synchronzaton du to th ndrct couplng wll also bcom of a gnralzd typ. Th combnaton of drct and ndrct couplng stll lads to ampltud dath. As an xampl consdr two coupld Rösslr systms. In Eq. (), w kp th paramtrs of on systm fxd and vary th paramtr c of th othr systm. W fnd that for suffcnt strngth of couplng th systms go to th ampltud dath stat vn for larg dvatons n c. Also, th ampltud dath stat occurs whn th ndvdual nonntractng systms ar n dffrnt dynamcal rgms. ThssshownnFg.3, whrtmsrsforagnralzd x 11, x 1 x 11, x 1 x 11, x 1 1-1 1-1 1-1 6 1 Tm FIG. 3: Tm srs of th frst varabls x 1, = 1, of two coupld non-dntcal Rösslr systms. Hr, th paramtrs a and b of th two Rösslr systms ar chosn to b sam (a = b =.1), whl th paramtrs c of th systms ar chosn such that on of th systms s n a chaotc stat (c = 18) and th othr n a prodc stat (c = 4). Th dampng paramtr of th nvronmnt s chosn to b κ = 1. Gnralzd synchronzatonfor (ε d,) = (1.,.). Antphas-synchronzaton for (ε d,) = (.,1.). (c) Ampltud dath for (ε d,) = (1.,1.). synchronzd stat [Fg. 3], ant-phas synchronzd stat [Fg. 3] and an ampltud dath stat [Fg. 3(c)] ar shown for two non-dntcal Rösslr systms. 3. LINEAR STABILITY ANALYSIS W prsnt an analyss of th stablty of th stady stat of two systms coupld va th schm gvn n (c)
4 Eq. (1). For ths, w wrt th varatonal quatons formd by lnarzng Eq. (1) as ξ 1 = f (x 1 )ξ 1 +ε d β(ξ ξ 1 )+ γz, ξ = f (x )ξ +ε d β(ξ 1 ξ )+ γz, ż = κz γt (ξ 1 +ξ ), (4) whr ξ 1, ξ and z ar small dvatons from th rspctv valus. W dnot synchronzng and antsynchronzng tndncs through th varabls ξ s and ξ a rspctvly as gvn by Thn Eq. (4) can b wrttn as ξ s = f (x 1 )+f (x ) ξ s = ξ 1 ξ, ξ a = ξ 1 +ξ. (5) ξ s + f (x 1 ) f (x ) ξ a ε d βξ s, ξ a = f (x 1 ) f (x ) ż = κz γt ξ a. (6) ξ s + f (x 1 )+f (x ) ξ a + γz, For stablty, all th Lyapunov xponnts obtand from Eq. (6) should b ngatv. In gnral, t s not asy to analyz th stablty of th synchronzd stat from Eq. (6). Howvr, consdrabl progrss can b mad f w assum that th tm avrag valus of f (x 1 ) and f (x ) ar approxmatly th sam and can b rplacd by an ffctv constant valu µ. In ths approxmaton w trat ξ 1 and ξ as scalars. Ths approxmaton smplfs th problm such that only th rlvant faturs rman and s xpctd to gv faturs nar th transton. Ths typ of approxmaton was usd n Rfs. [33, 34] and t was notd that t dscrbs th ovrall faturs of th phas dagram rasonably wll. Thus, Eq. (4) bcoms ξ s = µξ s ε d ξ s, ξ a = µξ a + z, (7a) (7b) ż = κz ξ a. (7c) W not that Eqs. (7b) and (7c) ar coupld whl Eq.(7a) s ndpndnt of th othr two. Th synchronzng tndncy s gvn by Eq. (7a) and th corrspondng Lyapunov xponnt s λ 1 = µ ε d. (8) Th ant-synchronzng tndncy s gvn by Eqs. (7b) and (7c). Th corrspondng Jacoban s ( ) µ ε J = / κ and th gnvalus ar λ,3 = (µ κ)± (µ κ) 4(ε µκ). (9) As notd n th prcdng scton, ampltud dath s obtand whn both synchronzng and ant-synchronzng tndncs ar prsnt and th corrspondng couplng constants ar gratr than th crtcal valus rqurd for th rspctv phnomna. Th synchronzng and ant-synchronzng tndncs bcom ffctv whn th corrspondng Lyapunov xponnts,.., th ral parts of th gnvalus, ar ngatv. From Eq. (8) w obtan th condton ε d > µ/, (1) whl from Eq. (9) w gt th followng condtons: (1) If (µ κ) < 4(ε µκ), λ,3 ar complx and th condton of stablty s κ > µ, (11) () If (µ κ) > 4(ε µκ), λ,3 ar ral and th stablty condton bcoms κ > µ and ε > µκ. (1) If Eqs. (1) and (11) or (1) ar smultanously satsfd, th oscllatons can not occur and th systms stablz to a stady stat of ampltud dath. For a gvn κ and µ, th transton to ampltud dath occurs at crtcal couplng strngths ε dc and c that ar ndpndnt of ach othr. That s, and ε dc = const, (13) c = const. (14) Ths gnral stablty crtra ar numrcally vrfd for dffrnt systms n th followng scton. W can also analyz th stablty of ampltud dath by notng that th ampltud dath corrsponds to a fxd pont of th coupld systm. Thus, th condton for th stablty of ampltud dath s that all Lyapunov xponnts of th fxd pont ar ngatv. Ths can b don for dffrnt systms numrcally and s dscussd n th nxt scton. 4. NUMERICAL ANALYSIS W apply our schm to two chaotc systms, Rösslr and Lornz. 4.1. Coupld Rösslr systms Now, w apply th schm of couplng ntroducd n Eq. (1) to th cas of two chaotc Rösslr systms. Th occurrnc of ampltud dath n ths cas s llustratd n Fg. 1(c). Ths s furthr confrmd by calculatng th Lyapunov xponnts [35] also. Whn th systms ar n
5..6 I II II ε d.1.3 I.45.55.65 1 3 κ FIG. 4: Transton from rgon of oscllatons (I) to rgon of ampltud dath (II) n th paramtr plan -ε d for coupld Rösslr systms. Numrcal smulatons ar don wth κ = 1. Th ponts mark th paramtr valus (c,ε dc ) at whch th transton to ampltud dath occurs. Sold trangls show th transton to ampltud dath as ε d s ncrasd for a constant. Th horzontal ln formd by ths trangls confrms th stablty condton Eq. (13). Smlarly, th crcls corrspond totranston tothampltuddathstatas s ncrasdfor a constant ε d and confrm th stablty condton of Eq. (14). th ampltud dath stat, all th Lyapunov xponnts of th coupld systm ar found to b ngatv. Fgur 9 shows th largst Lyapunov xponnt of th coupld systm as a functon of couplng strngth. W study th transton to dath by dntfyng rgons of ampltud dath n th paramtr plan of couplng strngths ε d for a chosn valu of κ. To charactrz th stat of ampltud dath, w us an ndx A, dfnd as th dffrnc btwn th global maxmum and global mnmum valus of th tm srs of th systm ovr a suffcntly long ntrval. Th cas whr A = rprsnts th stat of ampltud dath, whl A ndcats oscllatory dynamcs. Th paramtr valu at whch A bcoms s thus dntfd as th thrshold for onst of stablty of ampltud dath stats. Usng ths ndx, th transton curvs n th paramtr plan ε d ar plottd n Fg. 4. W not that th ponts obtand from numrcal smulatons agr wth th stablty crtra Eqs. (13) and (14) obtand n th prcdng scton. W also vrfy numrcally th crtra for transton to ampltud dath gvn n Eq. (1). For ths, th numrcally obtand valus of ε c ar plottd aganst κ n Fg. 5. Th ln corrsponds to th stablty condton Eq. (1) and th ponts ar obtand from numrcal smulatons. It s sn that th agrmnt s good for largr valus of κ. Howvr, for small valus of κ, th ponts FIG. 5: Transton from rgon of oscllatons (I) to rgon of ampltud dath (II) n th paramtr plan κ ε for th coupld Rösslr systm. Th ponts of ampltud dath ar obtand numrcally whn th ndx A bcoms zro. Th sold curv s a lnar ft corrspondng to th stablty condton Eq. (1), wth th ffctv µ =.1. Th dvaton from straght ln bhavor for small valus of κ s dscussd n th txt. dvat from straght ln bhavor. Th rason can b sn from Eq. (11) whch gvs th lowr lmt on κ. As κ dcrass, th dampng of th nvronmnt varabl y s rducd. Howvr, ths dampng s ssntal for th ant-synchronzng tndncy arsng from th couplng to th nvronmnt. Ths lads to th dvatons for small valus of κ. Fortwocoupld Rösslrsystmsasgvnn Eq.(), w study th complt phas dagram n th paramtr plan of couplng strngths, dntfyng th rgons of dffrnt dynamc stats such as ampltud dath, complt synchronzaton and ant-synchronzaton. Ampltud dath stats ar dntfd usng th ndx A as mntond abov. To dntfy synchronzd or ant-synchronzd stats w us th asymptotc corrlaton valus as th ndx, calculatd usng th quaton C = < (x 11 (t) < x 11 (t) >)(x 1 (t) < x 1 (t) >) > < (x11 (t) < x 11 (t) >) >< (x 1 (t) < x 1 (t) >) >. (15) Th phas dagram thus obtand for th coupld Rösslr systm s shown n Fg. 6. Whn th couplng strngths ε d and ar small, th systms ar not synchronzd (wht rgon). For small valus of, whn ε d s ncrasd, th systms synchronz (lght-gray rgon). Whn s ncrasd, th systms bcom antsynchronzd (dark-gray rgon). Whn both th couplng strngths ar abov a crtan thrshold as gvn by th stablty condtons Eqs. (13) and (14), th systms
6.3 4 ε d.15 A.3.6.4.43.46 FIG. 6: Rgons of dffrnt dynamcal stats n th paramtrplanofcouplngstrngthsε ε d nthcas oftwocoupld Rösslr systms. Th ndcs usd to dntfy th dffrnt rgons ar th avrag corrlaton C [ Eq. (15) ] and th ndx A. Th black rgon corrsponds to th stat of ampltud dath (A ), th lght-gray rgon to th synchronzd stat (C 1), th dark-gray rgon to th ant-synchronzd stat (C 1) and th wht rgon corrspond to th stat whr C < 1. Hr, th paramtrs ar th sam as n Fg. 1. FIG. 7: Th ndx A as a functon of for a fxd valu of ε d =. for two coupld Rösslr systms. As s ncrasd, w obsrv a contnuous transton to th stat of ampltud dath (c.45). x 3 stablz to th stat of ampltud dath (black rgon). W also not that th transton from complt synchronzaton to ant-synchronzaton corrsponds to a phas transton whr th avrag phas dffrnc btwn th oscllators changs from to narly π. Ths s smlar to th phas-flp bfurcaton rportd n th contxt of tmdlay coupld systms [36, 37]. Th natur of th transtons to th stat of ampltud dath s furthr charactrzd by fxng on of th paramtrs or ε d and ncrasng th othr. Ths s shown n Fg. 7, whr th ndx A s plottd for ncrasng for a chosn valu of ε d. Hr, th transton from th oscllatory to th ampltud dath stat s contnuous such that, as th couplng strngth s ncrasd, th ampltud of th oscllatons gradually dcrass to zro. A smlar transton s obsrvd for th cas whr s kpt fxd and ε d s ncrasd. W also notc from th tm srs and phas spac plot that, as th couplng strngth ncrass (ε d or ), th Rösslr systms undrgo th full rvrs proddoublng squnc to th on-cycl stat bfor gong to th ampltud dath stat. Thn, th transton to th stat of ampltud dath occurs va a suprcrtcal Hopf bfurcaton. Th bfurcaton dagram for ths transton s shown n Fg. 8. Ths s furthr confrmd by computaton of th fxd ponts of th coupld systm and thr stablty nar th transton rgon. Numrcal smulatons show that th coupld Rösslr systms n Eq.() stablz to th stady stat corrspondng to synchronzd stats of th subsystms. Ths synchronzd.5.5 FIG. 8: Bfurcaton dagram obtand by plottng th maxma of x 11 (x ) for suffcntly long prod by ncrasng th couplng strngth for ε d =. and κ = 1. Hr, th systms go through a rvrs prod-doublng bfurcaton as ncrass, ladng to a on-cycl (at.35) bfor ampltud dath occurs (c.45). stady stats ar obtand from Eq. () as x 1 = (c± c 4abκ/(κ ε a))/, x = x 1/a, x 3 = b/(x 1 c), y = x 1/κ. (16) Of th two fxd ponts, th on wth th plus sgn n th scond trm of th x 1 quaton s unstabl and th on wth th mnus sgn n th scond trm of x 1 quaton bcom stabl n th ampltud dath stat. Th natur of th transton to th stabl fxd pont s dtrmnd by th gnvalus of th corrspondng Jacoban, and w fnd that at th transton, ral parts of th complx conjugat pars of gnvalus bcom ngatv
7 Egn valu.8 -.8.4.43.46 Lyapunov xponnt. -. -.6.4.43.46 ε d 8 4 FIG. 9: Ral parts of th largst gnvalus of th Jacoban obtand from Eq. (), valuatd at th fxd pont gvn n Eq. (16), for ncrasng and fxd valus of ε d =. and κ = 1. At th transton (c.435), th ral parts of on par of complx conjugat gnvalus cross zro. Largst Lyapunov xponnt of th coupld Rösslr systm gvn n Eq. () for ncrasng for fxd valus of ε d =. and κ = 1. Th zro crossng of th largst Lyapunov xponnt (c.435) ndcats th transton to th ampltud dath stat. In both fgurs, zro s shown as a dottd ln. 8 16 [Fg. 9], ndcatng a suprcrtcal Hopf bfurcaton as dscrbd n Rf.[38]. In th ampltud dath rgon, all th Lyapunov xponnts of th coupld systm [ gvn n Eq. ] ar found to b ngatv. Th largst Lyapunov xponnt of th coupld systm crosss zro at th transton, and ths s shown n Fg. 9. Th natur of th transton s found to b th sam whn s kpt fxd and ε d s ncrasd. Th abov numrcal rsults ar prsntd for on st of paramtrs of th Rösslr systm. W hav vard th paramtrs and vrfd that th mthod works for othr valus of th paramtrs. 4.. Coupld Lornz systms W rpat th sam study n th cas of two coupld Lornz systms. It s ntrstng to not that, n ths cas, th coupld systms stablz to a fxd pont that corrsponds to ant-synchronzd stats for th subsystms (x 11 = x 1,x 1 = x,x 13 = x 3 ) as shown arlr n Fg. (c). Th rgons of dffrnt dynamcal stats n th paramtr plan of couplng strngths n ths cas ar shown n Fg. 1. Whn both ε d and ar small, th systms ar not synchronzd (wht rgon). For vry small valus of ε d and larg, th systms ar ant-synchronzd (dark-gray), and whn ε d s ncrasd from ths stat, th systms go to th ampltud dath stat (black rgon). For small valus of and larg ε d, th systms ar synchronzd (lght-gray). As ncrass, th systms frst los synchronzaton, and for largrvalusof, thystablztothstatofampltud dath (black). In th d-synchronzd stat bfor th ampltud dath stat, th attractor n th phas spac s hghly dstortd and th systm gos through a stat of frustraton, tryng to stablz to th ant-synchronzd stat from th synchronzd stat bfor dath occurs. Ths s llustratd n Fg. 11, whr th tm srs of th synchronzaton rror btwn th two Lornz systms s FIG. 1: Rgons of dffrnt dynamcal stats n th paramtr plan of couplng strngths ε d n th cas of Lornz systms. Th black rgon corrsponds to th stat of ampltud dath (A ), th lght-gray rgon corrsponds to th synchronzd stat (C 1), th dark-gray rgon corrspond to th ant-synchronzd stat (C 1) and th wht rgon corrspond to th stat whr C < 1. x 11 -x 1 15-15 11 15 19 Tm FIG. 11: Synchronzaton rror (x 11 x 1) as a functon of tm n th cas of Lornz systms for ε d = 5 and = 4.665. shown. Nar ths transton rgon, som ntal condtons rman n a chaotc transnt stat for a long tm bfor bcomng stablzd to th fxd pont. Th phnomna of mult-stablty and hystrss ar also obsrvd n ths rgon. Th natur of th transton to ampltud dath n ths cas s shown n Fg. 1. Unlk th cas of coupld Rösslr systms, hr w s that th ampltud of oscllatons drops suddnly at a crtcal strngth of couplng. Thus, th transton s drctly from th chaotc to th ampltud dath stat. A smlar typ of transton from th chaotc to th ampltud dath stat n th cas of tm-dlay coupld Lornz systms has bn rportd n Rf. [6]. W furthr charactrz ths transton by computng th fxd ponts of th coupld systm gvn n Eq. (3) and valuatng thr stablty nar th transton rgon. Numrcal smulatons show that th coupld Lornz systms n Eq. (3) stablz to th stady stat corrspondng to th ant-synchronzd stat of th subsys-
8 A 4 Egn valu.1 -.1 -. 7 7.5 8 8.5 Lyapunov xponnt 1-1 7 7.5 8 8.5 7.5 8 8.5 FIG. 1: Th ndx A as a functon of for a fxd valu of ε d = 5. for two coupld Lornz systms. Not that th transton to ampltud dath s suddn (c 7.88) as opposd th cas of coupld Rösslr systms whr th transton s contnuous. tms. Ths stady stats ar obtand from Eq. (3) as x ((r 1)σ ε d )b 11 = ±, σ +ε d x 1 = (σ +ε d) x σ 11, x 13 = (σ +ε d) x 11 σb, x 1 = x 11, x = x 1, x 3 = x 13, y =. (17) For both th solutons, w fnd that at th transton, th ral parts of th complx conjugat pars of gnvalus of th corrspondng Jacoban bcom ngatv [Fg. 13]. For ε d > ε dc, numrcally an unstabl lmt cycl s found to coxst wth th stabl stat of ampltud dath for crtan ntal valus. As thr s no stabl lmt cycl bfor ampltud dath, and an unstabl fxd pont bcoms stabl, t sms that ths s a sub-crtcal Hopf bfurcaton. All th Lyapunov xponnts of th systm [Eq. 3] ar found to b ngatv at th ampltud dath stat. At th transton, th largst Lyapunov xponnt of th coupld systm bcoms ngatv as shown n Fg. 13. FIG. 13: Ral parts of th largst gnvalus of th Jacoban obtand from Eq. (3), valuatd at th fxd pont gvn n Eq. (17), for ncrasng and fxd valus of ε d = 5 andκ = 1. Atthtranston (c 7.31), thral parts ofon par of complx conjugat gnvalus cross zro. Largst Lyapunov xponnt of th coupld Lornz systm gvn n Eq. (3)for ncrasng for fxd valus of ε d = 5 and κ = 1. Th zro crossng of th largst Lyapunov xponnt (c 7.43) ndcats transton to th ampltud dath stat. In both fgurs, zro s shown as a dottd ln. a varty of systms and fnd that, for all cass whch ar synchronzabl, th mthod works. W undrstand that xtnsv numrcal smulatons may not analytcally stablsh th gnralty of th mchansm. Howvr, th rsults of our numrcal smulatons appars to ndcat that ths mthod s qut gnral. W also gv an ntutv physcal argumnt to support our clam of gnralty of th mthod. Our mchansm conssts of havng two typs of couplng. Th frst s drct couplng, whch lads to synchronzaton. W know that, f th couplng constant s suffcntly larg, th synchronzaton condton nsurs that th largst Lyapunov xponnt transvrs to th synchronzaton manfold s ngatv. Consdrng th spac of coupld oscllators as a product of th ndvdual systms and a ntwork of two nods, th synchronzaton manfold corrsponds to th drcton s = (1,1) T n th ntwork coordnats. Smlarly, th couplng through th nvronmnt whch nsurs ant-synchronzaton lads to th condton that th largst Lyapunov xponnt transvrs to th drcton a = (1, 1) T s ngatv. Snc w hav coupld only two systms, nsurng that th largst xponnts transvrs to both synchronzng and ant-synchronzng drctons,.. s and a, ar ngatv, mpls that all th Lyapunov xponnts ar ngatv. Thus th systm must convrg to a fxd pont. In ths scton, w prsnt th rsults of applyng ths mthod to prodc, tm-dlay, hyprchaotc, and drvn systms and dffrnt schms of drct couplng. 5. AMPLITUDE DEATH IN OTHER CASES W hav prsntd a mchansm for nducng ampltud dath n coupld systms du to th comptng ffcts of synchronzng and ant-synchronzng tndncs. W hav shown ths n th contxt of two chaotc systms, namly, coupld Rösslr and coupld Lornz systms. Totst th gnraltyofth mthod, w applyt to 5.1. Ampltud dath n prodc systms W study two standard lmt cycl oscllators, namly Landau-Stuart and van dr Pol oscllators, coupld usng th schm gvn n Eq. (1). Th Landau-Stuart systm s a nonlnar lmt cycl oscllator, whch has bn prvously usd as a modl systm for studyng th phnomnon of ampltud dath [, 5]. In our cas, th
9 dynamcs of two coupld Landau-Stuart systms s gvn by th followng st of quatons: ẋ 1 = (x 1 +x )x 1 ωx +ε d (x j1 x 1 )+ y, ẋ = (x 1 +x )x +ωx 1, ẏ = κy ε x 1. (18) From numrcal analyss of th abov quatons wth ω =, w s that, for small valus of and ε d, th systms ar synchronzd. For small valus of ε d and larg valus of, th systms ar n th ant-synchronzd stat. Whn th strngths of both ε d and ar suffcntly larg, w obsrv ampltud dath. W not that hr th transton to th stat of ampltud dath s suddn, smlar to th cas of Lornz systms. Howvr, n ths cas, for a gvn strngth of couplng, th stablty of ampltud dath stat dpnds on th ntal condtons ndcatng mult-stablty, that s, som ntal condtons go to th ampltud dath stat, whl som othr ntal condtons rman n th oscllatory stat. Such a multstablty has also bn rportd for ampltud dath phnomna n th cas of Landau-Stuart oscllators usng conjugat couplng []. For som ntal valus, ampltud dath occurs vn n th absnc of drct couplng ( ε d = ). A possbl xplanaton s that th µ of th ndvdual systm s ngatv or zro such that, th stablty condton κ > µ gvn n Eq. (11) s always satsfd. It s sn that th ara of th basn of ampltud dath ncrass as ε d s ncrasd. W rpat th sam study for th cas of two coupld prodc van dr Pol systms gvn by th followng quatons: ẋ 1 = x +ε d (x j1 x 1 )+ y, ẋ = α(1 x 1 )x x 1, ẏ = κy ε x 1. (19) W choos th paramtr α = 1 such that th systm has a stabl lmt cycl whn both th couplngs ar absnt (.., = ε d = ). For a sutabl strngth of th drct couplng ε d, th systms synchronz and ampltud dath s obsrvd whn both drct and ndrct couplngs ar abov thr rspctv thrsholds. As far as th natur of th transton to th stat of ampltud dath s concrnd, th van dr Pol oscllator shows a dffrnt bhavor from that of othr systms. W fnd that th natur of th transton to ampltud dath dpnds on th typ of couplng paramtr. Fxng ε d and ncrasng, w s a smooth transton smlar to that n th cas of Rösslr systms and by fxng and ncrasng ε d, w gt a suddn transton, as n th cas of Lornz systms. 5.. Ampltud dath n tm dlay systms Th Macky-Glass tm-dlay systm s wll studd as a modl xhbtng hyprchaos. Stablzaton to fxd pont, or ampltud dath n such systms has bn rportd by us of stablzaton mthods such as convntonal fdback, trackng fltrs and dlayd fdback[39]. Hr w consdr two Macky-Glass systms coupld va both drct and ndrct couplngs as gvn by ẋ = αx + βx τ ẏ = κy 1+x m τ +ε d (x j x )+ y, x j, () j=1, whrx rprsntsthmacky-glasssystm[4]andx τ s th valu of th varabl x at a dlayd tm x (t τ). Th paramtrs of th Macky-Glass systms ar chosn such that, th ndvdual systms ar n th hyprchaotc rgm ( α = 1, β =, τ =.5, m = 1). For vry small valusofε d and, thtwosystmsarnotsynchronzd. Th systms synchronz as ε d s ncrasd. For larg valus of ε d, f s ncrasd, th systms go to a stat of ampltud dath. Whn ε d s small and s larg, th systms ar n an ant-synchronzd stat. As w ncras ε d, th systms go to th ampltud dath stat. Howvr, for larg valus of ε d and, th Macky-Glass systms appars to show a dffrnt bhavor from th Rösslr or Lornz systms. W obsrv a r-ntrant bhavor to rhythmognss, both as ncrass and as ε d ncrass. Ths transton also satsfs our stablty condtons Eqs. (13) and (14). W fnd that th transton to th stat of ampltud dath n th cas of two coupld Macky-Glass systms s contnuous and that th systms go through a rvrs prod-doublng squnc rachng a lmt cycl bfor th ampltud dath occurs. Ths s smlar to th cas of Rösslr systms dscussd arlr. 5.3. Ampltud dath wth altrnat schms for drct couplng So far, w hav studd drct couplng of th dffusv typ. Synchronzaton s also possbl wth drct couplng of dffrnt typs. Hr, w now study two such typs of couplng. 5.3.1. Lornz systms wth rplacmnt couplng Hr w consdr Lornz systms coupld usng a dffrnt schm of couplng, namly rplacmnt couplng,
1 as gvn by th followng quatons ẋ 1 = σ(x j x 1 )+ y, ẋ = (r x 3 )x 1 x, ẋ 3 = x 1 x bx 3, ẏ = κy ε x 1. (1) ε d 1.5 Hr, th drct couplng s of th rplacmnt typ, such that th x varabl n th frst functon of th frst systm s that of th scond systm, and vc vrsa. Ths typ of couplng lads to synchronzaton as rportd n Rf. [41]. W ntroduc ndrct couplng through th varabl y. W fnd that, for sutabl valus of th couplng strngth, th systms stablz to a stat of ampltud dath. 4 5.3.. Synaptcally coupld Hndmarsh-Ros modl of nurons Th Hndmarsh-Ros systm s a modl of nurons whch shows th spkng and burstng bhavor of th mmbran potntal of a sngl nuron [4]. W tak two nurons wth xctatory synaptc couplng [43] btwn thm and ntroduc an ndrct couplng as gvn by th followng quatons: ẋ 1 = x x 3 1 +ax 1 x 3 +I + V r x 1 y +ε d 1+xp( λ(x j1 θ)), ẋ = 1 bx 1 x, ẋ 3 = ρ(s(x 1 +χ) x 3 ), ẏ = κy ε x 1. () Hr, th varabl x 1 rprsnts th mmbran potntal of a nuron and th varabls x and x 3 ar rlatd to on currnts across th mmbran. W choos th paramtrs of th systm such that th ndvdual nurons ar n th chaotc burstng stat. Whn =, and ε d s suffcntly larg, th bursts of both nurons bcom synchronzd. For largr valus of ε d and, w obsrv patchs of ampltud dath (shown n Fg. 14). Thus as w kp on of th couplng paramtrs fxd ( or ε d ), and ncras th othr, w obsrv a transton to th ampltud dath stat and agan r-ntrant bhavor to spks. 5.4. Ampltud dath wth altrnat dynamcs for th nvronmnt In th prvous sctons, w hav takn th ntrnsc dynamcs of th nvronmnt to b that of an ovrdampd harmonc oscllator. Hr, w show that ampltud dath FIG. 14: Rgon of ampltud dath (black) n th paramtr plan of couplng strngths ε d for two coupld Hndmarsh- Ros nurons. W choos th paramtrs of th systm to b a = 3,b = 5, ρ =.6, s = 4, χ = 1.6, and I = 3.. Th paramtrs n th couplng trms ar chosn to b V r = 3, θ =.5, λ = 1, and κ = 1. For th nuronal systms, th crtron usd n th txt for dntfyng th ampltud dath stat,.., th ndx A, nds to b changd. Hr, w us th maxmum voltag from both nurons (x m) as th ndx to dntfy th dath stat. Th crtra for synchronzd and ant-synchronzd stats also nds modfcaton. Hnc, th othr rgons ar not shown xplctly n th fgur. s possbl wth othr ntrnsc dynamcs for th nvronmnt also. For ths, w consdr th cas of two Rösslr systms coupld wth a dynamc nvronmnt, whr th ntrnsc dynamcs of th nvronmnt s that of a dampd harmonc oscllator. Th dynamcs of th coupld systm n ths cas s gvn by ẋ 1 = x x 3 +ε d (x j1 x 1 )+ y 1, ẋ = x 1 +ax, ẋ 3 = b+x 3 (x 1 c), ẏ 1 = y ε x 1, ẏ = y κy 1. (3) Hr, th varablsy 1 and y rprsntatwo-dmnsonal nvronmnt togthr formng an undrdampd harmonc oscllator. For vry wak couplng (ε d, ), th two Rösslr systms ar not synchronzd. Whn th couplng strngth ε d s ncrasd whl s kpt fxd at zro, th systms bcom synchronzd. On th othr hand, whn s ncrasd whl ε d s kpt fxd at zro, th systms bcom ant-phas-synchronzd. Whn both and ε d ar abov thr rspctv thrsholds, ampltud dath s obsrvd. W rpat th sam study by takng th ntrnsc dynamcs of th nvronmnt as that of an ovrdampd
11 Duffng oscllator. Th quatons n ths cas ar ẋ 1 = x x 3 +ε d (x j1 x 1 )+ y, ẋ = x 1 +ax, ẋ 3 = b+x 3 (x 1 c), ẏ = y κy 3 ε x 1. (4) For small valus of and ε d, w s that th systms ar not synchronzd. Thy bcom synchronzd as ε d s ncrasd from ths stat. For small valus of ε d and larg valus of, th systms ar n an ant-phassynchronzdstat. Whn thstrngthsofbothε d and ar suffcntly larg, w obsrv ampltud dath. Th phas dagram n ths cas s qualtatvly smlar to that gvn n Fg. 6. 5.5. Ampltud dath n hyprchaotc Rösslr systms W also consdr th cas of two hyprchaotc Rösslr systms as gvn by th followng quatons ẋ 1 = x x 3 + y 1, ε d cosθ(snθ(x 1 x j1 )+cosθ(x 3 x j3 )), ẋ = x 1 +ax +x 4 + y, ẋ 3 = b+x 3 x 1 ε d snθ(snθ(x 1 x j1 )+cosθ(x 3 x j3 )), ẋ 4 = cx 3 +σx 4 + y 3, ẏ 1 = κy 1 ẏ = κy ẏ 3 = κy 3 x 1, x, x 4, (5) whr,j = 1, and j. W choos th paramtrs of th systm such that, th ntrnsc dynamcs of th systms s hyprchaotc. For ths systm, th mthod of tm-dlay couplng s found nffctv for producng ampltud dath[3]. In th absnc of couplng va th nvronmnt, th drct couplng va a scalar sgnal rsults n synchronzaton of th two hyprchaotcsystms for sutabl valus of paramtrsε d and θ, as rportd n Rf. [44]. W tak th nvronmnt to b thr-dmnsonal n ths cas. W fnd that ampltud dath occurs n th coupld systm for sutabl valus of th couplng strngths and κ. Ths s shown n Fg. 15. 5.6. Small oscllatons n drvn systms W apply th schm dscrbd n ths papr to drvn systms such as drvn van dr Pol and Duffng systms. x 1 1-1 15 3 Tm FIG. 15: Tm srs of th frst varabls x 1, = 1, of th two coupld hyprchaotc Rösslr systms gvn n Eq. (5) n th ampltud dath stat. Paramtrs of th systms ar chosn to b a =.5, b = 3, c =.5, and σ =.5. Th paramtrs usd n th drct couplng trm ar ε d =.5, and θ = π/3 and n th ndrct couplng trms ar κ = 1 and = 4. In such drvn systms, th fxd pont s not a soluton for th ndvdual or coupld systms. So th ntrprtaton of ampltud dath as n othr systms nds to b changd. Hr, w ntrprt th ampltud dath stat as th stat of vry small ampltud oscllatons. Drvn van dr Pol systms wth drct dffusv couplng and ndrct couplng through th nvronmnt can b wrttn as ẋ 1 = x +ε d (x j1 x 1 )+ y, ẋ = α(1 x 1)x x 1 +βcos(ωt), ẏ = κy ε x j1. (6) j W fnd that, whn = and ε d ncrass, th two systms bcom synchronzd. In ths stat, f w start ncrasng, w gt a stat of small oscllatons or ampltud dath. Fgur 16 shows th tm srs for such a stat. Fgur 17 plots th ndx A as a functon of. W frst s a transton from a lmt cycl to two dffrnt lmt cycls for th two systms. Ths stat subsquntly gos to th ampltud dath stat contnuously as ncrass furthr. On th othr hand, f w kp fxd and ncras ε d, w fnd a suddn transton to th ampltud dath stat. As ε d ncrass furthr, w fnd a contnuous transton to a stat of ncrasng ampltud oscllatons. Ths s shown n Fg. 17. Th sam study s rpatd for th cas of two coupld Duffng systms gvn by th followng quatons: ẋ 1 = x +ε d (x j1 x 1 )+ y, ẋ = αx +x 1 x 3 1 +βcos(ωt), ẏ = κy ε x j1. (7) j W fnd that a rgm of small oscllatons s possbl n ths cas also. Ths s shown n Fg. 16. Thus, w hav llustratd that th mthod for nducng ampltud dath n coupld systms ntroducd n ths papr works for prodc, tm-dlay, hyprchaotc, and
1 1.5 4 x 11, x 1 A -1.5 1. 3 4 5.75 3.5 x 11, x 1 4 A -1. 5 1 Tm FIG. 16: Tm srs of th frst varabl x 1 of th coupld drvn systms showng small oscllatons. Hr, w ntrprt th small oscllatons as th stat of ampltud dath (s txt). Drvn van dr Pol systm for (ε d,) = (1.,3.5). Th paramtrs of th ndvdual systms ar takn to b α = 8.53, ω =.63, and β = 1. Drvn Duffng systm for (ε d,) = (1.,4.). Th paramtrs of th ndvdual systms ar takn to b α =.5, ω = 1, and β =.3 drvn systms. It s ffctv n qunchng th dynamcs vn wth dffrnt forms of drct couplng, such as rplacmnt couplng and synaptc couplng of nurons. Th contxt of coupld nurons prsnts an mportant cas of ampltud dath whch could xplan th mchansm of dsrupton or supprsson of synaptc sgnals n th cas of nuronal dsordrs, such as Alzhmr s dsas, as bng du to nducd actvty and fdback through a protn calld amylod bta (Aβ). Usng numrcal studs on a fw nuronal modls, w hav shown that th comptng ffcts of synaptc actvty and th ndrct ntracton mdatd by th protn Aβ lad to sub-thrshold actvty and synaptc slncng [45]. 6. DISCUSSION In ths papr, w show that ndrct couplng through a dynamc nvronmnt n addton to drct couplng can lad to ampltud dath n chaotc systms such as Rösslr and Lornz systms. Th approxmat stablty analyss dvlopd for gnral cass gvs th transton rgon n paramtr spac whch s furthr supportd by.5 1 1.5 ε d FIG. 17: Transton to ampltud dath stat n two coupld drvn van dr Pol oscllators. Hr th dath stat corrsponds to a stat of small oscllatons snc th fxd pont s not a stabl stat. Th ndx A rmans fnt although vry small n th ampltud dath stat. Th ndx A as a functon of for fxd ε d = 1.. Th transton s contnuous (c 3.8). Thr s a suddn transton du to a ptchfork bfurcaton of th lmt cycl obsrvd at.55. Th ndx A as a functon of ε d for fxd = 3.5. Th transton to th ampltud dath stat s suddn (ε dc.7). W also obsrv a r-ntrant contnuous transton to prodc oscllatons wth ncrasng ampltud at ε d 1.14. drct numrcal smulatons. Th natur of th transton to ampltud dath s found to b of two typcal typs, on contnuous and th othr dscontnuous. In concluson, th mthod for dtrmnng ampltud dath ntroducd n ths papr s qut gnral and works for dffrnt typs of systms such as prodc, chaotc, hyprchaotc and tm-dlay systms and also wth dffrnt typs of drct couplng, for xampl, dffusv, rplacmnt couplng, couplng va scalar sgnals and synaptc couplng. W hav dmonstratd that our mthod works n hyprchaotc Rösslr systms whr tm-dlay couplng s nffctv n nducng ampltud dath. W hav also prsntd a physcal argumnt for th gnralty of th mthod. Although th thory s dvlopd for th cas of dntcal systms, w hav vrfd that ths mthod of nducng ampltud dath works
13 n th cas of non-dntcal systms as wll. In fact, w dd not fnd any xcpton to our schm provdd th coupld systms ar synchronzabl. W also not that th mthod ntroducd n ths papr can b mplmntd n practcal cass. What s ndd s th dsgn of a sutabl nvronmnt whch can ntroduc th approprat ndrct couplng btwn th systms. Morovr, n many natural systms, th nvronmnt or xtrnal mdum xsts and can b nstrumntal n causng supprsson of th dynamcs. Th mchansm prsntd hr provds an xplanaton for ths phnomnon n such cass. Acknowldgmnts W thank th rfrs for thr constructv commnts and suggstons. [1] A. S. Pkovsky, M. G. Rosnblum, and J. Kurths, Synchronzaton: A Unvrsal Concpt n Nonlnar Scncs, Cambrdg Nonlnar Scnc Srs, (Cambrdg Unvrsty Prss, London, 3). [] K. Kanko, Thory and Applcatons of Coupld Map Lattcs (Wly, Nw York, 1993). [3] E. Ott, Chaos n Dynamcal Systms (Cambrdg Unvrsty Prss, Cambrdg, 1993). [4] M. Dolnk and M. Mark, J. Phys. Chm. 9, 45 (1988). [5] M. F. Crowly and I. R. Epstn, J. Phys. Chm. 93, 496 (1989). [6] K. Bar-El, Physca D 14, 4 (1985). [7] M. Dolnk and I. R. Epstn, Phys. Rv. E 54, 3361 (1996). [8] K. Tsanva-Atanasova, C. L. Zmlk, R. Brtram, and A. Shrman, Bophys. J. 9, 3434 (6). [9] A. Kosska, E. Volkov, and J. Kurths, Chaos, 313 (1). [1] I. Ozdn, S. Vnkataraman, M. A. Long, B. W. Connors, and A. V. Nurmkko, Phys. Rv. Ltt. 93, 1581 (4). [11] R. Hrrro, M. Fguras, J. Rus, F. P, and G. Orrols, Phys. Rv. Ltt 84, 531 (). [1] M. W and J. Lun, Appl. Phys. Ltt. 91, 6111 (7). [13] J. Bnford, H. Sz, W. Woo, R. R. Smth, and B. Hartnck, Phys. Rv. Ltt. 6, 969 (1989). [14] M. Y. Km, R. Roy, J. L. Aron, T. W. Carr, and I. B. Schwartz, Phys. Rv. Ltt. 94, 8811 (5). [15] P. Kumar, A. Prasad, and R. Ghosh, J. Phys. B 41, 1354 (8). [16] D. J. Slko, Ann. N. Y. Acad. Sc. 94, 17 (). [17] R. E. Tanz, Nat Nurosc. 8 977 (5). [18] B. Caughy, and P.T. Lansbury Jr., Annu. Rv. Nurosc. 6, 67 (3). [19] R. E. Mrollo and S. H. Strogatz, J. Stat. Phys. 6, 45 (199). [] B. Ermntrout, Physca D 41, 19 (199). [1] Y. Zha, I. Z. Kss, and J. L. Hudson, Phys. Rv. E 69, 68 (4). [] R. Karnatak, R. Ramaswamy and A. Prasad, Phys. Rv. E 76, 351 (7). [3] M. Dasgupta, M. Rvra, and P. Parmananda, Chaos, 316 (1). [4] K. Konsh, Phys. Rv. E 68, 67 (3). [5] D. V. Ramana Rddy, A. Sn, and G. L. Johnston, Phys. Rv. Ltt. 8, 519 (1998). [6] A. Prasad, Phys. Rv. E 7, 564 (5). [7] C. U. Cho, V. Flunkrt, P. Hövl, H. Bnnr, and E. Schöll, Phys. Rv. E 75, 466 (7). [8] R. Dodla, A. Sn, and G. L. Johnston, Phys. Rv. E 69, 5617 (4). [9] F. M. Atay, Phys. Rv. Ltt. 91, 9411 (3). [3] Y. Chn, J. Xao, W. Lu, L. L, and Y. Yang, Phys. Rv. E 8, 466 (9). [31] K. Konsh, Phys. Rv. E 7, 661 (4). [3] K. Konsh, Phys. Ltt. A 341, 41 (5). [33] V. Rsm, G. Ambka, and R. E. Amrtkar, Phys. Rv. E 81, 4616 (1). [34] G. Ambka and R. E. Amrtkar, Phys. Rv. E 79, 566 (9). [35] A. Wolf, J. B. Swft, H. L. Swnny, and J. A. Vastano, Physca D 16, 85 (1985). [36] A. Prasad, S. K. Dana, R. Karnatak, J. Kurths, B. Blasus, and R. Ramaswamy, Chaos 18, 3111 (8). [37] R. Karnatak, N. Puntha, A. Prasad, and R. Ramaswamy, Phys. Rv. E 8, 4619 (1). [38] S. H. Strogatz, Nonlnar dynamcs and chaos: wth applcatons to physcs, bology, chmstry, and ngnrng, ( Wstvw Prss, 1994). [39] A. Namajūnas, K. Pyragas, and A. Tamašvčus, Phys. Ltt. A 4, 55 (1995). [4] M. C. Macky, and L. Glass, Scnc 197, 87 (1977). [41] J. Güémz and M. A. Matías, Phys. Rv. E 5, R145 (1995). [4] J. L. Hndmarsh and R. M. Ros, Proc. R. Soc. London, Sr. B 1, 87 (1984). [43] I. Blykh, E. d Lang, and M. Haslr, Phys Rv Ltt 94, 18811 (5). [44] J. H. Png, E. J. Dng, M. Dng, and W. Yang, Phys. Rv. Ltt. 76, 94 (1996). [45] V. Rsm,G.Ambka, R.E.Amrtkar and G.Rangarajan, (unpublshd).