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: cygnus x-2 has long been considered a typical low - mass x - ray binary containing a neutron star ( e.g. , smale 1998 ) with a relatively long orbital period ( @xmath2d ; cowley , crampton & hutchings 1979 ) . recent spectroscopic observations by casares , charles & kuulkers ( 1998 ) and a study of the ellipsoidal light variation by orosz & kuulkers ( 1998 ) yielded accurate values for the mass ratio ( @xmath3 ) , the individual masses ( @xmath4 and @xmath5 for the neutron star and the companion , respectively ) , and the radius of the companion ( @xmath6 ; unpublished ) . these observations confirm the _ present _ low mass of the companion .    rather surprisingly , however , casares et al . ( 1998 ) also showed that the spectral type of the companion is , to within two subtypes , a9iii and that _ the spectral type does not vary with orbital phase_. using the binary parameters above and adopting a temperature of @xmath7k for this spectral type ( straiys & kuriliene 1981 ) , we obtain a luminosity @xmath8 for the companion star . on the other hand , we can also estimate the expected temperature and luminosity for a low - mass subgiant that is consistent with the binary parameters of cyg x-2 : using the calibrated stellar models by han ( 1995 ) , we find an expected temperature of 4500@xmath9k and an expected luminosity of 18@xmath10 . this much higher temperature and luminosity of the companion star could in principle be caused by x - ray heating ( also see cowley et al . however , the fact that the spectral type does not vary with orbital phase despite an inclination @xmath11 ( orosz & kuulkers 1998 ) would imply extremely efficient redistribution of the irradiation energy around the star by irradiation - induced circulation currents , much higher than is , for example , observed in the companion of her x-1 ( kippenhahn & thomas 1979 ) . the spectra of casares & charles ( 1998 ) also show a helium absorption line whose strength varies with orbital phase . the very presence of this line on the companion indicates regions of much higher temperature than found on an a9 star . since the strength of this line is strongest on the illuminated side of the secondary , this suggests that it is caused by external illumination and heating , while the non - varying component of spectral type a9 reflects the true intrinsic spectral type of the companion . since a spectral type a9 is not consistent with a low - mass subgiant , this immediately suggests that the original mass of the companion must have been much higher and that most of the mass of the companion must have been lost from the system . indeed , in  2 we will show that the present parameters of cyg x-2 can be understood if , at the beginning of the mass - transfer phase , the companion was a somewhat evolved star of mass @xmath12 that lost most of its envelope as a result of highly non - conservative mass transfer . cyg x-2 may therefore be related to a hitherto little studied class of intermediate - mass x - ray binaries ( imxbs ) . in  3 we discuss some of the implications of this conclusion and how it may affect our understanding of x - ray binaries in general . we note that a similar suggestion concerning the evolutionary status of cyg x-2 was recently made independently by king & ritter ( 1999 ) . however , as we will show in  2 , our best model ( case ab ) differs substantially from their suggested model ( early case b ) , both in its details and its predictions . in order to explore the possibility that cyg x-2 is the descendant of an intermediate - mass x - ray binary , we performed a series of binary calculations using an up - to - date , standard , henyey - type stellar evolution code ( kippenhahn , weigert & hofmeister 1967 ) . our calculations use solar metallicity ( @xmath13 ) , an initial hydrogen abundance of @xmath14 , a mixing - length parameter @xmath15 , and use a standard prescription for the mass - transfer rate @xmath16 ( see , e.g. , ritter 1988 ) . for each binary - evolution sequence , one needs to specify what fraction , @xmath17 , of the mass lost by the donor is accreted by the neutron star , and the specific angular momentum of any matter that is lost from the system . we somewhat arbitrarily choose @xmath17 to be 1/2 , and limit the maximum accretion rate to the eddington accretion rate , taken to be @xmath18 and kept constant throughout each run . we further assume that the mass lost from the system carries away the specific angular momentum of the accreting neutron star which has an initial mass of 1.4@xmath19 . we include angular - momentum loss due to gravitational radiation ( although it is of no real importance here ) , but do _ not _ consider magnetic braking , since during all slow evolutionary phases , the envelopes of our donor stars tend to be fully radiative ( for a general background review see , e.g. , ritter 1996 ) . we have also performed some calculations that include magnetic braking and found that this would not change any of the main conclusions in this paper . the evolution of a binary is very sensitive to the evolutionary state of the donor star at the beginning of the mass - transfer phase . for example , if the donor is initially relatively unevolved ( early case a ; kippenhahn & weigert 1967 ) , it will subsequently mimic a single star of the same actual mass rather than its initial mass ( e.g. , hellings 1983 ) , provided that is is not too far out of thermal equilibrium . since the companion of cyg x-2 does not presently resemble a single star of the same mass , this type of evolution can not be applicable . to illustrate this evolution , we show in figure  1 the evolutionary track in the hertzsprung - russell ( h - r ) diagram of a 3.5@xmath19donor star that is somewhat evolved ( marked ` case a ' ) . at the beginning of the mass - transfer phase , it has used up only half of its initial hydrogen supply in the center . after a short , rapid mass - transfer phase ( dashed portion ) , the star spends most of the remaining main - sequence phase in the low - mass region of the main sequence ( dot - dashed portion ) and then evolves like a low - mass star up the hayashi line ( solid portion ) . it never resembles a star like the companion of cyg x-2 ( the parameters for cyg x-2 are shown as a boxed region in fig .  1 ) . in order to explain the presently overluminous companion , it must have been relatively evolved at the beginning of the mass - transfer phase . this constraint allows two possibilities : either the secondary was near the end of the main sequence and had already developed a core structure typical of a post - main - sequence intermediate - mass star ( case ab ) or it had already left the main sequence and filled its roche lobe while evolving through the hertzsprung gap ( case b ) . the curves marked ` case ab ' in figure  1 and figure  2 illustrate the evolution of a 3.5@xmath19  star that fills its roche lobe for the first time near the end of the main sequence , when its central hydrogen abundance has been reduced to @xmath20 . the exact value of @xmath21 is not crucial , but should probably be less than @xmath22 . the panels in figure  2 show the radius and roche - lobe radius ( panel a ) , the orbital period ( panel b ) , the masses of the two components and the secondary s core mass ( panel c ) , and the mass - loss rate from the donor ( panel d ) since the beginning of the mass - transfer phase . as is most evident from the last panel , one can clearly distinguish three separate phases , one very rapid phase and two slower phases . _ the rapid initial phase _    because of the large initial mass ratio , the initial mass - transfer rate is of order @xmath23 , and most of the mass lost from the companion has to be ejected from the system . the high @xmath16 implies a mass - loss time scale that is short compared to the thermal time scale of the donor , and the donor star is therefore driven significantly out of thermal equilibrium ( i.e. , it is undersized and underluminous for its mass ) . the mass - transfer rate only starts to decrease significantly after the mass ratio has been reversed and the system , and hence the roche lobe , begin to expand . at this stage , mass transfer is driven entirely by the thermal expansion of the companion . the rapid phase ends once the companion has re - established thermal equilibrium ( after @xmath24yr , of order the thermal time scale of the star ) . by this time , the secondary s mass has decreased to 0.900@xmath19 . since the material that is now exposed at the surface has undergone partial cno burning , its composition shows the signature of cno processing ( enhanced nitrogen , decreased carbon and oxygen ) and the surface hydrogen mass fraction is reduced to 0.55 ( from 0.7 ) .    _ the hydrogen core - burning phase _    once the star has returned to thermal equilibrium , the further evolution is driven by the nuclear evolution of the core ( hydrogen core burning ) . during this phase , which in this example lasts @xmath25yr , the mass - transfer rate is of order @xmath26 and the mass decreases to 0.876@xmath19 . the phase ends when hydrogen has been exhausted in the core and the system becomes temporarily detached , accompanied by a small hook in the h - r diagram ( fig .  1 ) . _ the hydrogen shell - burning phase _    after the exhaustion of central hydrogen , the star expands again and soon starts to fill its roche lobe for a second time . initially the evolution of the secondary resembles that of a low - mass star evolving off the main sequence . however , its behavior changes qualitatively once it has lost all of the material that was outside the convective core at the beginning of the mass - transfer phase ( which had a mass of @xmath1 ) . at this point , the surface hydrogen abundance drops to 0.1 ( the central hydrogen abundance in the core at the beginning of mass transfer ) , and the star now has many of the characteristics of a non - degenerate helium star , except for the fact that the evolution , and hence mass transfer , are driven by hydrogen burning in an outward - moving shell . unlike the case of early case b ( see   2.2 ) , the star has now no tendency to become a giant , which would lead to mass transfer on a thermal time scale ( we have tested this by evolving such a star without further mass loss ) . the mass - transfer rate is therefore determined by the time scale for hydrogen - shell burning and lies in the range of @xmath27 to @xmath28 ( see the inset in fig . this phase ends when the hydrogen - rich envelope has almost been exhausted and hydrogen - burning stops . the star still continues to evolve , now as a detached star , burns helium in the core as an o subdwarf ( with a luminosity @xmath29 and temperature @xmath30k ) and eventually ends its evolution as a co white dwarf with an unusually low mass of @xmath31 . the final neutron - star mass in this calculation is @xmath32 .    _ the appearance of the secondary _    in the initial rapid phase ( dashed portion of curve ab in fig .  1 ) , most of the transferred mass has to be ejected from the system . since the mass - loss rate exceeds the eddington rate by several orders of magnitude , the system is unlikely to look like a typical x - ray binary . since this phase is rather short - lived ( @xmath33 ) , relatively few systems should be found in this phase . one possibly related , somewhat more massive observed counterpart is the famous system ss433 , which appears to eject most of the transferred mass by means of two relativistic jets ( e.g. , margon 1984 ) . on the other hand , in the two nuclear - evolution - driven phases , the system would in most respects resemble a typical low - mass x - ray binary ( lmxb ) and would almost certainly be observationally classified as a low - luminosity lmxb in the core - burning phase , and a high - luminosity lmxb in the shell - burning phase . both phases are also relatively long - lived ( @xmath34 and @xmath35 ) and hence relatively more systems should be found in these phases than in the rapid phase ( see   3 for discussion ) .    in this particular model , the star spends a large fraction of the shell - burning phase near a location in the h - r diagram where cyg x-2 lies ( the dashed portion of track ab in fig .  1 ) . indeed during the phase when the orbital period increases from @xmath36 to @xmath37d , the model provides an excellent description of the main observed properties of cyg x-2 ( see   1 ; the mass of the secondary is in the range of 0.45 and 0.5@xmath19 ) . an important prediction of this model is that the surface composition should be the same as the _ core composition _ of the secondary at the beginning of the first mass - transfer phase , i.e. , be severely hydrogen depleted ( in this particular model , @xmath38 ) and show strong evidence of cno processing . ( we caution , however , that there are a sufficient number of uncertainties in the binary - evolution model , in particular associated with the angular - momentum loss from the system , that somewhat different initial binary parameters can almost certainly also reproduce the observed properties of cyg x-2 . )      an alternative model for cyg x-2 is that it started to fill its roche lobe in the hertzsprung gap ( early case b ; also see kolb 1998 and , in particular , king & ritter 1999 ) . we have calculated this type of evolution for a 3.5@xmath19  star that has just left the main sequence ( dashed curve marked ` case b ' in fig .  1 ) . the initial evolution is similar to the previous case . the mass - transfer rate ( dashed curve in fig . 2(d ) ) reaches a peak of @xmath39 and starts to drop once the mass ratio has been reversed . in the subsequent somewhat , but not much slower phase , mass - transfer is driven by the thermal evolution of the envelope ( the star is trying to evolve across the hertzsprung gap to become a giant ) . this phase ends and mass transfer stops completely once most of the hydrogen - envelope has been lost and the hydrogen - burning shell has been extinguished . the secondary continues to evolve as a detached star , burning helium as an o subdwarf and finally becomes a co white dwarf of mass @xmath40 ( the neutron star has only a slightly increased final mass of @xmath41 ) . as is clear from figure  2(d ) , apart from the initial phase , the mass - transfer history is very different from the case ab scenario above . in the slower phase , mass transfer is driven by the thermal expansion of the secondary rather than the nuclear evolution of the core or shell . since the thermal time scale is much shorter than the evolutionary time scale during the hydrogen - shell burning phase in the case ab scenario , the mass - transfer rate in the case b scenario is @xmath42 , an order of magnitude above the eddington rate . it is not clear whether such a high mass - transfer rate would be consistent with the observed properties of cyg x-2 . this problem is somewhat reduced if the initial secondary mass was somewhat lower ( @xmath43 ) . however , all case b scenarios suffer from the generic problems that the mass - transfer rate is significantly super - eddington and that the characteristic lifetime of this phase is very short ( of order the thermal time scale of the star , @xmath33 ) . in all the case b models we have calculated , the secondaries spend very little time in the h - r diagram in the general neighborhood of cyg x-2 . for even lower masses , the thermal time scale would be longer , but the star would have a tendency to evolve towards the hayashi line , which again would be inconsistent with the inferred parameters of the secondary in cyg x-2 . we have performed a whole series of early - case - b binary sequences . none of them was able to reproduce all of the observationally inferred parameters of cyg x-2 simultaneously ( i.e. , the position in the h - r diagram , the constituent masses , and @xmath16 ) . while we can not strongly rule out an early - case - b scenario for cyg x-2 , keeping in mind the general uncertainties in modeling non - conservative mass transfer , we conclude that it is more likely that cyg x-2 is at present in a hydrogen - shell burning phase in a case ab scenario . an important implication of our modeling of cyg x-2 is that it suggests that there may be a whole class of x - ray binaries which either have intermediate - mass secondaries now or had them in the past . such systems have received hardly any attention in the past ( see , however , pylyser & savonije 1988 , 1989 ; kolb 1998 ) . to some degree , this is the result of observational findings , since previously only one such system ( her x-1/hz her ) had been unambiguously identified . in addition , the theoretically predicted large initial mass - transfer rates are much higher than those typically inferred for observed x - ray binaries . however , as our calculation show , this rapid phase is relatively short - lived and very few systems should be found in this phase . it is also not clear what they would look like . a further theoretical uncertainty is how the mass is lost in this rapid phase . as noted earlier , ss433 may provide an observational answer to both of these questions ( the appearance initially and the mass - loss mechanism ) . on the other hand , from an evolutionary point of view , imxbs should be relatively common , since they are much easier to form and do not require the same of amount of fine - tuning as lmxbs ( see , for example , bhattacharya & van den heuvel 1991 ) . the observational properties of imxbs would vary significantly , depending on the initial mass of the secondary and the evolutionary phase at the beginning of the mass - transfer phase . we have initiated a systematic investigation of imxbs and found that , in the longer - lasting evolutionary phases , they typically resemble ` classical ' lmxbs ( podsiadlowski & rappaport 1999 ) . this suggests that perhaps a significant fraction of x - ray binaries , presently classified as lmxbs , are actually imxbs or their descendants . this may have important implications for our understanding of x - ray binaries . for example , it may help to explain some of the systematic differences between galactic and globular - cluster lmxbs , the observed period and luminosity distributions , and may help to resolve the apparent discrepancy between the number of millisecond pulsars and their putative progenitors ( e.g. , kulkarni & narayan 1988 ) . there is a group of about two dozen radio pulsars that are found in nearly circular orbits with low - mass white - dwarf companions ( see , e.g. , taylor , manchester & lyne 1993 ; rappaport et al . the orbital periods of these systems range from less than 1 day to 1232 days . it has been conjectured that this type of system evolves from a low - mass ( @xmath44 ) subgiant or giant donor star which transfers its envelope to a neutron star via stable roche - lobe overflow ( see , e.g. , joss , rappaport & lewis 1987 ; rappaport et al . what remains is a low - mass white dwarf ( the core of the giant ) in a wide , nearly circular orbit about the neutron star which has been spun up by the accretion process . joss , rappaport , & lewis ( 1987 ) and rappaport et al .  ( 1995 ) showed that there is a nearly unique relation between the final orbital period , @xmath45 , of such a `` recycled '' binary pulsar and the mass of the white dwarf , @xmath46 ( see also savonije 1987 , and refsdal & weigert 1970 , 1971 for related discussions ) . in short , this relation results from the fact that the radius of a low - mass giant is a nearly unique function of its core mass . it then follows for a roche - lobe filling star in a binary that @xmath45 , at the end of mass transfer , is a function of only the core mass of the donor . thus , the binary remains as a fossil relic of the giant donor star and its degenerate he core . the theoretical relation between @xmath45 and @xmath46 is found to be roughly consistent with the inferred masses of the white dwarf companions in these systems . unfortunately , this model can not be confirmed to the level of confidence that one would like . first , for most of these systems , only the mass function is measured , and thus the uncertainty in the white dwarf mass is rather large . second , there are a number of these systems for which the orbital period is really too short to be explained by this scenario ( i.e. , @xmath47d ) . finally , there are a few systems where the lower limit on the mass of the white dwarf is only marginally consistent with the theoretical relationship . the last of these difficulties could be mitigated if some of these systems formed with donors of intermediate mass . we have shown in this work that neutron stars with intermediate - mass donor stars can attain a final evolutionary state which is very much like that of the binary radio pulsars discussed here . in particular , the system which we use to model cyg x-2 would end its evolution with a 0.42@xmath19  co white dwarf in a nearly circular 12-day orbit about a spun - up neutron star . since the initial mass transfer in this system occurred when the 3.5@xmath19  donor star was late in its main - sequence phase ( late case a ) , the core was neither completely hydrogen depleted nor degenerate . this combination enables the final orbital period to remain relatively short for a white - dwarf mass as large as 0.42@xmath19 . by contrast , in the more standard low - mass giant scenario the corresponding final orbital period with this same white dwarf mass would have been @xmath48 days . put another way , in this same scenario , a system with a 12-day orbital period would be expected to have @xmath49 . in the future , deep spectroscopy of the companion white dwarfs ( see , e.g. , van kerkwijk & kulkarni 1997 ) , as well as the possible detection of the shapiro delay in the arrival times of the radio pulses ( see , e.g. , kaspi , taylor , & ryba 1994 ) , could substantially improve the determination of the white dwarf masses and thereby discriminate between low - mass and intermediate - mass progenitors . it is also possible that neutron stars accreting from intermediate mass donor stars could attain a final evolutionary state consisting of a very low - mass white dwarf in a short orbital period binary ( e.g. , @xmath50d ) , and help explain binary radio pulsars with such properties . we plan to explore such scenarios in future work . thus , evolutionary scenarios and population synthesis studies which are employed to reproduce the distribution of orbital periods and white dwarf masses in binary radio pulsars , must include consideration of intermediate mass donor stars . we have shown that the observationally inferred parameters of the x - ray binary cyg x-2 can be understood if the secondary had an initial mass of @xmath12 and started to transfer mass near the end of the main sequence ( or less likely just after leaving the main sequence ) . cyg x-2 is therefore related to a class of intermediate - mass x - ray binaries that has been little studied before . our favored model for cyg x-2 ( case ab ) predicts that the secondary should , at present , be severely hydrogen depleted . we suggest that perhaps a significant fraction of x - ray binaries , that are presently classified as lmxbs , could be related to this class . this could have far - reaching implications for our understanding of x - ray binaries . detailed future observations of the secondaries in x - ray binaries as well as systematic binary population synthesis studies should help to assess the importance of this relatively unexplored evolutionary channel . bhattacharya , d. , van den heuvel , e. p. j. 1991 , phys.rep . , 203 , 1 casares , j. , charles , p. , & kuulkers , e. 1998 , , 493 , l39 casares , j. , & charles , p. 1998 , private communication cowley , a. p. , crampton , d. , & hutchings , j. b. 1979 , , 231 , 539 han , z. 1995 , ph.d . thesis , cambridge hellings , p. 1983 , , 96 , 37 joss , p. c. , rappaport , s. , & lewis , w. 1987 , apj , 319 , 180 . kaspi , v. m. , taylor , j.h . , & ryba , m.f . 1994 , apj , 428 , 713 . king , a. r. , & ritter , h. 1999 , , submitted kippenhahn , r. , & thomas , h .- c . 1979 , a&a , 75 , 281 kippenhahn , r. , & weigert , a. 1967 , z. astrophys . , 65 , 251 kippenhahn , r. , weigert , a. , & hofmeister , e. 1967 , in methods in computational physics , vol .  7 , ed .  b. alder , s. fernbach , & m. rothenberg ( new york : academic press ) , 129 kolb , u. 1998 , , 297 , 419 kulkarni , s. r. , & narayan , r. 1988 , apj , 355 , 755 margon , b. 1984 , ara&a , 22 , 507 orosz , j. a. , & kuulkers , e. 1998 , , submitted podsiadlowski , ph . , & rappaport , s. 1999 , in preparation pylyser , e. h. p. , savonije , g. j. 1988 , a&a , 191 , 57 pylyser , e. h. p. , savonije , g. j. 1989 , a&a , 208 , 52 rappaport , s. , podsiadlowski , ph . , joss , p. c. , distefano , r. , han , z. 1995 , , 273 , 731 refsdal , s. , & weigert , a. 1970 , a&a , 6 , 426 refsdal , s. , & weigert , a. 1971 , a&a , 13 , 367 . ritter , h. 1988 , a&a , 202 , 93 ritter , h. 1996 , in evolutionary processes in binary stars , ed . r. a. m. j. wijers , m. b. davies , & c. a. tout ( dordrecht : kluwer ) , 223 savonije , g. j. 1987 , nature , 325 , 416 smale , a. p. 1998 , , 498 , l141 straiys , v. , & kuriliene , g. 1981 , , 80 , 353 taylor , j. h. , manchester , r. n. , & lyne , a. g. 1993 , apjs , 88 , 529 van kerkwijk , m. , & kulkarni , s. 1997 , in iau symp . 188 , the hot universe , kyoto , japan , in press      * figure 1 . * evolutionary tracks of the secondaries in three binary calculations in the hertzsprung - russell diagram . the secondary has an initial mass of @xmath0 and the primary , assumed to be a neutron star , an initial mass of @xmath51 in all calculations . the dotted curve shows the track of a @xmath0 star without mass loss . the mass - loss tracks , labelled case a , ab and b , start at different evolutionary phases of the secondary ( ` case a ' : the middle of the main sequence ; ` case ab ' : the end of the main sequence ; ` case b ' : just after the main - sequence ) . the dashed portions in each track indicate the rapid initial mass - transfer phase , the dot - dashed and solid portions the slow phases where mass - transfer is driven by hydrogen core burning and hydrogen shell burning , respectively ( only in case a and case ab ) . the beginning and end points of the various phases are marked by solid bullets , the small figures next to them give the mass of the secondary at these points . the boxed region labelled ` cyg ' indicates the observationally determined parameter region for the secondary in cyg x-2 . the tracks after mass transfer has ceased are not shown . * figure 2 . * key binary parameters for the case ab binary calculation as a function of time ( with arbitrary offset ) . panel ( a ) : radius ( solid curve ) and roche - lobe radius ( dot - dashed curve ) of the secondary ; panel ( b ) : the orbital period ( solid curve ) ; panel ( c ) : the mass of the secondary ( solid curve ) , of its hydrogen - exhausted core ( dotted curve ) , and of the primary ( dot - dashed curve ) ; panel ( d ) : the mass - loss rate from the secondary ( solid curve ) ; the inset shows a blow - up of the second slow mass - transfer phase ( hydrogen shell burning ) . the dashed curves in panels ( b ) and ( d ) show the orbital period and mass - transfer rate for the case b calculation for comparison .
Abstract:

Answer: the x - ray binary cygnus x-2 ( cyg x-2 ) has recently been shown to contain a secondary that is much more luminous and hotter than is appropriate for a low - mass subgiant . we present detailed binary - evolution calculations which demonstrate that the present evolutionary state of cyg x-2 can be understood if the secondary had an initial mass of around @xmath0 and started to transfer mass near the end of its main - sequence phase ( or , somewhat less likely , just after leaving the main sequence ) . most of the mass of the secondary must have been ejected from the system during an earlier rapid mass - transfer phase . in the present phase , the secondary has a mass of around @xmath1 with a non - degenerate helium core . it is burning hydrogen in a shell , and mass transfer is driven by the advancement of the burning shell . cyg x-2 therefore is related to a previously little studied class of intermediate - mass x - ray binaries ( imxbs ) . we suggest that perhaps a significant fraction of x - ray binaries presently classified as low - mass x - ray binaries may be descendants of imxbs and discuss some of the implications .


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 properties of crystalline substances often crucially depend on the structure and dynamics of imperfections of the crystal lattice . in particular , point defects such as interstitials and vacancies play a pivotal role in determining the stability , transport properties , growth characteristics , and mechanical behavior of materials . recent impressive experimental advances , such as optical tweezers and confocal microscopy @xcite , now permit to study the fundamental properties of point defects in condensed matter systems with `` atomistic '' space and time resolution . recently , a number of experimental studies have focused on the structure and dynamics of point defects in two - dimensional assemblies of micrometer sized colloidal particles @xcite and , in particular , on their effective interactions @xcite . in studying such defect interactions the question arises to which degree they can be rationalized in terms of continuum elastic theory . as a first step towards answering this question , in this article we investigate numerically the disturbances caused by isolated interstitial particles and compare the results with the predictions of continuum theory . in carrying out such a comparison , it proves crucial that in solving the equations of elasticity theory boundary conditions are used that match those of the simulations . for the periodic boundary conditions usually applied in computer simulations , the displacement fields of single defects can be determined using the technique of ewald summation familiar from electrostatics @xcite . while elasticity theory properly describes the disturbances and interactions created by lattice imperfections on a larger scale , discrete lattice effects dominate on spatial scales of the order of few lattice constants . the remainder of this paper is organizes as follows . in sec . [ sec : model ] we define the model and describe the numerical methods . the treatment of point defects in a two - dimensional elastic continuum is discussed in sec . [ sec : elasticity ] and comparison with the numerical results is discussed in sec . [ sec : results ] . for certain defect configurations one observes an exponential rather than algebraic decay of the displacement fields . this behavior can be understood in terms of a simple bead - spring model introduced in sec . [ sec : spring ] with parameters related to the elastic constants of the material . some concluding remarks are provided in sec . [ sec : discussion ] . in this paper we study a two - dimensional crystal of soft particle interacting via the gaussian potential @xcite @xmath0 where @xmath1 is the inter - particle distance and @xmath2 and @xmath3 set the energy and length scales , respectively . in the following , energies are measured in units of @xmath2 and distances in units of @xmath3 . this so - called gaussian core model , used here as a generic model for a system of soft spheres , is a realistic description for the short - ranged effective interactions between polymer coils in solution @xcite . in three dimensions , the gaussian core model can exist as a fluid , a bcc- and an fcc - solid depending on temperature and density @xcite . in two dimensions , the perfect triangular lattice is the lowest energy structure of gaussian core particles at all densities @xcite . computer simulations indicate that also in this system of purely repulsive particles point defects such as interstitials , vacancies or impurity particles of different size display attractive ( as well as repulsive ) interactions both in two and three dimensions @xcite .    to study the displacement field of a single interstitial numerically , we prepare a configuration of particles arranged on the sites of a perfect lattice configuration and insert an extra particle of the same species . after insertion , the system is relaxed to a new minimum energy configuration by steepest descent minimization , i.e. , we study the defect structure at @xmath4 . typically , @xmath5 steepest descent steps are carried out . the system we study here consists of @xmath6 gaussian core particles ( without the extra particle ) at a number density of @xmath7 corresponding to a lattice constant @xmath8 . periodic boundary conditions apply to the simulation box of length @xmath9 and height @xmath10 . the aspect ratio of the almost square simulation box is @xmath11 . we quantify the perturbation caused by the defect in terms of the displacement field @xcite @xmath12 here , @xmath13 and @xmath14 denote the position of particle @xmath15 with and without the defect , respectively . as we will see in the following sections , simple point defects generate remarkably intricate displacement patterns that can be understood in terms of elasticity theory only on large length scales .    at @xmath4 , the elastic constants describing the macroscopic response of the system to perturbations can be calculated as a function of density from simple lattice sums . for a density of @xmath16 , the lam coefficients ( see sec . [ sec : elasticity ] ) of the perfect triangular lattice have values @xmath17 and @xmath18 . at this density , the pressure is @xmath19 and the energy density is @xmath20 corresponding to an energy per particle of @xmath21 . the bulk modulus , which in two dimensions is related to the the lam coefficients by @xmath22 , has a value of @xmath23 . while close to a point defect the displacement field is highly anisotropic and strongly dependent on the atomistic details of the interactions , for large distances elasticity theory is expected to be valid . the differential equations describing the equilibrium of an elastic continuum are usually expressed in terms of the strain tensor @xcite @xmath24 where @xmath25 denotes the @xmath15-component of the displacement @xmath26 and @xmath27 the @xmath15-th component of the position @xmath28 . for a given external volume force @xmath29 with components @xmath30 acting on an isotropic system such as a crystal on a triangular lattice , hook s law leads to the equilibrium condition for the strain : @xmath31 here , @xmath32 and @xmath33 are the so - called lam coefficients and summation over repeated indices is implied . solving this equation for a singular force yields the green s function from which the response of the elastic continuum to an arbitrary force can be obtained by integration . to model the displacement field caused by the introduction of point defects using linear continuum elasticity theory , we determine the displacement field caused by two pairs of opposing forces , one pair acting along the @xmath34-axis and the other one along the @xmath35-axis @xcite . this idealized model of a defect is equivalent to inserting a small circular inclusion into a hole of different size @xcite . each force of this pair is of equal magnitude @xmath36 but with opposite sign acting on two points separated by a small distance @xmath37 . such a force pair exerts a zero net force on the material . in the limit @xmath38 where the force @xmath39 in a way such that @xmath40 remains constant , the equilibrium condition for the displacement can be written as @xmath41 assuming that the displacement can be written as the derivative of a potential , @xmath42 one obtains @xmath43 this equation is the poisson equation of electrostatics with a singular disturbance . since , as noted above , @xmath44 is a solution of @xmath45 ( see , for instance , ref . @xcite ) , we obtain the green s function @xmath46 from which the displacement field @xmath47 follows by differentiation according to equ . ( [ equ : displacement_phi ] ) , @xmath48 in comparing the results of particle simulations with those of elasticity theory it is important to realize that the displacement fields predicted by continuum theory are of a long - range nature . therefore , it is crucial that corresponding boundary conditions are used in both cases . all simulations discussed in this paper are done with periodic boundary conditions in order to minimize surface effects and preserve the translational invariance of the perfect lattice . hence , also the continuum calculations need to be carried out with periodic boundary conditions . since the defect fields for the infinite material are long - ranged , the displacement field in the periodic system can not be obtained by simply summing up the contributions of the periodic images . in fact , such a naive summation of the contribution of all image defects diverges . a more appropriate treatment that avoids this problem consists in determining the green s function of the poisson equation ( [ equ : poisson ] ) for periodic boundary conditions . in this case , the solution of this equation in two dimensions , known from electrostatics @xcite , can be written as ewald sum of a logarithmic potential embedded in a neutralizing background , @xmath49 \right . \nonumber\\       & &     \hspace{1.5 cm } \left . - \frac{2 \pi}{a } \sum_{{\bf k } \neq 0 } \frac{e^{-k^2/4\eta^2}}{k^2 } \cos({\bf k}\cdot { \bf r } )   + \frac{\pi}{2 \eta^2 a } \right\}.\end{aligned}\ ] ] here , @xmath50 is the exponential integral . the first sum is over all lattice vectors @xmath51 in real space and the second sum is over all reciprocal vectors @xmath52 in fourier space . the adjustable parameter @xmath53 , set to a value of @xmath54 here , determines the rate of convergence of the two sums and @xmath55 is the area of the rectangular simulation cell . the fourier space sum can be evaluated accurately using about 2,500 reciprocal space vectors . from equ . ( [ equ : ewaldpotential ] ) for the scalar function @xmath56 the displacement field of a point defect in a system with periodic boundary conditions is found by differentiation , @xmath57 for the systems considered in this paper , the real space sum may be truncated after the first term . since the value of @xmath58 is undetermined , the parameter @xmath59 is treated as a fit parameter in the following . the ewald sums of the above equations describe the effects of `` image defects '' at the center of the periodically replicated domains . first , we study the displacement field of a single interstitial . to generate such a defect , we insert an extra particle of the same species into a perfect 2d - crystal on a triangular lattice . after insertion , the system is relaxed to a new minimum energy configuration by steepest descent minimization , i.e. , we study the defect structure at @xmath4 . typically , @xmath5 steepest descent steps are carried out . in each step each particle is moved in the direction of the force acting on the particle where the absolute value of the displacement in chosen to be small enough to ensure that the energy of the system decreases in each step . ( a )     ( a ) , the @xmath60 ( b ) and the @xmath61 ( c ) configurations . the length of the arrows representing the displacements of the particles from their position in the perfect lattice are exaggerated for better visibility . in the figures at the bottom the blue spheres represent the particles and the small gray spheres indicate the position of the lattice sites of the perfect crystal.,width=142 ]    ( b )     ( a ) , the @xmath60 ( b ) and the @xmath61 ( c ) configurations . the length of the arrows representing the displacements of the particles from their position in the perfect lattice are exaggerated for better visibility . in the figures at the bottom the blue spheres represent the particles and the small gray spheres indicate the position of the lattice sites of the perfect crystal.,width=142 ]    ( c )     ( a ) , the @xmath60 ( b ) and the @xmath61 ( c ) configurations . the length of the arrows representing the displacements of the particles from their position in the perfect lattice are exaggerated for better visibility . in the figures at the bottom the blue spheres represent the particles and the small gray spheres indicate the position of the lattice sites of the perfect crystal.,width=142 ]    the extra particle can deform the crystal in different ways @xcite and produces displacement fields of different symmetries ( see fig . [ fig : displacementinter ] ) . in one configuration , called @xmath62 interstitial or crowdion and shown in fig . [ fig : displacementinter]a , the additional particle pushes one particular particle of the crystal out of its equilibrium position . both the original particle and the additional particle arrange themselves at equal distance around the lattice position of the original particle . the displacement pattern arising for this type of interstitial has two - fold symmetry and , of course , occurs in all three low - index lattice directions with equal probability . one may suspect that this defect configuration , with a symmetry that differs from the symmetry of the underlying triangular lattice , is caused by the rectangular periodic boundary conditions that are applied to the system . to rule out this possibility , we have repeated the calculation with hexagonal periodic boundary conditions obtaining the same result . another low - energy defect configuration is the @xmath60 interstitial with three - fold symmetry ( see fig . [ fig : displacementinter]b ) . in this case the interstitial particle is located at the center of a basic lattice triangle and pushes its neighbors outward from their original positions . a third important interstitial configuration is the @xmath61 interstitial or dumbbell interstitial shown in fig . [ fig : displacementinter]c . in the 2d gaussian - core model under the conditions studied here the @xmath62 pattern has a slightly lower energy than the @xmath60 interstitial and the @xmath61 interstitial . the energy difference between a @xmath62 and a @xmath60 interstitial is @xmath63 and the difference between @xmath62 and @xmath61 is @xmath64 . all three displacement patterns are important for the diffusion of interstitials . an @xmath62 interstitial is very mobile in the direction of its main axis . the @xmath60 and @xmath61 forms are visited as intermediate configurations when the @xmath62 interstitial changes the orientation of its main axis and hence its direction of motion @xcite . next , we compare the displacement fields determined numerically with the predictions of elasticity theory . in particular , we verify to which extent the @xmath65-behavior modulated by the periodic boundary conditions and embodied in equ . ( [ equ : ewalddisplacement ] ) is realized in the particle system . the complex displacement patterns of the various interstitial configurations shown in fig . [ fig : displacementinter ] obviously differ from this expectation , at least near the defect , and indicate that continuum theory is not applicable in this region . far away from the defect , however , the perturbation caused by the defect is small and the response of the material should be described accurately by linear elasticity theory .     as a function of distance from the defect @xmath1 for the @xmath62 interstitial . each red dot corresponds to one particle . the solid line represents the @xmath66 behavior . here , @xmath67 was used as this value yields the best fit of the results obtained vie ewald summation to the results of the particle simulations in the far field . inset : angle @xmath68 between the displacement vector @xmath26 the position vectors @xmath28 as a function of the distance @xmath1 from the defect site . , width=302 ]    the magnitude @xmath69 of the displacement vector @xmath47 is shown as a function of the distance from the interstitial in fig . [ fig : disp_theta ] for the @xmath62 defect configuration . each point in the figure corresponds to one individual particle . for short distances , the displacement magnitude is not a unique function of the distance @xmath1 reflecting the anisotropic nature of the defect . for larger distances , however , the displacement magnitude is mostly determined by the distance @xmath1 . eventually , however , the periodic boundary conditions lead to a spread of the displacement magnitude for even larger distances and a splitting into two branches corresponding to the @xmath34- and @xmath35-directions and the directions along the diagonals , respectively . in the regime where @xmath47 behaves isotropically , the displacement follows the approximately @xmath65-form predicted by elasticity theory for a point defect in an infinitely extended medium . the orientation of the displacement vector @xmath47 , depicted in the inset of fig . [ fig : disp_theta ] , behaves in an analogous way . the angle @xmath68 between @xmath47 and the position vector @xmath28 , shown as a function of the distance @xmath1 from the defect , is not a unique function of @xmath1 near the defect . for larger @xmath1 , @xmath68 vanishes indicating that in this distance regime the displacement vector points straight away from the defect . at even larger distances , the periodic boundary conditions imposed on the system eventually cause the angle @xmath68 to spread again .     and @xmath70 of the @xmath71 interstitial as s function the of the distance along the @xmath34-axis and @xmath35-axis , respectively ( solid lines ) . here , the direction of largest displacement of the @xmath71 defect is oriented in @xmath34-direction . also plotted is the displacement computed from continuum theory according to equ . ( [ equ : ewalddisplacement ] ) ( dashed line ) , simple @xmath65-behavior ( dotted line ) and the displacement obtained for the simple mechanical model described in the main text ( dash - dotted line ) . the inset shows the region close to the defect location . as in fig . [ fig : disp_theta ] a defect strength of @xmath72 was used for the evaluation of the displacement from elasticity theory . , width=302 ]    the displacement fields calculated according to equ . ( [ equ : ewalddisplacement ] ) and numerically for an interstitial in the @xmath62 configuration are compared in fig.[fig : absdisplacement ] . in this figure , the displacement components @xmath73 and @xmath70 are depicted as a function of the distance from the defect along the @xmath34-axis and @xmath35-axis , respectively . the prediction of continuum theory , calculated using the ewald summation of equ . ( [ equ : ewalddisplacement ] ) , agrees well with the displacement field of the particle system for distances larger than about 10 lattice constants . near the defect , the predictions of continuum theory differ from the simulation results . the deviation is particularly pronounced in the direction of the main axis of distortion of the @xmath62 defect , in which the displacement appears to decay exponentially up to a distance of about @xmath74 . this unexpected exponential behavior can be understood in terms of a simple model with harmonic interactions . this model consists of a one - dimensional chain of particles in which each particle is connected to its two neighbors with springs of force constant @xmath75 ( except the first and last particle , which are coupled only to their neighbors on the right and left , respectively ) . in addition , each particle is attached to a fixed lattice position with another spring of force constant @xmath76 . the hamiltonian of this system is @xmath77 where @xmath78 is the number of particles , @xmath79 is the position of particle @xmath80 and @xmath81 is the equilibrium distance between two neighboring particles . in the minimum energy configuration of this chain , the particles are arranged such that @xmath82 . we now imagine that particle @xmath83 is pushed to the right by a distance of @xmath84 while particle @xmath85 is kept fixed at @xmath86 . if the system is then relaxed to a new energy minimum , all other particles will be displaced from their original positions too . for this simple model , the response of the system to the displacement of the first particle can be calculated analytically by direct matrix inversion ( see [ ap : springmodel ] ) . in the large @xmath87 limit , one finds that the displacement of the particles from their original position decays exponentially with their position , @xmath88 where @xmath89 is the displacement of particle @xmath80 due to the forced displacement @xmath84 of the first particle . the decay constant @xmath90 is related to the force constants of the model by @xmath91 to compare the prediction of this simple model with the simulation results we have to determine the force constants @xmath75 and @xmath76 felt by the particles in the main axis of the defect . while the force constant @xmath75 arises from interactions within this main axis , the force constant @xmath76 is related to interactions of the particles in the main axis with those from adjacent rows . accordingly , we determine @xmath75 by calculating numerically the energy change caused by slightly displacing one single particle in a one - dimensional row of otherwise fixed gaussian core particles without the presence of the neighboring rows . the distance of the particles in the row is chosen to be equal to the lattice constant at the density @xmath92 considered throughout the paper . from the energy as a function of the displacement one obtains a force constant of @xmath93 . to determine the force constant @xmath76 we calculate the energy change caused by translating a whole row of particles in the perfect crystal . the particles in the row are fixed with respect to each other and the remaining particles are kept at their lattice positions . from the energy change per moved particle a force constant of @xmath94 follows . the decay constant of @xmath95 calculated according to equ . ( [ equ : lambdamechanical ] ) with these force constants is in perfect agreement with the computer simulation results shown in fig . [ fig : absdisplacement ] .    for a system in which only nearest neighbor interactions are important , the force constants @xmath75 and @xmath76 can be simply related to the bulk modulus @xmath96 and the shear modulus @xmath33 . then , the force constant @xmath75 is given by @xmath97 where @xmath98 is the pair potential at distance @xmath99 . since in this case the elastic moduli are given by @xmath100 and @xmath101 one obtains @xmath102 to the extent that the response of the system to shear is determined by the interaction of neighboring parallel rows of particles , the energy density caused by shifting a whole row of atoms between two fixed ones is the same as that of a shear of appropriate magnitude . accordingly , the force constant @xmath76 is related to the shear modulus by @xmath103 this expression remains also valid if interactions beyond nearest neighbors are included between adjacent rows of particles . in terms of the elastic constants , the constant @xmath90 describing the exponential decay of the displacement field along the principal axis can be expressed as @xmath104 or , in terms of the poisson ratio @xmath105 , @xmath106 for a density of @xmath92 , inserting the poisson ratio of @xmath107 determined from a simple lattice sum yields @xmath108 , only slightly different from the correct value @xmath109 . this deviation occurs , because in the gaussian core model at the density @xmath92 interactions between non - nearest neighbor particles are important in determining the elastic constants ( in fact , considering only nearest neighbors would produce a negative shearing modulus @xmath33 in this case ) . for systems , in which only nearest neighbor interactions are relevant , the above expression is expected to hold accurately . point defects in two - dimensional crystals , such as interstitials and vacancies , can assume configurations with symmetries that vary from the symmetry of the underlying triangular lattice . while close to the defect the displacement field is highly anisotropic and strongly dependent on the atomistic details of the interactions , for large distances elasticity theory , which predicts isotropic behavior , is valid . for the particular i@xmath110 interstitial configuration , the displacement decreases exponentially with distance along the main defect axis . the decay constant is simply related to the material properties via the poisson ratio , which measures the ratio between transversal and axial strain upon stretching . in comparing the displacement fields computed from particles simulations with those obtained with continuum elasticity theory it is crucial to use equivalent boundary conditions in both cases . since particle simulations are usually carried out with periodic boundary conditions , also the differential equations of elasticity theory need to be solved for a periodic system . we have shown here that ewald summation , a technique routinely used in computer simulations to determine the electrostatic interactions of charges and dipoles , can be used for this purpose . in this method the sum over all interactions with periodic image defects is split into two sums in real space and reciprocal space , respectively . this particular treatment of the long - ranged nature of displacement fields effectively introduces a neutralizing background that leads to convergent sums . note that exactly the same expression apply also to a system that is enclosed in a rigid container . outside a core region near the defect , displacement patterns determined using such ewald summation agree perfectly with those calculated in particle simulations . the authors would like to thank christos likos , martin neumann , and andreas trster and for useful discussions . this research was supported by the university of vienna through the university focus research area _ materials science _ ( project `` multi - scale simulations of materials properties and processes in materials '' ) . 99                                                    in this appendix we calculate the response at @xmath4 of the one - dimensional bead - spring model of sec . [ sec : spring ] to a forced displacement of the first particle in the chain . the potential energy of the @xmath78 particles , located at positions @xmath79 , is given by @xmath111 where @xmath81 is the equilibrium distance and @xmath75 and @xmath76 are force constants . the vector @xmath112 includes the positions of all particles . minimizing the potential energy with respect to the particle positions @xmath79 by requiring that @xmath113 for all @xmath80 , one finds that at the potential energy minimum the particle positions are @xmath114 . we now displace particle @xmath83 by an amount @xmath84 from its original position @xmath115 and keep particle @xmath85 fixed at position @xmath116 . if we hold particle @xmath83 at this new position while minimizing the potential energy , all particles from @xmath117 to @xmath87 will move to new equilibrium positions . thus , the minimum energy configuration of the system is a function of the displacement @xmath84 of particle @xmath83 , which may be viewed as a parameter that is controlled externally and perturbs the system . to make this distinction between the displacement of particle @xmath83 and that of all other particles more explicit , we denote @xmath84 with an extra symbol , @xmath118 . the displacement @xmath89 of the particles @xmath119 is then a function of @xmath120 , @xmath121 where @xmath122 and @xmath123 denote the particle position in the minimum energy configuration with and without perturbation , respectively . in the following , we will calculate the displacements @xmath124 as a function of the perturbation strength @xmath120 . since condition ( [ equ : minpot ] ) defines the position of the energy minimum as a function of the perturbation strength @xmath120 , its derivative with respect to @xmath120 must vanish , @xmath125 application of the chain rule then leads to @xmath126 this condition must hold for all @xmath15 . defining @xmath127 and @xmath128 we can rewrite equ . ( [ equ : minpotderivationwithrespecttoa ] ) as @xmath129 inversion of the matrix @xmath130 then yields the vector @xmath131 , @xmath132 once @xmath133 is known , @xmath122 can be obtained by integration .    for the bead - spring model considered here , the first and second derivatives of the potential energy with respect to the particle coordinates are given by @xmath134 @xmath135 and @xmath136 to solve equ . ( [ equ : matrixnotationinverse ] ) we have to invert the symmetric tridiagonal matrix @xmath137 . for this particular matrix , the inverse matrix is known analytically @xcite , @xmath138}{\sinh(\alpha ) \sinh[(n+1)\alpha ] } \right . \nonumber \\ & & \hspace{1.0 cm } \left.-\frac{\cosh[(n + 1 - i - j)\alpha]}{\sinh(\alpha ) \sinh[(n+1)\alpha ] } \right\ } , \end{aligned}\ ] ] where @xmath139 since for our model @xmath140 , we obtain @xmath141 and hence @xmath142-\cosh[(n - i)\alpha]}{2\sinh(\alpha ) \sinh[(n+1)\alpha]}.\ ] ] for large @xmath87 , this equation simplifies to @xmath143 integration with respect to @xmath120 then yields @xmath144 where the integration constant is given by @xmath145 . thus , the displacement of particle @xmath80 is proportional to the displacement of particle @xmath83 and decays exponentially with the distance from the origin , @xmath146 with a decay constant @xmath90 that depends on the force constants @xmath75 and @xmath76 only .
Abstract:

Answer: point defects such as interstitials , vacancies , and impurities in otherwise perfect crystals induce complex displacement fields that are of long - range nature . in the present paper we study numerically the response of a two - dimensional colloidal crystal on a triangular lattice to the introduction of an interstitial particle . while far from the defect position the resulting displacement field is accurately described by linear elasticity theory , lattice effects dominate in the vicinity of the defect . in comparing the results of particle based simulations with continuum theory , it is crucial to employ corresponding boundary conditions in both cases . for the periodic boundary condition used here , the equations of elasticity theory can be solved in a consistent way with the technique of ewald summation familiar from the electrostatics of periodically replicated systems of charges and dipoles . very good agreement of the displacement fields calculated in this way with those determined in particle simulations is observed for distances of more than about 10 lattice constants . closer to the interstitial , strongly anisotropic displacement fields with exponential behavior can occur for certain defect configurations . here we rationalize this behavior with a simple bead - spring that relates the exponential decay constant to the elastic constants of the crystal .


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 dynamics of nonlinear systems is largely nontrivial and numerical simulations , involving sophisticated mathematical techniques , usually need to be made to fully understand their temporal evolution . the situation is even more complicated in spatially extended systems , like in nonlinear optics , in which ( multi)dimensional variables are present in the dynamics . in these cases an analytical ( or semi - analytical ) approach that allows us to gain physical insight about the system is only possible if we consider universal models @xcite , which formally capture the dynamics of nonlinear systems close to critical points ( like the threshold of emission in a laser ) . these models , also known as order parameter equations ( ope ) provide a simplified yet powerful description of the system . their universal character is due to the fact that very different systems ( biological , chemical , physical , etc . ) are described by the same ope , the only difference being in the meaning of variables and parameters under this approach . the symmetries of the system play a central role in the final form of the ope , those equations are able to capture qualitatively the essential elements of the dynamics . these ope can be real or complex depending on the fundamental variable ( electric field , temperature ) that is considered . in nonlinear optics , complex ope are commonly used since a complex ginzburg landau ( cgle ) for lasers for finite positive detuning was derived @xcite and have proven being very helpful for understanding a variety of systems @xcite . another equation , a complex swift - hohenberg ( cshe ) , which is valid for small detunings ( positive and negative ) was later obtained @xcite . one fundamental symmetry for complex ope describing spatially extended nonlinear systems is the phase symmetry , as it determines the nature of patterns which are possible in them . in systems with continuous phase symmetry , the dynamics is usually turbulent and the presence of vortices and spiral waves is common . when the systems only allows a finite number of phases ( by means a @xmath0 forcing @xcite for instance ) the dynamics become more ordered and controllable . furthermore , when a @xmath1 ( parametric ) forcing is applied to a systems with continuous phase symmetry , the dynamics of the system allows only two phases in it , which leads to the appearance of novel structures in the system like phase domains , domain walls and localized structures ( cavity solitons ) . the proper ope for the system must reflect this change and a generic equation can be derived for systems with @xmath0 forcing @xcite . for parametric forcing a family of ope , like the parametric complex ginzburg landau ( pcgle ) @xcite are obtained .    optical systems ( like lasers ) are usually insensitive to parametric forcing as the nonlinear response of these systems to high frequencies ( twice the natural frequency ) is negligible . in the last years a novel technique involving a new kind of forcing , known as `` rocking '' , was introduced @xcite . this is a @xmath2 forcing ( the frequency of the injection is close to the frequency of the system ) so it is appropriate for optical systems . the key factor is that the injection oscillates in time ( temporal rocking ) @xcite or space ( spatial rocking ) @xcite with a given frequency . this modulation modifies the dynamics of the system which becomes phase bistable . it can be shown that a pcgle describes the universal dynamics of a laser close to threshold under rocking when the undriven system exhibits an homogeneous hopf bifurcation . rocking has been successfully applied ( theoretically and experimentally ) to a wide range of systems @xcite . this has led to the derivation of opes describing those systems under certain limits and that provide relevant information of the dynamics . numerical simulations also prove that the behaviour predicted by these opes extends ( qualitatively ) far beyond the conditions imposed in their derivation , which increases the utility of this universal description . a complex swift - hohenberg equation with parametric gain , which describes a photorefractive oscillator ( pro ) under the injection of rocking was derived in @xcite . here we show that a very similar equation describes lasers with small detuning when rocking forcing is present , as well as we generalize the result of  @xcite to more general setups . the structure of the article is as follows . in section ii we present the equation which is object of our study and obtain their homogeneous solutions . in section iii a linear stability analysis is performed to both trivial and homogeneous solutions to study instabilities against perturbations of wavevector @xmath3 . in section iv we present some numerical simulations of our model for different values of parameters and show the patterns that can be found . we conclude in section v. in the appendices we present the derivation of the equation for lasers and pro and a comparison with the real swift - hohenberg equation . we consider the spatio - temporal dynamics of extended systems close to a bifurcation to traveling waves of long wavelength , when the system is forced in time close to its natural frequency ( 1:1 resonance ) and uniformly in space . however the type of forcing considered here is not the classic periodic one but has an amplitude which varies on time , specifically whose sign alternates periodically in time . the kind of forcing which we refer to can be thus expressed as:@xmath4 where @xmath5 is a @xmath6-periodic function of time of period @xmath7 , and @xmath8 is a high frequency , almost resonant with the natural frequency of the undriven system . a simple realization of @xmath9 is the function:@xmath10 where @xmath11 and @xmath12 is an amplitude .  according to ( [ forcing ] ) and ( [ periodic ] ) the forcing phase ( sign ) alternates in time . as shown in the appendices , where two relevant optical systems are analyzed under the action of such forcing , the state of the system can be expressed as:@xmath13 where @xmath14 is a @xmath6-periodic , spatially uniform contribution merely following the forcing , plus a spatially 2-dimensional field @xmath15 whose dynamics is governed by the following parametrically driven swift - hohenberg equation ( written in dimensionless form):@xmath16 the first line of this equation contains the usual complex swift - hohenberg equation , which models pattern formation arising from a finite wave number instability to traveling waves close to threshold @xcite , whereas in the second line we find the additional terms which appear when rocking is considered @xcite . the first extra term , @xmath17 , merely moves the reference frame to the central frequency of the rocking @xmath8 ( see appendices ) . the last extra term , @xmath18 , is the actual novelty as it breaks the phase invariance of the system , which becomes phase bistable as the equation only has the discrete symmetry @xmath19 . as for the parameters , @xmath20 measures the distance from threshold ( it can be removed from the equation by simple scaling but we keep it not to overwhelm the notation ) , @xmath21 controls the nonlinear frequency shift of the system , @xmath22 controls diffraction / dispersion . other parameters are @xmath23 , the detuning of the cavity from the natural frequency of the unforced system in the optical case , and @xmath24 , which is the detuning of the forcing from the natural frequency of the system . finally the `` rocking parameter '' @xmath25 is proportional to the squared amplitude of rocking @xmath26 and also depends on its frequency @xmath27 in a way whose exact form depends on the system considered ( see appendices ) ; note that when @xmath28 the effect of rocking is lost and ( [ swe1 ] ) becomes a usual complex swift - hohenberg equation .    in the following we will take @xmath29 as to not increase the degrees of freedom of the problem without need and also because in the analyzed cases that parameter happens to be zero . on the other hand @xmath30 in the pro case , while @xmath31 in the case of laser the spatially uniform nontrivial solutions of ( [ swe1 ] ) , or `` rocked states '' , can be expressed as @xmath32 where :    @xmath33       the state @xmath34 is always unstable as follows from a standard linear stability analysis , so we will not consider it in the following . on the other hand as @xmath35 can take two values ( differing by @xmath36 ) this produces two phase - locked states ( with same amplitude ) and phase differing by @xmath36 ( bistable phase locking ) . from now on we rename @xmath37 ( @xmath38 ) . the existence of these states requires @xmath39 and @xmath40 ; see ( [ rock ] ) . moreover , they exist only if @xmath41 , where :    @xmath42{c}\left\vert \theta\right\vert , \text { if } \left\vert \theta\right\vert < ( \mu-\delta^{2})/2\\ \gamma_{-}\text { otherwise}\end{array } \right . \text { \ \ \ \ \ \ } \nonumber\\ \text{\ } \gamma_{\pm }   &   = \frac{2(\mu-\theta^{2})}{3}\text{\ \ } \pm\frac { 1}{3}\sqrt{(\mu-\theta^{2})^{2}-3\delta^{2}}.\end{aligned}\ ] ]    on the plane @xmath43 the existence region is a closed domain ( `` rocking balloon '' ) as we can see in fig . 1 . outside that region we have the trivial solution and phase unlocked solutions ( traveling waves or patterns ) . [ ptb ]    fig1.jpg    we have performed a linear stability analysis of the trivial and uniform solutions against perturbations with wavevector @xmath3 . the eigenvalue with the largest real part ( for the stability of the trivial solution we just set @xmath44 ) reads:@xmath45   } \label{eigenval}\ ] ]    as we have two different nonlocal terms in the swift - hohenberg ( [ swe1 ] ) , @xmath46 will have , in general , two local maxima . one of them @xmath47 is associated with a real eigenvalue and it will give rise to static patterns through a patten - forming bifurcation ; the another one @xmath48 corresponds to a complex eigenvalue and it will produce oscillatory solutions ( homogeneous or traveling waves ) through a hopf bifurcation from the trivial solution . we could not obtain an exact analytical expression for @xmath49 but an approximated one by considering ( as it is observed numerically ) that it is close to @xmath50 . writing @xmath51 in ( [ eigenval ] ) and expanding @xmath52 to second order in @xmath53 , the resulting quadratic expression can  be maximized and solved for @xmath53 , leading to@xmath54 where @xmath55 for the trivial solution and @xmath56 in the case of the rocked states . the expression for @xmath57 is just @xmath58 so for @xmath59 we will have homogenous oscillations while for @xmath60 we will obtain traveling waves . [ ptb ]    fig2modi.jpg       all these results are summarized in fig . 2 ( positive @xmath23 ) and fig . 3 ( negative @xmath23 ) for @xmath61 ( laser case ) and @xmath62 . for positive @xmath63 if @xmath25 large , the trivial solution , which always exists , is stable . as we decrease the rocking parameter @xmath25 , we see bifurcations to other solutions like static patterns ( for the trivial state , @xmath64 ) and homogeneous oscillations ( for @xmath65 @xmath66 @xmath67 . the transition to uniform rocked states is supercritical in the upper bound ( the amplitude @xmath68 becomes @xmath69 along that line ) . regarding the lower bound ( two straight lines , where @xmath70 ) there is a saddle - node bifurcation which connects with oscillations by means of a complex eigenvalue for positive @xmath24 . for negative @xmath24 the rocked states become unstable close to the left edge of the balloon to instabilities ( real eigenvalue ) of wavenumber  @xmath71 . in fig.4 we can see an example of the temporal dynamics of spatially homogenous states close to the bifurcation : as we go closer to this , the period of the oscillations becomes larger and is infinite at the bifurcation . [ ptb ]    fig3modi.jpg    [ ptb ]    fig4.jpg    the analysis for negative @xmath23 is richer , as expected , because in this regime the swift - hohenberg equation ( [ swe1 ] ) without rocking ( @xmath28 ) already displays traveling waves of wavenumber @xmath72 ) . for large positive @xmath24 we still have a hopf bifurcation connecting trivial and traveling wave solutions ( for @xmath73 , @xmath74 , @xmath75 ) . for negative , and small positive @xmath24 ( @xmath76 ) the trivial solution destabilizes to a pattern with leading wavenumber @xmath71 at a certain @xmath77 . additionally , for small @xmath25 we have traveling waves as in the previous case . regarding the rocked states , these are destabilized in two ways : ( i ) through a complex eigenvalue with most unstable wavevector @xmath75 , which happens symmetrically around @xmath24 ( see fig . 3 ) and ( ii ) through a real eigenvalue with most unstable wavevector equal to @xmath78 , which happens , as in the case of positive detuning , close to the left side of the balloon ( here the unstable region is bigger ) . the presence of  these two unstable wavenumbers can be only seen in the simulations in the transient development of the instabilities ( fig . 5 ) as only one of the two spatial frequencies eventually survives . close to the lower bound we find ( i ) and close to the upper bound we get ( ii ) . in a range of values of @xmath79 we find that both instabilities arise simultaneously ( the two peaks of @xmath80 at @xmath81 @xmath82and @xmath83 becomes both positive as in fig . 6 : co - dimension-2 point ) , this happens for a small region in the @xmath43 plane ( fig.3 ) . this behaviour also appears ( for negative @xmath24 ) outside of the region where rocked states exist along a line which separates static patterns from traveling waves ( fig . 3 ) .    [ ptb ]    dobleanillo.jpg    [ ptb ]    fig6.jpg we studied the spatial patterns which appear outside the tongue where rocked ( phase - bistable ) states exist . in the former case , they basically confirm the previous analysis as we were able to obtain static spatial patterns which arise from real eigenvalue instabilities , selecting a particular spatial frequency as we have seen in the previous analysis . then the ( slow ) dynamics finally leads to a pattern where two spatial modes ( of opposite wavenumber ) are dominant , leading to roll - like patterns ( as in the right column of fig . 5 ) . additionally , we observe phase - unlocked patterns as ( slightly modulated ) traveling waves for negative @xmath23 ( fig . 7 ) , as a result of an instability governed by a complex eigenvalue . [ ptb ]    traveling.jpg    on the other hand , numerical experiments of ( [ swe1 ] ) confirm the existence of a variety of spatial patterns due to the phase bistability and the instabilities studied in the previous section . in the simulations we fix @xmath61 ( remember that @xmath30 for photorefractive oscillators and  @xmath31 for lasers ) as it can be shown that , after proper rescaling , the parameter @xmath22 does not affect the dynamics of the system and its value can be set an arbitrarily . different values of @xmath22 just change the temporal scale and the spatial scale associated to the real eigenvalue instability but what is relevant to the dynamics is the ratio between that scale and the one associated to the complex eigenvalue , which is determined by @xmath23 . [ ptb ]    dominiostran.jpg    traveling waves ( negative @xmath23 ) and homogenous oscillations ( positive @xmath23 ) are found in the phase - unlocked regions . inside the `` rocked '' balloon , two uniform rocked states of opposite phase can be connected through domain walls , generating phase domains ( fig . in two spatial dimensions , these domains are always a transient state before one phase dominates @xcite . close ( but still inside the stability region ) to the edge where the uniform rocked states lose their stability with real eigenvalue ( see section iii ) the walls become unstable due to curvature effects @xcite giving rise to the appearance of labyrinth patterns ( fig . 9 ) . before reaching the threshold where these labyrinths appear , we find dark - ring cavity solitons @xcite , which can be written / erased as it is shown in fig . this happens for both positive and negative @xmath84    [ ptb ]    labtrans.jpg    [ ptb ]    solitons.jpg we have proposed and studied , analytically and numerically , a complex swift - hohenberg equation with a parametric term that breaks the phase invariance , giving rise to phase - bistable patterns . this equation models diverse nonlinear optical cavities and can be thought of as a universal equation for rocked , spatially extended systems close to a homogeneous hopf bifurcation . in particular we have derived such an equation for two optical cavities :  the two - level laser and the photorefractive oscillator . the structure of the equation , with two nonlocal terms , produces two kinds of instabilities as revealed by a linear stability analysis of the homogeneous states . these two differ in the character of the eigenvalue governing the instability ( real or complex ) and in the spatial scale generating a complex variety of spatial patterns , phase locked and phase unlocked . the existence of extended patterns like ising domain walls and labyrinths is confirmed with numerical solutions as well as of dark - ring ( phase ) solitons . our starting point is the set of maxwell - boch equations with a bichromatical injected signal @xcite:@xmath85   + i\nabla^{2}e+f\cos(\omega t)e^{i\theta t}\label{mb}\\ \partial_{t}p   &   = -(1-i\delta)p+(r - n)e\nonumber\\ \partial_{t}n   &   = b\left [   -n+\frac{1}{2}\left (   e^{\ast}p+p^{\ast}e\right ) \right ] \nonumber\end{aligned}\ ] ]    the complex fields @xmath86 and @xmath87 are the scaled envelopes of the electric field and polarization , @xmath88 is proportional to the difference between the population inversion and its steady value in the absence of lasing . @xmath89 and @xmath90 , where @xmath91 , @xmath92 , and @xmath93 are , respectively , the decay rates of @xmath86 , @xmath87 , and @xmath94 . the transverse laplacian @xmath95 , where the spatial coordinates @xmath96 have been normalized so as to make unity the diffraction coefficient , and @xmath97 where @xmath6 is the physical time . @xmath98 is the pump parameter and the detuning @xmath99 , being @xmath100 ( @xmath101 ) the cavity ( atomic ) frequency . ( [ mb ] ) have been written in the frequency frame @xmath102 of the on - axis , or plane - wave @xmath103 , lasing solution in the absence of injected signal . in particular , this means that the actual injected field at the entrance face of the amplifying medium , @xmath104 , is of the form:@xmath105    we will consider two types of scales : fast and slow .      regarding the scales , we will consider that the detuning @xmath23 and laplacian @xmath106 are @xmath107 @xmath108as we consider class c lasers , @xmath109 and @xmath110 are @xmath111 additionally , @xmath112 and @xmath113 time scales are @xmath114 and @xmath115 , the pump is @xmath116 ( we are close to the threshold)@xmath108 the variables will be written in the form:@xmath117    @xmath118    at this order  @xmath119    @xmath107    this is the first nontrivial order and reads:@xmath120 where@xmath121{cc}-\sigma & \sigma\\ r_{0 } & -1 \end{array } \right )   , \ ] ]    and we have introduced the notation@xmath122{c}e_{i}\\ p_{i}\end{array } \right )   , \quad i=1,2,3,\ldots\label{vi}\ ] ] eq . ( [ l0v1c1 ] ) can be easily solved with the help of the left eigenvectors of matrix @xmath123:@xmath124 projecting eq . ( [ l0v1c1 ] ) onto @xmath125 we obtain @xmath126 , and projecting onto @xmath127 we obtain@xmath128 hence@xmath129{c}1\\ 1 \end{array } \right )   . \label{v1c1}\ ] ]    @xmath130    making use of eqs . ( ) and ( ) , we obtain@xmath131    where  @xmath132{c}i\nabla^{2}e_{1}-i\sigma\delta e_{1}+f_{1}\cos\left (   \omega t_{1}\right ) e^{i\theta t_{2}}\\ i\delta e_{1}\end{array } \right)\ ] ]    projecting eq . ( [ dt1v1c1 ] ) onto @xmath125 and making use of eq . ( [ p1c1 ] ) we obtain@xmath133 we solve@xmath134    projecting onto @xmath127 and making use of eq . ( [ dt1e1c1 ] ) we obtain@xmath135    @xmath136    @xmath137    where@xmath138{c}0\\ \left (   r_{2}-\left\vert e_{1}\right\vert ^{2}\right )   e_{1}+i\delta p_{2}\end{array } \right)\ ] ]    projecting@xmath139   = \nonumber\\ -\left (   1+\sigma\right )   \frac{\partial e_{1}}{\partial t_{2}}+\sigma\left ( r_{2}-\left\vert e_{1}\right\vert ^{2}\right )   e_{1}\ ] ]    we can rewrite this as following@xmath140    where@xmath141    the solution can be written formally as:@xmath142 to ensure convergence it must be true that:@xmath143    taking into account that if @xmath144 is homogeneous then it does not depend on @xmath145 the following condition must be fulfilled ( we make use of ( [ ea1])):@xmath146    in order to avoid divergences ( in the case @xmath147 )    we finally compute the time derivative @xmath148 up to third order in @xmath53 which , making use of ( [ ea2 ] ) and ( [ ea1 ] ) , setting @xmath149 and rescaling to the original variables , we obtain:@xmath150      the scales in this case will be : @xmath151    @xmath152 @xmath153@xmath154{cc}-\sigma & \sigma\\ 1 & -1 \end{array } \right)\ ] ]    the left - eigenvectors are @xmath155 and @xmath156 with eigenvalues @xmath69 and @xmath157 respectively . we can use this to write:@xmath158{c}e_{1}\\ p_{1}\end{array } \right)\end{aligned}\ ] ]    we obtain @xmath159    @xmath160@xmath161 @xmath132{c}i\nabla^{2}e_{1}-i\sigma\delta e_{1}\\ i\delta e_{1}\end{array } \right)\ ] ]    using the same procedure as above , we obtain ( if we set @xmath162 ) : @xmath163    @xmath164@xmath165{c}f\cos\left (   \omega t_{1}\right )   e^{i\theta t_{2}}\\ ( r_{2}-\left\vert e_{1}\right\vert ^{2}e_{1}+i\delta p_{2}\end{array } \right)\end{gathered}\ ] ]    now we obtain the following equation:@xmath166 using ( [ p2 ] ) and undoing the scaling:@xmath167    we can write:@xmath168    setting @xmath169 the previous equation becomes@xmath170    if the frequency @xmath27 is high compared with the dynamics of the system we can set : @xmath171 @xmath172 with @xmath173 and as it is done in @xcite separate the slow dynamics from the fast dynamics :    specifically @xmath174@xmath175    no we solve at different orders    @xmath176@xmath177    @xmath178    the equation reads at this order:@xmath179    this can be written as @xmath180 which can be solved : @xmath181 . the boundness of this requires:@xmath182    therein , we obtain the following solvability condition:@xmath183    this condition is exactly the same equation as [ fast]@xmath184so we recover the result that we obtained considering fast scales by just assuming that @xmath185 and @xmath27 are large . so , independently of the set of scales we consider , the final equation is consistent . setting @xmath186 , @xmath187 , @xmath188 , @xmath189 , @xmath190 and defining @xmath191 we can write ( removing the commas for simplicity):@xmath192 our starting point is the set of equations as in @xcite in which we consider a bichromatical injection and we make the change ( @xmath193 for convenience . after removing the commas , we get:@xmath194    where @xmath195  is the slowly varying envelope of the intracavity field , @xmath196 is the complex amplitude of the photorefractive nonlinear grating , * * r**=(x , y ) are the transverse coordinates , @xmath197 is the product of the cavity linewidth @xmath198 with the photorefractive response time @xmath199 ( @xmath200 under typical conditions , and @xmath201 s ) , @xmath202 is time measured in units of @xmath199 , the detuning ( @xmath203 ( @xmath204 and @xmath100 are the frequencies of the pump and its nearest longitudinal mode , respectively ) , @xmath205 is the ( real ) gain parameter that depends of the crystal parameters and the geometry of the interaction . the scales will be @xmath206 ( results are independent of the specific scale provided that @xmath109 very large ) and @xmath207 which implies @xmath208 . we also assume that @xmath205 is close to threshold , @xmath209 and that @xmath210 we will have `` two times '' , fast and slow : @xmath211 so that @xmath212we also set:@xmath213                      we solve:@xmath225   + \nonumber\\ \left [   -\partial_{t_{2}}\varphi_{2}+(g_{2}-1)\varphi_{1}-(-i\delta + i\nabla^{2})^{2}\varphi_{1}+i\nabla^{2}\varphi_{2}-\left\vert \varphi _ { 1}\right\vert ^{2}\varphi_{1}\right ]   -\nonumber\\ \left [   \left\vert e_{1}^{\omega}\right\vert ^{2}e_{1}^{\omega}+2e_{1}^{\omega}\left\vert \varphi_{1}\right\vert ^{2}+(e_{1}^{\omega})^{\ast}\varphi_{1}^{2}+2\left\vert e_{1}^{\omega}\right\vert ^{2}\varphi _ { 1}+\left\vert e_{1}^{\omega}\right\vert ^{2}\varphi_{1}^{\ast}\right ] \nonumber\end{gathered}\ ] ]
Abstract:

Answer: we study pattern formation in a complex swift hohenberg equation with phase - sensitive ( parametric ) gain . such an equation serves as a universal order parameter equation describing the onset of spontaneous oscillations in extended systems submitted to a kind of forcing dubbed rocking when the instability is towards long wavelengths . applications include two - level lasers and photorefractive oscillators . under rocking , the original continuous phase symmetry of the system is replaced by a discrete one , so that phase bistability emerges . this leads to the spontaneous formation of phase - locked spatial structures like phase domains and dark - ring ( phase- ) cavity solitons . stability of the homogeneous solutions is studied and numerical simulations are made covering all the dynamical regimes of the model , which turn out to be very rich . formal derivations of the rocked complex swift - hohenberg equation , using multiple scale techniques , are given for the two - level laser and the photorefractive oscillator