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: recent experiments on giant plasma membrane vesicles ( gpmvs ) isolated from living cells and model membranes ( e.g. , artificially prepared lipid bilayers made of lipids and cholesterol ) display lateral inhomogeneities over length scales much larger than the typical sizes of the constituent lipid and protein molecules and stay in coexisting liquid - ordered ( @xmath0 ) and liquid - disordred ( @xmath1 ) phases @xcite . the universal scaling exponents that characterize the miscibility transition in gpmvs or model lipid bilayers are experimentally found to be close to those of the two - dimensional ( 2d ) ising model @xcite . ref . @xcite , by using lattice - based monte - carlo simulations , studies phase separations and critical fluctuations in two - component lipid membranes , and discusses the possibility of both quasi - abrupt transitions ( first order transitions ) and second order transitions , depending upon the details of lipid compositions . plasma membranes of living cells are however asymmetric @xcite . this asymmetry may affect the macroscopic properties at the critical point @xcite , although the knowledge about it is still at a preliminary level @xcite . hence , in the absence of detailed experimental results on the possible universal scaling properties of phase transitions in asymmetric membranes , studies of phase transitions in asymmetric membranes theoretically by using models with definite symmetries should be welcome . in particular , studies on how a generic curvature - composition nonlinearity , allowed by specific symmetry considerations , can affect the resulting phase behavior of the model , should be of interest .    in this article we ask whether asymmetry in an inhomogeneous fluid membrane can significantly enhance fluctuations near phase transitions . since any purported effects of asymmetry are likely to depend on how asymmetry affects the system free energy , we theoretically propose and study a model belonging to a particular symmetry for phase transitions in both symmetric and asymmetric inhomogeneous membranes . our model is a reduced model , in which the bilayer nature is not kept explicitly ; instead it is replaced effectively by an inhomogeneous single layer for simplicity . our main results are that depending upon the bending modulus @xmath2 the composition field that describes the inhomogeneity of the asymmetric membrane in our model may display ( i ) both first and second order phase transition through a critical point ( cp ) and ( ii ) a tricritical point ( tp ) . in addition , our model predicts _ generic nontrivial enhancement _ of fluctuations of asymmetric inhomogeneous membranes near tp and cp with a variance that depends on the system size @xmath3 in a power law fashion in two dimensions . further in our model , it is necessarily of second order belonging to the 2d ising universality class for symmetric membranes . we provide a mean - field ( mf ) analysis of the phase diagram of the model and calculate the universal scaling exponents that characterize the cp and tp in our model by using perturbative renormalization group ( rg ) methods . unlike more commonly used linear curvature - composition coupling @xcite , our model involves non - linear curvature - composition coupling belonging to a particular symmetry ( see below ) . thus it should be useful in understanding the general implications of non - linear curvature - composition relationships in experiments on membranes @xcite . validity of our model and the results that follow may be tested by measuring fluctuations in asymmetric inhomogeneous membranes and used to contrast with symmetric membranes . from perspectives of theoretical physics , our model is an example of a fluid membrane with internal degrees of freedom that displays rich critical behavior.the rest of the article is organized as follows : in sec . [ model ] we construct our model . we then discuss the mean - field analysis of our model in sec . [ mft ] , following by enumeration of the scaling exponents in sec . [ scale ] . we conclude and summarize our results in sec . [ conclu ] . a real asymmetric lipid bilayer generally consists of two lipid monolayers , which are constitutionally different ( i.e. , non - identical ) , e.g. , cell membranes of eukaryotic cells  @xcite . there are microscopic processes ( e.g. , flip - flop  @xcite ) which lead to exchange of the constituent molecules between the two monolayers of a bilayer . however , these exchange processes are typically very slow . in addition , there are biochemical processes which may favor maintaining compositional differences across the leaflets  @xcite . thus , the compositions of the two leaflets over relevant physical and biological time scales may be very different . coarse - grained models for an asymmetric lipid bilayer thus require two sets of scalar fields , one each for each of the monolayers , representing the compositional degree of freedom in the corresponding monolayer , see , e.g. , refs . @xcite for details . we , however , use a simpler approach , where a single field @xmath4 is used together with a height field to mathematically describe an asymmetric membrane . here , the asymmetry is no longer described by the composition differences between the two monolayers , but rather by the breakdown of the inversion symmetry of the membrane . thus , effectively , only a single layer is considered instead of a bilayer structure . there is already a substantial body of work that uses a single composition field degree of freedom , together with a height field , to model an asymmetric inhomogeneous membrane ; see , e.g. , refs . our model is in the same spirit as them . here are the details of our model : we use a coarse - grained description to build our model . for simplicity , we ignore the bilayer structure , and model it implicitly by a single inhomogeneous fluid membrane , microscopically made up of two components ( say , two types of lipids ) , @xmath5 and @xmath6 , with local concentrations @xmath7 and @xmath8 , with vanishing surface tension @xcite , described by a height field @xmath9 in the monge gauge @xcite and the concentration difference field @xmath10 $ ] , as the relevant compositional degree of freedom , since we are interested in phase transitions in the membrane @xcite , and @xmath11 . the detailed form of the free energy functional depends sensitively on the underlying symmetry of the model system . we choose @xmath4 to have the _ ising - symmetry _ , interacting via local ferromagnetic ising - like interactions  @xcite . further , the system being asymmetric is _ not _ invariant under the inversion symmetry @xmath12 . our imposition of these symmetries for an asymmetric membrane , dictates that the system and hence the corresponding free energy functional @xmath13 should be invariant under @xmath14 ( hereafter symi ) . thus in our model , one may have symmetry - allowed _ curvature - composition coupling _ of the form @xmath15 . in contrast , symmetric membrane models are typically invariant separately under @xmath16 and @xmath12 ( symii ) . therefore , taking everything into consideration @xmath17 is & & = d^dx [ r 2 ^ 2 + b 2()^2 + u 4!^4 + v 6!^6 + & + & 2(^2h)^2 + ^2 ^ 2h + _ 1 ^2 ( ^2h)^2 ] , [ free ] where @xmath18 with @xmath19 as the mf critical temperature for @xmath4 , and taking @xmath4 ( having the dimension of concentration ) to scale as @xmath20 , @xmath21 ( in 2d ) is a measure of energy - scale associated with composition fluctuations , where @xmath22 is the boltzmann constant and @xmath23 is ( microscopic ) correlation length of the order of a few angstroms . coupling constants @xmath24 and @xmath25 determine the strengths of lipid - lipid interactions ; @xmath26-term is added for thermodynamic stability ( see below ) . the @xmath27-term is related to a local _ spontaneous curvature _ @xmath28 @xcite . alternatively , it may be viewed as a _ local fluctuation _ in @xmath19 , such that local critical temperature @xmath29 . coupling constant @xmath27 embodies asymmetry in the model and implicitly contains information about geometry ( e.g. , shape , packing ) and asymmetric distribution of lipid molecules in the system . on physical grounds , following the arguments in ref . @xcite , we set @xmath30 , where @xmath31 is a mean concentration and @xmath32 is a typical ( molecular ) spontaneous curvature . this yields @xmath33 . since @xmath23 , the microscopic correlation length , should scale with the linear size of the constituent lipid molecules , @xmath27 may be experimentally varied by considering samples made of lipids of different molecular sizes . of course , for a pure fluid membrane @xmath34  @xcite . notice that our model is not the most generic in that a symi violating asymmetric linear curvature - composition coupling term may be present . effects of such a coupling , i.e. , @xmath35 being proportional to @xmath4 with the system invariant under @xmath36 ( hereafter symiii ) , have been extensively studied @xcite . our choice of curvature - composition coupling is made for the purpose of illustration  @xcite and complementary to that in refs . furthermore , the @xmath37-term may be interpreted as a contribution to @xmath38 or to a local ( effective ) bending modulus @xmath39 . here , we have omitted the geometric nonlinearities , which are subdominant to the existing nonlinear terms ( the @xmath40- and @xmath27-terms )  @xcite . note that the absolute sign of @xmath27 is arbitrary and may be switched by @xmath12 . the @xmath27-term violates the inversion symmetry @xmath12 ; where as for a symmetric membrane one has @xmath41 , leading to invariance under symii . lastly , since our model effectively assumes the bilayer to be a single incompressible sheet , it is unable to capture the effects of coupling of the local bending and densities of the two monolayers @xcite , which yields a new slow mode related to fluctuations in the density difference in the two monolayers . notice that in constructing our model we have been guided solely by general symmetry considerations , without any reference to microscopic or molecular level details . we believe it is useful to theoretically construct models with specific symmetry using this approach and study the ensuing macroscopic consequences , in view of the relative lack of quantitative results on the universal scaling properties of phase transitions in asymmetric membranes , in comparison with their symmetric counterparts . it is instructive to begin with a mf analysis for the phase transitions in the model , ignoring the @xmath37-term ( it is subdominant to the existing terms ; see below ) . the mf is constructed by ignoring all spatial correlations and in terms of constant values of the order parameter field @xmath42 and the curvature @xmath43 . here , @xmath44 and @xmath45 may both be zero or non - zero depending upon the phase concerned . the geometric nonlinearities are neglected . first consider a symmetric membrane , for which within mf , @xmath46 for all @xmath47 and @xmath48 for @xmath49 and @xmath50 for @xmath51 , respectively , thus yielding ising - like mf behavior with a second - order transition at @xmath19 . in contrast , for an asymmetric membrane at the mf level for @xmath49 one obtains , @xmath52 . however , for @xmath51 , one has @xmath53 , thus giving @xmath54 , implying a non - zero radius of curvature for an asymmetric membrane in the ordered phase . ignoring the spatial correlations of @xmath4 in the spirit of mf analyses , an effective landau free energy @xmath55 for a constant @xmath42 may now be constructed by replacing @xmath44 . it has the form @xmath56 where @xmath57 ( @xmath58 is a @xmath59 numerical coefficient ) can be @xmath60 positive , negative or zero . within the mean field picture , the effective landau free energy ( [ effre ] ) may be analysed in a standard way  @xcite : one finds the following possibilities : ( i ) a cp at @xmath61 with @xmath62 , ( ii ) a first order phase transition at @xmath63 for @xmath64 , and ( iii ) a tp at @xmath65 , i.e. , @xmath66 as schematically shown in fig .  [ first ] ; this phase diagram is identical to that observed in the normal - superfluid transition in liquid helium mixtures  @xcite . evidently , from its definition as given above @xmath62 or @xmath67 for large or small @xmath2 , respectively , implying lipid domain formations proceed through second or first order transitions for stiff ( large @xmath2 ) or soft ( small @xmath2 ) membranes , respectively . we now critical exponents @xmath68 and @xmath69 through the relations @xmath70 where @xmath71 and @xmath72 at both cp and tp . the presence of the first order transition in this model clearly underlines the requirement of a @xmath73 term in @xmath13 for thermodynamic stability . we now complement our mf analysis above by analyzing the properties of the fluid membrane in the spirit of mf . we start with the partition function corresponding to the free energy functional ( [ free ] ) : @xmath74.\label{part}\ ] ] by expanding perturbatively and integrating @xmath4  @xcite ( we ignore the @xmath40-term here ) , one may further obtain for effective bending rigidity @xmath75 to the lowest order ( setting @xmath76 ) , where @xmath77 is a wavevector . we shall see below that @xmath78 becomes scale - dependent , i.e. , it depends on the local scale @xmath79 . since @xmath80 is positive definite , @xmath81 generally , opening up the possibility of @xmath78 being zero or negative . since @xmath78 is generically smaller than @xmath2 , the asymmetric membrane will be generically softer or more flexible than a pure fluid membrane , a feature that should also be displayed by models belonging to symiii  @xcite , which use bilinear curvature - composition coupling . one may further define a persistence length @xmath82 , such that @xmath83 . noticing that for a system with linear size @xmath3 , the correction to @xmath78 scales as @xmath84 in 2d at the mf @xmath19 given by @xmath61 , we find @xmath85^{1/2}$ ] , thus @xmath86 for fixed @xmath2 and @xmath47 . physically , this essentially gives a linear scale over which the membrane remains flat . clearly , in our model , for an experimentally accessible asymmetric membrane with nearly flat conformation , ( bare ) @xmath27 should be sufficiently small , such that @xmath87 , linear dimension of the system . however , for quantitative predictions on @xmath82 , one requires to have a numerical estimate of @xmath27 , which at present is lacking . construction of atomistic models and numerical studies on them linking microscopic structures with macroscopic properties or parameters ( e.g. , @xmath88 ) would be useful  @xcite , which are beyond the scope of this work . for @xmath89 , our model shows instabilities akin to refs . notice however that the instabilities here for @xmath89 are due to a nonlinear curvature - composition coupling , unlike for models in refs . @xcite , where bilinear curvature - composition couplings are responsible for the instabilities . in the remaining part of the article below , we only consider @xmath90 . furthermore , considering that @xmath78 has got a @xmath91 ( perturbative ) correction to its bare value @xmath2 coming from the curvature - composition interaction , thence at 2d , with @xmath92 as the local normal to the membrane ( to the leading order in smallness ) , the variance of @xmath93 has a ( perturbatively obtained , ignoring logarithms ) contribution @xmath94 in addition to the contribution from ( bare ) @xmath2 . thus , @xmath95 formally diverges as @xmath96 . this is significantly stronger than the well - known @xmath97 dependence of @xmath98 that ensues for a pure symmetric fluid membranes  @xcite . even for an inhomogeneous symmetric membrane , @xmath99 must have a weaker @xmath3-dependence than @xmath84 , since the dominant nonlinearity responsible for the @xmath84-behavior ( i.e. , the @xmath27-term ) in our model is always absent for a symmetric membrane . this is thus a feature of an asymmetric membrane described by our model that may be directly testable in experiments and can be contrasted with the same from symmetric membranes . plane , depicting line of critical points , a tp and a first order transition.,height=226 ]    one may also present an equivalent phase diagram in the @xmath100 plane : see fig . [ second ] for a schematic phase diagram in the @xmath100 plane . plane , depicting lines denoting first ( broken line ) and second ( thick continuous line ) transitions and a tp.,height=226 ]    compare now with linear curvature - composition models @xcite having a curvature - composition coupling of the form @xmath101 ( with @xmath41 ) . it can be shown in that ( a ) @xmath102 does not contribute to the fluctuation corrections of @xmath4 ; with @xmath40 remaining as the most relevant nonlinearity , large - scale properties of @xmath4 is identical to that of the ising model , and ( b ) due to the absence of any nonlinear term there are no non - trivial fluctuation corrections to @xmath2 from @xmath102 , and hence , as long as @xmath103 , @xmath104 is identical to that of a pure fluid membrane . in contrast , our model for asymmetric membranes displays generic non - ising behavior of @xmath4 including a first order transition and scale - dependent diverging fluctuations of @xmath104 . furthermore , as already pointed out above , for a symmetric membrane with symii at 2d , @xmath105  @xcite , much weaker than its @xmath84-dependence for an asymmetric membrane in our model  @xcite . we confirm below the existence of cp and tp , extract the corresponding universal scaling exponents and show the divergence of the variance of asymmetric membrane fluctuations in a perturbative rg calculation  @xcite . for a symmetric membrane ( @xmath41 ) the remaining geometric non - linearities are irrelevant at 2d in the presence of the @xmath40-term with upper critical dimension @xmath106 , and hence are ignored . effectively , thus in our model there are no _ relevant _ coupling between @xmath4 and @xmath104 in a symmetric membrane . consequently the critical behavior of @xmath4 is identical to the ising model , thus belonging to the 2d ising universality class , already studied extensively in a rg framework @xcite : one obtains @xmath107 with @xmath108 ( the @xmath109 interaction is irrelevant in a rg sense ) . this is in agreement with experimental results on lipid bilayers by several groups  @xcite . further , fluctuation statistics of @xmath104 is identical to that of a pure fluid membrane . assuming a nearly flat conformation ( see discussions in the concluding section ) , predictions from our model is more complicated for an asymmetric membrane for which @xmath110 with @xmath106 , along with @xmath40 are relevant couplings ( in a rg sense ) . thus the ensuing universal scaling behavior at the physically relevant dimension @xmath111 is expected to be different from both the mean - field scaling behavior for an asymmetric membrane discussed before and 2d ising behavior of a symmetric membrane . in order to analyze this , we perform a systematic rg analysis of the model following standard procedure well - documented in the literature  @xcite . in particular , the renormalization @xmath112-factors for the different vertex functions , which are introduced to absorb ultra - violet ( uv ) divergences in the theory , are calculated by using a dimensional regularization together with a minimal subtraction scheme in terms of an arbitrary momentum scale @xmath113  @xcite . the vertex functions in the present model , which have primitive divergences and hence require renormalization , are given by ( the right hand sides of the expressions below are the bare values of the vertex functions in our model that can be easily read off from the free energy functional ( [ free ] ) ) . @xmath114 where @xmath115 is the vertex generating functional and is the legendre transform of the free energy corresponding to the partition function ( [ part ] ) ; @xmath116 are fourier wavevectors . the resulting rg flow equation yields non - trivial scaling of the correlation functions and thermodynamic quantities at the rg fixed points ( fp ) . the relevant feynman diagrams are shown in the appendix below . for calculational convenience we define an implicitly scale - dependent bare coupling constant @xmath117 by @xmath118 , such that the inhomogeneous part of the fluctuation corrections to @xmath40 is absorbed and define renormalized coupling @xmath119 . we similarly define renormalized @xmath120 . defining the rg @xmath121-function by @xmath122 , where @xmath123 , we find _ u & = & u_r[-+ 72u_r - 48_r^2 ] , [ betau ] + _ & = & _ r[-2 + 24u_r - 6_r^2 + 192_r^4 ] . [ betalambda ] the fps are determined by the zeros of the @xmath121-functions above : ( i ) gaussian ( or trivial ) fixed point given by @xmath124 , ( ii ) ising fp given by @xmath125 , ( iii ) tricritical fp given by @xmath126 and ( iv ) nontrivial fixed point ( nfp ) given by @xmath127 . very small value of @xmath128 at nfp suggests larger @xmath82 and hence easier experimental accessibility at nfp . each of these fps represents a different universality class , characterized by a different set of values of the critical exponents  @xcite . for example , we find for    * gaussian fp : @xmath129 , @xmath130 , * ising fp : identical to the symmetric membrane case , * tp : @xmath131 , @xmath132 , and * nfp : @xmath133 .    thus , @xmath134 at all the fps , in agreement with our qualitative discussions above on @xmath81 . from the scaling exponents calculated above , we find that    * @xmath135 at the tp , * @xmath136 at the ising fp , * @xmath137 is a spatial separation vector in @xmath138-dimension . thus , @xmath139 varies most slowly with @xmath140 at tp . vanishing @xmath68 at this order at tp is fortuitous and does not imply mf result ; @xmath69 picks up correction over its mf value already at this order . let us find out how the correlation of local normal @xmath141 to the membrane surface behaves . noting that @xmath142 , we find @xmath143 scales as @xmath144 at the nfp and @xmath145 at the tp . hence ,    * at tp , @xmath146 , * at nfp , @xmath147 .    compare this with a symmetric inhomogeneous membrane or a pure fluid membrane . disregarding the geometric nonlinearities ( which are irrelevant in a rg sense in the presence of @xmath40 and @xmath27 ) , for a symmetric membrane variance @xmath148 varies as @xmath97  @xcite . thus , @xmath95 depends on @xmath3 much more strongly than @xmath99 . for example , measuring lengths in the unit of molecular cut - off @xmath149 and taking @xmath150 ( typical size of an eukaryotic cell ) , we find @xmath151 at tp and 2 at nfp . finally , we consider linear stability of the different fps . we find ( i ) the gaussian fp @xmath152 is unstable along both @xmath153 and @xmath154 directions , ( ii ) the ising fp is unstable along the @xmath154 direction , but stable along the @xmath153-direction , ( iii ) the tp is unstable along the @xmath153 direction but stable along the @xmath154 direction , and finally ( iv ) the nfp is stable along both the directions . not surprisingly , the tp can be reached only by proper tuning of the parameters . the rg flow lines are schematically described in fig . [ flow ] .     plane . arrows indicate stable directions ( see text).,height=207 ] in this article , we have proposed and studied a model , belonging to a particular symmetry , for asymmetric inhomogeneous membranes . we have analyzed the nature of phase transitions and fluctuations in the model . notice that in our model , the local spontaneous curvature @xmath28 . hence , the average spontaneous curvature @xmath155 in mf , and is zero in the high temperature disordered phase ( @xmath156 ) and non - zero in the low temperature ordered phase ( @xmath51 ) . however , if we go beyond mf , then @xmath157 is non - zero , due to fluctuations , for both @xmath156 and @xmath51 . thus , when fluctuations are taken into account , the model continues to have a non - zero global spontaneous curvature @xmath158 , patrametrized by the ratio @xmath159 , at all @xmath47 . the value of @xmath158 generally depends on @xmath47 , and is different in the disordered and ordered phases , with a divergence ( for a formally infinitely large system ) at @xmath19 . thus , our model provides a description for a system with an overall curvature , which changes with @xmath47 . for the validity of the calculational scheme employed here , one should have @xmath160 . our model displays generic enhancement of fluctuations of asymmetric membranes near its second order phase transition , in contrast to symmetric membranes . for the latter , critical fluctuations of @xmath4 are independent of @xmath104 and belongs to the 2d ising universality class , a feature observed in experiments on lipid bilayers and gpmvs @xcite . our results for miscibility transitions in asymmetric inhomogeneous membranes are much more dramatic displaying generic non-2d ising - like behavior , with both second and first order transitions and a tp . further , asymmetric membranes are found to have dramatically enhanced @xmath3-dependent ( formally diverging in the thermodynamic limit ) fluctuations at tp and cp in 2d , in contrast to symmetric membranes . thus measurements of membranes fluctuations near cp or tp are important , and should yield signatures of asymmetry and composition - curvature coupling . first order transitions in our model is noteworthy , since a similar possibility has been discussed in ref . experiments on carefully prepared asymmetric membranes ( see , e.g. , ref . @xcite ) by standard , e.g. , fluorescence , methods ( see , e.g. , ref . @xcite ) should be useful , although more complicated curvature - composition coupling @xcite may be needed for quantitative reproduction of experimental data . at a technical level , for a generic nonlinear coupling of the form @xmath161 , where @xmath162 is an arbitrary function of @xmath4 , ( in our case @xmath163 ) , there is a general possibility of nontrivial renormalization of @xmath2 via nonlinear contributions to the self - energy @xmath164  @xcite at various orders of the underlying perturbation expansion . our results here is a specific example of such a possibility . as mentioned above , asymmetric inhomogeneous membrane models belonging to symiii also display reduction in @xmath78 due to composition fluctuations and thus are qualitatively similar to our results . in fact , far away from @xmath19 , models belonging to both symi and symiii yield qualitatively similar behavior for membrane fluctuations ( they should , however , predict different phase diagrams for @xmath4 ) . but their predictions for membrane fluctuations differ significantly as @xmath165 ( i.e. , near the critical zone ) , since ( considering @xmath166 ) , in the absence of any nonlinear term that involves @xmath104 in models belonging to symiii , there are no nontrivial ( scale - dependent ) fluctuation corrections to @xmath2 . hence , a linear curvature composition coupling model can not lead to any change in the scaling of the asymmetric membrane fluctuations ( although the magnitude can be enhanced ) . in contrast , due to the non - trivial renormalization of @xmath78 , asymmetric membrane fluctuations in our model displays enhanced fluctuations near @xmath19 , with _ non - trivial corrections _ to the scale dependences of the fluctuations . in a recent study on symmetric membranes by ayton _ @xcite , the curvature composition coupling is of the form @xmath167 , which breaks the ising symmetry . it can be shown , by integrating the composition field , that @xmath168 in the model of ref . @xcite , in the model leading to enhancement of membrane fluctuations , qualitatively similar to our results here . thus , our study may be viewed as an asymmetric membrane analog of ref . @xcite , where the ising symmetry is kept , but the inversion symmetry is broken .    direct comparisons of our results with experiments on inhomogeneous asymmetric bilayer is difficult , mainly due to the reduced nature of our model . nonetheless , whether or not mechanisms as discussed in our model is operative in a specific experimental set up , e.g. , our predictions on the connection of @xmath2 with the order of transitions and large fluctuations of asymmetric membranes near tp or cp , may be tested by measuring membrane fluctuations . complementing our coarse - grained modeling by numerical studies on microscopic models linking microscopic structures with macroscopic properties ( e.g. , @xmath169 ) would be useful in the present context  @xcite . since the signature and magnitude of the effective coupling @xmath60 depends directly on @xmath170 , our system may be tuned by controlling the membrane bending stiffness @xmath2 . alternatively , as discussed earlier , since @xmath171 ( taking @xmath172 , although @xmath173 generally ) , @xmath27 can be tuned by considering lipid membranes with lipid molecules having different linear sizes , and hence with different @xmath23 . considering the strong dependence of @xmath27 on @xmath23 , this should be a promising route . asymmetric membranes may be prepared and their phase behavior investigated by combining the langmuir - blodgett / schaefer method  @xcite with fluorescence - based imaging . membrane fluctuation measurements may be done by spectroscopic methods @xcite . bending modulus @xmath2 may be tuned experimentally , e.g. , by cholesterol or bar proteins  @xcite . accessing tp experimentally is expected to be considerably more difficult due to the additional tuning required . since @xmath174 generally , for a sufficiently large membrane a first order transition may take place , depending upon the value of @xmath2 . . since _ in - vivo _ membranes ( e.g. , red blood cell membranes ) tend to have small but finite shear modulii , it will be interesting to theoretically investigate effects of interactions between in - plane displacements and compositional degrees of freedom . since the divergence of asymmetric membrane fluctuations occurs exactly at cp or tp , it will also be interesting to study the possibility of budding transitions in our model and if the diverging fluctuations are signatures of the budding dynamics , as discussed in ref . lastly , as we mentioned above , a limitation of our model is the lack of bilayer structure in it . this , as elucidated in ref . @xcite , may affect the macroscopic behavior . this can be incorporated by generalizing our single compositional field description to two such fields . both of these fields should be mutually interacting , representing the interactions between the monolayers , in addition to their couplings with @xmath104 , which should be of the form allowed in symi . analyses of this generalized model , although more complicated , can be done in a straight forward way as here . qualitative features found in our model should be preserved , in addition to possible emergence of new features , e.g. , multicritical points and additional slow mode  @xcite etc . we hope our results will stimulate further theoretical and experimental works along these directions . ab thanks the max - planck - gesellschaft ( germany ) and the department of science and technology / indo - german science and technology centre ( india ) for partial financial support through the partner group programme ( 2009 ) . here we show all the relevant one- and two - loop feynmann diagrams for our model . the internal continuous line denote @xmath175 correlator . the internal broken line denote @xmath176 correlator . the external continuous line stands for @xmath4 and the external broken line stands for @xmath104 . 99 see for recent reviews : f. a. heberle and g. w. feigenson , _ cold spring harb perspect biol _ * 3 * , 004630 ( 2011 ) ; k. simons and j. l. sampaio , _ ibid . _ , p. 4697 ; e. l. elson _ et al _ , _ annu . biophys . _ * 39 * , 207 ( 2010 ) . t. baumgart _ et al _ , _ proc . ( usa ) _ * 104 * , 3165 ( 2007 ) . s. l. veatch and s. l. keller , _ biophys . j _ * 85 * , 3074 ( 2003 ) ; a. r. honerkamp - smith , s. l. veatch and s. l. keller , _ biochimica et biophysica acta _ * 1788 * , 53 ( 2009 ) . s. l. veatch _ et al _ , _ proc . a. _ * 104 * , 17650 ( 2007 ) ; s. l veatch , _ acs chem . biol . _ * 3 * , 287 ( 2008 ) . a. r. honerkamp - smith _ et al _ , _ biophys . j _ * 95 * , 238 ( 2008 ) ; a. r. honerkamp - smith , b. b. machta and s. l. keller , _ phys . . lett . _ * 108 * , 265702 ( 2012 ) . et al _ , _ new j. ph . _ * 13 * , 045019 ( 2011 ) . b. alberts , d. bray , j. lewis , m. raff , k. roberts , and j.d . watson , _ molecular biology of the cell _ garland , new york ( 1994 ) . allender and m. schick , _ biophys . j. _ * 91 * 2928 ( 2006 ) ; a. j. wagner _ et al _ , _ biophys . j. _ * 93 * , 4268 ( 2007 ) . since artificial lipid bilayers are generally symmetric and some asymmetry is believed to be lost during preparation of gpmvs , experimental results on them generally correspond to symmetric membranes ; see , e.g. , ref . @xcite above . it is still debatable at present whether plasma membranes of living cells undergo any phase transition . see for asymmetric membrane models with linear curvature - composition relation ( i.e. , invariance under symiii ) w. helfrich , _ z. naturforsch . _ * 28 * , 693 ( 1973 ) ; s leibler , _ j. physique _ * 47 * , 507 ( 1986 ) ; t. taniguchi , _ phys _ * 76 * , 4444 ( 1996 ) ; y. jiang _ et al _ , _ phys e _ * 61 * , r57 ( 2000 ) ; e. j. wallace et al , biophys j , 88 , 4072 ( 2005 ) ; j. gmez , f. sagus and r. reigada , _ phys . e _ * 80 * , 011920 ( 2009 ) ; t. baumgart et al , annu . chem , 62 , 483 ( 2011 ) ; c. zhu et al , biophys . j. , 102 , 1837 ( 2012 ) . in the limit of equal compositions ( which generally does not happen ) , the single field model used here may be retrieved from a bilayer model with two fields when the monolayers are _ locked _ , i.e. , in the limit of strong inter - monolayer interactions . this is experimentally realizable  @xcite . c. zhu , s. l. das and t. baumgart , _ biophys . j. _ * 102 * , 1837 ( 2012 ) . r. a. moss and s. bhattacharya , _ soc . _ * 117 * , 8688 ( 1995 ) . d. allender and m. schick , _ biophys . j. _ * 91 * , 2828 ( 2006 ) ; g. g. putzel and m. schick , _ biophys . j. _ * 94 * , 869 ( 2007 ) ; a. j. wagner , s. loew and s may _ biophys . j. _ * 93 * , 4268 ( 2007 ) ; y. hirose , s. komura and d. andelman , _ chemphyschem _ * 10 * , 2839 ( 2009 ) . the surface tension @xmath177 of an interface saturated by surfactant molecules is expected to vanish ( see , e.g. , s. a. safran , statistical thermodynamics of surfaces , interfaces and membranes ( addison- wesley , reading , mass . 1994 ) , and thus @xmath177 for a fluid membrane fluctuating freely in a solvent is usually assumed to be zero . it may be noted that biomembranes may experience non - zero @xmath178 , which however , becomes irrelevant on length scales @xmath179 . , edited by d. nelson , t. piran , and s. weinberg world scientific , singapore ( 1989 ) . see , e.g. , s. leibler and d. andelman , _ j. physique _ * 48 * , 2013 ( 1987 ) , for similar approaches . w. helfrich , _ j. phys . ( france ) _ * 46 * , 1263 ( 1985 ) ; h.j . deuling and w. helfrich , ibid . * 37 * , 1335 ( 1976 ) . m. haataja , _ phys . e _ * 80 * , 020902 ( r ) ( 2009 ) . see , e.g. , ref . @xcite and u. seifert , _ adv . phys . _ * 46 * , 13 ( 1997 ) for detailed discussions on fluid membranes . dropping the ising symmetry ( symi ) and not insisting upon any other symmetry ( e.g. , symiii ) allows us to have a @xmath180-term even for a symmetric membrane . since present available experimental results point towards 2d ising like behavior for symmetric membranes , these considerations as above have led us to use the ising symmetry to build our model that uses a nonlinear curvature - composition coupling . note that nonlinear curvature composition coupling of the form @xmath167 has been considered for symmetric membranes , see , e.g. , ref . this violates the ising symmetry . in contrast , we have retained the ising symmetry , but broken the inversion symmetry of the membrane . these geometric nonlinearities originate from ( i ) expansion of the area element with respect to the monge gauge base plane for small height fluctuations and ( ii ) nonlinear corrections to the form of the mean curvature in the monge gauge . there may be an additional symmetry - allowed inversion symmetry breaking geometric nonlinearity of the form @xmath181 , which may arise in an asymmetric membrane from a linear curvature term in @xmath13 : @xmath182 , since @xmath183 does not contribute to the free energy , it being a total derivative . while this term may be relevant ( same as the @xmath27- or @xmath40-terms ) we did not consider this term here for simplicity . u. seifert and s. a. langer , _ europhys . * 23 * , 71 ( 1993 ) ; u. seifert and s. a. langer , _ biophys . _ * 49 * , 13 ( 1994 ) e. h. graf _ et al _ , _ phys . _ * 19 * , 417 ( 1967 ) . p. m. chaikin and t. c. lubensky , _ principles of condensed matter physics _ , foundation books , new delhi ( 1998 ) . formally , integrating over @xmath4 gets problemtic near @xmath19 . this can be circumvent in a rg framework . i. r. cooke , k. kremer and m. deserno , _ phys . e _ * 72 * , 011506 ( 2005 ) . even if the geometric nonlinearity @xmath184 is considered , _ renormalized _ or effetive scale - dependent @xmath2 decreases as @xmath185  @xcite . thus @xmath186 continues to scale as @xmath97 . the other possible geometric nonlinearities are all subleading to the @xmath27- and @xmath40-nonlinearities considered here . therefore , although for a symmetric membrane ( @xmath41 ) all the geometric nonlearities are present ( having invariance under symii ) , their contributions to @xmath99 will not generate a @xmath84 dependence found in @xmath95 for an asymmetric membrane in our model . since in our model the @xmath4-field , being ising like , has a lower critical dimension @xmath187 and thus displays the standard order - disorder transition at 2d , with a diverging correlation length near @xmath19 , formal and rigorous enumeration of the scaling exponents requires the machinery of rg calculations . d. j. amit , _ field theory , the renormalisation group , and critical phenomena _ ( academic press , ny , 1978 ) . g. hooft and m. veltman , _ nuc . b _ , * 44 * , 189 ( 1972 ) ; c. de dominicis , e. brzin and j. zinn - justin , _ phys . b _ , * 12 * , 4945 ( 1975 ) ; e. brzin , j. c. l. guillon , and j. zinn - justin , _ phase transitions and critical phenomena _ , vol . 6 ( academic press , ny , 1976 ) ; j. zinn - justin , _ quantum field theory and critical phenomena _ ( oxford science publication , oxford 2010 ) . fluctuation corrections to coupling constants @xmath40 and @xmath27 are evaluated up to the one - loop order . if higher order loop corrections are included , then corrections to the exponents higher order in @xmath188 should arise . however , we do not forsee appearance of any new divergence that can not be absorbed by renormalizing the existing model parameters and fields . see appendix for the relevant feynmann diagrams . m. d. collins and s. l. keller , _ proc . usa _ * 105 * , 124 ( 2008 ) . g. s. ayton _ et al _ , _ biophys . j. _ * 88 * , 3855 ( 2005 ) . j. liu and j. c. conboy , _ langmuir _ * 21 * , 9091 ( 2005 ) ; j. yuan , c. hao , m. chen , p. berini and s. zou , _ langmuir _ * 29 * , 221 ( 2013 ) . t. betz and c. sykes , _ soft matter _ * 8 * , 5317 ( 2012 ) . see , e.g. , j. song and r. e. waugh , _ biophys . j _ * 64 * , 1967 ( 1993 ) ; l. r. arriaga _ et al _ , biophys . j. * 96 * , 3629 ( 2009 ) ; c. l. armstrong _ et al _ , eur . j , ( 2012 ) for studies on how cholesterol may be used to tune membrane bending rigidity . see also a. yaghmur _ et al _ , _ plos one _ , e479 ( 2007 ) for studies on tuning curvature by designer lipid - like peptide surfactants . additionally , bar proteins are known to rigidify membranes , see , e.g. , j. gilden and m. f. krummel , _ cytoskeleton ( hoboken ) _ * 67 * , 477 ( 2010 ) . p. b. s. kumar _ et al _ , _ phys . lett . _ * 86 * , 3911 ( 2001 ) .
Abstract:

Answer: we propose and analyze a model for phase transitions in an inhomogeneous fluid membrane , that couples local composition with curvature nonlinearly . for asymmetric membranes , our model shows generic non - ising behavior and the ensuing phase diagram displays either a first- or a second - order phase transition through a critical point ( cp ) or a tricritical point ( tp ) , depending upon the bending modulus . it predicts generic nontrivial enhancement in fluctuations of asymmetric membranes that scales with system size in a power law fashion at the cp and tp in two dimensions , not observed in symmetric membranes . it also yields two - dimensional ising universality class for symmetric membranes , in agreement with experimental results .


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 local polaron model describes a localized electronic state which is coupled to a boson field . the local state is then assumed to be hybridized with a continuum of delocalized states and are thus not an eigenstate of the electronic part of the hamiltonian . one of the first applications of the model was the coupling of plasmons to core holes@xcite and valence holes@xcite in metals . the boson field then represents the plasmons which can be excited by the introduction of a structureless core - hole or a valence - hole which may be hybridized with metallic states . plasmon excitation spectra are typically measured using electron energy loss spectroscopy ( eels ) where the energy loss of highly energetic electrons are measured after transmission through a metallic film . a similar application is that of deep - level spectroscopy@xcite where the boson field represents a phonon system which can be excited by the introduction of a core - hole . hybridization of the core - hole is then introduced to capture the degeneracy of the core - hole with a continuum of states with no core - hole but a high energy auger electron present . a somewhat different line of application is that of certain rare earth compounds which are known to give rise to mixed valence states@xcite . these states are characterized by an alternating valence in an otherwise periodic lattice which can result in unusual thermodynamic properties . the reason is that the difference in valence results in a difference in ionic radii and the extra valence electron thus have a strong coupling to the phonon system . the model designed to capture the effect consists of a localized @xmath0-state ( the extra valence electron ) coupled to a continuum of delocalized electrons and a phonon field coupled to the @xmath0-state . allthough , there is orders of magnitude differences between typical plasmon and phonon energies the physics in the models are very similar and only the model parameters differ .    finally the local polaron model has been applied to the problem of resonant tunneling@xcite in the context of electronic transport , and the very similar problem of hot electron femtochemistry at surfaces@xcite ( hefats ) . the idea of hefats is that an adsorbate system on a metal surface can have unoccupied electronic states which obtain a broadening due to interaction with the metallic states . if a hot electron ( an electron above the fermi level ) is generated in the metal it may interact with the unoccupied state and induce a chemical reaction on the surface . as an example we can think of a single molecule on a metal surface with one unoccupied electronic state well above the fermi level . a hot electron with an energy that matches the unoccupied orbital has the possibility of tunnelling from the metal to the molecule resulting in a transient occupation of the orbital . if the molecule was initially in an equilibrium position the electron will assert a force on the internal molecular degrees of freedom and can excite vibrational modes of the molecule before it tunnels back into the conductor . the molecule may acquire enough energy in this process to undergo a chemical reaction or a desorption event . a clever method to produce hot electrons is based on a metal - insulator - metal ( mim ) heterostructure as suggested by gadzuk@xcite . with an ideal mim device it is possible to tune hot electrons to any desired resonance of an adsorbate system and the approach thereby suggests the highly attractive possibility of performing selective chemistry at surfaces . such devices have been constructed and characterized@xcite and comprise a promising candidate for advanced hefats experiments . common to all these applications is that a linear bosonic coupling term has been assumed . it is by no means obvious that linear coupling captures the possibly complicated interaction of bosons and electrons , although it is probably often a good approximation . to examine the local polaron model beyond linear coupling we calculate the consequences of substituting the linear coupling term with a quadratic coupling term . in principle we should add the quadratic coupling on top of the linear , but this renders the model somewhat tedious to work with and the physics of quadratic coupling become hidden in complicated expressions . in contrast , having only quadratic coupling allow us to obtain inelastic scattering amplitudes very similar to those with linear coupling and the comparison is very instructive . in terms of bosonic potentials , a linear coupling term corresponds to a shift in the potential minimum whereas a quadratic coupling term corresponds to a shift in the frequency of the potential . the paper is organized as follows : in section [ models ] we present the local polaron model with a general coupling function and no bosonic dispersion . the electronic part is briefly reviewed and the wide band limit which is imposed in the remainder of the paper is defined . we then present the well known spectral function and inelastic scattering probabilities of the model with linear coupling and compare with a quadratically coupled model calculated in the present work . it is shown that for inelastic scattering involving a large number of bosons , quadratic coupling can give rise to much larger scattering probabilities . in section [ application ] we apply the theory to hot electron mediated excitation of the different modes of co adsorbed on cu(100 ) . the model parameters are calculated using density functional theory and the delta self - consistent field method and we find that linear coupling dominates desorption probabilities for the normal modes along the molecular axis , but vanishes for the frustrated rotations where quadratic coupling has to be taken into account . in appendix [ decoupling ] we derive a path integral representation of the newns - anderson retarded green function and show that the special properties of the wide band limit allow us to decouple bosonic and electronic degrees of freedom . appendix [ linear ] and appendix [ quadratic ] present the details of the calculations leading to the spectral functions and inelastic scattering probabilities associated with linear and quadratic coupling . in appendix [ combined ] , we show how a linear transformation of creation operators makes it possible to obtain the exact green functions including both linear and quadratic coupling . the general model we are concerned with is composed of a newns - anderson type hamiltonian@xcite coupled to a dispersionless ( single frequency ) boson field through a single electronic state . a dispersionless boson field naturally corresponds to a single mode of oscillation in an adsorbate system , whereas we can think of the dispersionless model as describing an einstein band if the boson field represents a phonon system . thus it is a model of non - interacting metallic electrons @xmath1 , a localized resonant state @xmath2 , and a harmonic oscillator described by the coordinate @xmath3 or equivalently the bosonic creation and annihilation operators @xmath4 and @xmath5 . the hamiltonian is given by @xmath6 where @xmath7 creates an electron in the state @xmath2 and @xmath8 creates an electron in the state @xmath1 . the function @xmath9 couples the resonant electron to the oscillator degrees of freedom . if one considers a hole coupled to the bosons instead of an electron , the order of @xmath7 and @xmath10 should be exchanged . it is natural to taylor expand the coupling function in the vicinity of the ground state minimum @xmath11 . including only the zeroth order term @xmath12 results in the newns - anderson model . since it is quadratic in the electronic creation and annihilation operators one could in principle formally diagonalize it . however , it is much more useful to investigate the resonant state @xmath2 ( which is not an eigenstate of the hamiltonian ) and we are thus led to consider the retarded green function : @xmath13 where @xmath14 is an electronic vacuum state and @xmath15 it is easily calculated in the energy domain using the dyson equation and the result is @xmath16 with @xmath17 and @xmath18 assuming the hopping matrix elements @xmath19 to be constant , @xmath20 becomes proportional to the metal density of states . if we furthermore assume the metal density of states to be wider than the resonance energy we can write @xmath21 and @xmath22 . this is the wide band limit which will be imposed in the present paper . it allow us to separate electronic and bosonic degrees of freedom in the general case and we can calculate green functions corresponding to linear and quadratic coupling exactly .    in the wide band limit the electronic retarded green function is @xmath23 and the spectral function is a lorentzian with full width at half maximum given by @xmath24 . when boson coupling terms are included ( first and second order taylor expansions of @xmath9 in ) the spectral function changes and inelastic scattering on the resonance becomes possible . suppose the resonance is initially unoccupied and the oscillator is in the state @xmath25 . the differential probability that an incoming particle ( hole or electron ) with energy @xmath26 will scatter through the resonance into a state of energy @xmath27 is given by the inelastic scattering matrix which can be expressed in terms of a two - particle green function@xcite as @xmath28 where @xmath29 the probability of transferring a given amount of energy to the bosons can thus be obtained by integrating the inelastic scattering matrix over the relevant values of @xmath27 . in this paper we focus on inelastic scattering by electrons , since this is the relevant quantity in the context of hefats and eels . however , the two - particle green function also appear in the calculation of optical transition amplitudes of an adsorbed molecule@xcite and knowing @xmath30 allows one to calculate a variety of observable quantities . finally , we note that the lifetime of the electron is independent of the boson coupling in the wide band limit . the probability that the state @xmath2 is unoccupied and that the oscillator is in any state at time @xmath31 given that the state was occupied and the oscillator was in the state @xmath32 at @xmath33 is @xmath34 which is proved in appendix [ decoupling ] . and @xmath35 and the coupling function @xmath36 . the vertical distance between the two potentials at the ground state minimum is @xmath37 . ] consider the state @xmath38 with the oscillator at @xmath3 and an electron occupying the resonance . the expectation value of the hamiltonian on such a state will depend on the value of @xmath3 due to the coupling @xmath9 and if we could calculate the electronic energy for all values of @xmath3 we would obtain an excited state potential @xmath39 . doing the same for the state with no electron in the resonance @xmath40 would result in a different potential @xmath35 and the coupling function should then be given by @xmath41 which is illustrated in figure [ fig : potential ] . in the model we have implicitly assumed that the potential @xmath35 is quadratic , but in general it could have any form . the potentials @xmath42 and @xmath35 are called born - oppenheimer surfaces and are obtained by moving the oscillator adiabatically in the electronic environment .      . the strength of the coupling is proportional to the derivative of the excited state at @xmath11 . ] we now taylor expand the coupling function @xmath9 to first order and express the boson coordinate in terms of creation and annihilation operators . this gives an interaction term : @xmath43 with @xmath44 where @xmath45 and @xmath46 are the mass and frequency of the oscillator . this model corresponds to the potentials and coupling function shown in figure [ fig : linear_pot ] .    as shown in appendix [ decoupling ] , the bosonic degrees of freedom decouple from electronic degrees of freedom in the wide band limit and the retarded green function thus becomes a product of an electronic part given by and an bosonic part . since the interacting term is linear in the oscillator coordinate the oscillator part of the hamiltonian can be diagonalized by  completing the square  or , equivalently , performing a canonical transformation which shifts the boson coordinate an amount proportional to @xmath47 : @xmath48 the retarded green function can then be calculated exactly for the @xmath25th excited state giving @xmath49,\end{aligned}\ ] ] where @xmath50 is the @xmath25th laguerre polynomial and @xmath51 . in this paper @xmath52 denotes a dimensionless effective coupling constant and @xmath53 denotes the _ exact _ green function corresponding to a coupling term @xmath54 and _ not _ the contribution from an @xmath25th order perturbative calculation as is sometimes custom . the spectral function is given by @xmath55 and for the ground state we obtain @xmath56 the spectral function is thus a sum of lorentzians of width @xmath24 and an internal spacing of @xmath46 and the amplitude of the @xmath45th peak follows a poisson distribution . it should be noted that the peaks do not represent excited states of the oscillator . it is the spectral function of the resonant electron with the oscillator in the ground state and the different peaks show that the coupling term mixes the eigenstates of the isolated oscillator . the real part of the self energy is always negative and given by @xmath57 and all physical observables are invariant to @xmath58 since linear coupling corresponds to a shifted harmonic oscillator and the direction of the shift is irrelevant . the two - particle green function and inelastic scattering matrix can also be calculated exactly@xcite and the probability that an incoming electron scatters on the resonance and excites the oscillator from the ground state to the @xmath25th excited state is @xmath59 with @xmath60 the probability of exciting the @xmath25th vibrational state thus essentially conserves the poisson distribution , but the lorentzians are replaced by the interference factor @xmath61 . the results and can also be obtained using a disentangling theorem@xcite as shown in appendix [ linear ] .      . ]    we will now consider an quadratic excited state potential energy surface @xmath42 which has a minimum that coincides with the ground state minimum , but has a different harmonic evolution . the potentials and coupling function corresponding to this is shown in figure [ fig : quadratic_pot ] . alternatively we could regard this model as a second order taylor expansion of the phonon coupling function @xmath9 when the first order contribution vanishes . the interaction term in the hamiltonian becomes @xmath62 with @xmath63 in this work we will only consider bound excited state potentials of the form @xmath64 and we can then write @xmath65 in the wide band limit the electronic and bosonic degrees of freedom decouple and the boson propagator can be evaluated using a generalization of the baker - campbell - hausdorff formula@xcite . as a result the retarded green function and spectral function can be calculated exactly . the derivation is shown in appendix [ quadratic ] and gives for the ground state : @xmath66 where @xmath67 this result is valid for @xmath68 which implies that @xmath69 . again , the spectral function is a sum of lorentzians , but with the @xmath45th peak damped by a factor of @xmath70 instead of a poisson distribution . the real part of the self energy is now given by half the frequency shift . the internal spacing between the peaks is @xmath71 and we see that the quadratic coupling only mixes the oscillator ground state with the even excited states of @xmath42 . this is due to the mirror symmetry of @xmath9 which implies that only oscillator states with equal parity mix . in the quadratic case , the effective dimensionless coupling @xmath72 is not simply given by @xmath73 as may have been anticipated , and thus the @xmath45th term in does not correspond to a @xmath45th order perturbative calculation of @xmath74 in @xmath75 . the calculation leading to the exact result is very different from the perturbative approach and we have checked that the second order taylor expansion of indeed gives the result obtained from second order perturbation theory . allthough the spectral function shows a series of peaks spaced by @xmath71 it is only possible to excite an integer number of @xmath46 through inelastic scattering . the reason is of course that the boson field is completely decoupled from the electronic states in the asymptotics of a scattering event and will thus be observed in an free oscillator eigenstate . again , the symmetry of the quadratic coupling means that transitions involving an uneven number of bosons are forbidden . the two - particle green function and inelastic scattering matrix are calculated in appendix [ quadratic ] and the probability for an incoming hot electron to excite @xmath76 quanta of oscillation when initially in the ground state is : @xmath77 with @xmath78 the structure is very similar to the case of linear coupling . with linear coupling the probability for an electron to create @xmath25 bosons are proportional to the @xmath25th order taylor expansion of @xmath79 and normalized by @xmath80 whereas in the quadratic case the probability to create @xmath76 bosons are proportional to the @xmath25th order taylor expansion of @xmath81 and normalized by @xmath82 . in the context of eels and plasmon excitations , one would now observe a series of peaks spaced by @xmath83 . if the plasma frequency is not known the spacing itself can not give clues to whether linear or quadratic coupling governs the transitions , but one could use the relative amplitude between peaks since these follow a poisson distribution if linear coupling dominates and the distribution @xmath84 if quadratic coupling dominates . if both linear and quadratic coupling is present one would observe a coupling dependent enhancement of every second peak . in a model with linear coupling , the probability of exciting @xmath76 vibrational quanta is proportional to @xmath85 whereas it is proportional to @xmath86 in a quadratic coupled model . this implies that if @xmath87 a quadratic coupling term will give rise to larger inelastic scattering probabilities than a linear term . even with @xmath88 a quadratic coupling term may have stronger effect for large @xmath25 since the expansion coefficients of @xmath89 decay slower than those of @xmath90 . this is illustrated in figure [ fig : prob ] , where the probability of transferring @xmath25 vibrational quanta to the ground state is shown for linear and quadratic coupling .     through resonant inelastic scattering with linear and quadratic coupling . the parameters are @xmath91 , @xmath92 , and @xmath93 @xmath94 . even though @xmath88 the quadratic coupling becomes dominating for large @xmath25 due to the slowly decaying expansion coefficients . one should also note the spacing between peaks which is @xmath46 for linear coupling and @xmath71 for quadratic coupling . the centers of the probability distributions are approximately shifted by @xmath95 for the linear coupling and @xmath96 for the quadratic coupling relative to the bare resonance energy @xmath37 , since this is where the binomial coefficients in and have their maxima . ] as an example of a system where the dynamics can be approximated by a local polaron model , we consider the problem of hot electron mediated energy transfer on a metal surface . such an energy transfer can lead to desorption of adsorbed molecules @xcite or induce chemical reactions which can not proceed by thermal heating.@xcite . the conceptual picture of the process is the following : hot electrons are generated in the metal by means of an mim device or a femtosecond laser . the hot electrons may then interact with a chemisorped molecule by tunneling from the metal to an unoccupied molecular state and excite vibrational states in the molecule . if enough energy is transferred to the molecule either by a single or multiple scattering events , the molecule may eventually desorp or break an internal chemical bond . as a particular example we will calculate transition probabilities for co adsorbed on cu(100 ) .    to calculate inelastic scattering probabilities within the local polaron model we need to obtain the coupling function @xmath9 . as described in section [ adiabatic ] , we can fix the molecule at different positions and calculate the potential energy surfaces @xmath42 and @xmath35 at each point and @xmath41 . the model does not directly contain coulomb interactions between electrons , but these are included in the calculation of @xmath9 which thus becomes an effective coupling that is supposed to contain all the electronic interaction associated with the excited state of the molecule . the potential energies @xmath42 and @xmath35 has been obtained using the code ` gpaw`@xcite which is a real - space density functional theory ( dft ) code that uses the projector augmented wave method.@xcite in all our calculations we used the revised perdew - burke - ernzerhof ( rpbe ) exchange - correlation functional @xcite since this has been designed to perform well for molecules adsorbed on surfaces , and has been shown to perform better than the original pbe functional@xcite for adsorbed molecules . we set up a cu(100 ) surface consisting of three atomic layers with the top layer being relaxed . 10   of vacuum has then been introduced above the slab and 0.50 monolayer of adsorbate molecules relaxed at top sites which is the preferred adsorption site . both molecules adsorb with their molecular axis perpendicular to the surface with o pointing away from the surface . we then did a normal mode analysis and mapped out the three ground state potential energy functions @xmath97 corresponding to the two normal modes that involve the perpendicular degrees of freedom and a frustrated rotation . the perpendicular modes roughly correspond to an internal stretch @xmath98 and center of mass @xmath99 . we do not include the three remaining molecular modes since one is another frustrated rotation with identical properties to the one considered , and the two frustrated translations are only weakly coupled to the resonant electron and are not expected to play a significant role in the femtochemistry . in all calculations we use a p(2x2 ) cell , sample 12 irreducible k - points in the surface plane , and use a grid spacing of 0.2  .     is a generalized coordinate representing the deviation from equilibrium . @xmath100 corresponds to a @xmath101 angular deviation from the perpendicular position . ]    to find the excited state potential energies @xmath102 corresponding to the three normal modes of interest , we have used the method of linear expansion @xmath103scf which has been published in a previous work@xcite and implemented in ` gpaw ` . in the previous publication we have tested the method against inverse photo - emission spectroscopy , and found that it performed well for molecules chemisorped on surfaces.@xcite in each step of the self consistency cycle an electron is removed from the fermi level , the density of an excited state is added to the total density , and the band energy of this state is added to the total energy . to get the band energy right we need to expand the excited state on the kohn - sham orbitals found in each iteration . the method is thus a generalization of the usual @xmath103scf where occupation numbers are changed . instead of changing occupation numbers we occupy an orbital which is not an eigenstate of the kohn - sham hamiltonian but a superposition of eigenstates in such a way that the state is as close as possible to the original molecular state . in the present case the excited state is the anti - bonding @xmath104 orbital of co. in figure [ fig : pes ] we show the ground and excited state potential energy surfaces corresponding to the frustrated rotation along with @xmath9 . it is clear that the excited state potential is not exactly a quadratic potential and the parameter @xmath75 which we need to calculate transition probabilities will depend on how we fit this potential to a quadratic form . however the width of the gaussian ground state vibrational wavefunction corresponds to @xmath105  and for low lying excitations we can thus use this region of the potential which is rather flat . in fact , a closer look at the excited state potential reveals that the ground state minimum geometry actually has an unstable extremum in the excited state , but since the curvature is rather small we will simply approximate it by a constant potential . for both perpendicular modes we find that @xmath106 and quadratic coupling can thus be neglected . in contrast , due to symmetry the excited state potential energy of frustrated rotation is invariant to @xmath107 and the linear coupling term thus vanishes . we have calculated the excitation energy to @xmath108 and the resonance width is estimated from the kohn - sham projected density of states to @xmath109 . in table [ tab ] , we display the calculated parameters corresponding to the three modes . .parameters for co adsorbed on cu(100 ) . all number are @xmath110 . note that while the quadratic coupling for the two perpendicular modes are very small and thus neglectable , the linear coupling of frustrated rotation vanishes exactly due to symmetry . [ cols="^,^,^,^",options="header " , ]     we note that when calculating transition probabilities we should include all modes in the model , because even if the modes are not coupled directly they have an indirect coupling since they all interact with the resonance . it is possible to obtain expressions for the scattering matrix including more than one mode , but these are rather complicated to handle and for weakly coupled systems the physics can usually be extracted from three one - mode models.@xcite in figure [ fig : co_prob ] we show the calculated probabilities for a hot electron to excite the different modes of co adsorbed on cu(100 ) . the internal stretch and and frustrated rotation show transition probabilities on the same order of magnitude whereas the center of mass vibrations are very unlikely to get excited . this is in accord with calculations of the electronic friction coefficients of this system@xcite which is very closely related to the coupling function @xmath9.@xcite        the frequency of internal vibration is five times larger than both the center of mass and frustrated rotation frequencies and as previously shown@xcite the stretch mode will completely dominate the total energy transfer . thus , in a simple model where hot electron mediated desorption@xcite is reduced to calculating the probability of transferring the chemisorption energy to the adsorbate , the internal mode governs the desorption probability . nevertheless , our estimate of @xmath24 is based on the kohn - sham density of states which may give a poor description of the electronic spectral function @xmath111 . if @xmath24 is significantly smaller than our estimate , the quadratically coupled frustrated rotation will play an important role in hot electron mediated desorption for this system . we have calculated the spectral function and inelastic scattering amplitudes in a local polaron model with quadratic coupling to bosons . the probability of exciting @xmath25 bosons is found to be damped by a distribution function given by the @xmath25th taylor expansion of @xmath112 which decays much slower than the poisson distribution appearing in a linearly coupled model . hence for comparable values of linear and quadratic coupling constants , a quadratic term will dominate inelastic scattering probabilities involving a large number of bosonic excitations .    as an application we have considered the problem of hot electron mediated vibrational excitations of molecules adsorbed on metal surfaces . the coupling constants were calculated from the excitation energy along the molecular normal modes using delta self - consistent field dft . it was found that quadratic coupling is important for exciting the frustrated rotations since this mode does not couple linearly due to symmetry . a major approximation in the model is the quadratic assumption for the ground state potential . in our numerical example with hefats it is clear from figure [ fig : pes ] that the potentials is not exactly quadratic . for the center of mass mode the situation is even worse and a morse potential is much better suited to describe this mode . the anharmonic deviations are likely to have a significant effect on high lying excited states , but renders the model much more complicated . in fact , since the coupling to the internal stretch mode seems to govern the rate of energy transfer,@xcite one has to assume that the energy is readily redistributed to other degrees of freedom and anharmonic coupling is thus expected to play a vital role in the actual desorption process . the wide band limit has been essential in the derivation of scattering amplitudes , and we do not have the means to solve the model ( with linear or quadratic coupling ) exactly beyond this approximation . however , it would be very interesting to do perturbation theory with the general retarded green function to examine the effect of energy dependence in the electronic self energy . we are grateful to karsten jacobsen for his advice and help with the path integral representation of the newns - anderson model and to jakob schitz for helpful comments and guidance into the exciting field of hot electron femtochemistry . the center for individual nanoparticle functionality ( cinf ) is sponsored by the danish national research foundation . this work was supported by the danish center for scientific computing . the path integral representation of propagators often give a renewed insight into the underlying physics , allthough the mathematical complexity can be somewhat larger . by writing the newns - anderson retarded green function as a sum over paths , we see that each path can be understood as sequence of jumps from the resonance to the metallic band and we then have to sum over all possible time intervals between each jump . in the wide band limit the time spend in the metal band goes to zero and the electron thus spends all the time of propagation in the resonant state . the newns - anderson model is given by with @xmath113 : @xmath114 and @xmath115 the path integral representation is derived by dividing the time interval @xmath31 in @xmath116 intervals of length @xmath117 . when @xmath116 becomes sufficiently large we can take @xmath118 . we then insert @xmath119 complete sets of states @xmath32 such that the green function becomes a sum over @xmath116-fold products of matrix elements @xmath120 assuming that @xmath121 the states @xmath32 can either be @xmath2 or @xmath1 and the matrix elements @xmath122 , @xmath123 and @xmath124 represent propagation in the resonance , propagation in the band and a jump from band to the resonance respectively . when we take the limit @xmath125 , becomes formally exact and the jumps between band and resonance become instantaneous . it is then most convenient to order the terms in according to the number of jumps . since the endpoints of the time interval is at the resonance , a jump into the band has to be accompanied by a jump back into the resonance and each such `` band excursion '' comes with a factor of @xmath126 where @xmath127 is the time spend in the @xmath128th excursion into band . it is also clear that @xmath129 excursions into the resonance has to be accompanied by @xmath130 resonant propagation factors @xmath131 where @xmath132 is the @xmath128th time interval in the resonant state . finally , for a given number of band excursion we have to integrate over all possible band and resonance time intervals and the retarded green function becomes @xmath133 where the delta function has been introduced to ensure that the time intervals sum to @xmath31 and the theta function has become redundant . we can use the delta function to eliminate the @xmath132 integration variables and get @xmath134 with @xmath135 by using that @xmath136 it is now possible to evaluate and recover the result . in the wide band limit , @xmath137 , which implies that the electron does not spend any time in the band and the retarded green function becomes a sum over paths which are composed of _ instantaneous _ excursions into the band . we use the notation @xmath138 to represent a position in state space and @xmath139 as a formal expression representing the sum over all paths weighted by the newns - anderson action @xmath140 . we now proceed with the full hamiltonian . introducing the bosonic coordinate @xmath3 the full retarded green function ( with the boson field in the state @xmath11 ) can be written @xmath141 where @xmath142 is the newns - anderson action , @xmath143 is the free bosonic action corresponding to the hamiltonian @xmath144 , and @xmath145 is the interaction part of the action corresponding to the hamiltonian @xmath146 . however , there is a much nicer way to handle the coupling to the boson field . one can think of the bosons as an environment influencing the paths of the resonant electron and it can be shown that the green function can be written@xcite @xmath147 with the environment time evolution operator @xmath148 where @xmath149 is the environment hamiltonian evaluated on an electronic state fixed at @xmath138 and @xmath150 denotes time ordering . thus , the price we pay in separating bosonic and electronic degrees of freedom is an explicit path dependence in the environment part of the propagator . in general it is not possible to evaluate @xmath151 on all possible paths , but in the wide band limit it is particularly simple . the reason is that the resonant electron stays on the resonance in all possible paths and the environment hamiltonian is therefore _ independent _ of the electronic path . in fact , the environment propagator becomes @xmath152 and the electronic and bosonic degrees of freedom completely decouple in the retarded green function : @xmath153 the situation is very similar for the two - particle green function : @xmath154 the resonant electron is first propagated forward in time from @xmath155 to @xmath31 and then backward in time from @xmath156 to @xmath157 . the interaction vanishes between the @xmath158 and @xmath159 because the resonant state is unoccupied here . again , the full green function can be written in terms of a bosonic influence propagator and in the wide band limit the electronic and bosonic degrees of freedom decouple so : @xmath160 with @xmath161 this can be seen by applying the arguments above to both of the time evolution operators and the fact that @xmath162 . we observe that for any coupling function @xmath9 the above form of the two - particle green function implies that @xmath163 . this means that in the wide band limit , the resonant lifetime is unaffected by the phonon coupling since the probability of finding the electron in the state @xmath2 at time @xmath31 ( and the oscillator in any state @xmath164 ) given that it was there at @xmath33 ( where the oscillator was in the state @xmath32 ) is @xmath165 so , in the wide band limit the resonant state always has a well defined lifetime given by @xmath166    the problem of calculating the inelastic scattering matrix has now been reduced to evaluating the phonon propagator @xmath167 . in general this is not an easy task , but we will show that in the case of linear and quadratic coupling terms we can use a disentangling theorem@xcite to write the exponential operators in a form that allows a direct evaluation of the expectation value . the theorem is a generalization of the baker - campbell - hausdorff theorem and states that if @xmath168 , @xmath169 , and @xmath170 are three operators with a closed commutator algebra then @xmath171 , where @xmath172 , @xmath173 , @xmath174 , and @xmath175 are known functions of @xmath176 , and the commutation parameters . the model with linear coupling function given by @xmath177 is well - known and the one - particle green functions can be obtained exactly in the wide band limit by a canonical transformation.@xcite here we derive it using the formalism above and the disentangling theorem.@xcite to obtain the one - particle green function we need to evaluate the boson propagator @xmath178 using the disentangling theorem on the second exponential operator leads directly to the expression @xmath179,\qquad g_1=\big(\frac{\lambda_1}{\omega_0}\big)^2\end{aligned}\ ] ] where @xmath180 is the @xmath25th laguerre polynomial . to obtain the ground state spectral function we taylor expand @xmath181 , perform a fourier transformation , and take the imaginary part . the procedure can also be used to obtain the two - particle green function . the object of interest is now the two - particle boson propagator which we write @xmath182 where we used two different forms of the disentangling theorem to move @xmath183 and @xmath184 to the left and right respectively . we then use the theorem again to move @xmath185 to the left and obtain @xmath186 finally , we can use the baker - campbell - hausdorff theorem to collect all lowering operators at the the right . the evaluation of the remainder gives a laguerre polynomial and the result is @xmath187,\end{aligned}\ ] ] with @xmath188 the inelastic scattering matrix has previously been derived@xcite and we will not repeat the calculation here . we now consider a quadratic coupling function of the form @xmath189 to obtain the retarded green function we would like to calculate the boson propagator @xmath190 we proceed by disentangling the exponential operator : @xmath191 where @xmath192 and @xmath193 . this is valid for a bound excited state potential with positive second derivative in which case the argument of the square root is positive . in the case of an unbound excited state potential , the functions @xmath0 and @xmath194 involve real hyperbolic functions and the spectral function acquires a qualitatively different structure . acting with the operator @xmath195 on a state @xmath32 gives @xmath196}\frac{\alpha^l}{l!}\big(\frac{n!}{(n-2l)!}\big)^\frac{1}{2}|n-2l\rangle,\end{aligned}\ ] ] where @xmath197 $ ] means the integer part of @xmath198 . collecting it all and noting that @xmath199 , gives the retarded green function in the wide band limit : @xmath200}\frac{h^{2l}}{(l!)^2}\frac{n!}{(n-2l)!}\end{aligned}\ ] ] with @xmath201 to find the spectral function of the oscillator ground state we taylor expand the square root and obtain @xmath202 with @xmath203 fourier transforming and taking the imaginary part then gives @xmath204 where we also used that @xmath205 . note that the condition of @xmath68 implies that @xmath69 . we now need the propagator @xmath206 we restrict the calculation to the ground state two - particle green function which using the disentangled expression becomes @xmath207 it is also possible to obtain a closed expression that does not involve the infinite sum , since instead of inserting a complete set in we could have brought all lowering operators to the left by repeated use of the disentangling theorem . however , to calculate the inelastic scattering matrix we need to integrate over @xmath156 which is more convenient in the present form . performing the @xmath156 integral and using the resulting delta function to replace @xmath208 with @xmath209 , leaves the two remaining integrals as complex conjugates . we note that @xmath210 and write @xmath211 with @xmath212 the reason we have excluded the @xmath213 term is that it does not give rise to inelastic scattering and the elastic part of the scattering matrix have an additional term that we do not consider here.@xcite this expression implies that quadratic coupling can only give rise to inelastic scattering events involving an even number of vibrational quanta . this is also true if the initial state is not the ground state , since from we see that in general @xmath214 has to be even . to evaluate the inelastic scattering matrix we note that @xmath215 where we used that @xmath216 . the taylor expansion of the square root gives @xmath217 leading to @xmath218 and @xmath219 it is in principle straightforward to generalize the expressions above to the case of a linear _ and _ a quadratic coupling term in the hamiltonian . the linear term can be transformed away by noting that @xmath220 with @xmath221=1.\end{aligned}\ ] ] since the commutator algebra of @xmath222 and @xmath223 is identical to that of @xmath5 and @xmath4 we can immediately write down the one - particle boson propagator in its disentangled form @xmath224 using that @xmath225 we can evaluate the propagator in the ground state and obtain @xmath226 the expression clearly reduces to @xmath227 and @xmath228 in the limits @xmath229 and @xmath230 respectively . it should be straightforward to obtain the spectral function by fourier transforming this expression after a taylor expansion of the exponentials . however , the result becomes rather involved and we will not attempt to do the calculation here .
Abstract:

Answer: we calculate the inelastic scattering probabilities in the wide band limit of a local polaron model with quadratic coupling to bosons . the central object is a two - particle green function which is calculated exactly using a purely algebraic approach . compared with the usual linear interaction term a quadratic interaction term gives higher probabilities for inelastic scattering involving a large number of bosons . as an application we consider the problem hot electron mediated energy transfer at surfaces and use the delta self - consistent field extension of density functional theory to calculate and compare coupling parameters and probabilities for exciting different vibrational modes of co adsorbed on a cu(100 ) surface