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 lepton flavor violating ( lfv ) neutrinoless conversion of a bound 1s - muon to electron in the field of a nucleus @xcite-@xcite @xmath10 is known as one of the best probes to search for the hypothetical muon and electron flavor non - conservation @xcite . so far , the experiments @xcite-@xcite seeking for @xmath11 conversion events have only succeeded to put upper bounds on the branching ratio @xmath12 i.e. the ratio of the muon - electron conversion rate relative to the total rate of the ordinary muon capture . the best upper limits have been extracted at psi by the sindrum ii experiments in the values ( at 90% confidence level ) @xmath13 @xmath14 these limits are improvements over the previous limits set at triumf @xcite a decade ago using the same targets .    at present , two @xmath1 conversion experiments are launched , the ongoing experiment at psi using @xmath3ti target @xcite , and the planned meco experiment at brookhaven @xcite using @xmath5al target . the expected sensitivity on @xmath2 in the psi experiment is @xmath15 while in the brookhaven experiment it will be roughly @xmath16 which implies an improvement over the existing limits by about four orders of magnitude . the meco experiment is going to be conducted in a new @xmath17 beam line at the ags , where the muons are produced using a pulsed proton beam @xcite . the proton energy will be chosen in the range of 8 - 20 gev to optimize the @xmath18 flux per unit time . furthermore , the number of electrons with energy @xmath19 mev , equal to the energy for the coherent peak in @xmath5al , is very much suppressed . this is in contrast with the case of @xmath20 , where the electron flux from @xmath21 decay is peaked at the energy of the electrons from @xmath22 .    for such experiments the knowledge of nuclear transition matrix elements for all accessible @xmath23 channels of the targets employed are of significant importance @xcite-@xcite . in this work we use the transition matrix elements calculated for the aforementioned isotopes in the coherent mode to constrain the lepton flavor violating parameters of various lagrangians predicting this exotic process ( e.g.  scalar and vector couplings , neutrino mixing angles and masses @xcite-@xcite , supersymmetric r - parity violating couplings @xcite-@xcite etc . ) . to this aim the recent experimental data and the expected sensitivity of the meco experiment . as is well known , only the coherent rate can be measured because it is free from background events from bound muon decay and radiative muon capture followed by a fully asymmetric @xmath24 pair creation @xcite . on the other hand , previous studies of @xmath1 conversion rates @xcite-@xcite have shown that for all mechanisms the coherent mode dominates the process ( [ reaction ] ) which means that this is the most important channel . the incoherent reaction leading to excited nuclear states is suppressed due to pauli blocking effects and it is much harder to calculate but its knowledge is also useful in order to determine the experimentally interesting quantity of the ratio of the coherent to the total @xmath6 rate ( see refs . @xcite-@xcite ) . the family quantum numbers @xmath25 , @xmath26 , @xmath27 are conserved within the standard model ( sm ) in all orders of perturbation theory . however this is an accidental consequence of the sm field content and gauge invariance . physics beyond the sm can easily spoil this property . processes like @xmath11 conversion , which is forbidden in the standard model by muon and electron quantum number conservation , play an important role in the study of flavor changing neutral currents and possible physics beyond the sm .    on the particle physics side @xcite , there are many mechanisms of the @xmath11 conversion constructed in the literature ( see @xcite-@xcite and references therein ) . all these mechanisms fall into two different categories : photonic and non - photonic . mechanisms from different classes significantly differ from the point of view of the nucleon and nuclear structure calculations . this stems from the fact that they proceed at different distances and , therefore , involve different details of the structure . the long - distance photonic mechanisms are mediated by the photon exchange between the quark and the @xmath11-lepton currents . these mechanisms resort to the lepton - flavor non - diagonal electromagnetic vertex which is presumably induced by the non - standard model physics at the loop - level . the @xmath11 lepton - flavor violating loop can appear as the @xmath28 loop [ fig.1(a ) ] with the massive neutrinos @xmath29 and the loop with the supersymmetric particles such as the neutralino(chargino)-slepton(sneutrino ) [ fig.1(d ) ] . in the r - parity violating susy models there are also lepton - slepton and quark - squark loops generated by the superpotential couplings @xmath30 and @xmath31 respectively . the short - distance non - photonic mechanisms contain heavy particles in intermediate states and can be realized at tree level [ fig.2 ] , at 1-loop - level [ fig.1(a , b , d , e ) ] or at the level of box diagrams [ fig.1(c ) ] . the tree - level diagrams can be constructed in the r - parity violating susy models with the virtual z - boson , squarks @xmath32 and sneutrinos @xmath33 [ fig.2 ] . the 1-loop diagrams of the non - photonic mechanisms include the diagrams similar to those for the photonic mechanisms but with the z - boson instead of the photon [ fig.1(a , d ) ] as well as additional z - boson couplings to the neutrinos and neutralinos [ fig.1(b , e ) ] . the box diagrams are constructed of the w - bosons and massive neutrinos [ fig.1(c ) ] as well as similar boxes with neutralinos and sleptons or charginos and sneutrinos . our purpose is to calculate the contribution of the above - described mechanisms to the @xmath6 conversion branching ratio ( [ def1 ] ) . at the first step we construct the effective @xmath6 conversion hamiltonian for these mechanisms in terms of nucleon degrees of freedom of a nucleus involved in the process this will allow us to accomplish the calculation of @xmath2 by applying the conventional nuclear structure methods using non - relativistic impulse approximation . for the photonic diagrams of fig . 1(a , d ) the hadronic vertex represents the usual nucleon electromagnetic current @xmath34 @xmath35 is the nucleon isospin doublet @xmath36 with @xmath37 , @xmath38 being the proton and neutron spinors .    in the case of the non - photonic mechanisms of fig . 1 the hadronic currents are either the sm neutral currents coupled to z - boson [ fig . 1(a , b , d , e ) ] or the effective currents derived after an appropriate fierz transformation from the effective operators of the box diagrams [ fig . 1(c ) ] involving heavy particles . these currents can be written in terms of the nucleon doublet field as @xmath39 n \qquad ( non - photonic ) \label{3.1.4}\ ] ] where @xmath40 , @xmath41 represent the vector , axial vector static nucleon form factors ( @xmath42 ) and the parameters @xmath43 , @xmath44 , @xmath45 , for the models adopted take the values given in ref . the parameter @xmath43 is defined as @xmath46 , with @xmath47 @xmath48 the isoscalar ( isovector ) couplings at the quark level . we should note that , in general , the parameters @xmath43 , @xmath49 and @xmath50 are functions of @xmath51 . for example , in the case of z - exchange we have @xmath52 . the corresponding leptonic currents for the photonic and non - photonic mechanism of figs . 1 are given in ref . @xcite .    in the minimal susy model with a most general form of r - parity violation ( for review see , for instance , @xcite ) all the possible tree - level diagrams for the @xmath6 conversion process are shown in fig . 2 @xcite,@xcite . to write down the low - energy @xmath11 conversion lagrangian at quark level one starts from the well known r - parity violating superpotential and integrates out heavy intermediate fields . for the diagrams of fig . 2 the corresponding 4-fermion low - energy effective lagrangian at the quark level takes the form @xcite ( first order of perturbation theory ) @xmath53 \nonumber \\ & & + \frac{g_f}{\sqrt{2}}\left [ \bar\eta_{r}^{di}j_{dr(i)}j_{l } + \bar\eta_{l}^{di}j_{dl(i)}j_{r}\right ] . \label{eff - q}\end{aligned}\ ] ] ( @xmath54 runs over generations so that @xmath55 with @xmath56 and @xmath57 ) where the coefficients @xmath58 contain the @xmath59 susy parameters @xcite . the color singlet currents @xmath60 and @xmath61 at the quark level are written as @xmath62 where @xmath63 . the leptonic currents are written as @xmath64    at the next step we rewrite eq . ( [ eff - q ] ) , specified at the quark level , in terms of the nucleon degrees of freedom . this is usually achieved by utilizing the on - mass - shell matching condition @xcite and gives @xmath65 . \label{nucl1}\ ] ] where now the hadronic ( nucleon ) currents are @xmath66 \,n ,   \label{nucur } \\ j^{\pm } & = & \bar{n}\,\,\left [ ( \alpha _ { \pm s}^{(0)}+\alpha _ { \pm s}^{(3)}\tau _ { 3})+(\alpha _ { \pm p}^{(0)}+\alpha _ { \pm p}^{(3)}\tau _ { 3})\gamma _ { 5}\right ] \,n ,   \nonumber\end{aligned}\ ] ] the coefficients @xmath67 , with @xmath68 , in eq . ( [ nucur ] ) include the nucleon form factors ( functions of the momentum transfer @xmath69 ) for scalar , vector , axial vector and pseudoscalar contributions , respectively ( for their definition see ref . @xcite ) . since , however , the maximum momentum transfer @xmath69 in @xmath1 conversion is much smaller than the typical scale of nucleon structure ( @xmath70 , with @xmath71 , the muon mass ) , we can safely neglect the @xmath69-dependence of the nucleon form factors . one of the most interesting quantities in @xmath6 conversion , both from theoretical and experimental points of view , is the branching ratio @xmath2 defined in eq . ( [ def1 ] ) . the expression which gives @xmath2 in the case of the dominant coherent channel has been written as @xcite @xmath72 where @xmath73 is nearly independent of nuclear physics  @xcite and contains the flavor - violating parameters mentioned before . thus , e.g. for photon - exchange @xmath73 is given by @xmath74 this expression contains the electromagnetic form factors @xmath75 , @xmath76 , @xmath77 , @xmath78 parametrized in a specific elementary model  @xcite . the function @xmath79 of eq . ( [ iii.2 ] ) includes about all the nuclear information . by assuming that the total rate of the ordinary muon capture rate is described by the goulard - primakoff function @xmath80 , @xmath79 is defined as @xmath81 where @xmath82 . thus , the nuclear aspects of the @xmath23 branching ratio @xmath2 are mainly included in the matrix elements @xmath83 @xcite , which in the proton - neutron representation , are written as @xmath84    in general , the transition matrix elements @xmath85 in eq . ( [ me2 ] ) depend on the final nuclear state populated during the @xmath11 conversion . for ground state to ground state transitions ( @xmath86 ) in spherically symmetric nuclei the following integral representation is valid @xmath87 where @xmath88 the zero order spherical bessel function and @xmath89 the proton ( p ) , neutron ( n ) nuclear density normalized to the atomic number @xmath90 and neutron number @xmath35 of the participating nucleus , respectively . the space dependent part of the muon wave function @xmath91 , a spherically symmetric function , can be obtained by solving numerically the schrdinger and dirac equations with the coulomb potential .    for the coherent rate in light nuclei the factorization approximation ( see ref . @xcite ) is very good and @xmath92 can be expressed in terms of the nuclear form factors @xmath93 ( for protons ) and @xmath94 ( for neutrons ) which are easily estimated ( for this channel only the ground state wave function of the studied nucleus is required ) . these form factors are defined as @xmath95 @xmath96 and contain the single - particle orbit occupancies @xmath97 for the evaluation of which one must use a nuclear model in the proton - neutron representation , e.g. qrpa ( for @xmath3ti and @xmath4pb see ref . @xcite ) , shell - model ( for @xmath5al see ref . @xcite ) etc .    by using the lagrangian ( [ nucl1 ] ) one can deduce a similar expression to that of eq . ( [ iii.2 ] ) for the @xmath11 conversion branching ratio as @xmath98 where @xmath99 \ } \label{rme.1}\end{aligned}\ ] ] with @xmath100 the quantity @xmath101 in eq . ( [ rme.1 ] ) depends weakly on the nuclear parameters determining the factor @xmath102 . in fact , the terms depending on @xmath102 are small and in practice the nuclear dependence of @xmath101 can be neglected . thus , the corresponding upper bounds on @xmath101 determine the sensitivity to the susy signals arriving at the detector in various @xmath11 conversion experiments ( see results below ) .    a usual approximate expression for the ratio @xmath102 of eq . ( [ phi ] ) , containing the nuclear parameters @xmath103 , the atomic weight , and @xmath90 the total charge of the nucleus , is the following @xmath104 the latter equation can be obtained assuming that @xmath105 and @xmath106 where @xmath107 represents the mean value of the muon wave function @xmath108 . for light and medium nuclei @xmath109 and eq . ( [ phitil ] ) is a good approximation . for isoscalar nuclei , i.e. @xmath110 , @xmath111 ( see sect . the pure nuclear physics calculations needed for the @xmath11 conversion studies involve mainly the integrals of eq . ( [ v.1 ] ) . for the currently interesting nuclei @xmath112 , @xmath113 and @xmath114 the results for @xmath115 and @xmath116 are shown in table 1 . they have been calculated using proton densities @xmath117 from the electron scattering data @xcite and neutron densities @xmath118 from the analysis of pionic atom data @xcite . we have employed an analytical form for the muon wave function obtained by solving the schrdinger equation using the coulomb potential produced by the charge densities discussed before . in this way the nucleon finite size was taken into consideration . moreover , we included vacuum polarization corrections as in ref . @xcite . in solving the schrdinger equation we have used modern neural networks techniques @xcite which give the wave function @xmath119 in the analytic form of a sum over sigmoid functions . thus , in eq . ( [ v.1 ] ) only a simple numerical integration is finally required . to estimate the influence of the non - relativistic approximation on the muon wave function @xmath120 , we have also determined it by solving the dirac equation . the results do not significantly differ from those of the schrdinger picture .    in table 1 we also show the muon binding energy @xmath121 ( obtained as output of the dirac and schrdinger solution ) and the experimental values for the total rate of the ordinary muon capture @xmath122 taken from ref . @xcite .    using values for @xmath115 and @xmath116 for a set of nuclei throughout the periodic table @xcite one can estimate the nuclear structure dependence of the quantities @xmath73 and @xmath101 . in table 2 we show the results obtained for the ratio @xmath102 of eq . ( [ phi ] ) and its approximate expression @xmath123 of eq . ( [ phitil ] ) . we see that @xmath102 is very small ( @xmath124 ) and , because the terms of the quantity @xmath125 depending on @xmath102 are also small , we conclude that the nuclear dependence of @xmath125 can be ignored as we have discussed in sect . 2 .    since there is at maximum @xmath126 difference between @xmath102 and @xmath123 and because the isovector terms of the quantity @xmath125 are small , especially for light nuclear systems , one can also use @xmath127 in the expression ( a , z ) dependence of @xmath125 . the results of table 1 can be exploited for setting constraints on the parameters of a specific gauge model predicting the @xmath11 process . in table 3 we quote the upper bounds for the quantities @xmath73 and @xmath125 derived by using the recent experimental data on the branching ratio @xmath2 given in eqs . ( [ rme - ti ] ) , ( [ rme - pb ] ) and the expected experimental sensitivity of the brookhaven experiment , @xmath128 [ see eq . ( [ rme - al ] ) ] . the limits of @xmath73 and @xmath101 quoted in table 3 are improvements by about four orders of magnitude over the previous ones . we should stress that the limits on the quantities @xmath73 of eq . ( [ rho ] ) and @xmath101 of eqs . ( [ rme.1 ] ) , are the only constraints imposed by the @xmath11 conversion on its underlying elementary particle physics . one can extract upper limits on the individual lepton flavor violation parameters ( couplings of scalar , vector currents , neutrino masses etc . @xcite ) under certain assumptions like the commonly assumed dominance of only one component of the @xmath11 conversion lagrangian which is equivalent to constrain one parameter or product of the parameters at a time . for example , in the case of susy mechanisms , in order to learn about the size and the regularities of possible violation of r - parity , one needs some information on the parameters @xmath129 and the sneutrino vacuum expectation values @xmath130 @xcite . using the upper limits for @xmath101 given in table 3 we can derive under the above assumptions the constraints on @xmath131 of eq . ( [ rme.1 ] ) and the products of various parameters . thus , the bounds obtained for the scalar current couplings @xmath132 in the @xmath133-parity violating lagrangian for the @xmath5al target @xcite are @xmath134 . the limit for @xmath132 obtained with the data of ti target @xcite is @xmath135 , i.e. more than two orders of magnitude weaker than the limit of @xmath5al .    with these limits it is straightforward to derive constraints on the parameters of the initial lfv lagrangian . in table 4 we list the most stringent constraints on the products of the trilinear couplings which we obtained from the experimental limit on @xmath11 conversion in @xmath136 [ see eq . ( [ rme - ti ] ) ] and from the expected experimental sensitivity of meco detector using @xmath5al as a target material [ see eq . ( [ rme - al ] ) ] . the corresponding constraints for @xmath137 are significantly weaker and they are not presented in table 4 . in this table @xmath138 denotes a scaling factor defined as @xmath139 which can be used for reconstructing the limits for the other experimental upper limits on the branching ratio @xmath140 . as seen from table 4 , the @xmath11-conversion limits on the products @xmath141 are significantly more stringent than those previously known in the literature @xcite and given in the 2nd column . in refs . @xcite it was shown that except few cases , the constraints on @xmath142 , @xmath143 and @xmath144 obtained from @xmath11 conversion data are better than those derived from any other process . as we have mentioned at the beginning , significantly better improvement on these limits is expected from the ongoing experiments at psi @xcite and even better from the meco experiment at brookhaven @xcite . this would make the @xmath11 conversion constraints better than those from the other processes in all the cases .    before closing we should note that the last four limits for @xmath145 in table 4 originate from the contribution of the strange nucleon sea . these limits are comparable to the other @xmath11 constraints derived from the valence quarks contributions . the transition matrix elements of the flavor violating @xmath6 conversion are of notable importance in computing accurately the corresponding rates for each accessible channel of this exotic process . such calculations provide useful nuclear - physics inputs for the expected new data from the psi and meco experiments to put more severe bounds on the muon - number - changing parameters ( isoscalar couplings , etc . ) determining the effective currents in various models that predict the exotic @xmath6 process . in the case of the r - parity violating interactions discussed here we have studied all the possible tree - level contributions to the @xmath11 conversion in nuclei taking into account the nucleon and nuclear structure effects . we found new important contribution to @xmath146 originating from the strange quark sea in the nucleon which is comparable with the usual contribution of the valence @xmath147 quarks .    from the existing data on @xmath2 in @xmath136 and @xmath137 and the expected sensitivity of the designed meco experiment @xcite we obtained stringent upper limits on the quantities @xmath73 and @xmath125 introduced in eqs . ( [ rho ] ) and ( [ rme.1 ] ) . they can be considered as theoretical sensitivities of a particular experiment to the @xmath6 conversion signals for various targets employed . thus , these quantities are helpful for comparing different @xmath11 conversion experiments . we also extracted the upper limits on the products of the trilinear parameters of the type @xmath141 which are significantly more severe than those existing in the literature . as to the other products , the following observation was formulated in refs . @xcite . if the ongoing experiments at psi @xcite and brookhaven @xcite will have reached the quoted sensitivities in the branching ratio @xmath2 then the @xmath11 constraints on all the products of the parameters @xmath148 , @xmath145 , @xmath149 measurable in @xmath11 conversion will become more stringent than those from any other processes . this is especially important in view of the fact that no comparable improvements of the other experiments testing these couplings are expected in the near future . marciano and a.i . sanda , phys . lett . * 38 * , 1512 ( 1977 ) ; phys . * b 67 * , 303 ( 1977 ) . vergados , phys . reports * 133 * , 1 ( 1986 ) . w. molzon , the improved tests of muon and electron flavor symmetry in muon processes , spring . , in press . kosmas , g.k . leontaris and j.d . vergados , prog . part . . phys . * 33 * , 397 ( 1994 ) . c. dohmen _ _ , ( sindrum ii collaboration ) , phys . lett . * b 317 * ( 1993 ) 631 ; h.k . walter , phys . * 61 * ( 1998 ) 1253 . w. honecker _ et al . _ , ( sindrum ii collaboration ) , phys . rev . lett . * 76 * ( 1996 ) 200 . s. ahmad _ et al . _ , ( triumf collaboration ) , phys . rev * 59 * ( 1987 ) 970 . w. molzon , the meco experiment : a search for @xmath150 with sensitivity below @xmath151 , invited talk at int . conf . on `` symmetries in physics at intermediate and high energies and applications '' , ioannina - greece , sept . 30 - oct . 5 , 1998 . kosmas and j.d . vergados , nucl . phys . * a 510 * ( 1990 ) 641 . t.s . kosmas , j.d . vergados , and a. faessler , phys . * 61 * ( 1998 ) 1161 . kosmas , a. faessler , and j.d . vergados , j. phys . * g 23 * ( 1997 ) 693 ; h.c . chiang , e. oset , t.s . kosmas , a. faessler , and j.d . vergados , nucl . phys . * a 559 * ( 1993 ) 526 . j. schwieger , t.s . kosmas , and a. faessler , phys . lett . * b 443 * ( 1998 ) 7 . r. barbieri and l. hall , phys . * b 338 * ( 1994 ) 212 ; r. barbieri , l. hall , and a. strumia , nucl . * b 445 * ( 1995 ) 219 . kim , p. ko , and d .- lee , phys . rev . * d 56 * ( 1997 ) 100 , hep - ph/9701381 . k. huitu , j. maalampi , m. raidal , and a. santamaria , phys . lett . * b 430 * ( 1998 ) 355 , hep - ph/9712249 . a. faessler , t.s . kosmas , s. kovalenko , and j.d . vergados , hep - ph/9904335 . g. bhattacharyya , nucl . 52 a ( 1997 ) 83 , ifup - th 43/97 , hep - ph/9709395 ; h. dreiner , in proc . of _ perspectives on supersymmetry _ , ed . kane , world scientific , 462 ( 1999 ) ; p. roy , plenary talk given at the _ pacific particle physics phenomenology conference _ , seoul , korea , oct . 31 - nov . 2 , 1997 , hep - ph/9712520 ; r. barbier et al . report of the gdr working group on the r - parity violation , hep - ph/9810232 . a. faessler , s. kovalenko , f. imkovic , and j. schwieger , phys . * 78 * ( 1997 ) 183 ; a. faessler , s. kovalenko , f. simkovic , phys . rev . * d 57 * , 55004 ( 1998 ) . t. siiskonen , j. suhonen , and t.s . kosmas , phys . rev . * c 60 * ( rapid communication ) , 62501 ( 1999 ) . h. de vries , c.w . de jager and c. de vries , atomic data and nuclear data tables , * 36 * ( 1987 ) 495 . gibbs and b.f . gibson , ann . nucl . part . * 37 * ( 1987 ) 411 ; c. garcia - recio , j. nieves and e. oset , nucl . phys . * a 547 * ( 1992 ) 473 . i.e. lagaris , a. likas , and d.i . fotiadis , comp . * 104 * ( 1997 ) 1 . t. suzuki , _ , phys . rev . * c 35 * ( 1987 ) 2212 . bednyakov , a. faessler , and s. kovalenko , phys . lett . * b 442 * ( 1998 ) 203 . t. siiskonen , t.s . kosmas , and j. suhonen , phys . rev . * c * , to be submitted . .transition matrix elements ( muon - nucleus overlap integrals @xmath152 of eq . ( [ v.1 ] ) ) evaluated by using the exact muon wave function obtained via neural networks techniques . other useful quantities ( see text ) are also included . [ cols=">,^,^,^,^,^ " , ]
Abstract:

Answer: the exotic neutrinoless @xmath0 nuclear conversion is studied within the conventional extensions of the standard model as well as in the minimal supersymmetric ( susy ) models with and without r - parity conservation . the dependence of the @xmath1 conversion rates on the nucleon and nuclear structure is consistently taken into account . using our calculated transition matrix elements and the available experimental data on the branching ratio @xmath2 for @xmath3ti and @xmath4pb as well as the expected experimental sensitivity for @xmath5al employed as a target in the planned at brookhaven @xmath6 conversion ( meco ) experiment , we extract very severe constraints for the flavor violation parameters . we especially emphasize on the constraints resulting for susy r - parity violating parameters .    6.5 in    # 1    # 10= -.015em0 - 0 .03em0 - 0 -.015em.0233em0    * lepton flavor violation : constraints from exotic muon to electron conversion * +    t.s . kosmas @xmath7 and sergey kovalenko @xmath8 + @xmath7__division of theoretical physics , university of ioannina gr-45110 ioannina , greece _ _ + @xmath9__departamento de fsica , universidad tcnica federico santa mara , casilla 110-v , valparaso , chile _ _


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: in recent years there has been a lot of interest in the field of plasmonics which seems to be the only viable path towards realization of nanophotonics : control of light at scales substantially smaller than the wavelength @xcite . however , plasmonic materials ( most notably metals ) suffer from large losses in the frequency regimes of interest , which resulted in a wide search for better materials @xcite . lots of attention has recently been given to plasmonics in graphene @xcite , which is a single two - dimensional ( 2d ) plane of carbon atoms arranged in a honeycomb lattice @xcite . one exciting point of interest of 2d materials is that they are tunable . for example , graphene can be doped to high values of electron or hole concentrations by applying gate voltage @xcite , much like in field effect transistors . furthermore , graphene can be be produced in very clean samples with large mobilities ( demonstrated by dc transport measurements ) @xcite . the dc scattering rates would imply small plasmon losses in graphene , however , it is still not clear how the scattering rates change with frequency , particularly in the infrared ( ir ) region . recent nano - imaging measurements @xcite have demonstrated somewhat increased plasmon losses at ir compared to the estimate based on dc transport measurements . measurements of optical transmission through graphene nano - ribbons @xcite have demonstrated strong increase of plasmon linewidth with frequency and losses that are much larger than the dc estimates . however , since the ribbon width in these experiments is very small ( 10 - 100 nm ) edge scattering can significantly increase the losses . nevertheless , a similar experiment @xcite with graphene nano - rings has demonstrated plasmon linewidths that approximately agree with the dc estimate . finally electron energy loss experiments ( eels ) @xcite on graphene sheets have demonstrated huge plasmon linewidths that increase linearly with plasmon momentum ; however , the ( dc ) transport measurements were not reported so it is not clear what was the actual quality of the graphene films . it is also interesting to note that similar results @xcite were obtained with eels on the mono - atomic silver film which could imply a common origin of plasmon damping in these two 2d systems . on one hand , metallic monolayers might be even more interesting from the point of view of plasmonics since they have abundance of free electrons even in the intrinsic case , while graphene has to be doped with electrons since it is a zero band gap semiconductor . on the other hand , graphene has superior mechanical properties and was demonstrated in a free standing ( suspended ) samples while metallic monolayers have only been observed on a substrate . instead of calculating plasmon linewidth , we will focus on a directly related problem of optical absorption , which is easier to analyze . in that respect , it was shown experimentally @xcite that suspended graphene absorbs around 2.3@xmath0 of normal incident light in a broad range of frequencies . however , if graphene is doped with electrons , then pauli principle blocks some of these transitions and there should be a sudden decrease of absorption below a certain threshold , which should theoretically occur at twice the fermi energy . nevertheless , optical spectroscopy experiments @xcite have shown that there is still lots of absorption even below this threshold . this absorption is much larger than the estimate based on dc measurements . a great deal of theoretical work addressed this problem @xcite , but to our knowledge , the experimental results have quantitatively not been explained yet .    in this paper , we focus on optical absorption mediated by charged impurity scattering . as we have already stated , the motivation for studying this problem follows from the fact that typical graphene samples can have large mobilities ( @xmath1 @xmath2vs ) in spite of a huge concentration of charged impurities @xcite ( @xmath3 @xmath4 ) , which is actually comparable to the typical concentration of electrons . the reason one can have such a large mobility is screening @xcite . in fact , if one assumes that electrons scatter from bare charged impurities described with the coulomb potential @xmath5 , then the resulting mobility is almost two orders of magnitude lower than the measured value @xcite . the only way to reconcile the experiment and theory is to say that the actual scattering potential is screened to @xmath6 , where @xmath7 is the static dielectric function . however , in the dynamical case , at finite frequency , screening is not so effective and @xmath7 should be replaced with the dynamic dielectric function @xmath8 . this will certainly influence the single particle excitations where an incident photon excites an electron hole pair through impurity scattering . moreover , at finite frequency one can have @xmath9 ( at the plasmon dispersion ) so there exists an additional decay channel where an incident photon excites a plasmon of the same energy , through impurity scattering . in other words , impurities break the translational symmetry ( momentum does not need to be conserved ) , which allows the photon to couple directly to a plasmon mode . very recently another group also calculated this process in graphene but only in the small frequency limit @xcite . here we give the result for the arbitrary frequency ( both for metallic monolayers and graphene ) which can be very different from the small frequency limit . more specifically , we calculate the optical absorption in the 2d electron systems with the randomly arranged charged impurities . first , we discuss the case of metallic monolayers which have a parabolic electron dispersion , and then the case of graphene with dirac electron dispersion . we focus on a decay channel where the incident photon emits a plasmon through impurity scattering , but we also discuss a case where the incident photon emits the plasmon and a surface optical phonon of the substrate . for graphene on @xmath10 substrate , the resulting optical absorption is very small compared to the experimental results @xcite , and not enough to reconcile the difference between the theory @xcite and the experiment@xcite . on the other hand we predict large optical absorption by plasmon emission via impurity scattering in suspended graphene . thus we believe that these ideas can be tested in suspended graphene . finally we note that for suspended graphene ( metallic monolayers ) the small frequency limit @xcite gives an order of magnitude lower ( larger ) result than the more exact rpa calculation . the case of the optical absorption in a bulk 3d system with parabolic electron dispersion and randomly arranged impurities was already studied by hopfield @xcite . it is straightforward to extend his result to a 2d system and here we provide only a brief description of the calculation . we study a system described by the hamiltonian @xmath11 , where @xmath12 represents kinetic energy of free electrons , @xmath13 describes electron - electron interaction which is conveniently represented through the screening effect , @xmath14 describes scattering with light , and @xmath15 scattering with impurities . electrons in a metallic monolayer can be described with a parabolic dispersion : @xmath16 , where @xmath17 is the electron momentum , and @xmath18 is effective mass of the electron . next , let us introduce a monochromatic light beam of frequency @xmath19 which is described by the electric field @xmath20 . this wave is incident normally on a 2d electron gas , that is @xmath21 is in the plane of the gas . if we are only interested in a linear response with respect to this electric field , then interaction of electrons with light takes a particularly simple expression : @xmath22 , where we have introduced electron charge ( @xmath23 ) . further on , since momentum is a good quantum number even in an interacting electron system , light scattering ( @xmath14 ) will not change the many - body eigenstates of @xmath24 , but only the eigenvalues , see ref . @xcite . then , one only needs to do the perturbation theory in the impurity scattering . unfortunately this trick ( due to hopfield ) works only in the systems with parabolic electron dispersion , while in the case of dirac electrons , like those found in graphene , one needs to do the perturbation theory both in the light scattering and in the impurity scattering , which is a much more tedious task . we can write the hamiltonian for impurity scattering as a fourier sum over wavevectors @xmath25 : @xmath26 , where @xmath27 is the total area of our 2d system . by calculating the induced current to the second order in @xmath28 one can find the real part of conductivity @xcite : @xmath29 note that this quantity ( @xmath30 ) determines the optical absorption in our system . here @xmath8 stands for a dielectric function of the electron gas and @xmath31 is the fourier transform of the coulomb potential between two electrons in 2d layer embedded between two dielectrics of average relative permitivity @xmath32 . we have assumed without loss of generality that the external field points in the @xmath33 direction ( @xmath34 ) and is parallel to the plane of our 2d electron gas .    in the case of randomly assembled impurities at positions @xmath35 , one can write for the scattering potential @xmath36 . note that we are assuming positively charged ( @xmath37 ) impurities embedded in a see of negative electrons ( @xmath23 ) . then by averaging over random impurity positions one has @xmath38 , where @xmath39 is the number of impurities @xcite . equation ( [ absorp2d ] ) depends on the loss function @xmath40 which generally contains contribution from single particle excitations and collective ( plasmon ) excitations . in this paper we focus solely on the plasmon contribution in which case one can write @xcite : @xmath41 where @xmath42 is the plasmon frequency determined by the zero of the dielectric function : @xmath43 . this term then represents the process where an incident photon excites plasmon of the same energy , through impurity scattering . the @xmath44-function from equation ( [ lossf_plasmon_omega ] ) extracts only a single wavevector from the sum in equation ( [ absorp2d ] ) , which corresponds to the plasmon wavevector at the given frequency @xmath19 . then one is left with integration over the angle @xmath45 which is straightforward to perform since @xmath46 . finally , we plot the conductivity from expression ( [ absorp2d ] ) in figure [ fig1 ] by using the dielectric function @xmath8 within random phase approximation ( rpa ) given in ref . @xcite . to represent the experiment @xcite , which studied silver monolayer on a silicon substrate , we choose @xmath47 , @xmath48 , the effective mass @xmath49 , where @xmath50 is the free electron mass , electron concentration @xmath51 @xmath4 , and we assume the impurity concentration @xmath52 @xmath4 . it is also convenient to look at the small frequency limit @xmath53 in which case only long wavelength @xmath54 plasmons contribute to the scattering . here @xmath55 and @xmath56 stand for fermi energy and fermi momentum , respectively . in this limit , one can use a simple drude model to obtain the dielectric function : @xmath57 in this case , plasmon dispersion is simply @xmath58 and one can easily evaluate equations ( [ absorp2d ] ) and ( [ lossf_plasmon_omega ] ) to obtain the conductivity : @xmath59 here we have introduced the thomas - fermi wavevector : @xmath60 , while @xmath61 stands for the impurity density . from fig . [ fig1 ] b we see that in the case of metallic monolayers the small frequency limit ( dashed line ) significantly overestimates the more exact rpa result ( solid line ) . unfortunately , the trick that hopfield used in the case of the parabolic dispersion does not work for dirac dispersion so one has to do the perturbation theory both in impurity scattering and in light scattering , while including the screening effect in every order of the perturbation theory . this is straightforward , but very tedious task , so we give the derivation of the optical absorption in the appendix . here we only write the final result : @xmath62 where we have assumed general impurity scattering hamiltonian : @xmath63 ( see the appendix for more details ) . in the case of charged impurities one has @xmath38 after averaging over random impurity positions . then , to find the contribution of plasmon emission process one can use equation ( [ lossf_plasmon_omega ] ) and the dielectric function which is calculated in ref . @xcite within the rpa . the resulting optical absorption , is plotted in figure [ fig2 ] . to resemble parameters from the experiment @xcite we choose electron concentration @xmath64 @xmath4 , and impurity concentration @xmath52 @xmath4 . furthermore , we plot the case of graphene sitting on the @xmath10 substrate where @xmath65 , but also the case of suspended graphene where @xmath66 . it is also convenient to look at the small frequency limit @xmath53 in which case only long wavelength @xmath54 plasmons contribute to the scattering . then , one can use a simple drude model to obtain the dielectric function in graphene @xcite : @xmath67 in this case the function @xmath68 takes a particularly simple expression ( see the appendix for more details ) : @xmath69 , and it is straightforward to evaluate expression ( [ absorp_graphene ] ) to obtain : @xmath70 note that this is the same result as in the case of metallic monolayers . this is expected because in the small frequency ( long wavelength ) limit , one does not expect to see specific details of the band structure . of course , in the graphene case , the thomas - fermi wavevector is given by a different expression : @xmath71 . we would like to note that the small frequency limit in the case of graphene was also recently obtained by another group @xcite . however , from figure [ fig3 ] one can see that the small frequency limit can be very different from the more general rpa result . if we now compare our results [ figure [ fig2 ] ( b ) ] with experiment @xcite , we see that this effect of plasmon emission is relatively small ( @xmath72 ) compared to the experimental results ( @xmath73 ) in this regime . one might ask what are the other potentially strong scattering mechanisms ? for example in experiment @xcite , graphene is sitting on @xmath10 , which is a polar substrate , so there is a strong interaction of electrons with the surface polar phonons at energy @xmath74 ev . this is described by the hamiltonian @xmath75 , where @xmath76 is the phonon creation operator . for the square of the scattering potential we can write @xcite : @xmath77 . we use parameters from ref . @xcite for @xmath78 scattering : @xmath79 , @xmath80 , and we assume that the van der waals distance between graphene and the substrate is @xmath81 nm . if we neglect the frequency dependence of @xmath82 , one can make an estimate of absorption simply by replacing @xmath83 with @xmath84 in relation ( [ absorp_graphene ] ) . strictly speaking this is valid only at large frequencies when @xmath85 , but it should give a reasonable estimate in the regime @xmath86 which is the relevant regime in experiment @xcite . the resulting absorption is still extremely small ( @xmath87 ) in the regime of interest ( @xmath88 ) . even though our analysis suggests that these loss mechanisms can not be distinguished from other loss mechanisms in current experiments involving graphene on a @xmath10 substrate , our calculations point out that they should be observable in suspended graphene ( see figure [ fig3 ] ) . suspended graphene is a much cleaner system as one can eliminate all the scattering mechanisms that originate from the interaction with the substrate . moreover , in optical transmission measurements on suspended graphene ( sketched in figure [ fig3 ] ( c ) ) one does not need to consider optical absorption of the substrate . suspended graphene can be doped by depositing electron - donor atoms like sodium or lithium . in that case one is left with impurity ions with the same number as the number of injected electrons . in figure [ fig3 ] ( b ) we plot optical absorption in suspended graphene for identical impurity and electron concentrations @xmath89 @xmath4 . one can see that there is a huge optical absorption through the plasmon emission channel as the real part of conductivity reaches @xmath73 . this would correspond to the @xmath90 reduction in the intensity of transmitted light , which could easily be observed as the @xmath91 reduction is already visible by naked eye @xcite . finally we note that the small frequency limit ( equation ( [ absorp2d_drude_graphene ] ) ) underestimates the more exact rpa calculation ( equation ( [ absorp_graphene ] ) ) by an order of magnitude . in conclusion , we have studied optical absorption of 2d electron gas in graphene and metallic monolayers with random distribution of charge impurities . this formalism can also treat other 2d electron systems like those found in heterostructures , single layer boron - nitride , or single layer molybdenum - disulphide where we expect similar behavior . specifically , we have focused on a decay channel where an incident photon excites a plasmon through impurity scattering . for the graphene sitting on a @xmath10 substrate , we have also studied a decay channel where an incident photon excites a plasmon and an optical phonon of the polar substrate . the resulting optical absorption is more than one order of magnitude lower than the experimental results @xcite , and not enough to reconcile the difference between the theory @xcite and the experiment @xcite . on the other hand we predict large optical absorption by plasmon emission via impurity scattering in suspended graphene . thus we believe that these ideas can be tested in suspended graphene . finally we note that for suspended graphene ( metallic monolayers ) the small frequency limit @xcite gives an order of magnitude lower ( larger ) result than the more exact rpa calculation . we use single particle density matrix ( spdm ) approach which is a convenient way to take into account both temperature and the pauli principle . equation of motion for spdm @xmath92 is given by @xcite : @xmath93 , \label{eq_of_motion}\ ] ] where the hamiltonian is given by @xmath94 here @xmath12 represents kinetic energy of free electrons , @xmath14 describes scattering with light , @xmath15 scattering with impurities , and @xmath95 describes electron - electron interactions which we only take in the form of a self - consistent screening field . in the case of graphene , electrons are described by dirac dispersion @xcite : @xmath96 where @xmath97 m / s is fermi velocity , @xmath98 is electron wavevector , @xmath99 , and @xmath100 are the pauli spin matrices . let us denote by @xmath101 eigenstates of @xmath12 , where @xmath102 stands for the conduction band , and @xmath103 for the valence band . then the eigenvalues of @xmath12 are given by dirac cones : @xmath104 . if we now introduce a light source described by the electric field @xmath105 , then scattering with light is determined by the hamiltonian : @xmath106 where @xmath23 is the electron charge . furthermore , we can write the hamiltonian for impurity scattering as a fourier sum over wavevectors @xmath25 : @xmath107 where @xmath27 is total area of our graphene flake , @xmath108 is the position operator , and @xmath83 is the fourier transform of the scattering potential . here we assume a general scattering potential and only later we will specify @xmath83 for the case of charged impurity scattering and surface polar phonon scattering . finally , one can also write the screening field as a fourier sum : @xmath109 but one has to keep in mind that different orders of the perturbation expansion will have different time dependence ( frequencies ) . here , the screening field is taken as a self - consistent electrostatic field that the electrons induce on themselves , so one can write @xmath110 , where @xmath111 is the fourier transform of the electron density , and @xmath112 is the fourier transform of the coulomb potential between two electrons . for a 2d electron gas embedded between two dielectrics of relative permitivity @xmath32 , one can write : @xmath113 . note that this is valid only in the electrostatic limit @xmath114 which is the relevant regime for our case . furthermore , since @xmath115 , one can write for the screening field : @xmath116 where we have taken into account 2 spin and 2 valley degeneracies . we are now interested in calculating the current response up to the linear order in the external electric field @xmath21 . since the electric field is uniform in the graphene plane , we are only interested in the @xmath117 term , and the current density operator is given by @xmath118 . the induced current will have only the @xmath33 component , since the electric field points in the @xmath33 direction . finally , the induced current density is given by @xmath119 , so we can write : @xmath120 to include also impurity scattering , we need to calculate the induced current up to the second order in @xmath28 . in other words we need to do a perturbation expansion of spdm : @xmath121 where @xmath122 is the equilibrium solution to equation ( [ eq_of_motion ] ) for independent dirac electrons in the absence of impurity scattering and light scattering , @xmath123 is solution of equation ( [ eq_of_motion ] ) correct up to a linear order in light scattering , @xmath124 is solution up to the linear order in impurity scattering , @xmath125 is solution up to linear order both in light scattering and impurity scattering , and @xmath126 is solution up to the linear order in light scattering and quadratic in impurity scattering . using equation ( [ eq_of_motion ] ) , we can now write the equation of motion for every order of spdm expansion : @xmath127 , \label{eq_of_motion_0}\ ] ] @xmath128 + \left [ h_l+h_l^s,\rho_0 \right ] , \label{eq_of_motion_l}\ ] ] @xmath129+\left [ h_i+h_i^s,\rho_0 \right ] , \label{eq_of_motion_i}\ ] ] @xmath130 + \left [ h_i+h_i^s,\rho_l \right]+ \left [ h_l+h_l^s,\rho_i \right]+ \left [ h_{li}^s,\rho_0 \right ] , \ ; \mbox{and } \label{eq_of_motion_li}\ ] ] @xmath131 + \left [ h_i+h_i^s,\rho_{li } \right]+ \left [ h_{li}^s,\rho_i \right]+ \left [ h_{lii}^s,\rho_0 \right ] . \label{eq_of_motion_lii}\ ] ] the equilibrium solution of equation ( [ eq_of_motion_0 ] ) describes the free electrons and is given by : @xmath132 where @xmath133 is the kronecker delta symbol and @xmath134^{-1}$ ] is the fermi - dirac distribution at temperature @xmath135 and fermi energy @xmath55 . using relation ( [ rho_0 ] ) we can write the solution of equation ( [ eq_of_motion_l ] ) as : @xmath136 which is a stady - state solution of spdm that oscillates at frequency @xmath19 . here we have used the following relation : @xmath137 . we have neglected the screening field @xmath138 in equation ( [ eq_of_motion_l ] ) since the 2d electron gas can not screen the uniform electric field . this can be seen below from equation ( [ epsilon_approx ] ) which gives the dielectric function of graphene in the long wavelenght limit . one can immediately see that @xmath139 which means that there is no screening in the @xmath140 limit . let us now focus on equation ( [ eq_of_motion_i ] ) . we can introduce a self - consistent hamiltonian @xmath141 , and write : @xmath142 , where @xmath143 is a self - consistent scattering potential that consists of a bare impurity scattering potential @xmath144 , and a screening field @xmath145 . by solving equations ( [ screening_field_tr ] ) and ( [ eq_of_motion_i ] ) in a self - consistent way one can show that @xmath146 , where @xmath7 is the static dielectric function . the dynamic dielectric function of graphene is generally @xmath147 and one can simply check that @xmath148 . finally , the solution to equation ( [ eq_of_motion_i ] ) can be written as @xmath149 to solve the next order of perturbation theory @xmath150 we need to include the screening field described by a hamiltonian @xmath151 . one can then solve equation ( [ eq_of_motion_li ] ) by using results ( [ rho_0 ] ) , ( [ rho_l ] ) and ( [ rho_i ] ) to obtain : @xmath152 next , one can use relation ( [ screening_field_tr ] ) to obtain the screening field in a self - consistent way : @xmath153 where @xmath8 is the dynamic dielectric function given in ( [ dielectric_function ] ) . note that the terms containing @xmath154 and @xmath155 have disappeared after summation over @xmath156 and @xmath98 . one can also demonstrate the following important property : @xmath157 . finally , one can use equation ( [ eq_of_motion_lii ] ) to find @xmath158 , and equation ( [ current_tr ] ) to find the induced current up to the first order in light scattering , and the second order in impurity scattering : @xmath159 note that we have neglected the screening field @xmath160 since we need only the @xmath140 component of @xmath158 to obtain @xmath161 , and there is no screening in the 2d electron gas in the @xmath140 case . also note that we have skipped the lower orders in the induced current since one can generally show that @xmath162 . on the other hand @xmath163 but we are here interested in the optical absorption below interband threshold @xmath164 , where @xmath165 . finally , to evaluate the current component @xmath161 from expression ( [ current_lii ] ) we need to use expression ( [ rho_i ] ) for @xmath166 and expression ( [ rho_li ] ) for @xmath150 . the resulting conductivity is : @xmath167 where the functions @xmath68 , @xmath168 and @xmath169 are given by the following expressions : @xmath170 @xmath171 @xmath172 however , if we are interested only in the contribution from the collective excitations , we can neglect the single particle excitations to obtain : @xmath173 note that this is the complete expression for the real part of conductivity , i.e. @xmath174 since @xmath175 in this regime , and generally @xmath176 . then , since we are only interested in the plasmon contribution one can write the loss function as @xmath177 where @xmath42 is plasmon frequency at a given wavevector @xmath178 , and @xmath179 is plasmon wavevector at a given frequency @xmath19 , which is determined by the zero of the dielectric function : @xmath180 . note that @xmath44-function from equation ( [ loss_function_plasmon ] ) extracts only a single wavevector from the integral @xmath181 in equation ( [ absorb_plasmon_contrib ] ) . moreover , one can explicitly perform the remaining integral @xmath182 . to demonstrate this we start by writing the expression for the dirac wave function in coordinate representation : @xmath183 it is straightforward to calculate the following matrix elements : @xmath184 @xmath185 @xmath186 furthermore , the product of the last three terms can be written as : @xmath187 + i \frac{\sin\varphi}{2 }   \left [ n '' e^{i\varphi_{\bf q } } - n e^{-i\varphi_{\bf q } } \right ] \right ) , \label{product_matrix}\end{aligned}\ ] ] where @xmath188 . finally one can show that : @xmath189 where @xmath190 depends only on the magnitude of the wavevector @xmath178 and is given by the following expression : @xmath191 now one can indeed see that the that integration over @xmath192 in equation ( [ absorb_plasmon_contrib ] ) simply contributes with the following factor : @xmath193 . finally , equation ( [ absorb_plasmon_contrib ] ) is reduced to the following expression : @xmath194 where @xmath178 is the plasmon wavevector at the frequency @xmath19 . to evaluate this expression one needs to calculate the double integral @xmath195 to evaluate the function @xmath190 . this can be further simplified at zero temperature when the fermi - dirac distrubution is a step function . in that case , we can group @xmath196 and @xmath197 terms in equation ( [ f_skalped ] ) to obtain : @xmath198 where @xmath199 stands for the fermi - dirac distribution of the conduction band , and we have assumed electron doping i.e. @xmath200 . we perform a numerical integration to evaluate the function @xmath190 ; however , one can obtain a closed expression in the small frequency limit when @xmath201 . in that case , only intraband transitions contribute and one can set @xmath202 in equation ( [ f_skalped_simple ] ) . furthermore , in that case plasmon wavevector @xmath178 is much smaller than the fermi wavevector @xmath56 so one can use the long wavelength expansions : @xmath203 @xmath204 next , it is straightforward to perform integration in equation ( [ f_skalped_simple ] ) to obtain long wavelength ( small frequency ) approximation : @xmath205 in a similar manner , from equation ( [ dielectric_function ] ) , one can obtain the dielectric function in this approximation : @xmath206 which is just the drude model for the dielectric function in graphene @xcite . this work was supported in part by the unity through knowledge fund ( ukf grant no 5/13 ) . the work of m. soljai was supported in part by the mit s3tec energy research frontier center of the department of energy under grant de - sc0001299 . this work was also supported in part by the army research office through the institute for soldier nanotechnologies under contract w911nf-13-d-0001 . z. fei , a.s . rodin , g.o . andreev , w. bao , a.s . mcleod , m. wagner , l.m . zhang , z. zhao , m. thiemens , g. dominguez , m.m . fogler , a.h . castro - neto , c.n . lau , f. keilmann , d.n . basov , _ gate - tuning of graphene plasmons revealed by infrared nano - imaging , _ nature 487 , 82 ( 2012 )      z. fang , s. thongrattanasiri , a. schlater , z. liu , l. ma , y. wang , p. m. ajayan , p. nordlander , n. j. halas , f.j . garcia de abajo , _ gated tunability and hybridization of localized plasmons in nanostructured graphene _ , acs nano * 7 * , 2388 ( 2013 )
Abstract:

Answer: screening is one of the fundamental concepts in solid state physics . it has a great impact on the electronic properties of graphene where huge mobilities were observed in spite of the large concentration of charged impurities . while static screening has successfully explained dc mobilities , screening properties can be significantly changed at infrared or optical frequencies . in this paper we discuss the influence of dynamical screening on the optical absorption of graphene and other 2d electron systems like metallic monolayers . this research is motivated by recent experimental results which pointed out that graphene plasmon linewidths and optical scattering rates can be much larger than scattering rates determined by dc mobilities . specifically we discuss a process where a photon incident on a graphene plane can excite a plasmon by scattering from an impurity , or surface optical phonon of the substrate .


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: numerous complex systems in nature and society can be simplified and abstracted by describing them as networks , in which nodes and edges represent constituent elements and their interactions , respectively . extensive studies @xcite have revealed common statistical features of real - world complex networks , such as the small - world property @xcite , the scale - free property @xcite , community structures @xcite , and degree - degree correlations @xcite . in order to clarify the origin of these features and/or properties of various dynamics on such networks , so many network models have been proposed so far . in most of previous network models , the number of nodes @xmath4 , namely the network size , is treated as an a priori given parameter . the network size can then take any value , and , as is often the case , the limit of infinite @xmath4 is taken in order to simplify the analysis . thus , these models implicitly assumes that networks are stably present no matter how large networks grow . this assumption is , however , not always valid in real - world systems . in an ecological network representing a closed ecosystem , for example , too many species destabilize the ecosystem and the number of species ( nodes in the ecological network ) can not increase unboundedly @xcite . also in a trading network , too many firms make the network fragile because of unstable cartels @xcite , increase of financial complexity @xcite , possible large - scale chain - bankruptcy , and other risks @xcite . as in these examples , due to intrinsic instability of large - scale networks , some sort of networks have their own limit in sizes only below which they can be stable @xcite . it is important to study limit sizes of connected networks and find a way to control them . such information for _ functional networks _ is particularly crucial , because unstable functional networks are directly connected to the instability of our modern society supported by them . functions provided by functional networks are guaranteed by global connectivity and normal operation of each node ( or edge ) . however , the larger a network grows , the lower the probability that all the nodes operate normally becomes . when failures are caused by _ overloads _ , even a few failures can spread to the entire network through a cascading process , which leads the fatal breakdown of the network @xcite . it is then significant to investigate how large functional networks can grow , while maintaining its global connectivity , by overcoming cascading overload failures . in this paper , we evaluate the upper limit of connected network size above which the network becomes unstable by cascading overload failures and discuss how the maximum stable size can be controlled . we examine uncorrelated random networks with poisson and power - law degree distributions by employing the model of cascading failures triggered by temporally fluctuating loads @xcite . the rest of this paper is organized as follows . in sec . [ sec : model_method ] , we briefly review the model of cascading overload failures proposed by ref . @xcite and explain how we analyze the stability of networks against cascading overload failures . our analytical results and their numerical confirmations are presented in sec . [ sec : results ] . section [ sec : sammary ] is devoted to the summary and some remarks . in this section , we outline the model of cascading failures induced by temporally fluctuating loads @xcite . in functional networks such as power grids , the internet , or trading networks , some sort of  flow " ( electric current in a power grid , packet flow in the internet , and money flow in a trading network ) realizes their functions . and flow , at the same time , plays a role of  loads " in these networks . the load on a node usually fluctuates temporally and the node fails if the instantaneous value of the load exceeds the node capacity . since flux fluctuations at a node exhibit the same scaling behavior with fluctuations of the number of non - interacting random walkers on the node @xcite , kishore _ modeled fluctuating loads by random walkers moving on a network and calculated the overload probability that the number of walkers exceeds the range allowed for a node @xcite . the model of cascading failures employed in the present work utilizes this overload probability .    in a connected and undirected network with @xmath5 edges , the probability @xmath6 that @xmath7 random walkers ( loads ) exist on a node of degree @xmath8 is presented by @xmath9 where @xmath10 is the total number of walkers and @xmath11 is the stationary probability to find a random walker on a node of degree @xmath8 @xcite . this leads a natural definition of the capacity @xmath12 as @xmath13 where @xmath14 is a real positive parameter characterizing the node tolerance , and @xmath15 and @xmath16 are the average and the standard deviation of @xmath6 , which are given by @xmath17 and @xmath18 , respectively . since the overload probability @xmath19 is the probability of @xmath7 to exceed @xmath12 , we have @xcite @xmath20 where @xmath21 is the regularized incomplete beta function @xcite and @xmath22 denotes the largest integer not greater than @xmath23 .    using the above overload probability , the cascade process of overload failures is defined as follows @xcite :    \(i ) prepare an initial connected , uncorrelated , and undirected network @xmath24 with @xmath4 nodes and @xmath5 edges , in which totally @xmath10 random walkers exist , and determine the capacity @xmath12 of each node according to eq . ( [ eq : capacity ] ) . @xmath10 is set as @xmath25 , where the parameter @xmath26 is the load carried by a single edge . \(ii ) at each cascade step @xmath27 , reassign @xmath28 walkers to the network @xmath29 at step @xmath27 , where the total load @xmath28 is given by @xmath30 here , @xmath31 is the total number of edges in the network @xmath29 and @xmath1 is a real positive parameter . \(iii ) for every node in @xmath29 , calculate the overload probability given by @xmath32 where @xmath33 and @xmath8 are the initial degree and the degree of the node at cascade step @xmath27 , and remove nodes from @xmath29 with this probability . \(iv ) repeat ( ii ) and ( iii ) until no node is removed in the procedure ( iii ) . the reduction of the total load in the procedure ( ii ) corresponds to realistic situations in which the total load is reduced to some extent during a cascade process to prevent the fatal breakdown of the network function . when a company goes bankrupt on a trading network , for example , a large - scale chain bankruptcy would be prevented by the reduction of the total debt ( loads ) realized by financial bailout measures . the exponent @xmath1 characterizes how quickly the total load decreases with decreasing the network size , which is called the load reduction parameter .    during the cascade , the network @xmath29 might be disconnected even though the initial network @xmath24 is connected . in such a case , a walker on a connected component can not jump to other components . therefore , the amount of walkers on each component is conserved in the random walk process . the overload probability then becomes dependent on how the total load is distributed to disconnected components . thus the overload probability deviates from eq . ( [ eq : ftau ] ) . this deviation is , however , small and the effect of disconnected components can be approximately neglected as argued in details in ref . @xcite . the validity of this approximation will be confirmed in the next section by numerical simulations in which walkers distributed proportionally to the number of edges in each component can not move to other components . we examine the stability of a network under cascading overload failures described above by analyzing the size @xmath34 of the largest connected component in the network after completed the cascading process . if @xmath34 is very small , the initial network is considered to be unstable . the quantity @xmath34 obviously depends on the initial network size @xmath4 , and the maximum value @xmath0 of @xmath34 with respect to @xmath4 provides the upper limit of the size of stable connected networks in a given cascading condition . the surviving component of size @xmath0 may experience further cascading failures after a long time , but simultaneously the component can grow during this period . in the competition between the growth and decay processes , the component smaller than @xmath0 can , in substance , stably grow up to @xmath0 . therefore , the connected network size fluctuates around this maximum size @xmath0 . in order to calculate the maximum stable size @xmath34 , we construct a master equation for the probability @xmath35 that a randomly chosen node has the degree @xmath8 at cascade step @xmath27 and the initial degree @xmath36 . it is convenient to introduce another probability @xmath37 of a node adjacent to a randomly chosen node of degree @xmath8 to experience an overload failure at cascade step @xmath27 . this probability is independent of @xmath8 for uncorrelated networks and is given by @xcite @xmath38 where @xmath39 is the average degree of @xmath29 . we then formulate the master equation for @xmath35 as    @xmath40+\delta_{k0}f_{w_{\tau-1}}(k_{0},k')\biggr\}. \label{eq : pikk}\ ] ]    in this equation , we do not remove overloaded nodes actually but leave them in the system as zero - degree nodes , which makes the theoretical treatment easier . the right - hand side of this equation represents the probability that a degree-@xmath41 node in @xmath42 becomes a node of degree @xmath8 at cascade step @xmath27 . the first term describes the situation that the degree-@xmath41 node does not experience an overload failure and @xmath43 nodes adjacent to this node fail . the second term stands for the case that the degree-@xmath41 node itself fails and becomes a zero - degree node . solving numerically eq . ( [ eq : pikk ] ) with the aid of eq . ( [ eq : phi2 ] ) , we can calculate @xmath35 iteratively starting from @xmath44 , where @xmath45 is the degree distribution function of @xmath24 . according to the procedure ( iv ) , we stop this iterative calculation at step @xmath46 satisfying the condition @xmath47 which implies that the expectation number of overloaded nodes becomes less than unity . we can obtain the largest connected component size @xmath34 at the final cascade step @xmath46 from the degree distribution @xmath48 of @xmath49 which is given by @xmath50 . employing the generating function formalism , @xmath34 is calculated by @xcite @xmath51 , \label{eq : ntau}\ ] ] where @xmath52 is the smallest non - negative solution of @xmath53 and @xmath54 is the generating function of the remaining degree distribution , which is defined by @xmath55 it should be noted that eq . ( [ eq : ntau ] ) does not mean that @xmath34 is proportional to @xmath4 because @xmath56 depends on @xmath4 . first , we calculated @xmath34 for the erds - rnyi random graph ( errg ) as an initial network @xmath24 . in this case , the binomial degree distribution function for @xmath24 is given by @xmath57 where @xmath58 . in this work , we fix the initial average degree as @xmath59 . although initial networks having this average degree are not completely connected with isolated nodes at a very low rate , this does not affect our conclusion . figure [ fig1 ] shows @xmath34 as a function of the initial network size @xmath4 for various values of the load reduction parameter @xmath1 . for these results , the node tolerance parameter @xmath14 and the load carried by a single edge , @xmath26 , are chosen as @xmath60 and @xmath61 . the lines in fig . [ fig1 ] represent @xmath34 calculated by eq . ( [ eq : ntau ] ) and the symbols indicate the results obtained by numerical simulations performing faithfully the cascade process from ( i ) to ( iv ) described in sec . [ subsec : model ] . in the numerical simulation , the overload probability at cascade step @xmath27 is calculated under the condition that random walkers can not jump to other components . namely , instead of eq . ( [ eq : ftau ] ) , we adopt the overload probability of a node in the @xmath62-th component given by @xmath63 where @xmath64 is the number of edges in the @xmath62-th component of @xmath29 and @xmath65 . the remarkable agreement between the symbols and the lines suggests that our approximation by eq . ( [ eq : ftau ] ) is quite accurate .     after completing cascading overload failures as a function of the size @xmath4 of initial errgs , for various values of @xmath1 . the thick lines indicate @xmath34 calculated by eq . ( [ eq : ntau ] ) and filled and open circles on thick lines show the results obtained by the numerical simulation described in the main text . the thin straight line is a guide to the eyes for @xmath66 . ( b ) relative largest component size @xmath67 as a function of @xmath4 . lines have the same meanings as those in ( a ) . all the results are calculated for @xmath68 , @xmath60 , and @xmath61.,scaledwidth=48.0% ]     as a function of the load reduction parameter @xmath1 . three lines represent the results for the errg ( solid line ) and sf networks with the degree distributions @xmath45 given by eq . ( [ cauchy ] ) with @xmath69 ( dashed line ) and @xmath70 ( dotted line ) . the minimum degree @xmath71 in eq . ( [ cauchy ] ) is set as @xmath72 . all the results are calculated for @xmath68 , @xmath60 , and @xmath61.,scaledwidth=48.0% ]    the quantity @xmath34 shown in fig . [ fig1 ] is exactly equal to @xmath4 as long as @xmath73 , regardless of the value of @xmath1 . this implies that the network never experiences overload failures until the network grows up to @xmath74 . thus , @xmath74 is determined by @xmath75 which does not depend on @xmath1 . when the network grows larger than @xmath74 , it starts to decay by initial failures and subsequent avalanche of failures . the largest component size @xmath34 after the cascade then becomes smaller than @xmath4 , but still increases with @xmath4 , at least unless @xmath4 is much larger than @xmath74 . for @xmath76 , when @xmath4 exceeds a certain value @xmath77 , @xmath34 rapidly decreases with @xmath4 . therefore , @xmath34 becomes maximum at @xmath78 . this maximum value of @xmath34 is nothing but @xmath0 mentioned at the beginning of sec . [ subsec : largest ] . figure [ fig1 ] clearly shows that the maximum stable size @xmath0 is an increasing function of @xmath1 . this is because a large value of @xmath1 , namely a rapid decrease of the total load @xmath28 during the cascade , prevents large - scale cascading failures in our model .    the @xmath1 dependence of @xmath0 is closely related to the percolation transition by cascading overload failures . as pointed out by ref . @xcite , there exists a critical value @xmath79 above which the largest component size diverges in proportion to @xmath4 in the thermodynamic limit . thus , @xmath34 goes to infinity as @xmath80 for @xmath2 , which implies the absence ( divergence ) of @xmath0 . on the other hand , for @xmath81 , @xmath34 is finite and varies with @xmath4 to be maximized at @xmath82 as mentioned above . as a consequence , @xmath0 increases with @xmath1 for @xmath81 and diverges at @xmath83 . a finite - size scaling analysis @xcite predicts that @xmath0 for @xmath81 behaves as @xmath84 if @xmath1 is close enough to @xmath79 , where the correlation volume exponent @xmath85 and the order parameter exponent @xmath86 characterize @xmath82 and @xmath87 as @xmath88 for large enough @xmath82 and @xmath89 for large enough @xmath4 , respectively . such a behavior is demonstrated by the solid line in fig . [ fig2 ] for the errg . the result suggests that the load control during a cascade is crucial to realize large functional networks .     for the errg with @xmath59 as a function of the load reduction parameter @xmath1 . four lines represent the results for different combinations of the node tolerance parameter @xmath14 and the load carried by a single edge @xmath26.,scaledwidth=48.0% ]    how is the maximum size @xmath0 affected by properties of the initial network @xmath24 ? to clarify the influence of the degree inhomogeneity in @xmath24 , we first calculate @xmath0 for scale - free ( sf ) networks with the degree distribution given by @xmath90 where @xmath91 , @xmath92 , @xmath93 , and the normalization constant @xmath94 are real positive constants . the degree distribution has asymptotically a power - law form , i.e. , @xmath95 for @xmath96 . the average degree @xmath97 can be controlled by @xmath91 for a specific value of @xmath92 . the dashed and dotted lines in fig .  [ fig2 ] represent the results for @xmath69 and @xmath98 , respectively . although these values of @xmath92 are not sufficiently small due to the restriction of computational resources to calculate @xmath99 iteratively and obtain @xmath34 in eq . ( [ eq : ntau ] ) , the maximum stable sizes for the sf networks are obviously greater than that for the errg . this implies that an sf network is more stable than the errg with the same average degree , which is consistent with the previous result showing the robustness of sf networks against cascading overload failures @xcite . next , we calculate @xmath0 for several combinations of the node tolerance parameter @xmath14 and the load carried by a single edge @xmath26 . the results depicted in fig . [ fig3 ] clearly indicate that @xmath0 increases with both @xmath14 and @xmath26 . it is obvious that the larger the node tolerance parameter @xmath14 , the more stable the network consisting of tolerant nodes becomes . the monotonous increase of @xmath0 with @xmath26 is due to the fact that @xmath100 given by eq . ( [ eq : ftau ] ) is a decreasing function of @xmath26 . it is interesting to notice that the maximum stable size @xmath0 shown in fig . 1(a ) is roughly proportional to @xmath82 . in fact , the ratio @xmath101 is about @xmath102 to @xmath103 independently of @xmath1 , unless @xmath1 is close to @xmath79 @xcite . this ratio is also insensitive to the parameters @xmath14 and @xmath26 , as well as the degree inhomogeneity . furthermore , as shown by fig . 1(b ) , the ratio @xmath67 slowly decreases with @xmath4 if @xmath104 but turns to a rapid decrease when @xmath67 becomes less than @xmath105 . these empirical facts give us significant information about the stability of the network . if the size of the largest component after cascading overload failures falls close to @xmath106 to @xmath107 of the size of the network before the cascade , the network is in immediate danger of a fatal breakdown . in order to accomplish further stable growth of the network , we need to raise the load reduction parameter @xmath1 . of course , the value of @xmath101 is peculiar to the present cascade model . however , it has been found that qualitative properties of our model are robust against changes in details of the model as long as failures are induced by temporally fluctuating loads @xcite . therefore , even for a real - world functional network , the ratio @xmath67 is supposed to decrease drastically with @xmath4 when this ratio falls below a certain value . our results suggest that we must take measures to prevent a fatal breakdown of a functional network if the decreasing rate of @xmath67 with increasing the network size becomes higher than its ordinary value . we have studied how large a functional network exposed to cascading overload failures can grow stably and evaluated the maximum stable size @xmath0 above which the network would face the crisis of a fatal breakdown . to this end , we employed the model of cascading overload failures triggered by fluctuating loads @xcite which is described by random walkers moving on the network @xcite . in this model , how quickly the total load is reduced during the cascade to prevent the fatal breakdown is quantified by the load reduction parameter @xmath1 . the maximum stable size @xmath0 was calculated by using the generating function technique and solving the master equation for the probability @xmath35 that a randomly chosen node has the degree @xmath8 at cascade step @xmath27 and the initial degree @xmath36 . our results show that @xmath0 is an increasing function of @xmath1 and diverges at a certain value @xmath79 . this implies that the faster the total load is reduced during the cascade , the larger the network can grow , and we can realize an arbitrarily large network if the total load is sufficiently quickly reduced ( @xmath2 ) . it has been also clarified that the degree inhomogeneity improves stability of the network . more precisely , for a given @xmath108 , @xmath0 for a scale - free network is larger than that for the erds - rnyi random graph with the same average degree . furthermore , from the empirical relation between @xmath0 and the network size @xmath82 giving @xmath0 , we argued how one detects and avoids the crisis of the network breakdown . the present results suggest that a certain relative size of the largest component size after cascading failures could be a signature for the impending network collapse . for further stable growth of the network , a more rapid reduction of the total load is required during a cascade of overload failures .    in this paper , we investigated only uncorrelated networks , while most of real - world functional networks have correlations between nearest neighbor degrees . for a correlated network , the probability @xmath109 of a node adjacent to a randomly chosen node of degree @xmath8 to experience an overload failure at cascade step @xmath27 depends on @xmath8 . this is in contrast to the case of uncorrelated networks , where @xmath109 given by eq . ( [ eq : phi2 ] ) is independent of @xmath8 . thus , the analysis becomes much more complicated for correlated networks than the present study . however , we suppose that our conclusion does not change qualitatively though @xmath0 depends on the strength of the degree correlation . a positive ( negative ) degree correlation simply makes a network robust ( fragile ) against various types of failures @xcite . thus , it seems plausible that also for our failure dynamics the degree correlation merely shifts the value of @xmath0 upward or downward . nevertheless , the relation between the stability of a functional network and its size must be strongly affected by the model of cascading failures . we hope that the problem of spontaneous instability in largely grown networks will be studied more extensively in diverse ways and models . 99 r.  albert and a .- barabsi , rev . phys . * 74 * , 47 ( 2002 ) . s.  boccaletti , v.  latora , y.  moreno , m.  chavez , and d .- u . hwang , phys . rep . * 424 * , 175 ( 2006 ) . s.  n.  dorogovtsev , a.  v.  goltsev , and j.  f.  f.  mendes , rev . phys . * 80 * , 1275 ( 2008 ) . d.  j.  watts and s.  h.  strogatz , nature * 393 * , 440 ( 1998 ) . barabsi , and r.  alberl , science * 286 * , 509 ( 1999 ) . m.  girvan and m.  e.  j.  newman , proc . usa * 99 * , 7821 ( 2002 ) . m.  e.  j.  newman , phys lett . * 89 * 026126 ( 2003 ) . m.  r.  gardner and w.  r.  ashby , nature * 228 * , 784 ( 1970 ) . r.  m.  may , nature * 238 * , 413 ( 1972 ) . r.  h.  porter , bell j. econ , * 14 * , 301 ( 1983 ) . a.  g.  haldane and r.  m.  may , nature * 469 * , 351 ( 2011 ) . r.  h.  heiberger , physica a * 393 * , 376 ( 2014 ) . j.  e.  meade , _ alternative systems of business organization and of workers renumeration _ ( routledge , london , 2010 ) . t.  shimada , sci . rep . , 4082 ( 2014 ) . a.  watanabe , s.  mizutaka , and k.  yakubo , j. phys . jpn . , * 84 * , 114003 ( 2015 ) . a.  e.  motter and y .- c . lai , phys . e * 66 * , 065102(r ) ( 2002 ) . p.  holme , phys . e * 66 * , 036119 ( 2002 ) . p.  crucitti , v.  latora , and m.  marchiori , phys . e * 69 * , 045104(r ) ( 2004 ) . s.  v.  buldyrev , r.  parshani , g.  paul , h.  e.  stanley , and s.  havlin , nature ( london ) * 464 * , 1025 ( 2010 ) . j.  gao , s.  v.  buldyrev , s.  havlin , and h.  e.  stanley , phys . lett . * 107 * 195701 ( 2011 ) . d.  zhou , h.  e.  stanley , g.   dagostino , and a.  scala , phys . e * 86 * , 066103 ( 2012 ) . s.  mizutaka and k.  yakubo , phys . e * 92 * , 012814 ( 2015 ) . m.  a.  de  menezes and a .- barabsi , phys . lett . * 92 * , 028701 ( 2004 ) . s.  meloni , j.  gmez - gardees , v.  latora , and y.  moreno , phys . lett . * 100 * , 208701 ( 2008 ) . v.  kishore , m.  s.  santhanam , and r.  e.  amritkar , phys . * 106 * , 188701 ( 2011 ) . j.  d.  noh and h.  rieger , phys . * 92 * , 118701 ( 2004 ) . m.  abramowitz and i.  a.  stegun , _ handbook of mathematical functions _ ( dover , new york , 1964 ) . m.  e.  j.  newman , s.  h.  strogatz , and d.  j.  watts , phys . e * 64 * , 026118 ( 2001 ) . d.  stauffer and a.  aharony , _ introduction to percolation theory _ ( taylor & francis , london , 1991 ) . if @xmath1 is close enough to @xmath79 , @xmath101 varies with @xmath82 according to @xmath110 . a.  v.  goltsev , s.  n.  dorogovtsev , and j.  f.  f.  mendes , phys . e * 78 * , 051105 ( 2008 ) . y.  shiraki and y.  kabashima , phys . e * 82 * , 036101 ( 2010 ) . m.  ostilli , a.  l.  ferreira , and j.  f.  f.  mendes , phys . e * 83 * , 061149 ( 2011 ) . t.  tanizawa , s.  havlin , and h.  e.  stanley , phys . e * 85 * , 046109 ( 2012 ) . s .- tan , j.  wu , m .- j . li , and x.  lu , epl * 114 * , 58002 ( 2016 ) .
Abstract:

Answer: we study how large functional networks can grow stably under possible cascading overload failures and evaluated the maximum stable network size above which even a small - scale failure would cause a fatal breakdown of the network . employing a model of cascading failures induced by temporally fluctuating loads , the maximum stable size @xmath0 has been calculated as a function of the load reduction parameter @xmath1 that characterizes how quickly the total load is reduced during the cascade . if we reduce the total load sufficiently fast ( @xmath2 ) , the network can grow infinitely . otherwise , @xmath0 is finite and increases with @xmath1 . for a fixed @xmath3 , @xmath0 for a scale - free network is larger than that for an exponential network with the same average degree . we also discuss how one detects and avoids the crisis of a fatal breakdown of the network from the relation between the sizes of the initial network and the largest component after an ordinarily occurring cascading failure .