attractors for dynamical systems are traditionally viewed as being asymptotically stable invariant sets which have a neighbourhood that absorbs all sufficiently close initial conditions . as a consequence
the basins of such attractors contain open sets .
whereas this concept works very well for simple attractors of ordinary differential equations , for more complicated attractors the smallest asymptotically stable attractor may contain much more than the asymptotic dynamics of ` typical ' initial conditions .
partly in response to this problem , milnor introduced a concept of measure attractor ( now usually called _ milnor attractor _ ) @xcite .
this is a compact invariant set whose basins of attraction has positive measure in phase space .
examples of riddled basins @xcite show that , even for smooth invertible dynamics , one can find milnor attractors with riddled basins ; i.e. such that any open set that intersects the basin has positive measure in the basin of a different attractor ; see @xcite for a recent review of riddled basins .    in the absence of smoothness or invertibility , it is clear that more exotic attractors may appear .
recent work @xcite
on globally coupled networks of oscillator with delayed pulse interactions @xcite indicates that extreme cases of riddled basin attractors , _ unstable attractors _ , can appear , where there is a neighbourhood @xmath1 of the attractor such that almost all points exit from @xmath1 under the dynamics .
even more surprisingly , numerical simulations @xcite seem to indicate that these attractors can appear in a robust way as long as there are sufficiently many oscillators in the system .
the main aim of this paper is to give a rigorous explanation for the appearance and robustness of unstable attractors in such systems .
this paper is organized as follows .
section  [ secunstableatts ] examines the general problem of unstable attractors in a systematic way .
we give some necessary conditions for the appearance of unstable attractors and some motivating examples .
subsequently , we address the relationship between networks of unstable attractors and robust heteroclinic cycles , and present a simple recipe of how to perturb a smooth flow with a robust heteroclinic cycle to a smooth semiflow with a network of unstable attractors .
in section  [ secpulsecos ] we consider the system of globally pulse coupled oscillators with delayed interactions studied in @xcite , including a generalization of the set of possible response functions .
because the system combines continuous time evolution with discrete time events , it is strictly speaking a hybrid system @xcite with delay .
the main result of this section is theorem  [ thmfinite ] which shows that the dynamics reduces , after a finite time , to a finite dimensional hybrid system .
it gives an explicit upper bound on the number of remaining dimensions in the system .
in section  [ secuasforpulsecos ] we prove that the system discussed in section  [ secpulsecos ] exhibits unstable attractors ( theorem  [ thm : uasexist ] ) and that these persist for an open set of parameter values ( theorem  [ thm : uasrobust ] ) .
the proof has two parts ; firstly we show ( by following an appropriate open set of initial conditions ) that a certain periodic orbit is a milnor attractor ; secondly we show that this orbit is a saddle .
we discuss this mechanism as being an effect of ` dimension jump ' in the system .
section  [ secdiscuss ] discusses obstacles to proving the existence and robustness of unstable attractors in a more general setting . for convenience some of the details of the proofs have been placed in appendices : appendix  [ app : k=1 ] describes the system in the case that the delay is short enough that only one delayed pulse can influence the future dynamics of the system .
appendix  [ app : unstable ] derives the return map on a section transverse to the unstable attracting periodic orbit .
consider a dynamical system defined by a semiflow @xmath2 on a finite dimensional manifold @xmath3 for @xmath4 ( recall that a semiflow is a family of maps such that @xmath5 whenever @xmath6 , @xmath7 , and such that @xmath8 ) . in
what follows we will consider continuous ( @xmath9 ) or discrete ( @xmath10 ) dynamics .
our examples arise from cases where @xmath11 is not invertible for all @xmath12 : one can not extend the semiflow to a flow .
for convenience we will assume that @xmath3 is compact and let @xmath13 denote lebesgue / riemann measure on @xmath3 .
recall that the @xmath14-limit set of @xmath15 is defined by @xmath16 a _ milnor attractor _ for @xmath17 is a compact invariant subset @xmath18 such that ( i ) the basin of attraction @xmath19 has positive measure in @xmath3 , i.e. @xmath20 , and ( ii ) any compact invariant proper subset of @xmath0 has a basin with strictly smaller measure @xcite .
it is a _ minimal milnor attractor _ if any compact invariant proper subset of @xmath0 has a basin with zero measure .
all explicit examples of milnor attractors given in this paper are fixed points or periodic orbits and thus they are minimal milnor attractors because they do not contain any compact invariant proper subset .
however , we will not explicitly discuss minimality in these cases .
given any subset @xmath21 we define the _ lingering subset _ of @xmath1 to be @xmath22 [ def : lingeringset ]    the lingering set is a ( forward ) invariant set consisting of points in @xmath1 which do not leave @xmath1 in the future ; note that if @xmath0 is an invariant set within @xmath1 then @xmath23 is non - empty as it contains @xmath0 .
[ def : ua ] we say a milnor attractor @xmath0 is an _ unstable attractor _ if there is a neighbourhood @xmath1 of @xmath0 with @xmath24    requirement ( [ eq : nocontractingbasin ] ) means that almost all trajectories in a neighbourhood of the attractor must leave this neighbourhood .
we distinguish between two classes of unstable attractors .
[ def : uapositivelocalbasin ] an unstable attractor @xmath0 is an _ unstable attractor with positive measure local basin _ ( resp .
[ def : uazerolocalbasin ] ) . ] if @xmath25 for all neighbourhoods @xmath1 of @xmath0 .
[ def : uazerolocalbasin ] an unstable attractor @xmath0 is an _ unstable attractor with zero measure local basin _ if there is a neighbourhood @xmath1 of @xmath0 with @xmath26    for an unstable attractor with positive local basin , ( [ eq : nocontractingbasin ] ) means that almost all trajectories starting in a sufficiently small neighborhood @xmath1 will first leave @xmath1 .
a positive measure will eventually return and be asymptotic to @xmath0 ( contrasting with with @xcite where there is a positive measure that remains in @xmath1 and is asymptotic to @xmath0 ) .
for an unstable attractor @xmath0 with zero local basin almost all initial conditions in a sufficiently small neighbourhood @xmath1 of @xmath0 will leave @xmath1 and never limit to @xmath0 .
has positive measure if @xmath27 ; it has full measure ( in @xmath1 ) if @xmath28 and @xmath29 . ]    clearly , linearly stable fixed points or periodic orbit attractors for a smooth invertible system can not be unstable attractors ; in such cases one can construct a neighbourhood that by the hartman - grobman theorem is contained within the basin of attraction . if @xmath0 is asymptotically stable then trivially it can not be an unstable attractor .
a more general observation is the following proposition characterizing a class of attractors that can not be unstable with zero local basin .
suppose that @xmath0 is a milnor attractor whose basin @xmath30 contains an open set , and suppose that @xmath17 is almost everywhere a local homeomorphism
. then @xmath0 can not be an unstable attractor with zero local basin .
[ prop : nounsts ]    * proof : *  consider a point @xmath31 in the interior of the basin of @xmath0 such that @xmath11 is a local homeomorphism in a neighbourhood of @xmath31 , and choose an open neighbourhood @xmath1 of @xmath0 . there is a finite time @xmath32 after which @xmath33 . pick a small neighbourhood @xmath34 of @xmath35 in @xmath1 ;
continuity of @xmath36 means that there is a neighbourhood @xmath37 of @xmath31 within the basin of attraction such that @xmath38 .
continuity of the local inverse means that @xmath39 is also a neighbourhood of @xmath35 .
however all points in @xmath37 are in the basin of @xmath0 and so @xmath39 is an open set in the basin of @xmath0 that is within @xmath1 .
hence @xmath40 for any @xmath1 that is a neighbourhood of @xmath0 .
@xmath41    in particular , this proposition ensures that continuous invertible flows with continuous inverse ( homeomorphisms ) do not exhibit unstable attractors with zero local basin .
we have so far been unable to determine whether the assumptions of proposition  [ prop : nounsts ] exclude the possibility that @xmath0 is an unstable attractor with positive measure local basin , though clearly under additional assumptions ( such as linear stability of @xmath0 ) we can exclude this case .
as we will see below , unstable attractors occur robustly for certain classes of noninvertible dynamical systems .
suppose @xmath42 is a space of dynamical systems under some suitable topology , such that one can obtain an unstable attractor for some @xmath43 .
if there are unstable attractors for all @xmath44 in a neighbourhood of @xmath45 then we say the unstable attractors for @xmath45 are _ robust _ ; this means that they are stable to perturbations of the system within this @xmath42 . in section  [ secuasforpulsecos ]
we show robustness of unstable attractors with respect to perturbations of a pulse response function defining a class of pulse coupled oscillator systems .
there is no simple relationship between the lingering set for a given neighbourhood of an attractor and the basin of attraction ; in particular there may be points that are in the lingering set and not in the basin , and vice - versa .
for example , near a centre in a planar system the lingering set can be an open set but the basin just one point .
similarly the flow of @xmath46 on the unit circle @xmath47 has a fixed point @xmath48 .
the basin of the attractor @xmath49 is the whole of phase space , but the lingering sets of small open neighbourhoods of @xmath48 will not contain nearby points to the right of @xmath48 .
unstable attractors are a related but different phenomena to riddled basin attractors @xmath0 @xcite .
recall that an attractor @xmath0 has _ riddled basin _ if @xmath30 is such that @xmath50 for all open sets @xmath1 that intersect @xmath30 .
in particular , any neighbourhood @xmath1 that contains @xmath0 will intersect @xmath30 in a set of positive measure ; this means that a positive measure set in the basin leaves @xmath1 . for unstable attractors
almost all points in a small enough @xmath1 will leave @xmath1 ; by the comment above , this does not have any direct implication for the basin of attraction of @xmath0 , though it does imply that any ` local basin of attraction ' relative to @xmath1 will have zero measure .
examples of riddled basin attractors so far all possess irregular or chaotic dynamics to provide attraction local to some points and repulsion local to others .
by contrast the dynamics observed on unstable attractors can be very simple , even equilibrium .      to illustrate the possible types of unstable attractors we consider the dynamical systems defined by the piecewise continuous mapping @xmath51\mapsto[0,1]$ ] according to @xmath52 and @xmath53\rightarrow[0,1]$ ] where @xmath54 except on @xmath55 in which case @xmath56 .
these two maps are illustrated in figure  [ figeg1 ] .
one can verify that the only invariant sets for @xmath45 and @xmath44 are the two fixed points @xmath48 and @xmath57 .
these are both minimal milnor attractors ; for @xmath45 their basins of attraction are respectively @xmath58}.\label{eq : basinf(0)}\ ] ] whereas for @xmath44 they are @xmath59.\label{eq : basing(0)}\ ] ] in both cases the disjoint union of the two basins is the interval @xmath60 $ ] . on the other hand , any open set @xmath61 that contains @xmath62
has lingering set @xmath23 that is just the point @xmath63 .
hence in both cases this point is an unstable attractor . considering the definitions ( [ eq : ualocalbasin ] ) and ( [ eq : uanolocalbasin ] ) above
, we see that the attractor @xmath64 is an unstable attractor with zero local basin for @xmath45 and an unstable attractor with positive local basin for @xmath44 .     and @xmath44 exhibit an unstable attractor @xmath65 with ( a ) zero and ( b ) positive local basin .
[ figeg1],width=302 ]      as noted in @xcite , unstable attractors can form networks such that a positive measure of sufficiently small perturbations from the unstable attractor @xmath66 are in the basin of attraction of some other unstable attractor @xmath67 .
we show that this is analogous to robust heteroclinic networks .    a set of unstable attractors @xmath66 , @xmath68 forms a _ network _ given by the directed graph with vertices @xmath66 .
there are directed edges from @xmath66 to @xmath67 if and only if for any neighbourhood @xmath69 of @xmath66 we have @xmath70    if in addition @xmath71 we say the network is _
closed_. all trajectories near a closed network can not leave a neighbourhood of the network .
in such a case the dynamics of the system is well described by the network even in the presence of small perturbations to the system .
we illustrate such a network schematically in figure  [ fignetworkofuas ] .
if a network is not closed , we only have @xmath72 and such a network of unstable attractors @xmath73 may still manifest itself in transients .
if there is a proper inequality in ( [ eq : notclosed ] ) then a positive measure of initial conditions starting near the network are not asymptotic to states that are in the network .
a network is called _ transitive _ if there is a directed path between any two attractors .
an arbitrary closed network will contain subnetworks that are transitive networks of unstable attractors .     may form a network where a link from @xmath73 to @xmath74 implies that the basin @xmath75 intersects every neighbourhood of @xmath73 in a set of positive measure . in
the example sketched above , there is a transitive subnetwork @xmath76 .
[ fignetworkofuas ] , width=151 ]      we briefly demonstrate that one can obtain unstable attractors from smooth semiflows with saddles by including small perturbations that induce collapse onto their stable manifolds .
consider a hyperbolic saddle equilibrium @xmath77 for a ( non - invertible ) semiflow @xmath2 on a compact @xmath3 with stable manifold @xmath78 .
we say @xmath17 exhibits _ collapse onto the stable manifold _ of @xmath79
if there is an open set @xmath80 and a @xmath81 such that @xmath82 . since @xmath83 is a saddle , and has thus an unstable manifold , such an @xmath1 can not contain @xmath83 .
similarly , one can define _ collapse onto the unstable manifold _ of a saddle @xmath83 .
this can occur on a neighbourhood that includes @xmath83 and in this case @xmath1 may contain @xmath83 .
an example of an unstable attractor constructed by using a pair of saddle equilibria is shown in figure  [ fig : twosaddles ] ; this is essentially the mechanism that causes the appearance of unstable attractors in the system we discuss in section  [ secpulsecos ] .
if there is a connecting orbit from @xmath84 to @xmath83 then collapse onto the unstable manifold of @xmath84 is simultaneously a collapse onto the stable manifold of @xmath83 ; the points in the open set @xmath34 intersecting the unstable manifold of @xmath85 tend to @xmath79 .
unstable attractors in this sort of scheme can be robust to perturbations of the system if the connecting orbits are in invariant subspaces .
is an unstable attractor .
the set @xmath86 is a heteroclinic connection from @xmath84 to @xmath83 .
we assume that the open set @xmath34 in the neighbourhood @xmath87 of another equilibrium @xmath84 collapses in finite time onto the stable manifold @xmath88 of @xmath79 .
hence @xmath89 ; moreover , almost all perturbations to @xmath83 within @xmath90 will leave @xmath90 and so @xmath83 is an unstable attractor . [
fig : twosaddles ] , width=302 ]      smooth dynamical systems with symmetries can exhibit heteroclinic cycles or homoclinic cycles as robust milnor attractors , if the connecting orbits lie within symmetry - forced invariant subspaces ; see for example @xcite .
we indicate briefly how one can use these to construct a large class of systems that have networks of unstable attractors .
consider a smooth semiflow @xmath2 on @xmath3 with an absorbing open region @xmath91 such that the only milnor attractor in @xmath92 consists of the set @xmath93 for hyperbolic saddle equilibria @xmath94 and the unstable manifolds @xmath95 .
we assume that the unstable manifolds are one dimensional and write the connecting orbits from @xmath94 to @xmath96 by @xmath97 given such an attractor @xmath98 , we can construct piecewise smooth semiflows @xmath99 such that the set of @xmath94 forms a network of unstable attractors , and such that @xmath17 and @xmath100 are identical except on a set of arbitrarily small measure .    to do this , pick a section @xmath101 to each connection @xmath102 and a neighbourhood @xmath103 of @xmath101 of arbitrarily small measure .
define @xmath100 to be a semiflow that is equal to @xmath17 when the trajectories are outside of the @xmath103 .
on entering the @xmath69 we require that the semiflow collapses onto the connecting orbit in finite time , as illustrated in figure  [ figuasfromhets ] .
the illustration depicts this for one - dimensional connections but there can clearly be higher - dimensional connections with qualitatively the same features .
can be locally perturbed in sets of arbitrarily small measure to obtain ( b ) a network of unstable attractors with the same connection topology . here
, trajectories entering the balls hit one of the connecting orbits after a finite time and continue out after a uniformly bounded time .
[ figuasfromhets ] , width=302 ]    note that robust networks of unstable attractors can appear under the same conditions as for robust heteroclinic attractors @xcite by adding collapse onto _ unstable _ manifolds because in this case unstable manifolds are contained within stable manifolds of other saddles .
this can arise in cases where there is a discontinuity in the semiflow , but we note that it is not necessary ; one could construct an unstable attractor similarly to figure  [ fig : twosaddles ] for a smooth semiflow , as long as this maps some open set onto a zero measure subset of the stable manifold of @xmath85 .
in this section we consider the system of oscillators with delayed pulse interactions investigated in @xcite .
first we will show that this system , although formally infinite dimensional , can be reduced under general assumptions to a finite dimensional system .
the number of dimensions remaining depends on the system parameters .
the system of oscillators we consider is technically a hybrid dynamical system .
a hybrid dynamical system is defined by a smooth flow in some regions of phase space with boundaries ; when the flow hits one of the boundaries it is then mapped instantaneously by a deterministic mapping to some other point on the boundary . a trajectory from a hybrid dynamical system
hence consists of intervals of continuous evolution interrupted by events that may be discontinuous .
an evolution of the system is a _ zeno state _
@xcite ( also referred to as a _ chattering state _ ) if there can be an infinite number of applications of the map in some finite time interval ; in such cases one needs to have extra rules to permit evolution beyond the first accumulation point in time of an infinite number of applications of the map .
we will show that there are no zeno states in the system considered below .
[ def : phasehistory ] we define @xmath104 to be the set of piecewise smooth ` phase history ' functions @xmath105\rightarrow[0,1)$ ] , @xmath106 such that :    * _ well defined limit from below .
_ for some discrete ( possibly finite ) , monotone increasing set of @xmath107 , @xmath108 we have @xmath109 is smooth for @xmath110 . *
_ no zeno states .
_ if there are infinitely many discontinuities @xmath111 then @xmath112 .
note that ( i ) means that the limit from below @xmath113 is defined for all @xmath114 . condition ( ii ) means that we allow no zeno states in the past .
the system @xcite of @xmath115 pulse - coupled oscillators with delay @xmath116 , where the phase histories are @xmath117 , is defined in terms of a semiflow on @xmath118 for all @xmath119 @xmath120 such that @xmath121 for all @xmath122 . on the interval
@xmath123 $ ] the function @xmath124 is defined below .
the boundaries of the intervals on which the phase history function @xmath125_i$ ] of oscillator @xmath126 is smooth are now denoted by @xmath127 , @xmath128 .
if @xmath129 we say oscillator @xmath126 _ fires _ at time @xmath130 .
given that the system has evolved a time @xmath131 from some initial condition @xmath132 at time @xmath133 , we list the past firings @xmath134 of oscillator @xmath126 as @xmath135 such that the @xmath136 are the times since the @xmath137th last firing of the @xmath126th oscillator before or at time @xmath131 .
note that the @xmath138 , @xmath139 count the firing events whereas the @xmath140 , @xmath139 count all discontinuities that may occur .
we write the phase of the oscillator @xmath126 as @xmath141 , namely as the most recent point in the phase history of that oscillator , and the previous firing times as @xmath138 .
the future dynamics will be determined only by the present phases and the times of the past firings ; formally there is a map @xmath142 such that @xmath143    thus the map @xmath144 directly gives the current phases @xmath145 and the past firing times @xmath136 , @xmath108 , of the oscillators @xmath146 .
similarly , it will be useful to define a map onto the phases and the @xmath147 past firing times : @xmath148      the semiflow @xmath149 in ( [ eq : semiflowforoscs ] ) is given as a hybrid system with continuous time evolution @xmath150 for all @xmath126 and @xmath147 , except when one of a number of possible discrete events occur .
these events are determined by a constant delay @xmath151 between one oscillator being reset and ` sending ' a pulse ( ` firing ' ) and its ` reception ' by another oscillator .
these events ` firing ' and ` reception ' are the only types of discrete events interrupting the continuous time evolution .
we say that at a time @xmath32    * _ oscillator @xmath126 fires _ , + denoted @xmath152 , if @xmath153 . *
_ oscillator @xmath126 receives a pulse from oscillator @xmath137 _ , + denoted @xmath154 , if @xmath155 for some @xmath147 and @xmath156 . * _ oscillator @xmath126 fires induced by reception of a pulse from oscillator @xmath137 _ , + denoted @xmath157 , if @xmath154 and this results in an immediate firing of oscillator @xmath126 ( see below ) .
observe that for any @xmath126 it is not possible that @xmath158 and @xmath154 ( for any @xmath137 ) occur simultaneously , whereas if @xmath154 occurs it is possible that @xmath159 , @xmath160 also occurs at the same time .
the phases @xmath161 and firing times @xmath162 , and thus the evolution @xmath163 occurs as ( [ eq : phisdotone ] ) until one of @xmath158 , @xmath164 or @xmath165 occur at time @xmath166 . by definition  [ def : phasehistory ] , we take @xmath167 and @xmath168 for all @xmath169 and all @xmath131 . at such an event
we define the new phases using the list below .
one can verify that all variables stay in their respective ranges if they are started in that range and evolving in this way .
* _ for each @xmath126 such that @xmath158 occurs _ + we set @xmath170 * _ for each @xmath126 such that @xmath154 occurs
_ + we set @xmath171 where @xmath172 is defined below , @xmath173 , @xmath174 , and @xmath175 is the number of ` input pulses ' that arrive at time @xmath176 .
* once all the @xmath154 have been checked , we say there is _ firing @xmath165 induced by pulse reception _ if @xmath154 occurs for some @xmath177 and results in @xmath178 . in this case
we set @xmath179 ,  @xmath180 and renumber @xmath181 for all @xmath182 , as if @xmath152 occurred .
one can verify that the order of ` processing ' simultaneous events with different @xmath126 does not affect the final outcome ; hence the map is well - defined .
the condition @xmath183 in ( [ eq : kinputsnow ] ) excludes self - interaction as in @xcite which we will consider exclusively in the following .
the response function @xmath184 expresses the response of an oscillator at phase @xmath185 to an input @xmath186 .
the model is specified once the delay @xmath187 , the coupling strength @xmath188 and the function @xmath34 are known . observe that given an initial @xmath189 and @xmath190 for @xmath191 we can define @xmath192 the phase history at time @xmath81 . this way we complete the definition of @xmath193 ( [ eq : semiflowforoscs ] ) for times @xmath123 $ ] .    for continuous functions @xmath194\rightarrow
{ \mathbb{r}}$ ] we define their _ difference increments _    @xmath195    we say @xmath34 is _ uniformly strictly monotonic increasing _ on @xmath60 $ ] if @xmath196 for all @xmath197 ( note that this implies that @xmath172 is strictly monotonic increasing in @xmath185 ) .
[ response functions with increasing response ] [ def : cv ] we define @xmath198 be the set of response functions @xmath199\times [ 0,c ] \rightarrow [ 0,\infty)$ ] for some @xmath200 with @xmath201 satisfying    * _ continuity : _
@xmath172 is jointly continuous in @xmath185 and @xmath202 .
* _ zero response to zero input : _
@xmath203 for all @xmath185 .
* _ uniform strict monotonicity : _ for any fixed @xmath174 the function @xmath204 is uniformly strictly monotonic increasing in @xmath185 on @xmath60 $ ] . *
_ upper bound linear in @xmath202 : _ for any fixed @xmath174 , monotonicity of @xmath34 means there is a unique @xmath205 that solves @xmath206 which gives the maximum response ; we assume this is bounded above by a linear function of @xmath202 .
more precisely we assume there is a @xmath207 such that @xmath208 for all @xmath209
.    condition ( iv ) allows one to obtain a bound on the total change in phase over one period .
note that @xmath210 for all @xmath211 and that for @xmath212 , the response @xmath172 is such that the updated phase ( [ eq : epshat ] ) is supra - threshold , resulting in an induced firing ; hence @xmath213 is the largest response to an input of magnitude @xmath202 .
we say that response functions @xmath214 give _ increasing response _ because the response increases    with @xmath185 . in particular , concave ( down ) potential functions give rise to this form of response ( see lemma  [ lemma : uv ] below ) .
we give this space the structure of a metric space with norm @xmath215\times[0,c ] }    @xmath216 .
closeness of @xmath34 and @xmath217 in this norm imply that both the functions and their differences increments are close , by the following lemma .
[ lemma : deltabound ] suppose that @xmath218 ; then @xmath219 .
* proof : *  suppose that @xmath218 , fix @xmath202 and drop the second argument of @xmath172 for convenience .
then @xmath220 for all @xmath185 .
the triangle inequality implies that @xmath221 for all @xmath222 and @xmath223 . for a given value of @xmath185
there is a sequence @xmath224 such that @xmath225 as @xmath226 . by considering @xmath227 in ( [ eq : dv2delta ] )
we can conclude that @xmath228 .
similarly by considering @xmath229 such that @xmath230 as @xmath226 we have using ( [ eq : dv2delta ] ) that @xmath231 .
this implies @xmath232 for all @xmath185 and @xmath202 .
@xmath41    the resetting rule @xmath233 with a twice continuously differentiable potential function @xmath1 was considered previously @xcite to describe a phase jump induced by a reception ( @xmath154 ) of a pulse of strength @xmath202 .
this ` mirollo - strogatz ' case is a special case of the more general reception rule we use here , if we substitute @xmath234 we now show that , if @xmath1 is increasing concave downwards , this fits into the more general reception scheme of definition  [ def : cv ] . in figure
[ fig : vanduplots ] we illustrate a typical choice of @xmath1 and its consequence for @xmath34 for given @xmath185 and @xmath202 .
[ lemma : uv ] consider a response defined by @xmath235 in ( [ eq : phi_uphi],[eq : vfromu ] ) where @xmath236 and @xmath1 is an unbounded twice continuously differentiable function @xmath237 that satisfies @xmath238 , @xmath239 , @xmath240 and @xmath241 for @xmath242 $ ] where @xmath243 ( @xmath244 ) are positive constants .
then @xmath245 for small enough @xmath246 .
note that the details of the definition of @xmath1 are only critical in the range @xmath60 $ ] .
this is because any time that @xmath247 we will get resetting ; the definition on @xmath248 is simply for expositional convenience .
* proof : *  one can directly verify that the conditions ( i)(ii ) in the definition of @xmath198 are satisfied . to show ( iii ) , note that @xmath1 monotonic increasing and surjective on its range means it has a unique inverse @xmath249 that is also monotonic increasing . to see that ( iii ) also holds , first observe that @xmath250 the condition @xmath251 means that @xmath1 is strictly monotonic increasing .
the condition @xmath252 means that @xmath253 for all @xmath15 satisfying @xmath254 . together with the monotonicity of @xmath249 we obtain , substituting @xmath255
, @xmath256 similarly , @xmath257 means that @xmath258 is monotonic decreasing ; @xmath259 implies @xmath260 for all @xmath185 and thus , using ( [ eq : ujumpbounded ] ) , @xmath261 substituting into ( [ eq : vprime ] ) this leads to @xmath262 hence for fixed @xmath174 we have @xmath263 for all
@xmath264 , and so ( iii ) holds . finally , to show ( iv ) , note that for @xmath265 eq .
( [ eq : phim ] ) , which determines @xmath205 , becomes @xmath266 so that by definition ( [ eq : phim ] ) , the maximum change at a reset is @xmath267 clearly , @xmath268 for all @xmath174 .
moreover , @xmath269 by definition ( eq .  ( [ eq : vonems ] ) ) , @xmath270 and @xmath271 ^ 3}<0\ ] ]    because @xmath251 and @xmath257 in the relevant range .
thus @xmath272 is increasing and concave and so we have the linear upper bound @xmath273 where @xmath274
. @xmath41     satisfying the conditions for lemma  [ lemma : uv ] that hence giving a resetting function @xmath172 with increasing response .
@xmath275 in this case is defined in ( [ eq : ub ] ) with @xmath276 .
[ fig : vanduplots ] , width=226 ]      we proceed to show that the system of pulse coupled oscillators with increasing response will reduce to a finite dimensional system for an open set of parameter values .
theorem  [ thmfinite ] gives an estimate @xmath147 for the largest number of ` past firings ' one needs to consider in the limit of large time .
we can effectively reduce to a system of dimension at most @xmath277 after a finite time .
the proof works by considering the maximum number of resets that can occur to an oscillator in any time interval @xmath278 $ ] .
this depends on the number of firings that occurred in the previous time interval of the same length , giving the recurrence inequality in lemma  [ lemestimate ] .
we reduce the dynamics from a hybrid delay differential equation with infinite dimensional phase space to a finite dimensional system by restricting to phase histories of the form : @xmath279 namely , those that have at most @xmath147 firings within the delay time @xmath116 .
[ lem : finite ] suppose that @xmath189 and @xmath280 are in @xmath281 and @xmath282 .
then @xmath283 for all @xmath7 .    *
proof : *  this result states that any two initial states in @xmath281 whose past @xmath147 firings and phases are the same will have precisely the same phases ( and hence firings ) at all times in the future .
this is because by definition all firings beyond the @xmath147th in the past will not affect the future phases .
@xmath41    the previous lemma does not necessarily imply that @xmath284 .
[ thmfinite ] suppose that @xmath285 . then there is a @xmath286 ( depending on @xmath202 , @xmath116 and @xmath34 ) such that for any @xmath287 there is a finite @xmath81 ( depending on @xmath189 , @xmath202 and @xmath34 ) such that @xmath288 for all @xmath289 .
moreover , after time @xmath32 the number of firings of an oscillator that can appear in any time interval of length @xmath116 is at most @xmath290 where @xmath291 is the smallest integer less than or equal to @xmath15 .    for any @xmath147 such that the conclusion of theorem  [ thmfinite ] holds we say the dynamics _ reduces _ to @xmath292 .
we prove theorem  [ thmfinite ] using the following lemma . for any initial condition @xmath132 let @xmath293 ( @xmath294 ) and suppose that the number of firings of @xmath295 in the interval @xmath296 is @xmath297 .
[ lemestimate ] for all @xmath298 we have @xmath299    * proof : *  suppose that there have been at most @xmath300 firing events of each of the @xmath115 oscillators during the time interval @xmath296 .
consider the events that affect the phase @xmath301 of an arbitrary oscillator @xmath302 in the future time interval @xmath303 ; there will be @xmath304 events in @xmath303 where at each event @xmath305 , @xmath306 pulses of the oscillators @xmath307 , @xmath177 are simultaneously received by @xmath126 .
the total number of received pulses is @xmath308 .
the total change @xmath309 in phase of oscillator @xmath126 up to time @xmath116 is hence at most @xmath310 hence oscillator @xmath126 will have fired at most @xmath311 times in the interval @xmath303 . since @xmath126 is arbitrary , we have in general at most @xmath312 firing events of each oscillator in the interval @xmath313 $ ] and the lemma follows .
@xmath41    * proof : *  * ( of theorem  [ thmfinite ] ) * since @xmath314 applying lemma  [ lemestimate ] inductively we have @xmath315 and so because @xmath316 , after some finite time @xmath317 ( which can be estimated from ( [ eq : timeestimate ] ) given @xmath318 ) we have @xmath319 the remainder of the theorem is a consequence of applying lemma  [ lem : finite ] .
@xmath41    the estimate of @xmath147 will typically be a very poor upper bound .
if the number of oscillators @xmath115 is small one can improve on these bounds by noting that the pulses are ` quantized ' into multiples of @xmath320 .
in particular , one can replace the upper bound ( [ eq : vone ] ) by requiring only that there is a @xmath207 such that @xmath321 for all @xmath322 . for larger @xmath115 this approaches ( [ eq : vone ] ) .    as a consequence of lemma  [ lemestimate ] we can now justify defining the dynamics as a semiflow on @xmath323 .
a system with @xmath214 that satisfies @xmath324 will remain free of zeno states for any @xmath174 and @xmath151 .
this is because , although the number of firings in an interval @xmath278 $ ] can increase if @xmath325 is large , the maximum rate is given in lemma  [ lemestimate ] and there can not be infinitely many firings in any time interval of finite length .
another consequence of the proof of theorem  [ thmfinite ] is that the semiflow on @xmath326 has a restriction to @xmath327 that is a well - defined semiflow .      in the case
that @xmath330 we can apply theorem  [ thmfinite ] to see that @xmath331 , and so @xmath329 .
in fact , this is not a very tight estimate , and in many cases where ( [ eq : k=1 ] ) does not hold the dynamics still reduces to @xmath329 . for a particular initial condition @xmath189 one can practically check the dimension it reduces to by finding the @xmath332 of the number of firings of one oscillator in any time interval of length @xmath116 for a neighbourhood of that initial condition .
if at most one previous firing is needed to determine the phases uniquely in the future , the dynamics reduces to a semiflow on a @xmath328 dimensional space @xmath333 where we write @xmath334 in shorthand . for the oscillators we consider , the semiflow @xmath335 for @xmath119 is well defined for @xmath202 and @xmath116 small enough ; appendix  [ app : k=1 ] gives an explicit form on this reduced space .
the main result in this section is a proof that unstable attractors exist and are robust in the system of delay pulse coupled oscillators discussed in section  [ secpulsecos ] .
for this section we consider the potential function @xmath336 used in @xcite .
computing the associated response function @xmath34 defined by ( [ eq : vfromu ] ) we find an affine relationship between @xmath185 and @xmath34 : @xmath337 and one can check that @xmath214 ; this is a particular case of lemma  [ lemma : uv ] ; the relationship between @xmath1 and @xmath34 is illustrated in figure  [ fig : vanduplots ] .
[ thm : uasexist ] consider the system of oscillators defined in section [ sect : defsemiflow ] with @xmath1 defined by ( [ eq : ub ] ) with @xmath338 , @xmath339 and @xmath340 and @xmath276 .
this system has a periodic orbit that is an unstable attractor with zero measure local basin .
* proof : *  by applying theorem  [ thmfinite ] one can directly verify that any initial condition @xmath341 will reduce to the @xmath342 dimensional space @xmath343 after a finite time , i.e. such that only the phase and last four firings of each oscillator contribute to the dynamics at any point in the future .
now consider the open set of initial conditions defined by @xmath344 where @xmath345^{8}\subset \pi_{1 } \left({\mathcal{p}}^4\right).\ ] ] as a slight abuse of notation , we will also write @xmath346 as meaning @xmath347 if the context is clear .
note that all the @xmath348 , meaning that @xmath349 .
theorem  [ thmfinite ] means that at points in the future we may in principle arrive at states that are in @xmath350 for @xmath351 but not for @xmath352 .
however , we will verify below that @xmath353 for all @xmath12 .
note that @xmath354 is an open set , as is its image under @xmath355 is open for any @xmath356 ; we will show that it is in the basin of attraction of a saddle periodic orbit . in what follows
we write @xmath357 to denote the @xmath147th intersection of the trajectories in @xmath358 with one or more of the ` events ' @xmath152 or @xmath154 .
since the information about the original initial condition of the @xmath359 has been lost we list only the first four components .
we compute that @xmath360\times[0.749,0.751]\times[0.349,0.351]\times\{0\ } \label{eq : c1}\ ] ] due to the first event being @xmath361 .
the next event is @xmath362 which will give @xmath363~~\}\label{eq : c2}\ ] ] and so we have reduced to a one dimensional subspace of phase space at this point .
a very tedious direct numerical verification listed in table  [ tab : c0traj ] shows that _ all _ of the points starting in @xmath364 after a further 32 events reach a single point @xmath365 where @xmath366 and @xmath367 can be computed exactly ( see appendix  b ) .
we present in figure  [ fig_return ] the values of the coordinates depending on choice of @xmath368 over a larger range than in @xmath364 ; observe that a one dimensional set of possible initial values converges to exactly the same point .
the point @xmath369 ( [ eq : c22p ] ) is a cross section to a periodic orbit where two pairs of oscillators remain synchronized .
since the periodic orbit through @xmath369 contains the open set @xmath346 in its basin it is a milnor attractor for any reduction of the system to finite dimensions , for example to @xmath370 .
@xmath371     dependence of dynamics on the initial phases @xmath372 .
panel ( a ) shows the phases ( @xmath373(@xmath374 ) , @xmath375(@xmath376 ) , @xmath377 ( @xmath41 ) , @xmath378(@xmath379 ) ) after @xmath380 further intersections with discontinuities of the semiflow ( ` events ' ) in dependence of @xmath368 .
observe that there is a large range of @xmath368 for which the dynamics arrives at precisely the same point .
panel ( b ) shows the actual attractor reached from the above initial conditions ( shown after the @xmath381 reset of oscillator @xmath382 .
surprisingly , the same attractor is reached by several initial conditions whose trajectory after 34 events differs substantially .
, width=340 ]    now we show that the orbit through @xmath369 is unstable . in order to do this
we show that almost all perturbations of the first two coordinates grow to finite size under the action of the local return map .
a complete stability analysis is presented in the appendix [ app : unstable ] .
in particular , as is shown below , this proves instability in general , not only linear instability .    in this analysis
we consider the return map from the plane @xmath383 to itself .
the pulses @xmath359 are kept locked to the phases @xmath301 such that the nonlinear return map ( [ eq : f(delta ) ] ) @xmath384 is three - dimensional and maps a perturbation @xmath385 to the periodic orbit started at ( [ eq : c22p ] ) .
@xmath386    to the perturbation @xmath387 when the plane @xmath388 is reached the next time .
if we take a general perturbation ( satisfying @xmath389 ) just before arriving at @xmath390 , observe that the @xmath391 are all larger than @xmath116 and hence will not affect the future of the phases .
starting with such a perturbation , when it hits the next event we can arrive at an arbitrary perturbation satisfying @xmath392 such that immediately after the perturbation , events occur in the order @xmath393 . note that the order of the events @xmath394 and @xmath361 does not play a role as they become synchronized at their first firing .    as demonstrated explicitly in appendix b
, there is a small @xmath395 such that the return map ( [ eq : fdelta ] ) defined in the three - dimensional region @xmath396 is given by a smooth map with @xmath397 ( and such that @xmath398 is a line that intersects @xmath399 ) . in particular , given the map @xmath400 in terms of an arbitrary function @xmath214 ( see appendix b ) we analyse of the fully nonlinear local dynamics of a sufficiently small perturbation @xmath401 , @xmath402 , to the periodic orbit . for the following estimate ,
we write @xmath403 for @xmath404 in shorthand notation .
first note that requirement ( iii ) in def .
[ def : cv ] implies that there is a @xmath405 such that for all @xmath406 , @xmath407 for all @xmath174 , in particular , @xmath408 is monotonic increasing in @xmath185 .
this implies that    @xmath409    hence we find that @xmath410 implies that any perturbation @xmath411 satisfying @xmath412 increases ( e.g. in maximum norm ) . because the function @xmath413 is non - decreasing in @xmath414 , we see that there is a @xmath395 for which any sufficiently small perturbations ( satisfying @xmath415 ) leave the region @xmath399 in finite time @xmath416 , bounded by @xmath417    by the permutation symmetry of the system , the analysis in the case @xmath418 is analogous .
this means that almost all sufficiently small perturbations @xmath419 leave a neighbourhood of the attractor in finite time , proving that the attracting periodic orbit considered is also unstable .
the case @xmath420 can be treated analogously by the symmetry of the globally coupled network .
taken together there is a neighbourhood @xmath421 of @xmath422 such that all points starting in @xmath421 with @xmath423 will iterate away until they leave this neighbourhood .
this proves instability of the orbit through @xmath369 ( [ eq : c22p ] ) .
@xmath41    figure  [ fig : xc ] shows the dynamics starting from one initial condition ( [ eq : c0 ] ) used in the proof of the theorem . the convergence in finite time to the attractor
a small random perturbation to the phases of all oscillators is expanded by the instability , yielding convergence towards another attractor .    ,
@xmath340 , @xmath339 , @xmath424 , @xmath276 ) exhibiting an unstable attractor . at time zero , initial phases in @xmath346 are shown .
the time axis is discrete and shows the phase when oscillator @xmath425 fires . just after the firing event @xmath426 of oscillator @xmath425 at times @xmath427 ,
all phases @xmath428 are plotted ( @xmath373(@xmath374 ) , @xmath375(@xmath376 ) , @xmath377 ( @xmath41 ) , @xmath378(@xmath379 ) ) .
there is convergence of the trajectory in state space in finite time towards the periodic orbit passing through ( [ eq : c22p ] ) .
an instantaneous small random perturbation ( @xmath429 ) to the phases at time @xmath430 leads to a switching to another attractor , demonstrating instability of the original periodic orbit .
[ fig : xc],width=264 ]    [ thm : uasrobust ] the unstable attractors of theorem  [ thm : uasexist ] are robust in the sense that there is an open neighbourhood in @xmath198 of @xmath431 given by ( [ eq : vfromu],[eq : ub ] ) and an open neighbourhood of @xmath432 such that a system defined by any @xmath217 , @xmath433 , and @xmath434 in these neighbourhoods , has an unstable attracting periodic orbit that remains close to that of theorem  [ thm : uasexist ] .
* proof : *  we refer to the return map @xmath435 ( [ eq : f(delta ) ] ) derived in the appendix  b and already discussed in the proof of theorem  [ thm : uasexist ] .
note that in the derivation of the return map @xmath435 , the form of the transfer function @xmath436 and the precise values of the parameters @xmath116 and @xmath202 are not important , because @xmath435 continuously depends on @xmath436 , @xmath116 , and @xmath202 and changes its form only if the sequence of events occurring along the trajectory on or near the periodic orbit is altered .
thus , an equivalent periodic orbit exists if the dynamical system is slightly perturbed , i.e. if @xmath436 , @xmath116 , and @xmath202 , are slightly changed from their original values . in particular , one can choose @xmath437 , where @xmath34 is sufficiently close to the original .
since the mechanism of dimensional reduction and attractivity depends only on the supra - threshold synchronization mechanism given by the induced firing events @xmath165 , the slightly perturbed periodic orbit will also be an attractor .
moreover , for a sufficiently small perturbation to @xmath34 , @xmath116 , and @xmath202 we still obtain @xmath438 which implies instability .
the proof of the previous theorem ) . again
, since the return map @xmath435 does not change form if the parameters @xmath436 , @xmath116 , and @xmath202 are slightly perturbed , also the proof of instability is independent of these parameters of the system .
this completes the proof .
@xmath41    this result implies that qualitatively the same dynamics is obtained by any flow defined by a @xmath217 nearby the original @xmath34 ( [ eq : vfromu ] ) ; in particular it does not have to be defined in terms of a any potential function @xmath1 by ( [ eq : vfromu ] ) .
since the proof of theorem  [ thm : uasexist ] is rather involved , we discuss some intuition in terms of what we term a dimension jump of the semiflow . given a semiflow @xmath439 that is not invertible , one can nevertheless define , for all @xmath440 , a set - valued inverse for any @xmath12 by @xmath441 and the semiflow can be extended uniquely to an invertible flow precisely when @xmath442 is a single point for all @xmath15 and @xmath12 .
given sufficient regularity on @xmath443 , for example sufficient to ensure that @xmath442 is a union of manifolds , we can compute its dimension .
if all @xmath442 are non - empty then @xmath444 will be a non - decreasing function of @xmath131 .
the _ dimension jump _ @xmath445 at @xmath15 is of particular interest and we refer to @xmath446 as the _ dimension jump set _ , shown schematically in figure  [ fig_jumpset ] .
such behaviour can be found in systems exhibiting sliding modes @xcite , and we suggest that such systems may display unstable attractors .     that has a unique inverse @xmath447 for @xmath448 .
if we define @xmath449 then @xmath31 is in @xmath450 the jump set ( in this case @xmath451 ) and for any @xmath452 the dimension of the set @xmath453 is greater than zero .
[ fig_jumpset],width=377 ]    observe that trajectories converging to an invariant set that pass through @xmath450 may cause locally unstable invariant sets to be unstable attractors .
for the system of delay coupled oscillators with @xmath329 there are two events on which dimension jump can occur :    1 .
for one of the firing events @xmath158 or @xmath165 we lose dependence on the delayed phase @xmath454 .
2 .   at a resetting @xmath455 that induces firing in two oscillators @xmath456 and @xmath457
we synchronize the oscillators @xmath458 and @xmath459 beyond that point ; this also causes a dimension jump .    in out system ,
these events have ` jump sets ' that are codimension one and typically they cause a loss of one dimension ; however when the oscillators become synchronized , several such events can occur at once resulting in a loss of more dimensions .
the proof of theorem  [ thm : uasexist ] works by choosing an open set of initial conditions that reduces dimension down to one dimension after two further discrete events .
the final loss of dimension that reduces the open set down to one point in the section occurs only a while later , and this is because there is a global connection in phase space that we need to follow .
in this paper we give a rigorous definition for a novel type of milnor attractor , the _ unstable attractor _ , introduced in @xcite .
we classify these attractors in terms of presence or absence of local basin .
these attractors are surprising in many ways ; notably they have ` local repelling ' properties in that they have neighbourhoods such that almost all points leave the neighbourhood after a finite time .
we show that unstable attractors not only exist , but even appear robustly in certain classes of hybrid dynamical systems of interest as physical and biological models ; namely in networks of pulse - coupled oscillators .
the robustness that we have demonstrated is limited in the sense that we only consider perturbations that keep the time delay identical for all connections and the coupling function identical for all oscillators .
we do not consider perturbations that break the identical nature of the oscillators or the symmetry of the system .
on introduction of such general perturbations we expect no robustness for such unstable attractors , although one may obtain attractors that remain ` nearby ' to where a network of unstable attractors used to exist , i.e. close to a collection of unstable periodic orbits and their heteroclinic connections .    for perturbations that break only part of the symmetry
, unstable attractors may still persist .
for instance there is numerical evidence that unstable attractors can still exist even the network is not coupled all - to - all , cf .
the exact consequences of non - identical oscillators and more general broken symmetry remains to be studied .
if we perturb the system in such a way that it becomes smooth and invertible , the unstable attractors will disappear ; any saddle periodic orbit will have a zero measure basin simply by existence of a global stable manifold of dimension less than that of phase space ; however as indicated in section  [ secunstableatts ] one should be able to make the system arbitrarily differentiable and still obtain unstable attractors , as long as non - invertibility is retained .
in proposition  [ prop : nounsts ] we have shown that certain conditions on a general dynamical system exclude the possibility of unstable attractors with zero local basin . to generalize this finding to also exclude unstable attractors with positive local basin , we believe it will be necessary to strengthen these conditions
. generalization of our result is not straightforward because a ( local ) basin may be highly irregularly , as for riddled basin attractors .
our proofs of the main results are not very elegant ; this is because they rely on global features of the hybrid dynamical system , and we have not found them easily reducible to simpler arguments .
we believe however that it should be possible for exact numerical methods @xcite to identify when unstable attractors are present in such systems .
the numerical results in @xcite suggest that networks of unstable attractors exist in a wide range of systems of pulse - coupled oscillators .
preliminary studies have shown that the networks of unstable attractor may often be transitive for a wide range of parameters and system sizes .
revealing the mechanisms underlying the existence , structure and prevalence of such networks of unstable attractors is a challenging task for future research .
* acknowledgements : * pa thanks the max - planck - institut fr dynamik und selbstorganisation ( formerly : mpi fr strmungsforschung ) for their hospitality and financial support , as well as a grant from the leverhulme foundation for partial support .
mt thanks leonid bunimovich , stefan grokinsky and raoul - martin memmesheimer for useful and interesting conversations relating to this work as well as the department of mathematical science , university of exeter for their hospitality and financial support .
we give the semiflow of globally coupled oscillator networks explicitly for the case @xmath329 by the following algorithm .
consider @xmath460 where the @xmath126th component represent the @xmath126th phases and the @xmath461th component represents the time since the @xmath126th oscillator fired .
suppose that @xmath462 for all @xmath7 .
we note that we can represent the flow on @xmath463 by @xmath464 for all @xmath465 such that @xmath466 and @xmath467 .
after time @xmath468 we hit the first event and apply a mapping according to the following algorithm :    * set @xmath469 .
* for each @xmath126 such that @xmath470 we say @xmath455 occurs for all @xmath471 and set @xmath472 for @xmath186 defined as in ( [ eq : epshat ] ) .
* if there is an @xmath126 such that @xmath178 then we say @xmath158 occurs and set @xmath179 and @xmath473 .    after applying this algorithm we are assured that @xmath474 , as @xmath475 and @xmath476 for @xmath477 . thus we can take this state as a new initial condition and rerun the algorithm sequentially from event to event .
for the main results of this paper ( theorems 2 and 3 ) we use a nonlinear return map @xmath400 of perturbations @xmath419 to the unstable attracting periodic orbits in a network of @xmath338 oscillators using the semiflow detailed in appendix  [ app : k=1 ] .
we detail the derivation of this return map and analyze the periodic orbit that is the unstable attractor .
we choose @xmath478 where @xmath479 , and consider cases where the initial @xmath480 .
first we analyse the dynamics on the periodic orbit . to establish instability we then consider , without loss of generality , perturbations to the periodic orbit phases at a point on the orbit where the time since last pulse is greater than @xmath116 for all oscillators .
this means that the future depends only on the current phases and we can effectively consider @xmath481 for all @xmath482 and all @xmath131 satisfying @xmath483 where @xmath484 defines the time of the last event such that @xmath485 and @xmath486 the time of the next event analogous to ( [ eq : algort ] ) . within this convention
, we consider the return map from the three - dimensional plane defined by @xmath388 to itself and establish instability for this return map to prove instability of the periodic orbit .
the following analysis applies for the parameters @xmath340 , @xmath339 , and the function @xmath424 ( [ eq : ub ] ) where @xmath276 .
it will turn out that it also applies for a neighbourhood thereof , showing robustness of the instability of the periodic orbit . for compact notation
, we introduce the _ transfer function _ @xmath487 this transfer function mediates the action of an incoming pulse of strength @xmath202 on the phase @xmath185 .
we obtained the sequence of discrete events and the numerical values in the following analysis for the specific choice ( [ eq : vfromu],[eq : ub ] ) .
first we determine the dynamics of the unperturbed orbit ( table  [ tab : unperturbed_po ] ) . throughout the table
, the phases after each event are labelled @xmath488 where the first index labels the oscillator @xmath482 and the second index counts the events @xmath147 occurring for a _ perturbed _ trajectory ( see table [ tab : perturbed_po ] below ) starting with @xmath489 for the initial condition ( first row of table [ tab : unperturbed_po ] ) .
since the second oscillators @xmath425 and @xmath490 in the two synchronous groups have the same phases as the oscillators @xmath382 and @xmath491 , respectively , we define @xmath492 for all @xmath147 .
we show that the unperturbed orbit indeed is periodic with period @xmath493 given by the time @xmath131 in the last row of table [ tab : unperturbed_po ] where@xmath494 is also read off the table . here
we use the abbreviation @xmath495 for the strength of an individual pulse .    .
events occurring during one period of one of the unstable attracting periodic orbits in the system of four coupled oscillators ( see text for explanation ) .
[ tab : unperturbed_po ] [ cols="<,<,^,^",options="header " , ]     table [ tab : perturbed_po ] displays the time evolution of a general perturbation event by event . as for the unperturbed orbit described above
, the left column gives the sequence of events and the second left column gives the time of these events .
the right four columns give the phases of oscillators @xmath496 right after an event has occurred .
the first row gives the initial condition at @xmath176 , just after the signal has been sent , @xmath426 , by oscillator @xmath425 .
before , at time @xmath497 , the signal of oscillator @xmath382 has been sent , denoted @xmath498 ; this signal will travel until time @xmath499 . the initial phases ( [ eq : phi(0 ) ] ) are displayed in the first column . here
@xmath500 by definition .
throughout the table , the phases after each event are labelled @xmath501 where the first index labels the oscillator @xmath126 and the second index counts the events @xmath147 starting with @xmath489 for the initial condition as was also used in table  [ tab : unperturbed_po ] above .
for instance , the first event ( after the initial condition ) , @xmath329 , at time @xmath499 is the reception @xmath502 of the signal sent by oscillator @xmath382 .
since there are no self - interactions , the phase of oscillator @xmath382 is only shifted in time , @xmath503 the phases of all other oscillators @xmath504 additionally jump because of the incoming subthreshold signal , @xmath505    again , as for the unperturbed dynamics , going through the table row by row , the phases are calculated event by event .
the phases @xmath501 for @xmath496 right after the events are approximated by the unperturbed phases @xmath488 defined in table [ tab : unperturbed_po ] , @xmath506 for all @xmath147 for which the @xmath488 are defined , @xmath507 . here
@xmath508 denotes a function of all @xmath509 , @xmath510 , that approaches zero at least linearly if all @xmath511 , i.e.@xmath512 .
we thus obtain the return map@xmath513 from the last row of table [ tab : perturbed_po ] . here
the @xmath514 again depend on the original perturbation vector @xmath419 . from this identification
we obtain the components @xmath515 note that the component functions @xmath516 in ( [ eq : a : f_i ] ) are of order @xmath517 ; moreover the dynamics clearly depends only on the initial perturbation @xmath518 to this orbit .
it is important to note that although @xmath519 and @xmath520 in general , this group is resynchronized at the next return , @xmath521 .
10 p. antsaklis , w. kohn , a. nerode and s. sastry ( eds ) _ hybrid systems .
papers from the @xmath522 workshop held at cornell university , ithaca , new york ( 1994 ) .
lecture notes in computer science , 999 .
springer , berlin , ( 1995 ) .
p. ashwin , riddled basins and coupled dynamical systems . to appear in `` dynamics of coupled map lattices and of related spatially extended systems '' ,
chazottes and b. fernandez , springer lecture notes in physics ( 2005 ) .