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: to present a global panorama of the dynamics of a system by plotting poincar maps has a long tradition for two degrees of freedom ( 2-dof ) systems where these maps live on a 2 dimensional domain ( for a good explanation of the idea of poincar maps see chapter 6 in @xcite ) . unfortunately things are not so easy for more degrees of freedom . already for 3-dof the domain of the poincar map has dimension 4 and a graphical representation of the map is impossible . however , such 4 dimensional maps usually have the property that there are 2 dimensional invariant surfaces in the domain , i.e. surfaces with the property that initial conditions on these lower dimensional surfaces lead to trajectories which lie completely in this lower dimensional surface . and then it is easy to present graphically the restriction of the map to any one of such surfaces . such restricted 2 dimensional maps are a great help to understand important properties of the full 4 dimensional map .    to be useful such 2 dimensional invariant surfaces must have some kind of robustness property , i.e. they must survive parameter changes such that we can follow the perturbation scenario in the restricted map . there are two types of surfaces with this property . first , the so called normally hyperbolic invariant manifolds ( abbreviated nhims , for their general properties see @xcite ) , however , it seems that these objects do not play any important role in our present galaxy model . second , lower dimensional surfaces may be invariant because of symmetry reasons , and in the present galaxy model there are various surfaces of this type which will be very useful in the following . a further contribution to the understanding of 4 dimensional maps comes from a possible partially integrable limit case of the full system and from the corresponding possibility to build up the 4 dimensional map as stack of 2 dimensional reduced maps . because of its importance for the present work let us give a short summary of the stack construction in a form appropriate for the case found in the dynamics of the galactic potential . assume a hamiltonian 3-dof system where the unperturbed case has a rotational symmetry around the @xmath0-axis and where correspondingly @xmath1 , the @xmath0-component of the angular momentum , is conserved . then the unperturbed system can be reduced to a 2-dof system which depends parametrically on the value of @xmath1 . this leads to a 1 parameter family of reduced poincar maps each one acting on a 2 dimensional domain . from the reduced maps we obtain the non - reduced poincar map acting on its 4 dimensional domain in a two step process first , we pile up the continuum of 2 dimensional reduced maps to a 3 dimensional stack where @xmath1 acts as stack parameter . second , we form the cartesian product of this pile with a circle representing the cyclic angle canonically conjugate to the conserved quantity @xmath1 . the result is a 4 dimensional construction . the full 4 dimensional poincar map ( still for the unperturbed case ) acts on the resulting 4 dimensional domain as reduced map in each invariant horizontal leaf of constant @xmath1 and in addition applies a rotation to each copy of the circle . if we now add a perturbation which destroys the rotational symmetry around the @xmath0 axis then the invariant foliation into the horizontal leaves is destroyed . however , if the perturbation conserves some discrete symmetry then the 4 dimensional map will have corresponding invariant 2 dimensional surfaces @xmath2 also after the perturbation . these 2 dimensional surfaces do not coincide with any one of the invariant 2 dimensional horizontal leaves of the unperturbed stack , they are transverse to the stack structure . most important , the dynamics restricted to @xmath2 is sensitive to the perturbation and can be exploited to obtain an understanding of the perturbation scenario for the full 3-dof system . the restricted map on @xmath2 visualizes the decay of the stack structure under the perturbation . the idea works the same if the role of @xmath2 is taken over by a nhim . for all details of this approach and for examples of the stack construction see @xcite . for the analysis of a 3-dof molecular system taking advantage of this procedure see @xcite . the present article shows that the method is equally successful for the investigation of the dynamics of a 3-dof galaxy model .    in order to explore the orbital dynamics of galaxies one should build first suitable models describing sufficiently and realistically the properties of the galaxy . usually observations provide the necessary information on the construction of the dynamical models . a galactic model can be characterized as successful and realistic only if the derived results agree with the corresponding observational data . in most of the cases , the galaxy models are either spherical or axially symmetric . for instance , in a spherically symmetric potential all three components of the angular momentum and of course the total angular momentum are conserved . therefore , the motion of the stars is plane and takes place in the plane perpendicular to the vector of the total angular momentum . spherically symmetric models for galaxies were studied by @xcite . in an axially symmetric potential on the other hand , only the @xmath0 component of the total angular momentum is conserved . many previous papers are devoted to the distinction between order and chaos in axially symmetric potentials ( see e.g. , @xcite ) . furthermore , axially symmetric galaxy models were presented and examined in @xcite . another interesting category of galactic potentials are the so - called composite galaxy models . in those models the potential is multi - component and each component describes a different part of the system . composite axially symmetric galaxy models describing the motion of stars in our galaxy were also studied by @xcite . in these models the gravitational potential is composed by three superposed disks : one representing the gas layer , one the thin disk and one for the thick disk . recently in @xcite , a class of realistic triaxial models for galaxies was provided . in particular , the authors extended an earlier method proposed by @xcite to three - dimensional systems by replacing the radial with an ellipsoidal symmetry in the total mass density . moreover , triaxial galaxy models were also constructed by @xcite and @xcite . the layout of the paper is as follows : section [ galmod ] , contains a detailed presentation of the structure and the properties of our galactic gravitational model . in section [ numres ] , we construct poincar maps in order to investigate the orbital properties of the 6 dimensional phase space . the paper ends with section [ disc ] , where the discussion and the main conclusions of our numerical analysis are presented . the total gravitational potential @xmath3 is three - dimensional and it consists of three components : a central , spherical component @xmath4 , a flat disk @xmath5 and a triaxial dark matter halo potential @xmath6 . the spherically symmetric nucleus is modeled by a plummer potential @xcite @xmath7 here @xmath8 is the gravitational constant , while @xmath9 and @xmath10 are the mass and the scale length of the nucleus , respectively . here we must point out , that potential ( [ vn ] ) is not intended to represent a black hole nor any other compact object , but a dense and massive bulge . therefore , we do nt include any relativistic effects .    in order to model the disk we use the miyamoto - nagai potential @xcite @xmath11 where @xmath12 is the mass of the disk , while @xmath13 and @xmath14 are the horizontal and vertical scale lengths of the disk . for the description of the properties of the dark matter halo we use the logarithmic potential @xmath15 where @xmath16 and @xmath17 are the flattening parameters along the @xmath18 and @xmath0 axes , respectively , @xmath19 is the scale length of the halo , while the parameter @xmath20 is used for the consistency of the galactic units . the choice for the logarithmic potential was motivated for several reasons : ( i ) it can model a wide variety of shapes of galactic haloes by suitably choosing the parameter @xmath17 . in particular , when @xmath21 the dark matter halo is prolate , when @xmath22 is spherical , while when @xmath23 is oblate ; ( ii ) it is appropriate for the description of motion in a dark matter halo as it produces a flat rotation curve at large radii ( see fig . [ rotvel ] ) ; ( iii ) it allows for the investigation of flattened configurations of the galactic halo at low computational costs ; ( iv ) the relatively small number of input parameters of eq . ( [ vh ] ) is an advantage concerning the performance and speed of the numerical model ; and ( v ) the flattened logarithmic potential was utilized successfully in previous works to model a dark matter halo component ( e.g. , @xcite ) . we use a system of galactic units , where the unit of length is 1 kpc , the unit of mass is @xmath24 ( solar masses ) and the unit of time is @xmath25 yr ( about 100 myr ) . the velocity units is 10 km / s , the unit of angular momentum ( per unit mass ) is 10 km kpc s@xmath26 , while @xmath8 is equal to unity . the energy unit ( per unit mass ) is 100 km@xmath27s@xmath28 . in these units , the values of the involved parameters are : @xmath29 @xmath30 , @xmath31 , @xmath32 , @xmath33 , @xmath34 and @xmath35 , while @xmath16 and @xmath17 are treated as variable parameters in the intervals @xmath36 and @xmath37 . the values of the disk and the nucleus were chosen with a milky way - type galaxy in mind ( e.g. , @xcite ) and they secure positive density everywhere . here we must point out that our total gravitational potential is truncated at @xmath38 kpc , otherwise the total mass of the galaxy modeled by this potential is infinite which of course lacks of physical meaning .        for @xmath39 the total gravitational potential @xmath3 becomes axially symmetric @xmath40 , where @xmath41 and the @xmath0-component of the angular momentum is conserved . one of the most important quantities in disk galaxies is the circular velocity in the galactic plane which is defined as , @xmath42 a plot of @xmath43 for our galactic model , so the rotation curve is the same regardless of the particular value of @xmath17 . ] is presented in fig . [ rotvel ] as a black curve . moreover , in the same plot , the red line shows the contribution from the spherical nucleus , the green curve is the contribution from the disk , while the blue line corresponds to the contribution form the dark matter halo . it is seen that each contribution prevails in different distances form the galactic center . in particular , at small distances when @xmath44 kpc , the contribution from the spherical nucleus dominates , while at mediocre distances , @xmath45 kpc , the disk contribution is the dominant factor . on the other hand , at large enough galactocentric distances , @xmath46 kpc , we see that the contribution from the dark halo prevails , thus forcing the rotation curve to remain flat with increasing distance from the center . we also observe the characteristic local minimum of the rotation curve due to the massive nucleus , which appears when fitting observed data to a galactic model ( e.g. , @xcite ) . the equations of motion for a test particle with a unit mass @xmath47 are @xmath48 where , as usual , the dot indicates derivative with respect to the time . the hamiltonian for the total potential @xmath3 reads @xmath49 where @xmath50 , @xmath51 and @xmath52 are the momenta per unit mass conjugate to @xmath53 , @xmath18 and @xmath0 , respectively , while @xmath54 is the numerical value of the hamiltonian , which is conserved . it is rather obvious what is the appropriate partially integrable limiting case of our galaxy model . it is the rotationally symmetric case with @xmath39 . then in cylindrical coordinates @xmath55 the potential is independent of @xmath56 and correspondingly the canonically conjugate action @xmath1 , the @xmath0 component of the total angular momentum , is conserved . and we have a foliation of the system into leaves of constant @xmath1 . first we treat the case of @xmath57 and @xmath58 in all detail and we give short remarks on other values of these two quantities later . here , we have to point out that the energy level controls the size of the maps and particularly @xmath59 which is the maximum possible value of the @xmath53 coordinate . we chose that energy level which yields @xmath60 kpc .    from the functional form of the hamiltonian in eq . ( [ ham ] ) it is also easy to see that the surfaces @xmath61 given by @xmath62 , @xmath63 given by @xmath64 , and @xmath65 given by @xmath66 are invariant . any initial condition starting in one of these surfaces leads to a trajectory which lies entirely in the respective surface . now we come to the question which intersection condition to use for the poincar map . it must be an intersection condition which is appropriate to investigate the restricted dynamics in the surfaces @xmath61 , @xmath63 and @xmath65 and it must also be appropriate for the stack construction in the symmetric case . the best intersection condition we found fulfilling all these side conditions is the condition that the cylindrical radius @xmath67 is maximal , or equivalently that its conjugate momentum @xmath68 changes from positive to negative . the 4 dimensional poincar map generated with this intersection condition will be called @xmath69 in the following . the restrictions of this map to the invariant surfaces @xmath61 , @xmath63 and @xmath65 will be called @xmath70 , @xmath71 and @xmath72 , respectively . as in most cases for realistic systems the transversality of the intersection condition is violated for a few important trajectories . in our case these trajectories include two simple ones . first the trajectory @xmath73 which oscillates up and down along the z - axis with @xmath67 constantly zero . and second , in the symmetric case the trajectory @xmath74 which encircles the origin in a circle with constant value of @xmath67 and in the plane @xmath75 . as we will see later these periodic trajectories appear in some restricted maps as energetic boundaries . next we explain the stack idea in detail . consider the symmetric case with @xmath39 , where in cylindrical coordinates the coordinate @xmath56 does not appear in the hamiltonian and where @xmath1 is conserved and correspondingly can be treated as if it would be a parameter . because of the rotational symmetry the angle @xmath56 will be ignored for the moment . then we have a foliation of the phase space and also of the domain of the poincar map into invariant leaves belonging to the various values of @xmath1 and each leaf can be treated as independent dynamical system . the reduced system is 2-dof ( @xmath67 and @xmath0 degrees of freedom ) and its reduced poincar map ( called @xmath76 in the following ) acts on the 2 dimensional domain with the coordinates @xmath0 and @xmath52 . in fig . [ map1 ] we show the numerical plots of @xmath76 for the examples of the @xmath1 values 0 , 30 , 60 , 90 , 120 , and 150 in parts ( a ) , ( b ) , ( c ) , ( d ) , ( e ) , and ( f ) , respectively . we observe that with increasing value of @xmath1 the domain of the reduced map shrinks and the map becomes more regular . for the value 160 of the total energy @xmath54 the maximal possible value of @xmath1 is found as @xmath77 and when @xmath1 converges to this value from below then the domain of the reduced map shrinks to a single point , the origin . and in this limit case this single point represents the periodic trajectory @xmath74 . fig . [ chaos ] shows the relative fraction of area covered by chaotic trajectories as function of @xmath1 . it demonstrates in a quantitative form how this fraction converges to 0 for increasing value of @xmath1 . for the determination of the chaotic percentage we used an image processing numerical method which is explained in the appendix of course , there are also corresponding leaves with negative values of @xmath1 and the opposite orientation of rotation around the @xmath0-axis . however , because the hamiltonian is even in @xmath1 these cases are nothing new and only are mirror images of the cases with positive values of @xmath1 . therefore , there is no need to show numerical examples for negative values of @xmath1 .     when @xmath39 and @xmath58 . ]     as a function of the angular momentum @xmath1 . ]    for the case of @xmath78 shown in part ( a ) of fig . [ map1 ] any trajectory moves in a 2 dimensional plane in position space which contains the @xmath0-axis . in addition for @xmath39 the motion is independent of the angle @xmath56 of this plane . therefore for @xmath39 and @xmath78 @xmath76 coincides with the maps @xmath70 and @xmath71 . next we explain how we pile up the stack of reduced systems to represent the 4 dimensional poincar map . first , imagine that we construct the reduced map for all values of @xmath1 in the interval @xmath79 $ ] and pile up this continuum of 2 dimensional maps to a 3 dimensional stack where the vertical coordinate ( the stack parameter ) is @xmath1 . second , we form the cartesian product of this stack with a circle representing the cyclic angle @xmath56 . thereby we arrive at a 4 dimensional object . the natural coordinates of this object are @xmath0 , @xmath52 , @xmath56 and @xmath1 , i.e. exactly the 4 coordinates of the domain of @xmath69 , i.e. this 4 dimensional object can be identified with the domain of the 4 dimensional poincar map @xmath69 . finally , we have to understand how @xmath69 acts on this construction . in each copy of the @xmath0 - @xmath52 plane the map acts according to @xmath76 for the respective value of @xmath1 . @xmath1 itself is conserved . and each copy of the circle suffers a rigid rotation . [ rot ] shows @xmath80 as function of @xmath1 for the trajectory moving in the equatorial plane @xmath75 . by choosing the equatorial plane we obtain exactly that twist curve which we need in a moment to understand the perturbed map @xmath72 . it is only plotted up to a value of @xmath81 since higher values are irrelevant for the real galaxy . however , the function continues up to @xmath82 where the rotation angle reaches a value around -1.2 . for negative values of @xmath1 the function should be continued antisymmetrically . after having reached the global understanding of the 4 dimensional map for @xmath39 we have to study the perturbation scenario of this map , i.e. its behaviour under changes of @xmath16 , or even better we should use @xmath83 as perturbation parameter . the perturbation has a strong effect , if the @xmath56 dependence of the perturbation comes in resonance with the rotation angle @xmath80 of the unperturbed system . in our case the perturbation appears in the argument of the logarithm of the halo potential only and in this argument we set @xmath84 in this sense the perturbation consists of a fourier component of order 2 only and has most effect when the rotation angle @xmath80 has values around @xmath85 . [ rot ] shows , this happens for @xmath1 values around 0 . therefore the perturbation will have strong effects around @xmath86 also for small values of @xmath87 . higher order resonance effects are expected to be small for small values of @xmath87 and to become efficient for large values of @xmath87 only . the next important resonances are a 1:3 resonance for @xmath88 , a 1:4 resonance for @xmath89 and a 3:8 resonance for @xmath90 . corresponding resonances exist for negative values of @xmath1 . the perturbation scenario is presented in the most evident form by the restricted map @xmath72 on its domain with coordinates @xmath56 and @xmath1 . these are the canonical coordinates of the @xmath56 degree of freedom and this degree of freedom is exactly the one which is decoupled from the rest of the system in the symmetric case . in the symmetric case @xmath72 is very trivial since @xmath1 is conserved and the domain is foliated into the invariant lines of constant @xmath1 ( such lines are circles of the variable @xmath56 , one for each value of @xmath1 ) . each of these circles is just rigidly shifted by the rotation angle @xmath91 . in this sense the map @xmath72 for @xmath39 is a pure twist map . under perturbations any twist maps follows the universal perturbation scenario of the standard map @xcite . for useful information on the interpretation of structures found in poincar maps see also chapter 4 in @xcite .    in fig . [ map2 ] we show 4 examples for the perturbed map @xmath72 . only the part of the domain up to @xmath81 is included in the plot since the model is only realistic up to this value . in principle the domain goes up to @xmath82 and the natural boundary of the domain is the trajectory @xmath74 . the plots should be supplemented symmetrically for negative values of @xmath1 .    for the moderate perturbation strength @xmath92 in part ( a ) of the figure the 1:2 resonance at @xmath78 has already a very large effect , it causes two large secondary islands and a fine chaos strip which reaches up to @xmath1 values around 45 . the centres of the secondary island are the two points @xmath93 and @xmath94 of an elliptic orbit of period two , and the centres of the chaos strips are the two points @xmath95 and @xmath96 of a hyperbolic orbit of period two . the fine chaos strip still comes very close to a separatrix curve and it indicates the position of the initial segments of the stable and unstable manifolds of the unstable orbit of period two . the secondary structure of the 1:3 resonance at @xmath1 values around 110 consists of two distinct island chains of period 3 . these two different island chains are transformed into each other by the discrete symmetry @xmath97 . the corresponding chaos strip is extremely fine and not displayed in the figure . all other secondary structures are still too small to become visible at this perturbation strength .    for @xmath98 in part ( b ) of the figure the secondary structure of the 1:2 resonance has grown to a width in @xmath1 of approximately 70 and its chaos strip has become wider so that it is clearly distinct from a separatrix curve . further , the fine chaos strip of the 1:3 resonance has become visible and the island chain of the 3:8 resonance has become visible . however , in total the map still is dominated by regular motion .    for the perturbation @xmath99 in part ( c ) of the figure the interior of the 1:2 resonance islands has become chaotic on a large fraction of the area . and we see the fine chaos strip belonging to the 3:8 resonance . finally in part ( d ) we present the case @xmath100 . here the 1:2 resonance island has already dissolved into a large chaotic sea with only a few surviving islands inside of the chaotic sea , they come from higher levels of the hierarchy . in addition the chaos strips of the 3:8 and of the 1:3 resonances have become so wide that we can distinguish them very clearly from separatrix curves . at least for @xmath1 values below 100 the map is now dominated by instability , whereas for large values of @xmath1 it is still mainly integrable . in total we see in fig . [ map2 ] a perturbation scenario following the standard map as it is expected for any generic perturbed twist map . it should be emphasized that the dynamics in the surface @xmath65 is independent of @xmath17 and therefore figs . [ rot ] and [ map2 ] are equally valid for any other value of @xmath17 . it is seen in fig . [ map2 ] that the most prominent centres for the perturbation scenario are the points of period 2 on the line @xmath78 . the elliptic points at @xmath101 , @xmath78 and at @xmath102 , @xmath78 correspond to a trajectory in position space which oscillates along the @xmath53-axis . and the hyperbolic points at @xmath103 , @xmath78 correspond to a trajectory in position space which oscillates along the @xmath18-axis . in fig . [ orbs2 ] we show two trajectories in position space with initial conditions close to these two central trajectories . the outermost solid line is the zero velocity curve at the @xmath104 plane defined as @xmath105 . we choose the examples for the perturbation @xmath98 corresponding to part ( b ) of fig . [ map2 ] . also the nearby trajectories shown start with initial conditions @xmath75 and @xmath106 , they lie in the invariant plane @xmath65 and thereby represent points in the poincar plot of fig . the trajectory plotted in fig . [ orbs2]a has initial conditions @xmath107 , @xmath108 . because it represents the corresponding point in fig . [ map2]b it must have the additional initial condition @xmath109 and the initial value of @xmath67 is given by the energy value @xmath57 resulting in @xmath110 . [ orbs2]a is the resulting trajectory in the equatorial plane @xmath65 plotted in coordinates @xmath53 and @xmath18 . it represents a small kam curve in fig . [ map2]b encircling the elliptic point of period 2 . it is the smallest kam curve around the elliptic point included in fig . the direction of oscillation of this trajectory fluctuates quasi - periodically around the values @xmath101 or equivalently @xmath102 . the width of these fluctuations in @xmath56 is given by the horizontal width of the corresponding kam curve in fig . the trajectory plotted in fig . [ orbs2]b has initial conditions @xmath111 , @xmath112 . also it has the additional initial condition @xmath109 and the initial value of @xmath67 is given by the energy value @xmath57 resulting in @xmath113 . [ orbs2]b is the resulting trajectory in the equatorial plane with coordinates @xmath53 and @xmath18 . it represents motion in the small chaos strip in fig . [ map2]b containing the hyperbolic point of period 2 . the direction of oscillation of this trajectory starts in @xmath18 direction , i.e. around the values @xmath103 with the initial value 2 of @xmath1 . then the direction slowly rotates away from the @xmath18 direction passes through the @xmath53 direction quite rapidly and returns slowly to the @xmath18 direction before making the next rotation of the direction of oscillation . because it is a chaotic trajectory the rotations of the direction of oscillation occurs irregularly and also the orientation of rotation changes in the long run again and again in an irregular fashion . when the direction of oscillation lies in the @xmath53 direction this corresponds to a point in the small chaos strip of fig . [ map2]b at @xmath56 values close to 0 or @xmath85 and correspondingly the angular momentum of the trajectory lies around 70 at this moment . accordingly the trajectory is more round at this moment .    to understand the 4 dimensional map @xmath69 for @xmath114 it is helpful to have in addition the restricted maps @xmath70 and @xmath71 . first , in the surface @xmath63 the dynamics is independent of @xmath16 . therefore the plot of the map @xmath71 always coincides with the plot of the map @xmath76 for @xmath78 as shown in fig . [ map1]a . for @xmath39 also the map @xmath70 coincides with the map of fig . . however , for @xmath114 the map @xmath70 shows a weak dependence on @xmath16 . in fig . [ map3 ] we show the cases of @xmath70 for @xmath98 and @xmath100 in parts ( a ) and ( b ) respectively . by comparing figs . [ map1]a , [ map3]a and [ map3]b we see that there is a rather small dependence of @xmath70 on @xmath16 , but only in fine details . the energetic boundaries of fig [ map3 ] are given by the trajectory @xmath73 . for this trajectory along the @xmath0-axis all energy is in the @xmath0 motion and the dynamics along this line is independent of the value of @xmath16 . finally we intent an approximate global understanding of the 4 dimensional map @xmath69 for general values of @xmath16 . the plane @xmath61 can also be characterized as plane given by @xmath78 and @xmath103 and the plane @xmath63 is at the same time the subset characterized by @xmath78 and @xmath101 or @xmath102 . the 4 dimensional domain of @xmath69 can be considered as cartesian product of the @xmath56 - @xmath1 space with the @xmath0 - @xmath52 space . over each point of the @xmath56 - @xmath1 space we imagine sitting one copy of the @xmath0 - @xmath52 space . each one of the two copies of the @xmath0 - @xmath52 space sitting over the points @xmath103 , @xmath78 can be identified with @xmath61 . these fibers are invariant and the restriction of @xmath69 to these fibers coincides with @xmath70 . equally each one of the two copies of the @xmath0 - @xmath52 space sitting over the points @xmath101 , @xmath78 and @xmath102 , @xmath78 can be identified with @xmath63 . these fibers are invariant and the restriction of @xmath69 to these fibers coincides with @xmath71 . then in the neighbourhood of such invariant planes we can imagine as a first approximation for the complete map @xmath69 the product of the restricted map @xmath72 and @xmath70 or @xmath72 and @xmath71 , respectively . for general points in the 4 dimensional domain this is a rather crude approximation but it helps to get some first imagination of the approximate action of the full map @xmath69 in 4 dimensions . so far we have considered the case @xmath58 only . in figs . [ map4 ] and [ map5 ] we show the plots necessary to understand the cases @xmath22 and @xmath115 . first we have to remember that the plots of figs . [ rot ] and [ map2 ] are independent of @xmath17 . in figs . [ map4 ] and [ map5 ] we show the plots of the map @xmath70 for the various combinations of @xmath16 values of 1 and 1.5 and @xmath17 values of 1 and 1.5 . remember that the plot of @xmath71 always coincides with the corresponding @xmath70 plot for @xmath39 . thereby we have the necessary information to understand also the cases @xmath22 and @xmath115 .    before closing this section we would like to present in fig . [ orbs3](a - c ) some examples of each of the basic types of regular resonant orbits . the orbit shown in fig . [ orbs3]a is a 0:1:1 resonant periodic orbit with initial conditions @xmath116 . a 0:2:1 resonant periodic orbit with initial conditions @xmath117 is given in fig . [ orbs3]b , while the initial conditions of the 0:3:2 periodic orbit of fig . [ orbs3]c are @xmath118 . the @xmath119 notation we use for the regular orbits is according to @xcite and @xcite , where the ratio of those integers corresponds to the ratio of the main frequencies of the orbit , where main frequency is the frequency of greatest amplitude in each coordinate . main amplitudes , when having a rational ratio , define the resonances of an orbit . for the periodic orbits of fig . [ orbs3](a - c ) we have @xmath120 because the trajectories lie in this invariant plane , while from the intersection condition we know that @xmath121 since for this particular situation @xmath51 coincides with @xmath68 and our intersection condition is @xmath109 . finally the initial value of @xmath18 , which coincides with @xmath67 in this case , has to be calculated by the total energy integral . the values of the parameters @xmath16 and @xmath17 are as in fig . the purpose of the whole work is to obtain knowledge for which parameter values ( here we think of parameters @xmath16 and @xmath17 that control the geometry of the dark matter halo ) and for which regions of phase space we have stability and instability and which types of trajectories play an important role . an inspection of figs . [ map1 ] , [ map2 ] , [ map3 ] , [ map4 ] , [ map5 ] demonstrates that the 1:2 resonance in the invariant plane @xmath65 ( most clearly shown in fig . [ map2 ] ) plays the most prominent role . it depends strongly on @xmath16 but is completely independent of @xmath17 . this importance is easy to understand since the functional form of the perturbation enforces a 1:2 coupling with the @xmath56 degree of freedom , i.e. the degree of freedom which is decoupled from the rest of the dynamics in the partially integrable limit case . because the resonance occurs around @xmath78 it mainly affects trajectories passing close to the center of the galaxy . as can be seen from the functional form of the halo potential , the @xmath18 direction ( or the angle values @xmath103 as seen from eq . ( [ pert ] ) ) is the one most affected by the perturbation , i.e. by the breaking of the rotational symmetry . with increasing @xmath16 the width of oscillations in @xmath18 direction shrinks . as seen in fig . [ map2 ] for moderate values of @xmath16 initial conditions close to @xmath78 and @xmath101 or @xmath102 lead to stable oscillations in @xmath53 direction . on the other hand , initial conditions close to @xmath78 and @xmath103 lead to unstable behaviour . the trajectory first oscillates in @xmath18 direction , then its direction of oscillation rotates and at the same time the angular momentum increases and becomes maximal when the trajectory oscillates in @xmath53 direction . then the direction of @xmath67 oscillation keeps on rotating while the angular momentum decreases again until it returns to @xmath86 when the direction of @xmath67 oscillation returns to the @xmath18 direction . this is the behaviour of trajectories in the chaos strip starting from the hyperbolic points @xmath78 , @xmath122 in plot shown in fig . [ map2 ] . for larger perturbation around @xmath123 the oscillation in @xmath53 direction becomes inverse hyperbolic , see the transition from fig . [ map2]b over [ map2]c to [ map2]d . then most motion with small values of @xmath1 becomes unstable .    on the other hand , the parts of fig . [ map1 ] for larger values of @xmath1 and fig . [ map2 ] show that motion with large values of @xmath1 always stays stable . this indicates that the approximately circular motion of the disc around the centre of the galaxy is always stable . for larger values of @xmath16 the secondary stable structures of the 1:3 resonance and the 3:8 resonance become moderate size . this leads to stable bundles of trajectories in the neighbourhood of the corresponding periodic orbits . influences of the parameter @xmath17 can only exist for trajectories which go far outside of the galactic plane @xmath75 . the most extreme case is the trajectory @xmath73 . [ map4 ] and [ map5 ] suggest the following : remember that the intersection condition is @xmath109 . then points in figs . [ map4 ] and [ map5 ] with large values of @xmath52 are trajectories which pass through the maximal value of @xmath67 when at the same time they have a large momentum in @xmath0 direction . and these are trajectories having their major part of the energy in the @xmath0 motion . figs . [ map4 ] and [ map5 ] indicate that such trajectories are regular for large values of @xmath17 . the situation is different for small @xmath17 as is evident from fig . here we have stable islands of moderate size only . and in particular , we see chaos close to the whole boundary . this indicates that the trajectory @xmath73 has become unstable . this is understandable : a small @xmath17 allows very wide oscillations in @xmath0 direction , wider than in @xmath53 or @xmath18 directions . this leads to overfocusing of the trajectory @xmath73 and therefore to inverse hyperbolic behaviour ( e.g. , @xcite ) .    in order to have a more complete view on how energy influences the orbital structure of the dynamical system we constructed additional maps for lower and higher values of the standard @xmath57 level . our numerical calculations suggest that for low energies , that is the case of local motion around the nucleus , the trajectories do not go out into the dark matter halo and therefore the values of the halo parameters are irrelevant . whereas for large values of the energy the trajectories go out far and are strongly influenced by geometrical details of the halo . this work has been supported by dgapa under grant number ig-101113 and by conacyt under grant number 79988 . the authors would also like to thank the anonymous referee for the careful reading of the manuscript and for all the apt suggestions and comments which allowed us to improve significantly both the quality and the clarity of the paper . allen c. , santilln a. , 1991 , rmaa , 22 , 255    bailin j. , simon j.d . , bolatto a.d . , gibson b.k . , power c. , 2007 , apj , 667 , 191    binney j. , 2012 , mnras , 426 , 1328    binney j. , tremaine s. , 2008 , galactic dynamics , princeton univ . press , princeton , usa    carpintero d.d . , aguilar l.a . , 1998 , mnras , 298 , 21    chirikov b.v . , 1979 , physics reports , 52 , 263    cretton n. , de zeeuw p.t . , van der marel r.p . , rix h.w . , 1999 , apj suppl . series , 124 , 383    dehnen w.d . , 1993 , mnras , 265 , 250    gmez f. , helmi a. , brown a. , li y. , 2010 , mnras , 408 , 935    gonzalez f. , jung c. , 2012 , j. phys . a : math . theor . , 45 , 265102    gonzalez f. , drotos g. , jung c. , 2014 , j. phys . a : math . theor . , 47 , 045101    helmi a. , 2004a , pasa , 21 , 212    helmi a. , 2004b , apj , 610 , l97    irrgang a. , wilcox b. , tucker e. , schiefelbein l. , 2013 , a&a , 549 , a137    jackson e.a . , 1990 , perspectives of nonlinear dynamics , cambridge university press    jung c. , merlo o. , seligman t.h . , zapfe w.p.k . , 2010 , new j. phys . ,  12 , 103021    keller s.c . , da costa g.s . , prior s.l . , 2009 , mnras , 394 , 1045    keller s.c . , murphy s. , prior s.l . , da costa g.s . , schmidt b. , 2008 , apj , 678 , 851    law d.r . , majewski s.r . , johnston k.v . , 2009 , apj let . , 703 , l67    lichtenberg a.j . , lieberman m.a . , 1983 , regular and stochastic motion , new york : springer    lin y - d . , reichl l.e . , jung c. , 2015 , j. chem . phys . , 142 , 124304    miyamoto w. , nagai r. , 1975 , pasj 27 , 533    moore b. , kazantzidis s. , diemand j. , stadel j. , 2004 , mnras , 354 , 522    prior s.l . , da costa g.s . , keller s.c . , 2009 , apj , 704 , 1327    rindler - daller t. , dejonghe h. , zeilinger w.w . , 2005 , mnras , 356 , 1403    ruika a. , palou j. , theis c. , 2007 , a&a , 46 , 155    terzi b. , graham a.w . , 2005 , mnras , 362 , 197    terzi b. , sprague b.j . , 2007 , mnras , 377 , 855    wiggins s. , 1994 , normally hyperbolic invariant manifolds in dynamical systems , berlin : springer    wolfram s. , 2003 , the mathematica book , wolfram media , inc    zhao h.s . , 1996 , mnras , 278 , 488    zotos e.e . , 2011 , new astronomy , 16 , 391    zotos e.e . , 2012 , new astronomy , 17 , 576    zotos e.e . , 2014 , a&a , 563 , a19    zotos e.e . , carpintero d.d . , 2013 , cemda , 116 , 417    zotos e.e . , caranicolas n.d . , 2014 , nonlinaer dynamics , 76 , 323 when constructing a poincar map it would be of particular interest to know how much area on the phase plane is covered by initial conditions that correspond to chaotic orbits . here we are going to present a simple but yet very effective numerical method for determining the chaotic percentage in poincar maps . to begin with it is important to have a poincar map constructed in fine detail . the density of points must be high enough that the following two conditions are fulfilled : first , in chaotic seas almost all pixels must be occupied . second , kam lines must appear as continuous curves without interruptions . after obtaining the data points from the numerical integration routine we load them into the software program _ @xcite and we plot them thus creating the map . this software program contains a plethora of image - processing tools . by using the appropriate functions we locate all the white holes in the map which correspond to stability islands of regular orbits . then the map is divided into two parts : ( i ) the mask of the full available region on the phase plane and ( ii ) the mask of the regular regions . for both parts we measure how many pixels each mask covers . therefore , if we assume that the mask of the full region occupies @xmath124 pixels and the mask of the regular regions occupies @xmath125 pixels then the chaotic area on the phase plane is @xmath126 . in fig . [ img ] we present an illuminating example combining the detected stability islands with the full poincar map of fig . [ map1]a in an attempt to to point out the great efficiency of the numerical method .
Abstract:

Answer: we explore the orbital dynamics of a realistic three dimensional model describing the properties of a disk galaxy with a spherically symmetric central dense nucleus and a triaxial dark matter halo component . regions of phase space with regular and chaotic motion are identified depending on the parameter values for triaxiality of the dark matter halo and for breaking the rotational symmetry . the four dimensional poincar map of the three degrees of freedom system is analysed by a study of its restriction to various two dimensional invariant subsets of its domain . = 1    hamiltonian systems  ordered and chaotic motion symplectic map of dimension four


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: slowly pulsating b - type supergiants ( spbsg ) have emerged as a new class of pulsating variables after saio et al . ( 2006 ) have detected pulsational frequencies in the most satellite data of the blue supergiant hd 163899 ( b2 ib / ii , klare & neckel 1977 , schmidt & carruthers 1996 ) . saio et al . ( 2006 ) identified 48 frequency peaks in the most light curve and attributed them to g- and p - mode pulsations . because so far only one object of this type is found , the spbsg class is not yet well defined as for the range of effective temperature and luminosity . also the evolutionary status of these pulsators remains unrevealed ; they may be either in the shell hydrogen burning phase as well as after core helium ignition . the finding of saio et al . ( 2006 ) has prompted a few groups ( godart et al . 2009 , daszyska - daszkiewicz , ostrowski , pamyatnykh . 2013 ; hereafter d2013 ) to reanalyse pulsation stability in models of b - type stars after the terminal age main sequence ( tams ) and to further studies of spbsg variables . the presence of g - mode pulsations in b - type post - main sequence stars has been explained by a partial reflection of some modes at an intermediate convective zone ( icz ) related to the hydrogen - burning shell or at a chemical gradient zone surrounding the radiative core . however , all studies of these objects published so far are based on the assumption that hd 163899 has not reached the phase of core helium ignition , i.e. , it is in the phase of shell hydrogen burning . this assumption might not necessary be fulfilled , because the blue loop can reach temperatures of early b spectral types . in this paper we calculate models which undergo core helium burning on the blue loop and compare them with models in the hydrogen shell burning phase . the structure of the paper is as follows . in section [ domains ] , we present the new instability domains which include pulsational instability on the blue loops . the effects of opacity , metallicity , overshooting , element diffusion and mass loss on the pulsational instability areas are studied . propagation diagrams and properties of instability parameter and kinetic energy of modes for representative models are described in section [ models ] . in section [ identification ] we construct photometric diagnostic diagrams for the b - type supergiant models and discuss a prospect for mode identification . an attempt to interpret the oscillation spectrum of hd 163899 is presented in section [ hd163899 ] . the last section contains conclusions . the evolutionary models were calculated with mesa evolution code ( modules for experiments in stellar astrophysics , paxton et al . 2011 , paxton et al . 2013 ) . mesa allows to follow the evolution of massive stars to the pre - supernova phase . our previous calculations ( d2013 ) were limited to the phase before core helium ignition due to the lack of helium burning in warsaw - new jersey evolution code we used . we considered the mass range of 13@xmath318 @xmath0 . all computed models have an initial hydrogen abundance of @xmath4 and agss09 metal mixture ( asplund et al . the initial metal abundance varies from @xmath5 to @xmath6 . opal ( iglesias & rogers 1996 ) and op opacity tables ( seaton 2005 ) were used . all effects of rotation were neglected . non - adiabatic pulsation analysis was performed using the code of dziembowski ( 1977 ) . in fig.[fig1 ] we present evolutionary tracks including various effects : opacity , metallicity , overshooting , mass loss and element diffusion . they will be discussed in detail in the following subsections . the convective zones are determined by the ledoux criterion . it seems to be more justified than the schwarzschild criterion because of a direct sensitivity of the chemical gradient to evolution . lebreton et al . ( 2009 ) argued that the use of the ledoux criterion leads to a narrower icz and hence to the smaller number of unstable modes . in case of our models , mesa produces the wide icz with both criterions and this effect is not significant . for the envelope , we adopted the parameter of the mixing length theory ( mlt ) of @xmath7 which is the higher value than usually used for massive stars . during the main sequence ( ms ) and hydrogen - shell burning phases , the value of @xmath8 has a negligible effect on structure of massive stars but on the red giant branch ( rgb ) a higher efficiency of convection is needed and hence we opted for the higher value of @xmath8 . we took into account convective overshooting from the hydrogen and helium core and inward overshooting from non - burning convective zones , using the exponential formula ( herwig 2000 ) : @xmath9 where @xmath10 is the diffusion coefficient derived from mlt at a user - defined location near the boundary of a convective zone , @xmath11 is the pressure scale height at that location , @xmath12 is the distance in the radiative layer away from that location , and @xmath13 is an adjustable parameter , which we set to @xmath14 for most of calculated models . this value of @xmath13 mimics the behaviour of step overshooting with @xmath15 .    in our models the core helium ignition commences on the way towards the rgb but at this stage the energy produced from the helium burning is much lower than the energy produced from the hydrogen shell burning . the shell dominates the evolution up to the tip of rgb . at this point the energy produced in a helium core becomes comparable to the energy from the shell and the luminosity of a star drops . a star might eventually move towards higher effective temperatures and form a blue loop . the behaviour of the blue loops is still very poorly known and it is difficult to predict the effect of different parameters on their properties . we found that the convective overshooting has the most visible effect on their emergence , especially the inward overshooting from the non - burning zone . without that overshooting it is impossible to obtain blue loops in our models independently of metallicity , opacity , diffusion , etc . this effect is shown in fig.[fig2 ] , where we depicted two evolutionary tracks for the mass of 15 @xmath0 calculated with and without inward overshooting from the non - burning zone . the blue loop emerges if the @xmath16-gradient left by the evolution during ms is erased by the outer convective zone on the rgb . the range of this convective zone is increased by the inward overshooting . on the other hand , the overshooting from the hydrogen core generally acts against emerging of the blue loops . we show the effect of different values of the overshooting parameter on the evolutionary tracks in the central - left panel ( @xmath17 ) and in the bottom - left panel ( @xmath18 ) of fig.[fig1 ] . in both cases the blue loops are present but for the lower value of @xmath13 they have a larger extension . if we increase the overshooting only from the hydrogen core to @xmath18 and we leave @xmath17 for the helium core overshooting and the inward overshooting from non - burning zone , the blue loops still emerge but their extension is lower . this is caused mainly by the reduced overshooting from the helium core which decreases the amount of helium available to burn and hence leads to smaller extension of the blue loop .     with ( solid line ) and without ( dashed line ) inward overshooting from the non - burning convective zone . opal opacity tables and metallicity @xmath19 were used.,width=325 ]    the presented examples show that properties of mixing , especially the convective overshooting , have a huge impact on the behaviour of the blue loops . it is difficult to discriminate between different effects but generally it seems that the inward overshooting from the outer convective zone on the rgb is indispensable to obtain blue loops for the models in the studied range of masses . the overshooting from the hydrogen core suspends the emerging of the blue loops and the overshooting from the helium core has a slight effect on their extension .      for the models considered in section2.1 , we computed pulsational instability for modes with the degrees @xmath20 . they are shown as thick lines in in fig.[fig1 ] . the top panels show models calculated with @xmath21 and opal ( the top - left panel ) and op opacity tables ( the top - right panel ) . the central panels depicts models calculated for @xmath22 and opal tables without and with mass loss ( the central - left and central - right panels , respectively ) . the bottom - left panel shows the models calculated with @xmath22 , opal tables and convective overshooting parameter @xmath18 . the bottom - right panel depicts models calculated with element diffusion . details of the used mass - loss models and implementation of diffusion are described later in this section . similarly to our previous results ( d2013 ) there is the p- and g - mode instability beyond tams for the models that undergo hydrogen shell burning . most of the unstable modes are non - radial , but the radial fundamental mode is also unstable in hotter models beyond main sequence . the main difference between our current and all previous calculations is the presence of unstable non - radial modes on the blue loops for more massive models ( @xmath23 and @xmath24 ) . the radial modes are stable in this evolutionary stage in all calculated models . the existence of pulsation instability on the blue loop , as well as the blue loops themselves , depend critically on the metallicity , @xmath25 . with lowering the metallicity the blue loops reach the higher effective temperatures which has a huge impact on pulsation instability . in models with the lower value of @xmath26 the driving layer related to the metal opacity bump is located too deep in the star to be effective . on the other hand , the lower value of metallicity the lower the @xmath25 opacity bump . this crucially weakens instability for the ms and h - shell burning models but it seems to have a negligible effect on the instability on the blue loops , especially when compared to the effect of the effective temperature . in models with @xmath27 the instability in the phase before core helium ignition is entirely quenched whereas there are still a lot of unstable modes in the models with very well developed blue loop . the pulsational instability on the blue loops vanishes at @xmath28 . the reason for that is a difference in the internal structure resulting from a presence of an inward overshooting from the nonburning zone in the blue - loop models . the main effect of this overshooting is visible during the evolution on the rgb , but there is also a thin zone of overshooting below the icz . that leads to enrichment of hydrogen in layers just below the icz , which causes a steep increase of the @xmath16-gradient in a thin zone . as a consequence the brunt - visl frequency , @xmath29 , increases as well , what makes trapping of pulsation modes more effective and some of them can be driven by the @xmath1-mechanism operating in the z - bump despite low metallicity . the opposite situation occurs for models with higher metallicities . for example with @xmath6 there is a large instability area beyond ms but there are no unstable modes in the blue loop models , because of too low temperature the blue loops reach . the exceptions are the highest considered masses - @xmath30 for opal opacities and @xmath31 for op opacity tables . for those masses there are a few models on the hottest part of the blue loops that have unstable modes .    for every calculated mass and independently of metallicity or opacity tables we obtained unstable radial modes in the cooler models on the red giant branch or its vicinity . this result is similar to the one from saio et al . ( 2013 ) and it was expected because in this area of hr diagram a classical instability strip is located . those models are high - mass classical cepheids , in which the pulsations are driven in a region of the second ionization of helium . it is important to note that there are two effects that make the pulsational code we use unsuitable for such stars : the linear approximation and the lack of convection - pulsation interaction . we use the convective flux freezing approximation which is not adequate in large outer convective zones . we do not study properties of these stars to any further extent in this paper . ( solid lines ) and and the lamb frequencies , @xmath32 , for @xmath33 ( dashed lines ) and @xmath34 ( dotted lines ) , as a function of @xmath35 , for model1 ( top panel ) and model2 ( bottom panel ) . the scales of @xmath26 are different to represent the actual temperature profiles of the models.,title="fig:",width=332 ]   ( solid lines ) and and the lamb frequencies , @xmath32 , for @xmath33 ( dashed lines ) and @xmath34 ( dotted lines ) , as a function of @xmath35 , for model1 ( top panel ) and model2 ( bottom panel ) . the scales of @xmath26 are different to represent the actual temperature profiles of the models.,title="fig:",width=332 ]    as depicted in fig.[fig1 ] , for @xmath21 and @xmath36 there is a large instability strip that covers models both before and after core helium ignition . a comparison of the top- and central - left panels gives the effect of metallicity . for @xmath21 ( top - left panel ) , the range of effective temperatures that covers unstable models on the blue loop is a little wider than for @xmath22 ( central - left panel ) . this is a consequence of a larger extension of the loops towards the higher effective temperatures for lower metallicieties . as expected , an increase in @xmath25 results in wider instability area on ms and beyond tams . the upper panels show an effect of opacity tables on instability domains for @xmath21 , but the main conclusions are the same for every tested metallicity . there are significant differences on the ms phase , due to a long - known fact that for op tables instability starts earlier than for opal tables . during the hydrogen shell burning phase and on the blue loop the differences are very subtle . the loops and instability areas have slightly different shapes but generally they are similar for both sources of the opacity data .    in the central panels we compare instability domains for models with @xmath22 calculated without ( the central - left panel ) and with mass loss ( the central - right panel ) . for @xmath37 we used the wind s scheme of vink et al . ( 2001 ) whereas for @xmath38 a recipe of de jager et al . ( 1988 ) was adopted . we do not study very massive models and hence the amount of mass which is lost during the evolution is rather modest . for instance , a star with initial mass of @xmath39 loses only about @xmath40 during the evolution from zams to the point when helium is depleted in its core . the maximum value of the mass - loss rate is reached on the tip of the red giant branch and for the @xmath39 star it has the value of @xmath41 . mass loss seems to have a negligible effect on instability of models immediately beyond tams , whereas this effect on the instability on the blue loops is much more pronounced . with the reduced mass a star reaches lower luminosity during helium burning phase and also a blue loop reaches lower effective temperature than for models calculated without mass loss . both of these effects lead to smaller instability range during this phase of evolution . the description of wind which we used for lower effective temperatures ( de jager et al . ( 1988 ) ) gives relatively low values of mass loss when compared to other models we checked ( nieuwenhuijzen & de jager 1990 , van loon et al . 2005 , reimers 1975 ) . it means that the instability areas on the blue loops calculated with other mass - loss models would be even smaller or entirely vanished . also with sufficiently high mass loss the blue loops do not emerge at all . we will not discuss effects of mass loss anymore in this paper because they do not change any qualitative results presented here . despite this it has to be remembered that it is a common phenomenon and it should be enabled during construction of detailed asteroseismic models . the bottom - left panel depicts the instability areas calculated for higher value of convective overshooting , @xmath18 . the most important difference when compared to models with @xmath17 ( central - left panel ) is smaller extension of the blue loops and hence there are no unstable models present during this evolutionary phase .    in the bottom - right panel of fig.[fig1 ] we present models calculated with @xmath22 , opal opacity tables and element diffusion . mesa calculates element diffusion by solving burgler s equations using the method and diffusion coefficients of thoul et al . the effect of diffusion on the shape of evolution tracks and instability areas is negligible , both for models before and after core helium ignition . the small changes of abundances of chemical elements in the central regions of the star do not modify significantly a behaviour of pulsation modes . the minor differences of frequencies might be important during a comparison to the observations but they do not influence the qualitative results presented here . here , we compare pulsation properties of two representative models with similar positions in the h - r diagram but in two different evolutionary stages , one during the shell hydrogen burning phase ( @xmath39 , @xmath42 , @xmath43 , hereafter model1 ) and the second during core helium burning ( @xmath44 , @xmath45 , @xmath46 , hereafter model2 ) . they are calculated with opal opacity tables , agss09 metal mixture and @xmath19 . propagation of pulsational modes in the radial direction is confined to the regions determined by the two characteristic frequencies ( e.g. , unno et al . 1989 ) , i.e. , the lamb ( acoustic ) frequency ( @xmath47 ) and brunt - visl ( buoyancy ) frequency ( @xmath48 ) . the values of @xmath47 depend mostly on the mode degree , @xmath49 , whereas the brunt - visl frequency is very sensitive to the structure of a star ( e.g. , a presence of convective zones or chemical gradients ) .    in fig.[fig3 ] , we present propagation diagrams , i.e. , @xmath50 , @xmath51 , for models 1 and 2 ( the top and bottom panels , respectively ) . the main difference between these models is the presence of a convective core ( @xmath52 ) in the model on the blue loop , whereas the model that undergoes shell hydrogen burning has a radiative core . beyond the core , the propagation diagrams of these two models are qualitatively similar ; in particular , both have a fully developed intermediate convective zone . the two spikes in the run of brunt - visl frequency that occur in the vicinity of convective zones in model2 have different origin . the one at the border of the convective core is both a consequence of the very steep @xmath16-gradient in this region of star ( the main reason ) and also a numerical effect related to the extreme difficulties with precise determination of the border of helium core by evolution codes . due to a very high sensitivity of triple - alpha process to a change of temperature , the location of border changes slightly with every step of evolution . that results in an artificial increase of the @xmath16-gradient , which is clearly visible in behaviour of @xmath29 . the second spike , at the innermost border of icz is related to the presence of inward overshooting from this convective zone . that leads to a steep change of @xmath16 and a peak in the propagation diagram .    in fig.[fig4 ] , we present the instability parameter , @xmath53 , for modes with the harmonic spherical degrees @xmath54 , as a function of frequency for models 1 and 2 ( the top and bottom panels respectively ) . instability parameter measures the net energy gained by a mode during one pulsational cycle and is defined as ( stellingwerf 1978 ) :    @xmath55    where @xmath56 is the global work integral . this definition gives normalization @xmath57 $ ] . if @xmath58 a pulsation mode is unstable and if @xmath59 it is stable . there are two global maxima of @xmath53 in model1 , which is similar to the behaviour during the ms evolution ( d2013 ) . the first one , at low frequencies , is related to the very high - order g - modes , while the second one , at higher frequencies , is related to the low - order p- and g - modes . all non - radial modes in the second group are mixed modes and they behave like p - modes only in the outer envelope of the star . in the frequency range presented in fig.[fig4 ] , i.e. , from 0.04 to 3.5 d@xmath60 , the values of the radial order , @xmath61 , are in ranges of about @xmath62 for the @xmath33 modes and @xmath63 for the @xmath34 modes . in model1 the fundamental radial mode and the first overtone are unstable with the ratio of periods @xmath64 . the most interesting behaviour of @xmath53 is visible in the vicinity of the first maximum . there is a clearly visible pattern of consecutive local minima and maxima , which does not occur in ms models . this phenomenon is related to the partial trapping of pulsation modes in the radiative core and the envelope . the detailed description of this phenomenon can be found in our previous paper ( d2013 ) .    in model2 , which burns helium in its core , a situation is different . there are only a few very high - order low - frequency g - modes which are unstable . these maximum values of @xmath53 are shifted towards lower frequencies when compared to model1 . there is no pattern of local minima and maxima and the calculated frequency spectrum is much sparser than for model1 . this is typical for every blue - loop that a number of unstable modes is very low . model2 has one unstable dipole mode with frequencies of @xmath65 c / d and two unstable quadruple modes with frequencies @xmath66 c / d . radial modes are always stable in models during core helium burning phase , regardless of metallicity , opacity tables , etc . however , the ratio of periods between the first overtone and the fundamental mode is visibly lower than in model1 and has a value @xmath67 . in model2 , the radial orders of pulsation modes are much larger , in particular for the quadruple mode . for the considered range of @xmath68 , cf . fig.[fig4 ] , the values of @xmath61 are : @xmath69 for @xmath33 and @xmath70 for @xmath34 .    in fig.[fig5 ] , we compare the kinetic energy densities , @xmath71 ( top panels ) , and differential work integrals , @xmath72 ( bottom panels ) , for unstable and stable modes of the degree @xmath34 with very close frequencies . the left group of panels is for model1 ( before core helium ignition ) and the left group of panels for model2 ( after core helium ignition ) . the behaviour of presented modes is typical and representative for all unstable and stable modes during shell hydrogen burning and on the blue loops . the properties of pulsational modes in the shell hydrogen burning model have been already presented in d2013 but for comparative purposes we repeat them here .    in model1 , both unstable ( @xmath73 c / d , @xmath74 ) and stable ( @xmath75 , @xmath76 ) modes exhibit a very strong damping in the radiative helium core , which is visible in the behaviour of @xmath72 in the bottom panels of the left group . as shown in the upper panels of this group , most of the kinetic energy of the modes is confined in the central region . the most notable difference between the presented modes is that the unstable one has some energy deposited also in the radiative envelope , whereas the stable one has almost no energy in this region . this is due to the fact that the unstable mode has been partially reflected at the icz and is trapped in two cavities : the core and the outer layers . that energy , although its amount is small when compared to the amount trapped in the core , is sufficient for the @xmath1 mechanism , operating in the @xmath2bump , to efficiently drive the pulsations . this is visible in the behaviour of differential work integral - there is a maximum around @xmath77 and it dominates over the damping in the core . on the contrary , the stable mode has not been reflected at icz because it does not have a node of eigenfunction there . that means that it penetrates the core with high amplitude and is effectively damped there . the above behaviour of @xmath71 explains the pattern of local maxima and minima of @xmath53 described earlier on .    for model2 we also compare two very close modes : the unstable one with @xmath78 c / d and @xmath79 and the stable one with @xmath80 c / d and @xmath81 . their kinetic energy densities and differential work integrals are shown in the right group of panels in fig.[fig5 ] . the behavior of @xmath71 is different for these modes than for the modes from the model1 . the entire energy of the stable mode is confined to the radiative zone in - between the convective core and the icz ( there is strong damping in this region of the star ) , whereas the unstable mode has almost all of its energy trapped in the outer radiative zone . that means that the unstable modes of the blue loop models have to be almost entirely reflected at the icz . this fact might be the explanation why we do have far fewer unstable modes from models burning helium in its core than from models before core helium ignition and also the lack of the structure of @xmath53 observed in model1 . the behaviour of kinetic energy density depicted in fig.[fig5 ] is typical for all unstable and stable modes both during hydrogen shell burning and on the blue loops . in our previous paper ( d2013 ) , we have shown that one should not expect regular patterns in the oscillation spectra of blue supergiants which the asymptotic theory predicts for g - mode pulsations . however , it turned out that there is a prospect for identification of the degree , @xmath49 , for modes excited in the blue supergiant models before core helium ignition from photometric diagnostic diagrams ( d2013 ) . here , we will check whether the photometric amplitudes and phases of pulsations in models after core helium ignition have similar properties . the expression for the complex amplitudes of the light variations in a given passband within the linear and zero - rotation approximation can be found , e.g. , in daszyska - daszkiewicz et al . ( 2002 ) , and need not be repeated here . to compute the values of these observables , we used the nonadiabatic pulsational code of dziembowski ( 1977 ) and tables of the fluxes in photometric passbands and nonlinear limb - darkening coefficients computed by daszyska - daszkiewicz & szewczuk ( 2011 ) on the basis of the nlte atmosphere models of lanz & hubeny ( 2007 ) .    in fig.[fig6 ] , we show positions of unstable modes with the degree @xmath82 excited in the supergiant models in the photometric diagnostic diagrams employing the @xmath83 and @xmath84 johnson filters . these passbands are similar to those used on the bright - star target explorer constellation ( brite ) which is already operating . we used the opal data for the agss09 mixture and @xmath22 . in the left panel we marked positions of all unstable modes in models with masses @xmath85 whereas in the right panel we chose modes excited in models with @xmath86 and @xmath87 . these temperature range is a typical maximum error and larger than the error for hd 163899 ( saio et al . 2006 ) . with different symbols we discriminate between different values of @xmath49 and evolutionary stages . their meanings are described in the legend of the plots .    in the diagram for all studied masses , there is a perfect separation between radial and non - radial modes and fairly good separation between dipoles and quadruples . a discrimination between models undergoing shell hydrogen burning and models on the blue loops is not obvious . modes excited in models on these evolution stages tend to group in slightly different locations but in a plot including such a range of masses their distinction would be uncertain given observational errors in the photometric amplitudes and phases . the situation improves slightly when we are able to constrain the mass and effective temperature . after fixing the range of @xmath26 in the right panel of fig.[fig6 ] , the separation between phases before and after core helium ignition can be noticeable . fig.[fig7 ] presents how the photometric observables depend on the mode frequency considering again the br johnson passbands and the same models as in the right panel of fig.[fig6 ] . from this plot we can conclude that much better separation between modes of different values of @xmath49 , from both the amplitude ratios and phase differences , occurs for higher frequency modes which are unstable only in model1 . and @xmath88 , respectively . the models have been calculated with @xmath21 and opal opacity tables . the grey area depicts the error of @xmath26 for hd 163899 ( saio et al . 2006).,width=332 ]    the greatest uncertainty concerning spbsg stars is their evolutionary status . although the basic parameters ( @xmath26 , @xmath89 , @xmath25 , @xmath90 , etc . ) of hd 163899 are very poorly determined or unknown , our calculations still can give some insights about the star . in fig.[fig8 ] , we present a comparison of evolution times in the pulsational unstable region of hr diagram during shell hydrogen burning phase and on the blue loop , for models with masses @xmath91 and @xmath92 , respectively , and @xmath21 . the star spends only a very short time ( about 8 000 years ) in the instability strip during the shell - hydrogen burning phase compared to the time spent on the unstable region of the blue loop ( about 800 000 years ) . there should be much more stars in the region of hr diagram occupied by blue supergiants that undergo core helium burning than stars that moves towards red giant branch ( rgb ) . on the one hand , this is the main and only argument why spbsg stars should be on the blue loop . on the other hand , hd 163899 is the only know star of this type and we might be very lucky to catch it during the very short phase of pulsational instability before core helium ignition . other arguments also prefer this evolutionary stage . the instability area during a crossing towards rgb is wider in @xmath93 than on the blue loop . this is especially important when we take into account the effect of mass loss which is ubiquitous effect for massive stars . in such a case the instability area on the blue loops become very small . the most important argument is a number of excited modes ; there is much more unstable modes in models during the shell - hydrogen burning phase than in models on the blue loop . the last argument against the evolution on the blue loop is a complete lack of unstable high - frequency modes in this phase . such modes are present in the spectrum of hd 163899 and they are present in our models before core helium ignition .    in fig.[fig9 ] , we compare the oscillation spectrum of hd 163899 ( top panel ) and the theoretical one computed for a model in the shell - hydrogen burning phase ( @xmath53 vs. @xmath68 , bottom panel ) . pulsational modes with the degree up to @xmath94 were considered . this is one of models in which the larger number of pulsation modes are excited . the observed spectrum consists of 48 frequencies between @xmath95 to @xmath96 c / d ( c.f . table 1 in saio et al . this range covers both g - modes ( high - order ) and p - modes in post - ms models . there are also frequencies between those two groups which do not occur in pulsation calculations . we do not know @xmath90 for the star and hence we do not include rotation in our calculations . this is left for the future work . the estimated critical velocity for model1 , which has the stellar radius @xmath97 , is @xmath98 km / s which corresponds to the critical rotation frequency of about @xmath99 c / d . it is very unlikely that a supergiant would rotate at such a high rate , but even with lower rotation rate the very low dominant frequency , @xmath100 c / d , could be explained by retrograde modes . the rotational splitting could also solve the biggest problem of our pulsation models which is too low number of unstable g - modes . as for the observed range of frequencies which occur between the theoretical g- and p - mode maxima , it could be filled in with the high @xmath49 modes of various azimuthal orders , @xmath101 . our goal was to study pulsational properties of b - type supergiants with masses @xmath102 , especially to check the possibility whether spbsg stars can undergo core helium burning . in addition to previously known spbsg instability strip ( saio et al . 2006 , d2013 ) which covers the stars that move towards rgb after departure from ms , we found the instability strip for the blue - loop models for the temperatures of b spectral types . its width depends on the extension and shape of the blue loops and hence it varies with many parameters , especially with metallicity , overshooting and mass loss . one of the most important results of our paper is the fact that an inward overshooting from the outer convective zone on the rgb phase is indispensable for the emergence of blue loops at all . the lower values of metallicity in the main sequence and shell hydrogen burning models acts towards decreasing pulsational instability whereas in the case of the blue loops models the most important factor is the effective temperature . for lower values of @xmath25 the blue loops have larger extension and the driving layer is located at the proper place for pulsation to be excited . we found unstable modes in the blue loop models with metallicity as low as @xmath103 . there is a difference in behaviour of instability parameter , @xmath53 , between models before and after core helium ignition . as we previously shown ( d2013 ) , in models during shell - hydrogen burning phase there are two groups of unstable mods ; one related to the very high - order g - modes and the second related to the p - mode behaviour . on the contrary , in models on the blue loops only the very high - order g - modes are unstable and the number of unstable modes is much lower than in models that departed from tams and move towards rgb . that might be explained by the difference in the kinetic energy density distribution between those two evolutionary phases . in models that undergo shell hydrogen burning , only a partial reflection of pulsation mode at the icz is sufficient to efficiently drive the mode , whereas on the blue loop the mode has to be almost entirely reflected in order to be unstable .    comparing the properties of models before and after core helium ignition and frequencies detected in the light curve of hd 163899 , we drew a conclusion that spbsg stars should rather undergo shell hydrogen burning and not core helium burning on the blue loop . however , to fully support this statement a better determination of basic parameters of hd 163899 is needed and , in the next step , seismic model of the stars should be constructed taking into account the effects of rotation . we found also that mode identification from multicolour photometry is applicable to spbsg stars and hence multi - colour photometry for hd 163899 would be desirable . the analysis of such data could finally lead to constraining the evolutionary status of the star and help to understand pulsations of spbsg stars . this work was financially supported by the polish ncn grants 2013/09/n / st9/00611 , 2011/01/m / st9/05914 and 2011/01/b / st9/05448 . 99 asplund , m. , grevesse , n. , sauval , a. j. , scott , p. : 2009 , ara&a 47 , 481 daszyska - daszkiewicz j. , dziembowski w. a. , pamyatnykh a. a. , goupil m - j . , 2002 , a&a , 392 , 151 daszyska - daszkiewicz , j. , ostrowski , j. , pamyatnykh a.a . : 2013 , mnras 432 , 3153 ( d2013 ) daszyska - daszkiewicz j. , szewczuk w. , 2011 , apj , 728 , 17 dziembowski , w. : 1977 , aca 27 , 95 godart , m. , noels , a. , dupret , m .- a . , lebreton , y. , 2009 , mnras 396 , 1833 herwig , f. : 2000 , a&a , 360 , 952 iglesias , c. a. , rogers , f. j. : 1996 , apj 464 , 943 de jager , c. , nieuwenhuijzen , h. , van der hucht , k. a. : 1988 , a&as 72 , 259 klare , g. , neckel , t. : 1977 , a&as 27 , 215 lanz t. , hubeny i. 2007 , apjs , 169 , 83 lebreton , y. , montalbn , j. , godart , m. , et al . : 2009 , coast 158 , 277 van loon , j. th . , cioni , m .- l. , zijlstra , a. a. , loup , c. : 2005 , a&a 438 , 273 nieuwenhuijzen , h. , de jager , c. : 1990 , a&a 231 , 134 pamyatnykh , a. a. : 1999 , aca 49 , 119 pamyatnykh , a. a. , ziomek , w. : 2007 , coast 150 , 207 paxton , b. , bildsten , l. , dotter , a. , et al . : 2011 , apjs 192 , 3 paxton , b. , cantiello , m. , arras , p. , et al . : 2013 , apjs 208 , 4 reimers , d. : 1975 , mem . liege 8 , 369 saio , h. , kuschnig , r. , gautschy , a. , et al . : 2006 , apj 650 , 1111 saio , h. , georgy , c. , meynet , g. : 2013 , mnras 433 , 1246 schmidt , e. g. , carruthers , g. r. : 1996 , apjs 104 , 101 seaton , m. j. : 2005 , mnras 362 , l1 stellingwerf r. f. : 1978 , aj 83 , 1184 thoul , a. a. , bahcall , j. n. , loeb , a. : 1994 , apj 421 , 828 unno , w. , osaki , y. , ando , h. , saio , h. , shibahashi , h. : 1989 , nonradial oscillations of stars , tokyo , university of tokyo press vink , j. s. , de koter , a. , lamers , h. j. g. l. m. : 2001 , a&a 369 , 574
Abstract:

Answer: evolutionary tracks and pulsational analysis of models with masses of 13 - 18 @xmath0 are presented . we address two important questions . the first one deals with one of the most unresolved problems in astrophysics , i.e. , the existence of a blue loop after core helium ignition ; the so called `` to loop or not to loop '' problem . we show that inward overshooting from the outer convective zone in the red giant phase is prerequisite for the development of the blue loop . our second question concerns pulsational instability of models in the core helium burning phase . we present for the first time that models on the blue loop can have unstable modes driven by the @xmath1 mechanism operating in the @xmath2bump . contrary to post - main sequence models in the shell hydrogen burning phases , pulsational instability of the blue loop models depends mainly on effective temperature and metallicity is of secondary importance . finally , we try to interpret the oscillation spectrum of the blue supergiant hd 163899 , the only member of the spbsg class , and to get some clue on the evolutionary status of the star . [ firstpage ]    stars : early - type  stars : supergiants  stars : oscillations