Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: in 1967 - 73 , the four _ vela _ spacecraft ( named after the spanish verb _ velar _ , to keep watch ) , that where originally designed for verifying whether the former soviet union abided by the limited nuclear test ban treaty of 1963 , observed 16 peculiarly strong events ( klebesadel , olson and strong 1973 , bonnell and klebesadel 1996 ) . on the basis of arrival time differences , it was determined that they were related neither to the earth nor to the sun , but they were of cosmic origin . therefore they were named cosmic @xmath0-ray bursts ( grbs hereafter ) . grbs appear as brief flashes of cosmic high energy photons , emitting the bulk of their energy above @xmath1 0.1 mev ( fig . they are detected by instruments somewhat similar to those used by the particle physicists at their laboratories . the difference is that grb detectors have to be placed onboard balloons , rockets or satellites . in spite of the abundance of new observations of grbs , their energy source and emission mechanism remain highly speculative . the konus experiment on _ veneras 11 _ and _ 12 _ gave the first indication that grb sources were isotropically distributed in the sky ( mazets et al . 1981 , atteia et al . 1987 ) . based on a much larger sample , this result was nicely confirmed by batse on board the _ cgro _ ( meegan et al . 1992 ) . about 300 grbs occur annually in the full sky , but only few of them are localized accurately . the apparent isotropy was interpreted in terms of grbs produced at cosmological distances , although the possibility of a small fraction of the sources lying nearby , within a galactic disc scale of few hundred pc , or in the halo of the galaxy , could not be discarded . another result was that the time profiles of the bursts are very different , with some grbs lasting a few ms and others lasting for several minutes . in general , there was no evidence of periodicity in the time histories of grbs . however there was indication of a bimodal distribution of burst durations , with @xmath225% of bursts having durations around 0.2 s and @xmath275% with durations around 30 s. a extensive review of the observational characteristics can be found in fishman and meegan ( 1995 ) .    a deficiency of weak events was also noticed , and all these observational data led many researchers to believe that grbs are indeed at cosmological distances . in this case , the released energies could be as high as 10@xmath3 erg and models could involve coalescence of neutron stars in double systems , neutron star - black hole systems , accretion induced collapse in white dwarfs or @xmath4failed@xmath4 type i supernovae . it was also proposed the possibility of having a mixture of two populations : a cosmological plus a galactic one , but the latter , probably formed by accreting neutron stars , will account for only a very small fraction of the total population . see nemiroff ( 1984 ) for a review of the different theoretical models . it is well known that an important clue for resolving the grb puzzle is the detection of transient emission -at longer wavelengths- associated with the bursts . here i review all the efforts in the search for grb counterparts throughout the electromagnetic spectrum . i will first review the searches prior to 1997 , and afterwards i will discuss the important discoveries achieved last year . previous reviews can be seen in schaefer ( 1994 ) , hartmann ( 1995 ) , vrba ( 1996 ) , greiner ( 1996a ) and hurley ( 1998 ) . the observation of x - rays below 20 kev is important in order to understand the nature of the physical processes in the grb sources . most of the experiments devoted to grb study had lower cut - off energies of the order of 20 - 40 kev . however , a few of them were designed in order to cover part of the soft x - ray range .    already in 1972 , the apollo 16 spectrometer detected x - ray photons in the 2 - 8 kev range for grb 720427 ( trombka et al . wheaton et al . ( 1973 ) observed in grb 720514 the presence of an x - ray signal which lasted _ longer _ than the gamma - ray event . in 1984 , katoh et al . reported an x - ray tail for grb 811016 that lasted more than 30 s after the termination of the hard event . _ ginga _ also detected emission down to 1.5 kev , with soft x - ray tails found for 8 events ( murakami et al . 1991 , yoshida and murakami 1994 ) . moreover , a significant fraction ( @xmath2 10% ) of the grbs detected by watch on _ granat _ displayed x - ray tails in the 6 - 20 kev range ( castro - tirado et al . connors and mcconnell ( 1996 ) have discussed the proposal that the x - ray transient detected by _ heao 1 _ following grb 780506 could indeed be related to the burst itself . therefore it seemed clear that , at least for some bursts , there was evidence for x - ray emission lasting longer than the grb . grindlay et al . ( 1982 ) and pizzichini et al . 1986 initiated the searches for x - ray quiescent counterparts by means of _ einstein_. five grb error boxes were observed and , with one possible exception ( grb 781119 , grindlay et al . 1982 ) , no point source was detected in any of the observations . assuming that the bursts originated from galactic neutron stars , only upper limits to the surface temperatures were derived ( pizzichini et al . boer et al . ( 1988 ) observed the grb 781119 region , but failed to detect the source with _ exosat_. recently , _ asca _ and _ rosat _ `` rediscovered '' the x - ray source in grb 781119 , with a flux of f@xmath5 1.6 @xmath6 10@xmath7 erg @xmath8 s@xmath9 ( hurley et al . 1996 , boer et al . further evidence for an x - ray source associated with a grb was obtained by _ asca _ in the grb 920501 error box ( murakami et al . the x - ray spectrum is heavily absorbed at low energies , implying an absorption consistent with the galactic value in that direction , and suggesting that the x - ray source is _ quite _ distant . using the _ rosat _ all - sky - survey data , greiner et al . ( 1995 ) detected 27 x - ray sources in 45 small grb error boxes . on the basis of the optical identification of theses sources , no convincing association to grbs was found . moreover , the number of sources corresponded to the expected number of background sources . the upper limit to the x - ray flux was f@xmath10 10@xmath7 erg @xmath8 s@xmath9 . another field ( grb 940301 ) was observed by _ rosat _ four weeks later and revealed no fading counterpart ( greiner et al . 1996b ) . searches at euv wavelengths have only been performed since 1992 , following the launching of _ euve_. castro - tirado et al . ( 1998a ) has reported a serendipitous observation in the extreme ultraviolet window by _ euve _ , @xmath2 11 hours after the occurrence of grb 921013b . no source was detected in the euv range . a deep observation on the grb 920325 error box was performed 17 months after by _ euve_. no quiescent source was found , thus constraining a neutron star as a possible counterpart ( hurley et al . 1995 ) . the strategy behind this kind of search is the assumption that grbs repeat . some scientists have assumed that grbs could be a repetitive phenomenon , and they have used all the available information at the major plate archives in order to look for optical transient ( ot ) emission in some of the smaller grb error boxes . atteia et al . ( 1985 ) searched for about 1500 hours but no such transients were seen . schaefer discovered in 1981 the first ot in a harvard plate taken in 1928 ( ot 1928 ) , located inside the small grb 781119 error box . another ot was found in the grb 781006b error box in three different plates taken in 1966 . a large amplitude flaring dme star , probably unrelated to the burst , was proposed to be related to the ot ( greiner and motch 1995 ) . another good candidate is ot 1905 , found by hudec et al . ( 1994 ) slightly outside the interplanetary network error box for grb 910219 . although these ots seem to be real , another @xmath2 40 candidates were found , but most of them were rejected as they turned out to be plate defects . see a discussion in hudec ( 1993 , 1995 ) . however , there is no convincing proof that the ots found in plate archives and grbs are related to the same physical phenomenon . deep searches were carried out in the 80 s for the smallest grb error boxes provided by the interplanetary network : grb 781119 ( studied by pedersen et al . 1983 ) and grb 790406 ( chevalier et al . recently , with the advent of the newer big ccd detectors , deep searches have been carried out of grb 790329 and other error boxes ( vrba , hartmann and jennings 1995 ) . the most powerful telescopes , such as the 6-m in the caucasus and the hubble space telescope ( hst ) have also been used for observing the error boxes of grbs : grb 790613 ( sokolov et al . 1996 ) and grb 790325 and 920406 ( schaefer et al . 1997a ) . in any case , no quiescent grb counterpart has been firmly established . in general , the small grb error regions did not contain galaxies to fairly deep limits ( the so - called @xmath4no - host problem@xmath4 , fenimore et al . 1993 ) , and it was obvious that the bursts , if at cosmological distances , could not come from normal galaxies . either they had to occur in subluminous galaxies or their progenitors somehow had to be ejected from the parent galaxy . the first search for optical emission arising _ simultaneously _ with a grb was performed by grindlay , wright and mccrowsky ( 1974 ) . hudec et al . ( 1987 ) and hudec ( 1993 ) found five events for which time - correlated plates were available . greiner et al . ( 1996c ) found plates within @xmath11 12-h of the bursts , but none revealed interesting brightenings in the error box . for grb 920824 , a rather deep plate was taken simultaneously with the burst in dushanbe . they conclude that the optical emission of typical grbs is at level below @xmath12 @xmath2 2 _ during _ the burst and @xmath12 @xmath2 20 - 400 a few hours _ after _ the burst , with @xmath13 and @xmath14 being the peak gamma - ray ( and optical ) fluxes . stronger upper limits have been obtained by grosce ( lee et al . 1996 ) who have recorded sky images for more than 30 grb triggers , in many cases while the burst was still in progress . additional results were obtained by castro - tirado et al . ( 1994b ) regarding follow - up observations of eight bursts detected by watch . no transient optical emission was detected in the following hours down to 18th magnitude . similar observations were reported for grb 930131 by schaefer et al . ( 1994 ) and mcnamara and harrison ( 1994 ) . again , no firm counterpart was detected in their observations . barthelmy et al . ( 1994 ) obtained upper limits for any transient optical emission arising from three bursts , between 35 hr and 8 days following the burst . the only ot found nearly - simultaneously to a grb event , is the one reported by borovicka et al . they found a star - like spot at the edge of the grb 790929 error box on a plate taken only 7.1 hr after the high energy event . also , two additional objects were found at the same position on another plates . the objects are consistent with the star hde 249119 . further optical studies detected additional optical eruptions of hde 249119 , reporting ten optical brightenings with amplitudes greater than 0.5 mag ( stpn & hudec 1996 ) , but it is unlikely that this star would be related to grb 790929 . the first search for counterparts in the near - ir was initiated by schaefer et al . ( 1987 ) , who investigated seven grb error boxes at a wavelength of 2.2 @xmath15 m ( k - band ) . no convincing candidates were found , with a detection limit of k = 19 for two error boxes . an imaging survey of the six smallest grb error boxes from the third interplanetary network has revealed relatively bright galaxies in each of them ( larson , mclean and becklin , 1996 ) . although the statistical significance of their result is not overwhelming , it leaves the `` no - host problem '' as an open question . klose et al . ( 1996 ) obtained deep imaging of the grb 790418 error box , showing several galaxies as well . deep infrared imaging of the quiescent x - ray source located within the grb 920501 error box has been performed by blaes et al . one of the objects is a type 1 seyfert galaxy , and drinkwater et al . ( 1997 ) has identified the galaxy with the x - ray source seen by _ einstein _ , _ asca _ and _ rosat _ , strengthening the above - mentioned results . is this seyfert galaxy the true grb counterpart ? searches for _ quiescent _ counterparts have also been conducted in the past in the ir . schaefer et al . ( 1987 ) used the iras data base at wavelenghts of 12 , 25 , 60 and 100 @xmath15 m and looked for candidates within 23 well localized grb error boxes . no convincing counterparts were found . two simultaneous observations were obtained by the cobe satellite , when correlating grbs that occurred during an 8 month overlap period ( bontekoe et al . 1995 , stacy et al . 1996 ) . the 2@xmath16 upper limit quoted is 31000 jy at 53 ghz . a likelihood analysis was performed by ali et al . ( 1997 ) in order to look for a change in the level of millimetric emission from the location of 81 grbs during the first two years of the mission ( 1990 - 91 ) . no positive detection was achieved , with 95 % confidence level upper limits of 175 , 192 and 645 jy respectively at 31 , 53 and 90 ghz . another team at the fallbrook low - frequency immediate response telescope ( flirt ) searched for prompt radio emission at 74 mhz arising from grb sources . the delays were of the order of 2 to 12 min . preliminary results for 12 grbs indicated that no source was found above a @xmath2 100 jy detection limit ( balsano et al . 1996 ) . a serendipitous observation towards the error box of grb 920711 was performed at 151 mhz with the cambridge low - frequency synthesis telescope ( clfst ) from 2 to 4 days after the occurrence of grb 920711 . no source was detected down to 40 mjy ( koranyi et al . this instrument also observed grb 940301 within 1 hour after the burst , imposing a 200 jy upper limit , and 80 mjy between 1 - 36 days after the event . clfst observations were conducted for another two bursts , with flux limits of 35 jy for 2.5 hours after grb 950430 and 16 jy for 1 hour after grb 950706 ( green et al . 1996 ) . a search for a radio transient at 4.86 and 1.4 ghz within the grb 930131 failed to detect any source ( schaefer et al . 1994 ) . vla observations were performed 23 hours after grb 940217 , but only an upper limit of @xmath17 4 @xmath15jy / beam was given by palmer et al . no fading / flaring radio counterpart was found either at 1.4 ghz ( upper limit of 3.5 mjy ) or 0.4 ghz ( upper limit of 55 mjy ) within the grb 940301 error box ( frail et al . 1994 ) . it was suggested that if grbs were to be related to highly magnetized , rotating neutron stars , they would produce polarized radio emission during the quiescent state . in fact , vla maps of the grb 781119 error box showed two sources ( hjellming and ewald 1981 ) . one showed some degree of polarization , and the other coincided with the _ rosat _ x - ray source ( hurley et al . 1996 ) . schaefer et al . ( 1989 ) found four radio sources in 10 small grb error regions . none of them displayed unusual properties and the number of sources was the expected number of chance background sources . they concluded that none of them were associated with grbs . deep vla observations of three grbs were also performed by palmer et al . grb 920501 was observed 13 days and 9 months after the burst , at frequencies of 1.4 and 8.4 ghz , imposing upper limits of @xmath17 235 @xmath15jy / beam and @xmath17 77 @xmath15jy / beam respectively . a strong limit was derived for grb 930706 as @xmath17 1 @xmath15jy / beam after 9 days . photons with energies up to 18 gev have been detected in only one case , 1.5-hr after grb 940217 ( hurley et al . there are models that predict tev emission from grbs ( meszaros , rees and papathanassiou 1994 , halzen et al . a positive detection would constrain models with respect to the emission mechanism and the distance , because no tev photons are expected if they occur at cosmological distances due to the interaction of the photons with the microwave and infrared background . this process leads to pair production and might become important at e @xmath18 30 tev , as pointed out by wdowczyk , tkaczyk and wolfendale ( 1992 ) , mannheim , jartmann and funk ( 1997 ) and stecker and de jager ( 1997 ) . in fact , at e @xmath18 50 tev , there seems to be evidence of a detection from the dublin asa ( grb 910511 , plunkett et al . the rest of the current arrays and cerenkov telescopes have provided only upper limits at this energy threshold . a systematic search was conducted at aes - kgf ( bhat et al . 1996 ) for pev @xmath0-rays from the locations of about 39 high fluence grbs during apr 91-mar 93 . upper limits of 2 @xmath6 10@xmath19 photons @xmath8 s@xmath9 were obtained . another search was performed for 115 grbs at the casamia detector . flux limits of 6 @xmath6 10@xmath19 photons @xmath8 s@xmath9 at e @xmath18 160 tev were obtained ( catanese et al . 1996 ) . follow - up searches have also been carried out with hegra above 1 tev for more than 150 bursts , but no significant excess was detected by krawczynski et al . ( 1996 ) or padilla et al . in the 0.25 - 4 tev range , it has been suggested that the whipple observatory can detect , under favourable conditions , @xmath2 1 event per year ( only @xmath2 0.01 per year in the case of the hegra array ) . however , only upper limits have been set up by the whipple telescope ( connors and mcconnell 1996 ) and the eas - top asa experiment at e @xmath18 10 tev ( agglietta et al . 1995 ) . the search for neutrinos , antineutrinos and muons arising either from intrinsic grb mechanisms as proposed by mszaros and rees ( 1993 ) and paczynski and xu ( 1994 ) or earthrock neutrino - produced muons , has been performed at several places , but none of the experiments have detected neutrinos or muons in excess of the expected background in correlation with a grb . these are the mount blanc lsd telescope ( agglietta et al . 1995 ) , the irvine - michigan - brookhaven ( imb ) detector ( losecco 1994 , becker - szendy 1995 ) and the soudan-2 muon experiment ( de muth et al . upper limits are 6 @xmath6 10@xmath20  @xmath8 ( for @xmath21 , @xmath22 , @xmath23 at 20 @xmath17 e@xmath24 @xmath17 100 mev ) and 2.1 @xmath6 10@xmath25  @xmath8 ( for @xmath26 , @xmath27 , @xmath28 at 20 @xmath17 e@xmath24 @xmath17 100 mev ) ( agglietta et al . with the launching of _ bepposax _ , it was possible on 28 feb 1997 to detect the first _ clear _ evidence of a long x - ray tail -the x - ray afterglow- following grb 970228 . see figure 2 . the x - ray fluence was @xmath2 40 % of the gamma - ray fluence , as reported by costa et al . therefore , the x - ray afterglow is not only the low - energy tail of the grb , but also a significant channel of energy dissipation of the event on a completely different timescale . another important result was the non - thermal origin of the burst radiation and of the x - ray afterglow ( frontera et al . 1998 ) . further x - ray afterglows were observed by _ bepposax _ in grb 970402 ( piro et al . 1997a ) , grb 970508 ( piro et al . 1997b ) , grb 971214 ( antonelli et al . 1997 ) and probably in grb 970111 ( feroci et al . 1998 ) and grb 971227 ( sofitta et al . 1997 ) . another three were observed by _ rossixte _ in grb 970616 ( marshall et al . 1997a ) , grb 970815 ( smith et al . 1997 ) and grb 970828 ( marshall et al . 1997 , remillard et al . the latter was also detected by _ ( murakami et al . 1997 ) and _ rosat _ ( greiner et al . 1997 ) . the x - ray spectrum as observed by _ is strongly absorbed ( n@xmath29 = 10@xmath30 @xmath8 , murakami et al . 1998 ) , suggesting that the event occurred in a dense medium . bor et al . ( 1997b ) reported an observation by euve @xmath2 20 hr after grb 971214 . no source was detected at 100 @xmath31 . assuming n@xmath29 = 10@xmath32 @xmath8 , the upper limit to the flux is 1.7 @xmath6 10@xmath7 erg @xmath8 s@xmath9 . the improvement of grb localizations by _ bepposax _ ( in combination with the _ third interplanetary network _ of satellites ) has led to the detection of the first three optical counterparts in 1997 . the first optical transient associated to a grb was identified on 28 feb 1997 , 20 hr after the event , by groot et al . the ot was seen in earlier images taken by pedersen et al . ( 1997 ) and guarnieri et al . ( 1997 ) , in the rising phase of the light curve . the maximum was reached @xmath2 20 hr after the event ( v @xmath2 21.3 ) , and followed by a power - law decay as t@xmath33 ( galama et al . 1997 , masetti et al . an extended source was seen at the ot position since the very beginning by both ground - based and _ hst _ observations ( van paradijs et 1997 , sahu et al . 1997 ) . hst _ observations taken 6 months after the event were reported by fruchter et al . both the ot ( at v = 28 ) and the extended source ( v = 25.6 ) were seen . the extended source surrounding the point - source was interpreted as a galaxy , according to the similarities ( apparent size , magnitude ) with objects in the _ hst _ deep field . the second optical transient is associated with grb 970508 and was discovered by bond ( 1997 ) . it was observed only 3 hr after the burst in unfiltered images ( pedersen et al . 1998a ) and 4-hr after in both the u and r - bands ( fig . it displayed a strong uv excess ( castro - tirado et al . 1998 , galama et al . the optical light curve reached a peak two days after ( r = 19.7 , djorgovski et al . 1997 ) and was followed by a power - law decay as t@xmath33 ( fig . the long - lasting ot was observed at many ground - based observatories ( chevalier et al . 1997 , garcia et al . 1997 , mignoli et al . 1997 , schaefer et al . 1997b , sokolov et al . optical spectroscopy was performed at calar alto and la palma ( castro - tirado et al . 1997a ) and hawaii ( metzger et al . 1997a ) . the latter spectrum allowed a direct determination of the grb 970508 redshift ( @xmath34 ) and was the first proof that at least a fraction of the grb sources lie at cosmological distances ( metzger et al . unlike grb 970228 , no host galaxy has been seen so far , but the optical spectroscopy and the flattening of the decay in late august ( pedersen et al . 1998b ) suggest that the contribution of a constant brightness source -the host galaxy- has been detected . the luminosity of the galaxy is well below the knee of the galaxy luminosity function , @xmath35 , and the detection of deep mg i absorption and strong [ o ii ] emission -both indicative of a dense medium- suggests that the host could be a dwarf , blue , rapidly forming starburst galaxy ( pian et al . 1998 ) . the third optical transient is related to grb 971214 and was identified by halpern et al . ( 1997 ) as an i = 21.2 object that faded 1.4 mag in 1 day . independently , this object was also noticed as suspicious by itoh et al . further observations proved that the decay followed a power - law decline ( diercks et al . 1997 , castander et al . 1997 ) with t@xmath36 , similar to , but steeper than , the other two grbs with optical counterparts .    in spite of intensive searches , no optical counterparts were found within one day of the event for grb 970111 , grb 970402 or grb 970828 ( castro - tirado et al . 1997b , pedersen et al . 1998b , groot et al . 1998 ) . fading near  ir counterparts were observed for grb 970228 , grb 970508 and grb 971214 . for grb 970228 , the ir - flux declined from j = 21.0 ( k @xmath18 19.5 ) on mar 18 ( klose and tuffs 1997 ) to j = 23.5 ( k = 22.0 ) on mar 30 ( soifer et al . 1997 ) . in the case of grb 970508 , the ir flux decayed from k = 18.4 to k = 19.2 in three days ( morris et al . 1997 ) . grb 971214 was detected in the k - band 4 hr after the burst ( gorosabel et al . 1998 ) and decreasing from j = 20.3 on dec 15 to j = 21.5 on dec 16 ( tanvir et al . 1997 ) .      following the detection by bepposax of grb 970402 , a target - of - opportunity observation by _ iso _ was performed 55-hours after the event , in the isocam 8 - 15 @xmath15 m ( iras - equivalent bandpass ) and isophot 174 @xmath15 m band , in order to monitor the entire bepposax grb error box . a further observation was carried out 8 days later , with a similar set - up , in the x - ray afteglow error box . about 50 sources were detected in a field centered on the grb error box . seven of them lie within the smaller x - ray error circle . with the exception of the variable star bl cir , no other variable object was found in the entire grb error box . 5-@xmath16 upper limits for any new object are f@xmath37 @xmath17 0.14 mjy and f@xmath38 @xmath17 350 mjy ( castro - tirado et al . a negative result was also obtained in the same bandpass for grb 970508 ( hanlon et al . 1998 ) . grb 970508 was detected as a continuum point source at 86 ghz with the iram plateau de bure interferometer . a @xmath2 2 mjy flux was observed on 19 - 21 may 1997 . the object was not detected any more after 28 may . no significant intra - day variation was seen , and it was excluded that the flux was modulated significantly by interstellar scintillation ( bremer et al . 1998 ) . in the case of grb 971214 , no source was detected at 850 @xmath15 m at the james clerk maxwell telescope . the upper limit was 1 mjy on dec 17 - 22 ( smith and tilanus 1998 ) . after two unsuccesful searches for transient radio emission in grb 970111 and grb 970228 ( frail et al . 1997a ) , prompt vla observations of the grb 970508 error box allowed detection of a variable radio source within the error box , at 1.4 , 4.8 and 8.4 ghz . it coincided with the optical transient and displayed a bizarre behaviour , in the sense that it remained at a steady level of brightness , flaring occasionally by a factor of two or three ( frail et al . 1997b ) , supporting the idea of the host galaxy as an active galactic nucleus . vlbi observations did not resolve the source ( taylor et al . frail et al . proposed that the fluctuations could be the result of strong scattering by the irregularities in the ionized galactic interstellar gas . the damping of the fluctuations with time indicates that the source expanded to a significantly larger size . the source was also detected at 15 ghz ( pooley and green 1997 ) . all these observational efforts have provided the complete picture of a grb multiwavelength spectrum , in the case of grb 970508 ( fig . the observational characteristics of the three grb counterparts ( or _ hypernovae _ , a name proposed by paczynski ) can be accommodated in the framework of the fireball models ( mszros et al . 1994 , mszros and rees 1997 , tavani 1997 , vietri 1997 , waxman 1997 , wijers et al . 1997 ) , in which a compact source releases 10@xmath3 ergs of energy within dozens of seconds . the plasma accelerates to relativistic velocities ( the fireball ) . the blast wave is moving ahead of the fireball , and sweeps up the interstellar matter , producing an afterglow at frequencies gradually declining from x - rays to visible and radio wavelenghts . although some of the predicted behaviour has been observed , there are still many open questions ( see katz and piran 1998 ) . the existence of an x - ray afterglow in _ all _ bursts seems to be confirmed . the first optical / ir counterparts have been found in 1997 . grb 970508 was also detected in radio and mm , but prompt searches at other wavelengths failed to detect grb 970111 , grb 970402 and grb 970828 . _ bepposax _ and _ rossixte _ have opened a new window in the grb field and it is widely accepted now that most grbs , if not all , lie at cosmological distances . it is expected that _ bepposax _ , _ rossixte _ and _ cgro _ will facilitate the discoveries of other counterparts and , together with the new high - energy observatories ( _ axaf _ , _ spectrum x/@xmath39 _ , _ xmm _ , _ integral _ , _ hete 2 _ ) and other satellites of the future _ 4th interplanetary network _ , will definitively solve the long - standing gamma - ray burst mystery . i am grateful to many colleagues for very fruitful discussions , in particular to s. brandt , m. cervio , e. costa , m. feroci , g. fishman , f. frontera , s. golenetskii , j. greiner , j. gorosabel , s. guzyi , r. hudec , k. hurley , j. isern , c. kouveliotou , n. lund , l. metcalfe , h. pedersen , m. de santos , a. shlyapnikov , r. sunyaev , w. wenzel , b. wilson , c. wolf , w. wenzel , a. yoshida and many others . i should also mention here the continuous support of my wife , m. e. alcoholado - feltstrm , who always patiently awaits for my return when sometimes i run off to the nearest observatory in order to perform the grb follow - up observations and further 40-hr non - sleeping data analysis .
Abstract:

Answer: since their discovery in 1973 , gamma - ray bursts ( grbs ) have remained for many years one of the most elusive mysteries in high energy - astrophysics . the main problem regarding the nature of grbs has usually been the lack of knowledge of their distance scale . about 300 grbs are detected annually by batse in the full sky , but only a few of them can be localized accurately to less than half a degree . for many years , follow - up observations by other satellites and ground - based telescopes were conducted , but no counterparts at other wavelengths were found . the breakthrough took place in 1997 , thanks to the observation by _ bepposax _ and _ rossixte _ of the fading x - ray emission that follows the more energetic gamma - ray photons once the grb event has ended . this emission ( the afterglow ) extends at longer wavelengths , and the good accuracy in the position determination by _ bepposax _ has led to the discovery of the first optical counterparts -for grb 970228 , grb 970508 , and grb 971214- , greatly improving our understanding of these puzzling sources . now it is widely accepted that most bursts originate at cosmological distances but the final solution of the grb problem is still far away . + key  words : multiwavelength observations ; gamma - ray bursts .


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: it has been shown experimentally that photoelectron spectra obtained from strong - laser matter interaction contain a wealth of information about the driving laser field and the target @xcite . decoding this information and discriminating the various spectral features with respect to laser parameters , target geometry , many - electron effects etc . is often difficult . on one hand , experiments and _ ab initio _ solutions of the time - dependent schrdinger equation ( tdse ) yield the `` exact '' spectra but little insight . on the other hand , the widely used workhorse of strong - field physics , the strong field approximation ( sfa ) @xcite , allows for an intuitive interpretation when formulated in terms of interfering quantum orbits via the saddle - point approximation ( see , e.g. , @xcite ) . however , the sfa is rather inaccurate for the case of long - range coulomb potentials due to the neglect of the interaction between the photoelectron and the parent ion after the time of ionization . in fact , various spectral features originating from this coulomb interaction have been discovered recently , among them the so - called `` low - energy structure '' @xcite , asymmetries in elliptically polarized laser fields @xcite , and holographic side - lobes @xcite .    within the trajectory - based sfa , the matrix element for ionization is written as a coherent sum over ionization times for which the corresponding quantum orbits all lead to a given , final momentum @xmath0 at the detector ( see , e.g. , @xcite ) . if there are more than one such orbits of similar weight , interference comes into play . interfering orbits originating from subsequent laser cycles give rise to peaks in energy spectra that are separated by @xmath1 with @xmath2 being the laser frequency , the so - called `` above - threshold ionization '' ( ati ) peaks . we call this `` inter - cycle interference '' , in contrast to `` intra - cycle interference '' @xcite where the interfering orbits originate from within one laser cycle , leading to , e.g. , the holographic interference patterns observed in @xcite . intra- and inter - cycle interference in time is in analogy to a double slit and a grating in position space , respectively . once the sfa is formulated in terms of trajectories , one may coulomb - correct them @xcite . this is straight - forwardly done in the classically allowed region , i.e. , from the so - called `` tunnel exit '' to the detector . in the current paper we will extend the coulomb - correction of photoelectron spectra to the sub - barrier region where time and orbit are complex . the imaginary time method has since many years been employed for the derivation of coulomb - corrected tunneling ionization _ rates _ ( see @xcite for a review ) . in fact , the difference in the probability for tunneling through a triangular barrier of a short - range potential and through the corresponding coulomb barrier can be many orders of magnitude @xcite . however , to the best of our knowledge sub - barrier coulomb effects on the interference pattern of photoelectron _ spectra _ have not been investigated . the outline of our paper is as follows . in section [ sec : tcsfa ] we review the trajectory - based coulomb - corrected strong - field approximation ( tcsfa ) in real time and space before we introduce , in subsection [ sec : subcc ] , the sub - barrier coulomb correction ( sub - cc ) . in section [ sec : intracyc ] results on the intra - cycle interference are presented . in particular , the fringe shift in polarization direction due to the sub - cc is derived and discussed . there is no sub - cc exactly in perpendicular direction , as discussed in subsection [ sub : perpendicular_direction_phase_analysis ] . the nodal structure along low - order ati rings is investigated in section [ sec : ationring ] . we summarize in section [ sec : summ ] . atomic units @xmath3 are used , unless noted otherwise . a quantum orbit in plain sfa is a semi - classical trajectory of a photoelectron starting at a complex time @xmath4 , subsequently governed by the external light field , and finally detected with the asymptotic momentum @xmath0 . the transition amplitude between the bound state to the continuum state of asymptotic momentum @xmath0 is a coherent superposition of contributions from all relevant quantum orbits ( indexed by @xmath5 ) that lead to the asymptotic momentum @xmath0 , @xmath6 @xmath7 is a prefactor whose form is not of interest in this work ( see , e.g. , @xcite for more details ) . @xmath8dt\ ] ] is the action in which @xmath9 with the vector potential of the laser field @xmath10 in dipole approximation , and @xmath11 is the ionization potential . all the trajectories @xmath5 leading to an asymptotic electron momentum @xmath0 fulfill the equation of motion @xmath12 and start at the complex saddle - point ( i.e. , ionization ) time @xmath13 , which is determined by the saddle point equation ( spe ) @xmath14^{2}=-i_{p}.\label{eq : spe}\ ] ] with @xmath15 and @xmath0 real - valued , @xmath13 is complex . from ( [ eomsfa ] ) and ( [ vpa ] ) follows that the canonical momentum is conserved , i.e. , @xmath16 const . the amplitude ( [ spsum ] ) clearly demonstrates the multi - slit in time - nature of quantum mechanical ionization dynamics @xcite and lays the cornerstone for our further analysis . the key idea of the tcsfa is to insert the coulomb force due to the ion of nuclear charge @xmath17 into the electronic equation of motion , @xmath18 and the corresponding coulomb potential into the action , @xmath19dt,\label{eq : action_tcsfa}\ ] ] so that the amplitude ( [ spsum ] ) becomes @xmath20 here , we write @xmath21 to indicate that @xmath22const now . instead @xmath23 . the saddle - point times are also modified ( hence @xmath24 ) , as they satisfy @xmath25^{2}=-i_{p } \label{spedash}.\ ] ]     we integrate parallel to the imaginary time axis down to @xmath26 on the real axis , and then along the real axis to infinity . the corresponding situations in position space around the tunneling barrier are illustrated : at @xmath27 the electron is inside the atom in front of the barrier , at @xmath26 it leaves the barrier at the tunnel exit and subsequently moves in the classically allowed region . ( b ) real parts of the saddle point times for a monochromatic laser with the electric field @xmath28 and vector potential @xmath29 , @xmath30 . the superscripts of @xmath31 indicate the signs @xmath32 , @xmath33 of the asymptotic momentum @xmath34 ; the numbers @xmath35 and @xmath36 indicate whether the corresponding quantum trajectory is `` long '' or `` short , '' respectively . the negative vector potential ( red ) and the electric field ( blue ) in arbitrary units are indicated . panels ( c)(e ) illustrate the asymptotes of real trajectories in perpendicular direction @xmath37 ( c ) , for an angle @xmath38 ( d ) , and parallel to the polarization direction , @xmath39 ( e ) . the bold solid lines @xmath40 and @xmath41 represent trajectories in plain sfa where initial and asymptotic momentum @xmath0 are equal . thin solid lines represent the corresponding tcsfa trajectories , which start off like the dashed trajectories with a modified initial momentum but then get coulomb - distorted such that the asymptotic momentum is @xmath0 . [ fig : demo],scaledwidth=45.0% ]    as illustrated in fig . [ fig : demo](a ) , the integration in complex time in @xmath42 can be splitted into two parts , @xmath43 the sub - barrier action @xmath44 , where @xmath45 , may be calculated along the path parallel to the imaginary axis down to the real axis . provided the integrand is analytic , the result for the total action @xmath46 is independent of the integration path . the action @xmath47 in the classically allowed region is evaluated along the real axis . the well - established imaginary time method ( itm ) @xcite has been used to study quantum tunneling through ( time - dependent ) barriers in various contexts ( e.g. , in @xcite ) . at @xmath48 , the electron is in the bound state with dynamical variables in general being complex . at @xmath49 , the electron emerges at the outer turning point @xmath50 ( tunnel exit ) . the boundary conditions for the quantum orbits may be chosen such that from that time @xmath26 on all variables are purely real @xcite . let us first summarize the coulomb - correction to this real part of the propagation , as it has been employed in @xcite .      for @xmath51 , the real trajectories with the coulomb interaction taken into account are obtained by solving the equation of motion ( [ eomfull ] ) with the initial conditions @xmath52 ( tunnel exit ) where @xmath53 is the electron excursion , and @xmath54 the corresponding saddle - point times @xmath55 are determined via eq . ( [ spedash ] ) . the equations of motion are solved numerically , along with @xmath56 , up to the end of the laser pulse . the asymptotic momentum @xmath57 can then be determined via kepler s formulas . the values of @xmath57 , @xmath58 , and other entities of interest are saved in a data base . in the previous works @xcite the uncorrected @xmath59 ( i.e. , without coulomb potential ) of the plain sfa with @xmath60 was evaluated analytically . finally , photoelectron spectra @xmath61 can be calculated by adding coherently all the @xmath62 of the quantum orbits from the data base that have their asymptotic momentum @xmath57 in a bin around the final momentum of interest @xmath0 . the size of the bin determines the resolution of the photoelectron momentum spectrum . in the current work we used @xmath63 bins to sample the final momentum interval @xmath64\times[0,1]$ ] in the spectra of fig . [ fig : spec_2d_sin2_2um_1e14_nc4](b)(d ) below . thanks to the azimuthal symmetry with respect to the laser polarization direction @xmath65 it is sufficient to consider only , e.g. , @xmath66 . we covered randomly the initial momentum interval @xmath64\times[0,2]$ ] by 32 million points . for each initial momentum several saddle - points contribute . in the case of the 3-cycle pulse considered below we thus dealt with @xmath67 million trajectories . because of focusing along the polarization axis the randomly sampled initial momenta in lateral direction @xmath68 should cover a larger interval than the final @xmath69 of interest . on the other hand , it rarely occurs that @xmath70 so that the initial momentum interval in polarization direction may be chosen the same as the final one . if , for an initial momentum @xmath71 the final momentum @xmath69 turns out to be @xmath72 we simply reverse the signs of both @xmath68 and @xmath69 . the finer the binning and the higher the dynamic range of interest the more trajectories are required to ensure good enough statistics . it is also advisable to shoot randomly in the initial momentum plane , not in a regular , grid - like pattern , in order to avoid artificial structures in the final momentum spectra . the numerical effort to obtain tcsfa spectra like the one in fig . [ fig : spec_2d_sin2_2um_1e14_nc4](c ) by solving the equations of motion for 100 million trajectories is significant but trivially parallelizable @xcite . one may attempt to solve the equation of motion ( [ eomfull ] ) in complex space and time under the barrier as well , for instance , by propagating from the tunnel exit `` backwards '' in imaginary time from @xmath26 to @xmath27 . however , we found that the photoelectron spectra obtained in this manner are completely spoiled by artifacts such as caustics , destroying the interference pattern observed in the _ ab initio _ tdse results . if one wanted to propagate in forward direction from @xmath27 to @xmath26 it is not clear what the initial conditions should be . choosing , e.g. , @xmath73 one would have to iterate @xmath74 and the imaginary @xmath75 such that at the tunnel exit position and velocity become real , and @xmath76 because @xmath27 is determined by @xmath21 . in the derivation of tunneling ionization _ rates _ @xcite such problems are circumvented by matching the action @xmath77 , evaluated with the zeroth order ( i.e. , plain sfa ) quantum orbit with the highest statistical weight , to the exponent of the unperturbed , initial state wavefunction . this is justified , as for @xmath78 the coulomb potential certainly dominates the external laser field . moreover , the quantum orbits are almost linear @xcite and spatially confined in the polarization direction , pointing from the origin to the tunnel exit . we pursue a similar strategy in this work in so far that we consider the plain - sfa quantum orbits for the evaluation of @xmath77 only . of course , for the calculation of photoelectron _ spectra _ we can not restrict ourselves to just the dominating quantum orbits . the plain - sfa quantum orbits fulfill @xmath79 where @xmath80 is the randomly sampled momentum from the real - part propagation described in the previous subsection . the sub - barrier trajectory is @xmath81 which ensures that @xmath82 is purely real and @xmath83 is purely imaginary . we write the sub - barrier action as @xmath84 where @xmath85dt\label{eq : w_sub_0}\ ] ] is the plain - sfa action , and @xmath86 is the sub - cc action . here , we write on purpose @xmath87 instead of the non - analytic @xmath88 .    in the preceding studies @xcite only @xmath89 was taken into account whereas in this work , the action takes the full form @xmath90 including the sub - cc action @xmath91 . the action may diverge if the complex trajectory @xmath92 . however , thanks to the negligible probability of sampling momenta @xmath21 that lead to such singularities they turned out to be harmless in the actual numerical implementation of our sub - cc . in order to investigate the influence of the sub - cc on the interference pattern in photoelectron momentum spectra we first consider h(1s ) ( @xmath93 ) in a linearly polarized laser pulse . the vector potential of the laser pulse of carrier frequency @xmath2 is given ( in dipole approximation ) by @xmath94 for @xmath95 and zero otherwise . @xmath96 is the peak electric field strength , @xmath97 is the number of laser cycles . figure  [ fig : spec_2d_sin2_2um_1e14_nc4](a ) shows the photoelectron momentum spectrum , calculated from the exact numerical solution of the time - dependent schrdinger equation ( tdse ) using the qprop code @xcite . the laser parameters are @xmath98 , @xmath99 , @xmath100 ( corresponding to @xmath101 @xmath102 m and an intensity of @xmath103w/@xmath104 ) . here we focus on the almost vertically aligned interference fringes @xcite . these fringes stem from the `` intra - cycle interference '' @xcite of two trajectories `` born '' at times @xmath105 and @xmath106 at the respective tunnel exits within one laser cycle . the electric field has opposite signs at @xmath105 and @xmath106 , meaning that the tunnel exits point in opposite directions . in the context of the sfa , the two trajectories are typically denoted by `` short '' and `` long '' according to whether the signs of the projections of @xmath50 and @xmath80 onto the polarization direction are equal or opposite , respectively . in the former case the photoelectron drifts directly from the tunnel exit to the detector whereas in the latter case it starts off into the `` wrong '' direction before it is reversed by the laser field . the vertical fringes are also present in the sfa spectrum fig . [ fig : spec_2d_sin2_2um_1e14_nc4](b ) . however , the positions are shifted relative to the tdse result , as indicated by the vertical , dashed lines . in fact , when there is constructive interference in the tdse result , there is destructive interference in the sfa spectrum , and _ vice versa_. this indicates that there is an erroneous phase - shift between the long trajectory and the short trajectory of @xmath107 in the plain sfa as compared to the exact tdse result . it is expected that incorporating the coulomb - attraction between the outgoing electron and the parent ion should cure this disagreement . let us now `` switch - on '' the coulomb - attraction for the propagation in the classically allowed region ( i.e. , from the tunnel exit to the detector ) , as described in section [ sec : repart ] . the result is shown in fig . [ fig : spec_2d_sin2_2um_1e14_nc4](c ) . as discovered previously @xcite , side - lobes are reproduced in that way . however , the vertical interference fringes still have the same position as in the plain - sfa result , meaning that the erroneous phase - shift between the long and the short trajectory of @xmath107 is _ not _ cured by the coulomb - correction out of the coulomb - barrier . only if the sub - cc , as described in section [ sec : subcc ] , is incorporated as well , do the fringe positions coincide with the _ ab initio _ tdse result [ see fig .  [ fig : spec_2d_sin2_2um_1e14_nc4](d ) ] . hence , the sub - cc not only affects the statistical weight of a quantum orbit via the imaginary part of @xmath91 but also the interference pattern in photoelectron spectra via the real part of @xmath91 . this finding ( and its analytical derivation in the following ) is the main result of this paper . -@xmath108 plane ( @xmath109 along the polarization direction ) calculated using ( a ) tdse , ( b ) plain sfa , ( c ) tcsfa without sub - cc , and ( d ) tcsfa with sub - cc . the calculation is for h(1s ) in a 4-cycle , @xmath110-envelope , near - infrared laser pulse of wavelength @xmath111 m and intensity @xmath103 w/@xmath104 . four orders of magnitude in probability are shown . [ fig : spec_2d_sin2_2um_1e14_nc4 ] ]    let us suppose that only two quantum orbits are dominating the transition amplitude for a certain asymptotic momentum @xmath0 . given the two individual transition amplitudes @xmath112 and @xmath113 , the probability to find a photoelectron with asymptotic momentum @xmath0 is @xmath114^{*}\left[m_{\boldsymbol{p}}^{(1)}+m_{\boldsymbol{p}}^{(2)}\right]$ ] with @xmath115 . the action @xmath116 is of the form ( [ eq : action_w ] ) . interference with pronounced contrast requires the weights ( i.e. , absolute values ) of @xmath112 and @xmath113 being comparable , @xmath117 . the assumption that this is the case then yields an ionization probability @xmath118 . hence , the interference pattern is determined by the phase difference @xmath119 $ ] . for @xmath120 with @xmath121 even we have interference maxima , for @xmath121 odd interference minima , provided the prefactors @xmath7 in ( [ spsum ] ) do not introduce extra phase - differences @xcite . the separation of the action ( [ eq : action_w ] ) suggests that we split the phase @xmath122 accordingly , @xmath123 from the simulation results in fig . [ fig : spec_2d_sin2_2um_1e14_nc4 ] we know already that it is @xmath124 that is responsible for the intra - cycle fringe shifts towards the correct positions , as predicted by the _ ab initio _ tdse result . in the following we will derive an analytical expression for @xmath124 in the long - wavelength limit .      in the long - wavelength limit ( @xmath125 ) it is possible to evaluate the sub - cc phase difference @xmath124 in the polarization direction analytically . consider a monochromatic laser field @xmath126 where @xmath127 , and we assume @xmath128 . we will prove in this subsection that in the long - wavelength limit we obtain a momentum - independent sub - cc phase shift of @xmath129 in particular , for the case of fig . [ fig : spec_2d_sin2_2um_1e14_nc4 ] where @xmath130 we have @xmath131 , turning constructive interference maxima of the plain sfa into destructive interference minima in the tcsfa with sub - cc and _ vice versa_. as we know already from the results in fig . [ fig : spec_2d_sin2_2um_1e14_nc4 ] that the coulomb correction of the propagation in real space is _ not _ responsible for the fringe shift , we ignore its effect on @xmath40 and @xmath41 ( and thus also on @xmath132 and @xmath133 ) in this section and assume @xcite @xmath134 for a monochromatic laser field ( [ eq : mono_pulse ] ) , @xmath135 , and @xmath136 the saddle point equation ( [ eq : spe ] ) yields @xmath137 equation  ( [ eq : spe_reduced1 ] ) shows that the possible @xmath31 can be determined as the intersections between @xmath138 and the negative vector potential @xmath139 , as indicated in fig . [ fig : demo](b ) . the figure shows that a pair of saddle points per cycle , e.g. , @xmath140 and @xmath141 , exists for a specified @xmath34 ( dashed horizontal lines ) . the so - called intra - cycle interference is thus due to two saddle - points for which @xmath142 are fulfilled . in the polarization direction , @xmath143 , the sub - cc action is @xmath144 with @xmath145 along the integration path @xmath146 $ ] , as shown in fig . [ fig : demo](a ) . the superscript @xmath147 is suppressed here . for the sub - barrier trajectory in polarization direction calculated from ( [ eq : complex_traj ] ) @xmath148+p_{z}\left[t-\mathrm{re}\,t_{s}\right]\label{eq : complex_traj_mono_z}\ ] ] results . separating the real and imaginary parts of @xmath149 , @xmath150 one finds @xmath151 in the long wavelength limit @xmath125 , the functions @xmath152 and @xmath153 in @xmath154 and @xmath155 can be expanded about @xmath156 , leading to @xmath157 with @xmath158 let us consider emission times @xmath26 for positive @xmath34 . in this case @xmath159 and thus @xmath160 . as @xmath161 @xcite as well , @xmath162 . hence the integrand in @xmath91 in ( [ eq : wccsub ] ) may be simplified [ by choosing the positive root ] to @xmath163^{-1/2}=z/\tilde{z}(\tau)$ ] for @xmath164 and @xmath163^{-1/2}=-z/\tilde{z}(\tau)$ ] for @xmath165 so that it is sufficient to consider the simplified action @xmath166 the real part of this action , @xmath167 is relevant for interference patterns in the spectra . inserting ( [ eq : a_and_b_in_limit ] ) leads with ( [ eq : spe_reduced2 ] ) for @xmath168 and @xmath169 , @xmath170 to    @xmath171\right|_{t_i}^0\\ \stackrel{\delta\to 0}{\simeq }   \mp \frac{z}{\sqrt{2i_p } } \ , \frac{\pi}{2}. \label{eq : rewsubccpm}\ ] ]    it follows from ( [ eq : a_and_b_in_limit ] ) and ( [ eq : sfa_condition ] ) that , for a given @xmath34 , the two saddle - point times @xmath132 and @xmath133 are such that @xmath172 and @xmath173 . the opposite signs of @xmath154 have been considered already in our calculation of @xmath174 . as claimed in ( [ eq : claim ] ) .    for a hydrogen - like ion in a state with principal quantum number @xmath176 we have @xmath177 and thus @xmath178 . in particular , for ground states @xmath179 we have independently of @xmath17 that @xmath180 . this is the required phase shift for bringing the intra - cycle fringe pattern in agreement with the _ ab initio _ tdse result [ see fig .  [ fig : spec_2d_sin2_2um_1e14_nc4](d ) ] . we compared tdse and tcsfa results for other parameters @xmath17 , @xmath11 , and @xmath181 and find good agreement if the sub - cc is included , provided the minimum number of photons required for ionization @xmath182 and the keldysh parameter @xmath183 , where @xmath184 is the ponderomotive energy . moreover , the smaller @xmath17 the better is the approximation to use the zeroth - order , plain - sfa trajectory for the sub - cc . excited states @xmath185 may introduce extra phase - differences via the prefactors @xmath7 in ( [ spsum ] ) @xcite , not discussed in this work .    the reader may have noticed that the precursors of ati rings , visible in the exact tdse result of fig . [ fig : spec_2d_sin2_2um_1e14_nc4](a ) and in the sfa result ( b ) as circular structures around the origin , are actually less clear in the tcsfa ( c ) , let alone in the sub - cc tcsfa ( d ) . the reasons for this deterioration are the overemphasized classical boundaries and caustics that spoil for few - cycle pulses the low - energy part of the spectra by rendering the few trajectories that contribute to the ati - ring structure relatively less important . the situation improves for more - cycle pulses , as in fig . [ fig : spec_2d_const_800nm_1e14 ] below , where more trajectories contribute to inter - cycle features such as ati rings . we are also interested in effects from the sub - cc in the perpendicular direction @xmath186 . in the sfa , the short and long trajectories are asymptotically parallel , as shown in figs . [ fig : demo](c  e ) . in the perpendicular direction [ fig . [ fig : demo](c ) ] , the parallel asymptotes are vertical when two saddle points merge as @xmath187 in ( [ eq : spe_reduced1 ] ) . one easily verifies that in fig . [ fig : demo](b ) @xmath188 merges with @xmath189 and @xmath190 merges with @xmath191 . hence , in plain sfa the momenta and the saddle point times of long and short trajectories become degenerate in perpendicular direction [ @xmath40 and @xmath41 in fig .  [ fig : demo](c ) ] . however , with the long - range coulomb interaction included , this degeneracy is lifted . in order to have a common asymptotic momentum @xmath192 the initial momenta @xmath193 must be different , as indicated in fig . [ fig : demo](c ) .    in order to cancel the longitudinal shift induced by the long - range coulomb distortion , the two trajectories should satisfy @xmath194 the corresponding saddle points are easily found for these criteria from fig . [ fig : demo](b ) . taking , e.g. , the saddle point @xmath195 as the reference , the other saddle point is @xmath196 . therefore , the initial conditions for these trajectories are @xmath197    in the following , the phase contributions @xmath198 , @xmath199 and @xmath200 are deduced , and the interference pattern in the perpendicular direction is inferred . the subscript `` @xmath201 '' indicates that the phase difference between the trajectories has been generalized to allow for the real - space trajectory - distortion by the coulomb force . the interference pattern in perpendicular direction is determined by the symmetry properties of the interfering trajectories , while it is independent of the detailed , coulomb - modified dynamics because both trajectories of a pair collect equal phase distortions . this suggests that the sfa , tcsfa with or without sub - cc result in the same interference structure in the perpendicular direction , as will be discussed in the following .      we show first that the sub - cc phase difference @xmath198 does not contribute in ( [ eq : phase_angle ] ) , @xmath202=0.\label{eq : phi_c_subcc_perpendicular}\ ] ] the real part of the sub - cc action reads @xmath203 for a monochromatic laser pulse ( [ eq : mono_pulse ] ) , the zeroth - order sub - barrier trajectory ( @xmath204 as a function of @xmath205 can be cast into [ see ( [ eq : aandb ] ) , ( [ eq : a_and_b ] ) ] @xmath206 and @xmath207 substituting @xmath208 and @xmath209 into @xmath149 ( [ eq : z_tau_sub_mono ] ) , we find with the help of the conditions ( [ eq : perpendicular_conditions ] ) @xmath210 . for @xmath211 ( [ eq : x_tau_sub_mono ] ) clearly follows @xmath212 . as a consequence , the integrand in ( [ eq : re_w_subcc_components ] ) is identical for both trajectories . as the integration starts for both trajectories at @xmath213 we obtain @xmath214 , i.e. , the sub - cc phase difference vanishes , and ( [ eq : phi_c_subcc_perpendicular ] ) is proved .    in passing , we would like to mention that caution has to be exercised regarding which saddle - points are chosen for @xmath215 if zeroth - order , _ plain _ sfa were used , for which the degeneracy of the saddle points is _ not _ lifted . if the `` usual '' long and short trajectories @xmath195 and @xmath216 with @xmath217 and @xmath218 shown in fig . [ fig : demo](b ) are chosen , @xmath219 can be expressed by variables of @xmath195 as @xmath220 compared to ( [ eq : z_tau_sub_mono ] ) , this expression differs by a sign only in the first term . since @xmath221 , the imaginary parts of @xmath222 in the denominator of ( [ eq : re_w_subcc_components ] ) for @xmath195 and @xmath216 are opposite in signs . it seems a paradox occurs : in the exact perpendicular direction , integrands are identical since the imaginary part in ( [ eq : z2plus2 ] ) is zero , and thus no phase difference exists . but a tiny longitudinal momentum @xmath223 would lead to the abrupt appearance of a finite phase difference , which is unphysical . the problem does not arise in the above employed sub - cc , thanks to the coulomb distortion of the trajectories in real - space , which enforces @xmath224 , thus determining the relevant saddle - point times to be @xmath195 and @xmath196 , _ not _ @xmath195 and @xmath216 .      after having identified the relevant saddle - points @xmath195 and @xmath196 we now show that @xmath225=0.\label{eq : phi_c_sub0_perpendicular}\ ] ] substituting a monochromatic laser field ( [ eq : mono_pulse ] ) into ( [ eq : w_sub_0 ] ) yields @xmath226 where @xmath227 . the real part reads @xmath228 inserting @xmath209 and @xmath208 , which satisfy ( [ eq : perpendicular_conditions ] ) , we find @xmath229 confirming ( [ eq : phi_c_sub0_perpendicular ] ) .      the interference pattern in the perpendicular direction is entirely determined by the real - time trajectories and the corresponding phase difference @xmath230 given real trajectories @xmath231 , @xmath232 , @xmath233 , @xmath234 for @xmath235 $ ] , the action due to the propagation in real time is @xmath236\\   & -\frac{z}{\sqrt{z^2(t)+x^2(t)}}+i_{p}\biggr]dt . \end{split}\ ] ] the symmetry properties for two interfering trajectories in the perpendicular direction are @xmath237 moreover , @xmath238 . casting @xmath239 in terms of the variables of @xmath195 one finds @xmath240\\   & \quad   -\frac{z}{\sqrt{{z_{\perp}^{(1+)}}(t)^2+{x_{\perp}^{(1+)}}(t)^2}}+i_{p}\biggr\}\ , dt \\ & = - \int_{t_{r}^{(1+)}}^{t_{p}+\pi/\omega } i(t)\,d t. \end{split}\end{aligned}\ ] ] the phase difference @xmath241 is thus given by @xmath242 independent of the specific saddle - point pair under consideration . as @xmath243 the asymptotic position in polarization direction remains finite , @xmath244 , but @xmath245 , leading to a vanishing phase difference due to the coulomb potential . in the kinetic energy term of the integrand , @xmath246 and @xmath247 , so that @xmath248 using @xmath249 from the energy conservation for ati in the perpendicular direction with @xmath181 the number of absorbed photons and @xmath250 we find @xmath251 this result agrees with the findings in ref . @xcite that every other ati ring is destructively `` interfered away '' in perpendicular direction . because of the pre - factor @xmath7 in ( [ spsum ] ) the interference is destructive for odd @xmath181 in the case of even - parity initial states and _ vice versa_.    the interference concerning the ati rings in perpendicular direction predicted by ( [ eq : phi_c_re_perpendicular ] ) is confirmed by the results for ionization of h(1s ) in a 9-cycle trapezoidal - envelope laser pulse of 800 nm wavelength shown in fig .  [ fig : spec_2d_const_800nm_1e14 ] . the other laser parameters are given in the figure caption . ati peaks appear as rings of radius @xmath252 in the photoelectron momentum spectra . it is known from experiments that the ionization probability along the ati rings is not isotropic but exhibits nodal structures . these structures , also visible in fig . [ fig : spec_2d_const_800nm_1e14 ] , have been analyzed recently within the sfa @xcite . the radius of the ati rings is well reproduced already by the plain sfa , which is verified by comparing fig .  [ fig : spec_2d_const_800nm_1e14](b ) with the _ ab initio _ tdse result fig . [ fig : spec_2d_const_800nm_1e14](a ) . however , the positions of the nodes along low - order ati rings do not agree well . not even the _ number _ of minima and maxima agrees . take , e.g. , the second ati ring ( @xmath253 , absorption of 14 photons ) . the tdse result in ( a ) shows 9 maxima along that ring , the plain sfa result in ( b ) only 7 . the real - time coulomb - corrected sfa in ( c ) does not affect the number of maxima . instead , the sub - cc sfa in ( d ) does give the correct number . along the next ati - ring ( absorption of 15 photons ) tdse and plain sfa agree already ( 10 maxima ) , and tcsfa with or without sub - cc does not change the nodal structure . along the first ati ring the tdse shows 8 maxima , which is not properly reproduced , even with the sub - cc tcsfa , probably because our sub - barrier coulomb correction involving only the zeroth - order plain - sfa trajectories is insufficient . but for a 9-cycle , trapezoidal - envelope laser ( 2 , 5 , and 2 cycles for up - ramping , constant - envelope , and down - ramping , respectively ) . the wavelength is @xmath254 nm , the intensity is @xmath103 w/@xmath104 . in the spectrum calculated using the tcsfa with sub - cc ( d ) only trajectories of type t1 and t2 in the nomenclature of @xcite have been considered @xcite . [ fig : spec_2d_const_800nm_1e14 ] ]    figure  [ fig : on_ring_nodal ] shows the momentum distribution along the second ati ring for the _ ab initio _ tdse , the sfa , and the tcsfa with or without sub - cc . tcsfa with sub - cc agrees well with the _ ab initio _ tdse result . the results from the tcsfa without sub - cc and the plain sfa are similar . they agree well with the tdse around @xmath37 but deviate for smaller angles . as discussed above , there is one maximum less in the plain sfa and the tcsfa without sub - cc , as compared to the tdse and the tcsfa with sub - cc . hence we conclude that the nodal structure of low - order ati rings is affected by the sub - cc . higher - order ati rings may be effected by the side - lobes ( clearly visible in fig . [ fig : spec_2d_const_800nm_1e14](a ) and discussed in @xcite ) .     in fig . [ fig : spec_2d_const_800nm_1e14 ] according to the _ ab initio _ tdse , the ( plain ) sfa , the tcsfa ( without sub - cc ) , and the tcsfa with sub - cc result . the polar angle @xmath255 is in the range @xmath256 $ ] , @xmath39 corresponding to @xmath168 , and @xmath37 corresponding to @xmath215 . the plain sfa result is as noisy as the tcsfa because it was calculated using the same statistical sampling approach as for the tcsfa ( which is explained briefly in sec . [ sec : repart ] and in more detail in @xcite ) . [ fig : on_ring_nodal],scaledwidth=35.0% ] we extended the previously introduced trajectory - based coulomb - corrected strong field approximation towards the inclusion of sub - barrier coulomb effects . to that end we used the plain strong field approximation quantum orbits and inserted them into the coulomb - corrected sub - barrier action . the ( additional ) phase - difference between the so - called `` long '' and `` short '' trajectories introduced in that way has a pronounced qualitative effect on the photoelectron spectra , namely a phase shift of @xmath257 along the polarization axis for linear polarization in the long - wavelength limit . in the case of ionization from the ground state of atomic hydrogen or hydrogen - like ions , the phase shift is @xmath107 , turning constructive intra - cycle interference maxima into destructive interference minima and _ vice versa_. we benchmarked the predictions of the trajectory - based coulomb - corrected strong field approximation including sub - barrier effects by comparisons with _ ab initio _ solutions of the time - dependent schrdinger equation . we also showed that the sub - barrier coulomb correction affects the nodal structure along low - order above - threshold ionization rings . our semi - classical method is applicable to sufficiently long wavelengths ( @xmath182 ) . however , the nuclear charge @xmath17 ( and thus @xmath11 ) has to be small enough to render the zeroth - order orbit a good approximation for the sub - barrier coulomb correction . moreover , the trajectory - based coulomb - corrected strong field approximation displays too many unphysical caustic structures in the spectra if the keldysh parameter @xmath258 is not small enough ( i.e. , around @xmath35 or smaller ) . it thus turns out that our method is applicable in a regime where the direct _ ab initio _ solution of the time - dependent schrdinger equation in position space is stretched to its limits : long wavelengths and high ponderomotive energy . our method is , in principle , also applicable to more complex systems for which strong field approximations already have been or can be developed . the most straightforward extension is towards two coulomb - centers ( e.g. , h@xmath259 ) . also non - sequential ionization in effective two - electron systems is accessible . however , due to the higher - dimensional momentum space to be sampled in such systems it will be numerically very demanding to obtain differential photoelectron spectra with good enough statistics . the influence of many - electron polarization effects could be investigated by incorporating a polarization potential into the equations of motion for the quantum trajectories @xcite . we thank s.v . popruzhenko for fruitful discussions . the work was supported by the deutsche forschungsgemeinschaft ( sfb 652 ) . yan acknowledges support from the international max planck research school for quantum dynamics ( imprs - qd ) in heidelberg . references r.  gopal , k.  simeonidis , r.  moshammer , th . ergler , m.  drr , m.  kurka , k .- u . khnel , s.  tschuch , c .- d . schrter , d.  bauer , j.  ullrich , a.  rudenko , o.  herrwerth , th . uphues , m.  schultze , e.  goulielmakis , m.  uiberacker , m.  lezius , and m.f . kling , phys . lett .  * 103 * , 053001 ( 2009 ) . y.  huismans , a.  rouze , a.  gijsbertsen , j.  h.  jungmann , a.  s.  smolkowska , p.  s.  w.  m.  logman , f.  lpine , c.  cauchy , s.  zamith , t.  marchenko , j.  m.  bakker , g.  berden , b.  redlich , a.  f.  g.  van der meer , h.  g.  muller , w.  vermin , k.  j.  schafer , m.spanner , m.  y.  ivanov , o.  smirnova , d.  bauer , s.  v.  popruzhenko , and m.  j.  j.  vrakking , science * 331 * , 61 ( 2011 ) . y.  huismans , a.  gijsbertsen , a.  s.  smolkowska , j.  h.  jungmann , a.  rouze , p.  s.  w.  m.  logman , f.  lpine , c.  cauchy , s.  zamith , t.  marchenko , j.  m.  bakker , g.  berden , b.  redlich , a.  f.  g.  van der meer , m.  yu .  ivanov , t .- m .  yan , d.  bauer , o.  smirnova , and m.  j.  j.  vrakking , phys .  rev .  lett .  * 109 * , 013002 ( 2012 ) . exaggerated caustic structures appear that spoil the interference structure if t3 and t4 trajectories ( in the nomenclature of @xcite ) are included . the situation improves with increasing @xmath261 and decreasing keldysh parameter @xmath258 .
Abstract:

Answer: we use a quantum trajectory - based semi - classical method to account for coulomb interaction between the photoelectron and the parent ion in the classically forbidden , sub - barrier region during strong - field tunneling ionization processes . we show that  besides the well - known modification of the tunneling ionization probability  there is also an influence on the interference pattern in the photoelectron spectra . in the long - wavelength limit , the shift of the intra - cycle interference fringes caused by sub - barrier coulomb effects in the laser polarization direction can be derived analytically . we compare our results with _ ab initio _ solutions of the time - dependent schrdinger equation and find good agreement in the long - wavelength regime , whereas the standard strong field approximation fails . we show that the nodal structure along low - order above - threshold ionization rings is also affected by sub - barrier coulomb effects .