Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: almost all of us have already met someone far from home , who turns out to share a friend with us . milgram first put this social phenomenon on a firm basis , finding that any two individuals are separated , on average , by six acquaintances @xcite . furthermore , our society is composed of groups within which , if individual a is acquainted with persons b and c , then b is likely to know c. in the language of social networks , this means that the individuals a , b and c are arranged on a triangle . these two facts tell us that in the social universe , global distances are small , and local interconnectedness , i.e.  _ clustering _ , is high . a network that possesses both of these topological properties is called a _ small world_. in order to find a model that accounts for both of these properties , one could choose a regular lattice whose nodes are indeed locally highly interconnected . but since global distances are large in this type of network , it is an ineligible candidate for a small world model . in a random graph @xcite , that is a set of nodes ( e.g.  laid out on a virtual circle ) with connections between them established at random , a link is more likely to point to a far away node than to one close by . ( here we refer to the distance on the virtual circle s circumference . ) as a consequence , the presence of the very many long - range links account for the desired global property , but the degree of clustering is low . watts and strogatz combined these two insights and proposed a model that interpolates between a regular lattice and a random graph , thus capturing both the global and local topological properties mentioned above @xcite . but the small world property is not merely exhibited by social networks . yet , through the increased availability of data , it was found that the simultaneous occurence of high local _ and _ global interconnectedness is prominent to a much wider class of systems , notably the internet @xcite , the world wide web @xcite , metabolic networks @xcite and food webs @xcite . in addition to the small world property , there is another important fact when it comes to characterizing the topology of a complex network : not all the nodes have the same number of edges . the corresponding measure is the _ degree distribution _ @xmath0 which gives the probability that a randomly chosen node has degree @xmath1 , that is @xmath1 edges . for most of the just mentioned examples , this distribution was found to be @xmath2 , @xmath3 , implying the absence of a characteristic ( degree ) scale , hence the name _ scale - free _ network . this emergence of scaling can be understood , for example , in terms of _ growing _ networks . starting from a small core graph , at each time step a node is added together with a certain number of edges that are connected to existing nodes , the latter being chosen _ preferentially _ , that is a link is more likely established to a high degree node @xcite . this fat - tailed degree distribution implies the presence of _ hubs _ ( high degree nodes ) which hold together the network and play a crucial role in issues such as robustness or fragility @xcite . this explosion of research activity in the field of complex networks has also shed new light on epidemic spreading since the latter can be regarded as a dynamical process occuring on a complex network : a computer virus spreads on the internet , and hiv ( as a biological example ) propagates on top of the web of human sexual contacts . also these two examples fall into the category of scale - free networks @xcite , and it is this topological property that accounts for the absence of a finite epidemic threshold in the corresponding spreading phenomenon @xcite . the degree distribution is only a first way to characterize the degree related topology of a complex network . indeed , by analyzing scientific collaboration networks , researchers working with many others ( high degree nodes ) tend to collaborate with other `` hubs '' . this means that there exist _ degree correlations _ , and the just mentioned property has been called _ assortative mixing _ , holding generally for social networks @xcite . on the other hand , in the internet ( at the autonomous system level ) , high degree nodes are more likely connected to low degree ones , thus exhibiting _ disassortative _ mixing . the influence of such degree correlations on the spreading of an epidemic was investigated in detail @xcite , finding that neither assortative nor disassortative scale - free networks exhibit a finite epidemic threshold @xcite . the insights about the role of these degree related topological properties in epidemic spreading have been gained at the mean - field level .    besides degree correlations , triangles are ubiquitous in complex networks as outlined in the first paragraph , and more generally , many loops of short length were found in these systems @xcite . motivated by watts and strogatz model @xcite which uses a regular lattice , i.e.  an ordered network possessing many loops , as starting point , we shall also use the concept of _ local ordering _ when referring to the loop structure of a complex network . more information can be extracted if the triangles are sorted according to the degrees of their corners , finding that mainly the low degree nodes account for the high level of clustering @xcite . this suggests the presence of interesting modular organizations , and similar results were obtained from an analysis of loops of length 4 @xcite . local ordering plays a crucial role in the function of a metabolic network ( with scale - free topology ) . indeed , the loop structure is much richer than in the scale - free model based on growth and preferential attachment @xcite . this work also investigates the number of triangles as a function of the system size for this model , a result which was generalized to loops up to length 5 , yielding robust scaling relations @xcite . obviously the presence of loops has an effect on the spreading behavior since , with respect to a treelike topology , there exist many more paths along which the virus can propagate . different strategies in order to gain insights about the role of local ordering properties have been proposed . an interpretation of how clustering influences the stationary spreading behavior was obtained by mapping the epidemic process onto bond percolation @xcite . another approach is to abandon the mean - field level and take into account spatial correlations which govern the epidemic dynamics . matsuda et al . first used the ordinary pair approximation in order to study a population dynamical problem @xcite . this approximation , as its name anticipates , accounts for pair correlations and lies at the basis of improved pair models @xcite which uncover the role of local ordering in a rather indirect way : clustering enters by making a number of assumptions about the open ( @xmath4 ) and closed ( @xmath5 ) triple correlations . in cluster approximations , a time - dependent probability is assigned to each configuration of the fundamental cluster whose choice is guided by the network topology @xcite : for investigating the spreading dynamics on a triangular lattice one uses a triangle as fundamental cluster whereas the star is the appropriate choice for a random network . higher - order correlations are therefore embedded very explicitly . moreover the systematic improvability of this method makes it a powerful tool to study probabilistic systems .    in order to understand the role of local ordering properties , the exploration of temporal correlations seems to be an even more natural approach : for example within a two - step description , it matters whether the local topology is treelike or if loops of short length are present . the method is illustrated for the susceptible - infected - susceptible model ( see for example @xcite ) , homogeneous networks are used as starting point , and analytical estimates are obtained also for disordered graphs obeying @xmath6 , @xmath7 being arbitrary . the paper is organized as follows . section ii describes the adopted model consisting of the contact network as well as the local dynamics . in section iii , we introduce the formalism , from which the two - step description is derived . for completeness , the necessary ingredients in order to arrive at the one - step site approximation , i.e.  the common mean - field level , are also shown . section iv explores the implications of the double - step approach for networks of degree 4 , the generalization to arbitrary degree is done in the following section . the major conclusions are drawn in section vi . the dynamical laws that describe the spreading of an infectious disease are determined by the contact structure which underlies the population . we therefore model the epidemic as a dynamical process on top of a given network that does not change in time . the nodes of the network represent individuals , and the links correspond to relationships between individuals along which an infective agent can propagate . since the aim of this paper is the investigation of the role of loops of short length , we adopt a rather simple epidemiological model where the individuals can be only in two possible states , namely infected ( i ) or susceptible ( s ) . because the nodes repeatedly run through the cycle susceptible @xmath8 infected @xmath8 susceptible , it is called sis model . in the physics community , it has recently been formulated as follows @xcite : a node susceptible to the disease gets infected with probability @xmath9 if it is connected to at least one infected nearest neighbor . on the other hand , infected nodes recover spontaneously with probability @xmath10 . this version of the sis model is formally advantageous with respect to its conventional formulation , where infected nodes can infect neighboring susceptible vertices with probability @xmath9 @xcite . in the latter case , susceptible nodes become infected with probability @xmath11 , @xmath12 being the number of infected nearest neighbors . in this paper , we will use the former version . by rescaling the time unit , we can reduce the number of parameters to one : the time evolution is determined by the effective spreading rate @xmath13 , and the recovery rate is set to 1 . the quantitative details of the behavior of the system still depend on the choice of @xmath14 . in particular , the effect of the loops are of higher order in @xmath14 , such that their influence is not seen in the continuous - time limit ( @xmath15 ) . as long as @xmath16 , we set this quantity to 1 without lack of generality . the other model constituent concerns the underlying network . we shall not attempt to examine the combined effect of the degree distribution , degree correlations and the loop structure . while the role of these degree related connectivity patterns in epidemic spreading is rather well understood , little attention has been given to the effect of local ordering . that is why we focus on this topological property . we therefore use networks in which every node has the same number of nearest neighbors ( fixed degree ) . the adopted strategy is to start with strictly homogeneous graphs ( networks in which identical connectivity patterns are `` seen '' from every node ) differing in the detailed loop structure ( fig . [ fig : hom4_ex ] ) . this leads us to an understanding of the role of these distinct local ordering properties , even for networks that are nolonger strictly homogeneous but which still obey @xmath6 , k being the constant degree .    in the case @xmath17 , examples of purely homogeneous networks are the random homogeneous network ( fig . [ fig : hom4_ex]a ) constructed according to the molloy reed algorithm @xcite , the square lattice ( fig . [ fig : hom4_ex]b ) , the kagom lattice ( usually used in condensed matter physics as this geometry represents one of the most frustrated antiferromagnetic systems , fig . [ fig : hom4_ex]c ) and a ring where nodes two units apart from each other are also directly connected ( fig . [ fig : hom4_ex]d ) . for the random network , the homogeneity definition given above only holds approximately . while loops of short length do not occur in the limit of large network size , long ones of various lengths may exist . therefore different nodes may not `` see '' identical connectiviy patterns . also the cayley tree lacks completely in loops . but as the last generation of nodes has degree 1 , it does not fall into the class of homogeneous networks . indeed , since the number of nodes belonging to generation @xmath18 is @xmath19 , most of the nodes are even comprised within the last generation . clearly , all the networks depicted in fig . [ fig : hom4_ex ] are characterized by the degree distribution @xmath20 , but they differ in the way the second neighbors are arranged , that is the loop structure becomes richer going from ( a ) to ( d ) . while the random network has no short loops at all , the ones of length 4 are a fingerprint of the square lattice . in the kagom lattice , every node is a corner of two distinct triangles . finally , in the ring of degree 4 , triangles and two different kinds of loops of length 4 are found , namely plaquettes ( as in the square lattice ) and the so - called primary quadrilaterals ( two adjacent triangles ) , see section ivb .    in order to study the epidemic process on disordered `` homogeneous '' topologies , we subject the networks in fig .  [ fig : hom4_ex]b - d to the following rearrangement algorithm ( fig . [ fig : rewir ] ) @xcite :    * choose randomly two links ( link 1 connecting node @xmath21 with @xmath22 and connection 2 linking vertex @xmath23 with @xmath24 ) that do not share a common node . * remove these two links and establish two new connections between @xmath21 and @xmath24 as well as @xmath23 and @xmath22 . repeating this rewiring procedure a certain number of times leads to locally varying numbers of loops while the degree distribution remains unaltered . in this section , we shall introduce the formalism that is used in order to describe the dynamics on top of the networks specified in the previous section . for completeness , we then derive the ( single - step ) mean - field approximation which could also be obtained heuristically . in the last paragraph , we derive a two - step description that allows us to gain insight about the loop structure in epidemic spreading .    on an exact level , we shall describe the epidemic dynamics by assigning a probabilty @xmath25 to each configuration @xmath26 at every instant of time @xmath27 . the vector @xmath26 contains the states @xmath28 of all the nodes @xmath29 of the network , @xmath28 being either 0 ( susceptible ) or 1 ( infected ) . the system probabilities satisfy at any time @xmath27 @xmath30 and evolve in time according to @xmath31 the transition matrix of the system @xmath32 is obtained from the matrix @xmath33 which shall denote the probability that the state of the arbitrary site @xmath34 changes from @xmath35 to @xmath36 , through @xmath37 n being the total number of nodes in the network . the matrix elements representing the probabilites for the possible events at the site @xmath34 are given by our version of the sis model , namely @xmath38 \\ w_{1 \rightarrow 1}^l & = 0 \qquad w_{0 \rightarrow 0}^l = 1-\lambda\bigl[1-\prod_{j \text{nn } l } ( 1-y_j)\bigr ] , \end{split}\ ] ] or in a more compact form @xmath39.\ ] ] the products in the above expressions have to be taken over all the nearest neighbors @xmath40 of node @xmath34 . the factor @xmath41 is 1 if at least one @xmath42 and 0 otherwise .    before deriving our two - timestep description , we show how the conventional mean - field approximation is retrieved through this formalism . at that level , the sites are considered independently from each other , and we write for the system probability @xmath43 i.e.  the system is described by the single variable @xmath44 [ the probability of being susceptible is @xmath45 . its dynamics is obtained from eq . ( [ equ : exact1 ] ) by summing it over all possible configurations @xmath26 , @xmath46 held fixed @xmath47 where the node @xmath48 can be chosen in an arbitrary way . the left hand side of the above equation corresponds to the probability that node @xmath48 is in state @xmath46 at time @xmath49 . with the ansatz ( [ equ : hom_mix ] ) , the time evolution is @xmath50.\ ] ] with @xmath51 and @xmath52 , by consequence , we obtain for small values of @xmath53 @xmath54 the stationary - state condition ( @xmath55 ) leads to a value for the epidemic threshold @xmath56 therefore all the networks depicted in fig . [ fig : hom4_ex ] are treated identically at this level of description . the only topological property that determines the critical value of the effective spreading rate @xmath57 is the number of nearest neighbors @xmath7 . improvements upon the mean - field description can be obtained by taking into account spatial correlations . based on the ordinary pair approximation @xcite , the presence of triangles enters by establishing a number of hypotheses about the open and closed triple correlations @xcite . this approach therefore embeds this topological property in a rather implicit way . in cluster approximations @xcite , the local topology and its associated spatial correlations translate directly into the choice of the cluster and a set of probabilities for its possible configurations . another strategy that serves to incorporate local ordering properties is to take into account temporal correlations . thus by performing two timesteps exactly , we expect that the way the second neighbors are arranged , enters very naturally into the description . for example , the cases where two nearest neighbors of an arbitrary node are also directly connected ( presence of a triangle ) , where they are linked via a second neighbor ( giving rise to a loop of length 4 ) or where the only path goes through the original node ( treelike structure ) , lead to different results . in the remaining part of this section , we derive the general equation , special cases are then looked at within the following section . as outlined above , we iterate eq . ( [ equ : exact1 ] ) once @xmath58.\ ] ] we now again pass to a site approximation . by summing the above equation over all possible configurations @xmath26 , @xmath46 held fixed , we get @xmath59 , \end{split}\ ] ]        0 again being an arbitrarily chosen node . furthermore , the system probability @xmath60 is given by eq . ( [ equ : hom_mix ] ) , and the nodes @xmath61 denote the nearest neighbors of the arbitrarily chosen node 0 . as only @xmath62-states associated to the vertex @xmath48 , its nearest and second neighbors appear in the @xmath63-factors , the sum over the @xmath62-variables associated to nodes more than 2 links away from vertex 0 , is carried out trivially . a tour de force calculation leads to    @xmath64 . \label{equ : exact2}\end{gathered}\ ] ]        thereby the connectivity embedding factor @xmath65 \\ & = \begin{cases } 1 & \text{if $ z_\alpha=0 $ and at least one $ z_\sigma=1$,}\\ 0 & \text{otherwise . } \end{cases } \end{split}\ ] ] and the expectation value of a function of the states of node 0 , its nearest and second neighbors , these vertices collectively being denoted by @xmath66 , @xmath67\ ] ] were introduced for notational convenience . in the following , an expectation value of a product of @xmath68 @xmath69-factors will be referred to as a term of @xmath68-th order although it is proportional to @xmath70 . as is illustrated in detail in the next section , every term of eq . ( [ equ : exact2 ] ) corresponds to a subgraph of the graph composed of the nodes @xmath66 and whose links are according to the network under investigation . it can already be anticipated that the first term accounts for the degree distribution only , whereas the contributions of higher order will give insight about the role of the loop structure . in this section , we elaborate the implications of our two - step description , the analysis being restricted to the stationary state .    in the left part of fig . [ fig : hom4_reg ] , we show the simulation results for the graphs introduced in section ii . the ring - type network exhibits the largest epidemic threshold . for both the square and kagom lattices , the critical value is @xmath71 . we therefore anticipate that either four plaquettes or two triangles ( per node ) lead to the same effect in the regime of low prevalences . finally , the lowest epidemic threshold is found if the population is arranged on a homogeneous random network ( of degree 4 ) . the last result is very intuitive since in such a graph , global distances are small , making it more easy for a virus to spread . therefore , even if the effective spreading rate is rather low , a finite fraction of the population will be infected in the stationary state , hence the small value for the location of the onset of the epidemic . in summary , these results indicate that the poorer the loop structure , the lower the corresponding epidemic threshold .    in the right part of fig . [ fig : hom4_reg ] , the one - step and two - step site approximations are reported . the former corresponds to the steady - state solution of eq . ( [ equ : mf_1step ] ) for which @xmath72 at @xmath73 according to eq . ( [ equ:1step_thr ] ) . all the networks in question are therefore treated identically , the loop structure being ignored at this level of description . yet , the two - step solutions [ eq .  ( [ equ:2steps ] ) ] are diverse for the different graphs . going from right to left , the curves correspond to the ring , the square lattice , the kagom lattice and the molloy - reed network , that is they appear in the same sequence as at the level of simulation . furthermore , the curves corresponding to the kagom and square lattice also meet the @xmath74-axis at the same value of @xmath57 . these findings confirm our intuitive arguments given in the previous section . it has to be noted that the two - step estimates for the threshold values are still considerably inaccurate especially for the ring and lattices , but this just highlights the presence of higher - order spatiotemporal correlations . however , the important point is that the degeneracy associated to the one - step description disappers at the two - step level . we shall now evaluate all the terms of eq . ( [ equ : exact2 ] ) for a locally treelike topology . [ fig : tree_sub ] shows the subgraphs representing the terms in eq . ( [ equ : exact2 ] ) , in increasing order . thereby the correspondance is as follows : given the term @xmath75 , the nodes @xmath76 and @xmath77 are represented by filled circles whereas their nearest neighbors are drawn by empty circles . the links which enter at the level of the subgraph in question , are represented by solid lines whereas the ignored ones are dashed . if we denote the second neighbors of the central vertex 0 by @xmath78 for @xmath79 and follow eq . ( [ equ : expect ] ) , the first order contribution for @xmath80 ( subgraph in fig . [ fig : tree_sub]a ) is @xmath81}_{f_1}\bigr\ } \\ & \qquad = 4p + \mathcal o(p^2 ) \end{split}\ ] ] where the sum over the @xmath62-variables to which no circles are associated , has been carried out trivially . furthermore , we again have set @xmath82 in the third line ( this is also done below ) , and the time index was omitted since we are only interested in the steady state . this term appears with multiplicity 4 ( due to @xmath83 ) , giving the contribution @xmath84 to first order in @xmath85 . [ fig : tree_sub]b shows the subgraphs representing @xmath86 ( upper part ) and @xmath87 ( lower part ) . their contributions are @xmath88 occuring @xmath89 times and @xmath90 thus not giving a contribution to first order in @xmath85 . as a consequence , the second - order term ( that is the one proportional to @xmath91 ) is @xmath92 . as the procedure should now be clear , we only give the results for the remaining orders . the upper subgraph of fig . [ fig : tree_sub]c represents the term @xmath93 . its contribution is @xmath94 . the lower subgraph corresponds to @xmath95 , yielding @xmath96 . as the former term has multiplicity @xmath97 , the total third - order contribution is @xmath98 .    as far as the fourth order is concerned ( fig . [ fig : tree_sub]d ) , the subgraph involving node 0 as filled circle neither gives a contribution whereas the term @xmath99 having multiplicity 1 also gives @xmath94 , thus totally yielding @xmath100 .    collecting these findings , we obtain the following condition that determines the epidemic threshold for a treelike topology @xmath101 which is satisfied by @xmath102 . this is the value that can be extracted from the right panel of fig . [ fig : hom4_reg ] ( second curve from the left ) . it is easy to imagine that the preceeding analysis yields different results when triangles and loops of length 4 are present . caldarelli et al .  @xcite have classified loops of length 4 in a complex network into _ primary _ and _ secondary quadrilaterals _ [ fig : quadlat ] ) . in the former case , the external vertices which the loop is composed of are all nearest neighbors whereas secondary quadrilaterals are plaquettes , the external nodes being two nearest and a second neighbor . with these concepts , the loop structure of a strictly homogeneous graph can quantitatively be characterized as follows : by choosing an arbitrary node , the number of edges between its nearest neighbors is denoted by @xmath103 . @xmath104 and @xmath105 shall refer to the number of primary and secondary quadrilaterals . for the networks depicted in fig . [ fig : hom4_ex]b - d , we report the corresponding values in tab .  [ tab : simplenets_eq ] .
Abstract:

Answer: the spreading of an epidemic is determined by the connectiviy patterns which underlie the population . while it has been noted that a virus spreads more easily on a network in which global distances are small , it remains a great challenge to find approaches that unravel the precise role of local interconnectedness . such topological properties enter very naturally in the framework of our two - timestep description , also providing a novel approach to tract a probabilistic system . the method is elaborated for sis - type epidemic processes , leading to a quantitative interpretation of the role of loops up to length 4 in the onset of an epidemic .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: with the rapid growth in demand for higher data rates , wireless technology is shifting towards the systems with multiple antennas . multiple - input multiple - output ( mimo ) systems in wireless communications provide significant improvements in wireless link reliability and the achievable capacity @xcite . recently , the systems with large number of antennas ( tens to hundreds of antennas ) called as _ large - mimo _ systems are getting increased attention @xcite because they provide additional multiplexing and diversity gains . using spatial multiplexing in mimo systems , multiple data streams can be transmitted simultaneously from different transmit antennas . however , the detection of these data streams at the receiver end is challenging . to achieve minimum bit error rate ( ber ) performance in mimo systems , an exhaustive search over all the possible transmit vectors is required . this is called as the maximum likelihood ( ml ) search for symbol vector detection in mimo systems . however , as the number of antennas grow , the number of possible transmit vectors also increases exponentially and thus the ml search becomes computationally impractical . sphere decoder , a well known mimo detector achieves near ml performance but is practical only up to limited number of dimensions @xcite . + the traditional mimo detection techniques involve linear detection such as zero forcing ( zf ) detector and minimum mean squared error ( mmse ) detector . these detectors use linear transformation of the received vector in order to estimate the transmitted symbol vector . comparatively the mmse detector is superior over zf detector in terms of ber performance , but still their performance is far inferior compared to the ml performance . another type of detection involve non - linear detectors such as vertical bell labs layered space time architecture ( v - blast ) @xcite which utilizes ordered sic for symbol vector detection . some of the other mimo detection technique involve ant colony optimization based mimo detection @xcite , message passing based algorithm @xcite , channel hardening based algorithm @xcite , and lattice reduction aided mimo detection @xcite . as an alternative , in this study we employ sic based detection technique for mimo detection problem . + sic based mimo detection is a well known detection technique where symbols are detected sequentially @xcite . in sic , after detecting each symbol , its interference is canceled from the received vector in order to improve the instantaneous signal to interference plus noise ratio ( sinr ) for the remaining symbols . it is also known as layered detection where in each layer , one symbol is detected . however , it suffers from error propagation which occur due to wrong decisions in early stages of the algorithm @xcite . for improving the performance of sic , a well known technique known as v - blast @xcite technique is used which performs the following steps : ( 1 ) snr based ordering in the detection sequence , ( 2 ) using zf or mmse for nulling the interference among the data streams and ( 3 ) detecting a symbol and canceling its interference i.e. sic . further , to mitigate the effect of error propagation , one of the available algorithms in the literature include multi - branch sic ( mb - sic ) algorithm @xcite where the concept of multiple branch ( mb ) processing is utilized . each branch of mb - sic differs in ordering of the detection sequence thus resulting in a higher detection diversity over the conventional sic . recently , a multiple feedback ( mf ) strategy is proposed in @xcite for sic based mimo detection . multiple feedback sic ( mf - sic ) algorithm proposed in @xcite is based on the concept of shadow are constraint ( sac ) where if the decision fall in the shadow region ( unreliability region ) then multiple constellation points are used in the decision feedback loop . the best symbol from the multiple symbols used in the decision feedback is selected using the ml criteria . + in this article , we propose an improved multiple feedback successive interference cancellation ( imf - sic ) for symbol vector detection in mimo spatial multiplexing systems . in mf - sic , once a decision falls into the shadow region , conventional sic is used to find the best one symbol for the corresponding layer . however , it may happen that while computing the best symbol for a given layer , multiple decisions from the subsequent layers fall into the shadow region . this condition is not checked in mf - sic which sometimes results in error propagation and limits the ber performance . in improved mf - sic ( imf - sic ) , we overcome this limitation by checking the shadow region criteria recursively rather than using the conventional sic for searching the best candidate symbol . this results in significant reduction in error propagation and thus the ber performance can be improved . also , we propose an ordered imf - sic ( oimf - sic ) algorithm where we employ the log likelihood ratio ( llr ) based dynamic ordering in the detection sequence @xcite . by the term dynamic ordering we mean that after every successful decision about a symbol , the ordering of the detection sequence is updated based on new llr values . in oimf - sic , same ordering pattern is followed even for searching the best symbol when a decision is found unreliable . simulation results show that the proposed algorithms significantly outperform the conventional sic and the mf - sic based detection methods , and that they achieve a near ml performance . + rest of the article is organized as follows . in sect . [ sec2 ] , we present mathematical model of mimo system . an overview of some of the traditional mimo detection methods is given in sect . the proposed algorithms for symbol vector detection in mimo spatial multiplexing systems is discussed in sect . [ sec4 ] . in sect . [ sec5 ] , we show the simulation results on comparison of ber performance . finally , in sect . [ sec6 ] we conclude the article . in this section , we discuss the mathematical model of mimo system as shown in fig.[figs1 ] . we consider a point to point spatially multiplexed mimo link with @xmath0 transmit antennas and @xmath1 receive antennas . in spatial multiplexing , each antenna transmits a different information symbol @xmath2 for @xmath3 . let @xmath4^t$ ] be the @xmath5 dimensional transmit vector , @xmath6^t$ ] denote the transpose of a matrix . each entry @xmath2 of the transmit vector @xmath7 is taken from a signal constellation @xmath8 , for example @xmath9 for 4-qam signaling . the channel over which the symbol vector @xmath7 is transmitted is assumed to be a rayleigh distributed flat fading channel . let us consider an @xmath10 complex channel matrix @xmath11 with its elements @xmath12 for @xmath13 and @xmath14 assumed to be independent and identically distributed ( i.i.d . ) as complex normal with mean 0 and variance 1 i.e. @xmath15 . the received symbol vector @xmath16 , after demodulation and matched filtering can be written as @xmath17 where @xmath18 is the additive white gaussian noise ( awgn ) with its element @xmath19 for @xmath13 as i.i.d . the signal to noise ratio ( snr ) is defined as @xmath21 , where @xmath22 is the average energy per symbol defined as @xmath23)$ ] and @xmath24 is the element - wise noise variance . tr(@xmath25 ) denote the trace of a matrix , @xmath26 $ ] denote the expectation operation and @xmath27 denote the hermitian of a matrix . the channel state information ( csi ) is assumed to be perfectly known at the receiver but is unknown at the transmitter . the ml solution for a given received vector @xmath16 and a known channel matrix @xmath11 can be found as @xmath28 ml detection ( mld ) requires an exhaustive search over @xmath29 possible transmit symbol vectors where @xmath30 is cardinality of the constellation set @xmath8 . this search grows exponentially with increase in the modulation order ( @xmath30 ) or increase in the number of transmit antennas ( @xmath0 ) or both . in this section , we discuss some of the traditional mimo detection techniques such as linear detection techniques which involve zero forcing ( zf ) based mimo detection and minimum mean squared error ( mmse ) based mimo detection , and nonlinear detection techniques such as sic based mimo detection .      in zf based mimo detection , a linear transformation matrix @xmath31 is used in such a way that it removes the effect due to the channel impairments on the received vector @xmath16 i.e. @xmath32 which is the pseudo inverse of channel matrix @xmath11 . the linear transformation of the received vector using @xmath31 can be written as @xmath33 thus the zf estimate of the transmitted symbol can be found as @xmath34,\ ] ] where @xmath35 $ ] is the quantization operation which maps the soft values on to the nearest constellation point . the computational complexity of zf receiver is very less , however the ber performance is far inferior to the ml performance . the main drawback of zf based mimo detection is the problem of noise enhancement which degrades its performance . mmse based mimo detection is based on finding a suitable linear transformation matrix which minimizes the mean squared error between the transmit vector @xmath36 and the transformed received vector as @xmath37.\ ] ] the solution to eq . [ eq5 ] can be written as @xmath38 using the transformation matrix @xmath39 , the mmse solution can then be computed as @xmath40.\ ] ] mmse based mimo detection overcomes the problem of noise enhancement in zf detectors . thus it achieves better ber performance over zf but still the performance of mmse receiver is inferior as compared to the ml performance . one more drawback is that the mmse solution needs knowledge of the noise variance @xmath24 .      in sic based mimo detection , the symbols transmitted from different transmit antennas are detected in a sequential manner . after detecting a symbol for a particular transmit antenna , its interference is canceled from the received vector in order to reduce the interference between different data streams at the receiver . the received vector in eq . [ eq1 ] can be rewritten as @xmath41 using the knowledge of channel matrix @xmath11 a filter @xmath42 is defined as @xmath43 where @xmath44 is the matrix with its column taken from the @xmath45 to @xmath0 columns of the channel matrix @xmath11 . we use vector @xmath46 as the vector left after removing the interference due to @xmath47 as @xmath48 using @xmath42 , the first symbol @xmath2 is estimated as @xmath49=\mathcal{q}[\mathbf{w}_i^h\mathbf{y}_i],\ ] ] where @xmath50 is estimate of the transmitted symbol @xmath2 . @xmath42 , @xmath51 and @xmath46 are updated after every successful decision about a symbol . sic is an iterative detection technique , where in every iteration one symbol is detected . thus a total of @xmath0 iterations are performed to detect the complete symbol vector . a pseudo code of sic based mimo detection is given in algorithm [ algo1 ] . ; ; compute @xmath52 @xmath53 @xmath54 $ ] update @xmath55 update @xmath56 $ ]    however , sic suffers from the problem of error propagation which occurs due to the wrong decisions in early stage of the algorithm . once an error occurs , it propagates to the later stages and thus increases the number of errors . this degrades the ber performance of sic based mimo detection . in this section , first we discuss the mf - sic based algorithm for symbol vector detection in mimo systems . we will then discuss the imf - sic algorithm and oimf - sic algorithm . the conventional sic suffers from the problem of error propagation which degrades the ber performance . as a solution to mitigate the effect of error propagation , mf strategy has been proposed in @xcite for sic based mimo detection . the mf strategy is based on the concept of reliability of a soft decision in sic . let us consider the @xmath57 soft decision within sic as @xmath58 the constellation point near to the soft decision @xmath59 is nothing but its quantized value i.e. @xmath60 $ ] . let us define the distance between the soft decision and its quantized value as @xmath61\vert\ ] ] if the soft decision in sic is within a predefined threshold radius @xmath62 ( i.e. if @xmath63 ) around the nearest constellation point then the decision is said to be reliable and its quantized value is used in the decision feedback . but if the decision is not within the radius ( if @xmath64 ) then multiple neighboring ( say @xmath65 ) constellation points are used instead of the quantized value followed by the conventional sic . it can be thought of as running parallel streams one for each symbol and at last the best one is selected using the ml metric ( eq . [ eq2 ] ) . mf - sic performs superior than the conventional sic in terms of ber . a pseudo code of mf - sic is shown in algorithm [ algo2 ]    ; ; compute @xmath52 @xmath53 @xmath66\vert$ ] @xmath54 $ ] @xmath67 i.e. initialize by @xmath68 neighborhood of @xmath69 @xmath70 @xmath71 $ ] update @xmath72 update @xmath73 $ ] compute the solution vector @xmath74 @xmath75 @xmath76 update @xmath55 update @xmath56 $ ]      in imf - sic , we propose an improved mf strategy for sic based symbol vector detection in mimo spatial multiplexing systems . in mf - sic algorithm @xcite , multiple feedback loops are followed by the conventional sic to check for the optimal decision for a particular layer . however , while searching for the optimal decision , the shadow region ( unreliability ) criteria is not considered in the lower layers and hence there are chances of selecting a sub - optimal decision about a symbol which degrades the ber performance . in order to improve the accuracy about a decision among multiple symbols in the feedback loop , we check the shadow region criteria recursively in imf - sic algorithm . we therefore formulate a separate sub - routine ( function ) named as imf_sic in algorithm [ algo3 ] which checks the unreliability condition and optimally decide the symbol . sub - routine imf_sic is called from the main function which is given in algorithm [ algo4 ] and also from the sub - routine itself . we denote the number of calls of the sub - routine by @xmath77 in imf - sic algorithm which is nothing but the number of times unreliability condition is checked once a decision falls into the shadow region . the number of neighboring symbol which are used in the decision feedback is denoted by @xmath65 . +    ; ; compute @xmath52 @xmath53 @xmath66\vert$ ] @xmath54 $ ] @xmath78 imf_sic(@xmath79 ) update @xmath55 update @xmath56 $ ]    function [ @xmath50]=imf_sic(@xmath79 ) @xmath67 i.e. initialize by @xmath68 neighborhood of @xmath69 @xmath54 $ ] @xmath78 imf_sic(@xmath81 ) update @xmath72 update @xmath73 $ ] compute the solution vector @xmath74 @xmath75 @xmath76 return @xmath50    significant reduction in the effect of error propagation is noticed using imf - sic ( discussed in sect . [ sec5 ] ) which results in improvement in the ber performance over sic and the mf - sic . however , imf - sic does not include ordering in the detection sequence which plays a very important role in mitigating error propagation in sic based detection techniques . in the next section , we discuss ordered imf - sic ( oimf - sic ) algorithm for mimo detection .      in this section , we propose an ordered imf - sic ( oimf - sic ) algorithm which incorporates the llr based ordering in the detection sequence . error propagation can be mitigated successfully if the decisions in the early layers of sic are detected correctly . in several sic based detection algorithms , channel norm ordering is used such as in v - blast technique @xcite . however , channel norm ordering does not consider the instantaneous noise in the received vector which has a significant impact on the ber performance . thus a log likelihood ratio based ordering for multi - user detection in cdma systems @xcite has been proposed in @xcite for mf - sic based mimo detection to overcome the sub - optimal behavior of channel norm based ordering . we employ llr based dynamic ordering in the imf - sic algorithm to order the detection sequence . the term dynamic ordering is used in the sense that the ordering is updated after every successful symbol decision . however , the modification in the proposed method is that the ordering is also used in the oimf_sic sub - routine which help in a more precise decision . the llr value of a soft estimate @xmath82 can be written as @xcite @xmath83 where @xmath84 . the value of the argument @xmath85 which maximizes @xmath86 is arranged in descending order . let us consider a set @xmath87 be the set of indices of the symbols which are not detected so far . the set @xmath87 is updated after every successful decision . for ordering the detection sequence in every iteration we consider the indices which are in the set @xmath87 i.e. in eq . [ eq14 ] @xmath88 . let us also consider the ordered sequence @xmath89 which is updated in every iteration . in algorithm [ algo5 ] , we present the oimf - sic algorithm for mimo detection .    ; compute @xmath90 @xmath91 compute @xmath92 using the equation [ eq14 ] arrange the values @xmath93 in descending order update the detection order @xmath89 start detection using imf - sic given in algorithm [ algo3]with ordered sequence @xmath89 . in this section , we discuss the ber performance plots of the proposed algorithms and compare with the performance of the conventional sic and the mf - sic algorithm for @xmath94 , @xmath95 and @xmath96 mimo systems with 4-qam and 16-qam signaling . the results are generated using matlab and the ber is averaged over @xmath97 monte - carlo simulations . + in fig . [ figs2 ] , we compare the ber performance of the conventional sic and mf - sic with the proposed algorithms for @xmath94 mimo system with 4-qam . the value of threshold parameter @xmath62 used is 0.2 and the number of recursions @xmath98 . observation reveals that the proposed algorithms achieve a near ml performance and also the performance of imf - sic and oimf - sic algorithms is similar in terms of ber performance . at ber of @xmath99 , imf - sic results in approximately 1 db gain in snr over mf - sic . in fig . [ figs3 ] , the ber performance is compared for @xmath95 mimo system with 4-qam . a significant improvement in the ber performance of imf - sic and oimf - sic can be seen over sic and mf - sic based mimo detection techniques . an snr gain of around 2 db can be achieved by using imf - sic over mf - sic to achieve a target ber @xmath100 . due to near optimal performance of imf - sic for @xmath94 and @xmath95 mimo systems with 4-qam the advantage of dynamic ordering can not be observed . however , with increase in number of antennas or the modulation order , ordering plays an important role in mitigating the error propagation . [ figs4 ] presents the ber performance comparison for @xmath96 v - blast mimo system with 4-qam . a significant improvement in the performance of oimf - sic can be seen when compared with imf - sic for @xmath98 . however , to achieve the near ml performance the number of recursions are increased to @xmath101 and @xmath102 for oimf - sic . thus , the effect of increase in number of recursions on the ber performance can be observed . +     v - blast mimo system with 4-qam and the value used for @xmath103 , @xmath104 and @xmath105.,width=9,height=10 ]     v - blast mimo system with 4-qam and the value used for @xmath103 , @xmath104 and @xmath105.,width=9,height=10 ]     v - blast mimo system with 4-qam.,width=9,height=10 ]    the ber performance for @xmath94 and @xmath95 v - blast mimo systems with 16-qam signaling is shown in figs . [ figs5 ] and [ figs6 ] . in 16-qam the total number of constellation points is 16 and hence in order to get a reliable decision , we keep the number of neighboring contellation point in the feedback as @xmath106 . the threshold value @xmath62 is kept same i.e. @xmath103 . it can be observed that the ber performance of oimf - sic is close to within the 0.1 db of the ml performance . ordering in imf - sic results in almost 1 db gain in snr over the imf - sic without ordering . however , in fig . [ figs6 ] , to achieve the ber close to ml , the number of recursions is increased to @xmath101 in oimf - sic . this shows a snr gain of more than 2 db over imf - sic at a target ber @xmath99 . v - blast mimo system with 16-qam and the value used for @xmath103 , @xmath104 and @xmath107.,width=9,height=10 ]     v - blast mimo system with 16-qam.,width=9,height=10 ] we proposed an improved multiple feedback successive interference cancellation ( imf - sic ) algorithm for symbol vector detection in mimo systems . we also proposed an ordered imf - sic ( oimf - sic ) algorithm which employs log likelihood ratio ( llr ) based dynamic ordering in the detection sequence . the ordering of the detection sequence is updated after every successful decision and thus the effect of error propagation is mitigated . the proposed algorithms significantly improves the ber performance and outperform the conventional sic algorithm and the mf - sic algorithm , and achieve near ml performance . the authors would like to thank iit indore and ministry of human resource development ( mhrd ) india for their support . p.  w.  wolniansky , g.  j.  foschini , g.  d.  golden , and r.  a.  valenzuela , `` v - blast : an architecture for realizing very high data rates over the rich - scattering wireless channel , '' _ international symposium on signals , systems , and electronics _ 295 - 300 , 1998 . t.  lakshmi  narasimhan and a.  chockalingam , `` channel hardening - exploiting message passing ( chemp ) receiver in large - scale mimo systems , '' _ ieee jrnl . on sel . topics in signal proc . : special issue on signal processing for large - scale mimo communications _ , vol . 847 - 860 , 2014 . k.  a.  singhal , t.  datta , and a.  chockalingam , `` lattice reduction aided detection in large - mimo systems , '' _ ieee 14th workshop on signal processing advances in wireless commun . ( spawc ) _ , pp . 594 - 598 , 2013 . g.  golden , c.  foschini , r.  valenzuela , and p.  wolniansky , `` detection algorithm and initial laboratory results using v - blast space - time communication architecture , '' _ electronics letters _ 1 , pp . 14 , 1999 . r.  fa and r.  c.  de  lamare , `` multi - branch successive interference cancellation for mimo spatial multiplexing systems : design , analysis and adaptive implementation , '' _ iet commun . _ , vol . 5 , no . 484 - 494 , 2009 . p.  li , r.  c.  de  lamare , and r.  fa , `` multiple feedback successive interference cancellation detection for multiuser mimo systems , '' _ ieee trans . on wirel . 10 , no . 8 , pp . 2434 - 2439 , 2011 . l.  yang , `` receiver multiuser diversity aided multi - stage mmse multiuser detection for ds - cdma and sdma systems employing i - q modulation , '' in _ veh . fall ( vtc 2010-fall ) , 2010 ieee 72nd _ ieee , pp . 1 - 5 , 2010 .    t.  hesketh , p.  li , r.  c.  de  lamare , and s.  wales , `` multi - feedback successive interference cancellation with dynamic lo - likelihood - ratio based reliability ordering , '' in _ proc . of the tenth international symposium on wirel . systems ( iswcs 2013 ) _ ieee , pp . 1 - 5 , 2013 .
Abstract:

Answer: in this article , we propose an improved multiple feedback successive interference cancellation ( imf - sic ) algorithm for symbol vector detection in multiple - input multiple - output ( mimo ) spatial multiplexing systems . the multiple feedback ( mf ) strategy in successive interference cancellation ( sic ) is based on the concept of shadow area constraint ( sac ) where , if the decision falls in the shadow region multiple neighboring constellation points will be used in the decision feedback loop followed by the conventional sic . the best candidate symbol from multiple neighboring symbols is selected using the maximum likelihood ( ml ) criteria . however , while deciding the best symbol from multiple neighboring symbols , the sac condition may occur in subsequent layers which results in inaccurate decision . in order to overcome this limitation , in the proposed algorithm , sac criteria is checked recursively for each layer . this results in successful mitigation of error propagation thus significantly improving the bit error rate ( ber ) performance . further , we also propose an ordered imf - sic ( oimf - sic ) where we use log likelihood ratio ( llr ) based dynamic ordering of the detection sequence . in oimf - sic , we use the term dynamic ordering in the sense that the detection order is updated after every successful decision . simulation results show that the proposed algorithms outperform the existing detectors such as conventional sic and mf - sic in terms of ber , and achieves a near ml performance . mimo detection , spatial multiplexing , decision feedback , successive interference cancellation , bit error rate .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: many high energy experiments require pure electron beams . despite the steady improvement of the beam lines , contamination below a level of few % is very difficult to achieve . an example is the na64 experiment at cern in which it is mandatory to suppress hadron and muon contamination in the electron beam since such particles can generate irreducible background processes mimicking the experimental signature of a dark photon @xcite . na64 uses 100 gev electrons from the h4 sps beam line at cern which is one of the best existing beam lines at this energy in terms of beam purity @xcite . the standard technique of detecting cherenkov radiation to distinguish between electrons and heavy charged particles is very inefficient at energies of 100 gev @xcite . instead in na64 we detect the synchrotron radiation , produced by the electrons passing through a dipole bending magnet , using bgo crystals located downstream the magnet . na64 is a fixed - target experiment at the cern sps @xcite combining the active beam dump and missing energy techniques to search for rare events predicted in dark sector models @xcite . na64 uses 100 gev electrons from the h4 sps beam line at cern . the electron beam is dumped in an active target , an electromagnetic calorimeter ( ecal ) made of lead and scintillators in a sandwich geometry of shashlik type ( corresponding to about 40 radiation lengths ) @xcite . if dark sectors exist , bremsstrahlung photons generated in the target could produce dark photons ( @xmath3 ) via kinetic mixing @xcite . those could then decay in a pair of particles in the dark sector ( @xmath4 ) which will escape the setup undetected . the basic concept of the synchrotron radiation emitted from high energy electrons to reject heavy charged particles is to exploit the high suppression of the radiated power emitted by particles heavier than electrons passing through a magnetic field in order to discriminate them . the use of this technique is not new and detection of electrons or positrons in electrons beams with momenta ranging from 30 to 50 gev was reported earlier by @xcite and the achieved suppression factor was around 1% . however , here we present the results obtained with this method using the bgo detector employed in na64 . as will be shown , our results have a very good overall efficiency 95% and 3 orders of magnitude better rejection factor than previously reported using a similar scheme @xcite .    a charged particle in a magnetic field moves in a circular motion emitting photons along its trajectory due to the basic principles of electrodynamics . both quantum and classical theory of synchrotron radiation ( sr ) are well understood @xcite . in the range of interest for our experiment both treatments are equivalent and we can therefore use the classical approximation for our calculations . the total power @xmath5 emitted per unit length by a relativistic charged particle of energy @xmath6 with mass @xmath7 and with bending radius @xmath8 in a magnetic field @xmath9 perpendicular to its velocity is given by : @xmath10 where @xmath11 is the charge of the particle and @xmath12 the speed of light , since the emission angle of the synchrotron photons is proportional to the inverse of the lorentz factor @xmath13 , the photons are emitted tangentially to the particle trajectory . the total emitted power has an inverse scaling to the fourth power of the charged particle mass . therefore , heavy charged particles emit orders of magnitude less synchrotron radiation than light ones . for @xmath14 and @xmath15 with about 200 @xmath16 mass , one can estimate that they radiate @xmath17 times less than an electron . this would be the case if the particles propagate in an ideal vacuum . however , in a real experimental setup , vacuum windows , residual gas , beam counters such as scintillators and trackers result in interactions of the incoming particles with material . therefore , the suppression factor when crossing materials is limited by the emission of secondary electrons with enough kinetic energy ( several mev ) to leave a synchrotron - like signal in the detector . although most of the energy transfer due to ionization for heavy charged particles is only a few kev , rare high energy transfer is possible . the distribution of such secondary electrons with kinetic energy @xmath18 , where @xmath19 is the mean excitation energy of the atom / molecule , for a particle with velocity @xmath20 and charge @xmath21 passing through a material with atomic number @xmath22 , mass number @xmath23 and thickness @xmath24 is described by @xcite : @xmath25 the constant @xmath26 is defined as @xmath27 where n@xmath28 is the avogadro s number , @xmath29 is the classical electron radius and @xmath30 the electron mass . @xmath31 is a spin - dependent factor , which in our case for @xmath32 is very close to unity . w@xmath33 is the maximal energy transfer in a single collision to the electron : @xmath34 for a @xmath0 at 100 gev , @xmath35 is roughly 1 gev which covers completely the energy range where synchrotron radiation is emitted . eq.[eqn : knock - on ] is valid in the range @xmath36 .    the total emitted power scales with the energy of the charged particle to the fourth power . in principle , this scaling would make it possible to correlate the total energy emitted via synchrotron radiation with the momentum of the particle . however , the broad total energy distribution of the synchrotron radiation would result in a very low efficiency for rejecting 50 gev electrons in a 100 gev beam . furthermore , experimental issues limit this approach . in fact , electrons interacting with the material or the residual gas can generate bremsstrahlung photons that overlap with the synchrotron radiation spectrum and spoil the energy reconstruction capability . for this reason , in na64 such an approach alone is not enough to achieve the sensitivity required by the experiment and therefore a tracking system is mandatory to reject low energy electrons @xcite . the na64 experimental setup @xcite used for this measurement is shown in fig.[fig : na64setup ] . as shown in fig . [ fig : na64setup ] , a set of veto counters ( denoted as v2 in the drawing ) and 4 hadronic calorimeters ( hcal ) modules placed directly after the ecal ensures the hermeticity required to search for the missing momentum with the aimed sensitivity . the experimental signature is thus an energy deposition below a given threshold ( around 50 gev in na64 ) in the ecal , no activity in v2 and no energy deposition in the hcal . in fig . [ fig : ecal_energy ] the energy deposition in the ecal for a 100 gev electron beam is shown . in addition to the expected peak at 100 gev from electrons , one can see the presence of hadrons punching through the ecal and the peak close to zero from minimum ionising particles ( mip ) . since hadrons and muons are a source of background for the experiment @xcite , their contamination in the beam has to be suppressed at a level of @xmath1 . the synchrotron radiation detector ( srd ) described in this paper was designed to serve this purpose .            to reject low energy electrons and heavier charged particles , a synchrotron radiation detector ( srd ) requires good energy resolution ( see eq . [ eqn : power - emitted ] ) . moreover , in order to work efficiently at high rates of @xmath37 particles / s good timing properties are necessary . we chose bi@xmath38ge@xmath39o@xmath40 ( bgo ) crystals the properties of which are a good compromise to fulfil these requirements . the detector consists of 8 hexagonal crystals with an external diameter of 55 mm and a length of 200 mm . the crystals are grouped into two modules . each crystal is wrapped in teflon tape for efficient light collection and it is glued to an etl 9954 photomultipliers ( pmt ) . the gain of the pmts was set using a 100 gev @xmath0 beam impinging directly on the srd . the pions crossing the bgo have a mean energy deposition of about 60 mev as expected for a mip particle which is in the range of the sr radiation . the bgo has a density of 7.13 g/@xmath41 and because of the high atomic number of the bismuth component ( z=83 ) it has one of the largest probability per unit volume for photoelectric absorption of gamma rays @xcite . the light yield of about 8500 photons / mev coupled to the transportation losses and quantum efficiency of the pmt gives an energy resolution of about 17 % ( fwhm ) at 1.27 mev ( measured with a @xmath42na radioactive source ) . all the crystals were tested using cosmic radiation and their decay time was measured to be about 300 ns in good agreement with the value quoted in the literature @xcite . the bgo s signals are digitized by the msadc system @xcite which is connected to the daq of the experiment . the daq is downscaled version of compass ifdaq system @xcite . simulation of the expected sr signal was performed with the geant 4 package@xcite . the geometry of the na64 experiment was coded in geant 4 , including the 200 @xmath43 m mylar vacuum windows , the detailed composition of the trackers , scintillators and the residual gas was set at a level of @xmath44 mbar as in the measurements . the bgo saturation was taken into account using birks law with the constants taken from @xcite . the expected sr spectra for pions and electrons with energies of 50 gev and 100 gev are shown in fig . [ srspectrum ] . the plot shows the expected dependence on the incoming electron energy in the emission spectra for the realistic experimental conditions . moreover , the comparison between the sr spectra of pions and electrons illustrates clearly the principle of this technique that allows to discriminate between them by requiring an energy threshold in the synchrotron detector . for pions , one can see that the probability of detecting an event with energy above 1 mev ( the threshold in our detector ) is about @xmath45 . these sr - like signals originate from the interactions of the incoming pions with material which they ionise as predicted by eq . [ eqn : knock - on ] . furthermore , geant 4 reproduces the critical energy @xmath46 which divides the spectrum into two parts of equal power is : @xmath47 with the reduced plank constant @xmath48 and the bending radius @xmath8 . for 100 gev electrons in the @xmath49 t bending field this corresponds to @xmath5011.35 mev . the expected mean energy of a synchrotron photon @xmath51 mev is in very good agreement with simulation . the number of photons emitted per revolution in this energy range in the field of 7 t@xmath52 m is defined as : @xmath53 where @xmath54 is the fine structure constant . by scaling this equation for the fraction of the circle where the particles are inside the magnetic field , one obtains a mean number of emitted photon of about 24 . the srd geometrical acceptance is about one third , thus one can estimate that the sum of deposited energy is approximately 29.35 mev in good agreement with the results of the simulation as shown in fig . [ srspectrum ] . ( black dashed / solid line ) and 50/100 gev @xmath0 ( red dashed / solid line ) . ] the srd detector was tested during the na64 test beam run in july 2016 . the two bgo rows are parallel to the primary beam direction as shown in fig.[fig : newgeo ] . the dipole magnets installed in series ( see fig . [ fig : na64setup ] ) produce a total integrated magnetic field of 7 t@xmath52 m @xcite resulting in a nominal deviation for the incoming electrons at the srd / ecal positions of 31/34 cm from the undeflected beam axis . the two rows of bgo have about 9 cm distance from both the undeflected and the deflected beam axis ( fig.[fig : newgeo ] ) . this separation minimises the possibility for bremsstrahlung photons and neutrals produced in upstream beam particle s interactions with material ( vacuum windows , collimators , trackers ... ) and for particles in the beam halo ( the typical beam diameter used in the experiment was 2 cm ) to hit the srd . in fact , such interactions result in the saturation of the srd with a significant loss of efficiency due to the long decay time of the bgos .    the two crystals facing the beam ( labeled 3 and 7 in fig . [ fig : newgeo ] ) detect most of the energy emitted by the synchrotron radiation . we will refer to those as srd bgos from now on . the remaining six crystals are used to detect events with high energy deposition in the srd . in particular the last two crystals of each row ( labeled 0 and 4 in fig . [ fig : newgeo ] ) detect some energy only in the case of very energetic bremsstrahlung events and thus can be used as a veto ( see fig.[fig : newgeo ] ) . the six crystals after the srd bgos act also as a shield from backscattering particles coming from the ecal suppressing pions by an additional order of magnitude . finally in this new geometry it is possible to use the coincidence of the two srd bgo crystals to improve the tagging of synchrotron photons by rejecting knock - on electrons produced by incoming pions . in fact synchrotron radiation has a homogenous spectrum in the whole arc described by the primary and deflected beam and thus a signal is detected in both srd bgos . on the contrary , electrons generated by a @xmath0 undergoing hadron ionisation will mostly leave energy only in a single crystal as illustrated in fig . [ fig : newgeo ] . with the requirement of detecting in both srd bgos an energy deposition above a 1 mev the suppression factor is improved up to a level of @xmath1 . data with a 100 gev @xmath0 beam were taken to have a direct measurement of the suppression factor achievable for hadrons through synchrotron radiation measurements . the beam intensity was 5.3@xmath55 particles per spill . the trigger was given by the coincidence of the three plastic scintillator counters ( s1 , s2 and s3 shown in fig.[fig : na64setup ] ) . the additional requirement of an energy deposition below 60 gev in the ecal was applied in order to select an almost pure @xmath0 sample of @xmath56 collected events . the probability for electrons to punch through the 40 radiation length of the ecal was estimated to be at a level of @xmath57 @xcite . for the 100 gev electron beam run , a total of 220 spills were recorded with an intensity of 3.4@xmath58/spill . the same trigger used in the pion run was used for the electron data . in this case though , in order to reduce the pion contamination which is at a level of few % and obtain a pure sample of electrons only events with a total energy deposition in ecal + hcal above 90 gev but with less than 20 gev energy in the hcal were used . the energy spectra recorded by the srd bgo with electrons and pions are shown in fig.[fig : comp_spectra ] . the sr spectra obtained with the electron beam are used to perform the bgo calibration by comparison with the simulation . with this method a very good agreement of data and mc is achieved ( see fig.[fig : comp_spectra ] plot on the left ) . as a cross check using the obtained calibration constants , the data from the pion beam impinging directly on the srd are fit with a landau distribution . the obtained peak position of 60 mev is in good agreement with the prediction of the mc . time coincidence of signals above the energy threshold of 1 mev from both srd counters is required and high energy bremsstrahlung events are removed using the veto bgo . the suppression of synchrotron radiation emission detected for pions compared to electrons is clearly visible by comparing the two plots . for the electron spectrum , a 1% pileup beam events have been added to the simulation as predicted for the given spill intensity and with the known decay time of bgos . both spectra are in very good agreement with the simulation .            the efficiency for the electrons and the suppression factor for the pions are plotted in fig . [ fig : sup_mult ] as a function of the threshold on the energy deposited in the srd . we distinguish between two cases :    1 . the threshold is set on the total energy deposited in the srd . 2 .   both srd signals have to be in - time and above the threshold ( multiplicity requirement ) . one can see that applying the criterion 2 ) the efficiency only decreases slightly compared to 1 ) , while the suppression factor for pions is dramatically increased ( by two orders of magnitude ) with the requirement of having the two bgos in coincidence . this can be understood because the sr - like signal generated from secondary electron will leave a signal only in one of the two bgos while the sr from electrons is spread out uniformly as explained above . this is also nicely evidenced by tab . [ tab : hits ] where the fraction of events with different hit multiplicity in the srd bgo for both pion and electron runs are reported . .fraction of pion and electron events for different hit multiplicity in the srd from the data . [ cols="^,^,^,^",options="header " , ] in this paper we demonstrated that the detection of synchrotron radiation using bgo crystals is a very powerful method to tag 100 gev electrons bending in a 7 t@xmath52 m magnetic field . discriminating on the total energy deposited in the srd by setting a threshold of 10 mev , the contamination of pions in the beam can be suppressed down to a level of @xmath44 ( see fig . [ fig : sup_mult ] ) . moreover , we have shown that exploiting the granularity of the detector one can suppress the signal from secondary electrons generated via hadron ionization resulting in a suppression of hadron contamination down to a level @xmath1 maintaining an electron identification efficiency of @xmath295% . according to the geant 4 simulation validated with our measurements a further improvement by an additional order of magnitude is expected by increasing the detector granularity . in order not to degrade the efficiency due to pileup effects at higher beam intensity of @xmath59 electrons / s as achievable in the h4 beam line at the cern sps , one should use faster detectors such as , e.g. a pbsc sampling calorimeter or lyso crystals . these options will be tested by the na64 collaboration in the near future . we gratefully acknowledge the support of the cern management and staff and the technical staffs of the participating institutions for their vital contributions . this work was supported by the hiskp , university of bonn ( germany ) , jinr ( dubna ) , mon and ras ( russia ) , the russian federation program `` nauka '' ( contract no . 0.1764.gzb.2017 ) , eth zurich and snsf grant no . 169133 ( switzerland ) , and grants fondecyt 1140471 and 1150792 , pia / ring act1406 and pia / basal fb0821 conicyt ( chile ) . part of the work on mc simulations was supported by the rsf grant 14 - 12 - 01430 . we thank the compass daq group and the institute for hadronic structure and fundamental symmetries of tu munich for the technical support . a. andreas et al . , `` proposal for an experiment to search for light dark matter at the sps '' , cern - spsc-2013 - 034 ; spsc - p-348 ( 2013 ) , arxiv:1312.3309 [ hep - ex ] . d. banerjee et al . [ na64 collaboration ] , `` search for invisible decays of sub - gev dark photons in missing - energy events at the cern sps '' , phys . 118 , 011802 ( 2017 ) . j. alexander et al . , `` dark sectors 2016 workshop : community report '' , fermilab - conf-16 - 421 , arxiv:1608.08632 . l. b. okun , sov . jetp 56 ( 1982 ) 502 [ zh . eksp . 83 ( 1982 ) 892 ] . p. galison and a. manohar , phys . b136 , 279 ( 1984 ) . b. holdom , phys . b166 , 196 ( 1986 ) . h. wiedemann , `` synchrotron radiation '' , springer - verlag berlin heidelberg , 2003 . j. s. dworkin et al , `` electron identification using a synchrotron radiation detector '' , nim a247 , 412 ( 1986 ) . w. busza , m. chen and d. luckey , `` application of synchrotron radiation for the identification of electrons in high energy beams '' , nim 105 , 613 ( 1972 ) . r. crittenden and j. krider , `` synchrotron radiation detectors for tagging electrons '' , nim a 276 , 643 ( 1989 ) . j. allison et al . , nim a835 , 186 ( 2016 ) . g. f. knoll , `` radiation detection and measurement '' , wiley , 2002 . s. huber , j. friedrich , b. ketzer , i. konorov , m. kramer , a. mann , t. nagel , s. paul , ieee trans . 58 , 1719 ( 2011 ) . m. bodlak , v. frolov , v. jary , s. huber , i. konorov , d. levit , a. mann , j. novy , s. paul and m. virius , `` new data acquisition system for the compass experiment '' , jinst 8 , c02009 ( 2013 ) . v. avdeichikov et al . , nim a439 , 158 ( 2000 ) . k. a. olive et al . ( pdg ) , `` passage of particles through matter '' chin . c38 ( 2014 ) . d. banerjee , p. crivelli and a. rubbia , `` beam purity for light dark matter search in beam dump experiments '' , adv . high energy phys . , 105730 ( 2015 ) .
Abstract:

Answer: in high energy experiments such as active beam dump searches for rare decays and missing energy events , the beam purity is a crucial parameter . in this paper we present a technique to reject heavy charged particle contamination in the 100 gev electron beam of the h4 beam line at cern sps . the method is based on the detection with bgo scintillators of the synchrotron radiation emitted by the electrons passing through a bending dipole magnet . a 100 gev @xmath0 beam is used to test the method in the na64 experiment resulting in a suppression factor of @xmath1 while the efficiency for electron detection is @xmath295% . the spectra and the rejection factors are in very good agreement with the monte carlo simulation . the reported suppression factors are significantly better than previously achieved .    0.25 cm


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: it has been almost five decades , since the pioneer work of hnon and heiles [ 1 ] , who studied the regular and chaotic character of motion in a system of two coupled harmonic oscillators . during these years several methods have been presented in order to distinguish order from chaos . among them , it is worth mentioning , the method based on the fourier spectra of trajectories ( see [ 2,18 ] ) , the lyapunov characteristic exponents ( l.c.es)(see [ 19 ] ) and the ks entropy ( see [ 3 ] ) . apart from the above methods , in systems of two degrees of freedom , one can apply the classical method of the poincar phase plane , in order to distinguish regular from chaotic motion . on the contrary , in a 3d system the above method can not be applied . on this basis , it is of particular interest to locate the regions of initial conditions in 3d potentials generating regular or chaotic motion . the aim of the present article is to investigate the nature of motion in the 3d potential @xmath3 , \ \ \\end{aligned}\ ] ] where @xmath4 is the common frequency of oscillation along @xmath5 and @xmath6 axis , while @xmath7 is the perturbation strength . potential ( 1 ) represents three coupled harmonic oscillators in the 1:1:1 resonance . potentials of this type are also known as perturbed elliptic oscillators ( see [ 5,6 ] ) . the basic reason for the choice of potential ( 1 ) is that perturbed elliptic oscillators appear very often in galactic dynamics and atomic - particle physics ( see [ 5 ] and references therein ) . a second reason is that it displays exact periodic orbits , interesting sticky orbits together with large chaotic regions . therefore , it gives a good chance in order to check the reliability of the @xmath1 spectrum ( see [ 7,10 ] ) , which is used as an indicator to track the sticky orbits and to distinguish between regular and chaotic motion . the outcomes of this research are mainly based on the numerical integration of the equations of motion @xmath8 where the dot indicates derivatives with respect to the time . the hamiltonian to the potential ( 1 ) reads @xmath9 where @xmath10 and @xmath11 are the momenta per unit mass conjugate to @xmath5 and @xmath6 , while @xmath12 is the numerical value of the hamiltonian . the paper is organized as follows : in section 2 we present the analysis of the structure of the @xmath13 poincar phase plane of the 2d system for a set of values of energy @xmath12 . in the same section , we give some examples for the reader to understand the usefulness of the @xmath1 spectrum , since it has the ability to describe islandic and sticky motion . in section 3 we study the character of orbits in the 3d system . special interest is given to the evolution of the exact periodic orbits . particularly , we are interested to find the value of energy , where orbits starting near the exact periodic orbits change their character from regular to chaotic . we close our investigation with section 4 , where a discussion and the conclusions of this research are presented . in the following we shall study the character of motion in the corresponding 2d potential @xmath14 it is well known , that in this potential the two axes @xmath15 and @xmath16 and the straight lines @xmath17 are exact periodic orbits ( see [ 8 ] ) . the hamiltonian to potential ( 4 ) writes @xmath18 where @xmath19 is the numerical value of the hamiltonian . it is interesting to note , that potential ( 4 ) has not a finite energy of escape , so the corresponding zero velocity curves ( zvcs ) are always closed . in all numerical calculations , we adopt the values : @xmath20 , while the energy is treated as a parameter . figure 1a shows the @xmath13 , @xmath21 phase plane , when @xmath22 . one can see , that a large part of the phase plane is covered by chaotic orbits . the regular regions are occupied by invariant curves around the points @xmath23 and @xmath24 . the above points give the positions of two exact periodic orbits , the first is the @xmath16 axis while the second describes the @xmath17 straight line orbits on the @xmath13 phase plane . furthermore , one observes a large number of smaller islands of invariant curves produced by secondary resonances and some sticky regions as well . figure 1b is similar to fig . 1a but when @xmath25 . here we see that the majority of the phase plane is covered by a chaotic sea . the regular regions are confined only around the straight line periodic orbits , while the @xmath16 axis has now become unstable . some smaller islands of invariant curves are also present . the most interesting feature observed in this case , is the sticky regions around the chain of islands of invariant curves on the @xmath15 as well as on the @xmath26 axis . we shall come to this point later in this section . if we increase the value of energy to @xmath27 , we obtain the phase plane shown in figure 1c . in this case , the @xmath16 axis has returned to stability , while the straight line periodic orbits are unstable . almost all the phase plane is chaotic , except of a small region around the center and some small islands of invariant curves which are products of secondary resonances . a sticky region around the chain of the small islands of invariant curves near the center is also observed . figure 1d is similar to fig . 1a but when @xmath28 . this value of the energy is the value of the energy of escape of the 3d potential ( see next section ) . here almost all the phase plane is covered by chaotic orbits . there is a tiny regular region near the center . a careful observation also shows some very small regular islands of invariant curves embedded in the chaotic sea . 0.02 cm    one can see , that the fraction of the phase plane covered by chaotic orbits increases as the energy @xmath19 increases . figure 2a shows a plot of the percentage of the phase plane @xmath29 covered by chaotic orbits vs @xmath19 . we observe that @xmath29 increases rapidly , as @xmath12 increases . dots indicate values of @xmath29 found numerically , while the solid line is a fourth degree polynomial fitting curve . we would like to make clear , that @xmath29 is estimated on a completely empirical basis , by measuring the area on the @xmath13 phase plane occupied by chaotic orbits ( see [ 11 ] ) . in order to have an estimation of the degree of chaos from another point of view , we have plotted the maximum l.c.e vs @xmath19 . the results are shown in figure 2b . note that the l.c.e increases linearly with @xmath19 . here , we must point out that it is well known that the l.c.e has different values in each chaotic component ( see [ 20 ] ) . since we have regular regions in all cases and a large chaotic area with small sticky regions embedded in each chaotic sea , we calculate the average value of the l.c.e in each case by taking fifty orbits with different initial conditions @xmath30 , in each chaotic domain . in all cases , the obtained values of the l.c.es were different in the fourth decimal point , in the same chaotic area . 0.02 cm    0.02 cm    0.02 cm    0.02 cm    interesting information not only for the chaotic or regular character of orbits but also for islandic and sticky motion can be obtained using the dynamical parameter of the @xmath1 spectrum . the parameter @xmath31 is defined as @xmath32 where @xmath33 are the successive values of the @xmath34 orbital elements on the poincar @xmath13 phase plane . the dynamical spectrum of the parameter @xmath31 is its distribution function @xmath35 where @xmath36 are the number of the parameters @xmath31 in the interval @xmath37 , after @xmath38 iterations . this dynamical spectrum has been proved a very reliable tool in several cases ( see [ 7,10 ] for details ) . the nature of a 2d orbit can be identified by the shape of the @xmath1 spectrum . if the structure of the spectrum is a well formed @xmath39-type shape , then the corresponding orbit is regular . on the other hand , if the shape of the spectrum is complicated , asymmetric , with a lot of of small and large abrupt peaks , then this indicates that the orbit is chaotic . moreover , the @xmath1 spectrum can help us identify resonant orbits or orbits of higher multiplicity , as it produces as much @xmath39-type spectra as the total number of the islands of invariant curves in the @xmath13 phase plane . one additional advantage of this spectrum , is that it can be used in order to calculate the sticky period of a 2d sticky orbit and also to help us follow its time evolution towards the chaotic sea .    at this point , we should emphasize that there are some spectra , such as the @xmath40 spectrum , ( see [ 23 ] ) , which have no ability to identify islandic or sticky motion . furthermore , the @xmath40 spectrum needs the computation of two neighboring orbits and a large number of iterations ( about @xmath41 ) in order to provide reliable results regarding the character of a 2d orbit . using the @xmath1 spectrum we have managed to overcome these drawbacks . recently , in [ 24 ] we introduced the @xmath42 spectrum which proved to be a very sensitive and reliable spectrum . this new spectrum definition allow us to study more accurately the islandic motion of resonant orbits of higher multiplicity , as it produces as much spectra as the total number of the islands of invariant curves in the @xmath13 phase plane . the main advantage of the @xmath42 is that it has the ability to detect complicated resonant orbits especially in 2d hamiltonian systems , no matter how small are the corresponding islands of invariant curves . furthermore , for the study of orbits in 3d hamiltonian systems , we constructed and used in [ 25 ] the @xmath43 dynamical spectrum . this spectrum is an advanced form of the @xmath1 spectrum and it was especially designed for the study of 3d orbits , as it produces as much spectra as the total number of the 3d invariant tori in the @xmath44 phase subspace .    in what follows we present some examples in order to visualize the usefulness of the @xmath1 spectrum . figure 3a - b shows the @xmath1 spectrum and the corresponding l.c.e vs time of an orbit with initial conditions : @xmath45 and @xmath22 , while in all cases the initial value of @xmath46 is always found from the energy integral . we see a @xmath39 type spectrum indicating regular motion while the corresponding l.c.e tends asymptotically to zero . figure 4a - b is similar to fig . 3a - b for an orbit with initial conditions : @xmath47 and @xmath25 . here we see a spectrum with a lot of asymmetric small and large peaks indicating chaotic motion . the plot of the l.c.e vs time in fig . 4b verifies the results of the spectrum . the @xmath1 spectrum is useful in order to detect islandic and sticky motion . the results are shown in figure 5a - d . figure 5a shows an orbit with initial conditions : @xmath48 and @xmath25 . this orbit produces the two sets of five islands shown near the @xmath15-axis in fig . 1b . note that the corresponding spectrum produces ten individual regular spectra each one corresponding to an island . as the orbit produces two different sets of five islands on the positive and negative side of the @xmath15-axis , we observe two sets of five spectra almost symmetric to the origin . figure 5b shows the spectrum of a nearby orbit . the initial conditions are : @xmath49 and @xmath19=2 . this orbit produces the two sticky regions around each set of the five islands near the @xmath15-axis shown in fig . 1b . as there are two areas of sticky motion we see two separate sticky spectra with asymmetric large and small peaks . here the time is @xmath50 time units . figure 5c is similar to 5b but when @xmath51 time units . we see two different asymmetric spectra . this indicates that the test particle is still in the sticky region . figure 5d shows the spectrum when @xmath52 time units . in this case , we see that the two spectra have joined to produce a unified chaotic spectrum . this indicates that the test particle has left the sticky region and has gone into the chaotic sea . the l.c.e vs time of the orbit described above is shown in figure 6 . all l.c.es were computed for a time period of @xmath53 time units . the above results strongly suggest that the @xmath1 spectrum is a powerful dynamical parameter providing useful and reliable results for distinguishing between order and chaos . furthermore , the @xmath1 spectrum describes very satisfactorily islandic and sticky motion . in this section we shall investigate the regular or chaotic behavior of the orbits in the 3d potential ( 1 ) . working in a similar way as in caranicolas and varvoglis [ 9 ] , we find that this potential has an energy of escape given by @xmath54 for values of @xmath55 the zero velocity surface is open and the test particle is free to escape to infinity . in this paper we investigate only bounded motion that is we consider always @xmath56 . for the adopted values of the parameters @xmath20 we find that @xmath57 . figure 7 shows the surface @xmath58 . note that in the potential ( 1 ) the straight lines @xmath59 are exact periodic orbits . using the origin @xmath60 as a starting point we find from hamiltonian ( 3 ) that the starting value of the corresponding velocities are @xmath61 , while the value of @xmath46 is obtained from the energy integral . 0.02 cm    0.02 cm    0.02 cm    0.02 cm    in order to find the regular or chaotic nature of the 3d orbits we choose initial conditions @xmath62 , where @xmath63 is a point on the phase plane of the 2d system . this point lies inside the limiting curve @xmath64 which is the curve containing all the invariant curves of the 2d system . we use @xmath65 and the value of @xmath46 , for all orbits , is obtained from the energy integral ( 3 ) . the numerical calculations suggest , that orbits with initial conditions @xmath62 , such as @xmath63 is a point in the chaotic regions of fig . 1 a - d and for all permissible values of @xmath66 , produce chaotic orbits .    on the other hand , it would be also interesting to investigate the character of orbits with initial conditions @xmath62 , such as @xmath63 is a point in the regular regions of fig . 1 a - d . we shall use the @xmath1 spectrum and in order to verify the effectiveness of the @xmath1 spectrum in 3d dynamical systems , we shall compare the results with two other indicators , the classical method of the l.c.e and the spectral method @xmath2 , used by karanis and vozikis [ 18 ] . this method , uses the fast fourier transform ( fft ) of a series of time intervals , each one representing the time that elapsed between two successive points on the poincar @xmath13 phase plane , for 2d systems , while for 3d systems they take two successive points on the plane @xmath67 . note that the coupling of the third component , @xmath6 , carrying all the information about the 3d orbits , is hidden in the definition of the @xmath1 spectrum , but in any case it affects the values of @xmath68 and @xmath69 entering eq . thus , the @xmath1 spectrum provides implicit results for the 3d orbits . on the other hand , using the new @xmath43 spectrum the outcomes are explicit , since the coupling of the third component @xmath6 is located directly in its definition . figure 8 a - d shows ( a ) the 3d orbit , ( b ) the corresponding l.c.e . ( c ) the @xmath1 spectrum and ( d ) the @xmath2 indicator , of an orbit with initial conditions : @xmath70 . the value of energy is @xmath71 , while in all cases the value of @xmath46 is found from the energy integral . all dynamical indicators coincide that this orbit is regular . figure 9a - d is similar to fig . 8a - d but for a chaotic 3d orbit . initial conditions are : @xmath72 , @xmath71 . here one observes an asymmetric @xmath1 spectrum with a lot of small and large peaks indicating chaotic motion . the l.c.e in figure 9b and the @xmath2 indicator in figure 9d , coincide with the results of the spectrum . finally , we shall study the character of orbits starting near the 3d exact periodic orbits . figure 10 a - d is similar to fig . 8 a - d but for an orbit starting near the 3d exact periodic orbit . initial conditions are : @xmath73 , @xmath71 . the l.c.e , the @xmath1 spectrum and the @xmath2 indicator suggest for regular motion . figure 11 a - d is similar to fig . 9 a - d but when @xmath74 . here the dynamical indicators reveal the chaotic nature of this orbit . the numerical calculations indicate that orbits starting near the exact periodic orbits become chaotic when @xmath75 . figure 12 shows a plot of the l.c.e vs @xmath12 for the chaotic orbits in the 3d system . we see that the l.c.e increases linearly as @xmath12 increases . note that the values of the l.c.es of the 3d system are smaller than the corresponding values of the l.c.es of the 2d system . potentials composed of harmonic oscillators have been frequently used during the last years by several researchers [ 4 , 12 - 16 ] . in this article we have studied the regular or chaotic character of orbits in a 3d potential made up of perturbed elliptic oscillators . this potential displays exact periodic orbits , that are straight lines going through the origin , together with extended chaotic regions . our investigation begins from the 2d potential . in order to keep things simple we tried to study the regular or chaotic character of orbits for fixed values of the two other parameters and different values of the energy . our numerical calculations have shown that the extend of chaotic regions increases rapidly as the value of the energy increases . on the other hand , we have found that the degree of chaos also increases as the energy of increases . the numerical experiments show that the l.c.e increases linearly with @xmath19 . furthermore , the application of the @xmath1 spectrum shows that this dynamical parameter is a fast and reliable indicator in order to distinguish order from chaos . note that in order to obtain the @xmath1 spectrum we need only one orbit and 3000 - 10000 time units , while in order to obtain the l.c.e we use two orbits and much longer time periods . moreover , the @xmath1 spectrum describes in a very satisfactorily way islandic and sticky motion . on the other hand , the time needed in order to obtain the @xmath2 indicator , is of order of 5000 time units .    using the results of the 2d potential we proceeded to study the 3d model . first we tried to find the character of the 3d orbits with @xmath65 and initial conditions @xmath62 , such as @xmath63 is a point in ( i ) in the chaotic regions or ( ii ) in the regular regions of fig . the numerical results indicate that all orbits of type ( i ) are chaotic . the orbits of type ( ii ) display a different behavior . the numerical outcomes suggest that for small values of @xmath6 @xmath76 the orbits remain regular , while for larger values of @xmath6 @xmath77 the motion becomes chaotic . the regular or chaotic nature of the 3d orbits can be found using the @xmath1 spectrum . comparison with the l.c.e and the @xmath2 spectral method , shows that the @xmath1 spectrum is a good and fast indicator in order to distinguish regular from chaotic motion , in 3d dynamical systems . in the present research , we have choose to use dynamical parameters , such as the @xmath1 spectrum and the @xmath2 indicator , in order to characterize the nature of orbits , because the main advantage of these new methods is that they can provide reliable and conclusive results in short integration time periods , that is @xmath78 time units . here , we must point out that there are also other new dynamical indicators , such as the fli [ 17 ] , sali [ 21 ] and gali [ 22 ] . these indicators need at most 500 to 1000 time units of integration time , in order to reveal the true nature of an orbit in hamiltonian systems . the main disadvantage of the above methods is that we can not use them to identify resonant orbits of higher multiplicity which correspond to multiple islands of invariant curves in the phase plane . moreover , fli and gali are not very sensitive in the case of sticky motion . on the contrary , the sali method can provide reliable information about sticky orbits , using shorter integration time than the s(c ) spectrum .    note that an interesting characteristic of this investigations is that both , the 2d and 3d potentials display exact periodic orbits and large chaotic regions . the numerical results , in the 2d and the 3d model , show that , for small values of the energy , orbits starting near the exact periodic orbits remain regular , while for larger values of the energy become chaotic . the authors would like to thank the two anonymous referees for their very useful and also constructively suggestions and comments , which improved the quality of the present work . ferrer , s. , lara , m. , palacin , j. , juan , j.s . , viartola , a. , yanguas , p. : the hnon  heiles problem in three dimensions . ii . relative equilibria and bifurcations in the reduced problem . j. bifurc . chaos * 8 * , 1215 ( 1998 )
Abstract:

Answer: we study the nature of motion in a 3d potential composed of perturbed elliptic oscillators . our technique is to use the results obtained from the 2d potential in order to find the initial conditions generating regular or chaotic orbits in the 3d potential . both 2d and 3d potentials display exact periodic orbits together with extended chaotic regions . numerical experiments suggest , that the degree of chaos increases rapidly , as the energy of the test particle increases . about @xmath0 of the phase plane of the 2d system is covered by chaotic orbits for large energies . the regular or chaotic character of the 2d orbits is checked using the @xmath1 dynamical spectrum , while for the 3d potential we use the @xmath1 spectrum , along with the @xmath2 spectral method . comparison with other dynamical indicators shows that the @xmath1 spectrum gives fast and reliable information about the character of motion . = 1