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 granular character of many processes in nature has motivated the study of granular materials and in particular granular packings since kepler s time@xcite . one of such processes is the deposition of particles subject to gravity@xcite , which underlies the formation of sandstones , ceramics and in general , the emergence of porous materials . modelling porous media@xcite in a realistic way is important for instance to understand permeability of soils for petroleum engineering or in material sciences . however , a model for grain deposition should incorporate two ingredients not always easy to combine . first , empirical studies@xcite show that the shape of grains significantly deviates from the sphere with a non - negligible polydispersity in their sizes . until recently most computational studies dealt with round particles for simplicity , in solving specific problems such as obtaining the optimal packing@xcite , constructing space - filling configurations@xcite or studying elongated shapes@xcite . some progress was recently made when addressing the influence of particle shape on the macroscopic observables of the agglomerate , focusing on ellipsoidal particles@xcite and even in more general shapes@xcite , without accounting for the size polydispersity . second , deposition processes are fundamentally sequential and therefore should be modelled through sequential procedures . the procedure introduced by donev _ et al_@xcite for jammed configurations , captures some of the general features observed in granular deposits , such as their density and average coordination number , although it is based on a non - sequential procedure . however , deposition under gravity does not generally lead to a jammed state . recently , some sequential procedures were already carried out in two - dimensional deposits of mono - disperse elongated particles@xcite .    in this paper we focus on three - dimensional polydisperse deposits of general ellipsoids and study how macroscopic properties of such deposits depend on the shape and size polydispersity of particles . we show that the density and coordination number of deposits behaves similarly to what is observed in jammed systems . however , in contrast with randomly constructed jammed systems , deposits become strongly anisotropic when the shape of the particles deviate from the sphere . furthermore , we find that the size polydispersity plays a minor role when compared to the shape .    to this end , we develop an efficient algorithm for the calculation of the interactions forces at all contact points within a three - dimensional deposits of polydisperse ellipsoids and use it in a molecular dynamics@xcite method to perform our simulations . further , we introduce two shape parameters which describe any ellipsoid , enabling a systematic study of the role of particle shape in granular deposits . the outline of the paper is as follows . in sec . [ sec : model ] we describe the algorithm for calculating the contact forces between colliding ellipsoids with arbitrary size and shape as well as the simulation setup . in sec . [ sec : controlparam ] the control parameters describing the polydispersity in the size and the shape of a general ellipsoid are introduced . in sec . [ sec : macroprop ] the dependency of macroscopic observables taken from the entire deposit are studied as functions of the control parameters . we focus on the density , coordination number and anisotropy of the deposits . discussion and conclusions are given in sec . [ sec : conclusions ] . when dealing with large systems , calculating the inter - particle forces is by far the most time consuming part of the computations . this is due to the fact that , in general , every particle can interact with every other one . in granular materials , however , the inter - particle forces are short range making techniques such as verlet list and linked - cell algorithm indispensable in reducing the computation time . these technique are employed to eliminate the redundant calculations for pairs which are too far apart to have any interaction@xcite .    in time - driven molecular dynamics ( md ) methods of granular materials the system is evolved for one time step allowing the particles to overlap . then normal repulsion forces on the particles are calculated as a function of the amount of overlapping , which is defined usually by a line segment . here we refer to this line segment as the contact vector ( see fig .  [ fig : contact ] ) .    the contact vector between two ellipsoids can be defined in more than one way . the common method , introduced by perram and co - workers@xcite and developed by others@xcite , is based on finding two points , one on the surface of each ellipsoid by minimizing a potential through a variational problem with one geometric constraint . the two obtained points can be used to define the contact vector . however , the contact vector defined this way is not necessarily normal to the surface of the particles , particularly when particles are strongly asymmetrical , leading to non - physical motions of the particles . an alternative method consists in finding two points on the surfaces of the ellipsoids at which the corresponding surface normals point exactly in opposite directions . consequently , the vector connecting these two points is normal to the surfaces of both ellipsoids and can be used as a well - defined contact vector . this method is , however , much more computationally expensive than the former one . see ref . @xcite for more details . to retain the computational efficiency , we use the former method and introduce an additional step for correcting the direction of the contact vector .     used to define the normal repulsive force between two overlapping ellipsoids @xmath0 and @xmath1 . the contact point is chosed to be @xmath2 after properly deriving points @xmath3 and @xmath4 at the surfaces of the ellipsoids ( see text ) . the full procedure is sketched in fig . [ fig : diagram].,scaledwidth=45.0% ]    the geometric potential of the ellipsoid @xmath0 is defined as : @xmath5 where @xmath6 and @xmath7 are the cartesian coordinates of a point on the surface and the centroid of the ellipsoid , respectively . the equation @xmath8 describes the surface of ellipsoid @xmath0 . considering a pair of ellipsoids @xmath9 and @xmath10 , two minimum points @xmath3 and @xmath4 are found on the surfaces of @xmath0 and @xmath10 by minimizing @xmath11 and @xmath12 , respectively . these points are shown schematically in fig.[fig : contact ] . to obtain the minima of @xmath13 subject to the constrain of being on the surface of ellipsoid @xmath1 we minimize the following auxiliary potential : @xmath14,\ ] ] where @xmath15 is the lagrange multiplier associated with the constrain . setting to zero the gradient of this function with respect to @xmath6 , @xmath16 we obtain an equation for _ optimum _ points @xmath17 for @xmath0 : @xmath18    to derive an equation for @xmath15 we left - multiply eq . ( [ eq : lagrange ] ) by @xmath19 and use the fact that @xmath17 lies on @xmath1 and that the gradients of the potential functions of the ellipsoids at _ minimum _ point are in opposite directions : @xmath20 equations ( [ eq : min ] ) and ( [ eq : lambda ] ) can be solved iteratively to find @xmath21 up to the desired precision . the second point @xmath22 can be found similarly . the contact point @xmath2 where a contact force is exerted on both ellipsoids is defined as the midpoint between @xmath3 and @xmath4 . the contact vector for ellipsoid @xmath0 is defined as @xmath23 and for ellipsoid @xmath1 , @xmath24 ( see fig .  [ fig : contact ] ) .    although eq . ( [ eq : min ] ) involves inversion of a @xmath25-matrix , finding the minima can be very efficient by choosing the initial points properly . good approximations for minima are the minima from the last time step or , if the ellipsoids were disjoint in the last time step , the closest points on the surfaces of the ellipsoids which are known from collision detection procedure described later in this section . these points are then used as initial points , reducing the number of iterations significantly .    the overlapping vector calculated in this way , however , is not necessarily normal to the ellipsoids at the contact point . in fact , the normal directions to @xmath0 and @xmath1 at points @xmath3 and @xmath4 deviate significantly from each other when the particles are strongly aspherical . to overcome this problem we introduce an additional step for correcting the contact vector , as described next . if the calculated contact vectors deviate more than a given amount from the gradients of the ellipsoid s potential at the contact point @xmath2 , a new direction is calculated by averaging the gradients at @xmath2 . then , two new points are calculated from the intersection of the line along this direction passing through the contact point . these new points are used to define a new contact vector . although it is not generally possible to define a vector which is along the gradients of both ellipsoids at the contact point , this approximation gives satisfactory results in terms of energy conservation . having described the procedure for calculating the contact of two colliding ellipsoids , we next need an efficient way to determine when two previously disjoint ellipsoids collide . to this end , we use the technique proposed by wang _ et al_@xcite which is described briefly as follows .    in homogeneous coordinates the equation for the surface of a general ellipsoid is given by @xmath26 where @xmath6 is any point on the surface of the ellipsoid and @xmath27 is the @xmath28-matrix of the affine transformation which transforms the unit sphere centered at the origin into the ellipsoid .        given two ellipsoids @xmath29 and @xmath30 , their characteristic equation is defined as @xmath31 , whose roots are the eigenvalues of matrix @xmath32 . furthermore , it can be shown @xcite that ( i ) @xmath33 has at least two negative roots , ( ii ) @xmath33 has two distinct positive roots if and only if the two ellipsoids are disjoint , ( iii ) @xmath33 has a positive double root if and only if the two ellipsoids are externally touching . therefore , by examining the eigenvalues of @xmath34 it can be determined if two ellipsoids collide . when the ellipsoids are disjoint , the corresponding eigenvectors give the coordinates of the four vertices @xmath35 , @xmath36 , of a tetrahedron which is self - polar for both ellipsoids @xmath0 and @xmath1@xcite , i.e.  @xmath37 . two of the vertices lie outside both ellipsoids ( external vertices ) , while the other two are contained inside the ellipsoids , one in each ( internal vertices ) . in this case , a set of separation planes can be defined using three points exterior to both ellipsoids , two of them being the external vertices and the other a point on the line segment connecting the internal vertices . we choose the third point as the midpoint between the surfaces of two ellipsoids on this line segment .    having the separation plane between the two disjoint ellipsoids , in the subsequent time step we check if the ellipsoids intersect that plane . if not , the ellipsoids are still disjoint and the procedure stops . this helps to eliminate the need for checking for collision between the ellipsoids which is more computationally expensive . only when any of the ellipsoids intersects the separation plane it is necessary to check again for the collision between the ellipsoids . if the ellipsoids turn out to be still disjoint their separation plane and closest points are updated to be used in the next time step . if the ellipsoids do collide the contact vector is calculated and saved to be used for the calculation of the contact force . figure [ fig : diagram ] shows the flow - chart of the complete procedure which is performed for all potentially overlapping pairs of ellipsoids . it is worth noting that , to further speed up the procedure , only the pairs whose spherical envelopes intersect are evaluated . the contact force can be decomposed to normal and tangential forces . the tangential forces such as static and dynamical frictions are generally derived from the normal force . for calculating the magnitude of the normal force we use hertzian model@xcite : @xmath38 with @xmath39 being the length of the contact vector between the particles and @xmath40 being the dissipative constant depending on the material viscosity @xcite . since we are only concerned with the properties of the deposit at rest , we choose @xmath40 sufficiently large for the deposit to relax rapidly . this relaxation is enhanced by including a tangential force @xmath41 , @xmath42 being the dynamic friction coefficient and @xmath43 the relative tangent velocity of the particles .    finally , the procedure described above is used in a md simulation together with standard methods , namely a prediction - correction integrator for solving the equations of motions , quaternions for describing the orientation of the particles and linked - cell algorithm to eliminate unnecessary collision detections . for details of the md methods see e.g.  ref . @xcite and references therein . the size and shape of the particles are the main parameters whose effects on the properties of the system are to be studied . in the following we explain how these are defined and chosen when preparing the samples .     for different values of @xmath44 and @xmath45 as defined in eq . ( [ vv ] ) . the inset shows the corresponding size distributions @xmath46 ( see text).,scaledwidth=45.0% ]    we define the size of an ellipsoidal particle as the radius of the sphere with the same volume , i.e.  @xmath47 with @xmath48 , @xmath49 and @xmath50 being the three semi - axis radii of the ellipsoid . we refer to a sample as mono - sized or monodisperse in size when all its constituting particles have the same volume , i.e.  the same size @xmath51 . to introduce size polydispersity we adopt the approach by voivret _ et al_@xcite : instead of directly choosing the size distribution @xmath46 of the particles , we consider the distribution of the total volume @xmath52 of all particles with size @xmath51 .     ellipsoids with shape parameters @xmath53 and @xmath54 a polydispersity of @xmath55 , that is , @xmath56 and @xmath57 ( see text ) . colors are arbitrarily chosen for better visualization.,scaledwidth=45.0% ]    we chose @xmath58 as @xmath59 where @xmath44 and @xmath45 determine the shape of the distribution and @xmath60 is the normalization factor . figure [ fig : vr - nr ] shows @xmath58 for different values of @xmath44 and @xmath45 and the corresponding size distributions @xmath46 ( inset ) . for @xmath61 this function takes a symmetric shape with a peak at @xmath62 . for the samples studied in this work we have chosen @xmath63 which corresponds to a distribution ( solid line in fig . [ fig : vr - nr ] ) very close to a truncated gaussian . the width of the distribution is controlled by the polydispersity parameter @xmath64 which is defined as @xcite : @xmath65 for @xmath66 the particles are monodisperse while @xmath67 corresponds to infinite polydispersity . the shape of an ellipsoid is characterized by two parameters here defined as @xmath68 and @xmath69 , with @xmath70 . for prolate ( @xmath71 ) and oblate ( @xmath72 ) ellipsoids , the shape can be fully characterized by the aspect ratio @xmath73 . the particular case @xmath74 corresponds to the sphere . these shape parameters together with the size @xmath75 fully specify any ellipsoid . similar to the size polydispersity , two additional parameters @xmath76 and @xmath77 can be similarly defined for the polydispersities in the shape of the particles . in this study we consider systems of mono - shaped particles , describing each deposit by @xmath78 , @xmath79 and size polydispersity @xmath64 .    in this work , all the deposits are generated by releasing particles with randomly chosen positions and orientations and letting them fall in an open box , under gravity along the vertical direction . the box is limited from below by a rigid wall ( ground ) and is periodic in @xmath80 and @xmath81 directions . in order to reduce the boundary effects the first few bottom layers of particles are not considered in our analysis . figure [ fig : deposit ] shows a deposit of ellipsoidal particles with @xmath53 , @xmath54 and @xmath55 .     of monodisperse deposits as a function of the shape parameters @xmath78 and @xmath79 . the density is approximately symmetric with respect to the diagonal @xmath82 where it attains a maximum . * ( b ) * the density as a function of shape along the diagonal @xmath82 ( solid line ) together with the case of prolate ellipsoids ( @xmath83 , dashed line ) and of oblate ellipsoids ( @xmath84 , dashed - dotted line ) . in the inset the density is plotted as a function of @xmath82 for different polydispersities @xmath64 . @xmath85 depending on the amount of the size polydispersity.,title="fig:",scaledwidth=45.0% ]   of monodisperse deposits as a function of the shape parameters @xmath78 and @xmath79 . the density is approximately symmetric with respect to the diagonal @xmath82 where it attains a maximum . * ( b ) * the density as a function of shape along the diagonal @xmath82 ( solid line ) together with the case of prolate ellipsoids ( @xmath83 , dashed line ) and of oblate ellipsoids ( @xmath84 , dashed - dotted line ) . in the inset the density is plotted as a function of @xmath82 for different polydispersities @xmath64 . @xmath85 depending on the amount of the size polydispersity.,title="fig:",scaledwidth=45.0% ]     and @xmath79 . here @xmath86 . , scaledwidth=45.0% ] the macroscopic properties considered in our study are the packing density , the average coordination number and the anisotropy . figure [ fig : density ] shows the packing density in deposits of @xmath86 mono - disperse ellipsoids as a function of the shape . the immediate realization is that the density is to a large extend symmetric with respect to the axis @xmath82 . when the shape deviates from the spherical shape , i.e.  when @xmath87 or @xmath88 , the density first increases rapidly , then it attains a maximum and finally decreases slowly . from fig . [ fig : density]a one sees that typically , along lines of constant @xmath78 ( or constant @xmath79 ) the maximum is observed around @xmath82 . further , along @xmath82 a global maximum is attained at @xmath89 , as can be observed from fig . [ fig : density]b ( solid line ) . this absolute maximum curiously approximates the value of the golden ratio @xmath90 , which from the definition of @xmath78 and @xmath79 corresponds to when @xmath91 . whether such an observation is merely a coincidence or not needs further investigations which is out of the scope of our manuscript . these results are qualitatively similar to those obtained in refs . @xcite where a slightly higher density is observed , due to their method for generating the packings : instead of being generated via deposition , their samples are prepared in order to be in jammed state using a generalized form of lubachevsky - stillinger algorithm @xcite . the effect of polydispersity is examined by comparing the result for three different values of @xmath64 as shown in the inset of fig .  [ fig : density]b . the results indicate that the density seems to be insensitive to low and moderate size polydispersities . figure [ fig : coordnum ] shows the average coordination number @xmath92 as a function of the shape of the particles . for spheroidal particles , @xmath92 increases gradually with the aspect ratios toward a maximum . for a general ellipsoid the coordination number attains a maximum higher than the maximum observed for spheroids ( @xmath84 or @xmath83 ) , but decreases to values comparable to those of spheroids for larger @xmath78 and @xmath79 . this can be understood considering the fact that for @xmath93 , one of the semi - axis radii is significantly larger than the other two and therefore the particle can be regarded as a prolate ellipsoid .      to study the anisotropy of the deposits we investigate the orientational order of the contact vectors and the principal directions of the ellipsoids by calculating their angular distributions . in general two angles are needed to specify the orientation of a vector in three dimensional space , namely one azimuthal angle and one angle with the vertical direction . our deposits are symmetric by construction in azimuthal direction . this is also observed in our results ( not shown here for the sake of brevity ) . therefore , we will only consider the angle with the positive vertical direction .     between the contact vectors and the positive @xmath94-semiaxis for different shape parameters @xmath78 and @xmath79 , for the mono - sized case ( @xmath66 ) . * ( b ) * histogram of the angle @xmath95 for for a given shapen ( @xmath96 , @xmath54 ) but different size polidispersity . here @xmath85 depending on the amount of the size polydispersity . , title="fig:",scaledwidth=45.0% ]   between the contact vectors and the positive @xmath94-semiaxis for different shape parameters @xmath78 and @xmath79 , for the mono - sized case ( @xmath66 ) . * ( b ) * histogram of the angle @xmath95 for for a given shapen ( @xmath96 , @xmath54 ) but different size polidispersity . here @xmath85 depending on the amount of the size polydispersity . , title="fig:",scaledwidth=45.0% ]     between the largest semi - axis @xmath48 of each particle and the positive @xmath94-semi - axis . each curves corresponds to particular values of the shape parameters . here @xmath86 . , scaledwidth=45.0% ]    figure [ fig : contact - hist]a shows the results for deposits for different values of the shape parameters @xmath78 , @xmath79 where @xmath66 . the vertical axis , @xmath97 , represents the fraction of number of vectors within the solid angle swept by a conic shell whose axis is parallel to @xmath94 and its vertex angle is @xmath98 . note that the curves are symmetric with respect to @xmath99 because for each contact normal of a particle there is a contact normal with opposite direction from the other particle to which it is in contact . the results show that for deposits of spheres the distribution @xmath100 is independent of @xmath95 , and consequently the configuration of the contacts is isotropic and has no preferred direction . however , as the shape of the particles deviates from a sphere the number of side contacts between particles , i.e.  the number of contacts on the horizontal plane , decreases and the contact normals line up along the vertical direction . this is seen in fig . [ fig : contact - hist]a which shows for each shape there is a minimum at @xmath99 which corresponds to the number of contacts on the @xmath101-plane , while having a maximum around @xmath102 ( vertical direction ) . the minimum at @xmath99 together with the maximum at @xmath102 get more pronounced when increasing @xmath78 and @xmath79 . this means that for larger values of @xmath78 and @xmath79 which correspond to higher asphericities the particles tend to lie horizontally , for minimizing potential energy during deposition . consequently most of the contacts will occur with the particles beneath and above . this can have dramatic effect on the force chains and the response of under shear and compression . such a study is , however , beyond the scope the current work . figure  [ fig : contact - hist]a also shows two crossing points at a value slightly higher than @xmath103 , indicating this direction as an invariant directions while varying the shape of the depositing particles .    figure [ fig : contact - hist]b shows the comparison of the angular distribution of contact normals for deposits with @xmath104 , with different polydispersities @xmath64 . it is clear that the polydispersity has a minor effect .    to study the anisotropy of the orientation of the particles we do a similar analysis for the principal directions of each ellipsoid . figure [ fig : orient - hist ] shows the results for the largest semi - axis @xmath105 for the deposits . one can see that as the shape of the particles deviates from a sphere the largest semi - axis @xmath105 also tends to lie on the @xmath101-plane . an invariant orientation is also observed , now for @xmath106 . the polydispersities in the sizes of the particles has again no significant effects on the orientation anisotropy ( not shown ) . in this work , we studied systematically the static properties of deposits composed of ellipsoidal particles of different size and shape , focusing on their density and anisotropy . to this end , we developed an efficient molecular dynamics method for simulation of ellipsoidal particles which we used to generate deposits of such particles under gravity . for mono - sized particles , the density increases rapidly as the shape of the particles deviates from a sphere , reaching for particular values of the shape parameters a maximum significantly higher than the one observed in random close packing of spheres . with increasing further the shape parameters , that is , deviating further for a sphere density to decrease slowly . the densities observed in the deposits are generally lower than the ones obtained through procedures for optimal packings . thus , as expected , deposition is not an ideal process for achieving the largest densities . typically , porous media and media formed through deposition processes present also densities much lower than the optimal one . such observations point toward the possibility for using deposition of ellipsoids as a proper procedure for modeling porous media . this is the subject of future works .    in the case of polydisperse particles the density takes higher values with qualitatively similar behavior . a rather unexpected result is that the polydispersity in the sizes of the particles has almost no effect on the anisotropic behavior of the system . we also study the anisotropy of such deposits and show that as the shape of the particles deviates from a sphere stronger anisotropies are observed . such anisotropies appear due to the gravity and are absent in the systems without gravity . the authors thank bibhu biswal and jens harting for useful discussions and the _ fundao para a cincia e a tecnologia _ for financial support , under the fellowship _ sfrh / bpd/48974/2008 _ and contract _
Abstract:

Answer: we present a systematic overview of granular deposits composed of ellipsoidal particles with different particle shapes and size polydispersities . we study the density and anisotropy of such deposits as functions of size polydispersity and two shape parameters that fully describe the shape of a general ellipsoid . our results show that , while shape influences significantly the macroscopic properties of the deposits , polydispersity plays apparently a secondary role . the density attains a maximum for a particular family of non - symmetrical ellipsoids , larger than the density observed for prolate or oblate ellipsoids . as for anisotropy measures , the contact forces show are increasingly preferred along the vertical direction as the shape of the particles deviates for a sphere . the deposits are constructed by means of an efficient molecular dynamics method , where the contact forces are efficiently and accurately computed . the main results are discussed in the light of applications for porous media models and sedimentation processes .


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: there are mainly three kinds of electron emission processes from metal surface , they are : ( i ) thermionic emission , ( ii ) photoelectric emission and ( iii ) cold emission or field emission . thermionic emission is the external thermal energy induced flow of electrons through a metallic surface barrier into vacuum . this occurs because the thermal energy given to the electrons overcome the confining barrier potential . therefore , to have a thermo - electric current , the thermal energy or equivalently the temperature should have a lower critical value . this lower limit strongly depends on the work function of the material , which is the minimum energy needed to remove an electron with zero kinetic energy from the material surface . the photoelectric effect is a phenomenon in which electrons are emitted from matter as a result of absorption of electromagnetic energy in the form of quanta of the radiation incident on their surfaces . analogous to the thermionic emission , in the case of photo - emission , the incident electromagnetic waves or photons must have a lower critical frequency , below which there will be no photo - emission . this lower limit also strongly depends on the work function of the material . finally the field emission ( also known as cold emission ) is an electron emission process induced by strong external electric fields at zero or extremely low temperature . field emission can happen from solid and liquid surfaces , or from individual atoms . it has been noticed that the field emission from metals occurs in presence of high electric field : the gradients are typically higher than 1000 volts per micron and strongly dependent upon the work function . in modern science and technology , there are a lot of applications of field emissions , which includes construction of high - resolution and high luminosity electron microscopes and also used in modern nano science and technology @xcite . further , the thermal emission from metallic surface can be explained by classical boltzmann statistics , with only quantum correction needed is the inclusion of fermi energy . on the other hand , the photo - emission of electrons can not be explained unless the quantum nature of electromagnetic wave is considered , i.e. , old quantum mechanics or quantum theory is needed where photon concept of light is introduced . of course , the concept of work function associated with any kind of electron emission process is a quantum mechanical phenomenon . the field emission process on the other hand can only be explained by quantum tunneling of electrons . which has no counter classical explanation . however , for general type of surface barrier , this purely quantum mechanical problem can not be solved exactly , a semi - classical approximation , known as wkb is needed to get tunneling coefficients . now because of quantum fluctuation , electrons from the sea of conduction electrons ( electron gas ) always try to tunnel out through the metallic surface ( surface barrier ) . however , as soon as an electron comes out , it induces an image charge on the metal surface , which pulls it back and does not allow it to move far away in the atomic scale . but if a strong attractive electrostatic field is applied near the metallic surface , then depending on the fermi energy of electrons and the height of the surface barrier , it may overcome the effect of image charge and the electrons get liberated . since the external strong electric field is causing such emission and does not depend on the thermal properties of the metal , even the metal can be at zero temperature , it is called field emission or cold emission .    the theory of field emission from bulk metals was first proposed by fowler and nordheim in an epoch making paper in the proceedings of royal society of london in the year 1928 @xcite . fowler - nordheim tunneling is the wave - mechanical tunneling of electrons through a triangular type barrier produced at the surface of an electron conductor by applying a very high electric field . fowler nordheim tunneling current is found to have important implication in various scientific and technological applications , starting from cold cathode emission to vlsi . cold emission or field emission processes not only have significance in the terrestrial laboratories , but is found to be equally important in the electron emission processes from cold and compact stellar objects . for an example , in the case of a rotating neutron star , the existing large magnetic field , @xmath0 g for the conventional radio pulsars and @xmath1 g for the magnetars , produces a strong electric field at the poles , approximately given by @xmath2volt @xmath3 , parallel to @xmath4 at the poles @xcite , here @xmath5 is the time period of the neutron star in second . at the proximity of polar region the potential difference changes almost linearly with distance from the polar cap , which is a region very close to the magnetic pole . the repulsive surface barrier in combination with this attractive potential , forms a triangular type barrier at the poles . now , the study of plasma formation in a pulsar magnetosphere is a quite old but still an unresolved astrophysical issue . in the formation of magneto - spheric plasma , it is generally assumed that there must be an initial high energy electron flux from the magnetized neutron stars . at the poles of a neutron star the emitted charged particles flow only along the magnetic field lines . the flow of high energy electrons along the direction of magnetic field lines and their penetration through the light cylinder is pictured with the current carrying conductors . naturally , if the conductor is broken near the pulsar surface the entire potential difference will be developed across a thin gap , called polar gap . this is of course based on the assumption that above a critical height from the polar gap , because of high electrical conductivity of the plasma , the electric field @xmath6 , parallel with the magnetic field near the poles is quenched . further , a steady acceleration of electrons originating at the polar region of neutron stars , traveling along the field lines , will produce magnetically convertible curvature @xmath7-rays . if these curvature @xmath7-ray photons have energies @xmath8 ( with @xmath9 the electron rest mass and @xmath10 the velocity of light ) , then pairs of @xmath11 will be produced in enormous amount with very high efficiency near the polar gap . these produced @xmath11 pairs form what is known as the magneto - spheric plasma @xcite . the cold emission , therefore plays a significant role in magneto - spheric plasma formation . in turn , the motion of charged particles in the magnetosphere in presence of strong magnetic field causes pulsar emission in the form of synchrotron radiation . therefore the cold emission process indirectly affects the intensity of synchrotron radiation . although the field electron emission is a paradigm example of what we call the quantum mechanical tunneling , unfortunately , it is also a paradigm example of the intense mathematical difficulties associated with this tunneling process with general type barrier structure . exactly solvable models with simple type tunneling barrier lead to equations @xcite that underestimates the emission current density by a factor of 1000 or more . if a more realistic type barrier model is used by inserting an exact potential barrier in the simplest form of the schr@xmath12dinger equation , then a complicated mathematical problem with very little physical important arises over the resulting differential equation : it is known to be mathematically impossible in principle to solve this equation exactly in terms of the usual functions of mathematical physics , or in any simple way . to get even an approximate solution , it is necessary to use special approximate methods known in physics as semi - classical or wkb approximation . moreover , to the best of our knowledge , a relativistic version of cold emission model in presence of strong magnetic field , relevant for electron emission from the polar region of strongly magnetized neutron stars , even with a simple type potential barrier has not been reported in the literature . it has already been mentioned that the emission of high energy electrons and thereby formation of magnetosphere is caused by strong electric field at the polar region and this strong electric field at the poles is produced by the magnetic field of rotating neutron stars . we have organized the article in the following manner : in section-2 we have studied electron emission with a scalar type potential barrier . in section-3 we have studied the same thing replacing the scalar potential by a vector type potential barrier . in the last section we have given the conclusion and future prospect of this work . in the original paper by fowler and nordheim on cold emission , the conventional non - relativistic quantum mechanical equations , i.e. , the schr@xmath13dinger equation was used @xcite . in the relativistic version we shall use the dirac equation and introduce the triangular form of barrier potential at the surface of the system using the conventional form of relativistic quantum hadro - dynamics . in the simplified picture of quantum hadro - dynamic model the interaction can either be of scalar type or vector type or both . in this model the conventional form of dirac equation with scalar and vector type interactions is given by @xcite @xmath14\psi = 0\ ] ] where @xmath15 and @xmath16 are the vector and scalar coupling constants respectively , @xmath17 and @xmath18 are respectively the vector and the scalar fields and @xmath19s are the dirac @xmath7-matrices . in presence of strong neutron star crustal magnetic field , electrons can travel almost along the field direction . if the magnetic field strength @xmath20 , where @xmath21 is the typical strength of magnetic field at which the electron cyclotron quantum exceeds the corresponding rest mass energy , then in the relativistic scenario the landau levels of the electrons will be populated . in this situation , the electron momentum in the plane orthogonal to the direction of magnetic field , gets quantized and is given by @xmath22 , where @xmath23 are the landau quantum numbers . assuming that the magnetic field is along @xmath24-direction the corresponding electron momentum component changes continuously from @xmath25 to @xmath26 . it can be shown that the critical strength @xmath27 g . in this situation the dirac equation reduces to effectively one - dimensional in nature . further , we replace the scalar and vector part of the fields by the non - relativistic potentials at the interface between matter and vacuum with triangular in shapes . in the case of vector potential , only the zeroth part is retained whereas its spacial part vanishes because of rotational symmetry within the matter . in one - dimension , for the sake of convenience we use the form of dirac equation originally proposed by dirac , where @xmath7-matrices are replaced by @xmath28 and @xmath29 matrices . the triangular shape potential barrier is represented by the form @xmath30 , where following the notation of fowler and nordheim , @xmath31 is the surface barrier and @xmath6 ( absorbing the electron charge , we express @xmath32 by @xmath6 ) is the electric field at the polar region . the quantity @xmath6 acts as the driving force for electron emission . let us first assume that the potential barrier is of scalar type . then the dirac equation in one - dimension is given by @xmath33\psi = e\psi\ ] ] where @xmath34 is the energy eigen value . following the representation @xcite @xmath35 and defining @xmath36 as the effective electron mass , then with the spinor notation @xmath37 after eliminating @xmath38 from the eqns.(5)-(6 ) , we have @xmath39 to solve this equation analytically , we make the following transformation of coordinate : @xmath40 then we have @xmath41 this equation may be re - arranged in the following form @xmath42 which is the well known hermite differential equation , with @xmath43 , which gives @xmath44 and the solution is given by @xmath45 where @xmath46 is the normalization constant and @xmath47 is the well known hermite polynomial of order @xmath48 ( which is not necessarily integer ) . it is trivial to show that @xmath38 will also have the same kind of solution . this spinor solution ( @xmath49 or @xmath38 ) is for the electron which has already been liberated through the surface into vacuum under the influence of electric field @xmath6 .    now the one - dimensional form of free dirac equation for the electrons bound within the metal body ( in this case within the metallic crustal matter ) , i.e. , before emission from the system is given by @xmath50 then following the same technique as above , we have @xmath51 where @xmath52 ( exactly same type of equation can be obtained for @xmath38 ) . now following the notation of fowler and nordheim , we express the solution for free electrons within the system , confined by the surface barrier @xmath53 , in the form @xmath54= u_{ii}~~{\rm{(say)}}\ ] ] where @xmath55 is the probability amplitude for electrons moving along the positive direction of @xmath56 , whereas @xmath57 is the corresponding quantity for left moving waves . assuming the interface between the neutron star matter ( which in this case is the dense metallic iron crystal of neutron star crustal matter ) and vacuum is at @xmath58 , the wave function and their derivatives must be continuous at @xmath58 , i.e. , @xmath59 using the relation @xmath60 , we have @xmath61 and @xmath62\ ] ] where @xmath63 . these two conditions ( eqns.(16)-(17 ) ) may be rearranged in the following form @xmath64 where @xmath65 and @xmath66 are two real quantities . hence it is straight forward to verify that the transmission coefficient , given by @xmath67 vanishes exactly . in this section we therefore conclude that if the barrier in combination with the external electro - static driving force behaves like a scalar potential , electrons can not tunnel through the surface barrier whatever be their kinetic energies and the strength of external electric field . identical result will be obtained if @xmath38 is used instead of @xmath49 . this result has some analogy with the confinement of quarks inside hadronic bag with scalar type linear potential . in presence of a vector type potential barrier at the interface with the same kind of triangular shape as mentioned in the previous section , we have the conventional form of dirac equation @xmath68 decomposing the spinor @xmath69 into @xmath49 and @xmath38 as before and using the same representation for @xmath28 and @xmath29 , we have @xmath70 eliminating @xmath38 from these two equations ( eqns.(21)-(22 ) ) we have @xmath71+m^2u=(e - c+fx)^2u\ ] ] substituting @xmath72 , we have @xmath73 to obtain an analytical solution for @xmath49 , let us consider the non - linear transformation @xmath74 . then the above equation becomes @xmath75 to obtain a closed analytical solution , we neglect for the sake of simplicity the @xmath76 term . then we have @xmath77 where @xmath78 and @xmath79 . this is the standard differential equation for the well - known coulomb function with zero angular momentum @xcite . the solution is given by @xmath80=u_i ~~{\rm{(say)}}\ ] ] where @xmath81 .    ' '' ''    : @xmath82 _ where _ @xmath83    ' '' ''    to obtain tunneling co - efficient , we again consider the continuity of @xmath49 and its derivative at the interface ( i.e. , at @xmath58 ) . which means @xmath59 hence we get @xmath84 where @xmath85\ ] ] is the phase factor , and @xmath86\ ] ] where , as before @xmath87 . hence , from the definition of tunneling coefficient , we have @xmath88 ^ 2}\end{aligned}\ ] ] substituting @xmath89 , we finally get @xmath90 ^ 2}\ ] ] obviously , in this case , unlike the scalar type surface barrier , the transmission co - efficient is non - zero . further , unlike fowler nordheim non - relativistic model , here the tunneling coefficient dose not fall exponentially . in the above expression for tunneling coefficient , as before @xmath91 and @xmath92 , with @xmath93 and @xmath94 are the landau quantum numbers for the electrons . in presence of magnetic field , it has been observed that the work function increases with the increasing field strength . a numerically fitted form for magnetic field dependence is given by @xmath95 in ev , where @xmath96 is a real constant parameter @xmath97 @xcite . one can easily verify that the tunneling coefficient will vanish for @xmath98 and extremely large values for electron work function . to obtain electron cold emission current in the relativistic scenario in presence of strong magnetic field , we consider the motion of relativistic electrons along the direction of magnetic field and the system is assumed to be at zero temperature . it is well known that in presence of strong quantizing magnetic field @xmath1 g , which is the recently observed surface values of some strange type neutron stars , known as magnetars @xcite , the motion of the electrons at the crustal region of a neutron star is effectively one dimensional in nature ; motion is restricted along the direction of flux lines . in this case the density of states is given by @xcite @xmath100 in this expression we have taken planck constant @xmath101 , velocity of light @xmath102 and boltzmann constant @xmath103 ; known as natural unit . here @xmath48 is the landau quantum number and the factor @xmath104 takes care of singly degenerate zeroth landau level and doubly degenerate other levels with @xmath105 . the upper limit of landau quantum number is infinity at finite temperature , whereas at zero temperature it is given by @xmath106\ ] ] here @xmath107 is the electron chemical potential and @xmath108 $ ] indicates that @xmath109 is an integer but less than the actual number . now the number of electrons incident normally on unit area at the interface between matter and vacuum is given by @xmath110 + 1}\ ] ] where @xmath111 is the electron energy at @xmath48th landau level and @xmath112 is the lower limit of @xmath24 component of electron momentum . unlike thermionic emission and photo - emission , the cold emission is a purely quantum mechanical phenomenon ; electrons tunnel through the surface barrier and its emission in free space is driven by the external electric field . since the emission process is controlled by the tunneling probability of electrons through surface barrier , the actual expression for cold emission current will therefore be @xmath113 + 1}\times t_r\ ] ] in the case of neutron star , since the temperature is low enough compared to the chemical potential of electrons , we put @xmath114 in the fermi distribution in eqn.(38 ) , then we have after replacing the integration variables from @xmath115 to @xmath116 and omitting the index @xmath48 ( keeping in mind that @xmath117 ) in electron energy @xmath118 we have obtained cold emission current by evaluating numerically energy integral and the sum over landau quantum number . however , we have noticed that at the crust of magnetars the landau quantum number becomes exactly zero , i.e. , the electrons occupy only their zeroth landau level . in fig.(1 ) we have shown the variation of cold current with electric field and magnetic field by the solid curve and the dashed curve respectively . both the curves indicate the initial increase of cold current . since the work function increases with increasing magnetic field and electric field is generated by the rotating magnetic field of neutron stars , the cold current decreases abruptly after some magnetic field / electric field strength . in this article we have studied cold emission of electrons from strongly magnetized neutron star polar region using relativistic quantum mechanics . to get an analytical solution , we have considered triangular type potential barrier as was used by fowler and nordheim . we have noticed that a scalar type potential barrier does not allow electrons to tunnel out the interface , whatever be the strength of external electric field . with the vector type potential barrier , we have obtained an analytical solution for the tunneling coefficient . unlike the original work of fowler and nordheim the tunneling coefficient does not follow any exponential law . however , the variation of cold current with the electric field strength can be fitted numerically using a @xmath119-minimization program . the numerically fitted functional form is given by @xmath120 where @xmath121 and @xmath122 . the variation shows an exponential decay form is because of approximate linear dependence of electric field with the magnetic field strength of the rotatiting neutron star . we have obtained an analytical expression for the cold emission current which changes with the strength of magnetic field also . we found that the cold current initially increases almost linearly with the external electric field and finally when the work function becomes large enough because of high magnetic field strength , it drops abruptly and finally goes to zero . since the emission current becomes extremely small in presence of ultra - strong magnetic field , the magnetosphere of the magnetars will be of very low in charge density . as a consequence there will be very weak synchrotron emission in the radio wave band . 99 shi - dong liang and lu chen , phys . , 101 , 027602 ( 2008 ) . r.h . fowler and dr . l. nordheim , proceedings of the royal society of london 119 , 173 , ( 1928 ) . t.e . stern , b.s . gossling and r.h . fowler , proceedings of the royal society of london , 124 , 699 , ( 1929 ) . k. l. jensen , j. vac . technol . , b13 , 516 , ( 1995 ) ; see also r. g. forbes and j. h. b. deane , proceedings of the royal society of london , 463 , 2907 , ( 2007 ) . shapiro and s.a . teukolsky , black holes , white dwarfs , and neutron stars- the physics of compact objects , john wiley & sons , new york , 288 , ( 1983 ) . a. jessner , h. lesch and t. kunzl , apj , 547 , 959 , ( 2001 ) . m.a . ruderman and p.g . sutherland , apj , 196 , 51 , ( 1975 ) . see also , d.a . diver , a.a . da costa , e.w . laing , c.r . stark and l.f.a . teodoro , astro - ph/0909.3581 . s.l . shapiro and s.a . teukolsky , black holes , white dwarfs and neutron stars , john wiley and sons , new york , ( 1983 ) . michel , rev . phys . , 54 , 1 , ( 1982 ) michel , advances in space research , 33 , 542 , ( 2004 ) . a.k . harding and d. lai , rep . phys . , 69 , 2631 , ( 2006 ) . m. ruderman , phys . letts . , 27 , 1306 , ( 1971 ) . n. k. glendenning , compact stars ( second edition ) , springer , 162 , ( 2000 ) . r. jackiw and c. rebbi , phys . rev . , d13 , 3398 , ( 1976 ) ( see also r. s. bhalerao and b. ram , am . j. phys . , 69 , 817 , ( 2001 ) ) . m. abramowitz and i. a. stegun , handbook of mathematical functions , dover publications , inc . , new york , 538 , ( 1970 ) . s. mereghetti , astro - ph/0904.4880 . z.medin and d.lai , mnras , 382 , 1833 , ( 2007 ) ; z.medin and d.lai , phys . a74 , 062507 , ( 2006 ) ; a.harding and d.lai , astro - ph/0606674 . a. ghosh and s. chakrabarty ( submitted ) . k. huang , statistical mechanics , wiley estern private ltd . , new delhi , 237 , ( 1963 ) . a. ghosh and s. chakrabarty ( to be communicated ) .
Abstract:

Answer: a relativistic version of cold emission of electrons in presence of strong magnetic field , relevant for strongly magnetized neutron stars is obtained . it is found that in this scenario , a scalar type potential barrier does not allow quantum tunneling through the surface . whereas , in presence of a vector type surface barrier , the probability of electron emission is much larger compared to the original fowler - nordheim cold emission of electrons . it is found that the relativistic version in presence of strong magnetic field does not follow exponential decay .


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: quantum heat engines ( qhe ) have recently attracted great interest since it was originally proposed in @xcite . like their classical counterparts , these engines employ a quantum system as the working substance and run cyclically between two heat baths that are at thermal equilibrium at different temperatures @xmath0 and @xmath1 ( @xmath2 ) . it absorbs @xmath3 heat from the hot bath at a temperature @xmath4 , performs certain mechanical work @xmath5 , and rejects heat @xmath6 to the bath at a temperature @xmath7 . due to the quantum property of the working substance , the qhes can give rise to large work efficiency @xmath8 beyond the limit set for a classical heat engine ( viz . , for a classical carnot engine , such limit is @xmath9 ) . scully and his coworkers have shown that it is possible to extract work from a qhe using quantum coherence @xcite or quantum negentropy @xcite , with a large work efficiency , compared to their classical counterparts .    several proposals for implementing quantum heat engines have been made using spin systems @xcite , harmonic oscillators @xcite , and multi - level systems @xcite . in @xcite , an ensemble of non - interacting two - level systems , driven by a local external field and subject to frictional force , is coupled to the baths for _ finite _ time and the optimized time - scale of engine operation has been derived . in @xcite , the interacting spins have been considered and the adiabatic processes have been performed by changing the local magnetic fields that drive the spins . correspondence of the efficiency of such coupled - spin system to the entanglement has been studied for heisenberg interaction @xcite and dzyaloshinski - moriya interaction @xcite . comparisons between quantum versions of carnot engine and otto engine have been studied in @xcite for a two - level system as well as for harmonic oscillators . it is further shown @xcite that photons in an optical cavity can also be used as a working substance for a qhe in which a coherent cluster of atoms serve the role of ` fuel ' through their superradiance . however , in all these proposals , the system undergoes a heat cycle , during which it interacts with a hot bath and a cold bath ( both modelled as a classical system ) in an alternative fashion . it is assumed that one can switch off the interaction with the bath during a certain stage of the heat - cycle , the so - called reciprocating cycle . for a quantum system , however , such an interaction is always on @xcite and it is experimentally challenging to turn it off or on during the heat cycle . dynamical decoupling from the bath by applying a certain pulse sequences @xcite to the system may serve the purpose . but , such a technological requirement can be overwhelming so as to mimic the classical heat cycle in a quantum system .    in this paper , we propose an experimentally feasible model of a single - ion qhe , that works like a reciprocating heat cycle , but maintaining the always - on interaction with the bath , like in the continuous heat cycles @xcite . the electronic degree of freedom of the ion is chosen as the working substance s , while its vibrational degree of freedom plays the role of the cold bath , that also interacts with the ion throughout the entire cycle . the thermal environment here acts as the hot bath . we show that while the hot bath thermalizes the state of s , the heat transfer to the cold bath can be achieved by a projective measurement of the states of s. practically speaking , in our model , only a single thermal bath is considered , while the vibrational mode is only an ancillary system , mimicking a ` cold ' bath . we emphasize that the assumption , that interaction with none of the baths is switched off , makes our model easily feasible in experiments . we show that the efficiency of this engine can be made close to unity by manipulating a local magnetic field ( that works as a  piston " ) applied to the ion . note that the proposal of implementing qhe using a single ion exists @xcite , in which the frequency of the linear paul trap is changed during the isentropic processes of the heat engine , while during isochoric processes , the system is weakly coupled to the hot or cold bath to achieve the thermal equilibrium at the bath temperature . this model clearly requires alternative coupling to the baths . this proposal has been implemented using @xmath10ca@xmath11 ion with an efficiency 1.9% at the maximum power limit @xcite . on the contrary , in our model , one only requires to change the driving magnetic field adiabatically and a projective measurement of the electronic states , both of which can be routinely achieved in a trapped - ion set up @xcite . the paper is organized as follows . in sec . ii , we discuss the qhe model based on a single trapped ion . we present all the relevant hamiltonians and the achievable efficiency in this qhe . in sec . iii , we describe all the required stages of the heat cycles . we also present the conditions that are relevant to successfully implement these stages . we conclude the paper in sec . in this paper , we focus on implementing a quantum otto engine ( qoe ) . this qoe cycle consists of four stages : a ) in ignition stroke ( an isochoric process ) the system s gets thermalized by absorbing @xmath12 heat from a hot bath at a temperature @xmath0 . b ) in the next stroke ( the ` expansion ' stroke ) , the system undergoes an adiabatic process , thereby maintaining the thermal equilibrium , such that any kind of exchange of heat with the bath is inhibited . during this process , an external work @xmath5 is done by the system amounting to an adiabatic change of a local driving field ( the ` piston ' ) . c ) in the following stroke , the ` exhaust ' stroke ( an isochoric process ) , the system releases @xmath13 energy to a cold bath , at an equilibrium temperature @xmath1 . d ) in the last stroke ( the ` compression ' stroke ) , the system initializes itself through an adiabatic evolution , that also initializes the local driving field . in the following , we will describe how to implement all these strokes using a trapped ion . a single trapped ion in lamb - dicke limit can be considered as a two - qubit coupled system [ see fig . 1 ] , in which the ion is confined in its two lowest lying internal states @xmath14 and @xmath15 ( represented by the relevant pauli matrices @xmath16 ) and the lowest lying vibrational states @xmath17 and @xmath18 ( such that @xmath19 vanishes , where @xmath20 is the creation operator of the vibrational mode ) . the hamiltonian that describes the interaction between the internal and motional states of the ion can thus be written as ( in unit of planck s constant @xmath21 ) @xmath22 where @xmath23 here , @xmath24 is the hamiltonian for the internal states of the ion , @xmath25 represents the energy of the phonons , relevant to the vibrational degree of freedom , and @xmath26 defines the interaction between the internal and the vibrational degrees of freedom of the ion . we consider that the internal states are driven by a local electric field with rabi frequency @xmath27 and a magnetic field of strength @xmath28 , applied along the quantization axis . the vibrational frequency of the ion is chosen as @xmath29 and @xmath30 represents the interaction between the internal states and the motional states of the ion .    in context of the qoe , we consider the internal degree of freedom as the working substance s. the eigenvalues of the relevant hamiltonian @xmath24 are given by @xmath31 , with the respective eigenstates @xmath32\;,\\ \label{npm}n\mp & = & \sqrt{2\left[g^2+b^2\mp b\sqrt{g^2+b^2}\right]}\;.\end{aligned}\ ] ]      following kieu @xcite , we can now calculate the efficiency of a qoe based on the above two - level system s. the average energy of this system can be written as @xmath33 ( @xmath34 ) , where @xmath35 is the occupation probability of the energy eigenstate @xmath36 . comparing @xmath37 with the first law of classical thermodynamics , @xmath38 ( @xmath39 defines non - exact differential , see e.g. , @xcite ) , one can identify in the quantum regime the infinitesimal heat transfer as @xmath40 and the infinitesimal work done as @xmath41 . clearly , the heat transfer is associated with the change in the probability distribution for various eigenstates , while the work done is related to change in the energy eigenvalues , keeping the probability distribution the same . for the two - level system s , the heat absorbed @xmath12 during the ignition stroke can thus be written as @xmath42\;,\ ] ] where @xmath43 is the @xmath44th energy eigenvalue while the system interacts with the hot bath at the temperature @xmath0 , @xmath45 is the initial probability of the @xmath44th energy eigenstates and @xmath46 is that after thermalization . in a similar way , the heat released @xmath13 during the exhaust stroke can be written as @xmath47\;,\ ] ] where @xmath48 is the @xmath44th energy eigenvalue while the system interacts with the cold bath at the temperature @xmath1 . note that during the adiabatic processes , the probabilities @xmath46 and @xmath45 do not change . however the energy eigenvalues @xmath49 change , due to the adiabatic change in the local driving field . the work done during the heat cycle can thus be written as @xmath50\;,\ ] ] leading to the following expression of the work efficiency @xmath51 : @xmath52}{\sum_{n=1}^2e_n^h\left[p_n(t_h)-p_n(t_l)\right]}\;,\nonumber\\ & = & \frac{(e_1^h - e_2^h)-(e_1^l - e_2^l)}{e_1^h - e_2^h}\nonumber\\ & = & 1-\frac{e_1^l - e_2^l}{e_1^h - e_2^h}\end{aligned}\ ] ] where we have used the fact @xmath53 . if the local driving fields are changed from @xmath54 and @xmath55 to @xmath56 and @xmath57 , respectively , during the expansion stroke , the efficiency takes the following form : @xmath58 this suggests that in the limit of @xmath59 , one can achieve the efficiency close to unity . in the following , we assume that the local electric field that drives the ion is kept constant at @xmath60 . further , the adiabatic process is performed by adiabatic changes of the magnetic field from @xmath61 to @xmath62 during the expansion stroke and from @xmath63 to @xmath61 during the compression stroke . such a control of the efficiency by a local magnetic field has also been proposed in @xcite . as outlined before , we consider the internal degree of freedom as the working substance s. on the other hand , in the lamb - dicke limit , the ionic motion is confined to its two lowest lying states , while the higher excited states are not populated . in this way , the motional states can be considered as the relevant two - level cold bath , such that average phonon number in the vibrational degree of freedom @xmath64 . note that the system s always interacts with this effective cold bath through the hamiltonian @xmath26 . the thermal environment at an equilibrium temperature @xmath0 also interacts with the system s and the ionic motion . the qoe consists of four different strokes - two isochoric strokes and two adiabatic strokes . in the following , we describe how to implement all these strokes with the system s and the two baths identified as above .      during this isochoric process ( @xmath65 , 1 ) , the ion interacts with the hot bath and gets thermalized . this leads to the following mixed states of the system s and the phonons : @xmath66 where @xmath67 is the occupation probability of the @xmath68th eigenstate @xmath69 of the total hamiltonian @xmath70 for the ( s+phonon ) system . to identify the eigenvalue @xmath71 , we rewrite the hamiltonian @xmath70 in the joint basis @xmath72 of the internal states and the motional states in the following matrix form : @xmath73 the eigenvalues of this hamiltonian are given by @xmath74 where @xmath75    here @xmath69 can be written in terms of the joint basis of the internal and motional states as @xmath76 therefore , the reduced density matrix of the system s can be obtained by taking partial trace over the phonon states as @xmath77 where @xmath78 $ ] is a normalization constant . this can be written in terms of the energy eigenstates @xmath36 of the system hamiltonian through the inverse transformation of the eq . ( [ en ] ) : @xmath79    using ( [ rhos ] ) and ( [ geeigen ] ) , we can have the occupation probability @xmath80 of the eigenstate @xmath36 as @xmath81 and @xmath82    assuming that the system s is initially prepared in the state @xmath14 , the heat absorbed @xmath12 during this stroke can be calculated using eq . ( [ qh ] ) as @xmath83    it must be reminded that during this stroke , the magnetic field is kept constant at @xmath54 , leading to the eigenvalues @xmath43 of the system hamiltonian @xmath24 to remain constant .      during this stroke ( @xmath84 , see fig . 1 ) , the magnetic field is adiabatically changed from @xmath61 to @xmath63 . the system still maintains thermal equilibrium with the hot bath . this means that the occupation probabilities of the two eigenstates @xmath85 do not change ; however , the corresponding eigenvalues @xmath86 change to the values @xmath87 . this leads to the following work done during this stroke : @xmath88 where @xmath89 are the eigenvalues of @xmath24 before the stroke and @xmath46 are given by ( [ p1h ] ) and ( [ p2h ] ) . note that we change only the magnetic field , that does not change the internal state and only leads to a zeeman shift . this is an isochoric processs , durng whcih the system release @xmath13 heat to the cold bath ( @xmath90 , see fig . the initial state of the coupled system ( s+phonons ) at this stage can be written as @xmath91 where @xmath92 is the @xmath68-th eigenstate of the hamiltonian @xmath70 with a new value of magnetic field @xmath93 and @xmath94 is the corresponding occupation probability . note that , in the combined basis of internal and motional states , this eigenstate can be written as @xmath95\;.\ ] ] clearly , at thermal equilibrium , the system s is entangled with the phonon state . the occupation probabilities of the states of the system depend upon the temperature @xmath0 . the release of heat from the system to the cold bath is equivalent to cooling down the system . this means that the occupation probabilities of the higher excited states of the system would reduce by such heat release . to facilitate the release of heat from the system to the phonon bath , here we propose a projective measurement of the state of the system s. by measuring the state @xmath14 of the system , the density matrix @xmath96 gets factorized and can be written as @xmath97 in this way , the system releases heat to the cold bath and thereafter gets decoupled from the cold bath .    through this cooling process , the probability distribution of the eigenstates of the system hamiltonian @xmath24 change , keeping the corresponding eigenvalues the same , as local driving fields are kept constant during this stroke . the heat released from the system to the cold bath can be thus be calculated using eqs . ( [ qh ] ) and ( [ geeigen ] ) as @xmath98\;.\ ] ]      during this stroke ( @xmath99 , see fig . 1 ) , the system again undergoes through an adiabatic process , during which the magnetic field strength is now adiabatically changed from @xmath63 to @xmath61 . the system remains in contact with the the hot bath and the vibrational mode - however , as during the expansion stroke , the occupation probabilities of the energy eigenstates @xmath36 remain unaltered at the values @xmath100 . the eigenvalues change from @xmath48 to @xmath43 due to the change in the magnetic field . this leads to the following work done during this stroke : @xmath101      we find that the qoe absorbs heat @xmath12 [ eq . ( [ qh1 ] ) ] during the ignition stroke and releases heat @xmath13 [ eq . ( [ ql1 ] ) ] during the exhaust stroke , while it does certain work during the two adiabatic strokes . therefore , the work efficiency of the qoe can be calculated as @xmath102\{z_{gn}^2-p_n(t_h)\}}{\sum_{n=1}^{2}e_n^h\{p_n(t_h)-z_{gn}^2\}}\;.\nonumber\\\end{aligned}\ ] ]      the compression and expansion strokes in an otto engine are performed adiabatically , during which the system does not share heat with the bath . in the present case , for a quantum version of the otto engine , one needs to consider quantum adiabatic process . in such a process , the probabilities of the different eigenstates of the system hamiltonian remain unaltered , though the relevant eigenvalues are adiabatically changed . the adiabatic condition in quantum mechanics reads as @xmath103 @xcite . in the present case , the adiabatic strokes are performed by changing the magnetic field strength @xmath28 only . using eqs . ( [ en ] ) , we obtain the following condition for adiabaticity : @xmath104 where @xmath105 is the adiabaticity parameter and @xmath106 are given by eq . ( [ npm ] ) . we consider a linear variation of the magnetic field , as given by @xcite @xmath107 in fig . 2 , we show how the condition of adiabaticity is satisfied for different choices of the rate of change @xmath108 of the magnetic field . clearly , the adiabaticity is maintained , even for the evolution for a finite time @xmath109 . with respect to time @xmath110 . the parameters chosen are @xmath111 , @xmath112 , and @xmath113 . clearly @xmath105 remains much less than unity for all the times for larger @xmath109 ( and therefore for slower change of the magnetic field ) , thereby referring to an adiabatic evolution.,width=302,height=226 ]    in the case of @xmath114 for all the times @xmath110 , we can rewrite the adiabaticity parameter @xmath105 as @xmath115 for a change of magnetic field from @xmath61 to @xmath63 in a finite time interval @xmath109 , we obtain the following condition for adiabaticity : @xmath116 it must be reminded that to obtain larger efficiency , one needs to have @xmath117 adiabatically [ see eq . ( [ eff - fin ] ) ] . therefore , as clear from eq . ( [ cond3 ] ) , the corresponding time - interval becomes quite large . we have shown how a single trapped ion can be used to implement quantum otto engine . the electronic degree of freedom of the ion can be considered as the working fluid , that interacts with the thermal bath and the vibrational degree of freedom . in the lamb - dicke limit , the vibrational mode is confined to two lowest eigenstates and therefore can be considered as a cold bath . we show that by adiabatically changing the local magnetic field , the work efficiency can be made close to unity . the interaction with the hot and the cold bath is never switched off . in this framework , we show that the heat release to the cold bath can be mimicked by a projective measurement of the electronic state of the ion . we emphasize that the present proposal can be implemented using the current trapped - ion technology , routinely used for quantum computing . 50 h. e. d. scovil and e. o. schulz - dubois , phys . * 2 * , 262 ( 1959 ) . m. o. scully , m. s. zubairy , g. s. agarwal , and h. walther , science * 299 * , 862 ( 2003 ) . m. o. scully , phys . lett . * 88 * , 050602 ( 2002 ) . m. o. scully , phys . . lett . * 87 * , 220601 ( 2001 ) . t. feldmann and r. kosloff , phys . e * 61 * , 4774 ( 2000 ) ; * 68 * , 016101 ( 2003 ) ; * 70 * , 046110 ( 2004 ) . f. wu , l. g. chen , f. r. sun , c. wu , and q. li , phys . e * 73 * , 016103 ( 2006 ) . h. t. quan , y. x. liu , c. p. sun , and f. nori , phys . e * 76 * , 031105 ( 2007 ) ; h. t. quan , p. zhang , and c. p. sun , phys . e * 73 * , 036122 ( 2006 ) . f. atlintas and o. e. mustecaplioglu , phys . e * 92 * , 022142 ( 2015 ) . g. thomas and r. s. johal , phys . e * 83 * , 031135 ( 2011 ) . t. zhang , w .- t . liu , p .- x . chen , and c .- z . li , phys . a * 75 * , 062102 ( 2007 ) ; he j .- zhou , he xian , and z. jie , chin . b * 21 * , 050303 ( 2012 ) . e. albayrak , int . . inf . * 11 * , 1350021 ( 2013 ) . f. atlintas , a. u. c. hardal , and o. e. mustecaplioglu , phys . e * 90 * , 032102 ( 2014 ) . zhang , eur . j. d * 49 * , 123 ( 2008 ) . j. h. wang , j. z. he , and x. he , phys . e * 84 * , 041127 ( 2011 ) ; j. h. wang and j. z. he , j. appl . phys . * 11 * , 043505 ( 2012 ) . d. gelbwaser - klimovsky , r. alicki , and g. kurizki , phys . e * 87 * , 012140 ( 2013 ) ; d. gelbwaser - klimovsky , r. alicki , and g. kurizki , europhys . lett . * 103 * , 60005 ( 2013 ) ; d. gelbwaser - klimovsky and g. kurizki , phys . rev . e * 90 * , 022102 ( 2014 ) ; sci . rep . * 5 * , 7809 ( 2015 ) ; r. kosloff and a. levy , annu . phys . chem . * 65 * , 365 ( 2014 ) .    l. viola and s. lloyd , phys . a * 58 * , 2733 ( 1998 ) ; d. vitali and p. tombesi , phys . a * 59 * , 4178 ( 1999 ) . o. abah , j. ronagel , g. jacob , s. diffner , f. schmidekaler , k. singer , e. lutz , phys rev . lett . * 109 * , 203006 ( 2012 ) . j. ronagel _ _ , arxiv:1501.03681 . d. leibfried , r. blatt , c. monroe , and d. wineland , rev . 75 * , 281 ( 2003 ) . t. d. kieu , phys . 93 * , 140403 ( 2004 ) .
Abstract:

Answer: we demonstrate how a quantum otto engine ( qoe ) can be implemented using a single ion and an always - on thermal environment . the internal degree of freedom of the ion is chosen as the working fluid , while the motional degree of freedom can be used as the cold bath . we show , that by adiabatically changing the local magnetic field , the work efficiency can be asymptotically made unity . we propose a projective measurement of the internal state of the ion that mimics the release of heat into the cold bath during the engine cycle . in our proposal , the coupling to the hot and the cold baths need not be switched off and on in an alternative fashion during the engine cycle , unlike other existing proposals of qoe . this renders the proposal experimentally feasible using the available tapped - ion engineering technology .


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 identification of the high - redshift precursors of normal present - epoch galaxies , and specifically a galaxy in its first episode of star formation , has long been an elusive goal ( e.g. , partridge & peebles 1967 ) . in past decades , a wide variety of models have been suggested and techniques attempted to identify such proto - galaxies ( see koo 1986 and pritchet 1994 for reviews ) . the proto - galaxy candidate serendipitously discovered by yee et al . ( 1996 ; hereafter paper i ) may offer an excellent chance to study a very young galaxy in detail . the @xmath5 galaxy , designated as ms1512-cb58 , ( hereafter , cb58 ) , lies in the field of the @xmath6 galaxy cluster ms1512 + 36 , which was observed as part of the canadian network for observational cosmology ( cnoc ) cluster redshift survey ( see yee , ellingson & carlberg 1996 ) . the galaxy was found to have @xmath7 by identification of about a dozen strong absorption lines indicative of young stars . optical colors gunn @xmath1 , @xmath8 , gunn @xmath9 and johnson @xmath10 indicated that the stellar population of this galaxy is less than 400 myr old . however , the uncertain extinction correction made it difficult to constrain the galaxy age to better than a range of 10400 myr . p - cygni profiles in the  absorption lines suggested that the rest uv flux is dominated by a population as young as 10 myr old . in this paper we present data in three additional photometric bands in the infrared : @xmath11 , @xmath12 and @xmath2@xmath3 . these data allow stronger limits to be placed on both the extinction and the age of the stellar populations in this galaxy . we also explore the possibility of underlying older populations and discuss briefly the effects of non - standard initial mass functions and gravitational lensing on our conclusions . @xmath11 and @xmath2@xmath3 images were obtained on 1995 sep 28 ut at the nasa infrared - telescope facility ( irtf ) with the nsfcam imager , which contains a 256x256 insb array ( shure & rayner 1993 ) . the plate scale was 0.30@xmath13 per pixel . the data were obtained under photometric conditions with seeing 0.8  0.9@xmath13 . for @xmath2@xmath3 , a 25 point grid with integration times of 15 seconds , 4 co - adds each , was used , resulting in a total integration time of 1500 seconds . for @xmath11 , a 15 point grid was used with 20 second exposures , 3 co - adds each , resulting in 900 seconds of integration . the ukirt faint standards # 28 and 26 ( casali & hawarden 1992 ) were observed for calibration just before and just after the galaxy observations . the data were reduced using standard techniques . images in the @xmath12 band were obtained at the cfht 3.6 m using the redeye - wide infrared camera with a 256@xmath14256 pixel nicmos3 hgcdte array on the night of 1995 dec 30 ut , under photometric skies . the pixel size of the detector is 0.5@xmath13 per pixel . six sets of 9 dithered frames of 35 seconds exposure were obtained , providing a total integration time of 1890 seconds . the photometry was calibrated using observations of the ukirt faint standard star # 23 ( casali & hawarden 1992 ) taken before and after the galaxy observations . photometry in the optical bands of @xmath1 , @xmath8 , @xmath9 , and @xmath10 is taken from paper i. the photometry for the near ir observations was derived in the same way as the optical data ( see paper i ) , using an aperture of 2.2@xmath13 diameter for the object and 9@xmath13 for the outer diameter of the sky aperture . the 2.2@xmath13 aperture photometry is then corrected to the `` total '' magnitude at an aperture equivalent to the @xmath8 isophote of 24.1 mag arcsec@xmath15 , assuming zero color gradients . the photometric data along with their uncertainties are listed in table 1 . also listed are the magnitudes in the ab system , where corrections for the ir photometry were based on the 0 mag fluxes from wamsteker ( 1981 ) . in order to test the calibration of photometry gathered from several different sources , the spectral energy distributions ( seds ) of galaxies of known redshift near cb58 are first examined . figure 1 plots the data for the cd galaxy of the cluster and a bright companion galaxy located 6@xmath13 w and 18@xmath13 s of the cd ( designated # 100983 in the cnoc catalog , abraham et al . the cd galaxy has a spectrum of an old population of stars but with strong [ oii ] emission , and possibly a somewhat blue continuum both presumably associated with the cooling flow detected in this cluster ( donahue , stocke & gioia 1992 ) . figure 1 also shows gissel spectral synthesis models from bruzual & charlot ( 1993 ) . both models are for a 13 gyr old single burst ( coeval ) population with a salpeter initial mass function ( imf ) . this standard imf has an upper mass cutoff 125 , and a lower mass cutoff of 0.1 . the models fit the data reasonably well , indicating that there are probably no gross inconsistencies in the calibration of the seds . note that the model is slightly lower than the cd data bluewards of the 4000   break , as expected .    in paper i , it was pointed out that interpreting the age of a stellar population from rest uv data is complicated by dust extinction . the sed was modeled assuming extinction occurs in the rest frame , and a slightly modified version of the average lmc curve from fitzpatrick ( 1986 ) was used . both here and in paper i , the 2175  hump in the lmc extinction was interpolated over , based on the fact that there is no evidence for the expected prominent dip due to this feature in the optical spectroscopy data . note that many local star burst galaxies also do not show this bump in the extinction in their iue spectra ( kinney et al . an extinction law derived from uv observations of starburst galaxies ( calzetti et al . 1994 ) was also tested and found to yield qualitatively the same results . the extinction in this case has a significantly grayer slope than that from the lmc , and hence requires a larger amount of extinction to provide enough reddening to fit the sed . this in turn will increase the already high rest optical luminosity of this galaxy . hence , the more conservative lmc law is adopted . three single burst gissel models were presented in paper i to illustrate the range of parameters allowed by the optical data : 400 myr old with no extinction , 200 myr old with the equivalent of @xmath16 , and 10 myr old with @xmath17 . the oldest single burst model possible was 400 myr , and that the maximum dust extinction possible was about @xmath18 , since the intrinsic spectrum would become bluer than the hottest stars if more dust is assumed . in figure 2 , these three models are plotted again , this time including the ir data . the spectral range is now from 13006000   rest . each model is normalized to match the optical data , as in paper i , and assumes a salpeter imf with the mass limits noted above . it is clear that the only model which comes close to fitting the new ir data is the youngest model , of a 10 myr old single burst and @xmath17 . the older models clearly predict too much ir flux in the @xmath11 through @xmath19 bands to match the data . no single burst model fits the data within the 95% confidence level and we can conservatively rule out all models older than 20 myr at the 99% confidence level . note that the @xmath11 and @xmath12 bands straddle the 4000    and balmer breaks and can be used alone as an age indicator which is relatively insensitive to extinction . the blue slope of the spectrum in this region clearly indicates that the sed is dominated by a young population . continuous star formation models provide a more physically plausible scenario . for these models , the rest uv flux is always dominated by very young stars , and in paper i it was found that no age discrimination was possible using just the optical data . however , over time it is expected that the observed ir flux will increase . figure 3 shows continuous star formation models of ages 5 , 10 , 50 , and 500 myr , assuming a salpeter imf . here the data have been corrected for extinction equivalent to @xmath17 to match the observed optical data . the ir data clearly favor the 10 myr continuous star formation model , and renormalization of the models relative to the data also provides good fits to ages as high as 20 myr . models older than about 35 myr are ruled out at the 95% confidence level . both of these estimates are consistent with preliminary modeling in paper i of the p - cygni profile of the  absorption line , which suggests a lower limit of 10 myr for a continuous star formation scenario . the ir data show that the observed episode of star formation is relatively young , but does it represent the very first incidence of star formation in this galaxy ? limits on the size of an underlying older population can be estimated from these data but are relatively broad , because even at @xmath19 ( rest @xmath8 ) , the continuum light is still dominated by the young population . the sed is modeled as two episodes of star formation : the observed burst , which is modeled as a 10 myr old continuous episode of star formation ( referred to hereafter as the 10myr model ) , and an older population , which is modeled as a single burst model with varying age and mass relative to the 10myr model . both models assume the standard salpeter imf . for all combinations , the best fits to the data suggest no contribution from an older population . upper limits at the 95% confidence level were determined for the fraction of total stellar mass in the older population . note that as the population ages , its sed differs more in shape from the 10myr model . however , the rapid fading of the single burst sed with age makes it possible to hide a slightly increasing mass fraction underneath the younger distribution . thus , the constraints on the mass fraction of the older population are a weak function of its age , ranging from 83% of the total stellar mass for a 100 myr old population , to 90% for a 2 gyr old population . figure 3 also shows an example of a 1 gyr old single burst combined with the 10myr model , where the older population has a mass equal to 85% of the combined stellar mass . this model can be ruled out at the 95% confidence level . the high redshift of cb58 indicates an upper limit on the age of any older population of about 2 gyr , depending on the cosmological parameters assumed . thus , it is possible to conclude that there is no evidence for an older stellar population , and most conservatively that cb58 has formed at least 10% of its stars within the previous 10 - 20 myr . a lower limit on the total stellar mass can be estimated by noting that the @xmath19 band is approximately equivalent to rest @xmath8 at @xmath7 . the @xmath20 magnitude of 19.61 therefore implies an extinction - corrected rest @xmath8 absolute magnitude of 27.94 mag ( @xmath21 = 75 km s@xmath22 mpc@xmath22 , @xmath23=0.1 , a@xmath24=0.93 mag ) . assuming a 10 myr old single burst model with no underlying older population , the galaxy is expected to fade approximately 6 mag over the next 10 gyr ( the lookback time to @xmath7 ) . note that after 10 gyr the differences between an instantaneous burst and just 10 myr of continuous star formation would be extremely subtle , so the choice of the initial model is not important . thus , the observed star formation in cb58 would be the precursor of a galaxy at least as bright as 21.9 mag , or about 2.7 @xmath25 . any subsequent star formation or the presence of an additional older population would serve to make its present - epoch luminosity even greater , which argues against a significant older population .    the high luminosity of cb58 implies prodigious rates of star formation . paper i estimates star formation rates of 4700  yr@xmath22 based on the 1500  flux and continuous star formation models by leitherer , robert , & heckman ( 1995 ) . the homogeneous appearance of the galaxy suggests that this intense star formation is not localized , but must be occurring across the 30 kpc visible galaxy on close to dynamical timescales . at this rate , a @xmath25 galaxy would be created in a few 10 s of myr . it is not clear whether such a short , intense phase of star formation over the entire galaxy is dynamically possible . it is also unclear whether the implied supernova rates would allow the galaxy to retain any gas for further star formation , or whether supernova - driven winds would quickly strip the galaxy of its interstellar medium . the implied star formation rates and total stellar masses are highly dependent on the initial mass function assumed . variations in the upper mass cutoff do not affect the age of the young population , although they can change the star formation rate by up to a factor of about 5 ( paper i , leitherer et al . lower - mass cutoffs as high as 10  have been predicted for regions of intense star formation ( e.g. scalo 1990 ) . evidence for such a cutoff has been mixed for nearby star forming galaxies , with some galaxies showing near - normal imfs ( e.g. hunter et al . 1995 ) and others showing indications of a lower - mass cutoff ( e.g. rieke et al .  1993 ) . the available gissel models only allow fits to models up to a fairly small 2.5  cutoff . such a low - mass cutoff provides a slightly poorer fit to the observed sed , but is not significantly worse and can not be ruled out . more extreme cutoffs would lower the ir flux of the models , but since massive stars contribute most of the sed , even into the ir , the data are not sufficient to discern between a slightly younger population and a more extreme cutoff . thus , a lower mass cutoff at masses greater than 2.5  might allow a somewhat older age for the galaxy . a similar conclusion would be reached for a flatter imf . a low - mass cutoff would of course also significantly lower the inferred star formation rate and total stellar mass for cb58 . episodes of intense star formation with an extreme low - mass cutoff would temporarily increase the galaxy s luminosity without adding significantly to its stellar mass . in this case , the example of cb58 would not add greatly to our understanding of the origin of present - day stellar populations in normal galaxies . the object is clearly extremely different from many other proto - galaxy candidates in its stellar sed and high luminosity . in comparison with the galaxies at @xmath26 studied by steidel et al . ( 1996 ) , cb58 is about 40 times brighter ( assuming the same extinction ) , and may have slightly higher equivalent width stellar absorption lines . an important difference is that the @xmath27 color of cb58 is 0.8 mag , bluer than the average value of 1.3 mag for the steidel et al .  sample . this may indicate more dust in the fainter galaxies or that they have an older stellar population on average . if the color difference is from dust , an additional @xmath28 of 0.10.2 ( for the lmc and calzetti et al . extinction laws , respectively ) in the steidel et al . galaxies is necessary to explain the color difference . in this case , cb58 would be still a factor of 20 brighter in the rest @xmath8 band . ages for the steidel et al . galaxies of greater than 100 myr would also explain the color difference . if continuous star formation for 10 and 100 myr is assumed for cb58 and the steidel et al . sample , respectively , then cb58 is undergoing star formation at a rate of about 40100 times that of the steidel et al .  sample . an important possibility to consider is that the high luminosity of cb58 is in part due to gravitational lensing from the foreground cluster ms1512 + 36 . the galaxy is located only 6@xmath13 from the cluster cd and is a very strong candidate for lensing . however , the image is resolved on both major and minor axes , appears quite homogeneous and its surface brightness profile can even be fit with an exponential disk model ( paper i ) . williams & lewis ( 1996 ) have modeled the possible lensing of cb58 , assuming a subcritical cluster potential . they found that if the cluster mass of e1512 + 36 is well represented by its velocity dispersion ( 690 ; carlberg et al . 1996 ) , and abell richness ( class 0 ; abraham et al . 1996 ) , then the lensing magnification is at most about 2 magnitudes . however , if the cluster mass is a factor of 2 larger , and assuming a specific geometry , then up to a factor of 40 magnification can be obtained without obvious distortion in the ground - based images . this factor of 2 in mass seems difficult to come by , as it requires a 3@xmath29 error on the velocity dispersion determination ( carlberg et al . williams & lewis note that the x - ray luminosity of the cluster is higher than the velocity dispersion suggests , but the presence of a cooling flow in the cluster will tend to increase the observed x - ray luminosity . furthermore , the x - ray luminosity and the measured velocity dispersion of ms1512 + 36 are well within the scatter of correlation of these quantities found by edge & stewart ( 1991 ) . thus there seems to be only a small possibility that this object is highly lensed , although magnification by factors of a few would not be surprising . unless the gravitational lensing magnification is much larger than the cluster dynamics and image morphology suggest , the star formation rates in this galaxy are probably on the order of at least several hundred to a thousand solar masses per year and are distributed fairly homogeneously across the entire galaxy . clearly such rates can not persist for long , and thus we expect a subsequent decrease in star formation and luminosity of this galaxy in future times . this initial high rate of star formation and later decline suggests that cb58 may be a precursor of a 13@xmath25 early - type galaxy . if the lensing magnification is much larger , then the stellar mass of cb58 would be correspondingly smaller , but its young age and the fraction of stars in the observed 1020 myr old episode of star formation would be unchanged . in this case , cb58 may represent the initial stages of any of a wide variety of galaxy types . we thank the staff of the irtf , especially the telescope operator david griep , for making these observations possible . the cfht observations were made possible by a generous grant of director s discretionary time by pierre couturier . e.e . would like to thank peter conti for numerous useful discussions . wishes to thank cfht for their hospitality while this work was being done . is supported by an operating grant from nserc of canada . abraham , r.g . , yee , h.k.c . , ellingson , e. , gravel , p. , carlberg , r.g . , pritchet , c.j . 1996 , , preprint bruzual , g.a .  & charlot , s. 1993 , , 405 , 538 calzetti , d. , kinney , a.l . , & storchi - bergmann , t. , 1994 , , 429 , 582 casali , m. & hawarden t. 1992 , ukirt newsletter , aug 1992 , 33 donahue , m. stocke , j.t . , & gioia , i.m . 1992 , , 385 , 49 edge , a.c . & stewart , g.c . 1991 , mn , 252 , 418 fitzpatrick , e.l . 1986 , aj , 92 , 1068 hunter , d. a. , shaya , e. j. , holzman , j. a. , light , r. m. , oneil , e.j . , & lynds , r. 1995 , , 448 , 179 kinney , a.l . , bohline , r.c . , calzetti , d. , panagia , n. , & wyse , r.f.g . 1993 , , 86 , 5 leitherer , c. , robert , c. & heckman , t.k . 1995 , , 99 , 173 partridge , b.r . & peebles , p.j.e . 1967 , , 147 , 868 pritchet , c.j . 1994 , pasp , 106 , 1052 rieke , g.h . , loken , k. , rieke , m.j . , & tamblyn , p. 1993 , , 412 , 99 scalo , j.m . , 1990 , in windows on galaxies , ed . g. fabbiano ( dordrecht : kluwer ) , 125 shure , m. & rayner , j. 1993 , irtf newsletter , 8 , no .  2 , 4 steidel , c.c . , giavalisco , m. , pettini , m. , dickinson , m. , & adelberger , k. , 1996 , , 462 , l17 wamsteker , w. 1981 ,  , 97 , 329 williams , l.l.r . , & lewis , g.f . , 1996 , yee , h.k.c . , ellingson , e. , & carlberg , r.g . 1996 , , 102 , 269 yee , h.k.c . , ellingson , e. , bechtold , j. , carlberg , r.g . & cuillandre , j .- c . 1996 , , 111 , 1783    ccccc @xmath1 & @xmath30  4930 & 21.08 & 0.10 & 21.15 @xmath8 & @xmath30  5500 & 20.64 & 0.12 & 20.64 @xmath9 & @xmath30  6540 & 20.60 & 0.10 & 20.41 @xmath10 & @xmath30  8060 & 19.92 & 0.12 & 20.35 @xmath11 & 12400 & 19.12 & 0.13 & 19.95 @xmath12 & 16000 & 18.42 & 0.10 & 19.82 @xmath19 & 21400 & 17.83 & 0.12 & 19.61
Abstract:

Answer: the spectral energy distribution of the proto - galaxy candidate ms1512-cb58 at @xmath0 discovered by yee et al . ( 1996 ) is presented . photometry in seven bands ranging from @xmath1 to @xmath2@xmath3 ( 13006000   rest wavelength ) are fitted with population synthesis models from bruzual & charlot ( 1993 ) . the data confirm a very young age for this galaxy , in agreement with ages estimated from preliminary  p - cygni profile modeling . single - burst models with ages greater than about 20 myr can be discarded at the 99% confidence level , and continuous star formation models with ages greater than about 35 myr can be discarded at the 95% confidence level . the spectral energy distribution is most consistent with a continuous star formation model of about 1020 myr , with reddening of @xmath4 . no evidence for an older population of stars is seen , but the possibility of an older population with as much as 90% of the galaxy mass can not be ruled out . we discuss the possible ramifications of a non - standard imf and gravitational lensing on the galaxy s age and mass . = cmr10