homo and heteroclinic bifurcations constitute the core of our understanding of complicated recurrent behaviour in dynamical systems .
the history goes back to poincar in the late @xmath0 century , with major subsequent contributions by the schools of andronov , shilnikov , smale and palis .
these schools have been founded on a combination of analytical and geometrical tools which developed a quite good understanding of the qualitative behaviour of those dynamics .
differential equations modelling physical experiments frequently have parameters which appear in the differential equations , and it is known that qualitative changes may occur in the solution structure of these systems as the parameters vary .
the dynamical behaviour of systems can be strongly influenced by special geometric or analytic invariant structures appearing in the equations ( _ e.g. _ divergence - free , preserving or reversing symmetries ) .
one may ask which dynamical behaviour we would expect to see in the presence of a given invariant structure .
generically , this is a hard question to answer , but some questions could be partially answered , by considering local and global bifurcations of low codimension .      a bykov cycle on a three - dimensional manifold is a heteroclinic cycle between two hyperbolic saddle - foci of different morse index , where one of the connections is transverse and the other is structurally unstable .
there are two different possibilities for the geometry of the flow around the one - dimensional connection depending on the direction solutions turn around it .
all literature about the michelson system @xcite considered that in the neighbourhoods of the two saddle - foci , trajectories wind in the same direction .
an immediate question arises :    * what happens if trajectories wind with opposite directions near each node ?
the first author and duarte proved in  @xcite that the set of @xmath1-divergence - free vector fields defined in a compact three - dimensional riemannian manifold without boundary , has a @xmath1-residual set such that any vector field inside it is anosov or else , the flow associated to it has dense elliptic solutions in the phase space .
furthermore , in @xcite , also in this context , the authors proved that if the vector field is not anosov , then it can be @xmath1-approximated by another divergence - free vector field exhibiting homoclinic tangencies .
however , a conservative vector field whose flow has a persistent bykov cycle may lie outside these residual / dense subsets and the developed theory can not be applied for this degenerated class of systems . since bykov cycles can not be anosov ( due to the presence of equilibria ) , the dichotomies in @xcite suggest the following problem :    * could we perform a @xmath1-perturbation within the set of vector fields whose flow has a bykov cycle in such a way that the elliptic periodic solutions are dense and/or tangencies occur ?    in this paper , we partially answer the questions * ( q1)**(q2)*. we present a mechanism for the coexistence of hyperbolic and non - hyperbolic dynamics arising in a neighbourhood of a conservative bykov cycle where trajectories turn in opposite directions near the two saddle - foci - see figure [ conservativo2 ] . we show that in a @xmath1-open class of divergence - free vector fields , tangencies of the invariant manifolds of two hyperbolic saddle - foci densely occur .
we prove that the global dynamics is persistently dominated by heteroclinic tangencies and by the existence of infinitely many elliptic points coexisting with hyperbolic dynamics arising from @xmath1transversality .
tangencies are strongly connected with the cocoon bifurcations .
we describe an extended version of these bifurcations which can be detected numerically with the _ time delay function
_ introduced by lau @xcite .
see also @xcite .
we also find the relations between these four dynamical phenomena : chirality , tangencies , generalized cocoon bifurcations and elliptic solutions . in the conservative setting ,
cocoon bifurcations are connected with unfoldings of a hopf - zero singularity , which has been shown to occur in climatological models with seasonal forcing
.     represents the set of @xmath1-divergence free vector fields on @xmath2 . on the dashed line
, we find the vector fields whose flow has a given bykov cycle and its continuation .
the set @xmath3 is the open set ( with the induced topology ) in which theorem [ main1 ] is valid .
, height=113 ]      we recall some important results about heteroclinic tangencies for general and conservative systems .
the systematic study of bifurcations of tangencies was started by gavrilov and shilnikov @xcite in the seventies , for the case of two - dimensional dissipative diffeomorphisms .
the authors established that diffeomorphisms with homoclinic tangencies might separate systems with regular and chaotic dynamics , _
i.e. _ they belong to the boundary of morse - smale systems and the transition through this boundary corresponds to an @xmath4explosion .    the theorem on the cascade of periodic sinks / sources near a tangency plays an important role
see newhouse @xcite : in the dissipative case , for any one parameter family that unfolds generically a quadratic homoclinic tangency , there exists a sequence of intervals of values of the parameter such that the corresponding diffeomorphism has sinks / sources .
one of the fundamental results in homoclinic bifurcations was established in newhouse @xcite : the existence of regions of the space of two - dimensional diffeomorphisms where tangencies are dense .
these regions are called _ newhouse regions _ ; the result has been generalized for the multidimensional case by gonchenko @xcite .
dynamics of systems within these regions is complex and rich and , as claimed in @xcite , it is impossible to give the complete description of bifurcations of such systems .
the majority of these results was obtained for the broader case of general systems without any restriction on the preservation of some invariant structure .
they can not be directly applied to volume - preserving and/or reversible systems because they require special considerations .
nevertheless , the main geometric and analytical arguments can be also used for systems with additional structures , having special care . in this direction , important results on the birth of elliptic periodic points in area - preserving maps under bifurcations of tangencies were obtained in @xcite .
see also the paper by lamb and stenkin @xcite about _ mixed dynamics .
_    bifurcations of single - round periodic solutions were studied in @xcite for two - dimensional symplectic maps close to a map having a quadratic homoclinic tangency .
the corresponding first return maps have been derived , bifurcations of their fixed points have been studied and bifurcation diagrams for one - parameter general unfoldings have been constructed . using renormalization results , the existence of one - parameter cascades of elliptic single - round periodic solutions has been proved .
let @xmath2 be a three - dimensional closed and connected @xmath5 riemaniann manifold without boundary , endowed with a volume - form and let @xmath6 denote the lebesgue measure associated to it .
let @xmath7 be a @xmath8 vector field , with @xmath9 . by picard s theorem on the existence and uniqueness of solutions of differential equations
, @xmath10 integrates into a complete @xmath8 flow @xmath11 in a sense that @xmath12 we say that @xmath13 is _ divergence - free _ if its divergence is equal to zero _ i.e. _ if @xmath14 computed in local coordinates @xmath15 , see @xcite .
equivalently , @xmath10 is divergence - free if , accordingly to the liouville formula , the volume - measure @xmath6 is invariant for the associated flow . in this case
we say that the flow is _ conservative _ , _ volume - preserving _ or _ incompressible _ and is such that @xmath16 , for all @xmath17 and for all @xmath18 .
we denote by @xmath19 , @xmath9 , the space of @xmath8 divergence - free vector fields on @xmath2 and we endow this set with the usual @xmath8 whitney topology .
let also denote by @xmath20 the space of @xmath8 general vector fields on @xmath2 without no divergence - free constraint .
we say that an equilibrium @xmath21 of @xmath10 ( _ i.e. _ @xmath22 ) , is _ hyperbolic _ if @xmath23 has neither pure complex nor zero as an eigenvalue . a hyperbolic saddle - equilibrium with a pair of conjugated non - real eigenvalues
is called a _
saddle - focus_. for regular points , ( _ i.e. _ non - equilibrium points ) we define hyperbolicity with respect to the poincar map in a standard way
see e.g.  @xcite . given a hyperbolic equilibrium point @xmath24 with respect to a vector field @xmath10 , we denote by @xmath25 its unstable / stable global manifold .
local unstable / stable manifolds are denoted by @xmath26 . when there is no ambiguity we omit the letter @xmath10 in these notations .
open class of divergence - free vector fields defined on a three - dimensional compact manifold , tangencies of the invariant manifolds of two hyperbolic saddle - foci occur densely.,height=264 ]      in the present paper , all equilibria are hyperbolic saddle - foci .
the dimension of the unstable manifold of a saddle - focus will be called the _
morse index _ of the saddle . given
two equilibria @xmath27 and @xmath28 , a _ heteroclinic trajectory _ from @xmath27 to @xmath28 , denoted @xmath29 $ ] , is a solution of @xmath30 contained in @xmath31
. there may be more than one solution from @xmath27 to @xmath28 .
a heteroclinic cycle is a finite collection of equilibria together with a set of heteroclinic trajectories connecting the equilibria in a cyclic way .
more details in field @xcite .    in a three - dimensional smooth manifold , a _ bykov cycle _ is a heteroclinic cycle associated to two hyperbolic saddle - foci with different morse indices , in which the one - dimensional manifolds coincide and the two - dimensional invariant manifolds have a transverse intersection @xcite .
the terminology _ bykov cycle _ is because it was bykov @xcite who studied dynamical properties from the existence of similar cycles in a more general context than conservative ones .
this type of cycle , or rather its bifurcation point , is also called by _
t - point _ by glendinning and sparrow @xcite .
in this section , we state the main theorems and we relate them with other results in the literature .      our object of study is the dynamics around an autonomous divergence - free smooth vector field whose flow has a bykov cycle , for which we give an accurate description here .
specifically , we study a @xmath1vector field:@xmath32 whose flow has the following properties ( see figure [ conservativo2 ] ) :    1 .
[ h1 ] there are two hyperbolic saddle - foci equilibria @xmath27 and @xmath28 . assuming that , for @xmath33 , the morse index of @xmath34 is @xmath35 , the eigenvalues of @xmath23 are : 1 .
@xmath36 and @xmath37 , where @xmath38 and @xmath39 , for @xmath40 ; 2 .
@xmath41 and @xmath42 , where @xmath43 and @xmath44 , for @xmath45 .
[ h2 ] there is a heteroclinic cycle @xmath46 consisting of @xmath27 , @xmath28 and two one - dimensional heteroclinic connections : @xmath47 $ ] and @xmath48 $ ] .
the solution @xmath47 $ ] is called the _
fragile connection_. 3 .
[ h3 ] at the heteroclinic connection @xmath48 $ ] , the two - dimensional manifolds @xmath49 and @xmath50 meet transversely .
property * ( h3 ) * is generic .
indeed , the transverse intersection of @xmath49 and @xmath50 of * ( h3 ) * persists under @xmath1perturbations , whereas the fragile connection @xmath47 $ ] does not . although bykov cycles appear naturally in systems with some symmetry as in knobloch _
et al _ @xcite , michelson @xcite , rodrigues and labouriau @xcite , unless it is explicitly said , in this paper we state general results without any kind of restriction about the symmetry .
there are two possibilities for the geometry of the solutions around the cycle , depending on the direction they turn around the one - dimensional heteroclinic connection @xmath29 $ ] , in the neighbourhoods of the saddle - foci . under the context of hypotheses
* ( h1)**(h3 ) * , we introduce the following definition adapted from @xcite .
let @xmath51 and @xmath52 be small disjoint neighbourhoods of two saddle - foci @xmath27 and @xmath28 with disjoint boundaries @xmath53 and @xmath54 , respectively .
typical trajectories starting at @xmath55 near @xmath50 go into the interior of @xmath51 in positive time , then follow the connection @xmath56 $ ] , go inside @xmath52 , and then come out at @xmath54 .
let @xmath57 be a piece of the described trajectory from @xmath53 to @xmath54 .
now , join its starting point to its end point by a line segment , forming a closed curve .
this curve is called the _ loop of _ @xmath58 .
[ chi ] the two saddle - foci @xmath27 and @xmath28 in @xmath46 have the _ same chirality _ if the loop of every trajectory in @xmath53 is linked to @xmath46 in the sense that , for every @xmath58 close to @xmath59 , the loop of @xmath58 and @xmath46 can not be disconnected by an isotopy .
otherwise , we say that @xmath27 and @xmath28 have _ different chirality_.    in contrast with the findings of @xcite , if the two nodes have different chirality , then the rotations may cancel out .
this is a key idea for the proof of the main results , which will be formalized via hypothesis * ( h4)*.    1 .
[ h4 ] the saddle - foci @xmath27 and @xmath28 have different chirality .
a good overview of the previous hypothesis has been considered in @xcite .
it corresponds to the _ non - concatenation property _ stated in @xcite .
property * ( h4 ) * means that there are two open neighbourhoods of the equilibria such that for any trajectory going from the first to the second , the direction of its turning around the heteroclinic connection @xmath47 $ ] is different .
a direct corollary is the following :    the condition * ( h4 ) * is persistent under _ isotopies _ : if it holds for @xmath10 , then it holds for continuous one - parameter families containing it , as long as the fragile connection remains .    in order to prove our main results
, we need some additional assumptions on the set of parameters @xmath60 that determine the linear part of the vector field @xmath10 at the equilibria . for any @xmath61 , let @xmath62 be the subset of parameters given by : @xmath63    the constant @xmath64 is related with the construction of the global map around the cycle and will be clarified later .
observe that condition ( [ final_formulae ] ) in the definition of @xmath62 is satisfied by an open set of parameters @xmath65 , @xmath66 , @xmath67 , @xmath68 and does not involve the quantities @xmath37 and @xmath69 .
the set @xmath70 has full lebesgue measure within the set of parameters @xmath62 .
hereafter , let @xmath71 denote the set of vector fields on @xmath2 satisfying * * ( h1)***(h5 ) * , which is a baire space @xcite .
as depicted in figure [ scheme1 ] , the set @xmath72 may be seen as an open set in the space @xmath73 , endowed with the induced topology .
let @xmath74 .
the vector field @xmath75 in dumortier _
et al _ @xcite associated to the michelson system satisfies * ( h1)**(h3 ) * but condition * ( h4 ) * does not hold .
@xmath76-reversibility where @xmath77 is a natural obstacle to * ( h4)*. in the next section we show a way to construct theoretically a divergence - free vector field displaying a bykov cycle and satisfying * * ( h1)****(h5)**. this is what we call a _ conservative bykov cycle_.      [ bykovexists ] any closed , connected , riemannian three - dimensional manifold @xmath2 supports a bykov heteroclinic cycle satisfying properties * ( h1)**(h5 ) * associated to some divergence - free vector field @xmath78 , @xmath9 .
let us consider a riemannian closed , connected and three - dimensional manifold @xmath2 .
we let @xmath79 be two hyperbolic saddle - foci associated to linear divergence - free vector fields @xmath80 and @xmath81 and complying the hypotheses * * ( h1 ) * * and * ( h5 ) * , and moreover in such a way that * ( h4 ) * also holds .
now , we would like to plunge in @xmath2 the open sets containing @xmath27 and @xmath28 .    by the poincar - hopf theorem  @xcite ,
we know that the sum of the indices of the equilibria of a vector field equals the euler characteristic of @xmath2 .
so , we carefully choose a finite number @xmath82 of traceless equilibria , _
i.e. _ equilibria associated to divergence - free linear vector fields @xmath83 and fulfilling the poincar - hopf theorem .    using moser s charts @xcite we send ( locally ) the vector fields @xmath84 and
then , using the pasting lemma @xcite we extend to a divergence - free vector field in the whole manifold .
let @xmath85 be the obtained vector field and we still denote by @xmath86 be two hyperbolic saddle - foci with the stable manifold of @xmath34 having dimension @xmath35 , @xmath87 .
using @xcite we consider @xmath88 @xmath1-close to @xmath10 with a dense trajectory with initial condition @xmath89 , @xmath90 .
since @xmath90 is dense , it passes arbitrarily close to @xmath91 and to @xmath92 . by the @xmath1-connecting lemma
@xcite we consider @xmath93 @xmath1-close to @xmath94 such that @xmath95 intersects @xmath96 .
since these manifolds are two - dimensional and @xmath2 is three - dimensional we can assume that the intersection is transversal .
using again @xcite , we consider @xmath97 @xmath1-close to @xmath98 with some dense trajectory for the flow @xmath99 and with @xmath100 still intersecting @xmath101 transversely .
since the solution is dense , it passes arbitrarily close to @xmath102 and to @xmath103 .
again , by the @xmath1-connecting lemma , we consider @xmath104 @xmath1-close to @xmath105 such that @xmath106 intersects @xmath107 , non - transversely of course , and with @xmath108 still intersecting @xmath109 transversely .
this satisfies the hypotheses * ( h2 ) * and * ( h3 ) * and therefore we obtain a bykov heteroclinic cycle associated to the divergence - free vector field @xmath110 .    by technical reasons , in section  [ sectionlocaldynamics ] , our main results ask for at least @xmath111-regularity of the vector field @xcite .
nevertheless , if we start with an element in @xmath112 exhibiting a bykov cycle and since the hyperbolic equilibria are @xmath1-stable , similar arguments to those used in proposition  [ bykovexists ] can ensure that @xmath1-arbitrarily close to the initial vector field , there exists an element in @xmath113 satisfying * ( h1)**(h5)*. this is due to the fact that @xmath114 is @xmath1-dense in @xmath112 by  ( * ? ? ?
* theorem 2.2 ) .
we address the reader to the statements of  ( * ? ? ?
* theorem 3.1 and 3.2 ) and the toolboxes used in the @xmath1-closing and connecting lemmas which do not affect the regularity of the initial vector field .      by @xcite
, the existence of hyperbolic equilibria implies that the flow is not @xmath1near an anosov one and thus , @xmath1close to the flow , tangencies are expected . furthermore , by @xcite ,
arbitrarily @xmath1close to a conservative flow displaying a bykov cycle , we may find an incompressible flow with a dense set of elliptic periodic solutions .
nonetheless , these approximated flows might not have a bykov cycle anymore , due to the fragile connection @xmath115 $ ] .
we start with some definitions that make our results more precise .
let @xmath116 , @xmath117 be a small tubular neighbourhood of the heteroclinic cycle @xmath46 and @xmath118 be a cross section to @xmath48 $ ] .
we say that the two - dimensional invariant manifolds @xmath49 and @xmath50 have a tangency of order @xmath119 inside @xmath117 if the following conditions hold :    * the manifolds @xmath49 and @xmath50 meet tangencially ; * there is a solution contained in @xmath120 which lies entirely in @xmath117 and has exactly @xmath121 intersection points with @xmath118 .
a tangency of order 1 inside @xmath117 is also called by _
primary tangency_.    our first result is the following :    [ main1 ] there exists a @xmath1dense subset @xmath122 of the @xmath123open set @xmath3 such that for any @xmath124 and any tubular neighbourhood @xmath117 of the bykov cycle , the flow of ( [ general ] ) displays a primary tangency inside @xmath117 between the manifolds @xmath49 and @xmath50 .
the perturbations of theorem [ main1 ] are performed within a degenerate subclass of incompressible flows exhibiting bykov cycles .
the density of flows exhibiting tangencies is achieved without breaking the bykov cycle .
the proof of this result will be addressed in section [ sectiont - pointhetero ] , where we concentrate our attention on the geometric intersection between the invariant manifolds .
now , we consider a generalized version of the cocoon bifurcations described in @xcite , which can be seen as a series of global and local bifurcations that create an infinite number of heteroclinic connections .
the authors of @xcite called cocoon to these series of bifurcations because it forms a cocoon like structure in a section transverse to @xmath46 .
[ definition_cocoon ] let @xmath125 .
we say that @xmath10 exhibits a _
generalized cocooning cascade of heteroclinic tangencies _ if , for every @xmath126 large , there exists a closed two - dimensional torus @xmath127 such that :    1 .   for @xmath128 , @xmath129
the vector field @xmath10 can be @xmath1-approximated by @xmath130 whose
flow displays a tangency between @xmath49 and @xmath50 , intersecting twice @xmath127 and having length greater than @xmath131 within @xmath127
.    definition [ definition_cocoon ] is different to the general case given in @xcite in which the authors considered generic unfoldings of @xmath10 . in @xcite , when following a given sequence of vector fields converging to @xmath10 , a saddle - node bifurcation takes place , creating a pair of elliptical and hyperbolic periodic solutions , followed by an infinite cascade of period doubling bifurcations .
since there are infinitely many heteroclinic tangencies occurring near the primary tangency , in section [ cocoon_section ] we prove that :    [ main2 ] let @xmath3 be the @xmath1-open set of theorem [ main1 ] .
for any @xmath132 , and any tubular neighbourhood @xmath117 of the bykov cycle , the vector field @xmath10 exhibits a generalized cocooning cascade of heteroclinic tangencies .    the proof is based on @xcite and on the @xmath1-connecting lemma for volume - preserving flows  @xcite .
these heteroclinic bifurcations are responsible for the folding and for the fractal structure of the two - dimensional invariant manifolds in the cross section @xmath118 previously described by lay @xcite . the cocoon bifurcations might be observed numerically with the _ time delay function_. a discussion of the mechanism explaining the scape of points is given in ( * ? ? ?
* section 4 ) .    , emphasizing the region in which we apply the @xmath1-connecting lemma .
the effect of moving the invariant two - dimensional manifolds corresponds to the motion suggested by the arrow.,height=226 ]    in the context of dissipative diffeomorphisms containing homoclinic points , newhouse @xcite introduced the term _ wild attractor _ for hyperbolic sets whose invariant manifolds have a tangency .
newhouse described what happens in a one - parameter unfolding , when a tangency splits , and discovered nontrivial , transitive and hyperbolic sets whose stable and unstable invariant manifolds have persistent nondegenerate tangencies .
a given tangency may be removed by a small perturbation , however one can not avoid the appearance of new tangencies .
this theory can be extended to heteroclinic tangencies  see ( * ? ? ?
* section 3.2 ) .
[ elliptic ] let @xmath76 be a diffeomorphism defined on a compact section @xmath2 transverse to @xmath46 .
we say that @xmath24 is an _ elliptic periodic orbit _ of period @xmath133 for the diffeomorphism @xmath76 if the following conditions hold :    * @xmath134 and @xmath135 , for all @xmath136 and * the eigenvalues of the map @xmath137 are non - real and have modulus equal to 1 .
from now on , the map @xmath76 of definition [ elliptic ] can be thought as the first return map to @xmath118 ( induced by the flow ) . by the generalised stokes theorem @xcite , the return map of a volume - preserving vector field to any transverse section @xmath118 preserves an ` area - measure ' , and so it may be seen as an area - preserving diffeomorphism
see lemma [ area - preserving ] .
newhouse s results near tangencies have been extended to conservative maps in @xcite : any area - preserving map with a heteroclinic tangency can be @xmath138approximated by an open domain in the space of area - preserving maps exhibiting persistent tangencies .
generic diffeomorphisms on such a domain exhibit infinitely many elliptic islands , which take the role of attractors in the dissipative context .
let @xmath76 be the first return map to a section transverse to the cycle @xmath46 .
a 1-periodic solution of the flow of ( [ general ] ) is a periodic solution associated to a fixed point of @xmath76 ( also called by _ single round periodic solution _ ) .
we obtain the following result :    [ corollaryperioddoubling ] let @xmath132 .
the map @xmath76 defines a sequence of conservative horseshoes accumulating on @xmath46 .
there are persistent heteroclinic tangencies of the invariant manifolds associated to periodic solutions and infinitely many 1-periodic elliptic solutions nearby .
the realization of @xmath76 as return map associated to a divergence - free vector field @xmath1-close to the original follows the same lines as ( * ? ? ?
* section 3 ) .    * the method . *
our analysis follows a classic procedure .
we first construct a model of the dynamics near a bykov cycle in terms of poincar maps between neighbourhoods of the saddle - foci , and then we analyze the algebraic bifurcation equations . near each equilibrium
we derive local maps by assuming that the flow is governed by appropriate linearized vector fields in these regions .
the flow near the connections is approximated by linearizing about each one , which enables us to derive global maps .
appropriate composition of the local and global maps yields the desired return map .
our aim is to provide a direct and intuitive geometric picture and to uncover the scaling laws of the codimension - one bifurcations that occur near the cycle .
in this section , we establish local coordinates near the equilibria @xmath27 and @xmath28 and define some notation that will be used in sequel .
the starting point is an application of a quite useful conservative version of sternberg s theorem for @xmath1linearization proved by banyaga , de la llave and wayne @xcite to @xmath1linearize the vector field around the saddle - foci and to introduce cylindrical coordinates around them . for each saddle
, we obtain the expression of the local map that sends points in the border where the flow goes in , into points in the border where the flows goes out in positive time _
i.e. _ by forward iteration
. we also establish a convention for the transition maps from one neighbourhood to the other  see also ( * ? ? ?
* lemma 3.4 ) where it is proved the conservative flowbox theorem .
since the open and dense non - ressonance conditions of ( * ? ? ?
* theorem 1.2 ) are satisfied , then around the saddle - foci the vector field @xmath10 is @xmath1conjugate to its linear part . in cylindrical coordinates @xmath139 the linearizations at @xmath140 and @xmath141 take the form , respectively : @xmath142 and @xmath143    rescaling coordinates , we may consider cylindrical neighbourhoods of @xmath140 and @xmath141 in @xmath144 of radius @xmath145 and height @xmath146 that we denote by @xmath51 and @xmath52 , respectively , as illustrated in figure [ conservativo3 ] .
their borders have three components :    * the _ wall
_ parametrized by @xmath147 and @xmath148 with the usual cover : @xmath149 here @xmath58 represents the height of the cylinder and @xmath150 is the angular coordinate , measured from the point @xmath151 in the heteroclinic connection @xmath152 $ ] . * two disks , the top and the bottom of the cylinder .
we assume the fragile heteroclinic connection @xmath153 $ ] goes from the top of one cylinder to the top of the other , and we take a polar covering of the top disk : @xmath154 where @xmath155 and @xmath156 .    on these cross sections
, we will define the poincar maps to study the dynamics near the cycle .
consider the cylinder wall of @xmath51 , that meets the stable manifold of @xmath27 , denoted by @xmath157 , on the circle parametrized by @xmath158 .
the top part @xmath159 of the wall near @xmath27 will be denoted by @xmath160 . by construction , trajectories starting at interior points of @xmath161
go into the cylinder in positive time and come out at the cylinder top , denoted @xmath162 .
solutions starting at interior points of @xmath162 go outside the cylindrical neighbourhood in negative time . in these coordinates ,
the invariant manifold @xmath163 is the @xmath164axis , intersecting @xmath162 at the origin of the polar coordinates . reversing the time
, we get dual statements for @xmath28 .
after linearization , the one - dimensional invariant manifold of @xmath28 , @xmath165 , is the @xmath164axis , intersecting the top of the cylinder , @xmath166 , at the origin of its coordinates .
trajectories starting at interior points of @xmath166 go into @xmath52 .
trajectories starting at interior points of the cylinder wall @xmath167 go into @xmath52 in negative time .
the set @xmath168 is parametrized by @xmath158 .
trajectories that start at @xmath169 leave the cylinder @xmath52 at @xmath167 .     and @xmath51 near @xmath28 and @xmath27 , respectively .
the relative configuration of the saddle - foci is the same as in figure [ conservativo2].,height=309 ]      let @xmath89 be a regular point , @xmath170 and @xmath171 be the poincar map , where @xmath172 , @xmath173 , is the surface contained in @xmath2 whose tangent space at @xmath174 is generated by @xmath175 , @xmath173 , and @xmath176 is a small neighborhood in @xmath177 of @xmath150 .
we can guarantee the existence of a continuous time-@xmath178 arrival function @xmath179 from @xmath176 into @xmath180 .
if @xmath181 is sufficiently small then : @xmath182 is a thin flowbox with base @xmath183 and top @xmath184 . due to the conservative flow - box theorem
( see ( * ? ? ?
* lemma 3.4 ) ) , with a volume - preserving change of coordinates , we can trivialize @xmath185 into a cylinder and thus we may obtain an expression for the global map @xmath186 .
next lemma is the conservative version of a result stated in bykov @xcite and its proof directly follows by applying @xcite and @xcite :    [ perturb ] there are a volume - preserving change of coordinates in @xmath187 such that the linear part of the global map @xmath188 has the form : @xmath189 for some @xmath61 .
consider the transition maps : @xmath190    by lemma  [ perturb ] , the first order ( linear ) approximation of @xmath191 , up to a change of coordinates , it is equal to :    @xmath192    the transition @xmath193 can be seen as a rotation by an angle @xmath194 . as in @xcite , we assume that @xmath195 , simplifying our computations .
the flow is transverse to the previous cross sections and the boundaries of @xmath51 and of @xmath52 may be written as the closures of the disjoint unions @xmath196 and @xmath197 , respectively .
the trajectory of the point @xmath198 in @xmath199 leaves @xmath51 at @xmath162 at : @xmath200 where @xmath201 . in a dual way
, points @xmath202 in @xmath203 leave @xmath52 at @xmath167 at : @xmath204 where @xmath205 .
the minus sign in the equation @xmath206 of suggests the condition stated on hypothesis * ( h4 ) * about different chirality .
now we describe the geometry associated to the local dynamics near each equilibrium .
we start with some useful definitions adapted from @xcite .
a continuous real valued map defined on an interval @xmath207 $ ] is _ quasi - monotonic _ if it is monotonic for a subinterval @xmath208\subseteq [ 0,1]$ ] where @xmath209 .
we shall define three objects to be used in the sequel : segments , spirals and helices :    1 .
[ def_segmento ] a _ segment _ @xmath210 in @xmath160 or @xmath167 is a smooth regular parametrized curve of the type @xmath211\rightarrow in(\sigma_1 ) \qquad \text{or } \qquad \beta : ( 0,1]\rightarrow out(\sigma_2)\ ] ] that meets @xmath212 or @xmath213 transversely at the point @xmath214 only and such that , writing @xmath215 both @xmath150 and @xmath58 are quasi - monotonic and bounded functions of @xmath216 and @xmath217 is bounded .
[ def_spiral ] a _ spiral _ in @xmath162 ( resp . :
@xmath166 ) around a point @xmath218(resp . :
@xmath219 is a curve @xmath220\rightarrow out(\sigma_1 ) \qquad \text{or } \qquad \alpha : ( 0,1]\rightarrow in(\sigma_2)\ ] ] satisfying @xmath221 and such that , if @xmath222 are its expressions in polar coordinates @xmath223 around @xmath224 , then @xmath225 and @xmath226 are quasi - monotonic , with @xmath227 , i.e. , winds infinitely many times around the point @xmath224 .
3 .   consider a cylinder @xmath228 parametrized by a covering @xmath229 $ ] , with @xmath230 where @xmath231 is periodic .
helix _ in the cylinder @xmath228 _ accumulating on the circle _ @xmath232 is a curve @xmath233\rightarrow c$ ] such that its coordinates @xmath234 satisfy @xmath235 and the map @xmath236 is quasi - monotonic .
the next result summarizes some basic ideas on the spiralling geometry of the flow near the saddle - foci and has been proved in ( * ? ? ?
* section 6 ) . since in its original form
the authors assumed implicitly that the property * ( h4 ) * does not hold , we make the necessary reformulations .
[ structures ] a segment @xmath210 :    1 .   in @xmath160 is mapped by @xmath237 into a spiral in @xmath162 around @xmath238 ; 2 .   in @xmath167
is mapped by @xmath239 into a spiral in @xmath166 around @xmath240 ; 3 .
[ item3 ] in @xmath160 is mapped by @xmath237 into a spiral in @xmath162 around @xmath238 , that is mapped by the conservative transition @xmath241 into another spiral , around @xmath242 .
this new spiral suffers the effect of the hyperbolic saddle distortion given by lemma  [ perturb ] .
4 .   if * ( h4 ) * does not hold , the spiral defined on ( [ item3 ] ) is mapped by @xmath243 into a helix in @xmath167 accumulating uniformly on the circle @xmath244 .
the complete study where hypothesis * ( h4 ) * does not hold has been done in @xcite for dissipative systems .
we assume that @xmath245 and @xmath246 define vertical lines across the cylinder walls of @xmath51 and @xmath52 , @xmath160 and @xmath167 , respectively . let @xmath247 , @xmath248 $ ] ,
be a parametrization of @xmath245 , where @xmath249 are the coordinates of the point @xmath48
\cap in(\sigma_1)$ ] .
let @xmath250 be the coordinates of @xmath251 .
the set @xmath252 ( @xmath253 ) defines an oriented curve in @xmath167 . under hypothesis *
( h4 ) * , it changes the direction of its turning around the cylinder @xmath167 at infinitely many points where it has a vertical tangent . in
what follows , we obtain the expression for @xmath254 , @xmath253 .
the proof runs along the same lines as in @xcite .
[ x , y ] let @xmath255 , @xmath248 $ ] , be a segment in @xmath160 parametrized by @xmath216 . then , for @xmath253 :    @xmath256)+ \phi(\varphi)\\   \\
y_2(s ) =    s [ c(\varphi)]\\ \end{array } \right.\ ] ]    where @xmath257 and @xmath258 with the argument @xmath259 taken in the interval @xmath260 $ ] , @xmath261 , that contains @xmath57 .    from now on ,
denote by @xmath262 the following map @xmath263     ( height ) and ( b ) ; the coordinate @xmath264 ( angular coordinate ) , with @xmath265 , using maple .
parameters : @xmath266 , @xmath267 and @xmath268.,height=264 ]      if @xmath269 , then @xmath270 and @xmath271 .
moreover , @xmath272    if @xmath273 as in the reversible case reported in @xcite , then @xmath274 is constant and the segment @xmath275 is mapped under @xmath262 into a vertical segment .      as suggested in figure [ maple1](a ) , if @xmath265 , then the coordinate @xmath276 is not a monotonic function of @xmath216 but @xmath277 .
moreover , for any @xmath278 , we have : @xmath279 and @xmath280    under hypotheses * * ( h1)***(h5 ) * , we show that the coordinate map @xmath274 is not a monotonic function of @xmath216 , since the curve @xmath281 reverses the direction of its turning around @xmath167 infinitely many times  see figure [ maple1](b ) .
this is the notion suggested by the following definition .
we say that the vector field @xmath10 has the _ dense reversals property _ if for the vertical segment @xmath282,\ ] ] the projection into @xmath283 of the points where @xmath284 has a vertical tangent is dense in @xmath285 .
the dense reversals property is the key step in the proof of theorem [ main1 ] . in order to prove it
, we use the assumption * ( h5 ) * on the parameters @xmath286 that determine the linear part of the vector field @xmath10 at the equilibria .
, we observe the line @xmath49 reversing the orientation of its angular coordinate .
the set @xmath285 is diffeomorphic to a circle.,height=188 ]    the proof will be divided into five parts , * ( i ) *  * ( v)*. the first part of the proof runs along the same lines to that of @xcite , which gives an optimal condition for the existence of tangencies .
let @xmath132 . *
( i ) * we need to compute the coordinate @xmath264 at the points where @xmath287 has a vertical tangent . differentiating the expression of lemma [ x , y ] with respect to @xmath216 , we get : @xmath288 and then @xmath289.\ ] ]
therefore @xmath290 has solutions if and only if @xmath291 define now the @xmath292periodic map : @xmath293    * ( ii ) * now we are going to prove that the left part of ( [ eqdxw ] ) has a zero if and only if @xmath294 indeed , the expression @xmath291 is equivalent to @xmath295 using the definition of @xmath296 and @xmath297 , the last equality yields the expression : @xmath298 which is the same as @xmath299    therefore the differential equation @xmath290 has solutions if and only if @xmath300 corresponding to the condition that defines @xmath70 . if @xmath301 $ ] is a solution of ( [ eqdxw ] ) , then there are infinitely many given by @xmath302 where @xmath303 . since @xmath304 , then @xmath290 has solutions @xmath305 where @xmath306
. see figure [ reversions ] .
* ( iii ) * for any @xmath307 and @xmath308 , we have @xmath309    , we observe the @xmath49 reversing the orientation of its angular coordinate .
two consecutive reversion points lie in region of the graph with different concavities , diving rise to different types of tangencies ( and therefore to different types of dynamics ) .
, height=151 ]    using the expressions we get that @xmath310 also , if @xmath311 $ ] , then @xmath312 lies in the same interval and @xmath313 the result follows , using the expression for @xmath264 in lemma [ x , y ] . in part * ( iii ) * , we have shown that @xmath314 hence , if the genericity condition @xmath315 in holds , then the points @xmath316 are dense in the circle defined by @xmath317 and thus @xmath10 has the dense reversals property .
with respect to the coordinate @xmath274 , in the wall of the cylinder @xmath52 , two consecutive reversion points lie in region of the graph with different concavities , as illustrated on figure [ wall ] .
this phenomenon implies the existence of chaotic dynamics near the tangency as we proceed to explain .    *
( iv ) * if @xmath318 is the projection of a point where @xmath284 has a vertical tangent , then @xmath50 is tangent to @xmath49 at the corresponding point .
otherwise an arbitrarily small modification of the vector field @xmath10 in a neighbourhood of @xmath319 will move @xmath318 , creating the tangency  see figure [ cocoon ] .
suppose @xmath320 is a curve close to a vertical segment and that it is parametrized by @xmath321with @xmath322 then , changing @xmath76 near @xmath319 if necessary , we may assume that @xmath323 meets @xmath281 at a point where the last curve has a vertical tangent .
the perturbation is achieved in part * ( v)*.    * ( v ) * in  @xcite is it used an unfolding argument to obtain tangencies from the dense reversals property . since we have no guarantee that a generic unfolding process holds in the volume - preserving setting , we should tackle the problem using another approach .
the key step is the volume - preserving hayashi s connecting lemma ( * ? ? ?
* and theorem e ) , whose formulation may be adapted for flows as follows :    ( * ? ? ?
* adapted ) [ wex ] let @xmath324 be a non - periodic solution of @xmath325 . for any @xmath1-neighbourhood @xmath326 of @xmath10 ,
there exists @xmath327 and @xmath328 such that for all @xmath329 $ ] and any two points @xmath330 outside the tube @xmath331}\varphi(t,(b(z,\delta))$ ] , if the forward trajectory of @xmath224 intersects @xmath332 and the backward trajectory of @xmath333 intersects @xmath332 , then there exists @xmath334 such that @xmath335 outside @xmath331}\varphi(t,(b(z,\delta))$ ] and there exists a one - dimensional connection @xmath336 $ ] .
the previous result says that if two distinct initial conditions @xmath224 and @xmath333 have solutions that visit a given neighbourhood of a point @xmath337 and the initial conditions @xmath224 and @xmath333 are sufficiently far away from a piece of the backward solution of @xmath164 , then we can find a divergence free vector field for which @xmath224 and @xmath333 are heteroclinically related .
the support of the perturbation lies on @xmath338}\varphi(t,(b(z,\delta)).\ ] ] by * ( i)(iii ) * , we clearly have two points @xmath339 and @xmath340 in the conditions of theorem [ wex ] .
it is obvious that the solution associated of @xmath164 can not intersect the heteroclinic connection @xmath29 $ ] .
therefore , we can always choose @xmath341 sufficiently small in order that the sets @xmath342}\varphi(t,(b(z,\delta))$ ] and @xmath29 $ ] are disjoint .
then , we can perform the volume - preserving perturbation given by theorem  [ wex ] to obtain the tangencial intersection between @xmath49 and @xmath50 by making a connection near @xmath164 . if the intersection is tangential the proof is completed .
otherwise , we consider a smooth path @xmath343}$ ] of elements on @xmath3 such that @xmath344 and @xmath345 .
since @xmath346 displays no intersection between @xmath347 and @xmath348 , @xmath349 displays a transversal intersection between @xmath350 and @xmath351 , and the accumulation is of a quadratic type , then there must be @xmath352 such that @xmath353 displays a tangencial intersection between @xmath354 and @xmath355 .
the classic cocoon bifurcation introduced in @xcite ` begins ' when the two - dimensional invariant manifolds of @xmath27 and @xmath28 have a tangency .
after the first tangency , when moving slightly a given one - dimensional parameter , we obtain two structurally stable heteroclinic connections . it is the beginning of a sequence of heteroclinic bifurcations .
saddle - node bifurcations and elliptic solutions have been detected numerically in @xcite .
we would like to prove that , for any @xmath356 , and any tubular neighbourhood @xmath117 of the bykov cycle @xmath46 , for every @xmath126 , there exists a closed two - dimensional torus @xmath127 such that @xmath357 and @xmath10 can be @xmath1-approximated by @xmath130 whose flow displays a tangency between @xmath49 and @xmath50 , intersecting twice @xmath127 and having length greater than @xmath131 within @xmath127 .
theorem  [ main1 ] shows that the vector field of ( [ general ] ) can be @xmath1-approximated by another vector field such that the associated flow displays a primary tangency inside @xmath117 between the manifolds @xmath49 and @xmath50 . beyond the tangency , there are infinitely many pieces of @xmath49 accumulating on @xmath358
. therefore the set @xmath359 accumulates on @xmath360 . repeating recursively the argument of the proof of theorem [ main1 ] in the appropriate direction
, we get curves in @xmath167 with the dense reversals property . using again the @xmath1-connecting lemma as in * ( v ) * , we may construct vector field @xmath130 such that its flow exhibits tangencies of order @xmath361 , with @xmath362 large , and then we may define a torus @xmath127 , not containing the saddle - foci , inside which heteroclinic tangencies of order @xmath361 spend most of the time ( say larger than @xmath126 ) .
an illustration of this phenomenon has been given in figure [ cocoon ] .
the main tool for detecting the manifolds @xmath49 and @xmath50 is the time delay function , defined by : @xmath363 where @xmath364 corresponds to integrate along the trajectory @xmath365 the time delay function measures how long it takes for a given solution to escape from the sphere of radius @xmath366 .
trajectories starting at the regions between two consecutive lines of figure [ cocoon ] ( corresponding to heteroclinic tangencies ) will turn one more time around the equilibria @xmath28 - more details in ( * ? ? ?
* section 7 ) .
this is why the map @xmath364 acquires a castle shape with a higher level between pairs of equilibria of the time delay map .    in the present paper
, we can not use the arguments of @xcite in which lebesgue almost all solutions leave @xmath117 and hence will escape from a given sphere centered on the origin , with radius @xmath367 .
the later paper uses explicitly that hypothesis * ( h4 ) * does not hold . by detecting logarithmic equilibria of this map ,
one locates the invariant manifolds .
a typical time delay function has several singularities as the one depicted in ( * ? ? ?
11 ) , corresponding to the heteroclinic tangencies . the time delay function acts as a numerical mechanism , searching along lines in a cross section , for pieces of the invariant manifolds of the two nodes .
heteroclinic connections form a skeleton for the structure of periodic and bounded aperiodic orbits . according to @xcite , when elliptic periodic solutions are detected , the time delay function has flat tops .
numerical results have been presented in ( * ? ? ?
* section 4 ) .
the analytical existence of these solutions is the core of next section .
the organizing centre for the dynamics may be characterized by the presence of a hyperbolic set which is flow - invariant , indecomposable and topologically conjugate to a bernoulli shift space with finitely many symbols .
the union of all of these sets is not uniformly hyperbolic , and countably many of these sets are destroyed under small generic perturbations .
homoclinic cycles of shilnikov type and subsidiary cycles are expected to occur @xcite .      in this section
we study the existence of elliptic fixed points for the map @xmath76 .
we will make use of the references @xcite .
the next result will be useful to make the bridge between conservative vector fields and conservative first return maps .
we thanks santiago ibez for the hints on the next proof .
[ return_map ] [ area - preserving ] the return map of a volume - preserving vector field to any transverse section @xmath118 preserves the area @xmath368 ( for some density @xmath369 ) , and so is an area - preserving diffeomorphism .
let @xmath10 be a volume preserving vector field , _ i.e. _ such that @xmath370 . take a transverse section @xmath118 with a well defined poincar map @xmath371 , where @xmath372 .
now , take a simply connected subset @xmath373 ( assume it is a disk to make it simple ) .
consider the solid @xmath2 obtained by saturation of the flow from @xmath374 to @xmath375 .
the boundary of @xmath2 , @xmath376 , can be written as @xmath377 , where @xmath378 is the saturation by the flow of the one - dimensional boundary @xmath379 of @xmath373 until it reaches @xmath380 . by the ostragadsky - stokes theorem @xcite ,
the triple integral of @xmath381 on @xmath2 equals the surface integral of @xmath382 on @xmath376 .
more formally : @xmath383 since @xmath384 , the left part of the previous equality is zero .
moreover , @xmath385 because @xmath10 is tangent to @xmath378 at each point .
hence @xmath386    in the above expression , the normal vector on @xmath374 is the opposite to the normal vector on @xmath375 because the normal vector must have a continuous variation on @xmath376 . by taking the normal vector pointing in the same direction
the conclusion is that @xmath387 meaning that there exists ` an adapted ' measure that is preserved ( by @xmath388 ) .
we will borrow the arguments of @xcite to show the existence of suspended horseshoes near the heteroclinic tangency .
the latter authors compute the first return map @xmath76 and apply the rescaling lemma for studying bifurcations of fixed points .
the rescaling lemma says that we may reduce the study of bifurcations of the first return maps to the classic analysis of bifurcations of the two - dimensional conservative hnon - like maps .    by @xcite
, we may find @xmath389 and a union of strips @xmath390 , @xmath391 , accumulating on the segment @xmath392 as @xmath393 such that @xmath394 has a conservative horseshoe near the cycle , for @xmath395 . for small perturbations of @xmath76 ,
infinitely many bifurcations of horseshoes creation or destruction occur , including the birth and disappearance of elliptic periodic points .
these horseshoe bifurcations must have different scenarios depending on a type of the initial tangency .
the character of reciprocal position of the strips @xmath396 and their images @xmath397 are essentially defined by the signs of two parameters that governs the linear part of the global map ; the authors of @xcite selected six different cases of maps with quadratic homoclinic tangencies . for these parameters and for @xmath119 large
, the position of all the strips @xmath396 and @xmath397 have regular intersections .
in particular , if @xmath398 , then the first return map associated to any vector field @xmath1-arbitrarily close to @xmath10 has a conservative horseshoe accumulating on @xmath399
. moreover , there are persistent heteroclinic tangencies of the invariant manifolds of periodic solutions .    in a small section transverse to the cycle ( not containing it ) , since the number of symbols of the bernoulli shift approaches @xmath400 as the strip approaches the cycles , a cascade of bifurcations occurs , associated to the creation and annihilation of horseshoes @xcite .
there are no guarantees that such perturbations still have the bykov cycle in its flow .    if @xmath401 is a linear map such that the @xmath402 , then @xmath403 if and only if its eigenvalues are complex ( non - real ) and conjugated .
let @xmath356 .
any neighbourhood of the points where we observe dense reversals property contains nontrivial hyperbolic subsets including infinitely many horseshoes , accumulating on the tangency .
in particular , there are infinitely many 1-periodic solutions accumulating on the tangency .
let @xmath404 and @xmath405 be the first return map of a vector field on @xmath2 satisfying ( * h1*)(*h5 * ) .
the points @xmath406 for which we observe reversals satisfy @xmath407 observing that @xmath408 and @xmath409 is bounded , it follows that : @xmath410    near the initial conditions in @xmath160 where we observe the dense reversals property , condition ( [ tang_final ] ) holds and thus , we conclude that @xmath411    as illustrated in figure [ trace ] , for @xmath412 , if @xmath413 is a neighbourhood in @xmath167 of a point for which we observe a reversal , then there is a line where the trace is zero and that splits @xmath413 into two connected components where the trace is positive or negative . in particular , we may define a strip @xmath414 in which the trace belongs to the interval @xmath415 , where all 1-periodic solutions should be elliptic .
intersecting the set of 1-periodic solutions that accumulate on the tangencies with the strip @xmath118 , yields infinitely many elliptic 1-periodic solutions accumulating on the tangency .
yields infinitely many elliptic 1-periodic solutions accumulating on the tangency . ,
height=113 ]      dynamics near an elliptic point is reminiscent of the dynamics near a center type equilibrium : under regularity assumptions on the vector field , a generic elliptic point is surrounded by a family of closed kam - curves which occupy a set of almost full lebesgue measure in a small neighbourhood of the fixed point . on each kam - curve
the map @xmath76 is conjugate to an irrational rotation . as the distance to the fixed point increases , the rotation number changes in a monotonic way .
the kam - curves do not fill the neighbourhood of the elliptic point .
rational rotation numbers correspond to resonant birkhoff zones between the closed curves and the dynamics in these regions is typically chaotic .
this gives rises to chaotic dynamics between tori regions coexisting with hyperbolic horseshoes and the bykov cycle .
adding the reversibility to the problem , symmetric elliptic trajectory may be seen as the limit of other elliptic periodic points
see gonchenko _
et al _ @xcite
.    when the fragile connection is broken , near the bykov cycle there are a wide range of curious homoclinic bifurcations as reported in @xcite .
some of them do not depend on the chirality and are generated by the swirling pattern due to the existence of saddle - foci
. they can be viewed as spiralling structures in a parameter diagram .
conservative structures might be comprised of elliptic islands separated by saddles .
there is another self - similar organization of embedded saddles and centers on smaller scales .
the closer one approaches a center of the rings , the greater the number of twists the outgoing heteroclinic connection makes before returning to the saddle
see @xcite .
in this paper , trying to answer the questions * ( q1)**(q2 ) * , we have proved that the extended version of the cocoon bifurcations for conservative systems , the heteroclinic tangencies and the elliptic periodic solutions are strongly related .
more precisely , in the present paper we have considered divergence - free vector fields whose flows have bykov cycles , we have found heteroclinic tangencies leading to the description of elliptic periodic solutions for the first return map to a cross section defined near the tangency . within a codimension 1 set in @xmath73 ,
tangencies between the invariant manifolds of the two saddle - foci densely occur .    the proof of the existence of elliptic solutions near the horseshoes does not run along the same lines of gonchenko and shilnikov @xcite , in which the authors used renormalisation to the hnon map .
along the proof of our results , we use the topological notion of chirality , making our findings different from those of @xcite .
we asked for the @xmath1 volume - preserving connecting lemma @xcite to perform local perturbations .
similar heteroclinic bifurcations have been considered in @xcite , where generic flat perturbations do not preserve the cycle .
the non - wandering dynamics near the bykov cycle is dominated by hyperbolic horseshoes , conjugate to a full shift over a finite alphabet , that accumulate on the cycle .
hyperbolic behaviour can not be separated from the ellitic behaviour by isotopies .
our underlying hypothesis on the maintenance of the bykov cycle is inspired in the @xmath76-reversible case , with @xmath416 , where the mantainance holds _ per se_.
preserving the bykov cycle is the most important consequence from the anti - symmetry ; this is why we decided to reach general results without the reversibility hypothesis .
the divergence - free character is crucial because it guarantees the stability of elliptic periodic solutions in our three - dimensional context . in
higher dimensional case , the divergence freeness is no longer sufficient and we should consider hamiltonian vector fields in order to get the stability of elliptic periodic solutions .    in the present paper
, we are able to prove the existence of heteroclinic tangencies between @xmath49 and @xmath50 .
our findings are richer that those of @xcite in the sense that we can explicitly give the manifolds for which we observe tangency .
note that the works @xcite proved the existence of a @xmath1-dense / residual set inside which tangencies and elliptic solutions occur via a mechanism developed by mora and romero @xcite to create open sets having a dense set of maps with tangencies nearby .
different chirality of the nodes is a new topological concept and has never been studied in the conservative setting .
this topological property has profound effects on the dynamics near the bykov cycle . of course
, in agreement with kaloshin s theorem on the prevalence of kupka - smale systems , these tangencies can not occur for a full lebesgue measure on @xmath2 .
the density of tangencies near the cycle is markedly different from the results about the michelson system proved in @xcite where the tangencies ( near the bykov cycle ) are rare . if @xmath77 , the @xmath76-reversibility prevents different chirality of the nodes .
the context of the proof of theorem [ main1 ] is different from that described in knobloch _
et al _ @xcite using spirals , a work that has been motivated by the michelson system . instead of restricting the flow to the wall of the cylinder ( as we do in theorem [ main1 ] ) ,
their conclusions , like bykov s @xcite , are supported by the study of spirals intersections . in (
* lemmas 6.6 and 6.7 ) it is shown that two logarithmic spirals with a common centre can have at most two tangencies , assuming an irrational ratio of the real part of the complex eigenvalues
see figure [ spirals](a ) .
although these spirals are only the leading order terms of the bifurcation equations , the density of tangencies can not occur near the cycle , as shown in @xcite . without breaking the fragile connection , assuming the same chirality of the nodes , tangencies may occur but not near the cycle .
bykov implicitly assumed that the chiralities of two nodes are different in formulas ( 3.4),(3.5 ) of @xcite .
indeed , observe the same sign of the exponent of @xmath417 in the formulas that represent the evolution of the angular component of the spirals : @xmath418 where @xmath419 represents the angular coordinates in @xmath162 and @xmath166 and the nonreal eignevalues are given by @xmath420 , @xmath421 . the same happens in formula ( 3.1 ) of bykov @xcite .
after taking logarithms of both sides , from equality ( 3.3 ) on , bykov neglects the assumption about the chiralities of the nodes , assuming that they are different . therefore , in the cross sections , the spirals corresponding to the two - dimensional invariant manifolds of the saddle - foci are oriented in the same way .
this explains the orientation of the spirals of figure 2 of @xcite in contrast to those depicted in figure 11 of @xcite .
these spirals are revived in figure [ spirals ] .
bykov never comments on the chiralities of the nodes , assuming implicitly that they are different .     and
@xmath422 have : case ( a ) : an infinite number of transversal intersections and at most two tangential intersections in total @xcite .
case ( b ) : an infinite number of transversal intersections and , for small perturbations which do not break the fragile connection , an infinite tangential intersections.,height=170 ]      the cocoon bifurcation is a set of rich bifurcation phenomena numerically observed by lau @xcite in the michelson system .
the authors in kokubu _
et al _ @xcite proved the existence of the cocoon bifurcations for the michelson system , which arise as the travelling wave equation of the kuramoto - sivashinky equation and also as a part of the limit family of the unfolding of the three - dimensional nilpotent singularity of codimension 3  see @xcite .
numerical results show that , when a given parameter varies , the system exhibits an infinite sequence of heteroclinic bifurcations due to the tangency of the two - dimensional invariant manifolds of saddle - foci .
each bifurcation creates new transverse heteroclinic connections . in @xcite ,
these bifurcations have been studied from a general point of view ; they have been explained as a consequence of an organizing center called the _
cusp transverse heteroclinic chain_. fixing the fragile connection , the cocoon bifurcations play the same role as the newhouse bifurcations ( generic unfoldings for the homoclinic tangency ) in the dissipative category .
thick cantor sets appear near the tangency , which are @xmath8-persistent , @xmath423 . in the present paper , we prove the existence of cocoon bifurcations within a degenerate set of divergence - free vector fields characterized by different chirality of the nodes .
in contrast to the findings of @xcite , in which the authors could break the cycle , our perturbations have been performed without breaking the fragile connection .
furthermore , in @xcite , the reversibility is used to get saddle - nodes with both stable and unstable sets of dimension 1 and the volume - preserving character plays no role .
we may use our results to illustrate the expected behaviour of the simple bifurcations of symmetric periodic solutions , namely saddle - node and period doubling via the creation and annihilation of horseshoes in a similar phenomenon to that described in @xcite .
although we deal with a degenerate class of vector fields , our results hold for any generic unfolding of the cycle .
the greatest difficulty arises when we want to preserve the bykov cycle , a property that emerge in @xmath76-reversible systems with @xmath416 . in a broader way
, our results hold for any kind of cycles having a one - dimensional solution connecting two saddle - foci .
in fact , following the approach in @xcite , it should be possible to show that vector fields containing bykov cycles are always @xmath1-close to another with a tangency ( and even a cocconing cascade of tangencies ) , without making any assumption on the chirality .
nevertheless , there is no guarantee that the cycle is preserved .    although a degenerate class of systems satisfy * ( h1)(h5 ) * , the present work may be seen as the description of a simple and general mechanism to create non - hyperbolic behaviour in a robust fashion near a singular cycle .
since the existence of cycles involving equilibria is not a generic phenomenon , they may exist outside the residual sets given by the classical dichotomies for divergence - free vector fields @xcite .
several questions remain to be solved for this kind of cycles , in particular about the existence and sign of the lyapunov exponents .
the lebesgue measure of the points that remain for all time in the neighbourhood of the cycle is an interesting point to be considered . in the spirit of @xcite ,
conservativeness - breaking is a another direction that should be investigated .
a first attempt to understand this phenomenon may be found in dumortier _
et al _ ( * ? ? ?
* section 4 ) .
the present article should be seen as a starting point for further related studies ; the ergodic characterization of the dynamics near these cycles ( even in the dissipative case ) is still far for being completely understood .
gonchenko , s. v. , shilnikov , l. p. & turaev , d. m. ,  homoclinic tangencies of an arbitrary order in newhouse domains  , _ itogi nauki tekh . ,
sovrem . mat .
_ , * 67 * , ( 1999 ) 69128 [ english translation in j. math .
, 105:1738 - 1778 , 2001 ] .
labouriau , i .
s. ; & rodrigues , a. a. p. , `` partial symmetry breaking and heteroclinic tangencies '' , in : s. ibez , j.s .
prez del ro , a. pumario , j.a .
rodrguez ( eds . ) , progress and challenges in dynamical systems , in : springer proc . math .
stat . , springer - verlag , ( 2013 ) , 281299 .