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: given a discrete distribution @xmath0 on @xmath1 classes @xmath2 , we consider an index called _ polarization measure _ , defined by @xmath3    the polarization measure has been introduced in economics for real distributions by . the discrete version we consider here has been used in . the polarization measure has the following interpretation . let @xmath4 be i.i.d . @xmath5 and consider the indicator of _ exactly two equal _ @xmath6 then @xmath7 } =   3 \sum_{x=0}^n \pi_x^2(1-\pi_x ) = 3 { \operatorname{pol}\left(\pi\right)}$ ] . [ fig : polarization ]    the polarization measure on the classes @xmath8 , as shown in fig . [ fig : polarization ] , has an unstable critical point at the uniform distribution , it is zero in the case of concentration in one class , and reaches its maximum 1/4 on distributions on two classes with equal probabilities . polarization measure was devised to be an index of the distance of a distribution from the three cases of maximal polarization . in fig . [ fig : polarization ] the simplex is represented as an equilateral triangle . in the following we shall use different sets of coordinates to represent the probability simplex , e.g. see fig . [ fig : gradpol - a ] ( left ) . we want to study the dynamics of this index , i.e. to characterise evolutions that maximise or minimise the index . this study requires tools from information geometry ( ig ) e.g. , . however , the following presentation is actually largely self - contained . the recourse to ig is not dispensable because the ordinary gradient flow of @xmath9 , as shown in fig . [ fig : gradpol - b ] ( left ) , does not lead to the extrema of interest on the border of the probability simplex . consequently , one wants to turn to a different way to compute the gradient , i.e. to the so - called amari s natural gradient . we use elementary fact of the theory of dynamical systems to characterise critical points of the gradient flow and refer to . the basics of ig are discussed in sec . [ sec : natural gradient ] . the application of ig to the polarization measure is described in sec . [ sec : pol ] . the possibility of a generalisation of such an index is shortly discussed in sec . [ sec : pols ] , while the reduction of the problem to the study of an exponential family is presented in sec . [ sec : polase ] . we suggest a possible application in sec . [ sec : app ] . further material , not directly related with the measure of polarization , but suggested by the methodology , is presented in the appendixes . differential equations on the probability simplex are well known in applications other then descriptive statistics . we briefly discuss the relations between these applications and our one in app . [ sec : replicator ] . in app . [ sec : hessian ] some issues related to the second order calculus are briefly discussed . we denote by @xmath10 the simplex of the probability function @xmath0 on @xmath11 . the interior of the simplex , @xmath12 , is the set of the strictly positive probability functions , @xmath13    the border of the simplex is the union of all the faces of @xmath10 as a convex set . we recall that a face of maximal dimension @xmath14 is called facet . a facet is a simplex of dimension @xmath14 . we define @xmath15 to be the vector space of random variables @xmath16 that are @xmath0-centered , @xmath17 } = 0 $ ] . in the geometry of @xmath18 , @xmath15 is the plane through the origin , orthogonal to the vector @xmath19 .    1 . the _ tangent bundle _ of the open simplex @xmath20 is the set @xmath21 2 .   if @xmath22 is a one - dimensional statistical model , geometrically a curve , its _ score _ @xmath23 belongs to @xmath24 for all @xmath25 . as the score is a centered random variable , hence @xmath26 is a curve in the tangent bundle . in fact , @xmath27 is meant to represent a generic velocity vector through @xmath0 , see fig . [ fig : tangentbundle ] . the score is a representation of the velocity along a curve , because of a geometric interpretation of c. r. rao s classical computation : @xmath28= \frac d { dt } \sum_x u(x)p(x;t)= \sum_x u(x )   \frac d { dt } p(x;t)=\\ \sum_x u(x )   \frac d { dt } \log\left ( p(x;t)\right)p(x;t)= \sum_x \left(u(x)-e_t[u]\right )   \frac d { dt } \log\left ( p(x;t)\right)p(x;t)=\\ e_t\left[\left(u - e_t[u]\right ) \frac d { dt } \log\left ( p(t)\right ) \right]=   { \left\langleu - e_t[u],d p(t)\right\rangle_{p(t)}}\end{gathered}\ ] ] we observe that the scalar product above is the scalar product on @xmath24 . a curve on the simplex is a parametric model . the probability @xmath0 is represented by a vector from @xmath29 to the point whose coordinates are @xmath30 . in fig . [ fig : tangentbundle ] , the velocity vectors are represented by arrows ; they are orthogonal to the vectors @xmath19 . [ fig : tangentbundle ]    in our context , a _ vector field _ @xmath31 is a mapping such that @xmath32 , i.e. such that the couple @xmath33 belongs to the tangent bundle . @xmath34 for all @xmath0 . because of our geometrical construction , here we prefer to call @xmath35 a vector field , but we want to stress its statistical meaning of centered function of the distribution . a _ differential equation _ is an equation of the form @xmath36 . given a real function @xmath37 , its _ gradient _ is the vector field @xmath38 such that for all curves @xmath39 we have @xmath40 } , u , v \in b_\pi.\ ] ]    the rao s computation in eq . is the prototypical gradient computation . the _ gradient flow equation _ is the differential equation @xmath41    along a solution of the gradient flow equation the value of @xmath42 is increasing because @xmath43 . actually the solution of the gradient flow equation is the curve of _ steepest ascent_.    computations are usually performed in a _ @xmath44 @xmath45 being an open set in @xmath46 . the @xmath47-th coordinate curve is obtained by fixing the other @xmath14 components and moving @xmath48 only . the scores of the @xmath47-th coordinate curves are the random variables @xmath49    the sequence @xmath50 is a vector basis of the tangent space @xmath51 . the representation of the scalar product in such a basis is @xmath52 where the matrix @xmath53_{i , i=1}^n$ ] is the _ fisher information _ matrix . if @xmath54 is the expression in the parameters of a function @xmath55 , that is @xmath56 , and @xmath57 is the expression in the parameters of a generic curve @xmath58 , then the components of the gradient in are expressed in terms of the ordinary gradient by observing that @xmath59    as @xmath60 , we obtain from @xmath61    the _ natural gradient _ is a vector @xmath62 whose components are the coordinates of the gradient @xmath63 in its @xmath0-basis , that is @xmath64    by substitution of the expression in we obtain @xmath65    fig . [ fig : albero ] is an illustration for the function @xmath66 . the _ common parametrization _ of the ( flat ) simplex @xmath67 is the projection on the solid simplex @xmath68 , that is @xmath69 in which case @xmath70 , @xmath71 , is the random variable with values @xmath72 at @xmath73 , 1 at @xmath74 , 0 otherwise , hence @xmath75 and @xmath76    the element @xmath77 of the fisher information matrix is @xmath78 } = \\ \sum_x \pi(x,\bm\eta)^{-1}\left((x = j)(j = h)+(x = 0)\right ) = \eta_j^{-1}(j = h ) + \left(1 - \sum_k \eta_k\right)^{-1}\end{gathered}\ ] ] hence @xmath79_{i , j=1}^n.\ ] ]    as an example we consider @xmath80 . the fisher information matrix , its inverse and the determinant of the inverse are , respectively , @xmath81 @xmath82 @xmath83 note that the computation of the inverse of @xmath84 is an application of the sherman - morrison formula and the computation of the determinant of @xmath85 is an application of the matrix determinant lemma .    for general @xmath86 , we have the following proposition , whose interest stems from the definition of natural gradient , see eq . . [ prop : i-1 ]    1 . the inverse of the fisher information matrix is @xmath87 2 .   in particular , @xmath85 is zero on the vertexes of the simplex , only . the determinant of the fisher information matrix is @xmath88 4 . the determinant of @xmath85 is zero on the borders of the simplex , only . 5 .   on the interior of each facet , the rank of @xmath85 is @xmath14 and the @xmath14 liner independent column vectors generate the subspace parallel to the facet itself .    1 .   by direct computation , @xmath89 is the identity matrix . 2 .   the diagonal elements of @xmath85 are zero if @xmath90 or @xmath91 , for @xmath71 . if , for a given @xmath47 , @xmath90 , then the elements of @xmath85 are zero if @xmath92 , @xmath93 . the remaining case corresponds to @xmath91 for all @xmath47 . then @xmath94 on all the vertexes of the simplex . 3 .   it follows from matrix determinant lemma . the determinant factors in terms corresponding to the equations of the facets . 5 .   given @xmath95 , the conditions @xmath96 and @xmath97 for all @xmath98 , define the interior of the facet orthogonal to standard base vector @xmath99 . in this case the @xmath95-th row and the @xmath95-th column of @xmath100 are zero and the complement matrix corresponds to the inverse of a fisher information matrix in dimension @xmath14 with non zero determinant . it follows that the subspace generated by the columns has dimension @xmath14 and coincides with the space orthogonal to @xmath101 . consider the facet defined by @xmath102 , @xmath103 for all @xmath95 . for a given @xmath47 , the matrix without the @xmath47-th row and the @xmath47-th column has determinant @xmath104 . on the considered facet this determinant is different to zero and @xmath85 has rank @xmath14 and their columns are orthogonal to the constant vector . an other parametrization is the _ exponential parametrization _ based on the exponential family with sufficient statistics @xmath105 , @xmath71 , @xmath106 where @xmath107 some of the properties discussed in prop . [ prop : i-1 ] should actually be discussed under the exponential parametrization , see e.g. , but we do not do that here . we will discuss the exponential parametrization below in sec . [ sec : polase ] to show that polarization can be seen as an expectation with respect to an exponential family . we apply now the general theory of the natural gradient to the study of the dynamics of the polarization measure . our goal is to find the lines of the steepest ascent of the function @xmath9 .    in the common parametrization we have @xmath108 and , for @xmath109 , @xmath110 with gradient @xmath111 the inverse of the fisher information matrix is @xmath112 the natural gradient is @xmath113 see in fig.s [ fig : gradpol - a ] and [ fig : gradpol - b ] the gradient fields .    [ cols="^,^ " , ]     to study the flow in the fixed points we consider the sign of the eigenvalues of the jacobian of the natural gradient , calculated in the fixed points , see @xcite .    the jacobian of @xmath114 is @xmath115 $ ] where @xmath116 @xmath117    the jacobian calculated in the vertexes are @xmath118 and the two eigenvalues of the three jacobian are both positive . then , the fixed points on the vertexes repel flow locally . moreover , the natural gradient is zero on the midpoints of the borders , @xmath119 , @xmath120 and @xmath121 . the jacobian calculated in the midpoints are @xmath122 and the two eigenvalues of the three jacobian are both negative . then , these fixed points attract flow locally .    in standard cases the jacobian of a gradient is the hessian matrix . this is not true anymore in ig where we compute the jacobian of the natural gradient . actually , in ig there are various notions of hessian , each one based on a different connection on the tangent bundle . we consider now an inverse problem in the simple case @xmath109 . we want a third degree symmetric polynomial which could be used as polarization measure . computations were performed using sage @xcite . a generic polynomial is @xmath123 and its expression in the parameters is @xmath124    the components of the natural gradient are @xmath125 note that the case @xmath126 is @xmath127 and @xmath128 . because of the multiplication by @xmath129 , the natural gradient @xmath130 is zero in each of the simplex vertexes @xmath131 , @xmath132 , @xmath133 . because of the symmetry , the natural gradient @xmath130 is zero in each of the mid - point of the edges , i.e. @xmath134 , @xmath135 , @xmath136 and at the uniform probability , @xmath137 . we need to turn to the discussion of the jacobian matrix of the natural gradient of @xmath138 .    at the vertexes , the jacobian is @xmath139 so that we want @xmath140 .    at the uniform probability @xmath137 the jacobian is @xmath141 a _ necessary _ condition of non definiteness is @xmath142    at the mid - points of the three edges @xmath134 , @xmath135 , @xmath136 , the values of the jacobian matrices are , respectively , @xmath143 and , imposing the necessary condition @xmath144 of eq . , they are @xmath145 the two eigenvalues of each of these matrices are both negative if @xmath146 . notice that this condition of attracting flow on the midpoints is the same condition of repelling flow on the vertexes . summarising , the conditions on the coefficients of a third degree symmetric polynomial that represent a measure with the properties stated above , are @xmath147    in conclusion , it is possible to design other measures of polarization in the form of a symmetric polynomial of degree 3 , or more general forms . in particular , it would be interesting to have a measure in a form similar to the entropy . the analogy is suggested by the fact that the gradient flow of the entropy gives trajectories that move from the uniform probability to one of the vertexes of the simplex . we have already observed that the function @xmath9 defined in eq . [ eq : pol ] is an homogeneous polynomial of degree 3 in the indeterminates @xmath148 . the general class of indexes based on homogeneous polynomials that we have discussed in sec . [ sec : pols ] is of special interest because they reduce to an expectation with respect to an exponential family , as we discuss now in the case of three sample points .    if the random variables @xmath4 are i.i.d . with the distribution of @xmath149 supported by @xmath8 , and coded as an exponential family i.e. @xmath150 then the joint probability function of @xmath4 is itself an exponential family , @xmath151 with sufficient statistics @xmath152 , @xmath153 . the model for @xmath154 is @xmath155 where @xmath138 is the table count of @xmath154 . the random variables of this model are represented in tab . [ tab:1 ] . @xmath156    the marginal polytope , that is the convex set generated by the values of the sufficient statistics , is illustrated in fig . [ fig : mpolytope ] . we refer to @xcite for the relevant theory . m.65m.30 & @xmath157    if @xmath158 is the indicator function of the set @xmath159 then the value of @xmath160 , when expressed in the parameters @xmath161 , is @xmath162}$ ] . thus , the problem of the maximisation of the polarization measure is rephrased to the problem of finding the maximum of the expected value of a random variable on a given exponential family . this approach has been considered in combinatorial optimisation , see e.g. . it has a number of issues . first , any convergent evolution along the exponential family has either a limit internal to the exponential family itself of a limit outside the exponential family , supported by a face of the marginal polytope , see . second , the maximum of the function can be reached as a limit of expected values if , and only if , there exists , among the distributions obtained conditioning the exponential family to one face of the marginal polytope , a distribution such that the expected value of the function is equal to the maximum of the function itself . we do not enter here in a detailed discussion , see more information in . the expectation parameters of the exponential family are @xmath163 } & = 3 \frac{\partial}{\partial\theta_1 } \psi(\theta_1,\theta_2 ) = 3 \frac { { \mathrm{e}}^{\theta_1}}{1 + { \mathrm{e}}^{\theta_1 } + { \mathrm{e}}^{\theta_2 } } , \\    { { \mathbb e}_{\bm \theta}\left[t_2\right]}&= 3 \frac{\partial}{\partial\theta_2 } \psi(\theta_1,\theta_2 ) = 3 \frac { { \mathrm{e}}^{\theta_2}}{1 + { \mathrm{e}}^{\theta_1 } + { \mathrm{e}}^{\theta_2}}.\end{aligned}\ ] ] these parameters are related to the @xmath164 by @xmath165 } = { { \mathbb e}_{\bm\theta}\left[x_j\right ] } + { { \mathbb e}_{\bm\theta}\left[y_j\right ] } + { { \mathbb e}_{\bm\theta}\left[z_j\right ] } = 3\eta_j , \quad j = 1,2.\ ] ]    the inverse of the map @xmath166 can be computed explicitly as @xmath167    the cumulant function is expressed as a function of @xmath164 as @xmath168 and the exponential family in eq . is expressed in the parameter @xmath169 as @xmath170    when the probabilities are expressed in the form of eq . , we can actually compute values for border cases @xcite . we recover in a different way the result already known . @xmath171 : : :    eq . becomes    @xmath172 ,    which is zero but for @xmath173 . this distribution is    concentrated on    @xmath174 , with    @xmath175 . it follows    @xmath176 } = 3/4 $ ] . @xmath177 : : :    same as the previous case . the distribution is concentrated where    @xmath178 , @xmath179 , and    @xmath180 } = 3/4 $ ] . @xmath181 : : :    in this case the distribution is    @xmath182 with    support on the face where @xmath183 , probabilities    @xmath184 , and    @xmath185 } = 3/4 $ ] . we have shown on an example that the maximal value of the polarization measure is actually reachable as a limit of the expected value of a random variable @xmath186 on the exponential family , but this does not mean that the maximum value of the expected value is equal to the maximum value of the random variable itself . this would be true if the maximum of the random variable were reached on a face of the marginal polytope . this is discussed in the literature we have cited above . we have considered a statistical index different from indexes for concentration or uniformity , such as the discrete gini index or boltzmann - gibbs - shannon entropy . the polarization measure has been discussed in a dynamic way , by considering its variation and computing the directions of steepest variation . the study requires tools suitable to discuss differential equation on a differentiable manifold . this methodology suggests to implement the velocity of variation itself as a statistical index . consider a study of the evolution of an index in time such as . in the time series @xmath187 the evolution of the index e.g. , @xmath188 , could be misleading , because an increase in the index could be associated to a shift from a basin of attraction to a different basin of attraction . we suggest a more precise local study as follows . given a movement from @xmath189 to @xmath190 , we look for a comparison of an estimate @xmath191 of the velocity vector to the gradient field of the index , that is compute @xmath192 . an estimator of the velocity is a mapping from a couple of densities @xmath193 , @xmath194 to the tangent space at the initial density @xmath195 . the inverse of such a mapping is discussed under the name of _ retraction _ in . the simplest example here being @xmath196 , which is suggested by @xmath197 it is a common practice in engineering to use the initial velocity of the riemannian geodesic connecting @xmath189 to @xmath190 . this would require the computation of the geodesic itself , which is done using the computations schetched in app . [ sec : hessian ] . 29 [ 1]#1 [ 1]`#1 ` urlstyle [ 1]doi : # 1    r.  abraham , j.  e. marsden , and t.  ratiu . _ manifolds , tensor analysis , and applications _ , volume  75 of _ applied mathematical sciences_. springer - verlag , new york , second edition , 1988 . isbn 0 - 387 - 96790 - 7 . doi : 10.1007/978 - 1 - 4612 - 1029 - 0 . url http://dx.doi.org/10.1007/978-1-4612-1029-0 . absil , r.  mahony , and r.  sepulchre . _ optimization algorithms on matrix manifolds_. princeton university press , princeton , nj , 2008 . isbn 978 - 0 - 691 - 13298 - 3 . with a foreword by paul van dooren . shun - ichi amari . natural gradient works efficiently in learning . _ neural computation _ , 100 ( 2):0 251276 , feb 1998 . issn 0899 - 7667 . doi : 10.1162/089976698300017746 . url http://dx.doi.org/10.1162/089976698300017746 . shun - ichi amari and hiroshi nagaoka . _ methods of information geometry_. american mathematical society , providence , ri , 2000 . translated from the 1993 japanese original by daishi harada . vladimir  i. arnold . _ ordinary differential equations_. universitext . springer - verlag , berlin , 2006 . isbn 978 - 3 - 540 - 34563 - 3 ; 3 - 540 - 34563 - 9 . translated from the russian by roger cooke , second printing of the 1992 edition . nihat ay and ionas erb . on a notion of linear replicator equations . _ j. dynam . differential equations _ , 170 ( 2):0 427451 , 2005 . issn 1040 - 7294 . lawrence  d. brown . _ fundamentals of statistical exponential families with applications in statistical decision theory_. number  9 in ims lecture notes . monograph series . institute of mathematical statistics , 1986 . n.  n. encov . _ statistical decision rules and optimal inference _ , volume  53 of _ translations of mathematical monographs_. american mathematical society , providence , r.i . , 1982 . isbn 0 - 8218 - 4502 - 0 . translation from the russian edited by lev j. leifman . manfredo  perdigo do  carmo . riemannian geometry_. mathematics : theory & applications . birkhuser boston inc . , boston , ma , 1992 . isbn 0 - 8176 - 3490 - 8 . translated from the second portuguese edition by francis flaherty . joan esteban and debraj ray . on the measurement of polarization . _ econometrica _ , 620 ( 4):0 819851 , 1994 . issn 00129682 . doi : 10.2307/2951734 . url http://dx.doi.org/10.2307/2951734 . paolo gibilisco and giovanni pistone . connections on non - parametric statistical manifolds by orlicz space geometry . _ idaqp _ , 10 ( 2):0 325347 , 1998 . issn 0219 - 0257 . paolo gibilisco , eva riccomagno , maria  piera rogantin , and henry  p. wynn , editors . _ algebraic and geometric methods in statistics_. cambridge university press , cambridge , 2010 . isbn 978 - 0 - 521 - 89619 - 1 .    narendra  s. goel , samaresh  c. maitra , and elliott  w. montroll . on the volterra and other nonlinear models of interacting populations . modern phys . _ , 43:0 231276 , 1971 . issn 0034 - 6861 . josef hofbauer . on the occurrence of limit cycles in the volterra - lotka equation . _ nonlinear anal . _ , 50 ( 9):0 10031007 , 1981 . issn 0362 - 546x . doi : 10.1016/0362 - 546x(81)90059 - 6 . url http://dx.doi.org/10.1016/0362-546x(81)90059-6 . serge lang . _ differential and riemannian manifolds _ , volume 160 of _ graduate texts in mathematics_. springer - verlag , new york , third edition , 1995 . isbn 0 - 387 - 94338 - 2 . luigi malag . _ on the geometry of optimization based on the exponential family relaxation_. phd thesis , politecnico di milano , 2012 . luigi malag and giovanni pistone . a note on the border of an exponential family . arxiv:1012.0637v1 , 2010 . luigi malag and giovanni pistone . combinatorial optimization with information geometry : newton method . _ entropy _ , 16:0 42604289 , 2014 . luigi malag , matteo matteucci , and giovanni pistone . stochastic natural gradient descent by estimation of empirical covariances . in _ ieee congress on evolutionary computation _ , pages 949956 . ieee , 2011 . doi : 10.1109/cec.2011.5949720 . luigi malag , matteo matteucci , and giovanni pistone . natural gradient , fitness modelling and model selection : a unifying perspective . in _ ieee congress on evolutionary computation _ , pages 486493 . ieee , 2013 . y.  ollivier , l.  arnold , a.  auger , and n.  hansen . . arxiv:1106.3708 , 2011v1 ; 2013v2 . francisco  j. pino and jordi vidal - robert . habemus papam ? polarization and conflict in the papal states . http://www.afse-lagv.com/lagv/submissions/index.php/lagv2013/lagv12/paper/viewfile/362/152 , 2013 . giovanni pistone . algebraic varieties vs. differentiable manifolds in statistical models . in paolo gibilisco , eva riccomagno , maria piera rogantin , and henry  p. wynn , editors , _ algebraic and geometric methods in statistics _ , chapter  21 , pages 339363 . cambridge university press , 2009 . giovanni pistone . algebraic varieties vs differentiable manifolds in statistical models . in paolo gibilisco , eva riccomagno , maria  piera rogantin , and henry  p. wynn , editors , _ algebraic and geometric methods in statistics _ , pages 341365 . cambridge university press , cambridge , 2010 . giovanni pistone . nonparametric information geometry . in frank nielsen and frederic barbaresco , editors , _ geometric science of information _ , number 8085 in lncs , pages 536 , berlin heidelberg , 2013 . springer - verlag . first international conference , gsi 2013 paris , france , august 28 - 30 , 2013 proceedings . giovanni pistone and maria piera rogantin . the exponential statistical manifold : mean parameters , orthogonality and space transformations . _ bernoulli _ , 50 ( 4):0 721760 , august 1999 . issn 1350 - 7265 . johannes rauh , thomas kahle , and nihat ay . support sets in exponential families and oriented matroid theory . . j. approx . _ , 520 ( 5):0 613626 , 2011 . issn 0888 - 613x . doi : 10.1016/j.ijar.2011.01.013 . url http://dx.doi.org/10.1016/j.ijar.2011.01.013 . alessandro rinaldo , stephen  e. fienberg , and yi  zhou . on the geometry of discrete exponential families with application to exponential random graph models . _ electronic journal of statistics _ , 3:0 446484 , 2009 . stein et  al . _ sage mathematics software ( version 6.4.1)_. the sage development team , 2014 . from , with @xmath198 , @xmath199 , and @xmath200 , the lv equation @xmath201 has a stationary point @xmath202 with @xmath203 in the variables @xmath204 the equations are @xmath205    it follows that @xmath206 so that @xmath207 is constant . note that @xmath208 and @xmath209 this shows the existence of periodic orbits , see .      because of the periodic behaviour , we can not expect the dynamic project to the simplex @xmath210 . following @xcite , we can go up to @xmath211    we add a constant population of one individual @xmath212 and define the transformation @xmath213 @xmath214    note that @xmath215 . we have @xmath216    @xmath217    @xmath218    define @xmath219 then @xmath220 hence the differential equation for @xmath0 is the _ replicator equation _ the replicator equation eq . is nothing else then a particular class of differential equations in @xmath222 as a submanifold of of @xmath223 . let @xmath224 , where @xmath225 is viewed as a random variable on @xmath8 . then @xmath226}$ ] is a vector field of @xmath222 . the differential equation is @xmath227 , see . let @xmath228 be the simplex as a sub - variety of @xmath229 : @xmath230    we consider different parametrizations . @xmath231    the jacobian of @xmath232 is @xmath233      the exponential family @xmath234 gives @xmath235    as @xmath236 , @xmath199 , then @xmath237 ( the solid simplex parameters ) are the _ expectation parameters _ of the exponential family . the jacobian is @xmath238 notice that , if @xmath239 , then @xmath240      @xmath241    the jacobian is @xmath242      a differential equation @xmath243 in the simplex @xmath222 , considered as sub - variety of @xmath223 , computed componentwise , has the form : @xmath244 where @xmath245 has to be orthogonal to the constant vector , i.e. @xmath246 , or @xmath247 . because of this condition , it is enough to consider @xmath248    the differential equation , written in the parameter @xmath249 has the form @xmath250 , i.e. : @xmath251      in this case , the jacobian is and    we have : @xmath252 and eq . in matrix form become : @xmath253 which implies : @xmath254 and @xmath255 , i.e. eq . . the choice to leave out @xmath256 is of course arbitrary .      in this case the jacobian is , hence eq . becomes @xmath257    we need to compute the inverse of the lower block , i.e.    @xmath258    the differential equation in the parameters @xmath161 is    @xmath259    or , written in @xmath0 : @xmath260      eq . in matrix form becomes : @xmath261 whose solution in @xmath262 is : @xmath263 and written in @xmath0 : @xmath264    notice that , in each parametrization , the condition @xmath247 is verified . the differential equation of the replicator , see , is : @xmath265 on the simplex we have : @xmath266 where the expected value of @xmath138 is @xmath267 so that @xmath268 and eq . becomes : @xmath269 or @xmath270    replacing @xmath271 and @xmath272 in the differential equations with the different parametrizations , written as a function of @xmath273 , we have , in the exponential parametrization , @xmath274 so that the differential equations are @xmath275    in the projective parametrization @xmath276 hence the differential equations are    @xmath277 the expectation parameters are @xmath278 and the information matrix is @xmath279            let @xmath288 be an univariate statistical model . let @xmath289 and @xmath290 be centered pivotal quantities ( vector fields of the statistical model ) . the variation of @xmath291 in the time is : @xmath292 see @xcite .
Abstract:

Answer: the polarization measure is the probability that among 3 individuals chosen at random from a finite population exactly 2 come from the same class . this index is maximum at the midpoints of the edges of the probability simplex . we compute the gradient flow of this index that is the differential equation whose solutions are the curves of steepest ascent . tools from information geometry are extensively used . in a time series , a comparison of the estimated velocity of variation with the direction of the gradient field should be a better index than the simple variation of the index .


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: when a wave propagates with a low energy in a quasi - one - dimensional ( q1d ) wave guide with an uniform cross section , the guide is effectively one - dimensional ( 1d ) since only the first transverse mode plays a role . however , when a guide comprises discontinuities such as bends and branching junctions , higher transverse modes come into play at the vicinity of the discontinuities . therefore , in principle , a full - mode or full - wave analysis is required , which means a substantial increase in the amount of calculation and the physics is often obscured . knowing that the guide is essentially 1d in the uniform sections , such a full - wave analysis for the entire system is actually unnecessary . it is the purpose of this paper to present a scheme to eliminate those redundancies and simplify the calculations , yet have those low - energy properties faithfully preserved . historically , there has been a number of schemes related to this purpose . the earliest one was due to kuhn and dates back to 1949 @xcite . then it was made more well - known by griffith in 1953 @xcite , and henceforth it has been called the griffith boundary condition @xcite . the scheme contains a set of equations relating the wave functions and their first derivatives on the lines connected to a junction @xcite . the equations are simple equations that satisfy the unitarity condition at the junction . the scheme has been intuitively stated and has no undetermined parameters . however , recently it has been pointed out by the author @xmath0 @xcite that it is not clear what kind of realistic q1d guides the scheme describes . later than kuhn , there was another scheme by shapiro in 1983 @xcite . the scheme starts by an unitary matrix relating the amplitudes of the inward and outward waves @xcite . though this scheme contains free parameters , it also has not been mentioned how the parameters are related to the guides in realistic spaces .    to the present , approaches to the related problems are belong to either one of the mentioned two categories @xcite . and a common feature of most of the schemes @xcite is that the relation between the reduced systems and the original systems is not addressed , and this hinders the application to the study of realistic systems .    in 2006 , the author _ et al . _ proposed a scheme @xcite , which has been the first attempt to relate reduced and unreduced systems . the scheme resembles that by kuhn @xcite , but has an extra phenomenological term with a tuning parameter . the parameter is to be fixed by a comparison between results from the reduced and unreduced systems . the scheme was shown to be a substantial improvement , and many low - energy transport properties were shown to be captured . but still , it is pointed out later in this paper that this simple phenomenology can be inadequate and a more general scheme is needed . section [ form ] illustrates the derivation of the connecting equations to be used in reduced systems , for two typical component structures  the l - bend and the t - junction . then in sec . [ examples ] , these two structures are assembled into more elaborated structures , and the wave propagation in the structures are studied in and compared between , the reduced and unreduced systems . in sec . [ conc ] , a few concluding remarks are given . the scheme to be proposed can be summarized as in the following . for a particular discontinuity , the scattering matrix ( @xmath1-matrix ) is evaluated using a full - wave treatment , and then the matrix is truncated leaving only those elements relating the first transverse modes in the branch guides . the truncated @xmath1-matrix is then used to connect the 1d wave functions in reduced systems . the two - dimensional ( 2d ) time - independent schr@xmath2dinger equation ( tise ) @xmath3 [ \partial_x ^2 \psi + \partial_y ^2 \psi ] + v(x , y ) \psi = e \psi$ ] , with @xmath4 and hard wall boundaries is considered . for a wave with an energy @xmath5 in a guide having a width @xmath6 , the wave function can be written as a sum of direct products of transverse modes and longitudinal waves . labeling the guide by @xmath7 and defining a coordinate system @xmath8 in the guide , where @xmath9 is in the longitudinal direction and @xmath10 is in the transverse direction ( @xmath11 ) , the wave function can be written as @xmath12 where @xmath13 is the longitudinal wave number for the @xmath14-th transverse mode , which can be propagating or evanescent , and @xmath15 is a large enough integer . in a 1d space , the tise becomes @xmath16 \partial_x ^2 \psi + v(x ) \psi = e \psi$ ] . for a line labeled by @xmath7 , and with a coordinate @xmath9 defined on it , the wave function for @xmath17 is @xmath18 where @xmath19 is the longitudinal wave number given by @xmath20 . when results from 2d guides and 1d lines are compared , @xmath21 ( @xmath22 ) is compared with @xmath23 ( @xmath24 ) , and @xmath25 is compared with @xmath19 . two discontinuities in 2d wave guides are considered , the l - bend and t - junction as shown in fig . [ discontinuities](a ) and [ discontinuities](b ) respectively . the discontinuities are divided into regions which are labeled by 1 , 2 , 3 , and d in our discussion . coordinates are also defined in the branch guides , and wave functions in the branch guides are in the form given by eq . ( [ psi_2d ] ) . the @xmath1-matrices are to be presented in terms of these coordinates . the symbols used to denote the l - bend and t - junction in reduced systems are shown in fig . [ discontinuities](c ) and [ discontinuities](d ) respectively . we only sketch the evaluation of the @xmath1-matrices here , since the techniques are well - established and are detailed in the literatures @xcite . the wave function in region d [ see figs . [ discontinuities](a ) and [ discontinuities](b ) ] is also expanded in terms of undetermined amplitudes , and they are connected to wave functions in other regions by the conditions of continuities of wave functions and normal derivatives of wave functions at the boundaries between the regions . otherwise , one may also solve the tise in a discretized space , where the tise is a set of finite - difference ( fd ) equations .    for the t - junction shown in fig . [ discontinuities](b ) , we can get an equation for the amplitudes of the waves in the branch guides ( regions 1 , 2 , and 3 ) such as @xmath26 = \left [ \begin { array } { cccc } s_{11 } ^{\rm t } & s_{12 } ^{\rm t } & s_{13 } ^{\rm t } & ... \\ s_{21 } ^{\rm t } & s_{22 } ^{\rm t } & s_{23 } ^{\rm t } & ... \\ s_{31 } ^{\rm t } & s_{32 } ^{\rm t } & s_{33 } ^{\rm t } & ... \\ ... & ... & ... & ... \\ \end { array } \right ] \left [ \begin { array } { c } a^{(1)}_1 \\ a^{(1)}_2 \\ a^{(1)}_3 \\ ... \end { array } \right ] , \label { st_2d}\end{aligned}\ ] ] where @xmath27 ( @xmath28 ) is an amplitude for an inward ( outward ) wave in eq . ( [ psi_2d ] ) . truncating the matrices in eq . ( [ st_2d ] ) and retaining only terms related to the first transverse modes in the branch guides , we get @xmath29 = s _ { \rm t } \left [ \begin { array } { c } a_1 \\ a_2 \\ a_3 \end { array } \right ] , \label { teqs}\end{aligned}\ ] ] where @xmath23 and @xmath24 , @xmath30 and 3 are amplitudes in eq . ( [ psi_1d ] ) , and @xmath31 , \label { st_1d}\end{aligned}\ ] ] where the coordinates are chosen as toward the junction , and @xmath32 at the junction as shown in fig . [ discontinuities](d ) . for the scheme to work , it is necessary that the energy is such that only the first transverse mode is propagating , and the exponential tails of the evanescent waves of higher modes emanated from the discontinuities are shorter than the distances between the discontinuities . we can also write eq . ( [ teqs ] ) in terms of @xmath33 . using @xmath34 / 2 $ ] and @xmath35 / 2 $ ] at @xmath36 , the connecting equation can be rewritten as @xmath37 = \left ( 1 - s _ { \rm t } \right ) \left [ \begin { array } { c } \psi_1 \\ \psi_2 \\ \psi_3 \end { array } \right ] , \label { st_psi}\end{aligned}\ ] ] where the wave functions and their derivatives are evaluated at the discontinuity , and the directions of the coordinates are defined to be toward the discontinuity . note that in eq . ( [ st_psi ] ) , it is not necessary that the origins of the coordinates be located at the discontinuity .    likewise , the connecting equation for the l - bend shown in fig . [ discontinuities](a ) can be written as @xmath38 = \left ( 1 - s _ { \rm l } \right ) \left [ \begin { array } { c } \psi_1 \\ \psi_2 \end { array } \right ] , \label { sl_psi}\end{aligned}\ ] ] where the directions of the coordinates are defined to be toward the discontinuity as shown in fig . [ discontinuities](c ) . note that , the @xmath1-matrices are symmetric ( @xmath39 and @xmath40 ) and unitary ( @xmath41 and @xmath42 ) . the unitarity implies @xmath43 / [ 2 e^{ik x_\eta } ] $ ] and @xmath44 / [ 2 e^{-ik x_\eta } ] $ ] , the unitarity can be rephrased as a more intuitive equality @xmath45 , or there is no net inflow of probability current to the discontinuity . the numerical results for the magnitudes and arguments of the elements of @xmath46 are plotted in figs . [ smat](a ) and [ smat](b ) respectively , and the results for @xmath47 are plotted in figs . [ smat](c ) and [ smat](d ) respectively , versus a dimensionless longitudinal wave number @xmath48 defined by @xmath49 . for the reference of the readers , the cutoff of the second transverse mode is at @xmath50 . the results are obtained using discretized spaces with 20 sites across a width of @xmath6 , and they are found to be in congruence with results from 10 sites within windows of @xmath51 on the horizontal axes , and windows of @xmath52 and @xmath53 \simeq 0.005 $ ] on the vertical axes , which implies that the continuous space limit has been approached . since the l - bend and t - junction appear quite often in practical problems , it may be convenient to have their @xmath1-matrices in analytic forms . within a finite range of @xmath48 , it is possible to approximate a @xmath1-matrix by analytic functions . for the @xmath46 at @xmath54 , we may approximate the magnitudes and arguments of the elements by @xmath55 within @xmath54 , the approximation for @xmath56 in eq . ( [ sl11mag_approx ] ) and the approximation for @xmath57 in eq . ( [ sl11arg_approx ] ) approximate the numerical results in figs . [ smat](a ) and [ smat](b ) respectively up to @xmath58 and @xmath59 < 0.01 $ ] accuracy . the exact equalities eqs . ( [ sl12mag])-([sl_approx ] ) are due to the symmetry and unitarity of the @xmath1-matrix and the exchange of leads 1 and 2 .    for the @xmath47 at @xmath54 , it is found that the numerical results for the elements can be approximated by @xmath60 for @xmath54 , the above approximations for @xmath61 and @xmath62 approximate the results in figs . [ smat](c ) and [ smat](d ) respectively up to @xmath63 and @xmath64 < 0.01 $ ] accuracy . the exact equalities are due to the symmetry and unitarity of the @xmath1-matrix , and the exchange of leads 2 and 3 . in this section , three wave guides which are composites of the discussed l - bend and t - junction are analyzed . the scattering amplitudes from the original 2d structures and the reduced multiply - connected 1d structures are compared . cases of far apart and close discontinuities , different orientations of guides , and straight and smoothly curved guides are considered . the first example is a 2d square loop resonator with two leads as depicted in fig . [ loop1](a ) . the translational invariant sections have the same width @xmath6 . for simplicity , the distances between the discontinuities are chosen to be the same , and are denoted by @xmath65 as shown . in fig . [ loop1](b ) , a reduced version for the structure in fig . [ loop1](a ) is shown . the magnitude of the transmission scattering amplitude @xmath66 has been plotted versus @xmath48 for the unreduced and reduced systems , in figs . [ loop1](c ) and [ loop1](d ) , for @xmath67 and @xmath68 respectively . the @xmath69 for the 2d structure is the scattering amplitude from the first transverse mode in lead 2 to the first transverse mode in lead 1 , and the dimensionless wave number @xmath48 is now defined by @xmath49 or @xmath70 depending on the context . the procedure for a full - wave evaluation of @xmath66 for the 2d structure is standard @xcite and it is not to be repeated here , but only the result is given . the result here is obtained with a fd tise , and the number of sites across a width @xmath6 is equal to 20 .    in the reduced system shown in fig . [ loop1](b ) , a line labeled by @xmath7 and given a coordinate @xmath9 has a wave function in the form given in eq . ( [ psi_1d ] ) . for the coordinates defined in fig . [ loop1](b ) , the wave functions on the lines are connected by @xmath71 = \left ( 1 - s _ { \rm l } \right ) \left [ \begin { array } { c } \left . \psi_3 \right| _ { x_3=d } \\ \left . \psi_4 \right| _ { x_4=0 } \end { array } \right ] ,    \label { node1 } \\ \displaystyle { { 1 + s _ { \rm l } } \over { ik } } \left [ \begin { array } { c } \left . ~ \displaystyle { { d\psi_4 } \over { dx_4 } } \right| _ { x_4=d }   \\ \left . - \displaystyle { { d\psi_5 } \over { dx_5 } } \right| _ { x_5=0 }   \end { array } \right ] = \left ( 1 - s _ { \rm l } \right ) \left [ \begin { array } { c } \left . \psi_4 \right| _ { x_4=d } \\ \left . \psi_5 \right| _ { x_5=0 } \end { array } \right ] , \label { node2 } \\ \displaystyle { { 1 + s _ { \rm t } } \over { ik } } \left [ \begin { array } { c } \left . - \displaystyle { { d\psi_3 } \over { dx_3 } } \right| _ { x_3=0 }   \\ \left . ~ \displaystyle { { d\psi_1 } \over { dx_1 } } \right| _ { x_1=0 }   \\ \left . ~ \displaystyle { { d\psi_6 } \over { dx_6 } } \right| _ { x_6=d }   \end { array } \right ] = \left ( 1 - s _ { \rm t } \right ) \left [ \begin { array } { c } \left . \psi_3 \right| _ { x_3=0 } \\ \left . \psi_1 \right| _ { x_1=0 } \\ \psi_6 \right| _ { x_6=d } \end { array } \right ] ,   \label { node3 } \end{aligned}\ ] ] and @xmath72 = \left ( 1 - s _ { \rm t } \right ) \left [ \begin { array } { c } \left . \psi_5 \right| _ { x_5=d } \\ \left . \psi_2 \right| _ { x_2=0 } \\ \left . \psi_6 \right| _ { x_6=0 } \end { array } \right ] . \label { node4a}\end{aligned}\ ] ] this contains 10 equations with 10 unknowns , when @xmath73 and @xmath74 are given . the scattering amplitude @xmath69 is obtained as @xmath75 at @xmath76 and @xmath77 . the magnitude @xmath66 is plotted in figs .  [ loop1](c ) and [ loop1](d ) , for values of @xmath65 corresponding to the original 2d structure . in the calculation , @xmath46 and @xmath47 use the numerical values shown in fig . [ smat ] .    in the case of @xmath67 [ see fig . [ loop1](c ) ] , results from the original 2d and the reduced systems are nearly indistinguishable when seen in the size of the plot , while in the case of @xmath68 [ see fig .  [ loop1](d ) ] , the two results have a perceptible difference , especially when @xmath48 becomes large . a criterion for the applicability of the one - mode reduction scheme is @xmath78 , which means that the exponential tails of the evanescent waves for the higher transverse modes emanating from the discontinuities , is much shorter than the distances between the discontinuities . in this regime , the evanescent waves from neighboring discontinuities do not overlap each other , and the discontinuities communicate with each other only via the first transverse modes in the translational invariant sections . hence , the higher transverse modes in the sections are redundant , and a single - mode description is adequate . the length of an exponential tail of the second transverse mode is of the order of @xmath79 . the values of @xmath80 are approximately equal to 0.18 ( @xmath81 ) , 0.23 ( @xmath82 ) , 0.28 ( @xmath83 ) , 0.37 ( @xmath84 ) , 0.48 ( @xmath85 ) , 0.96 ( @xmath86 ) , and 1.56 ( @xmath87 ) , for the values of @xmath48 given in the brackets @xcite . it is seen that @xmath88 becomes large only when @xmath48 approaches @xmath89 . for @xmath67 , @xmath90 at @xmath83 , and that justifies the reduction scheme in the entire range of @xmath48 in fig . [ loop1](c ) [ and also fig . [ loop2](c ) later in this section ] . for @xmath68 , @xmath91 at @xmath92 , and @xmath93 is certainly not `` large '' for fig . [ loop1](d ) [ and also fig . [ loop2](d ) later in this section ] . in spite of this , the scheme might still perform up to certain precision , until it really starts to breakdown at @xmath94 as seen in fig . [ loop1](d ) [ see also fig .   [ loop2](d ) ] . however , its reliability in this regime of @xmath48 is uncontrolled in general . another 2d wave guide as shown in fig . [ loop2](a ) is also analyzed . this structure resembles the one in fig . [ loop1](a ) , except that one of the leads is rotated by 90-degree . the width of the uniform sections are also denoted by @xmath6 , and the distances between the discontinuities are also denoted by @xmath65 . a reduced structure for the guide is shown in fig . [ loop2](b ) , and the magnitude of the transmission scattering amplitude @xmath66 is also plotted for both of the reduced and unreduced structures in figs . [ loop2](c ) [ for @xmath67 ] and [ loop2](d ) [ for @xmath68 ] . the connecting equations for the reduced structure in fig . [ loop2](b ) are eqs . ( [ node1 ] ) , ( [ node2 ] ) , ( [ node3 ] ) , and @xmath97 = \left ( 1 - s _ { \rm t } \right ) \left [ \begin { array } { c } \left . \psi_6 \right| _ { x_6=0 } \\   \left . \psi_2 \right| _ { x_2=0 } \\ \left . \psi_5 \right| _ { x_5=d } \end { array } \right ] . \label { node4b}\end{aligned}\ ] ] the transmission scattering amplitude @xmath69 is found as in the previous example . likewise , it is seen in fig . [ loop2](c ) that the reduction scheme is guaranteed to work in the @xmath78 regime . also , it is seen in fig . [ loop2](d ) that the scheme might still work qualitatively or semi - quantitatively , when @xmath48 is departed from this regime .    comparing the result in fig . [ loop1](c ) with that in fig . [ loop2](c ) , and the result in fig . [ loop1](d ) with that in fig . [ loop2](d ) , it is seen that scattering amplitudes can depend significantly on the orientations of the branch guides at a discontinuity . this indicates that reduction schemes with symmetric branch lines such as those in refs . @xcite are not adequate for some cases . the third example is a 2d annulus structure with an inner and an outer radii of @xmath98 and @xmath99 respectively , and two mutually perpendicular leads of width @xmath6 radially connected to the annulus as shown in fig . [ ring](a ) . figure [ ring](b ) shows a reduced version of it . for the 2d annulus , we may follow a mode - matching full - wave treatment formulated by xia and li @xcite . in the annulus the wave function @xmath100 can be expanded by radial and angular modes , @xmath101 , where a radial mode is given by @xmath102 , and @xmath103 . the @xmath104 and @xmath105 are the radial and angular coordinates respectively ; and the @xmath106 and @xmath107 are the bessel functions of the first and second kinds respectively . at the inner radius , @xmath108 for any @xmath105 ; at the outer radius , @xmath109 when @xmath105 is away from the leads , and @xmath110 , when @xmath105 is in the range of lead @xmath7 . in addition , the radial derivative @xmath111 is equated with the longitudinal derivative @xmath112 when @xmath100 and @xmath113 meet at the outer arc of the annulus . the difference between the straight transverse cuts of the leads and the outer arcs of the annulus is neglected . the wave functions in the leads and the annulus are hence matched , and one can get a set of equations relating the coefficients of the modes in the different regions . the transmission scattering amplitude @xmath66 for the 2d annulus is plotted for @xmath114 [ fig . [ ring](c ) ] and @xmath115 [ fig .  [ ring](d ) ] .    to apply a reduced calculation to the reduced system in fig . [ ring](b ) , note that for smoothly curved guides with small curvatures , the back - scattering is small and the guides can be treated as reflectionless for most purposes . guides with constant curvatures are translational invariant and indeed reflectionless , though the lengths and widths may not be rigorously defined . the junctions are treated as t - junctions . for the left junction , @xmath116 = \left ( 1 - s _ { \rm t } \right ) \left [ \begin { array } { c } \left . \psi_1 \right| _ { x_1=0 } \\ \left . \psi_3 \right| _ { x_3=0 } \\ \left . \psi_4 \right| _ { x_4=d_4 } \end { array } \right ] ; % \label { } \end{aligned}\ ] ] for the right junction , @xmath117 = \left ( 1 - s _ { \rm t } \right ) \left [ \begin { array } { c } \left . \psi_2 \right| _ { x_2=0 } \\ \left . \psi_3 \right| _ { x_3=d_3 } \\ \left . \psi_4 \right| _ { x_4=0 } \end { array } \right ] . % \label { } \end{aligned}\ ] ] we let @xmath118 , where @xmath119 and @xmath120 , and the subtraction is to approximately exclude the regions in the junctions [ see region d in fig .  [ discontinuities](b ) ] . the radius @xmath114 gives @xmath121 and @xmath122 , and @xmath115 gives @xmath123 and @xmath124 . results for @xmath66 are plotted in figs . [ ring](c ) and [ ring](d ) . if one analyzes the reduced system using the griffith scheme @xcite for a junction with three branches , the left junction has @xmath125 and the right junction has @xmath126 note that the equations at a junction are symmetric with respect to an interchange of any two branches . using the same @xmath127 and @xmath128 , results for @xmath66 are also plotted in figs . [ ring](c ) and [ ring](d ) . comparing the results from our scheme and the griffith scheme with the result from the 2d calculation , it is seen that the griffith result is qualitatively different from the 2d result in general , while our scheme captures the essential features in the 2d result , especially in the case of longer guide sections [ fig .  [ ring](c ) ] . the lengths of the guide sections in fig . [ ring ] are about the sizes of those in figs . [ loop1 ] and [ loop2 ] , but the disagreement between the 2d and reduced calculation results are seen to be more severe in fig . [ ring ] , especially in the case of short guide sections [ fig . [ ring](d ) ] . this is mainly due to vaguer notions of the lengths and widths of the guide sections between the discontinuities , and a stronger distortion of the shapes of the junctions from a `` t '' [ fig . [ discontinuities](b ) ] , when the guide sections are curved and short , in addition to a violation of the criterion @xmath129 . comparing figs . [ ring](c ) and [ ring](d ) , it is also seen that our scheme performs better when the discontinuities are more apart . this study has shown that the higher transverse modes in the uniform - cross - sectional sections in wave guides can give only a minor effect on the low - energy properties of the guides . a scheme to remove these modes and a criterion for the energy range ( @xmath78 ) in which the scheme applies have been proposed . in the scheme , a reduced system and its corresponding system in a realistic space have a sound relationship , and therefore the scheme may find more practical use in the analyses of realistic wave guides than previously reported schemes @xcite . the electronic spin degrees of freedom may also be included into the scheme by expanding the current one - mode @xmath1-matrices to two - mode @xmath1-matrices , for the up and down spin channels , and in the same way as in sec . [ examples ] , the two - mode @xmath1-matrices are calibrated by full - wave calculations for systems in realistic spaces . likewise , the precision of the scheme at short guide sections can also be refined by including more transverse modes and using multi - mode @xmath1-matrices . k. ruedenberg and c. w. scherr , j. chem . phys . * 21 * , 1565 ( 1953 ) . p. exner and p. seba , rep . * 28 * , 7 ( 1989 ) . t. kottos and u. smilansky , ann . of phys . * 274 * , 76 ( 1999 ) . xia , phys . rev . b * 45 * , 3593 ( 1992 ) . j. m. mao , y. huang , and j. m. zhou , j. appl . * 73 * , 1853 ( 1993 ) . p. singha deo and a. m. jayannavar , phys . b * 50*,11629 ( 1994 ) . m. v. moskalets , low temp . * 23 * , 824 ( 1997 ) . ryu and s. y. cho , phys . b * 58 * , 3572 ( 1998 ) . s. bandopadhyay , p. singha deo , and a. m. jayannavar , phys b * 70 * , 75315 ( 2004 ) . d. bercioux , m. governale , v. cataudella , and v. m. ramaglia , phys . * 93 * , 56802 ( 2004 ) . p. foldi , b. molnar , m. g. benedict , and f. m. peeters , phys . b * 71 * , 33309 ( 2005 ) . u. aeberhand , k. wakabayashi , and m. sigrist , phys . b * 72 * , 75328 ( 2005 ) . x. f. wang and p. vasilopoulos , phys . b * 72 * , 165336 ( 2005 ) . o. kalman , p. foldi , m. g. benedict , and f. m. peeters , arxiv:0806.2734 ( unpublished ) .
Abstract:

Answer: a scheme to reduce translational noninvariant quasi - one - dimensional wave guides into singly or multiply connected one - dimensional ( 1d ) lines is proposed . it is meant to simplify the analysis of wave guides , with the low - energy properties of the guides preserved . guides comprising uniform - cross - sectional sections and discontinuities such as bends and branching junctions are considered . the uniform sections are treated as 1d lines , and the discontinuities are described by equations sets connecting the wave functions on the lines . the procedures to derive the equations and to solve reduced systems are illustrated by examples , and the scheme is found to apply when the discontinuities are distant and the energy is low . when the scheme applies , it may substantially simplify the analysis of a wave guide , and hence the scheme may find uses in the study of related problems , such as quantum wire networks . wave guide , discontinuity , quantum wire , one - dimensional ,    43.20.mv,73.23.ad,73.63.nm,84.40.az


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: fk comae berenices ( ) , with its extreme activity , is a key target to study in the solar - stellar connection context . it is a rapidly rotating and apparently single g5 ii giant , even though the spectral classification can go as late as g8  iii , due to the fast rotation rate and spectral oddities . this star has a rotational velocity of _ v  @xmath4  i _ =  162.5  @xmath5  3.5  km  s@xmath6 and radial velocity -24  @xmath5  3  km  s@xmath6 ( huenemoerder et al . @xcite hereafter refered to as hrbn ) . the photometric and rotational period is 2.400 day ( jetsu et al . this gives a projected stellar radius of 7.7  @xmath5  0.2  r@xmath7 ( hrbn ) . as a member of the old disk aggregate , its distance and mass are approximately known : d  @xmath8  200300  pc and m@xmath9  1.5  m@xmath7 ( eggen & iben @xcite ) . some spectral peculiarities of fk com were first observed by merrill ( @xcite ) . amongst the oddities present in the spectra of this star can be included a strong , broad and erratically variable h@xmath0  emission ( e.g. ramsey et al . @xcite ) , strong uv chromospheric emission lines reminiscent of those observed in rs cvn binaries ( bopp & stencel @xcite ) , a quasi - sinusoidal light curve ( holtzman & nations @xcite ) and high x - ray luminosity ( walter @xcite ) . some of these spectral characteristics led bopp & stencel ( @xcite ) to define a new class of variable stars with fk comae as the prototype . the fk comae - type variables were thus defined as rapidly rotating , single , g - k type giants with strong chromospheric and transition region emission . stars belonging to this class are extremely rare . the space density of these stars was estimated by collier ( @xcite ) as 2  @xmath1010@xmath11  pc@xmath12 , indicating , when compared with the space density of the rare hertzsprung - gap g giants , that the evolution through this part of the h - r diagram must be extremely rapid . typical values of rotation rate from g - type stars are @xmath13  10  km  s@xmath6 for single stars and 3050  km  s@xmath6 for stars in close binary systems . the extreme rotational velocity , near breakup , of fk comae raises interesting questions on the possible evolutionary status of this star . two scenarios were advanced to explain the observed rotational and activity characteristics : accretion from a very small unseen companion or product of the recent coalescence of a close binary system . based on radial velocity variations measurements , mccarthy & ramsey ( @xcite ) proved the first hypothesis very unlikely , by imposing severe limits to the mass of this unseen companion . the currently accepted theoretical scenario is that fk comae - type stars are recently coalesced binary systems . according to webbink ( @xcite ) , w uma systems , due to the thermal instabilities caused by the primary s departure from main - sequence , evolve into configurations that can only be stabilized by a big loss of angular momentum , coalescing into a single star in a way that is mostly independent of their evolution . the rapid rotation would arise from the angular momentum of the former binary and the activity level from a rotation - driven magnetic dynamo ( ramsey et al . @xcite ) . the descendents of these binaries would be indistinguishable from single stars . welty & ramsey ( @xcite ) presented evidence that conforms with such scenario , based on the analysis of the width and shape of the spectral lines .    in fk comae , the most striking spectral feature is the h@xmath0  emission line . this line profile is extremely broad , asymmetric and presents variations on different time scales . the same characteristics are present in higher balmer lines , even though appearing less extreme due to the different optical thickness of the lines . the phase modulation in these lines indicates that the emitting extended material is corotating with the star . based on the h@xmath0 to h@xmath1 flux ratios , hrbn have proposed that the emission arises in structures similar to solar prominences . we present here a model of the h@xmath0  and h@xmath1  line profiles for a quiet phase of the circumstellar environment , derived from a data set obtained in may 96 . the model includes a 3d gaussian spatial distribution of emitting matter , extending to several stellar radii , and an optically thick near - stellar shell of absorbing material . the data sets from may and june 97 represent a more active state , demanding a different modelling technique . the h@xmath1  line profiles present a complex structure that we model as the superposition of gaussian - shaped components . the same components were found to be present in the h@xmath0  line profiles , superposed on a more distributed `` bulk '' emission . we were able to track these components in the velocity vs phase diagram for approximately 3 rotational periods . a large flare event was detected in our may 97 data set , in the balmer and d3 lines , and we computed the corresponding energy budget . the multi - site continuous spectroscopy ( musicos ) project was developed with the aim of providing continuous spectroscopic coverage ( in a large wavelength range and with instruments as similar as possible ) for a wide range of stellar objects ( e.g. foing & catala @xcite ; foing et al . @xcite ; catala et al . @xcite , @xcite ; huang et al . @xcite ) . following these orientations , the first fiber - fed echelle cross - dispersed musicos spectrograph was designed and built in meudon - paris observatory ( baudrand & bhm @xcite ) and a replica was built at the esa space science department in estec and adapted for the int . the esa - musicos spectrograph was installed at the 2.5  m isaac newton telescope ( int ) , at the _ observatorio del roque de los muchachos _ , la palma , spain . the commissioning campaign of this instrument took place from april 29th to may 6th . our first data set was obtained in this campaign . the resolving power achieved was @xmath14  35000 . the second data set was obtained in may  97 also at the int with the same instrument . the ccd detector used in both runs was a tektronix , 1024  @xmath10  1024 pixels , with a readout noise of @xmath8  5  e@xmath15 . to reduce the int data , we used the midas nov . 94 and nov . 96 echelle reduction packages . the reduction procedure used was standard for echelle spectra , with some problems in the correction for the blaze function due to the low level of counts of the flat - field exposures in the bluer part of the spectra . in addition to these int observations , we obtained in june  97 more h@xmath0  spectra at the _ observatoire de haute provence _ ( ohp ) , in france . we used the aurlie spectrograph at the 1.52  m coud telescope , with resolving power of 30000 . the reduction of this data set was done using the on - line reduction facilities available on site .    in table 1 we give the summary information on all the spectra obtained . the ephemeris used for the phase computations is from chugainov ( @xcite ) : 2442192.345  +  2.400e .    .list of 1996 and 1997 fk comae observations [ cols="^,<,^,>,<,>",options="header " , ]     in table  2 we give the parameters for each of the gaussian components fitted to the balmer line profiles , h@xmath0  and h@xmath1 . in fig . [ vel_all ] , we show several plots that describe their behaviour : from top to bottom , the peak intensity vs velocity , and velocity vs phase for both lines . in order to highlight their phase behaviour , we also plotted ( in the two lower graphs ) relevant velocity curves . figure  [ vel_time ] shows the velocity evolution of the components in time . we have estimated a maximum error in the velocity determinations of about 50  km  s@xmath6 . the plots for these two balmer lines were analysed together in order to clarify the structures position and evolution .      as it can be seen in fig . [ vel_all ] , three structures could be followed in phase , corresponding to v@xmath16  =  245  @xmath17 , v@xmath18  =  285  @xmath19 and v@xmath20  =  487  @xmath19 , both in h@xmath0  and h@xmath1 . it was observed that most of the components were detected within approximately @xmath5  2  _ v  @xmath4  i _ from the rotational axis , in agreement with welty et al . ( @xcite ) that placed the emitting material at @xmath21  2  _ r@xmath22_.    similar component analysis was carried out with the h@xmath0  spectra obtained at the ohp in june  97 . two representative spectra were fitted and the velocity positions of the components were compared with what was obtained for may  97 . even though some profiles appear similar ( fig . [ all97 ] , for phase 1.971 and 16.974 ) , the previously fitted velocity curves do not seem to describe the phase behaviour now observed ( fig . [ vel_all ] middle graph ) . in particular , one can notice the appearance of a new component near - zero velocity at phase @xmath8  0.47 , compared with one month earlier . however , the poor phase coverage in june  97 makes a detailed comparison ambiguous . we also find an increase of h@xmath0 equivalent width of 2    in june  97 compared to the quiet level of may  97 , for the two phases shown . most of the components detected in h@xmath0  have their counterparts in h@xmath1 . still , some higher velocity h@xmath0  components are not seen in h@xmath1 . this is due to the better detection sensitivity in h@xmath0 as the signal - to - noise ratio in h@xmath1  is lower , but also to the presence of the `` bulk '' emission in h@xmath0  that extends to higher velocities ( as it was also seen in the quiet state in may  96 we modelled earlier ) . a slight mismatch can be observed in the velocity position of the components in both lines . intrinsic measurement uncertainties and physical phenomena combine to create this effect . the origin of the emission is obviously not in point sources , their vertical extent making difficult to determine the center of mass of the matter distribution , due to projection effects and emission gradients in the extended structures . this lack of exact agreement between the two lines is then justified given the simplified geometrical model adopted and the ignoring of radiative transfer effects . the d3 ( @xmath2 5876   ) line is very weak or absent in inactive stars spectra . the presence of this line in cool stars like fk com indicates the existence of non - radiative heating processes , since the line can not be excited at photospheric temperatures . in fk comae , this line has two components : an excess absorption and near - surface emission wings . a non - uniform shell of excited could explain this line profile , with the corona as source of ionizing radiation as it is observed in the sun and active stars . a similar multi - component study was made for the d3 line . all the spectra present an absorption at near - zero projected velocity , even though its strength is variable . emission wings were detected above the limb , confirming the existence of a slightly extended shell of excited . the variability observed in the spectra can represent a back - radiation of a nonuniform coronal emission onto the upper chromosphere . the modelling technique we used did not allow us to take advantage of the rotational phase information that might be present in the data set , due to the complexity of the profiles . a large flare event was detected in these two balmer lines and also in d3 . the flaring material seems to extend from the stellar surface up to at least 3 stellar radii . the two spectra where we detected an increase in the excess emission , what we call a pre - flare spectrum and a flare spectrum , ( respectively ut date 16.0 and 16.8 in fig . [ all97 ] ) are marked by `` _ f _ '' in figs . [ vel_all ] and [ vel_time ] . if these spectra represent the same flaring structure viewed at different phases , we can predict a velocity curve allowing to identify the post - flare structures at later phases . such a velocity curve is also plotted in fig . [ vel_time ] , corresponding to _ v@xmath23 _  @xmath8  245 kms@xmath6 ( r  =  1.50 _ r@xmath24  i _ ) and @xmath25  =  0.612 . at five other phase positions the same structure is identified . indeed , this structure was first detected on may 15.8 and it was still present on may 21.8 , when our run ended . when transiting in front of the stellar disk , near zero projected velocity , this structure ( extended vertically ) would have smaller projected dimension , and in the spectrum we should detect a narrow peaked component . this is in fact observed in the spectrum corresponding to ut date 19.8 in fig . [ components ] . this flare event detected in the balmer lines caused an increase in equivalent width of 7  @xmath5  1    in h@xmath0  and 1.2  @xmath50.2    in h@xmath1  ( fig . [ ew ] ) . in the case of the measurements were very difficult , the flare increase in equivalent width being the only variation detected above the errorbars . using the cousins - bessel magnitudes for fk com and assuming a distance of 215  pc ( hrbn ) some flare energetics were computed . the calculated energy released at flare maximum is 8  @xmath5  1  @xmath1010@xmath26  erg  s@xmath6 for h@xmath0  and 1.1  @xmath5  0.1  @xmath1010@xmath26  erg  s@xmath6 for h@xmath1 . the h@xmath0  to h@xmath1  flux ratio is about 7 , perhaps indicating much lower densities when compared with solar flares . again the interpretation of such global ratios is dubious in this type of analysis . the total energy released during the flare for h@xmath0  is 11  @xmath5  3  @xmath1010@xmath27  erg . the uncertainties in the intrinsic flare radiative losses for h@xmath0  come mainly from the choice of reference quiet spectrum and residual emission from active regions around the flare area . comparing this energy with the one reported on the giant h@xmath0  flare , @xmath8  4  @xmath1010@xmath28  erg ( foing et al . @xcite ) , this seems to be possibly the largest h@xmath0  flare reported on a cool star . what is also remarkable is the long duration of this flare event ( total equivalent flare duration approximately 2.5  @xmath5  0.5  @xmath1010@xmath29  s ) . the ratios of the continuum ubvri counterparts to h@xmath0 , for similar flare events , are reported to be from 120 ( avrett et al . @xcite ; machado et al . @xcite ; houdebine @xcite ) up to 340 ( foing et al . the continuum ubvri counterpart of the fk comae flare could therefore amount to losses of 1.3  @xmath1010@xmath30 erg ( for a flare ratio ubvri / h@xmath0  of 120 ) . this can be compared to total radiative budgets for the largest reported cool star flares , 1.2  @xmath1010@xmath31  erg on hr  1099 ( foing et al . @xcite ) or 1.8  @xmath1010@xmath30  erg on ( cutispoto et al . thus , we can speculate that we may have observed an exceptionally energetic flare , but we have no independent photometric measurements to quantify the continuum counterpart . an enhancement of emission was measured at pre - flare and flare maximum with equivalent width 0.08 and 0.19   respectively and at a velocity consistent with h@xmath0  and h@xmath1  components . it was not measurable during the gradual phase , and in particular not in the post - flare spectra of may  19.8 . this is similar with measurements in solar and red dwarf flares , which show a more impulsive behaviour of this line , while the balmer lines show enhancements both in impulsive and the gradual phase ( foing @xcite ) . it seems also that the d3 flare maximum was reached before the h@xmath0  flare maximum . the changes observed in the large flare component width might also be associated with large kinetics ( turbulence and flows ) related with the flare event ( foing @xcite , foing et al . once again , the analysis of the line profiles gives us abundant evidence on the complexity of the circumstellar environment of the peculiar giant star fk comae .      for the more complex structure of the 1997 profiles a different modelling technique is applied . we have successfully identified and described the phase behaviour of emission components . the h@xmath1  spectral line was described as the superposition of several ( gaussian ) emission components . the velocity vs phase diagram allows to trace the same emitting structures as they corotate with the star . in the case of h@xmath0 , the profile is interpreted as a quasi - stationary distributed `` bulk '' emission , to which discrete components , lasting several days , are superposed . the same components were detected in both balmer line profiles . a strong flare event was detected in these two balmer lines and also in d3 . the flaring material seems to extend from the stellar surface up to at least 3 stellar radii . the evolution of the flaring structure was described by a velocity curve corresponding to _ v@xmath23 _ @xmath8  245 kms@xmath6 or r  =  1.50  _ r@xmath24  i _ from the rotational axis . the total energy released during the flare in h@xmath0  is @xmath8  10@xmath3 erg , the largest h@xmath0  flare energy reported on a cool star . the data set of june  97 , as it presents longer time intervals continuously covered , provides some more information . short time - scale variations were observed , that can not be associated with rotational modulation . several variable features are visible ; they can be interpreted as emission bursts with a time scale of about 1 hour that can be related to the excitation and recombination of giant magnetic loop structures around the star ( oliveira et al . @xcite ) . the d3 line is always present in the spectra , with absorption and emission wings components , as a clear indicator of the activity in this star . an absorption near - zero projected velocity is always present with variable intensity . the profiles also present near - stellar emission wings , particularly enhanced during the described flare event . the need is clear for a multi - site continuous data set to disentangle the intrinsic variations from rotational modulation . we have organized such a campaign in february / march 1998 , involving the following telescopes / observatories : int ( with the esa - musicos spectrograph ) , telescope bernard lyot , france ( with the musicos spectropolarimeter ) , david dunlop observatory ( canada ) , russian special astrophysical observatory ( sao ) and mcdonald observatory . the phase coverage thus obtained represents a major improvement in the analysis of the circumstellar emission on fk comae . this will be reported in a subsequent paper ( oliveira et al . in preparation ) . oliveira s research work is being supported by the praxis xxi grant bd9577/96 from the _ fundao para a cincia e a tecnologia _ , portugal . the authors would like to thank the ing / int and ohp personnel for their support . we thank p.  sonnentrucker , p.  ehrenfreund , j.a . de  jong , h.c . stempels , r.  le  poole , p.  gondoin , t.  beaufort , j.  telting , n.  walton and j.th . van  loon for support observations and for discussion . we acknowledge the referee l.w . ramsey for useful comments . jacco , you know what ? you know ...    avrett  e.h . , machado  m.e . , kurucz  r.l . , 1984 , in proceedings of `` lower atmosphere in solar flares '' , neidig  d.f . ( ed . ) , sunspot baudrand  j. , bhm  t. , 1992 , a&a  259 , 711 bessel  m.s . , 1979 , pasp  91 , 589 bopp  b.w . , stencel  r.e . , 1981 , aj  247 , 131 catala  c. , foing  b.h . ( eds . ) , 1990 , proceedings of the 2nd musicos workshop , meudon obs . catala  c. , foing  b.h . , baudrand  j. et al . , 1993 , a&a  272 , 176 chugainov  p.f . , 1976 , krymsk . astrofiz 54 , 89 collier  a.c . , 1982 , mnras  200 , 489 cutispoto  g. , pagano  i. , rodono  m. , 1992 , a&a  263 , l3 eggen  o.j . , iben  i. jr , 1989 , aj  97 , 431 foing  b.h . , 1989 , soph  121 , 117 foing  b.h . , catala  c. , 1991 , adspr  11 , 159 foing  b.h . , char  s. , ayres  t. et al . , 1994a , a&a  92 , 543 foing  b.h . , catala  c. , baudrand  j. et al . , 1994b , in asp conf . 64 , the 8th cambridge workshop in cool stars , stellar systems and the sun , j .- p.caillaut ( ed . ) , 699 holtzman  j.a . , nations  h.l . , 1984 , aj  89 , 391 houdebine  e.r . , 1992 , j.  20 , 213 huang  l. , zhai  d .- s . , catala  c. , foing  b.h . ( eds . ) , 1995 , proceedings of the 4th musicos workshop , publ . beijing/  toulouse/  estec huenemoerder  d.p . , ramsey  l.w . , buzasi  d.l . , nations  h.l . , 1993 , apj  404 , 316 ( hrbn ) jetsu  l. , pelt  j. , tuominen  i. , 1993 , a&a  278 , 449 machado  m.e . , avrett  e.h . , vernazza  j.e . et al . , 1980 , apj  242 , 336 mccarthy  j.k . , ramsey  l.w . , 1984 , apj  283 , 200 merrill  p.w . , 1948 , pasp  60 , 382 oliveira  j.m . , foing  b.h . , gondoin  p. et al . , 1998a , in `` new perspectives on solar prominences '' , iau colloquium 167 proceedings , webb  d. , schmieder  b. , rust  d. ( eds . ) , 243 oliveira  j.m . , foing  b.h . , gondoin  p. et al . , 1998b , in asp conf . 154 , the 10th cambridge workshop in cool stars , stellar systems and the sun , donahue  r.a . , bookbinder  j.a . ( san francisco : asp ) , cd-1524 ramsey  l.w . , nations  h.l . , barden  s.c . , 1981 , apj  251 , 101 walter  f.m . , 1981 , apj  245 , 677 webbink  r.f , 1976 , apj  209 , 289 welty  a.d . , ramsey  l.w . , iyengar  m. , nations  h.l . , buzasi  d.l . , 1993 , pasp  105 , 1427 welty  a.d . , ramsey  l.w . , 1994 , apj  435 , 848
Abstract:

Answer: we present results on spectroscopic observations of the fast - rotating active giant , obtained mainly with the esa - musicos spectrograph at the isaac newton telescope ( int ) in 1996 and 1997 and also with the aurlie spectrograph at the observatoire de haute provence ( ohp ) in 1997 . the profiles analysed are those of the balmer , h@xmath0  and h@xmath1 , and d3 ( @xmath2 5876  ) lines . we analyse the balmer line variability and phase behaviour . we confirm these lines as highly variable , with excess emission that originates from extended structures and exhibits clear signs of rotational modulation . we have described the line profiles for two distinct states of activity , using different modelling approaches . similar techniques were applied to the d3 spectra . a large flare event lasting several days was detected in both balmer lines and in the d3 line . the energy released during this flare in h@xmath0  is of the order of 10@xmath3  erg , making it the largest h@xmath0  flare reported on a cool star . our results confirm the extreme complexity of the circumstellar environment of fk comae . the activity level of this star is quite variable demanding different approaches to the line profile analysis .