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: classical continuum theories lack any length scales and as such provide leading order approximations for a number of problems that contain microstructures . microstructures typically introduce characteristic length scales that may induce gradient dependences of various kinds . several continuum theories have been developed to deal with microstructures and their attendant phenomena . these include micropolar , micromorphic , strain - gradient , non - local , etc . see eringen @xcite and nowacki @xcite for catalogs of such phenomenological theories . physical theories for complex fluids such as liquid crystals have also been introduced , based on symmetry breaking phase transitions that yield statistically - defined order parameters as additional thermodynamic variables . see de gennes and prost @xcite for discussions of the fundamental principles of order parameter physics for liquid crystals . see holm @xcite for a variational description of order parameter theories of complex fluids . often these phenomenological theories are combined with geometrical discussions based on the theory of cosserat and cosserat @xcite . the mathematical theory of continuum mechanics for complex , or composite materials produced a number of interesting phenomenological models in the 1960 s . among these are the models of differential type introduced by rivlin and ericksen @xcite and the multipolar models of green and rivlin @xcite . the history of how these models were tested in comparison with experiments and refined by _ ab initio _ assumptions of thermodynamics is recounted , for example , in dunn and fosdick @xcite , eringen @xcite and fosdick and rajagopal @xcite . recently , bellout @xcite considered a fusion of the models due to rivlin and ericksen @xcite and those of green and rivlin @xcite . the resulting theory introduced higher order spatial velocity gradients into the energy that regularized the model solutions and endowed the model with promising stability characteristics . the present work specializes to a subclass of the rivlin - ericksen - green multipolar fluids treated in bellout @xcite that has energy density given by @xmath2 where @xmath3 is the mass density and @xmath4 is the strain rate tensor . materials whose energy density takes this form are called _ gradient fluids of degree @xmath1 _ , where @xmath1 is the order of the velocity gradients appearing in ( [ grad - n - erg ] ) . the case @xmath5 is the euler fluid ( no velocity gradient dependence ) , while the case @xmath6 coincides with the @xmath7 grade fluid @xcite , whose energy depends on the velocity gradient through the strain - rate , @xmath8 . the equations of motion for gradient fluids are obtained from the eulerian form of hamilton s principle introduced in holm , marsden and ratiu @xcite called the euler - poincar theory for continua with advected quantities . in vector notation , this is @xmath9 for the class of lagrangians we shall consider , for @xmath10 , @xmath11 and @xmath7 gradient fluids , one has @xmath12 the term in @xmath13 in this lagrangian boosts the gradient fluid flow into a frame rotating with angular frequency @xmath14 , while the term in @xmath15 imposes the constraint @xmath16 . hence @xmath17 , as implied by substituting @xmath16 into the continuity equation , @xmath18 many mathematical regularity properties are available for the class of gradient fluids , especially for the case that the lagrangian @xmath19 in ( [ lag - grade-2 - 3 ] ) provides a norm ( when evaluated on the constraint surface , @xmath16 ) . however , these regularity properties for gradient fluids will be discussed elsewhere , following foias @xcite . the objective of the current paper is to investigate the stability properties of cc solutions of the gradient fluid equations . for cc solutions in an unbounded domain , the fluid velocity is linear in the spatial coordinate and the pressure is quadratic . the cc solutions may be regarded as the first term in a taylor expansion in space , around a stagnation point of the gradient fluid flow in a moving frame . we shall use the theory of elliptic instability to investigate the exact nonlinear growth rates when cc solutions interact with a wave packet whose phase is frozen into the cc flow for gradient fluids of degree @xmath20 . ( the euler case @xmath5 was studied in the original work of craik and criminale @xcite . see also craik @xcite , miyazaki @xcite and kerswell @xcite for subsequent developments . see also lagnado and simmen @xcite and goddard and alam @xcite for similar analyses for an upper - convected maxwell fluid and granular media , respectively . ) [ [ outline . ] ] outline . + + + + + + + +    section [ ep - form - sec ] summarizes the properties of ideal gradient fluids that follow directly from their euler - poincar formulation . these properties include energy conservation , momentum balance and kelvin circulation preservation , all of which follow from noether s theorem . we then specialize to gradient fluids of degree @xmath20 . section [ cc - soln - sec ] introduces the cc solutions for @xmath11 and @xmath7 gradient fluids . section [ elliptic - instab - sec ] discusses their stability properties for both inviscid and viscous cc solutions . here we introduce viscosity as in the theory of @xmath7 grade fluids , to which the gradient fluids reduce when @xmath6 . section [ conclusion - sec ] summarizes our conclusions . [ [ hamiltons - principle - for - first - and - second - gradient - fluids . ] ] hamilton s principle for first and second gradient fluids . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the mathematical basis common to all ideal fluid motions is hamilton s principle @xmath21 where @xmath22 is the lagrangian for the system . we work in the eulerian representation of fluids , where the euler - lagrange equation is replaced by the euler - poincar equation . see ref . @xcite for a detailed discussion of euler - poincar theory . this paper focuses on the incompressible motion of first and second gradient fluids in a rotating frame . thus , the class of lagrangians we shall consider has the form @xcite @xmath23 where @xmath24 is the vector potential for the coriolis parameter , i.e. , @xmath25 . specifically , we shall take , @xmath26 here @xmath15 is pressure ( a lagrange multiplier ) , @xmath3 is mass density , and @xmath27 is fluid velocity . through the function @xmath28 , the first and second gradient fluids depend on @xmath8 , the symmetric strain - rate tensor , @xmath29 that is , first and second gradient fluids allow energy to depend upon strain - rate @xmath8 and gradient of strain - rate @xmath30 , respectively , @xcite . the higher gradient fluids will allow energies that depend upon higher - order gradients of strain - rate . we introduced the dependence on @xmath8 , @xmath30 , etc . , instead of @xmath31 , @xmath32 , etc . , in equation ( [ lag - grade2 - 3 ] ) , so that the lagrangian @xmath19 will be invariant under rotations . consequently , the resulting euler - poincar equations will admit an angular momentum balance relation and will satisfy the requirements of material frame indifference . [ [ variational - derivatives - and - natural - boundary - conditions . ] ] variational derivatives and natural boundary conditions . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the variational derivatives of the lagrangian ( [ lag - grade2 - 3 ] ) for first and second gradient fluids are given by @xmath33 where @xmath34 , and we sum over repeated indices . the quantity @xmath35 is the stress tensor , whose definition assures that it is symmetric , @xmath36 , @xmath37 the added natural boundary condition for second gradient fluids , @xmath38 arises from an integration by parts . another application of integration by parts and use of the symmetry of @xmath39 gives @xmath40 where @xmath41 . another natural boundary condition has been introduced and applied , @xmath42 this condition may be satisfied when the fluid velocity has no normal component at the boundary , by requiring that the normal stress have no tangential component , @xmath43    [ [ the - euler - poincar - motion - equation . ] ] the euler - poincar motion equation . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the euler - poincar motion equation is @xcite @xmath44 the momentum density @xmath45 is defined as the variational derivative of the lagrangian with respect to the fluid velocity @xmath27 . for the gradient fluid lagrangian ( [ lag - grade2 - 3 ] ) , we see from that this is @xmath46 we denote @xmath47 , and @xmath48 is the material derivative along @xmath27 . the incompressibility condition @xmath49 follows from the continuity equation @xmath50 , evaluated for @xmath16 , as imposed by the pressure constraint . consequently , the euler - poincar motion equation ( [ epmotion ] ) obtained from the gradient fluid lagrangian ( [ lag - grade2 - 3 ] ) is expressed as , @xmath51 together with @xmath49 . next , we use the vector identity @xmath52 together with the coriolis relation @xmath53 , and introduce the standard dissipation law for the first and second gradient fluids . consequently , the motion equation takes the familiar form , @xmath54      [ [ kelvin - noether - circulation - theorem . ] ] kelvin - noether circulation theorem . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    in the absence of forcing and dissipation , the euler - poincar theory for lagrangians in the class ( [ lag - multipolar ] ) provides a kelvin - noether circulation theorem @xcite @xmath55 which holds for integrations around any closed curve @xmath56 moving with the fluid . for the first and second gradient fluids considered here , this becomes , @xmath57 stokes theorem then provides , for relative vorticity @xmath58 , that @xmath59 for any surface @xmath60 whose boundary @xmath61 moves with the fluid . consequently , we find the helmholtz vortex dynamics equation for the total vorticity , in the absence of forcing and dissipation , as @xmath62 thus , the kelvin - noether circulation theorem in the euler - poincar framework implies that the total vorticity @xmath63 is frozen into the flow of a non - newtonian , first or second gradient fluid . hence , its total vorticity @xmath63 satisfies the helmholtz vortex dynamics equation ( [ helm - vort - dyn ] ) . [ [ energy - conservation . ] ] energy conservation . + + + + + + + + + + + + + + + + + + + +    from the euler - poincar theory , one may compute the hamiltonian from the lagrangian @xmath22 in equation ( [ lag - grade2 - 3 ] ) for first and second gradient fluids in a rotating frame by applying the legendre transformation , @xmath64 the corresponding conserved energy is found by evaluating this expression on the constraint manifold , @xmath16 , as @xmath65 the surface integrals in the last two equations vanish , upon applying the normal - stress boundary condition ( [ normal - stress - bc ] ) , for the situation in which the velocity @xmath27 on the surface has no normal component . as a consequence , the inner product of the fluid velocity @xmath27 with the motion equation ( [ motion - eqn ] ) yields , @xmath66 the surface integral vanishes , in this energy balance relation for first and second gradient fluids in a rotating frame , provided @xmath27 vanishes on the boundary .    [ [ momentum - conservation . ] ] momentum conservation . + + + + + + + + + + + + + + + + + + + + + +    we express the euler - poincar equation ( [ epmotion ] ) in components as @xmath67 observe that for a gradient fluid lagrangian ( [ lag - multipolar ] ) given by @xmath68 we have variational derivatives @xmath69 where the @xmath70 refer to any dependence of the lagrangian density @xmath71 on higher spatial derivatives of @xmath3 and @xmath72 . therefore , upon performing the indicated differentiations by parts , one eventually finds the local conservation law for momentum , @xmath73 and momentum - stress tensor @xmath74 defined by @xmath75 here we abbreviate , by using variational - derivative notation to denote , @xmath76    the momentum conservation form ( [ mom - form ] ) is guaranteed by the euler - poincar equation for any choice of lagrangian that does not depend explicitly on the spatial coordinate . the coriolis vector potential @xmath24 introduces explicit spatial dependence into the lagrangian . consequently , although not all components of the momentum will be conserved , we may still write the motion equation ( [ motion - eqn ] ) as a * momentum balance relation * , @xmath77 where the momentum - stress tensor @xmath78 is given by ( [ stress - tens ] ) and @xmath79 is the completely antisymmetric tensor density , with @xmath80 .    equations ( [ gen - momentum ] ) for the momentum density @xmath81 and ( [ stress - tens ] ) for the momentum - stress tensor @xmath74 indicate how the derivation and analysis may proceed within the euler - poincar framework for gradient fluids of degree three , four , five , etc . these generalizations correspond to allowing the strain - rate @xmath82 in the lagrangian ( [ lag - grade2 - 3 ] ) to depend on higher gradients of the strain - rate @xmath8 . pursuing this direction further for @xmath83 degree gradient fluids is straightforward within the euler - poincar framework . however , the present paper stops at @xmath7 degree gradient fluids . [ [ choice - of - energy - density - wmathbfenablamathbfe - for - first - and - second - gradient - fluids . ] ] choice of energy density @xmath28 for first and second gradient fluids . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    in this paper , we will examine elliptic instability via exact nonlinear craik - criminale ( cc ) solutions for specific cases that apply for first and second gradient fluids . for this study , we shall choose the strain - rate dependence in the potential energy density as a norm , @xmath84 where @xmath85 and @xmath86 in tensor notation . see also bellout for a discussion of the role of this norm in proving the regularity properties of their rivlin - ericksen - green multipolar fluids . the case @xmath87 corresponds to the equations for second gradient fluids , and @xmath88 corresponds to the classic ns equations . for the choice in , we have @xmath89 upon defining @xmath90 , the motion equation ( [ motion - eqn ] ) takes the following form : @xmath91 for this choice of the energy density , the stress tensor in has the form : @xmath92 note , this stress tensor is not symmetric . the lagrangian in is also not invariant under rotations , when the coriolis vector potential @xmath24 is a fixed vector . in the absence of @xmath24 , this lagrangian regains invariance under rotations and the angular momentum in that case is conserved . however , the stress tensor in that case is still not symmetric . a solution to on an unbounded domain may be obtained , by taking velocity in the linear form , @xmath93 together with a pressure @xmath94 , which is quadratic in space . the matrix @xmath95 is a time dependent matrix such that @xmath96 and @xmath97 is the instantaneous velocity field at the origin . here , @xmath98 is a symmetric matrix defined as @xmath99 , where @xmath100 a typical approach is to choose a matrix @xmath95 for which the left hand side of is symmetric . then , the corresponding pressure @xmath101 is determined _ a posteriori _ by . we nondimensionalize the system using the variables @xmath102 , @xmath103 , @xmath104 , @xmath105 , @xmath106 , @xmath107 , where @xmath108 is a typical length scale and @xmath109 . the resulting equation with the prime notation suppressed is with @xmath110 replaced by @xmath111 . we construct a second solution to of the form @xmath112 with corresponding pressure @xmath113 . we refer to @xmath114 as the ` base ' flow and @xmath115 as the ` disturbance . ' the equations governing the disturbance are @xmath116 with @xmath117 , in which we mix tensor and vector notation , where @xmath118 and @xmath119 for @xmath120 , and @xmath121 . in the above equation , we have used the fact that @xmath114 is an exact solution to together with the vector identity @xmath122 for any two vectors @xmath123 . we choose the disturbance to be of the form @xmath124 @xmath125 , and @xmath126 and @xmath127 are scaling factors so that we can choose the initial conditions @xmath128 and @xmath129 . the unknown phase @xmath130 and the amplitudes @xmath131 , @xmath132 , and @xmath133 are to be determined . the incompressibility condition @xmath117 gives @xmath134 form this equation it follows that the nonlinear term @xmath135 in vanishes exactly . thus , in what follows , the sum @xmath112 is an exact solution to the nonlinear equations of motion in . by collecting on powers of @xmath136 and @xmath137 , the evolution equations for the amplitudes and phase are @xmath138 here , @xmath139 @xmath140 is the vorticity based ekman number , @xmath141 is the total vorticity of the system , and @xmath142 . note that the amplitude scaling @xmath126 is immaterial . without loss of generality , we set @xmath143 then taking the gradient of becomes @xmath144 we eliminate the pressure term by taking the dot product of with @xmath145 and by using @xmath146 , the first of which follows from and the second from : @xmath147}\cdot { { { \mathbf k } } }   }            + {   { { { \mathbf \pi}}\times{{\mathbf a}}}\cdot { { { \mathbf k } } }   } \big \}.\end{aligned}\ ] ] in summary , we have obtained a new exact incompressible solution to . the variables are amplitude @xmath131 and wave vector @xmath148 . once these are determined , the pressure terms follow from and . note that @xmath114 and @xmath112 are exact solutions to the nonlinear equations , but @xmath115 by itself is only a solution to linearized about @xmath114 . the exception is that in a rotating coordinate system ( @xmath149 ) , @xmath115 is an exact solution by itself since this scenario corresponds to @xmath150 in a non - rotating frame , where @xmath151 is rigid body rotation about the @xmath152axis ; cf . we emphasize that the operator @xmath153 acting on a vector represents the complete time derivative of that quantity in a lagrangian frame moving with @xmath114 . finally , the equation for @xmath131 is @xmath154 we examine the stability of a rotating column of fluid with elliptic streamlines whose foci lie on the @xmath155-axis : @xmath156 here , @xmath157 is the eccentricity of the ellipses , and the pressure is @xmath158 . with @xmath159 is analytically solvable : @xmath160^t\end{aligned}\ ] ] where @xmath161 and @xmath162 is the polar angle that @xmath145 makes with the axis of rotation . has the form @xmath163 where the elements of the matrix @xmath164 are periodic with period @xmath165 , the period of @xmath148 . therefore , the system can be analyzed numerically using floquet theory @xcite . we compute the monodromy matrix @xmath166 , that is , the fundamental solution matrix with identity initial condition evaluated at @xmath167 . will have exponentially growing solutions if @xmath168 , where @xmath169 are the eigenvalues of @xmath166 , with corresponding lyapunov - like growth rates given by @xmath170 thus , we can simulate numerically the solution to over one period and indisputably determine the exponential growth rates . we can be certain that at least one of the eigenvalues will always be unity because of the incompressibility condition and that the remaining two eigenvalues appear as complex conjugates on the unit circle or as real valued reciprocals of each other .      for flows with circular streamlines ( @xmath171 ) , the monodromy matrix can be analytically computed . it follows from that @xmath172 . then , @xmath173 is constant in time ( denoted by @xmath174 ) and has three linearly independent solutions : @xmath175 where @xmath176 , @xmath177^t$ ] and @xmath178^t$ ] are vectors orthogonal to @xmath145 , and @xmath179 is an arbitrary phase . clearly the first two solutions @xmath180 and @xmath181 satisfy . the monodromy matrix can be constructed from these three solutions : @xmath182 the three eigenvalues are @xmath183 . all of the eigenvalues lie on the unit circle , from which it follows that all solutions in the inviscid case for @xmath171 are stable . the values of @xmath184 for which @xmath185 are called ` critically stable ' and are given by @xmath186 , @xmath187 . at these parameter values an exponentially growing solution can appear ( together with an exponentially decaying one ) as @xmath188 increases from zero . bayly @xcite argues that the evenness of @xmath189 as a function of @xmath145 implies that the eigenvalues , if real and unequal , must be positive . this dismisses the odd choices of @xmath190 . furthermore , floquet theory is not applicable for the case @xmath191 . thus , the possible choices for critical stability are @xmath192 , @xmath193 . this corresponds to @xmath194 these are the critical parameter value at which @xmath131 suffers exponential growth as @xmath188 increases from zero . for the ns equations ( i.e. @xmath195 ) , only the @xmath6 choice ( called the ` principle finger ' ) is physically interesting . the other choices of @xmath1 are extremely thin fingers with growth rates ten orders of magnitude smaller than that of the principle finger @xcite . as @xmath196 and/or @xmath197 increase from zero , however , the fingers widen and the associated growth rate increases . finally , since @xmath198 , we conclude that there exists a band of stable eccentricities for @xmath199 additional understanding of this result emerges by following the analysis of waleffe @xcite and kerswell @xcite . by taking the dot product of with @xmath200 , we obtain ( for all @xmath188 and @xmath173 ) @xmath201 one can determine an exponential growth rate to leading order in @xmath188 by inserting the zeroth order solutions for @xmath145 and @xmath180 into the right hand side of this equation : @xmath202     + \frac{2\gamma(\upsilon_0 - 1)}{\upsilon_0}\sin^2\theta\sin(2t),\end{aligned}\ ] ] where @xmath203 . upon averaging over a period of @xmath180 , this quantity will vanish except when @xmath204 , corresponding to @xmath205 . compare this with . the maximum values for @xmath206 will occur at @xmath207 for @xmath204 , respectively , with growth rate @xmath208 valid for @xmath209 and @xmath210 not satisfying . thus , we see that the maximum growth rate increases as a function of @xmath196 and @xmath197 due to the @xmath211 dependence of the critical stability point up to a maximum of @xmath212 , after which a set of stable solutions emerges in a band of nonzero eccentricities . see fig . [ fig : sample ]    , @xmath213 , @xmath214 , @xmath215 computed on a @xmath216 grid for various values of @xmath210 : ( a ) -0.5 , ( b ) 0 , ( c ) 1.0 , and ( d ) 2.5 . note that the individual fingers touch the @xmath184 axis according to . [ fig : sample ] , width=480 ]    for nonzero values of @xmath188 , we must investigate the system numerically . we use the variable coefficient ordinary differential equation solver @xcite . the level surface of the growth rate for fixed @xmath217 is seen in fig . [ fig : bigfigs1 ] , and fig . [ fig : a1a2 ] shows the growth rate surface maximized over the @xmath218 plane as a function of @xmath219 .     for @xmath215 , @xmath220 , @xmath221 , @xmath222 , and various @xmath197 . figure ( a ) shows the neutral surface for @xmath223 and is an expansion of the boxed region in figure ( b ) . for @xmath224 , the critical stability point occurs at @xmath225 , which agrees with the classical ns results . the critical stability point shifts towards @xmath226 as @xmath197 increases according to @xmath227 . as @xmath228 exceeds unity , a stable band of rotating flows with nonzero eccentricities appears . the corresponding surfaces for @xmath224 and various @xmath196 are qualitatively the same , . see @xcite . [ fig : bigfigs1 ] , title="fig:",width=240 ]   for @xmath215 , @xmath220 , @xmath221 , @xmath222 , and various @xmath197 . figure ( a ) shows the neutral surface for @xmath223 and is an expansion of the boxed region in figure ( b ) . for @xmath224 , the critical stability point occurs at @xmath225 , which agrees with the classical ns results . the critical stability point shifts towards @xmath226 as @xmath197 increases according to @xmath227 . as @xmath228 exceeds unity , a stable band of rotating flows with nonzero eccentricities appears . the corresponding surfaces for @xmath224 and various @xmath196 are qualitatively the same , . see @xcite . [ fig : bigfigs1 ] , title="fig:",width=240 ]     maximized over the @xmath229 plane as a function of @xmath196 and @xmath197 for @xmath215 , @xmath230 , ( a ) @xmath222 and ( b ) @xmath231 . in figure ( a ) , the growth rate for @xmath232 ( corresponding to the classic ns case ) is @xmath233 . we see that the growth rate quickly increases to unity and then slowly decays . not shown is that @xmath234 as @xmath235 . [ fig : a1a2 ] , title="fig:",width=192 ]   maximized over the @xmath229 plane as a function of @xmath196 and @xmath197 for @xmath215 , @xmath230 , ( a ) @xmath222 and ( b ) @xmath231 . in figure ( a ) , the growth rate for @xmath232 ( corresponding to the classic ns case ) is @xmath233 . we see that the growth rate quickly increases to unity and then slowly decays . not shown is that @xmath234 as @xmath235 . [ fig : a1a2 ] , title="fig:",width=192 ]      the solutions to must be simulated numerically for @xmath236 . an interesting feature of this equation is that , unlike the ns equations , a change of variables will not remove viscosity from the problem . however , the qualitative results for ns hold true here . viscosity stabilizes the flow by lowering the maximum growth rate and introducing a stable band of eccentricities . this stabilization is slower than its ns counterpart , that is , the disspation in is of the form @xmath237 , not @xmath238 . see figs . [ fig : viscfixedbeta ] and [ fig : visfixedalpha ] .     for @xmath239 , @xmath230 , @xmath222 , as functions of @xmath196 for ( a ) @xmath240 , ( b ) @xmath241 , and ( c ) @xmath242 , . again , the corresponding surfaces for @xmath217 and various @xmath197 is qualitatively the same . [ fig : viscfixedbeta ] , title="fig:",width=192 ]   for @xmath239 , @xmath230 , @xmath222 , as functions of @xmath196 for ( a ) @xmath240 , ( b ) @xmath241 , and ( c ) @xmath242 , . again , the corresponding surfaces for @xmath217 and various @xmath197 is qualitatively the same . [ fig : viscfixedbeta ] , title="fig:",width=192 ]   for @xmath239 , @xmath230 , @xmath222 , as functions of @xmath196 for ( a ) @xmath240 , ( b ) @xmath241 , and ( c ) @xmath242 , . again , the corresponding surfaces for @xmath217 and various @xmath197 is qualitatively the same . [ fig : viscfixedbeta ] , title="fig:",width=192 ]     for @xmath222 , @xmath243 as a function of the wavenumber @xmath127 for ( a ) @xmath239 , @xmath244 and ( b ) @xmath217 , @xmath245 . for @xmath246 , the flow is inviscid . as @xmath127 increases , the leading order term in an asymptotic exapnsion of will be @xmath247 . since @xmath248 , viscous dissipation quickly takes over . [ fig : visfixedalpha ] , title="fig:",width=240 ]   for @xmath222 , @xmath243 as a function of the wavenumber @xmath127 for ( a ) @xmath239 , @xmath244 and ( b ) @xmath217 , @xmath245 . for @xmath246 , the flow is inviscid . as @xmath127 increases , the leading order term in an asymptotic exapnsion of will be @xmath247 . since @xmath248 , viscous dissipation quickly takes over . [ fig : visfixedalpha ] , title="fig:",width=240 ] future studies may investigate the roles of other aspects of nonlinear stress on elliptic instability , for example , in the rivlin - ericksen - green multipolar fluids analyzed in bellout @xcite . brf thanks the los alamos national laboratory for its hospitality , and gratefully acknowledges the financial support of the laboratory s turbulence working group and theoretical division . ddh is grateful for partial support by imperial college london and by us doe , under contract w-7405-eng-36 for los alamos national laboratory , and office of science ascar / ams / mics .
Abstract:

Answer: the energy of an @xmath0gradient fluid depends on its eulerian velocity gradients of order @xmath1 . a variational principle is introduced for the dynamics of @xmath0gradient fluids and their properties are reviewed in the context of noether s theorem . the stability properties of craik - criminale solutions for first and second gradient fluids are examined .


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: recently there has been increased interest in the possibility of using organic molecules as electronic components in nanoscale devices . one of the most important issues is the role of molecule - metal contact . the charge injection at the metal - molecule interface is usually modelled by the energy dependent probability for a substrate electron to tunnel to the molecular wire . the magnitude of the tunnelling probability depends on lineup of the molecular wire electronic states and the metal fermi level . when a closely packed array of molecular wires is formed on the metal electrode , the molecular states form bands dispersive in the lateral direction . this suggests that the tunnelling probability in addition to the energy of the injected electron depends also on the lateral component of the electron wavevector . the change of the metal work function due to a self - assembled monolayer ( sam ) may be employed in device applications to assist the charge injection into the molecular wire @xcite . the work function of a crystal surface has generally two contributions @xcite : ( 1 ) the electrostatic barrier due to the distortion in the charge distribution at the surface and ( 2 ) the many - body effects of the screened hole . it is important to know the relative contributions of the two components for the successful design of molecular electronic devices .    in this paper we present density functional theory ( dft ) calculations of the electronic structure of sams of phenylthiolate ( pt = -s - c@xmath0h@xmath1 ) on noble metal substrates cu(111 ) and au(111 ) . we find significant lateral dispersion of the molecular electronic states derived from the molecular orbitals of the sam . our calculations aim to address three essential questions : ( 1 ) what are the molecular states closest to metal fermi level and therefore the most important for transport properties ; ( 2 ) how large is the lateral dispersion of those states ; ( 3 ) what is the largest contribution to the workfunction change associated with sam formation , for which we find the electrostatic contribution from the polar molecules as determined by the dipole moment of the isolated array of molecules in addition to that of the bare metal surface . in addition we calculate two - photon photoemission ( 2ppe ) spectra using dft wavefunctions and the work function of the sam / cu(111 ) system to compare with the experiments by zhu _ et . we use the full potential linearized augmented plane wave ( fp - lapw ) method @xcite with local orbital extensions @xcite in the wien2k implementation @xcite . the gga @xcite exchange - correlation potential was used . well - converged basis sets were employed with a 4.7 ry plane wave cut off . five special k - points were used to sample the two dimensional brillouin zone . the surface plane was taken as the @xmath2 plane . we used slab calculations based on a super cell with three metal layers covered on one side by a pt monolayer for electronic structure calculations , and six metal layers with both sides covered for work function calculations . the pt monolayer is of the so called @xmath3 type @xcite , with one molecule per three substrate surface metal atoms . the former are arranged in a triangular lattice with side @xmath4 , where @xmath5 is the side of the triangular lattice formed by metal atoms @xmath6 in the ( 111 ) plane . for @xmath6=au and cu , the corresponding cross - section area for molecule is , respectively , 21.6 @xmath7 and 16.9 @xmath7 @xcite . the cross - sectional mean area of a pt molecule ( when directed normal to the substrate ) is 21.1 @xmath7 , slightly less than the above value based on the au substrate , thus indicating the possibility of a small degree of tilt @xcite , whereas no tilt is expected in the case of the cu substrate . the geometry of the molecular array has been fully relaxed on an unrelaxed substrate with metal atoms in ideal crystal positions . the monolayer on the au substrate was found tilted by 18@xmath8 , consistent with related experimental data @xcite . in the calculation , the pt molecule was constrained so as to maintain @xmath9 intramolecular symmetry @xcite . the lowest - energy binding site was the au - au bridge , analogous to the case of alkanethiols studied by selloni _ et . the s atom height is 2.1  . starting with alignment normal to the substrate , with the projection of the pt plane parallel to a nearest - neighbor au - au vector , the pt was tilted in the direction perpendicular to its plane ( tilting in the orthogonal direction caused a sharp rise in energy ) . the calculated tilt relaxation energy per molecule ( @xmath10 to @xmath11 ) was of the order of @xmath12 at room temperature , indicating the likelihood of appreciable thermal fluctuations , similar to the situation noted in @xcite .    in the case of the cu substrate we find zero tilt angle for the monolayer due to the close spacing controlled by the cu lattice constant , with the lowest energy when sulfur binds to the hollow site . the height of the s atom is 1.73  . the binding energy for the bridge site is higher by 170 mev . the frozen lattice approximation used here may cause a shift of the energy difference in binding energy of a few tens of mev for different adsorption sites , but it should not affect much the work function and the shape of the dispersion curves . the electronic structure along the high symmetry direction of the 2d brillouin zone is shown in fig . ( [ fig2]a ) for the pt sam on cu(111 ) . in order to distinguish which bands originate from the sam and which are due to the metallic substrate we calculate projected density of states of each band inside the spheres surrounding the atoms of the pt molecule @xcite . the cu surface states were identified by the similarities of the band dispersion of the covered and clean cu(111 ) three layer substrates shown on fig . ( [ fig2]b ) . at least seven molecular bands due to the sam can be identified . these molecular bands can be characterized by ( 1 ) the energy at the @xmath13-point , ( 2 ) the effective mass for the @xmath13-m direction in cases where a parabolic fit of the dispersion in the vicinity of the @xmath13 point is possible , ( 3 ) the bandwidth defined as the energy difference between the m and the @xmath13 point , and ( 4 ) the major atomic contributions . these band characteristics are summarized in table [ tab1 ] . .[tab1 ] analysis of the electronic bandstructure for the pt sam on cu(111 ) , sam on au(111 ) , and the free molecular arrays at the same geometries . the energies of the bands are given relative to the fermi level at the @xmath13-point ( @xmath14 ) . effective masses are reported for the bands where a parabolic fit ( along the @xmath13-m direction ) is possible . the bandwidth is reported for bands where the energy difference between the m and the @xmath13 point can be identified ( m denotes @xmath15 ) . the last column shows the band character for the sam / substrate systems , where c@xmath16 stands for side phenyl carbons and c@xmath17 for next to sulfur and top carbons . [ cols="^,>,>,>,>,>,>,^ " , ]      ( color online ) . band structure results ( energies relative to e@xmath18 ) . ( a ) pt sam on 3 layers of cu(111 ) . the size of the colored circles is proportional to the projected density of states of each band on a given atom of a molecule : sulfur - blue , c@xmath16 carbons - light blue , and c@xmath17 carbons - green . ( b ) 3 layer slab of cu(111 ) . from comparison of ( a ) and ( b ) four homo and three lumo bands can be identified due to the sam attachment.,width=303,height=226 ]    for the copper substrate the molecular levels closest to the fermi level are the two homo s ( bands 3 and 4 ) , which are of s ( 3p ) character . they have a very small dispersion ( of about 0.2 - 0.3 ev ) and the parabolic fits result in large effective masses ( -2.0 @xmath19 ) . the second occupied state ( ss@xmath20 , at -0.9 ev ) and the second unoccupied state ( ss@xmath21 , at 1.7 ev ) are surface states , which are at energies 0.0 ev and -1.6 ev , respectively , on the clean substrate in fig . ( [ fig2]b ) . the discrepancy with the experimental value of the -0.4 ev surface state energy at @xmath13 point is due to the finite size effect of our three layer slab geometry @xcite . the largest lateral dispersion is due to the bands 1 and 5 originating from the side phenyl carbons ( c@xmath22 for short ) . ( the minimal separation between such atoms on neighboring pt molecules is 2.61  . ) the molecular bands 2 , 6 , and 7 are mainly due to the next to sulfur and top carbon atoms ( c@xmath23 for short ) . these c@xmath23 atoms form the longest intermolecular contacts and the lateral dispersion of the c@xmath23-type bands is not significant compared to that of side carbon c@xmath22-type bands . there are certainly more electronic states in the band structure due to the sam , as is evident from the band count in fig . ( [ fig2 ] ) , but these are high energy states ( relative to fermi level ) and are not important for the conduction properties of the monolayer . these states are also well delocalized in space such that projected density of states inside the muffin tin spheres does not give a unique identification of the molecular bands . as will be clear in section [ sec4 ] , the band 6 at 2.8 ev ( the @xmath13 point ) and one at 6 ev have the largest contribution to the 2ppe intensity and thus are candidates for the peaks observed by zhu et . al . @xcite . ( color online ) . band structure results ( energies relative to e@xmath18 ) ( a ) pt sam on 3 layers of au(111 ) . the size of the colored circles is proportional to the projected density of states of each band on a given atom of a molecule : sulfur - blue , c@xmath16 carbons - light blue , and c@xmath17 carbons - green . ( b ) 3 layer slab of au(111).,width=302,height=217 ]    when the pt is attached to the gold substrate , the bandwidths of all molecular levels are narrower than for the cu case , because the nearest neighbor distance between the molecules is larger due to the larger lattice constant of the underlying substrate . the c@xmath22 atoms on neighboring molecules of the sam ( 18% tilted on au ) are separated by at least 3.13   from each other , whereas c@xmath23 carbons are 5.0   apart ( i.e. , @xmath24 ) . the electronic band structures of the sam / au(111 ) and the clean au(111 ) surfaces are shown in fig . ( [ fig3 ] ) . the band assignments are summarized in table ( [ tab1 ] ) . as in the case of the sam on cu , the two nearly degenerate sulfur - based homo s ( bands 3 and 4 ) are the molecular levels closest to the fermi level . as the lattice constant increases ( in proceeding from cu to au ) , the direct exchange mechanism responsible for the negative effective mass decays faster than the substrate mediated mechanism ( see below ) . this causes the sign of the effective mass of one of the sulfur bands to become positive on the au substrate . the most dispersive c@xmath22 bands 1 and 5 have bandwidth smaller by a factor of 2.5 than for the cu substrate . the occupied surface state ss@xmath20 is strongly hybridized with the sulfur 3p orbitals . ( color online ) . ( a ) band structure of the free standing pt array , based on the cu(111 ) lattice constant and ( b ) a pt array tilted by 18@xmath8 , based on the au(111 ) lattice constant ( energies relative to e@xmath18 ) . the size of the colored circles is proportional to the projected density of states of each band on a given atom of a molecule : sulfur - blue , c@xmath25 carbons - light blue , and c@xmath17 carbons - green.,width=313,height=189 ]    the lateral dispersion has two contributions the direct overlap between the molecules and the substrate mediated interaction . to estimate the relative importance of the two mechanisms for the band dispersion we compare the lateral dispersions obtained for the pt arrays in the absence of the substrates ( infinitely separated from the substrate ) , but maintaining the same geometry as determined with the substrate present . in this case the band structure shown in fig . [ fig4 ] has only a direct contribution . the signs of the dispersions for bands 3 and 4 are negative and the bandwidth is smaller for the au(111 ) case ( fig . ( [ fig4]b ) ) . in both cases the dispersion of the sulfur bands is larger in the absence of the substrate . this suggests that substrate - mediated interactions have an effect on the bandwidth opposite to that from the direct exchange . when the lattice constant is relatively large ( i.e. the au case ) , the substrate - mediated interaction even changes the sign of the band 4 dispersion . the coupling of the c@xmath16 atoms to the substrate is much smaller , so that substrate - mediated interaction in bands 1 and 5 is responsible for only 10% of the total bandwidth . the electronic structure of the organic sam can be directly probed in the 2ppe experiment @xcite . there are three electronic states involved in the 2ppe spectra : initial ( _ i _ ) , intermediate ( _ k _ ) , and the final ( _ f _ ) . if one of those states is localized and has the largest dipole matrix element in comparison with the rest of the states , then three scenarios are possible for the variation of kinetic energy ( @xmath26 ) with respect to the excitation photon energy @xmath27 , which is chosen to be less than the workfunction of the substrate : ( 1 ) @xmath28 ; ( 2 ) @xmath29 , and ( 3 ) invariance with respect to @xmath27 . for pt on cu(111 ) , two cases are reported @xcite : case(2 ) in which the intermediate state ( i ) is assigned as a lumo lying below the vacuum at 3.3 ev and case(3 ) , in which the final state ( f ) is assigned as a lumo lying above the vacuum at 6.4 ev ( relative to e@xmath18 ) . case ( 1 ) , corresponding to an initial homo state and was not observed in @xcite .    in the present work we use dft wavefunctions to calculate the 2ppe spectra at the @xmath13-point . we find two states to contribute most to the 2ppe intensity : a lumo ( band 6 ) at 2.8 ev ( of c@xmath17 carbon character ) and high energy bands at around 6.0 ev , whose character is complex due to hybridization of molecular and cu states . this assignment is consistent with the conclusion drawn in @xcite based on similarities of spectra taken for sams based on pt and alkanethiolates of different lengths . there are two possible 2ppe excitation mechanisms : direct and indirect . in the direct process with one - color pump and probe laser frequencies , the 2ppe intensity for @xmath30-light polarization is @xcite : @xmath31 where @xmath32 is the lifetime broadening of intermediate level @xmath33 and a denotes the cartesian component(x , y , or z ) . the sum over intermediate states in eq . ( [ ppd ] ) includes interference effects . the delta function may be replaced by a lorentzian to mimic the lifetime of the final state and the width of the laser pulse . the result for the @xmath13-point is shown in fig . [ fig6 ] where intensity contributions eq . ( [ ppd ] ) are summed over all initial states _ i_. the calculated curves were convoluted with a lorentzian function ( @xmath34 ev ) to account for the spectrometer resolution . we use intermediate and final level lifetimes @xmath35 ev . there is one pronounced peak at about 6.0 ev , whose intensity is resonant with the laser frequency at about 3.6 ev . as the excitation frequency increases the intensity transfers to higher energy final states ( @xmath36 ev ) . the resonance enhancement factor is sensitive to the choice of the damping parameters , and even for the large value ( 0.1 ev ) used in fig . ( [ fig6 ] ) , the intensity of the out - of - plane polarization varies by an order of magnitude . the resonance effect is related to the large dipole matrix element between the occupied surface state ss@xmath20 and the second lumo ( band 6 ) separated by 3.7 ev . for the in - plane polarization there is no strong laser frequency dependence of the 2ppe intensity , because there is no single pair of states which dominates the sum in eq . ( [ ppd ] ) . therefore the ratio of the out - of - plane and the in - plane polarizations in the direct mechanism model should be resonant with the laser frequency . prediction of the 2ppe spectra for the direct mechanism according to eq . ( ) for ( a ) out of plane and ( b ) in plane light polarizations.,width=285,height=209 ]    in the indirect mechanism the intermediate level can be populated by some incoherent process , for example , by the tunnelling of the photoexcited electron from the substrate to the molecule . in the second step , a photon promotes the photoexcited electron to the vacuum with probability given by the one photon absorption crossection : @xmath37 where fermi - dirac factors @xmath38 require that the intermediate and final states are empty and the intermediate state is less than one quantum of photon energy above the fermi level @xmath39 . it is assumed that all the intermediate states _ k _ have the equal lifetimes and probabilities to be excited in the first step . the result for the two polarizations and the same damping parameters as in direct process is shown in fig . there are no pronounced resonances and the ratio of intensities for the out - of - plane and the in - plane polarizations is relatively independent of the laser frequency . namely , the ratios of the integrated peak intensities at 6 ev for @xmath403.4 , 3.5 , 3.7 , and 3.9 are 1.8 , 4.6 , 3.8 , and 4.0 respectively . this is consistent with the experimentally reported value of @xmath41 at @xmath42 ev @xcite . comparison of our analysis with the available experimental data @xcite leads us to suggest that the most likely mechanism for the two photon photoemission process is indirect . an enhancement of 2ppe intensity with the laser frequency was observed for sch@xmath43 on ag(111 ) by the harris group @xcite . this is may be the same effect we find here for the direct mechanism . substrate assisted 2ppe spectra ( indirect mechanism ) calculated with eq . ( ) for ( a ) out of plane and ( b ) in plane light polarizations.,width=303,height=225 ] for the work function calculations we use a six cu layer substrate covered with a sam on both sides . the inversion symmetry of the supercell with two equivalent surfaces forces the net dipole and the electrostatic potential drop in the vacuum region to be zero . to separate the electrostatic and exchange correlation contributions to the work function we calculate the averaged surface dipole moment by integrating over the net charge density of the half unit cell : @xmath44dxdydz , \label{eq20}\end{aligned}\ ] ] where @xmath45 and @xmath46 are valence electron and ion charge densities respectively . the planes @xmath47 and @xmath48 contain the inversion centers of the supecell . the coulomb potential averaged in the lateral direction is shown on fig . [ fig8 ] along with that for the clean substrate and the two sam arrays in the absence of the substrate for the same geometry . the work function of the clean cu(111 ) is found to be 5.0 ev , in good agreement with the experimental value of 4.9 ev . the work function of the covered substrate is reduced to 2.7 ev compared with the experimental value of 3.7 ev @xcite . the isolated sam layer has surface dipole @xmath49 ( debye / molecule ) , and the expected coulomb potential drop of @xmath50 ev ( where @xmath51@xmath7 is the lateral area per molecule ) agrees well with the actual result @xmath52 ev ( see fig . ( [ fig8]b ) ) . analogously we define the electrostatic contribution to the clean cu surface workfunction to be @xmath53 ev ( see fig . ( [ fig8]a ) ) . this leaves a 0.4 ev contribution due to exchange - correlation effects . in the case of the sam on cu , the net surface dipole is approximately equal to the sum of dipoles of the clean cu surface and the isolated sam array @xmath54 debye . the difference between the simple estimate and the actual result ( @xmath55 debye ) is due to the charge transfer at the interface . the laterally averaged charge densities are shown in fig . ( [ fig9]a ) . the charge transfer electron density shown in fig . ( [ fig9]b ) is very small ( less then 0.1 @xmath56 per molecule ) . we calculate the workfunction change of the covered substrate to be -2.3 ev , where -2.8 ev is the electrostatic and 0.5 ev is the exchange - correlation contribution , respectively .    as we have shown the electrostatic contribution to the workfunction change can be immediately estimated from the dipole moment of the isolated sam array . however , the isolated sam dipole moment itself is not equal to that of isolated molecule . we have used the molecular nrlmol code @xcite and find the dipole moment of the s - c@xmath0h@xmath0 molecule to be 3.33 debye . this value is threefold larger than the dipole per molecule in the array , 1.19 debye . we have also performed constrained calculations by forcing occupations of two nearly degenerate homo orbitals : @xmath57 ( in the molecular plane ( xz ) ) and @xmath58 ( perpendicular to the molecular plane ) . the unconstrained spin restricted ( ie , pure doublet ) calculation results in 1.88 electrons in @xmath57 and 1.12 electrons in @xmath58 orbitals , with a net dipole moment of 3.33 debye . by forcing the @xmath58 orbital to be fully occupied the net dipole moment drops to 1.24 debye , very close to the sam value . this suggests that the charge redistribution between the two homo orbitals in the sam is the primary effect on the net dipole moment . ( color online ) . coulomb potential averaged in the lateral direction versus @xmath59-direction for ( a ) the clean cu(111 ) , ( b ) the single monolayer , and ( c ) the sam on cu(111 ) . the dipole moments according to eq . ( ) are shown by the arrows.,width=306,height=216 ]    ( color online ) . ( a ) the laterally average charges densities versus @xmath59-direction for the clean cu(111 ) ( blue ) - @xmath60 , isolated sam array ( red ) - @xmath61 and two - sided sam / cu(111 ) ( black ) - @xmath62 . ( b ) the laterally averaged charge density difference @xmath63 versus @xmath59-direction which integrates to the induced dipole p=0.08 ( d ) due to the charge transfer . , width=295,height=216 ] we have calculated the electronic structure of the phenylthiolate ( pt = s - c@xmath0h@xmath1 ) sam on cu(111 ) and au(111 ) substrates . we identified four homo s and three lumo s for each system . a sulfur - type homo ( band 4 ) is the closest to the fermi level for both substrates . it has two competing contributions to the lateral dispersion : ( 1 ) the direct exchange between the neighboring orbitals leads to a negative effective mass , whereas ( 2 ) the substrate - mediated interaction gives a positive effective mass contribution . the overall sign of the effective mass is changed from a negative to a positive value in going from the cu to the au substrate . the most dispersive molecular bands are of side carbon character ( c@xmath16 ) . the significant lateral dispersion of the electronic states derived from molecular orbitals of the sam suggests that the charge injection probability depends on the lateral component of the wavevector as well as the energy of the injected electron . this implies that the conductance of n molecules in an array is not simply n times the conductance of one , but further theoretical and experimental work is needed to quantify the lateral dispersion effect on transport . we have calculated the 2ppe spectra of a pt sam on cu(111 ) using the dft wavefunctions . in the direct mechanism we predict a very strong resonance of the 2ppe intensity for out - of - plane polarization . for the in - plane polarization only small intensity variation with laser frequency for the indirect mechanism we find a moderate intensity dependence on the excitation photon energy for both polarizations with a nearly constant ratio of the out - of - plane to the in - plane signals of about 3 - 4 , consistent with related experimental rsults @xcite . this fact suggests that the indirect mechanism contributes the most intensity to the 2ppe resonance observed in @xcite ( designated as resonance a ) . the energy positions of the observed molecular states at 6.4 ev and 3.4 ev are in reasonable correspondence with our calculated values of 6.0 ev and 2.8 ev , respectively . these molecular states are mainly due to lumos dominated by the pt carbon atoms adjacent to the s atoms .    calculation of the workfunction for the pt - sam on cu(111 ) yields a value of 2.7 ev , which may be compared with the experimental value of 3.7 ev . to reconcile theory and experiment one may assume a lower coverage of the sam in the experimental system . the larger effective area per molecule reduces the surface dipole moment and hence the coulomb contribution to the workfunction change . according to our calculations one has to assume @xmath64 @xmath7 to agree with experiment . this estimate is much larger than that for the @xmath3 structure , which implies @xmath51 @xmath7 . alternatively , the interpretation of the workfunction measurements in one photon photoemission spectra can be complicated if there is an exit barrier for substrate electrons to tunnel through the interface . then the onset of the one - photon photoemission signal may be due to the finite energy electrons with respect to the vacuum level . this in turn would imply that the agreement between the theory and the experiment for the energies of the lumo orbitals measured in 2ppe is accidental . indeed dft is known to underestimate the bandgaps in semiconductors . on the other hand the energy differences between the empty states are more reliable in dft and in the experiment the energy difference between the two lumo s is also independent of the workfunction value . we find an electrostatic contribution to the work function change of 2.8 ev ( associated with sam formation ) to be the dominant effect . therefore the dipole moment of the molecular sam in the absence of the metal substrate can be used to predict the the workfunction change of the covered substrate provided the charge transfer and correlation effects do not contribute much to the workfunction change as in the case studied here . the computations were performed on the bnl galaxy cluster . we are grateful to mike weinert , ben ocko , and x .- y . zhu for helpful discussions . this work was supported in part by doe grant no.de-ac-02-98ch10886 . t. vondrak , h. wang , p. winget , c. j. cramer , and x .- y . zhu , j. am . * 122 * , 4700 ( 2000 ) . d. j. singh , _ planewaves , pseudopotentials and the lapw method _ ( kluwer academic , boston , 1994 ) . d. j. singh , phys . b * 43 * , 6388 ( 1991 ) .          t. y. b. leung , p. v. schwartz , g. scoles , f. schreiber , a. ulman , surf . sci . * 458 * , 34 ( 2000 ) . the @xmath9 symmetry constrain leaves the following internal degrees of freedom for a ) the bondlengths : s - c@xmath65 = 1.792  , c@xmath65-c@xmath66 = 1.402  , c@xmath66-c@xmath43 = 1.396  , c@xmath43-c@xmath67 = 1.398  , c@xmath66-h = 1.090  , c@xmath43-h = 1.091  , c@xmath67-h = 1.090   and for b ) the bond angles : @xmath68 sc@xmath65c@xmath66=119.7@xmath8 , @xmath68 c@xmath65c@xmath66c@xmath43=119.4@xmath8 , @xmath68 c@xmath65c@xmath66h=119.9@xmath8 , @xmath68 c@xmath66c@xmath43c@xmath67=120.5@xmath8 , @xmath68 c@xmath66c@xmath43h=119.3@xmath8 , @xmath68 c@xmath43c@xmath67h=120.1@xmath8 . m. c. vargas , p. giannozzi , a. selloni , and g. scoles , j. phys . b * 105 * , 9509 ( 2001 ) . a. euceda , d. m. bylander , and l. kleinman , phys . rev . b * 28 * , 528 ( 1983 ) x.y . zhu , annu . phys . chem . * 53 * , 221 ( 2002 ) . m. r. pederson and k. a. jackson , phys . b. * 41 * , 7453 ( 1990 ) ; k. a. jackson and m. r. pederson , phys . b. * 42 * , 3276 ( 1990 ) ; m. r. pederson and k. a. jackson , phys . b. * 43 * , 7312 ( 1991 ) ; d. v. porezag and m. r. pederson , phys . b. * 54 * , 7830 ( 1996 ) ; a. briley , m. r. pederson , k. a. jackson , d. c. patton , and d. v. porezag , phys . b. * 58 * , 1786 ( 1998 ) .
Abstract:

Answer: we use first principles density functional theory to calculate the electronic structure of the phenylthiolate ( s - c@xmath0h@xmath1 ) self - assembled monolayer ( sam ) on cu(111 ) and au(111 ) substrates . we find significant lateral dispersion of the sam molecular states and discuss its implications for transport properties of the molecular wire array . we calculate the two photon photoemission spectra and the work function of the sam on cu(111 ) and compare them with the available experimental data . our results are used to discuss assignments of the observed spectral data and yield predictions for new electronic states due to the monolayer not yet accessed experimentally .


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: many adsorption systems exhibit either complete wetting below the triple point or triple point wetting transitions . these behaviors are recognized to arise from strongly attracting substrate potentials and/or small mismatch of substrate / adsorbate lattice constants @xcite . less attractive potentials exhibit incomplete or nonwetting behavior at the triple point . there is considerable interest in the nature of such weakly attractive situations near and above this temperature @xcite . general arguments and model calculations of cahn @xcite and of ebner and saam @xcite in 1977 implied that systems which are nonwetting at the triple point ought to exhibit transitions to wetting at some temperature below the critical point . experimental confirmation of this proposal was not found until appropriate systems were identified . in 1991 , it was argued that alkali metals provide the weakest adsorption potentials of any surfaces for he atoms and therefore wetting transitions were expected for such substrates at tow temperature @xcite . this hypothesis was subsequently confirmed and the corresponding transition has since been studied in many laboratories @xcite . more recently , similar phenomena have been predicted and/or seen for h@xmath0 , ne , and hg on various surfaces @xcite . key theoretical questions arising from these studies involve the reliability of both the assumed adsorption potentials and the statistical mechanical description of the transition , as well as the effects of surface heterogeneity @xcite present in any real experiment . in the course of these investigations , cheng _ @xcite ( ccst henceforth ) posited a simple model which interprets the transitions in terms of a balance between the surface tension cost of producing a thicker film and the energy gain associated with the film s interaction with the surface , v(z ) @xcite . there results an implicit equation for the wetting temperature ( t@xmath3 ) : @xmath4 @xmath5 where @xmath6 and @xmath7 are the number densities of the adsorbate liquid and vapor at coexistence , @xmath8 is the surface tension of the liquid , and @xmath9 is the equilibrium distance of the potential . if one assumes a simple 3 - 9 form of the adsorption potential @xmath10 then the transition condition becomes @xmath11 ccst used this model to predict wetting transitions for a large number of specific systems . more recently , the same general approach was used by chizmeshya _ et al . _ ( ccz ) @xcite , who based their predictions on new _ ab initio _ adsorption potentials . the latter potentials differ from earlier potentials in several respects ; especially important are differences in values of the coefficient c of the van der waals asymptotic potential , resulting in significantly larger well depths d. there is indication that these revised potentials are more consistent with experimental data than are their predecessors . particularly striking is the improved agreement for @xmath12he contact angle @xcite and wetting temperature on cs , the most studied wetting transition to date @xcite . these predictions have also been tested by recent experiments of hess _ _ @xcite and simulations of ne on alkali metals and mg @xcite . in this paper , we extend the simulation studies to a larger family of systems for which the hypothesis is relevant . we focus here on the following question : assuming that v(z ) is known , does the ccst model correctly predict whether a specific adsorption system exhibits wetting at and above the triple temperature ? the outline of this paper is the following . section [ label.3.sec:method ] summarizes our simulation methodology and presents detailed results for several systems which yield fairly distinct adsorption behaviors . these include complete wetting for all t ( kr / mg and xe / mg ) , an anomalous wetting transition near the triple point ( kr / pmg ) , a prewetting transition line ( ar / mg ) somewhat above the triple temperature , and prewetting transitions characterized by the formation of very thick films ( for ar , kr and xe on li ) .    here pmg is our own nomenclature for `` pseudo - magnesium '' , an ersatz surface which attracts kr atoms about 7% more strongly than mg itself @xcite . we also evaluate other hypothetical substrates in order to ascertain the triple point wetting threshold . in section [ label.3.sec:interpretation ] , we compare these findings with those obtained previously in simulations and/or experiments for inert gases on various surfaces @xcite . the simulation technique has been described in detail previously @xcite , so only a brief description is given here . for this study , isotherms for various model systems were calculated using the grand canonical monte carlo ( gcmc ) technique @xcite . for each isotherm , 6.8 million steps ( each step being an attempted creation , deletion , or displacement of a molecule ) were performed to reach equilibration and 4.5 million steps were performed in the data gathering phase . for some points near the transition regions , where the fluctuations are larger , the number of equilibration steps was increased to 20 million and 12 million steps were performed for data gathering .    as in our previous work , the model for the atom - atom interaction is a lennard - jones 12 - 6 potential and the atom surface interaction is modeled using the potential described by chizmeshya _ these potentials can be characterized by four parameters for each system studied . a well - depth @xmath13 and distance @xmath14 are needed to specify the adsorbate lj potential and a well - depth @xmath15 and the van der waals constant @xmath16 specify the atom - surface potential . the values of these parameters used in this study are listed in table [ label.3.table:param ] . the parameters correspond to ne , kr , ar , and xe adsorbed on alkali and alkali earth metal surfaces . the surface labeled pmg ( pseudo magnesium ) was obtained when we attempted to simulate xe on li , but used the @xmath13 and @xmath14 values for kr instead of xe . this potential is similar to that of kr on a mg surface ; however , the kr / pmg distance of closest approach is 20% smaller and the well depth is about 7% weaker than estimated for kr on mg . while this potential does not correspond to a `` real '' system , we shall see that this small decrease in interaction strength significantly affects the wetting properties of the substrate / adsorbate system .    in most cases the simulation cell had a height of h = 70  , with the metal surface forming the boundary at one end and a hard wall surface at the other . in xe / li and ar / li simulations we increased the height to 140 and 210  . the surface unit cell was taken to be @xmath17 giving a nominal monolayer coverage of roughly 100 atoms . periodic boundary conditions in these dimensions render the surface effectively infinite ; we used a fairly large cutoff ( 5@xmath14 ) to minimize long range corrections .    for each system , we simulated adsorption isotherms over relevant temperature ranges from the triple point to the critical temperature . the results are shown in figures [ label.3.fig:xeli1]-[label.3.fig:krmg ] . of the systems studied the least attractive surface of the surface interactions are for the ar / li , kr / li and xe / li cases . in all of these , we see evidence of prewetting transitions and , based on our results , we report wetting temperature . for ar / li we estimate t@xmath18 k , for kr / li , t@xmath19 k and t@xmath20 k for xe / li . all @xmath21 are reported in table [ label.3.table:param ] .    despite the fact that these systems are examples of prewetting transitions , there are some qualitative differences from the case ( say ne / mg ) where the wetting transition appears close to the triple temperature . figure [ label.3.fig:xeli1 ] shows a series of isotherms for xe / li from @xmath22 1.0 to 1.25 . the isotherm at 221 k ( the lowest @xmath23 ) shows nonwetting behavior , manifested as slight adsorption below saturated vapor pressure . the higher temperature isotherms however , show prewetting transitions below saturated vapor pressure , from which we estimate @xmath21 . the nature of the prewetting jump can be further characterized by the results in figures [ label.3.fig:xeli2 ] , [ label.3.fig:xeli3 ] , and [ label.3.fig:xeli4 ] , all at t=254 k. figure [ label.3.fig:xeli2 ] shows the dependence of the isotherm on the height @xmath24 [ @xmath25  ] of the simulation box . the rapid rise occurs at the same pressure , independent of box height . this is a clear indication that the transition is prewetting rather than capillary condensation . the detailed behavior is most easily understood by looking at the density profiles dependences on @xmath24 . figure [ label.3.fig:xeli3 ] shows these profiles at @xmath26 atm . note that the data for two largest boxes coincide but differ from those of the smallest box . evidently , this regime requires quite large values of @xmath24 in order to obtain reliable data ; this is not surprising because the compressibility of the system is particularly high at this transition as the slope of the isotherms suggest . the same phenomenon is even more dramatic for the case of a 1% higher pressure . as seen in figure [ label.3.fig:xeli4 ] , one must go to even larger @xmath24 ( @xmath27  ) in order to derive accurate adsorption data even though the liquid regime occupies only the region @xmath28  . while this sensitivity to @xmath24 and the details of the transition are quite interesting they are not the focus of this paper . suffice it to say that the density jumps discontinuously near @xmath29 atm , as seen in the different densities in figures [ label.3.fig:xeli3 ] and [ label.3.fig:xeli4 ] . we note one further point about these data which is that the film s density close to the surface is qualitatively similar to that on a relatively more attractive surface . there appears a quite large accumulation of adsorbate atoms in the region associated with the first layer . there is even a significant density maximum occurring where a dense second layer might occur . the latter is , however , preempted by the prewetting jump in the regime of temperature seen here . above the prewetting critical point , in contrast , one does expect such smooth film growth .    as the relative strength ( @xmath30 ) of the surface interaction increases , ( ar on mg , kr on pmg ) the isotherms give evidence of wetting transitions . the results for ar on mg , ( shown in figure [ label.3.fig:armg ] ) indicate prewetting transitions at about 10% above the triple point ( t@xmath31 = 81 k ) . from these data , we estimate the wetting temperature to be 90 k for ar / mg with a prewetting critical temperature of 95 k.    in figure [ label.3.fig:krpmg ] , the results for kr on pmg are shown . this system shows clear prewetting transitions below 128 k for all temperatures above the triple temperature ( estimated to be 116 k @xcite ) . although this is a fictitious system , it represents one of the first simulations of triple point wetting . the only other simulation of wetting in this region is the case of ne / mg which we reported in @xcite . finally , as the surface interaction strength further increases , one begins to see continuous growth for all t between the triple temperature and the critical point . this is evident in figure [ label.3.fig:krmg ] which shows isotherms for kr on mg . the isotherms were simulated at temperatures near the estimated triple temperature and demonstrate complete wetting on the mg surface by the adsorbate even at the lowest temperatures . the results for xe on mg ( not shown ) are qualitatively the same . to interpret our results we follow ccst , who evaluated their simple model in terms of the law of corresponding states . then we note that equation [ label.3.eq:criterion ] may be written in a universal form in terms of the reduced well depth @xmath32 and interaction range @xmath33 , defined as follows : @xmath34 the values of these parameters are reported in table [ label.3.table:param ] . the ccst model may then be written in terms of a dimensionless wetting parameter @xmath35\ ] ] the triple point experimental data summarized in ccst yield a value of 3.6 for the right side of this equation . hence the criterion is @xmath1 = 3.6 for a wetting transition precisely at the triple point . the latter value , however , presupposes a 3 - 9 adsorption potential , which we have not used in the present simulations . if , instead , we compute the integral of the actual potential used in the calculations , we find a value which is 8 to 10% higher than the result obtained with the 3 - 9 potential . taking this correction into account , the threshold for wetting at the triple point becomes somewhat lower for the realistic potential : @xmath36 this is shown in figure [ label.3.fig:threshold ] as a dashed curve . turning to the actual simulation results , we find that this prediction characterizes our data relatively well , insofar as the vicinity of the value 3.3 is a crossover regime of wetting behavior . for all gases exposed to the alkali metal surfaces the quantity @xmath1 is less than 2.3 , and the systems are nonwetting at the triple point . the case of ar / mg ( @xmath1 = 3.2 ) is one for which a wetting transition occurs about 10% above the triple temperature ; similar behavior was found in our previous study of ne / mg ( @xmath37= 3.4 ) . the case of kr / pmg involves a similar value ( @xmath1 = 3.2 ) , but we find anomalous triple point wetting transition behavior for that system . finally , the cases of kr / mg ( @xmath1 = 3.4 ) and xe / mg ( @xmath1 = 3.8 ) both correspond to complete wetting . for comparison with a more typical adsorption system , the case of ar / au ( @xmath1 = 6.7 ) is one which exhibits wetting at and above the triple temperature @xcite . we may add to this list of relevant systems the recently studied case of ar / co@xmath0 . _ @xcite found that it exhibits triple point wetting ; this finding agreed with their density functional calculation using their theoretical potential . this potential would be roughly consistent with the ccst model in that the latter would predict a wetting transition at a temperature only 3% above the triple point temperature , 83.8 k @xcite .    while the criterion @xmath1 = 3.3 successfully discriminates between the nonwetting and wetting data sets mentioned above , it is worthwhile to ascertain the validity of this criterion more accurately . hence we have performed extensive simulations at the triple temperature , varying the well depth near the values @xmath33 = 0.7 , 0.9 , and 1.4 in order to determine the wetting transition line . the results are shown in figure [ label.3.fig:threshold ] . one identifies this line as the point where the letter c is just above the letter w. note that these transition values ( dash - dot curve ) lie somewhat above the @xmath1 = 3.3 dashed line . this comparison indicates that the ccst model underestimates the well depth criterion for wetting . this finding is qualitatively consistent with results reported recently by ancilotto and toigo @xcite . using a density functional method , these workers found wetting transitions very close to the coexistence curve for several systems ( e.g. ne / li and na ) for which our simulations find nonwetting behavior . the different findings of the simulations may be due to the uncertainty in identifying wetting behavior close to saturation , especially at t approaching t@xmath38 . in these and other cases , the transitions occurred at ( typically @xmath39 10% ) higher temperature than is predicted by the ccst model . these findings are also consistent with one s realization that the ccst model predicts wetting transitions for all adsorption circumstances , however weakly attractive , while a more reliable theory yields drying transitions for _ ultraweakly _ adsorbing potentials @xcite .    to summarize , we have found the ccst equation [ label.3.eq:integral ] to be a semiquantitative guide to the occurrence of wetting transitions . it would be valuable to have additional , genuine experimental data to test this conclusion . we wish especially to encourage experimentalists to explore totally new systems , especially ar , kr , and xe adsorption on alkali metals . there are several reasons for doing so . one is that such experiments will provide tests of the weak adsorption potentials for these systems , for which no experimental data exist . a second is that triple point wetting can involve interesting behavior at and just below the triple temperature , as has been found in many experiments @xcite . we have not explored either these systems , or analogous ones , below the triple temperature . in order to reliably describe the crystalline phase , as is necessary there , such simulations would require a more careful handling of the unit cell dimensions and shape than has been employed here . in view of the many relevant systems , such a study seems well justified for the future . this research was supported by the national science foundation . it has benefited from discussion with f. ancilotto , v. bakaev , m. h. w. chan , j. z. larese , g. mistura , and f. toigo . we are especially grateful to _ fondazione aldo ing . gini _ for the generous support of ing . s. curtarolo .      c|ccccccc gas / surface & & rb & na & li & pmg & mg & co@xmath0 + ne & @xmath32 & 0.7 & 1.1 & 1.5 & & 2.8 & + @xmath40 k & @xmath33 & 1.6 & 1.4 & 1.4 & & 1.2 & + @xmath41 &behavior&d & n & n & & tp(22k ) & + & @xmath42 & 1.1 & 1.5 & 2.1 & & 3.4 & + ar & @xmath32 & 1.1 & 1.6 & 2.1 & & 3.5 & 3.8 + @xmath43 k & @xmath33 & 1.1 & 1.1 & 1.0 & & 0.9 & 1.0 + @xmath44 &behavior&n & n & w(130k ) & & w(90k ) & c + & @xmath42 & 1.2 & 1.7 & 2.1 & & 3.2 & 3.8 + kr & @xmath32 & & & 2.3 & 3.6 & 3.8 & + @xmath45 k & @xmath33 & & & 0.9 & 0.9 & 0.9 & + @xmath46 &behavior & & & w(175k)&tp(116k)&c & + & @xmath42 & & & 2.1 & 3.2 & 3.4 & + xe & @xmath32 & & & 2.8 & & 5.45 & + @xmath47 k & @xmath33 & & & 0.8 & & 0.7 & + @xmath48 &behavior & & & w(225k ) & & c & + & @xmath42 & & & 2.2 & & 3.8 & +                    j. rutledge and p. taborek , phys . . lett . * 69 * , 937 ( 1992 ) ; j. e. rutledge and p. taborek , j. low temp . phys . * 95 * , 405 ( 1994 ) ; d. ross , j. a. phillips , and p. taborek , j. low temp . phys . * 106 * , 81 ( 1997 ) ; d. reinelt , h. gau , s. herminghaus and p. leiderer , czech . j. phys . * 46 * , suppl . s1 , 431 ( 1996 ) . k. s. ketola , s. wang , and r. b. hallock , phys . lett . * 68 * , 201 ( 1992 ) ; r. b. hallock , j. low temp . phys . * 101 * , 31 ( 1995 ) ; a. f. g. wyatt , j. klier , and p. stefanyi , phys . 74 * , 1151 ( 1995 ) . d. ross , j. e. rutledge , and p. taborek , fluid phase equilibria * 150 - 151 * , 599 ( 1998 ) ; d. ross , j. e. rutledge , and p. taborek , science * 278 * , 664 ( 1997 ) ; j. a. phillips , p. taborek , and j. e. rutledge , j. low temp . phys . * 113 * , 829 ( 1998 ) ; j. a. phillips , d. ross , p. taborek , and j. e. rutledge , phys . b * 58 * , 3361 ( 1998 ) ; d. ross , p. taborek , and j. e. rutledge , phys . rev . b * 58 * , r4274 ( 1998 ) . the alert reader might notice that the adsorption potential parameters for kr / pmg coincide with those for xe / li ; this reveals the accidental origin of this serendipitous choice of fictitious substrate . prior to the recent work of mistura _ et al . _ , much theoretical attention had been devoted to the ar / co@xmath0 problem , but the assumed potential was less attractive than that implied by the recent study . @xcite . a. d. migone , j. g. dash , m. schick , and o. e. vilches , phys . b * 34 * , 6322 ( 1986 ) ; a. d. migone , a. hoffmann , j. g. dash , and o. e. vilches , phys . rev . b * 37 * , 5440 ( 1988 ) ; m. s. pettersen , m. j. lysek and d. l. goodstein , surf . sci . * 175 * , 141 ( 1986 ) ; l. bruschi and g. mistura , phys . b * 58 * , 1181 ( 1998 ) ; j. z. larese and m. y. m lee , phys . rev . * 79 * , 689 ( 1997 ) .
Abstract:

Answer: grand canonical simulations are used to calculate adsorption isotherms of various classical gases on alkali metal and mg surfaces . _ ab initio _ adsorption potentials and lennard - jones gas - gas interactions are used . depending on the system , the resulting behavior can be nonwetting for all temperatures studied , complete wetting , or ( in the intermediate case ) exhibit a wetting transition . an unusual variety of wetting transitions at the triple point is found in the case of a specific adsorption potential of intermediate strength . the general threshold for wetting near the triple point is found to be close to that predicted with a heuristic model of cheng _ et al . _ this same conclusion was drawn in a recent experimental and simulation study of ar on co@xmath0 by mistura _ et al . _ these results imply that a dimensionless wetting parameter @xmath1 is useful for predicting whether wetting behavior is present at and above the triple temperature . the nonwetting / wetting crossover value found here is @xmath2 3.3 .    _ pacs numbers : 64.70.fx , 68.35.rh , 68.45.gd , 82.20.wt_