Instruction: Write an abstract for the following acadamic paper.

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

Paper: it is widely believed that the quantum choromodynamics ( qcd ) exhibits a phase transition at some temperature @xmath1 , and quarks and gluons confined in the low temperature phase are liberated to form `` quark gluon plasma '' . at the beginning of 2000 , cern reported that the qgp state had been created in the heavy ion collision experiment@xcite . in this experiment , the @xmath0 suppression @xcite is regarded as a key signal of qgp formation . more reliable inspection will be performed in rhic project at bnl .    on the theoretical side , in spite of various approaches , we are still far from definite understanding of charmonium properties near the transition and the fate of the hadronic states in the plasma phase . we investigate this problem using lattice qcd , which enables us to incorporate the nonperturbative effect of qcd . the difficulties in studying the heavy quarkonium state at finite temperature can be categorized into two classes : ( i ) difficulties to obtain detailed information on the heavy quarkonium correlators at finite temperature , and ( ii ) analyses of the obtained correlators . although it is difficult to overcome the latter , the former is possible to be solved by employing the techniques developed in recent years .    at high temperature , sufficient information on the temporal correlator requires high resolution in the temporal direction . to achieve large temporal cutoff with limited computational resources , we employ the anisotropic lattice , on which the temporal lattice spacing is finer than the spatial one @xcite . calculations with a heavy quark need further implementation . the standard lattice action contains large systematic uncertainty for the quark mass larger than or comparable to the lattice cutoff . this problem was reexamined with the light of effective theoretical interpretation of the wilson ( and clover ) quark action by fnal group @xcite . their approach is well - matched with to the anisotropic lattice , since it naturally introduces the anisotropy in the spatial and the temporal hopping parameters . we employ this approach to describe the heavy quark on present lattices , which have not sufficiently large spatial cutoff compared with the charm quark mass . the latter problems , ( ii ) , are much more difficult to be overcome . this is because one is inevitably enforced to extract the spectral properties from the data at the short temporal distance , where high frequency component of the dynamics is significant . at present , there is no well - established way of lattice simulations to attack the spectroscopy at @xmath3 . for this reason , lattice studies of hadrons at finite temperature have been argued on the spatial correlation ( screening mass ) @xcite . recently , qcd - taro collaboration ( including three of us ) analyzed the space - time structure of the mesonic correlators in the temporal direction to examine the temperature effect on them @xcite . they caught a sign of the bound state even in the deconfined phase , as well as an indication of the chiral symmetry restoration . their strategy seems one of the best approaches also to the present subject , and is expected to give us novel information on the fate of charminium near the phase transition . therefore we proceed the analysis with the following steps .    at @xmath4 we construct good meson operators which has large overlap with the state of interest . temperature effects are observed in the meson correlator between the operators constructed in the previous step . we define the `` wave function '' in the coulomb gauge at finite temperature , and observe its @xmath5-dependence below and above @xmath1 to discuss the possibility of the bound state at @xmath6 . the quark action we use is in the following form . @xmath7 @xmath8 + where @xmath9 . the bare anisotropy @xmath10 of the fermion field appears in the ratio of the spatial and the temporal hopping parameters , @xmath11 . the wilson parameter in this form is @xmath12 , where @xmath13 is the renormalized anisotropy : @xmath14 . the hopping parameters are related to the bare quark mass as @xmath15 . we also define @xmath16 with @xmath17 as the parameter which plays the same role as @xmath16 on the isotropic lattice . as the field strength @xmath18 , we use the standard clover - leaf definition . at the tree level , the coefficients of the clover terms , @xmath19 and @xmath20 , are unity . this action is obtained by generalizing the form in @xcite to the anisotropic lattice , and by setting @xmath12 . this choice of the wilson term is adopted so that the action restores the axial interchange symmetry in the physical unit as @xmath21 . therefore we can use the same form in the light quark region as the @xmath22 improved wilson quark action @xcite . this is the reason that our form of the action is different from the form used in @xcite .    on the anisotropic lattice , we need additional tuning of the parameters called as `` calibration '' , to ensure the renormalized anisotropy is same for the gauge and the fermion fields @xcite . we define the fermionic anisotropy @xmath23 using the dispersion relation of the pseudoscalar and the vector mesons and tune the value of @xmath10 so that @xmath23 coincides with the gauge field anisotropy @xmath13 . we incorporate the mean - field improvement @xcite : @xmath24 ( @xmath251,2,3 ) and @xmath26 with the mean - field values of the spatial and the temporal link variables , @xmath27 and @xmath28 defined in the landau gauge . numerical calculation is done on the lattice of sizes @xmath29 where @xmath30 and @xmath31 , with the standard wilson gauge action at the coupling @xmath32 and the bare anisotropy @xmath33 , in the quenched approximation @xcite . on @xmath34 lattice , the renormalized anisotropy and the spatial lattice spacing from the string tension are determined as @xmath35 and @xmath36 gev . @xmath37 , @xmath38 and @xmath31 lattices correspond to the temperatures @xmath39 0.93 @xmath1 , 1.15 @xmath1 and 1.5 @xmath1 , respectively . we calculate the meson correlators at @xmath4 with the three sets of the parameters ( @xmath16 , @xmath10 ) listed in table  [ tab : params ] . .the quark parameters and the meson masses in the physical unit . the error of the meson masses include the statistical error of the scale @xmath40 . [ cols="^,^,^,^,^",options="header " , ]     each value of @xmath10 is determined by the calibration procedure mentioned above . observing the resultant meson masses summarized in table [ tab : params ] , ( @xmath16 , @xmath10 ) @xmath41 ( 0.1013 , 3.703 ) roughly corresponds to the charm quark . in the successive analysis , we use this set of parameters .    to optimize the meson correlator at @xmath4 , we employ the variational analysis . this is implemented by diagonalizing the correlator matrix @xmath42 the operator @xmath43 is defined with the smearing function @xmath44 in the coulomb gauge as @xmath45 where the @xmath46 matrix @xmath47 specifies the quantum number of the meson . as the smearing function @xmath48 , we use the eigenfunction obtained by solving the schrdinger equation for the @xmath49 state , @xmath50y(r ) & = & e y(r ) ,    \nonumber \\   y(r)&=&r\varphi(r),\end{aligned}\ ] ] where @xmath51 is the static quark potential obtained on the same lattice , and @xmath52 [ gev ] is used . [ fig : vanal ]    figure  [ fig : vanal ] shows the effective masses of the diagonalized correlators with two largest eigenvalues ( smallest masses ) . both correlators show plateaus beyond @xmath53 , though the statistical fluctuation is large for the correlator corresponding to the first excited state . we fit these data at @xmath54@xmath55 and find the meson masses in 2@xmath49 states as @xmath56 and @xmath57 in the lattice unit . although these give 1@xmath49 - 2@xmath49 splitting smaller than the experimental value , considering the large systematic uncertainty on this coarse lattice , they are not inconsistent results . we regard the operator corresponds to the ground state in this analysis as the optimized meson operator at @xmath4 . we start the argument at @xmath3 with the temperature dependence of the correlator between the optimized operators at @xmath4 difined in the previous section . figure  [ fig : ep ] shows the effective mass plot at finite temperatures in the vector channel . below @xmath1 ( @xmath58 ) , the plateaus of the effective mass seem to appear . although the effective mass of the ground state is slightly larger than that at @xmath4 , it is difficult to identify the plateau precisely and determine the mass quantitatively with present statistics . more detailed analysis with higher statistics may open a stage to discuss the potential mass shift of charmonium near to the @xmath1 @xcite .      above @xmath1 ( @xmath59 and @xmath31 ) , the correlators do not show any clear plateau in the effective mass plot . they decrease more rapidly than expected from the behavior at @xmath58 in both of the correlators correspond to the ground and the first excited states at @xmath4 . the ps channel shows similar feature , although the decrease of the effective mass is slightly milder . these behavior at least signal significant change in the nature of the correlators when the system crosses @xmath1 . in @xcite , the ( effective ) masses increase as @xmath60 in the pseudoscalar and the vector channels . the observed behavior in present work , however , shows qualitatively different nature of the correlators .    to discuss the spatial correlation between quark and antiquark at @xmath3 , we observe the `` wave function '' normalized at the spatial origin , @xmath61 @xmath62 if there is no bound state , like as with the free quark propagators , this wave function @xmath63 broadens with @xmath5 for any source smearing function . in contrast , the existence of the stable shape of @xmath63 give us a hint on the existence of the bound state . figure  [ fig : wavefunc ] shows the @xmath5-dependence of the wave function normalized at the spatial origin in the vector channel at @xmath6420 and 12 . the source function of these correlators are taken to be wider than the optimized one , to show clearly that @xmath63 keeps narrow shape . at the both @xmath65 , the wave function gradually decrease as @xmath5 , in contrast to the case with the free quark poropagators shown together . this indicates that the quark and antiquark tend to stay together even in the deconfined phase , up to 1.5 @xmath1 . these two observations are interesting , and at the same time puzzling . obviously , the phase transition causes significant change in the charmonium correlators , especially in the vector channel . on the other hand , the result of the wave function suggests the persistence of the bound states . one natural picture is that the mesonic spectral function still have a peak with rather large width above @xmath1 . this explains both of our observations . direct extraction of the spectral function from the lattice data @xcite may give us further information . for the definite understanding of the fate of chamonium at the phase transition , we need more systematic studies and development of procedures . the simulation has been done on intel paragon xp / s and nec hsp at insam , hiroshima university . this work is supported by the grant - in - aide for scientific research by monbusho , japan ( no.10640272,no.11440080 )
Abstract:

Answer: we study hadron properties near the deconfining transition in the finite temperature lattice qcd . this paper focus on the heavy quarkonium states , such as @xmath0 meson . we compare the meson correlators above and below @xmath1 and discuss the possibility of the @xmath2 bound state by observing the wave function .


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: the calculus involving arbitrary orders of derivatives and integrals is called the fractional calculus . recently , the fractional calculus has found many applications in several areas of science and engineering @xcite . the nonlocal property of the fractional derivatives and integrals is used to model the processes with memory effect @xcite . for example , the fractional derivatives are used to model more appropriate the dynamics of the non - conservative systems in the hamilton , lagrange and nambu mechanics @xcite . the continuous but non - differentiable functions admit the local fractional derivatives @xcite . the local fractional derivative give the measure on fractal sets . consequently , recently the @xmath0-calculus on the fractal subset of real line and fractal curves is built as a framework @xcite . fractal analysis is established by many researchers by using different methods @xcite . using @xmath0-calculus the newton , lagrange and hamilton mechanics were built on fractal sets @xcite . also , the schrdinger s equation on fractal curve was derived in @xcite . motivated by the above mentioned interesting results , in this work , we define the non - local derivative on fractal sets . these new derivatives can be successfully used to derive new mathematical models on fractal sets involving processes with memory . + we organize our manuscript as follows : + in section [ se-2 ] , we gave a brief exposition of @xmath0-calculus and defined fractal gamma and beta functions . in section [ se-3 ] we defined the non - local derivative on fractals as a generalized riemann - liouville and caputo fractional derivatives .  in section [ se-4 ] , mittag - leffler function and non - local laplace fractional on fractal sets are introduced . we solved the non - local differential equations on fractal using the suggested methods . section [ se-5 ] was devoted to our conclusion . in this section , we review the @xmath0-calculus @xcite .      the fractal geometry is the geometry of the real world @xcite . fractal shape is the object with the fractional dimension and the self similarity property @xcite . in the seminal paper , parvate and gangal have established a calculus on fractals which similar to riemann integration . the suggested framework become a mathematical model for many phenomena in fractal media @xcite . we recall that the triadic cantor set is a fractal that can be obtained by an iterative process . in figure [ cantor1f ] we show the triadic cantor set @xcite . the integral staircase function for the triadic cantor set is defined as @xcite . @xmath1 where @xmath2 is the @xmath3-dimension of triadic cantor set . in figure [ hddho12 ] we plot the integral staircase function for triadic cantor set . the definitions of @xmath0-limit , @xmath0-continuity and @xmath0-integration are given in the ref . the @xmath0-differentiation is denoted by @xmath4 and it is defined as @xmath5 if the limit exists @xcite . + * definition 1 . * the gamma function with the fractal support is defined as @xmath6 where @xmath7        * definition 2 . * the fractal beta function on the fractal set is defined as follows @xmath8 which is called two - parameter @xmath9 fractal integral , where @xmath10 and @xmath11 . + in the following we present some properties of fractal beta function . + 1 ) the fractal beta function has a symmetry @xmath12 . since , we have @xmath13 using the transformation @xmath14 , we conclude that @xmath15 2 ) using the transformation @xmath16 , we get following form for the fractal the beta function @xmath17 3 ) the beta fractal function is related to the fractal gamma function as @xmath18 * proof : * we have @xmath19 transforming to polar coordinates @xmath20 we obtain @xmath21 thus , the proof is completed . in this section , we define the non - local derivative for the functions with fractal support . + * definition 3 . * if @xmath22 $ ] ( @xmath2-order differentiable function on @xmath23 $ ] ) and @xmath24 then we have @xmath25 where if @xmath26 then we have fractal integral whose order is equal the dimension of the fractal , and @xmath27 are called the analogous left sided and the right sided riemann - liouville fractal integral of order @xmath28 . + * definition 4 . * let @xmath29 , then analogous left and right riemann - liouville fractal derivative are defined as follows : @xmath30 * definition 5 . * let @xmath31 $ ] , then the analogous left sided caputo fractal derivative is defined by @xmath32).\ ] ] also , the analogous right sided caputo fractal derivative has the form @xmath33 now , we give some important relations , namely @xmath34 * proof : * using the eq . ( [ xsaz125 ] ) we conclude @xmath35 let us consider @xmath36 therefore , @xmath37 while @xmath38 . as a result we obtain @xmath39 substituting eqs . ( [ bgthcxd963 ] ) and ( [ cdesx44 ] ) in eq . ( [ m258jyhnbgt ] ) we conclude that @xmath40 then , we have @xmath41 in view of eq . ( [ zaws ] ) we derive @xmath42 applying eq . ( [ vcxdzsax ] ) we get @xmath43 now , we consider following formula @xmath44 * proof : * by rewriting the eq . ( [ xsamjuza ] ) we get @xmath45 utilizing the eq . ( [ xsaza ] ) we conclude @xmath46 now , we write some important composition relations , namely @xmath47 * proof : * using the definitions we get @xmath48 applying , n - times integration by part it leads to @xmath49 the similar proof works for the following formulas @xmath50 in figures [ 5414 ] and [ 544 ] we compared the non - local standard derivative versus non - local fractal derivative and the generalized fractal integral .     and @xmath51 and their non - local derivative @xmath52 and @xmath53 , respectively . ]     and @xmath51 and their non - local integral @xmath54 and @xmath55 , respectively . ] in this section , we suggest the mathematical tools for solving the non - local fractal differential equations . now , we define the gamma function for the fractal calculus that will be used in non - local calculus on fractals . it is well known that the exponential function has important role in the theory of standard differential equation . the generalized exponential function is called the mittag - leffer function and plays an important role for fractional differential equations @xcite . + * definition 6 . * the generalized two parameter @xmath56 mittag - liffler function on fractal @xmath57 with @xmath2-dimension is defined as @xmath58 in the special case we have the following results , namely @xmath59      the laplace transformation is a very useful tool for solving standard linear differential equation with constant coefficients . the generalized laplace transformation is applied to solve the fractional differential equations . thus , in this section , we generalized the laplace transformation for the function with fractal support which is utilized to solve the non - local differential equation on the fractal set @xcite . + * definition 7 . * laplace transformation for the function @xmath60 is denoted by @xmath61 and it is defined as @xmath62=    \int_{s_{f}^{\alpha}(0)}^{s_{f}^{\alpha}(\infty)}f(x ) e^{-s_{f}^{\alpha}(s)s_{f}^{\alpha}(x)}d_{f}^{\alpha}x.\ ] ] now , we give the fractal laplace transformation of some functions . if we define the fractal convolution of two function @xmath60 and @xmath63 as follows : @xmath64 the fractal laplace transformation of power function of @xmath65 is @xmath66    \int_{s_{f}^{\alpha}(0)}^{s_{f}^{\alpha}(\infty)}s_{f}^{\alpha}(x)^{\beta}e^    { -s_{f}^{\alpha}(s){s_{f}^{\alpha}(x)}}d_{f}^{\alpha}x=    \frac{\gamma_{f}^{\alpha}(1+\beta)}{s^{\beta+1}}.\ ] ] * lemma : * the laplace transformation of the non - local fractal riemann - liouville integral is given by @xmath67=    \frac{\mathcal{f}_{f}^{\alpha}(s_{f}^{\alpha}(s ) ) } { s_{f}^{\alpha}(s)^{\beta}}.\ ] ] * proof : * the laplace transform of the fractal riemann - liouville integral is @xmath68 =      \mathcal{l}_{f}^{\alpha}\left [ \frac{1}{\gamma_{f}^{\alpha}(\beta)}\int_{s_{f}^{\alpha}(0)}^ { s_{f}^{\alpha}(x)}\frac{f(t)}{(s_{f}^{\alpha}(x)-s_{f}^{\alpha}(t))^{\alpha-\beta } } d_{f}^{\alpha}t\right].\ ] ] using the eqs.([zsewnja1485 ] ) and ( [ 14nhytgb ] ) we arrive at @xmath69     & = \frac{1}{\gamma_{f}^{\alpha}(\beta)}\mathcal{f}_{f}^{\alpha}(s_{f}^{\alpha}(s ) )    \mathcal{l}_{f}^{\alpha}[s_{f}^{\alpha}(x)^{\beta-1}],\nonumber    \\&=    \frac{1}{\gamma_{f}^{\alpha}(\beta)}\mathcal{f}_{f}^{\alpha}(s_{f}^{\alpha}(s ) )    \frac{\gamma_{f}^{\alpha}(\beta)}{s_{f}^{\alpha}(s)^\beta},\nonumber\\&=    \frac{\mathcal{f}_{f}^{\alpha}(s_{f}^{\alpha}(s ) ) } { s_{f}^{\alpha}(s)^{\beta}}.\end{aligned}\ ] ] the fractal laplace transform of the non - local fractal riemann - liouville derivative of order @xmath70 is given by @xmath71 where @xmath72 + 1 $ ] . the fractal laplace transform of the non - local fractal caputo derivative of order @xmath70 is given by @xmath73 where @xmath74)$ ] . in this section , we solve some illustrative examples . + * example 1 . * consider the following linear fractal equation @xmath75 with the initial condition @xmath76 where @xmath77 is cantor set dimension . by applying @xmath78 on the both sides of the eq . ( [ xz ] ) we obtain @xmath79    ) on the real - line and cantor set . ] * example 2 . * consider a linear fractal differential equation @xmath80 with initial condition as @xmath81 by applying @xmath82 on the both sides of the eq . ( [ n47712hytgb ] ) we arrive at @xmath83    ) on the real - line and cantor set . ]    in figures , [ 8945laqwqaqok44 ] and [ 8945lok4xxzsa4 ] we plot the solutions of eqs.([swz5xc9 ] ) and ( [ xde ] ) , respectively . + * example 3 . * consider a linear differential equation @xmath84 with the following initial condition , namely @xmath85 by inspection , the solution for the eq . ( [ juyh ] ) becomes @xmath86 in figure [ 8945lok44 ] we sketched the solution of eq . ( [ juyh ] ) on the cantor set and real - line .    ) on the real - line and cantor set . ]    * example 4 . * we examine the following non - local differential equation on a fractal subset of real - line , namely with the following initial condition @xmath87 for solving eq . ( [ bds128 ] ) we apply the fractal laplace transformation on both side of it and we get @xmath88 after some calculations we obtain @xmath89 by computing the inverse fractal laplace transform we conclude @xmath90 * remark : *  the figures [ 8945lok4xxzsa4 ] , [ 8945laqwqaqok44 ] , and [ 8945lok44 ] show that the solution of eqs . ( [ xde ] ) , ( [ swz5xc9 ] ) and ( [ juyh ] ) leads to the standard non - local fractional cases when @xmath91 , respectively . in this work , we defined new non - local derivatives on fractal sets . these new type of non - local derivatives can describe better the dynamics of complex systems which possess memory effect on a fractal set . four illustrative examples were solved in detail . finally , one can recover the standard non - local fractional cases when put @xmath91 . 999 vladimir v. uchaikin , _ fractional derivatives for physicists and engineers _ vol . 1 background and theory . application , springer , berlin , 2013 . d. baleanu , k. diethelm , e. scalas , juan j. trujillo , _ fractional calculus models and numerical methods . ser . on complexity , nonlinearity and chaos _ , world scientific , 2012 . stefan g. samko , anatoly a. kilbas , oleg i. marichev , _ fractional integrals and derivatives theory and applications _ , gordon and breach , new york , 1993 . i. podlubny , _ fractional differential equations_. academic , new york , 1999 . skwara urszula , et al . , applications of fractional calculus to epidemiological models , _ aip conference proceedings - american institute of physics_. * 2012 * , 1479 ( 1 ) , 1339 - 1342 . bruce j. west , m. bologna , p. grigolini , _ physics of fractal operators _ , new york , springer , 2003 alireza k. golmankhaneh , _ investigations in dynamics : with focus on fractional dynamics _ , lap lambert academic publishing , saarbrucken , 2012 . d. baleanu , alireza k. golmankhaneh , r. nigmatullin , ali k. golmankhaneh , fractional newtonian mechanics , _ cent . _ , * 2010 * , 8(1 ) , 120 - 125 . d. baleanu , alireza k. golmankhaneh , ali k. golmankhaneh , fractional nambu mechanics , _ int . j. theor _ , * 2009 * , 48(4 ) , 1044 - 1052 . d. baleanu , alireza k. golmankhaneh , ali k. golmankhaneh , and raoul r. nigmatullin , newtonian law with memory , _ nonlinear dyn . _ , * 2010 * , 60(1 - 2 ) , 81 - 86 . kiran m. kolwankar , anil d. gangal , fractional differentiability of nowhere differentiable functions and dimensions , chaos : an interdisciplinary _ j. nonlinear sci . _ , * 1996 * , 6(4 ) , 505 - 513 . a. parvate , anil d. gangal , calculus on fractal subsets of real - line i : formulation , _ fractals _ , * 2009 * , 17(01 ) , 53 - 81 . a. parvate , anil d. gangal , calculus on fractal subsets of real - line ii : conjugacy with ordinary calculus , _ fractals _ , * 2011 * , 19(03 ) , 271 - 290 . benoit b. mandelbrot , _ the fractal geometry of nature _ , freeman and company , 1977 . k. falconer , _ techniques in fractal geometry _ , john wiley and sons , 1997 . j. kigami , _ analysis on fractals _ , vol . 143 , cambridge university press , 2001 . xiao - jun yang , _ advanced local fractional calculus and its applications _ , world science , new york , ny , usa 2012 . alireza k. golmankhaneh , ali m. yengejeh , d. baleanu , on the fractional hamilton and lagrange mechanics , _ int . _ , * 2012 * 51(9 ) , 2909 - 2916 . alireza k. golmankhaneh , ali k. golmankhaneh , d. baleanu , lagrangian and hamiltonian mechanics on fractals subset of real - line , _ int . j. theor . phys . _ * 2013 * , 52(11 ) , 4210 - 4217 . alireza k. golmankhaneh , ali k. golmankhaneh , d. baleanu , about maxwell s equations on fractal subsets of @xmath92 , _ cent . j. phys . _ * 2013 * , 11 ( 6 ) , 863 - 867 . alireza k. golmankhaneh , ali k. golmankhaneh , d. baleanu , about schrodinger equation on fractals curves imbedding in @xmath92 , _ int . j. theor _ , * 2015 * , 54 ( 4 ) , 1275 - 1282 . hari m. srivastava , alireza k. golmankhaneh , d. baleanu and xiao - jun yang , local fractional sumudu transform with application to ivps on cantor sets , _ abstr . * 2014 * article i d 620529 , 1 - 7 .
Abstract:

Answer: in this manuscript we introduced the generalized fractional riemann - liouville and caputo like derivative for functions defined on fractal sets . the gamma , mittag - leffler and beta functions were defined on the fractal sets . the non - local laplace transformation is given and applied for solving linear and non - linear fractal equations . the advantage of using these new nonlocal derivatives on fractals subset of real - line lies in the fact that they are used for better modelling of processes with memory effect . * keyword : * fractal calculus ; triadic cantor set ; non - local laplace transformation ; memory processes ; generalized mittag - leffler function ; generalized gamma function ; generalized beta function


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: einstein s general relativity enables us to come up with a testable theory of the universe . hubble@xcite first found that the galaxies are receding from us not long after the birth of general relativity . the hubble s observation indicates that our universe is expanding . over the past decade , two groups @xcite observing supernovae reported that the luminosity distance can not be explained by a matter dominated universe . if one accepts the convincible assumption of homogeneity and isotropy of the universe which is approximately true on large scales , then the general relativity tells us that we now live in a dark energy dominated universe and the universe is accelerated expanding . dark energy , which has the property of negative pressure , is different from the classical matter and the found particles . a great amount of models have been proposed to study the possible candidates of dark energy and their dynamics @xcite . the most famous and acceptable candidate is the cosmological constant . however , the magnitude of cosmological constant @xmath0 is much smaller than the energy density of vacuum in quantum field theory . the theory of dark energy @xcite dominates the modern cosmology in the past decade . this situation is similar to the raise of the theory of ether , which has been considered as direct evidence of the aberration of starlight , an important astronomical effect known since eighteenth century . although the rapid progress in technology makes the astronomical observations more and more accurate , up to now , there is no direct evidence indicate what the dark energy it is . since the theory of dark energy is contrived , which requires fine tuning and apparently can not be tested in the laboratory or solar system , several modified theories of general relativity have been developed as the alternative source of cosmological acceleration @xcite . einstein first used the riemann geometry to describe the theory of gravitation . in this letter , we suppose that the spacetime in large scale may be described by other geometry instead of riemann geometry . to preserve redeeming feature of general relativity , this geometry must take riemann geometry as its special case . fortunately , paul finsler proposed a natural generation of riemann geometry - finsler geometry . finsler geometry as a more general geometry could provide new sight on modern physics . it is of great interest for physicists to investigated the violation of lorentz symmetry @xcite . an interesting case of lorentz violation , which was proposed by cohen and glashow@xcite , is the model of very special relativity ( vsr ) characterized by a reduced symmetry sim(2 ) . in fact , gibbons , gomis and pope@xcite showed that the finslerian line element @xmath1 is invariant under the transformations of the group disim@xmath2 . in the framework of finsler geometry , modified dispersion relation has been discussed@xcite . also , the model based on finsler geometry could explain the recent astronomical observations which einstein s gravity could not . a list includes : the flat rotation curves of spiral galaxies can be deduced naturally without invoking dark matter @xcite ; the anomalous acceleration@xcite in solar system observed by pioneer 10 and 11 spacecrafts should correspond to finsler - randers space @xcite ; the secular trend in the astronomical unit@xcite and the anomalous secular eccentricity variation of the moon s orbit@xcite should be correspond to the effect of the length change of unit circle in finsler geometry@xcite .    in this letter , we present a matter dominated navigation model . the influence of a possible drift ( wind ) in the navigation cosmolgical model makes the expanding universe  accelerated " . it is remarkable that the navigation cosmological model is described in the framework of finsler geometry . we find that the predictions of the navigation cosmological model could account for the observations of riess and perlmutter@xcite on distant supernovae . instead of defining an inner product structure over the tangent bundle in riemann geometry , finsler geometry is based on the so called finsler structure @xmath3 with the property @xmath4 for all @xmath5 , where @xmath6 represents position and @xmath7 represents velocity . the finsler metric is given as@xcite @xmath8 in finsler manifold , there exists a unique linear connection  -  the chern connection@xcite . it is torsion freeness and almost metric - compatibility , @xmath9 where @xmath10 is the formal christoffel symbols of the second kind with the same form of riemannian connection , @xmath11 is defined as @xmath12 and @xmath13 is the cartan tensor ( regarded as a measurement of deviation from the riemannian manifold ) . in terms of chern connection , the curvature of finsler space is given as @xmath14 where @xmath15 . the notion of ricci tensor in finsler geometry was first introduced by akbar - zadeh@xcite @xmath16 where @xmath17 . and the scalar curvature in finsler geometry is given as @xmath18 . in standard cosmology , following the cosmological principle , one gets the friedmann - robertson - walker ( frw ) metric@xcite . in another word , the spatial part of the universe is a constant sectional curvature space . here comes our major assumption , the gravity in large scale should be described in terms of finsler geometry . in light of the cosmological principle , the finsler structure of the universe should be written in such form @xmath19 where @xmath20 is scale factor with cosmic time @xmath21 , the structure @xmath3 is a constant flag curvature space . by making use of the geometrical terms we mentioned above , one obtains the components of ricci tensor @xmath22 where the dot denotes a derivative with respect to @xmath21 . the gravitational field equation in finsler space should reduce to the einstein s gravitational field equation while the finsler space reduce to the riemannian space . thus , the symmetrical tensor @xmath23 should involve in the gravitational field equation in finsler space . in general , finsler space is an anisotropic space . therefore , it has less killing vectors than the riemannian space , and it breaks the lorentz symmetry@xcite . it means the angular momentum in finsler space is not a conservative quantity . this fact implies that the energy momentum tensor is not symmetrical in finsler space . for these reasons , the gravitational field equation in finsler space should be taken in such form @xmath24 where @xmath25 is an asymmetrical tensor and @xmath26 are symmetrical part and asymmetrical part of energy momentum tensor respectively . this general form is agree with the result of asanov , its gravitational field equation contains the asymmetrical term in finsler space of landsberg type@xcite . dealing with field equation in the anisotropic and inhomogeneous cosmology is not a simple staff , and it is hard to find an exact solution of gravitational field equation while its energy momentum tensor involves the non diagonal terms@xcite . at first glance , we just deal with the symmetrical part of field equation ( [ field eq ] ) @xmath27 by making use of the equations ( [ ric0 ] ) and ( [ rici ] ) , the solution of equation ( [ field eqn ] ) is the same with the einstein s field equation deduced by frw metric . taking the energy - momentum tensor @xmath28 to be the form of perfect fluid , one can see that the evolution of scale factor @xmath20 is the same with the riemannian case . however , since the luminosity distant is related to the spatial geometry of the universe , one may expect that it is different from the riemannian luminosity distant . unlike riemann space , a complete classification of the constant flag curvature spaces remain an unsolved problem . however , by making use of the zermolo navigation on riemannian space , bao _ al_. @xcite gave a complete classification of finsler - randers space @xcite of constant flag curvature . the zermelo navigation problem @xcite aims to find the paths of shortest travel time in a riemannian manifold ( @xmath29 ) under the influence of a drift (  wind " ) represented by a vector field @xmath30 . in standard cosmology , our universe is very flat now . we may imagine that the universe is a flat riemannian manifold with flat friedmann metric @xmath31 and it influenced by the  wind " @xmath30 . the relation between the riemannian manifold ( @xmath29 ) and the randers metric @xmath3 is @xmath32 where @xmath33 and @xmath34 . one should notice that there is a map between the riemannian space which influence by the  wind " and the randers - finsler space@xcite . it means the effect of  wind " already accounted in the gravitational field equation in finsler space . the theorem@xcite of the classification supplies an interesting case where the randers metric @xmath3 has constant flag curvature @xmath35 : for @xmath36 and @xmath37 is flat . and @xmath38 satisfies the constraint @xmath39 , @xmath40 denotes lie differentiation . we set the vector field to be radial @xmath41 . then the flag curvature of randers metric is @xmath42 . after doing coordinate change and taking spherical coordinate , we get the randers metric as @xmath43 after scale change @xmath44 , the metric @xmath3 changes as @xmath45 and the cosmic scale factor changes as @xmath46 . taking the spatial part of the finsler structure ( [ finsler s ] ) to be of the form ( [ randers ] ) , we get metric of the universe in the framework of finsler geometry , and the space curvature of the universe is @xmath47 . one can deduce directly from such metric that the relation between the redshift @xmath48 and the scale factor @xmath20 is the same with the riemannian case @xmath49 where the subscript zero denotes the quantities given at the present epoch . since the non - radial part of the metric ( [ randers ] ) is the same with frw metric , the luminosity distant in the navigation cosmological model is given as @xmath50 . however , the proper distant is not the case , and it is given as @xmath51 the light traveling along the radial direction satisfies the geodesic equation @xmath52 then , we obtain @xmath53 the equations ( [ field eqn ] ) and ( [ dp ] ) can be solved . supposing the universe is matter dominated ( no more cosmological constant and dynamical dark energy ) , we obtain @xmath54 where @xmath55 is the hubble constant and @xmath56 is the density parameter for matter and satisfies @xmath57 substituting the equation ( [ dp1 ] ) into ( [ dp0 ] ) , we obtain the relation between luminosity distant and redshift in the navigation cosmological model @xmath58 versus the redshift @xmath48 for the navigational cosmological model . the data comes from riess _ et al_.@xcite and hubble space telescope ( hst)@xcite . and the value of hubble constant is set as @xmath59 . ]    here , we present numerical result for the relation between luminosity distant and redshift in the framework of the navigation cosmological model . the best fit curve is shown in fig.1 with the matter density parameter taken as @xmath60 . and the average value of chisquare is @xmath61 . thus , our prediction could account for experiment data given by riess _ et al_.@xcite . it should be noticed that the hubble constant @xmath55 we took is different from the hubble constant @xmath62 measured in flat space . the relation for @xmath63 and @xmath64 is @xmath65 by making use of the value of the hubble constant measured in flat space @xmath66@xcite , we have @xmath67 this geometrical parameter @xmath68 represents an  accelerated " effect provided by the vector field @xmath30 . our letter initiates an exploration of the possibility that the empirical success of the observations of type ia supernovae can be regarded as the influence of a navigational wind . the particles move on riemnnian manifold and influenced by a vector field ( the  wind " ) which proportion to the curvature @xmath69 of the space of the universe . its geodesic indeed is a finslerian geodesic . thus , the observations of riess _ et al_.@xcite on distant supernovae may be explained by the effect of the  wind " . our numerical result could account for astronomical observations . at last , we point out that the age of the universe is about 9.76gyr in our model . this is contradicted with the age ( @xmath70 ) of globular clusters in the milkyway@xcite . in the navigation cosmological model , we only involve the radial  wind " . the non - radial  wind " should be taken into account in future work , we expect that the effect of the non - radial  wind " may supply us a reasonable age of the universe . 999 e. hubble , _ proceedings of the national academy of sciences of the united states of america _ * 15 * , 168 ( 1929 ) . a. g. riess _ et al_. , astron . j. * 116 * , 1009 ( 1998 ) ; astron . j. * 117 * , 707 ( 1999 ) . s. perlmutter , _ et al_. , astrophys j. * 517 * , 565 ( 1999 ) . e. j. copeland , m. sami and s. tsujikawa , int . d * 15 * , 1753 ( 2006 ) . t. padmanabhan , phys . rept . * 380 * , 235 ( 2003 ) . s. bludman , arxiv : astro - ph/0605198 . d. colladay and v. a. kostelecky , phys . d55 , 6760 ( 1997 ) ; v. a. kostelecky , phys . d * 69 * , 105009 ( 2004 ) . a. g. cohen and s.l . glashow , phys . lett . * 97 * , 021601 ( 2006 ) . g. w. gibbons , j. gomis and c. n. pope , phys . d * 76 * , 081701 ( 2007 ) . f. girelli , s. liberati and l. sindoni , phys . d * 75 * , 064015 ( 2007 ) . z. chang and x. li , phys . b * 663 * , 103 ( 2008 ) . z. chang and x. li , phys . b * 668 * , 453 ( 2008 ) . j. d. anderson _ et al_. , phys . lett . * 81 * 2858 , ( 1998 ) ; phys . d * 65 * , 082004 ( 2002 ) ; mod . a * 17 * , 875 ( 2002 ) . x. li and z. chang , arxiv : gr - qc/0909.3713 . g. a. krasinsky and v. a. brumberg , celest . . astrn . * 90 * , 267 ( 2004 ) . e. m. standish , proc . iau colloq . * 196 * , 163 ( 2005 ) . j. g. williams and d. h. boggs , in proceedings of 16th international workshop on laser ranging , ed . s. schillak , ( space research centre , polish academy of sciences ) , 2009 . x. li and z. chang , arxiv : gr - qc/0911.1890 . d. bao , s. s. chern and z. shen , _ an introduction to riemann  finsler geometry _ , graduate texts in mathmatics * 200 * , springer , new york , 2000 . s. s. chern , sci . rep . tsing hua univ . a * 5 * , 95 ( 1948 ) ; or selected papers , vol . ii , 194 , springer 1989 . h. akbar - zadeh , acad . ( 5 ) * 74 * , 281 ( 1988 ) . s. weinberg , _ gravitation and cosmology : principles and applications of the general theory of relativity _ , wiley , new york , 1972 . x. li , z. chang and x. h. mo , arxiv : hep - th/1001.2667 . g. s. asanov , reports on mathematical physics * 59 * , 111 ( 2007 ) . c. w. misner , k , s , thorne and j. a. wheeler , _ gravitation _ , freeman , san francisco , 1973 . d. bao , c. robles and z. shen , j. diff . geom . * 66 * , 377 ( 2004 ) ; arxiv : math/0311233 . g. randers , phys . rev . * 59 * , 195 ( 1941 ) . e. zermelo , z. angew . mech . * 11*(2 ) , 114 ( 1931 ) . g. w. gibbons , _ et al . d * 79 * , 044022 ( 2009 ) . m. kowalski _ et al_. , astrophys . j. * 686 * , 749 ( 2008 ) . d. n. spergel _ et al_. , astrophys . * 170 * , 377 ( 2007 ) . r. jimenez , p. thejll , u. jorgensen , j. macdonald and b. pagel , mnras * 282 * , 926 ( 1996 ) .
Abstract:

Answer: we investigate a matter dominated navigation cosmological model . the influence of a possible drift ( wind ) in the navigation cosmological model makes the spacetime geometry change from riemannian to finslerian . the evolution of the finslerian universe is governed by the same gravitational field equation with the familiar friedmann - robertson - walker one . however , the change of space geometry from riemannian to finslerian supplies us a new relation between the luminosity distant and redshift . it is shown that the hubble diagram based on this new relation could account for the observations on distant type ia supernovae .


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: cosmological numerical cold dark matter ( cdm ) simulations predict a specific and universal shape for the dark matter density distributions of galaxies . this was investigated in detail by navarro , frenk & white ( 1996 , 1997 ) [ nfw ] who found that the mass - density distribution in the inner parts of simulated dark - matter - only cdm halos could best be described by a @xmath2 power law . this is however not what is typically observed : the mass - density distributions of , for example , dwarfs and late - type galaxies ( which do include baryons ) are best described by an almost constant density central core @xmath3 ( de blok et al . this disagreement has become known as the `` cusp / core problem '' .    there has been much discussion of the validity of both observations and simulations . my intention here is not to give an exhaustive review of this debate . rather i want to highlight some of the concerns and criticisms of both observers and simulators , and discuss them from an observer s point of view ( and in the spirit of the sydney panel discussion i will do this in a slightly provocative manner ) . the rotation curves of lsb galaxies are some of the primary evidence for the failure of the nfw model in the centres of galaxies . how reliable are these data ? an often heard complaint is that rotation curve observers can not seem to agree amongst themselves on their own data . unfortunately , this is occasionally used as an excuse to dismiss these data . this apparent disagreement is mostly a sociological effect . a large part of it can be traced back to spin ( where i use the non - astronomical definition of the word ) and over - reliance on a few discrepant rotation curves that have usually already been rejected by the observers . systematic effects due to pointing of the telescope or beam - smearing also feature prominently . in the latter case the few curves where beam - smearing was at the time a genuine problem have received a great deal of attention ( e.g.  swaters et al . 2000 ) , in contrast with the dozen or so where the problem was insignificant ( de blok et al . 2001a , de blok & bosma 2002 ) . beam - smearing was and is a red herring . similar conclusions apply to the pointing `` problem '' . here the observed shallow slopes are thought to be caused by the spectrograph slit not being precisely centered on the galaxy . this would then cause an nfw halo to mimic a halo with a shallow profile . the problem is that one would have to miss the centre of a galaxy with an nfw halo by an average of @xmath4 or so to observe an apparently shallow profile ( de blok , bosma & mcgaugh 2003 ) . in reality the centre of a galaxy can be acquired with an accuracy of @xmath5 , an order of magnitude better . in addition , if nfw halos are really prominently present in lsb galaxies one would have to miss the centres of these galaxies by a factor of 10 more than the pointing error for _ each _ and _ every _ of the over 80 galaxies observed . if pointing errors are not the problem , but rather offsets between the dynamical and photometric centres of galaxies , then the @xmath6 value translates into a physical offset of @xmath7 kpc ( for the galaxies in de blok et al . this would mean that in many lsb and dwarf galaxies optical and dynamical centres would be separated by roughly a disk scalelength !    in other words , pointing effects , though of course present in the data to a _ small _ degree , are by no means enough to challenge the observations . as a picture can say more than a thousand words , we illustrate this by comparing rotation curves of lsb galaxies observed by different independent groups at different times with different telescopes . there is now enough data in the literature that this can be done for over a dozen lsb galaxies . space does not allow us to present all galaxies here ( see de blok et al . 2003 for more comparisons ) , but we illustrate with two examples : lsb galaxies f583 - 1 and f568 - 3 ( fig .  1 ) . the latter has even been observed by three independent groups . to avoid confusion due to mixing of h@xmath8 and hi data ( as discussed below ) we compare the raw h@xmath8 data only . it is obvious that within the error bars the data agree beautifully . many of these curves are now starting to be confirmed by high - resolution 2-d velocity fields ( e.g. , chemin et al . , this volume ) . this would be quite a fluke if systematic errors really were dominant . in summary , _ the data are consistent . _ if the overwhelming majority of observational data are consistent with each other and therefore with dark matter cores , how does this explain the often - cited statement usually attributed to swaters et al . ( 2003 ) that lsb rotation curves are consistent with nfw halos ? the swaters et al . analysis actually does not claim this in such strong terms . its strongest statement in this respect is that 75% of the galaxies in their sample are `` consistent '' with a nfw model . what this means is made clear at the end of the swaters et al . paper : _ `` even though the majority of the galaxies in this sample seem consistent with steep inner slopes , none of the galaxies halos with @xmath9 , as most galaxies are equally well or better fitted by halos with shallower density profiles or even constant density cores . '' _ how one interprets this statement depends only on whether one prefers the mode of a distribution as representative ( without the need to invoke systematic effects ) , or the extremes of that distribution along with an epicyclic range of systematics . here occam s razor has often proved to be invaluable . the rotation curves in the literature span a large range in resolution , inclination , and general quality . it is therefore worth asking whether the data may in any way be affected by the effects of bars , non - minimum disc situations , resolution effects , etc . we thus isolate the highest quality rotation curves from the total sample from de blok et al . ( 2001a ) and de blok & bosma ( 2002 ) ( the open histogram in fig .  2 [ left panel ] ) . we have applied the following cuts : ( 1 ) we remove all galaxies with inclinations @xmath10 and @xmath11 ( single hatched histogram ) ; ( 2 ) we remove all galaxies with low - quality rotation curves , i.e.  with a small number of independent data points , with large error - bars , and large asymmetries ( double hatched histogram ) ; ( 3 ) galaxies which were clearly not dominated by dark matter ( e.g.  with bulges ) were removed ; ( 4 ) as the difference between cusp and core is most clearly visible in the innermost part of the galaxies , we demand that the inner part is resolved : we retain only those galaxies that have at least 2 independent data points in the inner 1 kpc ( grey filled histogram ) . this leaves us with a final sample of 19 galaxies , out of over 80 originally . a comparison of the different stages of pruning shows no systematic differences between the distributions of the initial full sample and final restricted sample . if anything , the peak of the distribution of inner slopes @xmath8 has shifted towards a more positive value . the tail towards steep negative slopes has disappeared .    as described above , an analysis by swaters et al . ( 2003 ) finds that the slopes in their sample cover the full range @xmath12 , with a preference for shallow slopes . below we show that a large fraction of the steep slopes are actually due to insufficient resolution . we have thus subjected their data to the same quality criteria that we used for ours . the results are over - plotted in fig .  2 ( left panel ; thick line ) . the number of remaining galaxies in the swaters et al . sample is small , but it is clear that the distribution is entirely consistent with our larger restricted sample . our restricted sample has three galaxies in common with their restricted sample . we find that for two of these galaxies the measured slopes are consistent at better than @xmath13 ( ugc 731 and ugc 4325 ) . for one galaxy ( ugc 11557 ) the slopes differ significantly : they find @xmath14 , while we find @xmath15 . an inspection of the density profiles in fig . 1 of de blok et al .  ( 2001b ) and fig .  5 in swaters et al . ( 2003 ) shows that this difference in slope is entirely due to the range over which the slope is measured . the data themselves are consistent . swaters et al.(2003 ) have chosen a larger fitting range , resulting in a much steeper slope . choosing a smaller range results in identical slopes . this comparison of restricted samples shows that once a proper pruning is made to only include minimum - disc , well - resolved galaxies , most of the differences are merely a matter of semantics . the observed _ shallow _ inner mass - density slopes of lsb and dwarf galaxies do therefore reflect these galaxies potentials .      as described in de blok et al . ( 2001b ) the inversion from rotation curve to mass - profile is only valid in the limit of spherical halos consisting only of dark matter . the last assumption is known as `` minimum disc '' and has been used in all rotation curve inversions described here . the minimum disc assumption ascribes all dynamics and all mass to the dark matter component , and is thus really a `` maximum halo '' model . as such it gives strong upper limits on the measured dark matter slopes . a proper treatment of non - minimum disc models should include the influence that dark and visible matter have on each other in the disc , but to first order we can derive the rotation curve under a non - minimum disc assumption and invert this curve . this should give an indication of the change in inferred dark matter slopes when going from a one - component to a two or three component mass model . we have done this for the rotation curves in de blok et al . ( 2001a ) for which photometry and hi data were available , and have derived the halo rotation curves assuming @xmath16 . the dark matter mass - density slopes found in this way are shown in fig .  2 ( right panel ) , where they are compared with the minimum - disc slopes for the same galaxies . the peak of the slope histogram is still found near @xmath17 , showing that lsb galaxies are indeed dark matter dominated . the bulk of the distribution has shifted to more positive slopes , making the discrepancy with the predicted @xmath18 value larger . it is important to realise that in more sophisticated models one has to take into account that the halo reacts to the presence of a disc by contracting slightly . this means that the original `` cosmological '' dark matter slopes must have been _ even less steep _ than the ones currently measured . despite the overwhelming evidence that measuring rotation curves is a reliable and repeatable exercise , claims do sometimes appear in the literature suggesting the opposite . in every one of these cases though , not all aspects of the data and their limitations are taken into account , which often leads to comparing apples with oranges .    a comparison of two datasets of ugc 4325 in hayahshi et al . ( 2003 ) suffers from this . one of them , from de blok & bosma ( 2002 ) , is a _ processed ( smoothed ) h@xmath8-only _ rotation curve . the second data set is taken from swaters et al . ( 2003 ) and consists of _ raw h@xmath8 data combined with processed hi data_. hayashi et al . ( 2003 ) compare both sets and point out the different nature of the error bars in both sets ( as one would expect given the different natures of the data sets ; however , it leads them to reject the use of @xmath19 in data analysis ) . they fail to take into account that the outer points of the swaters et al . rotation curve are determined by hi , not h@xmath8 , which leads to an entirely different weighting of the data . all this leads to an inappropriate and apparently unfavourable comparison .    comparing like with like , for example by comparing the raw swaters et al . h@xmath8 data with the raw de blok & bosma h@xmath8 data we see none of these glaring disagreements ( bosma , this volume ) . similarly , a comparison of the raw de blok & bosma h@xmath8 data with the raw swaters et al  hi data ( as shown in de blok & bosma 2002 ) also shows good agreement . the lesson to be drawn from this , is that data are not as black - and - white as theoretically motivated comparisons sometimes portray . data are not just `` good '' or `` bad '' ; they may be both , depending on the situation in which they are used . it is vital to be aware of the limitations . in every large collection of data it is always possible to find the odd one out ( which is precisely why one collects large data sets ) . for lsb rotation curves , when proper comparisons are made , agreement is almost always found . this may explain why virtually all observational papers agree on the inappropriateness of the nfw fitting formula for describing lsb galaxy dynamics ( contributions by , amongst others , bolatto , chemin , zachrisson , pizella , gentile , this volume ) .      in this context , it is instructive to compare the new high - resolution simulations by white and navarro ( this volume ) and hayashi et al.(2003 ) with observations . these authors present new high - resolution simulations of dark matter halos in a @xmath1cdm universe . these are now becoming just resolved enough to start probing the regime of dwarf and lsb galaxies , however , they also show the large gap that still exists between simulations and observations . observationally , we can plot the inner mass density slope @xmath8 against resolution @xmath20 of a rotation curve . such a diagram ( fig . 3 ) clearly shows that at high resolutions ( @xmath21 kpc ) the core and cusp model can be clearly distinguished from each other , but that that at low resolutions @xmath22 kpc the cusp model and the core model have identical slopes ( de blok et al . 2001b ) . using fig .  4 in hayashi et al . ( 2003 ) we have computed the values @xmath20 for their simulations and show the corresponding slopes in fig .  3 . the smallest scales these new simulations resolve are @xmath23 kpc , which is exactly the scale below which the core / cusp discrepancies become apparent . for now , the current models fall solidly in the range where the core and cusp models make identical predictions .    going back to the original nfw papers ( nfw 1996,1997 ) one sees that these older simulations also did not resolve scales smaller than 1 kpc . however , the discussion in these papers sometimes make it all too easy to forget this limitation . for example , nfw ( 1996 ) states ` the density profiles of cdm halos _ of all masses _ can be well fitted by an appropriate scaling of a `` universal '' profile _ with no free shape parameters _ ' ( emphasis added ) and the same paper presents an analysis of nfw rotation curves comparing them directly to those of dwarf galaxies . one can see how they would have provided an impetus for the cusp / core debate . simulations thus need to improve by another order of magnitude in resolution before the predictions become relevant for the observational regime currently probed . in summary , current simulations do not yet tell us anything about the nature of dark matter at scales @xmath24 kpc . the study of the distribution of dark matter within the visible discs of galaxies is now clearly driven by observations . increases in computing power will at some stage resolve the inner kpc in simulations as well , but given the `` coldness '' of cdm and its lack of cross - section it is unlikely that the slopes will decrease in the _ dark - matter - only _ simulations ( unless of course the properties of the dark matter particle are changed ) . it is thus interesting to note that the problems of cdm only become clear on length scales where the baryons start playing a role ( this applies to both the cusp / core problem as well as the missing dwarfs problem ) . this indicates that baryon physics is one of the missing pieces of the puzzle , and will very likely make a major contribution toward a solution ( see , for example , sancisi , this volume ) .    can observations of galaxy dynamics at scales larger than a kpc , where the simulations do have predictive power and baryons may be unimportant , help resolve the issue ? the answer to this is yes and no . cdm makes a robust prediction that the density of matter should drop off as @xmath25 at large radii in galaxies . however , finding such a drop - off is not so much a triumph for cdm , as it is a triumph for gravity . the @xmath25 behavior of dark matter at large scales is simply the behaviour of _ any _ pressure - less medium under the influence of gravity . the success of large - scale structure simulations has shown that at large scales this is exactly how dark matter behaves . an @xmath25 drop - off does _ not _ tell us how dark matter behaves at small scales and large densities . maybe the observed dark matter cores are telling us that dark matter has pressure at small scales . maybe dark matter is not cold at these scales . maybe they tell us something unexpected . the centers of galaxies are the only places where we can study dark matter under these conditions . it would thus be unwise to ignore the conclusions forced upon us by the data . * the available data overwhelmingly prefer soft cores . * there is a large range in inner slopes . observationally there is no evidence for a `` universal profile '' . * systematic effects play no significant rle . * observations probe a regime where simulations have no predictive power ( yet ) . * crucial ( baryon ) physics are still missing from the simulations . de blok , w.  j.  g. , bosma , a. , mcgaugh , s.  2003 , , 340 , 657 de blok , w.  j.  g.   bosma , a.  2002 , , 385 , 816 de blok , w.  j.  g. , mcgaugh , s.  s. , rubin , v.  c.  2001a , , 122 , 2396 de blok , w.  j.  g. , mcgaugh , s.  s. , bosma , a. , rubin , v.  c.  2001b , , 552 , l23 de blok , w.  j.  g.   mcgaugh , s.  s.  1997 , , 290 , 533 hayashi , w. , navarro , j.f . , power , c. , jenkins , a.r . , frenk , c.s . , white , s.d.m . , springel , v. , stadel , j. , quinn , t.r . , 2003 , mnras , submitted ( astro - ph/0310576 ) marchesini , d. , donghia , e. , ch incarini , g. , firmani , c. , conconi , p. , molinari , e. , zacchei , a. , 2002 , , 5 75 , 801 mcgaugh , s.  s. , rubin , v.  c. , de blok , w.  j.  g.  2001 , , 122 , 2381 navarro , j.f . white , s.d.m . , 1996 , , 462 , 563 navarro , j.f . , frenk , c.s . , white , s.d.m . , 1997 , , 490 , 493 pickering , t.  e. , navarro , j.  f. , rix , h .- w . , impey , c.  d.  1998 , asp conf . ser .  136 : galactic halos , 199 swaters , r.  a. , madore , b.  f. , van den bosch , f.  c. , balcells , m.  2003 , , 583 , 732 swaters , r.  a. , madore , b.  f. , trewhella , m.  2000 , , 531 , l107
Abstract:

Answer: low surface brightness ( lsb ) galaxies are dominated by dark matter . high - resolution rotation curves suggest that their total mass - density distributions are dominated by constant density cores rather than the steep and cuspy distributions found in cold dark matter ( cdm ) simulations . the data are best described by a model with a soft core with an inner power - law mass - density slope @xmath0 . however no single universal halo profile provides an adequate description of the data . the observed mass profiles appear to be inconsistent with @xmath1cdm . # 1_#1 _ # 1_#1 _ =    # 1 1.25 in .125 in .25 in


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 be able to fully describe hadrons from first principles in quantum chromodynamics ( qcd ) , one needs to be able to compute wave functions and spectra . progress toward making such a computation in the context of light - cone quantization has been steady  @xcite . two particular approaches will be discussed here . both make use of the discretized light - cone quantization ( dlcq ) technique pioneered by pauli and brodsky  @xcite for the numerical solution of quantum field theories . the approaches differ in how the theory is regulated ; one uses pauli  villars ( pv ) regularization  @xcite and the other supersymmetry . the use of pv regularization with dlcq in ( 3 + 1)-dimensional theories has been developed by brodsky , hiller , and mccartor  @xcite . the initial work involved simple many - body models  @xcite . new work discussed here and elsewhere  @xcite is for yukawa theory in a single - fermion truncation . the essential idea is to include pv particles in the dlcq basis . cancellation of ultraviolet infinities is then arranged by choosing imaginary couplings or an indefinite metric . in the case of yukawa theory there are choices to be made about the particle content of the pv sector . the work to be discussed here used three heavy scalars , two of which have negative norm . an alternative now under investigation is to use one heavy scalar and one heavy fermion , both with negative norm . this alternative was suggested by the work of paston _ et al_.  @xcite and has the advantage of being free of instantaneous fermion interactions . supersymmetric theories are simpler with respect to regularization but require greater care in the numerical discretization . a key insight to the correct discretization was made by matsumura , sakai , and sakai  @xcite . they noticed that by discretizing the supercharge @xmath0 and constructing a discrete light - cone energy @xmath1 from the superalgebra , one could retain supersymmetry in the discrete approximation , now called sdlcq  @xcite . ordinary dlcq discretizes @xmath1 directly and recovers supersymmetry only in the infinite resolution limit . the sdlcq technique has been refined and applied by pinsky and collaborators  @xcite , initially in 1 + 1 dimensions , but now also in 2 + 1 dimensions  @xcite . recent work on ( 2 + 1)-dimensional supersymmetric yang  mills theory is discussed here and elsewhere  @xcite . both approaches rely on dlcq . all light - cone momentum variables are discretized , with @xmath2 and @xmath3 , in terms of longitudinal and transverse length scales @xmath4 and @xmath5 . the integrals over wave functions that make up the mass eigenvalue problem @xmath6 are then approximated by the trapezoidal quadrature rule . this yields a matrix eigenvalue problem which is typically quite large but also quite sparse . lanczos techniques  @xcite are used to extract eigenvalues and eigenvectors for the lowest states , even in the case of an indefinite metric  @xcite .    because the longitudinal momentum is always positive , there exists a positive integer @xmath7 , called the ( harmonic ) resolution  @xcite , such that the total longitudinal momentum is @xmath8 and momentum fractions are given by @xmath9 . wave functions and the mass eigenvalue problem , where @xmath10 , are naturally expressed in terms of momentum fractions and the resolution @xmath7 . hence @xmath4 disappears , and @xmath7 effectively takes its place as the resolution scale . the transverse scale @xmath5 is set by a momentum cutoff and a transverse resolution . the remainder of this paper contains brief discussions of recent work on yukawa theory  @xcite in sec . [ sec : yukawa ] and super yang  mills theory  @xcite in sec . [ sec : sym ] . a summary and discussion of future work are given in sec . [ sec : future ] . when terms involving antifermions are eliminated , the yukawa light - cone hamiltonian becomes  @xcite @xmath11       b_{\underline{n}+\underline{m},-s}^\dagger b_{\underline{n},s }                                               a_{\underline{m } }   \nonumber \\    & & \rule{0.5in}{0 mm } + \mbox{h.c . } \nonumber \\     & & + \frac{mg}{\sqrt{8\pi}l_\perp }               \sum_{\underline{n}\underline{m}s}\frac{1}{\sqrt{m } }       \left[\frac{1}{n / k}\right .   \\ & &   + \left.\frac{1}{(n+m)/k}\right ]       b_{\underline{n}+\underline{m},s}^\dagger b_{\underline{n},s }                   a_{\underline{m } }     + \mbox{h.c . }   \nonumber \\     & & + \frac{g^2}{8\pi l_\perp^2}\sum_{\underline{n}\underline{m}\underline{m } 's }           \frac{1}{\sqrt{mm ' } }   \nonumber \\ & & \times   \left[b_{\underline{n}+\underline{m}+\underline{m}',s}^\dagger            b_{\underline{n},s } a_{\underline{m } ' } a_{\underline{m}}\frac{1}{(n+m)/k }                + \mbox{h.c . }   \right . \nonumber \\      & & + b_{\underline{n}+\underline{m}-\underline{m}',s}^\dagger              b_{\underline{n},s } a_{\underline{m}'}^\dagger a_{\underline{m } }   \nonumber \\    & &     \left.\times\left(\frac{1}{(n - m')/k}+\frac{1}{(n+m)/k}\right)\right ] \nonumber \end{aligned}\ ] ] where @xmath12            & = & \delta_{\underline{m},\underline{m}'}\ , , \nonumber\end{aligned}\ ] ] and @xmath13 the one - loop contribution to the fermion self energy is  @xcite @xmath14\ , , \nonumber \end{aligned}\ ] ] with @xmath15 a cutoff . for this to be finite and consistent with zero in the @xmath16 chiral limit , we need to perform three subtractions . therefore we add to the hamiltonian three heavy pv scalars  @xcite with masses @xmath17 and couplings @xmath18 determined by @xmath19 these constraints guarantee the desired one - loop subtractions if the norm of the i - th pv boson is @xmath20 . fermion self - induced inertia terms have not been included because they cancel once the heavy scalars are added . the state vector is @xmath21 with @xmath22 and normalization @xmath23    to cancel an infrared singularity in the instantaneous fermion term , we add an effective interaction that represents the contribution of the missing z graph . the added term is built diagrammatically from the pair production and annihilation terms @xmath24       b_{\underline{p},s}^\dagger d_{\underline{q}-\underline{p},s}^\dagger   \xi_i a_{i\underline{q } }       + \mbox{h.c . } \nonumber \\     & & + \frac{mg}{2l_\perp\sqrt{2l } }     \sum_{\underline{p}\underline{q}si }      \left[\frac{1}{p^+\sqrt{q^+ } }   \right . \\ & & -\left.\frac{1}{(q^+-p^+)\sqrt{q^+}}\right ]       b_{\underline{p},s}^\dagger d_{\underline{q}-\underline{p},-s}^\dagger   \xi_i a_{i\underline{q } }       + \mbox{h.c.}\ , ,   \nonumber\end{aligned}\ ] ] and the denominator for the intermediate state @xmath25    the bare parameters of the hamiltonian , @xmath26 and @xmath27 , are determined by input of `` data . '' the mass @xmath28 of the dressed single - fermion state is held fixed . this is imposed by rearranging the mass eigenvalue problem into an eigenvalue problem for @xmath27 : @xmath29 \tilde{\phi } } \\ & & -\int\prod_j dy'_j d^2q'_{\perp j}\sqrt{xx'}{\cal k}\tilde{\phi } ' = \delta m^2\tilde{\phi}\ , ,   \nonumber\end{aligned}\ ] ] where @xmath30 represents the original kernel and amplitudes are related by @xmath31 . to fix the coupling we use @xmath32 . from a numerical solution it can be computed fairly efficiently in a sum similar to the normalization sum @xmath33 the constraint on @xmath34 is satisfied by solving it simultaneously with the eigenvalue problem .    with the parameters fixed , we can compute various quantities . those considered include structure functions , the form factor slope at zero momentum transfer , average numbers of constituents , and average constituent momenta . a representative plot of the bosonic structure function @xmath35 is given in fig . [ fig : structfn ] . additional results can be found in ref . @xcite . for good results at stronger couplings , where first - order perturbation theory is insufficient , rather high resolution was required , with @xmath36 to 39 and as many as 15 transverse momentum points . this resolution was achieved by limiting the number of constituents to 3 , after verifying that the contribution from higher sectors was sufficiently small . the lagrangian for supersymmetric su(n ) yang  mills theory in 2 + 1 dimensions is @xmath37 where @xmath38 $ ] and @xmath39 $ ] . we work in light - cone gauge ( @xmath40 ) and the large-@xmath41 limit . the dynamical fields are @xmath42 and @xmath43 . the supercharge @xmath0 is @xmath44              + 2\psi\psi\right)\frac{1}{\partial_-}\psi\right]\,.\end{aligned}\ ] ] and @xmath1 is given by @xmath45 . by discretizing @xmath0 and computing @xmath1 from this anticommutator , supersymmetry is exactly preserved in the numerical approximation  @xcite . in addition to supersymmetry , we have transverse parity and the kutasov @xmath46 symmetry  @xcite . the mass eigenvalue problem can then be solved separately in each of the 8 symmetry sectors . in the largest calculation to date  @xcite , each sector contained roughly 230,000 basis states . a variety of results can be seen in ref . they include a number of studies of the coupling dependence of the mass values , as well as outcomes for structure functions . one striking feature is that , except for weak coupling , the average particle count for almost all eigenstates is at or near the maximum allowed by the resolution , even for the highest resolution considered . this means that for strong coupling the method does not capture all of the important pieces of the true eigenstate . such behavior is likely to be a consequence of dealing with massless constituents without introducing a mass scale through symmetry breaking or other means . another striking feature is that the average number of fermions in bosonic states has been observed to have a gap between 4 and 6 . this could only be checked for resolutions where a full diagonalization could be done . however , as can be seen in fig . [ fig:<nf > ] , the signal is quite clean for @xmath47 and 3 transverse momentum modes , where the average number is never above 4 unless it is precisely 6 , the maximum allowed ; below 4 there are no gaps at all . both ordinary dlcq with pauli  villars particles and sdlcq provide the means to compute masses and wave functions for eigenstates in multi - dimensional quantum field theories . these methods will continue to be explored in various contexts , leading eventually to consideration of ( supersymmetric ) qcd . in fact , paston _ et al_.  @xcite have already obtained a pv - like regularization of qcd that could , in principle , be solved by dlcq ; however , with present computing power the number of fields is probably too large for meaningful calculations . with qcd as a goal , work on the pv approach will next turn to the alternative regularization of yukawa theory , with one heavy scalar and one heavy fermion  @xcite , in both the one and two - fermion fock sectors . quantum electrodynamics will also be considered , as a first application to a gauge theory and as something of interest in its own right . the next step for work with sdlcq is to include a chern  simons term in super yang  mills theory and thereby give each constituent a nonzero mass . the extension of the method to 3 + 1 dimensions is also important , as is consideration of supersymmetry breaking . the work reported here was done in collaboration with s.j . brodsky and g. mccartor and with s.s . pinsky and u. trittmann and was supported in part by the department of energy , contract de - fg02 - 98er41087 , and by grants of computing time from the minnesota supercomputing institute . o. lunin and s. pinsky , in the proceedings of 11th international light - cone school and workshop : new directions in quantum chromodynamics and 12th nuclear physics summer school and symposium ( nuss 99 ) , seoul , korea , 26 may - 26 jun 1999 ( new york , aip , 1999 ) , p.  140 , hep - th/9910222 . f.  antonuccio , o.  lunin , and s.  pinsky , phys . d 59 ( 1999 ) 085001 , hep - th/9811083 ; p. haney , j. r. hiller , o. lunin , s. pinsky , and u. trittmann phys . d 62 ( 2000 ) 075002 , hep - th/9911243 ; j.r . hiller , s.s . pinsky , and u. trittmann , phys . d 63 ( 2001 ) 105017 , hep - th/0101120 .        c. lanczos , j. res . bur .  stand . * 45 * , 255 ( 1950 ) ; j.h . wilkinson , _ the algebraic eigenvalue problem _ ( clarendon , oxford , 1965 ) ; b.n . parlett , _ the symmetric eigenvalue problem _ ( prentice  hall , englewood cliffs , nj , 1980 ) ; j. cullum and r.a . willoughby , j. comput . * 44 * , 329 ( 1981 ) ; _ lanczos algorithms for large symmetric eigenvalue computations _ ( birkhauser , boston , 1985 ) , vol . i and ii ; g.h . golub and c.f . van loan , _ matrix computations _ ( johns hopkins university press , baltimore , 1983 ) .
Abstract:

Answer: applications of discretized light - cone quantization ( dlcq ) to ( 3 + 1)-dimensional yukawa theory with pauli  villars regulators and of supersymmetric dlcq ( sdlcq ) to ( 2 + 1)-dimensional super yang mills theory are discussed . the ability of these methods to provide wave functions as well as spectra is emphasized .


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: energy harvesting facilitates the battery charging to prolong the lifetime of wireless networks , especially under some extreme conditions , such as battle - field , underwater and body area networks @xcite . wherein , electromagnetic wave based wireless power transfer has received considerable research attentions from academic and industry due to two - fold reasons . first , it is a controllable power transfer mode . second , information and power can be simultaneously transferred in the form of electromagnetic wave @xcite @xcite . similar to information transmission , wireless power transfer may suffer from channel fading , resulting in low energy efficiency . in this end , energy beamforming is proposed by deploying multiple antennas at the power source @xcite @xcite . the impact of channel state information ( csi ) at the power source on the performance of wireless information and power transfer ( wipt ) is quantitatively analyzed in @xcite . recently , large scale multiple input multiple output ( ls - mimo ) techniques are also introduced to significantly improve the efficiency of wipt by exploiting the high spatial resolution @xcite . moreover , relaying technique is also proved as an effective way of improving the performance of wipt by shortening the transfer distance , since the distance has a great impact on both information and power transfer @xcite . a two - way relaying scheme is proposed in @xcite to offer a higher transmission rate with the harvested energy . in fact , through combing relaying schemes and multi - antenna techniques , especially the ls - mimo techniques , the performance of wipt can be improved significantly , even in the case of long - distance transfer . however , to the best of our knowledge , there is no work studying the problem of the ls - mimo relaying techniques for wipt . meanwhile , information transfer may encounter interception from the eavesdropper due to the broadcast nature of wireless channels . in recent years , as a supplementary of encryption techniques , physical layer security is used to realize secure communications , by exploiting wireless channel characteristics , i.e. , fading and noise . the performance of physical layer security is determined by the performance difference between the legitimate and eavesdropper channels @xcite . thus , the ls - mimo relaying technique is also an effective way of improving the secrecy performance @xcite . for secure wipt , long - distance transfer and short - distance interception are two challenging issues @xcite @xcite . to solve them , we introduce the amplify - and - forward ( af ) ls - mimo relaying technique into secure wipt . the contributions of this paper are two - fold :    1 . we derive an explicit expression of the secrecy outage capacity for such a secure wipt system in terms of transmit power , csi accuracy and transfer distance . we propose a power splitting scheme at the relay , so as to maximize the secrecy outage capacity . the rest of this paper is organized as follows . we first give an overview of the secure wipt system based on the ls - mimo af relaying scheme in section ii , and then derive the secrecy outage capacity under imperfect csi in section iii . in section iv , we present some simulation results to validate the effectiveness of the proposed scheme . finally , we conclude the whole paper in section v. we consider a time division duplex ( tdd ) ls - mimo af relaying system , where a single antenna source communicates with a single antenna destination aided by a multi - antenna relay , while a single antenna passive eavesdropper intends to intercept the message , as shown in fig.[fig1 ] . note that the number of antenna at the relay @xmath0 is quite large for such an ls - mimo relaying system , i.e. , @xmath1 or greater . the relay only has limited energy to maintain the active state , so it needs to harvest enough energy from the source for information transmission . considering the limited storage , the relay should be charged from the source slot by slot . the whole system is operated in slotted time of length @xmath2 , and the relay works in the half - duplex mode . then , the information transmission from the source to the destination via the aid of the relay requires two time slots . specifically , in the first time slot , the source sends the signal , and the relay splits the received signal into two streams , one for energy harvesting and the other for information processing . during the second time slot , the relay forwards the post - processing signal to the destination with the harvested energy . note that the direct link from the source to the destination is unavailable due to a long distance . we assume that the eavesdropper is far away from the source and is close to the relay , since it thought the signal comes from the relay . note that it is a common assumption in related literature @xcite @xcite , since it is difficult for the eavesdropper to overhear the signals from the source and the relay simultaneously . then , the eavesdropper only monitors the transmission from the relay to the destination . we use @xmath3 to denote the channel from @xmath4 to @xmath5 , where @xmath6 and @xmath7 with @xmath8 representing the source , the relay , the destination and the eavesdropper , respectively . @xmath9 is the distance - dependent path loss and @xmath10 is the small scale fading . in this paper , we model @xmath10 as gaussian distribution with zero mean and unit variance . @xmath9 remains constant during a relatively long period and @xmath10 fades independently slot by slot . thus , the received signal at the relay in the first time slot can be expressed as @xmath11 where @xmath12 is the normalized gaussian distributed transmit signal , @xmath13 is the transmit power at the source , and @xmath14 is the additive gaussian white noise with zero mean and unit variance at the relay . as mentioned above , the relay splits the received signal @xmath15 into energy harvesting stream with proportional factor @xmath16 and signal processing stream with @xmath17 . then , according to the law of energy conservation , the harvesting energy at the relay is given by @xmath18 where the constant parameter @xmath19 , scaling from 0 to 1 , is the efficiency ratio at the relay for converting the harvested energy to the electrical energy to be stored @xcite @xcite . with the harvested energy , the relay forwards the post - processing signal @xmath20 to the destination . then the received signals at the destination and the eavesdropper are given by @xmath21 and @xmath22 respectively , where @xmath23 and @xmath24 are the additive gaussian white noises with zero mean and unit variance at the destination and the eavesdropper . @xmath25 is the post - processing signal with @xmath26 being a transform matrix . we assume the relay has full csi @xmath27 by channel estimation , and gets partial csi @xmath28 via channel reciprocity in tdd systems . due to duplex delay between uplink and downlink , there is a certain degree of mismatch between the estimated csi @xmath29 and the real csi @xmath28 , whose relation can be expressed as @xcite @xmath30 where @xmath31 is the error noise vector with independent and identically distributed ( i.i.d . ) zero mean and unit variance complex gaussian entries . @xmath32 , scaling from 0 to 1 , is the correlation coefficient between @xmath29 and @xmath28 . a larger @xmath32 means better csi accuracy . if @xmath33 , the relay has full csi @xmath28 . additionally , due to the hidden property of the eavesdropper , the csi @xmath34 is unavailable . therefore , @xmath26 is designed only based on @xmath27 and @xmath29 , but is independent of @xmath34 . considering the better performance of maximum ratio combination ( mrc ) and maximum ratio transmission ( mrt ) in ls - mimo systems , we design @xmath26 by combining mrc and mrt . mathematically , the transform matrix is given by @xmath35 where @xmath36 is the power constraint factor . to fulfill the energy constraint at the relay , @xmath36 can be computed as @xmath37    based on the af relaying scheme , the signal - to - noise ratio ( snr ) at the destination is given by ( [ eqn9 ] ) at the top of the next page ,    @xmath38    where @xmath39 , @xmath40 and @xmath41 . similarly , the received snr at the eavesdropper is given by ( [ eqn10 ] ) at the top of this page ,    @xmath42    where @xmath43 and @xmath44 . letting @xmath45 and @xmath46 be the legitimate channel and the eavesdropper channel capacities , then the secrecy rate is given by @xmath47^+$ ] , where @xmath48^+=\max(x,0)$ ] @xcite . since there is no knowledge of the eavesdropper channel at the source and relay , it is impossible to maintain a steady secrecy rate over all realizations of fading channels . in this context , we take the secrecy outage capacity @xmath49 as the performance metric , which is defined as the maximum available rate while the outage probability that the transmission rate surpasses the secrecy rate is equal to a given value @xmath50 , namely @xmath51 in this section , we first analyze the secrecy outage capacity of such an ls - mimo af relaying system with energy harvesting , and then derive an optimal power splitting scheme to determine the proportional factor @xmath16 .      according to ( [ eqn11 ] ) , the secrecy outage capacity is jointly determined by the legitimate and eavesdropper channel capacities , so we first analyze the legitimate channel capacity . based on the received snr at the destination in ( [ eqn9 ] ) , we have the following theorem :    _ theorem 1 _ : the legitimate channel capacity in an ls - mimo af relaying system under imperfect csi can be approximated as @xmath52 , where @xmath53 is a half of the spectral bandwidth . please refer to appendix i.    it is found that the legitimate channel capacity is a constant due to channel hardening in such an ls - mimo af relaying system . then , according to the definition of the secrecy outage capacity , we have a theorem as below :    _ theorem 2 _ : given an outage probability bound by @xmath50 , the secrecy outage capacity for an ls - mimo af relaying scheme is @xmath54 . please refer to appendix ii .    while theorem 2 is useful to study the secure wireless information and power transfer system s secrecy outage capacity , the expression is in general too complex to gain insight . motivated by this , we carry out asymptotic analysis on the secrecy outage capacity at high transmit power regime , and derive the following theorem :    _ theorem 3 _ : at high transmit power regime , the secrecy outage capacity in this case is independent of transmit power @xmath13 . there exists a performance upper bound on the secrecy outage capacity .    please refer to appendix iii . _ remarks _ : the secrecy outage capacity will be saturated once the transmit power surpasses a threshold . this is because the af relaying system is noise limited due to noise amplification at the relay in this case . thus , it makes sense to find the minimum power that achieves the maximum secrecy outage capacity . according to theorem 2 , for a given transmit power , the secrecy outage capacity is a function of power splitting ratio @xmath16 . intuitively , a large @xmath16 leads to a high transmit power at the relay , but the received signal power decreases . moreover , from the view of wireless security , a high transmit power at the relay may also increases the interception probability . thus , it is necessary to optimize the power splitting ratio , so as to maximize the secrecy outage capacity . since @xmath55 is an increasing function , in order to maximize the secrecy outage capacity , it is equivalent to maximizing the term @xmath56 . thus , the optimal power splitting can be described as the following optimization problem : @xmath57    the objective function is not concave , so it is difficult to provide a closed - form expression for the optimal @xmath16 . however , because ( [ eqn13 ] ) is a one - dimensional function of @xmath16 , it is possible to get the optimal @xmath16 by numerical searching . specifically , by scaling @xmath16 from 0 to 1 with a fixed step , the optimal @xmath16 related to the maximum objective function is obtained . to examine the accuracy and effectiveness of the derived theoretical results for secure wipt in an ls - mimo af relaying system , we present several simulation results in the following scenarios : we set @xmath1 , @xmath58khz , @xmath59 , @xmath60 and @xmath61 without extra explanation . the relay is in the middle of a line between the source and the destination . we normalize the path loss as @xmath62 , and use @xmath63 to denote the relative interception path loss . specifically , if @xmath64 , the interception distance is shorter than the legitimate propagation distance , and then the interception becomes strong . in addition , we use snr@xmath65 to represent the transmit signal - to - noise ratio ( snr ) in db at the source . .,scaledwidth=50.0% ]    first , we validate the accuracy of the theoretical results with snr=10db . as shown in fig.[fig3 ] , the theoretical results coincide with the simulations nicely in the whole @xmath63 region under different requirements of outage probability . given a outage probability bound by @xmath50 , as @xmath63 increases , the secrecy outage capacity decreases gradually . this is because the interception capacity of the eavesdropper enhances due to the shorter interception distance . on the other hand , for a given @xmath63 , the secrecy outage capacity improves with the increase of @xmath50 , since the secrecy outage capacity is an increasing function of outage probability . , scaledwidth=50.0% ]    second , we investigate the impact of power transfer distance on the secrecy outage capacity with @xmath66 , @xmath67 and snr=0db for secure swip . note that optimal power splitting is adopted to improve the secrecy outage capacity . as seen in fig.[fig2 ] , even at a small @xmath68 , namely long transfer distance , there is a large secrecy outage capacity . as a simple example , at @xmath69 , the proposed scheme with @xmath70 can achieve @xmath71 kb / s . as @xmath72 increases , the secrecy outage capacity significantly improves . thus , the proposed scheme can solve the challenging problem of long - distance transfer for secure wipt . , scaledwidth=50.0% ]    then , we examine the effect of interception distance on the secrecy outage capacity with @xmath67 and snr=0db for secure swip . similarly , optimal power splitting is adopted to improve the secrecy outage capacity . as seen in fig.[fig6 ] , even at a large @xmath63 , namely short interception distance , there is a large secrecy outage capacity . as a simple example , at @xmath73 , the proposed scheme with @xmath70 can achieve @xmath74 kb / s . as @xmath72 increases , the secrecy outage capacity significantly improves . thus , the proposed scheme can solve the challenging problem of short - distance interception for secure wipt . next , we show the impact of snr on the secrecy outage capacity with @xmath66 . as seen in fig.[fig4 ] , the secrecy outage capacity is an increasing function of snr . however , as snr increases , the secrecy outage capacity will be saturated . this is because the af relaying system is noise limited due to noise amplification , which confirmed the claim in theorem 3 . thus , it makes sense to find the optimal snr to maximize the secrecy outage capacity in ls - mimo relaying systems with the minimum transmit power . finally , we check the effectiveness of the proposed optimal power splitting scheme with respect to a fixed scheme . specifically , the fixed scheme sets @xmath16 as 0.1 fixedly regardless of channel condition and transmit power . as seen in fig.[fig5 ] , the proposed scheme performs better obviously . especially , as snr increases , the performance gain becomes larger . note that the proposed scheme also will suffer from performance saturation . however , the performance bound can be lifted by adding antennas at the relay , which is a main advantage of the ls - mimo relaying scheme for secure wipt . a major contribution of this paper is the introduction of the ls - mimo relaying technique into secure wireless information and power transfer to significantly enhance wireless security and improve transmission rate . this paper derives a closed - form expression of the secrecy outage capacity in terms of transmit snr , power splitting ratio , antenna number and interception distance . furthermore , through maximizing the secrecy outage capacity , we present a power splitting scheme , which has a huge performance gain with respect to the fixed scheme . based on the snr @xmath75 at the destination , the legitimated channel capacity can be expressed as @xmath76 where @xmath53 is a half of the spectral bandwidth , since a complete transmission requires two time slots . @xmath77 denotes the real part of @xmath78 . @xmath28 is replaced by @xmath79 in ( [ app1 ] ) . ( [ app2 ] ) follows from the fact that @xmath80 scales with the order @xmath81 as @xmath82 while @xmath83 scales as the order @xmath84 , which can be negligible . ( [ app3 ] ) holds true because of @xmath85 and @xmath86 , namely channel hardening @xcite . therefore , we get the theorem 1 . according to ( [ eqn11 ] ) , for a given @xmath50 , we have @xmath87 where @xmath88 is the cumulative distribution function ( cdf ) of @xmath89 . in order to derive the secrecy outage capacity , the key is to get the cdf of @xmath89 . examining ( [ eqn10 ] ) , due to channel hardening , we have @xmath90 since @xmath91 is an isotropic unit vector and independent of @xmath34 , @xmath92 is @xmath93 distributed with 2 degrees of freedom . let @xmath94 , we can derive the cdf of @xmath89 as @xmath95 if @xmath96 , then we have @xmath97 since @xmath98 is impossible when @xmath99 , we have @xmath100 equivalently , we have @xmath101 hence , we get the theorem 2 . for the secrecy outage capacity in ( [ app9 ] ) , if the transmit power @xmath13 is sufficiently large , it can be approximated as @xmath102 where ( [ app10 ] ) holds true since the constant  1 " is negligible when @xmath13 is large enough . it is found that the secrecy outage capacity in this case is independent of transmit power @xmath13 . thus , we prove the theorem 3 .    1 f. zhang , s. a. hackworth , x. liu , h. chen , r. j. sclabassi , and m. sun ,  wireless energy transfer platform for medical sensor and implantable devices , " in _ proc . ieee embs 31st annual int . _ , pp . 1045 - 1048 , sep . 2009 . x. chen , c. yuen , and z. zhang ,  wireless energy and information transfer tradeoff for limited feedback multi - antenna systems with energy beamforming , " _ ieee trans . veh . 407 - 412 , jan . 2014 . x. chen , x. wang , and x. chen ,  energy - efficient optimization for wireless information and power transfer in large - scale mimo systems employing energy beamforming , " _ ieee wireless commun . 2 , no . 6 , pp . 667 - 670 , dec . 2013 . d. s. michalopoulos , h. a. suraweera , and r. schober ,  relay selection for simultaneous information transfer and wireless energy transfer : a tradeoff perspective , " [ online ] : http://arxiv.org/abs/1303.1647v1 . b. m. hochwald , t. l. marzetta , and v. tarokh ,  multiple - antenna channel hardening and its implications for rate - feedback and scheduling , " _ ieee tran . inf . theory _ , 1893 - 1909 , sept . 2004
Abstract:

Answer: in this paper , we address the problem of secure wireless information and power transfer in a large - scale multiple - input multiple - output ( ls - mimo ) amplify - and - forward ( af ) relaying system . the advantage of ls - mimo relay is exploited to enhance wireless security , transmission rate and energy efficiency . in particular , the challenging issues incurred by short interception distance and long transfer distance are well addressed simultaneously . under very practical assumptions , i.e. , no eavesdropper s channel state information ( csi ) and imperfect legitimate channel csi , this paper investigates the impact of imperfect csi , and obtains an explicit expression of the secrecy outage capacity in terms of transmit power and channel condition . then , we propose an optimal power splitting scheme at the relay to maximize the secrecy outage capacity . finally , our theoretical claims are validated by simulation results .