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 ultracompact binary stars v407 vul ( rx j1914.4 + 2456 , * ? ? ? * ; * ? ? ? * ) and hm cnc ( rx j0806.3 + 1527 , * ? ? ? * ; * ? ? ? * ) have been the subject of extensive study since their discovery . believed to be members of the interacting binary white dwarf ( am cvn ) class , each system displays only a single photometric period , presumed in both cases to be the orbital period . for hm cnc , @xmath0 min @xcite , and for v407 vul , @xmath1 @xcite . in a previous publication ( * ? ? ? * hereafter paper i ) we simulated the hydrodynamics of the accretion stream using the method of smoothed particle hydrodynamics ( sph ) . we found that the impact spot was small enough and thus the accretion stream ram pressure high enough for the specific kinetic energy of the flow to be advected and thermalized below the photosphere of the primary , such that it can be emitted `` downstream '' of the impact point in the equatorial flow that should develop . in paper i we discussed the physical nature of that accretion - driven flow and suggested the general features of any plausible photospheric temperature distribution that should result . the observational constraints for the systems hm cnc and v407 vul were discussed in detail in paper i ( also see * ? ? ? * ) . briefly , the systems are characterized by `` on / off '' x - ray light curves that trail the optical light curves by roughly 0.2 in phase ( * ? ? ? * and references therein ) , no hard x - rays are detected , the observed periods are decreasing with time at a rate consistent with gravitational radiation angular momentum loss @xcite , and flickering is minimal or absent @xcite . there are several models that have been proposed to explain the observations of these two systems : intermediate polar , polar , unipolar inductor , and direct impact . first , in the intermediate polar model originally suggested for v407 vul , the observed periods reflect a white dwarf spin period @xcite . however , the lack of observed flickering , emission features , or any other observed ( e.g. , orbital ) period , as well as the `` on / off '' character of the x - ray light curve , argue against this model . citing these issues @xcite proposed that v407 vul could be a helium analog of a polar ( am her ) system ( see also * ) , but the lack of observed he emission lines , polarization , and hard x - ray photons make it unlikely this model is correct for these systems ( e.g. , * ? ? ? the unipolar inductor ( ui ) model @xcite is the only model without a roche - lobe filling secondary . similar to the model proposed by @xcite for the io - jupiter system in which the conducting moon io generates an electromotive force ( emf ) as a result of the relative motion of the magnetic field of jupiter . the weak plasma environment allows currents to flow , driven by this emf , and results in observable bright spots where these currents are dissipated in the atmosphere of jupiter . the ui model was intriguing given initial reports by @xcite and @xcite that the systems v407 vul and hm cnc ( respectively ) displayed x - ray and optical light curves that were nearly anti - phased . however , @xcite were able to show that in fact timing errors had crept into the previous optical / x - ray analyses , and that for both systems the optical phase leads the x - ray phase by roughly 0.2 . a roughly - constant phase offset suggests a common physical cause for both systems , but the corrected result presents a severe problem for the ui model . in the ui model it is expected that the x - ray - bright footprints of the magnetic field lines that carry the currents should be nearly on the meridian that passes through the line of centres between the two stars . the expected ui geometry thus yields optical light curves anti - phased from the x - ray , assuming that the optical modulation comes from the heated face of the secondary star . the @xcite result that the optical leads the x - ray in phase by @xmath20.2 for both systems implies that the x - ray spots are rotated more than 90@xmath3 away from the line of centres _ for both systems_. first , it appears on the face of it that this can not plausibly be a stable minimum energy plasma current configuration in the co - rotating frame for both systems , but even if so , this geometry places the secondary star below the horizon of the x - ray spots thus removing the presumed heating source of the face of the secondary star . @xcite discuss these and further difficulties with the model , and conclude it does not apply to these systems . the final proposed model is the direct impact model . @xcite showed that if initial roche - lobe contact and the onset of mass transfer in the am cvn systems occurs at a short enough orbital period , the accretion stream can impact the surface of the primary white dwarf directly , and thus not form a disc . @xcite applied this to suggest that the observations for v407 vul could be explained by just such an algol - like model . marsh & steeghs ( 2002 ) proposed that the bulk of the accreted matter would thermalize below the photosphere , mix , and boil to the surface at a temperature cool enough to emit soft x - rays . the temperature excess should broaden in latitude as additional atmospheric material is entrained in the equatorial flow that would be established , and the peak emission energy should shift toward the optical with distance from the upwelling point .    because no hard x - rays are detected from these systems , for this model to be viable the accretion stream must have sufficient ram pressure to advect the specific kinetic energy of the accreted material below the photosphere before thermalizing . in paper i we simulated the system geometry @xcite proposed for v407 vul and found using sph simulations that this constraint is satisfied . the density in the accretion stream is roughly gaussian about the midpoint of the flow , with a 1-@xmath4 half width of 23 km in physical units at the impact point . the contours of mass flux for the impact region on the surface of the primary white dwarf are well approximated by a normal bivariate ( 2-d gaussian ) distribution . fitting the isoflux contours , we found @xmath5 km ( in longitude along the equator ) and @xmath6 km ( in latitude ) . the 3-@xmath4 ellipse encloses 99% of the mass flux , and has an area of @xmath7 or a fraction of the stellar surface of @xmath8 . the mean mass flux within this area is @xmath9 for @xmath10 , the mean density of the stream is @xmath11 , and the velocity of the accreted mass is roughly @xmath12 . within the direct impact model , the location of the accretion stream impact spot is fixed in the binary co - rotating frame ; the spot is on the accreting star s equator , and has an _ impact angle _ with respect to the line of centres that is a function of the mass ratio of the system and the radius ( mass ) of the primary white dwarf . for most physically - plausible system parameters yielding dynamically - stable mass transfer , the impact point has an impact angle roughly in the range 80 to 120@xmath3 @xcite . in this geometry , the altitude of the accretion stream as viewed from the impact point is @xmath13 , the bulk of the specific angular momentum is horizontal to the local stellar surface , and we expect that an equatorial flow will rapidly be established . the accreted material will penetrate far below the photosphere to a depth where the envelope pressure is of order the ram pressure . for the mass transfer rates typical of am cvn systems of such short orbital period ( @xmath14 ; * ? ? ? * ) , we expect @xmath15 , a pressure found at a depth of @xmath16 in a typical @xmath17 white dwarf of @xmath18  k @xcite . in a non - interacting white dwarf of this mass and temperature , the temperature , density , and degeneracy parameter at this depth are @xmath19 k , @xmath20 , and @xmath21 , respectively . the fractional mass of the layer exterior to this depth is @xmath22 , or roughly the mass accreted in a few days to a month at the @xmath23 rates quoted above . if we assume instead that the stream penetrates to a depth of where the envelope pressure is an order of magnitude higher than the ram pressure , @xmath24 , we find the depth to be @xmath25 km or equivalently @xmath26 in mass . the temperature , density , and degeneracy parameter at this depth in the unperturbed model are @xmath27 k , @xmath28 , and @xmath29 , respectively . the observed @xmath265 ev x - ray source temperature @xcite is thus consistent with upwelling of material from these depths . the thermalized temperature of the accreted gas ( @xmath25 kev ) will be orders of magnitude higher  and the density orders of magnitude lower  than the surrounding environment , leading to energetic convective overturn . in the original @xcite model , the authors proposed that the relative phasing of the optical and x - ray light curves might be explained by a temperature profile around the equatorial band characterized by a high temperature at the upwelling point and then a decreasing photospheric temperature as the angle from the impact point increases . because the equatorial flow is in the same direction as the orbital motion , a distant observer will see the characteristic temperature _ increasing _ with time , culminating with the passing of the x - ray - bright spot across the face of the star  this is what causes the optical light curves to lead the x - ray light curves in phase .    in the direct impact model the highest photospheric temperature should be some distance downstream of the accretion stream impact point where the upwelling occurs , and the range in latitude over which the photospheric temperatures exceed the nominal mean stellar photospheric temperatures should be of order but larger than the width of the accretion stream at impact . using the results obtained in paper i , the 3-@xmath4 profile enclosing 99% of the mass flux has a full width of @xmath2138 km , or equivalantly @xmath30 on a white dwarf of radius @xmath31 ( @xmath32 ) . the length of the 3-@xmath4 profile is 980 km , or equivalantly @xmath33 . we expect that the further downstream from the impact point , the broader the zone should be in which the temperature is above the nominal mean photosphere , and the smaller that temperature excess should be . we wrote a code in idl that models a temperature distribution on a cylindrical wall that represents the equatorial band of the accreting white dwarf . we take the latitude limits of this cylinder to be @xmath34 for the calculations here . a schematic representation of the computational layout is given in figure  [ fig : grid ] . the actual computational grid has 360 elements in longitude ( @xmath35 ) , and 100 elements in latitude ( @xmath36 ) . because of the symmetry with respect to @xmath37 , we do not need to explicitly calculate the contributions for @xmath38 since in our model @xmath39 . although integration over a spherical surface would yield slightly different results for very low inclination angles , it would not change the overall results that we present in this preliminary study .    in our idl code , we integrate numerically over the visible grid elements ( @xmath40 ) at each phase . we compute for each element the blackbody spectral energy distribution for that element s temperature , correct for projection angle ( including inclination ) and limb darkening , and then sum over all visible elements . finally , the contribution from the rest of the stellar surface is included , and the integrated spectral energy distribution is convolved with the sloan @xmath41 filter transmission functions for direct comparison with the observational results compiled by @xcite . for the x - ray band , we simply assume a window with 100% transmission from 20 to 60   roughly the sensitivity band of the _ chandra _ letg spectrograph used by @xcite in his study of hm cnc .    knowing that neither system is eclipsing , but otherwise uncertain about the system inclinations , we assume @xmath42 in the model . at this inclination the cylindrical approximation for the equatorial band is still valid . as the inclination decreases from @xmath43 , the fraction of the visible surface made up by the equatorial band declines and the temperature contrast or latitude limits of the band must increase to yield the observed light curve amplitudes . at very low inclinations the visibility drops for the hemisphere below the equator , and the lower latitudes of the equatorial band may be hidden ( see figures 2 and 3 ) , again requiring an increased temperature contrast or increased latitude limits to recover the observed optical light curve amplitudes . there are countless models for the equatorial - band temperature distribution that could be employed . for example , we initially implemented a model in which there was a simple linear temperature decline with distance from a compact high - temperature ( x - ray bright ) spot , but quickly discarded that model as it produced synthetic light curves with far too much x - ray luminosity over too great a range in longitude / phase . our next model included a single compact high - temperature ( x - ray emitting ) spot . downstream from this spot the temperature distribution is a decreasing function of both distance from the high - temperature spot , and from the equatorial symmetry plane . in this model the temperature falls off exponentially with longitude , and then falls off in latitude approximately as a gaussian function that has a width that increases with longitude . we use the standard kernel function of smoothed particle hydrodynamics to approximate the gaussian , and let our temperature distribution within the latitude limits of the equatorial band be given by the following expression :    @xmath44\ ] ]    where @xmath45 is the assumed background temperature of the equatorial band of the primary star , and @xmath46 is the peak temperature of the distribution ( @xmath47 ev ) . the temperature of the remainder of the star is @xmath48 . the function @xmath49 is the sph kernel function @xcite , which approximates a gaussian distribution but is identically zero beyond twice the smoothing length @xmath50 :    @xmath51 ^ 3 ,                                    & $ 1 \leq { \theta \over h } < 2$,\cr                                   \ \cr                           0 & $ { \theta \over h } \geq 2$.\cr                           } \ ] ]    physically we expect that the thermalized energy from the accretion stream should spread out in latitude as a function of distance from the upwelling point , and we assume the functional form    @xmath52    where in the results given below we use @xmath53 and @xmath54 . unlike our original linear - decline model , this one immediately gave results that reproduced the general observational results , including the on / off x - ray pulse shape and the optical / x - ray light curve phase offset . however , as we discuss below , a model with a single high - temperature spot does not reproduce the width of the x - ray pulses , and therefore we were forced to introduce a second x - ray spot @xmath265@xmath3 upstream of the original spot , characterised by temperature @xmath55 . so in practice we typically set @xmath56 and @xmath57 , where the results are not terribly sensitive to the choice for @xmath58 . finally , because the e - folding length @xmath59 required to recover the observed phase offset is such that the temperature of a grid element obtained using equation ( 1 ) extended beyond @xmath60 exceeds that for the same element directly , we set the temperature of each grid element to be the maximum of @xmath61 or @xmath62 . we show in figures [ fig : tempconthm ] and [ fig : tempcontv ] contour plots of the temperature distributions used below to produce synthetic lightcurves representing hm cnc and v407 vul , respectively . kk , and the contour levels are 100001 k , 100.5 kk , 110 kk , 150 kk , and 200 to 700 kk in steps of 100 kk . the white lines plotted over the second cycle show the function @xmath63 . the high - temperature feature at @xmath2@xmath64 and @xmath2@xmath65 is the impact point x - ray spot . ] kk , and the contour levels are 25010 k , 30 kk , 50 kk , and 100 to 700 kk in steps of 100 kk . the white lines plotted over the second cycle show the function @xmath63 . the high - temperature feature at @xmath2@xmath64 and @xmath2@xmath65 is the impact point x - ray spot . ] the phase - folded light curves compiled by @xcite are shown in figure  [ fig : bar07 ] . the authors reanalysed the _ chandra _ x - ray data for v407 vul obtained on 2003 february 19 ( solid line ) and 2003 november 24 ( dashed line ) . the x - ray data for hm cnc were digitized and fit from @xcite . barros et al .  ( 2007 ) obtained the optical light curves using the ultracam triple - beam ccd camera @xcite on the 4.2-m william hershel telescope in 2003 may and on the vlt in 2005 november . in the course of their analysis , the authors discovered a likely error in the @xcite v407 vul ephemeris , and in the optical timings of ( * ? ? ? * ; * ? ? ? * see barros et al . ( 2007 ) for full details ) .    in both v407 vul and hm cnc , the x - ray light curves are characterised by a pulse shape that rises steeply from zero detected flux in the `` off '' state to maximum flux , and then falls more gradually from peak light . additionally , in both systems the phase offset between the x - ray and optical data is roughly 0.2 in phase , with the optical peak flux leading the x - ray . contrasting the systems , the @xcite data show that , whereas the fractional amplitudes of the v407 vul optical light curves increase dramatically from the @xmath66 to @xmath67 bands , the hm cnc fractional amplitudes are relatively constant across all optical wavelengths . in addition , whereas the v407 vul optical light curves are very nearly pure sinusoids , the hm cnc light curves reveal a second - harmonic that has an amplitude 1520% that of the first harmonic . the optical spectrum of v407 vul is dominated by a late - g ( or early - k ) type stellar spectrum @xcite . barros et al .  ( 2007 ) suggest this star is either a line - of - sight coincidence in a dense field , or that it is the distant member of a hierarchical triple system , but that it can not be the direct cause of the optical or x - ray variations . whatever the ultimate determination , the longer - wavelength light curves are increasingly diluted from the flux from this star , causing at least in part the decreasing amplitude that is observed in longer - wavelength filters . we note that we omit the @xmath68-band light curves from the analysis below because of the increased noise in the case of hm cnc , and because of the dilution by the g star in the case of v407 vul . the primary goal of this project was to explore whether a plausible surface temperature distribution consistent with the direct impact model could give rise to the light curves observed for hm cnc and v407 vul , and in particular , whether this model could yield the phase offset and pulse shapes observed . the results we present in this section are representative of the model output , but we caution that the solutions are not unique . we begin by exploring the geometry and nature of the source of the x - ray light curves in both systems . in our model , we first assumed a small x - ray - bright region presumably just downstream of the impact point , and that the temperature distribution fell off exponentially as a function of longitude , while broadening in latitude . in figure [ fig : xraysonly ] we show a single cycle of the the x - ray phase - folded light curves for both systems with a simulated light curve also shown with and without limb darkening included . for v407 vul , we include only the 2003 february 19 data for clarity . the simulations shown here are of a single 65 ev ( 750 kk ) hot spot on the grid described above . the dotted line was calculated assuming no limb darkening , the solid line was calculated assuming a linear limb darkening law @xmath69 , and the short - dashed line was computed with a quadratic limb - darkening law @xmath70 , where @xmath71 and where we assume @xmath72 in both relations . even with this extreme limb darkening coefficient , it can be seen that for both systems the observed light curves rise more steeply than the synthetic curve . one possible effect not included that could give the observed steep slope is a stellar surface that is non - spherical . the hydrodynamic simulations of direct impact accretion between polytropic stars by ( * ? ? ? * see , e.g. , their figure 1 ) show just such an effect . ( solid line ) , and a quadratic limb darkening law again with @xmath72 ( short - dashed line ) . also shown is the light curve with the quadratic relation offset by 0.21 in phase ( long - dashed line ) , demonstrating that the final decline in the x - ray light curve of hm cnc is well matched by a second compact x - ray spot . ] figure [ fig : xraysonly ] also demonstrates that the pulse width in both systems is wider than would be expected for a single x - ray - bright spot rotating into and out of view . interestingly , both systems show a gentler slope during part of their decline , and that slope appears to be roughly the same for both systems . in hm cnc , the final decline is observed to be as steep as the initial rise , which is to say as steep as a single bright spot on the stellar surface . to demonstrate this we copied a portion of the single - spot fully - limb - darkened curve and shifted it by 0.21 in phase , and have shown this in figure [ fig : xraysonly ] as the long - dashed line . we find that to match the general shape of the x - ray light curves , we were forced to include a second compact x - ray spot some 60 to 65@xmath3 upstream of the primary x - ray spot and small enough that it does not contribute significantly to the optical flux . we speculate that the primary spot is the upwelling point , and the secondary spot is the region around the impact spot . while most of the specific kinetic energy of the accreted mass is advected and thermalized below the photosphere , the low - density outer sheath of the accretion stream ( @xmath73-@xmath4 from the midline ) will have @xmath74 and this should yield an x - ray emitting halo around the impact point . there will be a strong surface flow towards the impact spot in this region with speeds of order the sound speed as a result of the higher - density core of the accretion stream sweeping any and all atoms deep below the photosphere ; effectively , the high - density core of the impact spot is a hole in the surface of the star into which the surrounding material will flow . the direct impact hydrodynamic simulations of @xcite , while computed for a mass ratio of @xmath75 and thus higher than thought to apply to hm cnc and v407 vul , do indicate an upwelling point substantially displaced from the impact point . figure [ fig : hmcnc ] shows a sample model result for hm cnc plotted over the data . to determine suitable model parameters for each system we first find those that yield the best match `` by eye '' to the x - ray light curve , and then explore the remaining parameters to find a reasonable match to the optical light curves . the parameters for the hm cnc model shown in figure [ fig : hmcnc ] are given in table 1 and the contour plot of the temperature distribution is shown in figure [ fig : tempconthm ] , but again we emphasize that this set of parameters is almost certainly not unique and should be taken as representative , not definitive . as an example of this , no set of parameters that we found could cleanly reproduce the excess flux in the optical between phases 0.8 and 1.0 , but by artificially doubling the temperature of each grid element between longitudes 220@xmath3 and @xmath76 , we obtain a result which is a better approximation to the data . the physical source of this excess flux is not clear . if the broad outlines of the current model are correct , then the source can not be the heated face of the secondary star , as the phasing is is almost exactly wrong . for an assumed impact angle of roughly @xmath77 the upwelling point is at an angle of @xmath78 and the excess optical flux is at an angle of @xmath79 , or roughly when the back side of the secondary is towards the observer . the hydrodynamic simulations of @xcite show some evidence for large amplitude standing waves on the surface of the accreting star after the outer layers are spun up and perhaps this could be related to the observed excess flux . additional work will be required to reveal the source of this radiation . there are some physical conclusions we can make from the simulations that most closely approximate the data . the synthetic x - ray light curve is a reasonable approximation to the data , but the data rise more steeply to maximum , perhaps because the surface of the star is slightly non - spherical . our experiments clearly show that two distinct x - ray spots are required to reproduce the observed x - ray pulse shape  a single extended high - temperature ridge results in slopes of rise and fall that are shallower than those observed . observationally , the optical light curves all have nearly the same fractional amplitude regardless of filter , implying that the underlying source must have a very high temperature that is fluctuating , such that the optical passbands are far out on the rayleigh - jeans tail of the spectral energy distribution . @xcite suggest that compressional heating should result in these stars having temperatures @xmath80 kk , but we find fractional amplitudes of only @xmath20.04 if we assume the entire stellar surface has a temperature of 100 kk , and not just the equatorial band base temperature . .v407 vul , hm cnc model parameters . [ cols="<,<,<,<",options="header " , ]      figure [ fig : v407vul ] shows a sample model result for v407 vul , obtained using the parameters given in table 1 . the temperature distribution is shown in figure [ fig : tempcontv ] . the direct impact model provides a reasonable approximation to the observations . as discussed above , the optical spectrum of v407 vul is dominated by a g - type stellar spectrum , so the amplitude variations observed are not intrinsic to the direct impact scenario , and as a result are not reproduced by the simulations . @xcite note that the amplitudes are consistent with a constant g - star spectrum plus a variable blue component . for this model we show the @xmath81 and @xmath66 synthetic light curves as dotted lines . in order to recover the observed phase offset for this star , we must use an e - folding length of @xmath82 , but this then yields excess flux in the x - ray light curve in advance of the observed rise . light curves reasonably well ( solid lines ) . the g - star spectrum dilutes the observed amplitudes in longer passbands . ] we have shown that a variant of the direct impact model of @xcite can reproduce the general features observed from the ultracompact binary am cvn systems hm cnc and v407 vul . these preliminary results are promising , but should be taken as representative only . improved fits to the observations would be possible by iterating to a minimim @xmath83 using differential corrections to the model parameters as is implemented in the wilson - devinney code @xcite , or to allow the temperatures of the current grid points to vary freely ( with or without constraints ) and use a genetic algorithm or related technique ( e.g. , * ? ? ? * ) to find solutions for the temperature distributions that match the light curves of the individual systems . the solutions found in this way would not be unique as the number of grid points exceeds the number of data points , but additional physical insight might be gained through such a study . this is similar to the `` firefly '' approach used by @xcite and @xcite in their studies of eclipse mapping of polar accretion streams . the primary results of this work are    1 . the direct impact model of @xcite is a plausible mechanism to produce the observed x - ray and optical light curves observed in the ultracompact am cvn binaries hm cnc and v407 vul . using a reasonable guess for the general parameters of a temperature distribution model , the optical / x - ray phase offset is reproduced , as well as the general shapes of the light curves themselves . subtle differences between the model results and the observations indicate that the model is still incomplete . 2 .   in both systems the x - ray light curves rise significantly more steeply than would be predicted for a single compact hot spot on the surface of a spherical star , even assuming an extreme limb - darkening law and coefficient for the x - ray band . we speculate that the surface of the star may be non - spherical as a result of the ongoing accretion , and that the upwelling x - ray region in particular may have a photospheric radius slightly larger than than the average for the star . the hydrodynamic simulations of @xcite appear to support this scenario . the x - ray light curve from system hm cnc and v407 vul can only be reproduced by including a second x - ray - bright spot @xmath84 upstream of the primary ( brighter ) x - ray region . to fit the observed light curve , there must be two distinct spots  a high - temperature ridge spanning this range in longitude will not yield an x - ray light curve that reproduces the observed pulse shape . we speculate that the secondary spot is the actual impact point of the accretion stream , and that the primary spot is the region where the upwelling occurs ( again cf . the x - rays from the impact region most likely originate from the low density outer regions of the accretion stream for which @xmath74 . we might expect that the x - ray emission from the turbulent upwelling point might be highly variable  considerably moreso than the x - ray emission from the impact point  and a detailed statistical analysis of the existing x - ray time series might shed light on the viability of this proposed model . 4 .   for the accretion rates predicted for these two systems @xcite , the ram pressure is sufficient to penetrate to a depth of order 100 km below the photosphere or @xmath85 to @xmath86 in mass depth  i.e. , @xmath87 , where @xmath88 is the mass contained within radius @xmath89 . the temperatures at these depths in the envelope of a 20 kk white dwarf model are characteristic of the 65 ev temperatures measured from x - ray spectroscopic observations . the observed optical light curves of hm cnc all have roughly the same amplitude . to obtain this in the simulations , the base temperature of the equatorial band must be of order 100 kk , but the rest of the star must be at a lower temperature to reproduce the observed fractional amplitudes in the optical . @xcite suggest that compressional heating of the accretor should lead to mean photospheric temperatures of @xmath80 kk  marginally consistent with our simple model results . the amplitudes of the optical - band light curves from v407 vul increase dramatically from the @xmath68 band to the @xmath67 band . because the optical light is diluted with the g - star spectrum @xcite it is unknown how much of the amplitude trend in v407 vul is the result of this dilution , and how much might be intrinsic to the accreting star s temperature distribution . 6 .   in order to recover the observed phase offset between the optical and x - ray light curves in v407 vul , the model yields non - zero x - ray flux significantly in advance of the observed rise . this suggests that the actual temperature distribution resulting from direct impact accretion may be more complex than the simple functional form we have assumed for this preliminary study . while the current results are intriguing , what remains to be done is to model the hydrodynamics of the direct impact scenario directly as it applies to the envelope and atmosphere of the accreting white dwarf . direct hydrodynamical simulations in 3d should allow the estimation of the photospheric temperature on the surface of the primary star , and from that it is straightforward to calculate an estimate of the phased light curves that should result . the hydrocode pluto @xcite is well suited to such problems , and we are in the initial stages of implementing this project . we anticipate that explicit hydrodynamical simulations will not only reveal the true structure and dynamics of the outer layers of the accreting white dwarfs in hm cnc and v407 vul , but also the nature of the closely related system es cet @xcite where time - series spectroscopy reveals unusual features that may result from a truncated accretion disc , possibly resulting from a grazing impact of the accretion stream with the white dwarf atmosphere as the system makes the transition from a direct impact system to a normal am cvn accretion disc . thanks to susana barros for kindly sending an electronic copy of the observational phase - folded light curves . thanks also to gijs nelemans for insightful comments on a draft of this manuscript , and to danny steeghs and paul groot for useful conversations . thanks finally to the reviewer for comments that improved the manuscript . the author was on sabbatical from the florida institute of technology during the period of this work , and thanks radboud university and the netherlands organisation for scientific research ( nwo ) for support through visitor s grant 040.11.046 . this research has made use of nasa s astrophysics data system bibliographic services . espaillat c. , patterson j. , warner b. , woudt p. , 2005 , pasp , 117 , 189 gokhale v. , peng x.  m. , frank j. , 2007 , apj , 655 , 1010 goldberg d.  e. , 1989 , genetic algorithms in search , optimization and machine learning , addison - wesley , reading                marsh t.  r. , nelemans g. , steeghs d. , 2004 , mnras , 350 , 113 mignone a. , bodo g. , massaglia s. , matsakos t. , tesileanu o. , zanni c. , ferrari a. , 2007 , apjs , 170 , 228 monaghan j.  j. , lattanzio j.  c. , 1985 , a&a , 149 , 135                      steeghs d. , marsh t.  r. , barros s.  c.  c. , nelemans g. , groot p.  j. , roelofs g.  h.  a. , ramsay g. , cropper m. , 2006 , apj , 649 , 382 strohmayer t.  e. , 2003 , apj , 593 , l39 strohmayer t.  e. , 2004 , apj , 610 , 416 strohmayer t.  e. , 2005 , apj , 627 , 920 strohmayer t.  e. , 2008 , apj , 679 , l109 warner b. , robinson e.  l. , nather r.  e. , 1971 , mnras , 154 , 455
Abstract:

Answer: the interacting binary white dwarf ( am cvn ) systems hm cnc and v407 have orbital periods of 5.4 min and 9.5 min , respectively . the two systems are characterized by an `` on / off '' behaviour in the x - ray light curve , and optical light curves that are nearly sinusoidal and which lead the x - ray light curves in phase by about 0.2 in both systems . of the models that have been proposed to explain the observations , the one that seems to require the least fine tuning is the _ direct impact _ model of marsh & steeghs ( 2002 ) . in this model , the white dwarf primary is large enough relative to the semi - major axis that the accretion stream impacts the surface of the primary white dwarf directly without forming an accretion disc . marsh & steeghs ( 2002 ) proposed that in this situation there could be a flow set up around the equator with a decreasing surface temperature the further one measured from the impact point . in this study , we estimate the light curves that might result from such a temperature distribution , and find them to be reasonable approximations to the observations . one unexpected result is that two distinct x - ray spots must exist to match the shape of the x - ray light curves . [ firstpage ]    accretion , accretion discs  binaries : general , close  hydrodynamics  novae , cataclysmic variables  stars : individual(hm cnc ; v407 vul )  white dwarfs  x - rays : binaries  x - rays : stars .


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 @xmath0 scattering at hera interactions between photons of small virtuality @xmath2 and protons produce jets of high transverse energy @xmath3 . the presence of this hard scale @xmath3 allows the application of perturbative qcd to predict cross sections for the production of two or more high-@xmath3 jets which can be confronted with experimental data . this offers the opportunity to test qcd and to constrain the structure of the colliding particles .    at leading order ( lo ) qcd two distinct processes are responsible for the production of jets . in the direct photon process , the photon interacts as a point - like object with a parton in the proton , whereas in the resolved process the photon acts as a source of partons which then scatter with the partons coming from the proton . both lo processes are characterized by having two outgoing jets of large transverse energy . studies of dijet photoproduction at hera have shown that both classes of processes are present for the case of quasi - real photons , i.e. photons of extremely small virtuality @xmath4 @xcite . the comparison of theoretical predictions @xcite with these @xcite and more recent experimental data @xcite have given us some confidence that the parton distribution functions for the real photon available in the literature @xcite are consistent with the dijet production data . the parton content of photons with virtuality @xmath5 is reasonably constrained by data from deep inelastic scattering @xcite . unfortunately this is not the case for a photon target with non - zero , although small , virtuality @xmath2 . the only measurement for the virtual photon structure function available so far has been performed by the pluto collaboration at petra @xcite . they measured the structure function @xmath6 for @xmath7 gev@xmath8 and @xmath9 gev@xmath8 as a function of @xmath10 ( @xmath11 is the virtuality of the probing virtual photon , whereas @xmath2 always denotes the virtuality of the probed virtual target photon ) . more and better data should come from lep2 @xcite . on the theoretical side several models exist for describing the @xmath11-evolution equations of the parton distributions and the input distributions at @xmath12 with changing @xmath2 @xcite . these constructions use essentially the same methods as have been applied for the parton distributions of real photons . this way some smooth behavior towards @xmath5 , where previous results for the real photon should hold , is guaranteed . in @xcite this construction allows a calculation of the parton distribution functions ( pdf ) for virtual photons in lo and nlo . it incorporates a purely perturbative contribution and a non - pointlike hadronic contribution . unfortunately these pdf s for @xmath13 are not available in a form that parametrizes the @xmath11 evolution . such parametrized pdf s for virtual photons have been presented recently by two groups @xcite , but unfortunately only in lo . these two models differ somewhat in their method of extrapolation to @xmath14 , in the choice of the input distribution and the choice of the input scale @xmath12 . furthermore the pdf s of glck , reya and stratmann ( grs ) @xcite present distributions only for @xmath15 flavors , whereas schuler and sjstrand ( sas ) give pdf s for @xmath16 flavors @xcite . we expect these two structure function sets to be detailed enough so they can be tested in photoproduction experiments with virtual photons . first preliminary data have been presented by the zeus @xcite and the h1 @xcite collaborations . zeus studied the dijet cross section for @xmath17 gev in the range of @xmath18 . events with @xmath19 are assumed to be dominated by the direct process whereas events with @xmath20 give the resolved component . then the ratio of the direct - enriched and resolved - enriched cross section is measured as a function of @xmath2 for @xmath21 gev@xmath8 . in the h1 experiment the inclusive one - jet cross section is measured as a function of @xmath3 and rapidity @xmath22 for various @xmath2 bins @xcite . it is expected that the resolved component decreases relative to the contribution from the direct photon processes as the virtuality of the photon increases . some theoretical studies of the single inclusive and dijet inclusive cross section in lo have been presented recently @xcite . it is well - known that in nlo calculations the distinction between direct and resolved photoproduction becomes ambiguous . in this order , a large contribution in the nlo direct cross section is subtracted from the direct component and combined with the lo resolved term thus producing the scale ( @xmath23 ) dependence of the pdf s of the photon . therefore both components are related to each other through the factorization scale @xmath24 at the photon leg which determines the part of the nlo direct contributions to be absorbed into the resolved component . the @xmath24 dependence of the remaining nlo direct contribution cancels to a large extent against the dependence in the resolved cross section coming from the photon structure function . this connection has been worked out in detail @xcite and studied numerically for real photoproduction with @xmath25 @xcite . it is clear that this cancellation of the scale dependence must take place also for @xmath14 . so , for a consistent calculation up to nlo one needs to superimpose both contributions , the nlo direct and at least the lo resolved cross section , both computed with the same scale . the subtraction of the large contribution in the direct cross section which is shifted to the resolved term is defined only up to finite , non - singular terms in the limit @xmath26 . these finite terms may be fixed in a way that a smooth behavior towards the limiting case of real photons ( @xmath25 ) is achieved where these finite terms are usually defined in the @xmath27 subtraction scheme . for a complete nlo calculation we must include the nlo hard scattering parton - parton cross sections for the direct ( here one of the ingoing partons is the photon with virtuality @xmath14 ) and for the resolved process together with the two - loop evolved structure function of the proton and photon with virtuality @xmath14 . in this paper we shall work out the subtraction in the nlo direct contribution and superimpose the remaining direct contribution and the resolved cross section up to nlo . we study this cross section and the two contributions as a function of @xmath2 for various inclusive one- and two - jet cross sections . of particular interest is the behavior of the resolved component as a function of @xmath2 and the question at which @xmath2 the sum of direct and resolved cross sections approaches the unsubtracted direct cross section which one expects to give the correct description at sufficiently large @xmath2 . the outline is as follows . in section 2 we describe how to subtract the singular terms in the nlo direct cross section connected with the collinear singularity of the @xmath28 contribution . here we also define the finite terms which depend on the subtraction scheme used for @xmath25 . our result concerning the @xmath2 dependence of various one- and two - jet cross sections are presented in section 3 . in this section we also compare with the preliminary data from zeus @xcite . section 4 contains a summary and an outlook for future investigations . the nlo corrections to the direct process become singular in the limit @xmath25 . for real photoproduction these photon initial state singularities are usually evaluated , i.e. regularized , with the dimensional regularization method in which they are finite for @xmath29 . the singular contributions appear as poles in @xmath30 multiplied by the splitting function @xmath31 @xcite . these initial state singularities are absorbed into the pdf of the real photon ( @xmath5 ) . with no further finite terms subtracted ( for @xmath32 ) this defines the @xmath27 factorization scheme which must be consistently applied also for the nlo evolution of the photon pdf . for @xmath14 the corresponding contributions appear as terms proportional to @xmath33 , @xmath34 being the c.m . energy of the photon - parton subprocess . these terms are finite as long as @xmath14 and can be calculated with @xmath35 dimensions . since for small @xmath2 these terms are large they are absorbed as in the case @xmath5 into the pdf of the virtual photon which is present in the resolved cross section . concerning finite terms ( for @xmath26 ) which may also be subtracted together with the singular terms we have the same freedom as in the case @xmath5 . we fix these finite terms such that they agree with the @xmath27 factorization in the real photon case . to achieve this we must compare the contributions originating from the photon initial state singularities in the two cases @xmath5 and @xmath14 , which we shall do in the following . .classification of @xmath36 matrix elements [ cols="^,^,^",options="header " , ]     the procedure for virtual photons with virtuality @xmath2 is completely analogous . first one calculates the @xmath36 matrix elements , but now with @xmath14 and decomposes them into terms with the characteristic denominator from @xmath28 splitting which become singular in the limit @xmath26 . these terms after phase space integration up to a cut - off @xmath37 as in ( 1 ) have the same structure as ( 2)(5 ) with the lo parton - parton matrix elements factored out . the integration can be done with @xmath38 since @xmath14 . the phase space in ( 1 ) then contains the additional @xmath2-dependent factor @xmath39 @xmath40 which reduces to @xmath41 for @xmath42 . the factor @xmath43 in the equations ( 2)(5 ) takes the simple form @xmath44 which is singular for @xmath5 as to be expected . this singularity is absorbed into the pdf of the virtual photon with virtuality @xmath2 . instead of ( 8) the subtraction term is : @xmath45 after this subtraction the remaining finite term ( for @xmath26 ) in @xmath43 yields @xmath46 in addition to the singular term @xmath47 we have subtracted in ( 12 ) two finite terms in order to achieve in ( 13 ) the same result as in ( 9 ) for @xmath5 . therefore we call this form of factorization the @xmath27 factorization for @xmath14 . it is defined by the requirement that the remaining finite term is equal to the finite term in ( 2)(5 ) with @xmath43 replaced by @xmath48 as was obtained in the calculation for real photons . with this definition of factorization in the case @xmath14 we make sure that we obtain the same nlo corrections as in @xcite , where @xmath42 , when in the @xmath14 calculation we choose @xmath2 extremely small . so , for the actual calculations we apply the formulas ( 2)(5 ) where @xmath43 is now given by ( 13 ) . we note that for averaging over the spin of initial gluons and photons we apply a factor @xmath49 , that gives rise to some additional finite terms which should be included in ( 9 ) when expanded in @xmath30 . we have taken these terms into account also in our calculation . this completes the calculation of the contribution from the photon initial state singularity . it is clear that in nlo the pdf for the virtual photon must be given also in the @xmath27 factorization scheme . in ref . @xcite the pdf is constructed in the so - called dis@xmath50 scheme , which is defined as for real photons ( @xmath25 ) . this distribution function is related to the @xmath27 scheme pdf in the following way @xcite : @xmath51 where for @xmath52 or @xmath53 : @xmath54 \\      { \delta}f_{g/{\gamma}}(x , m_{\gamma}^2 ) & = & 0 \quad . \nonumber\end{aligned}\ ] ] if the pdf in this scheme is used to calculate the resolved cross section one must transform the nlo finite terms in the direct cross section . this produces a shift of @xmath55 as given in ( 13 ) by the same expression as ( 15 ) . then the relation is : @xmath56 \quad .\ ] ]    all other singular terms in the real corrections , i.e.  final state singularities and the contributions from parton initial state singularities have been calculated by graudenz in connection with the nlo corrections for jet production in deep inelastic ep scattering . they can be taken from his work @xcite together with the virtual corrections up to @xmath57 . ( for related work see @xcite . ) all these contributions are calculated in dimensional regularization . the appearing singularities in @xmath58 cancel when all singular terms are added . the remaining finite terms enter the nlo corrections for the jet cross sections . these finite terms depend on the phase space slicing parameter @xmath37 which is introduced to separate the singular regions of final and initial state infrared and collinear divergences . in this section , we present some characteristic numerical results for one- and two - jet cross sections as a function of the virtuality @xmath2 . we consider the contributions of the direct and resolved components and their sum . partly we shall compare the nlo results with lo predictions . the input for our calculation is the following . we have chosen the cteq3 m proton structure function @xcite which is a nlo parametrization with @xmath27 factorization and @xmath59 mev . this @xmath60 value is also used to calculate @xmath61 from the two - loop formula at the scale @xmath62 . the factorization scales are also @xmath63 . we also need the parton distribution of the virtual photon @xmath64 . for this we choose either the grs @xcite set or the sas1 m set of schuler and sjstrand @xcite . both sets are given in parametrized form for all scales @xmath65 so that they can be applied without repeating the computation of the evolution . unfortunately both sets are given only in lo , i.e. the boundary conditions for @xmath66 and the evolution equations are in lo . in @xcite pdf s for virtual photons have been constructed in lo and nlo . distinct differences occur for larger @xmath2 and @xmath67 which is mainly due to the different nlo perturbative boundary condition at @xmath66 , which does not exist for the real ( @xmath42 ) photon structure function . since neither of the two pdf s is constrained by empirical data from scattering on a virtual photon target we consider these lo distribution functions as sufficient for our exploratory studies on jet production and treat them as if they were obtained in nlo . in particular , we shall use the transformation formulae ( 14 ) for going from the dis@xmath50 to the @xmath27-scheme pdf , which makes sense only in nlo . then the parametrization @xcite is considered as the parametrization in the dis@xmath50 scheme . in the pdf of grs the input scale is @xmath68 gev and the restriction @xmath69 is implemented as to fulfill the condition @xmath70 . in addition the pdf of grs is constructed only for @xmath71 flavors . the production of the heavier @xmath72 and @xmath73 quarks is supposed to be added as predicted by perturbation theory of fixed order with no active @xmath72 and @xmath73 quarks in the proton and photon pdf s . in lo this amounts to adding the processes @xmath74 and @xmath75 to the cross section , keeping @xmath76 . since in this work we are primarily interested in studying the sum of the direct and resolved contributions and the influence of the consistent subtractions of the nlo direct part we refrain from adding the lo or nlo cross sections for direct heavy quark production as suggested in @xcite . so , our investigations in connection with the grs parametrization of the virtual photon pdf are for a model with three flavors only . for consistency we take also @xmath71 in the nlo corrections and in the two - loop formula for @xmath77 . to overcome this problem we studied the relevant cross sections also with the virtual photon pdf s of @xcite . here the @xmath72 quark is included as a massless flavor in the pdf which undergoes the usual evolution as the other massless quarks except for a shift of the starting scale @xmath12 . this @xmath78 pdf is considered only in the parametrization sas1 m with @xmath79 gev that is in the @xmath27 scheme . the cross sections we have computed are essentially for kinematical conditions as in the hera experiments . there , positrons of @xmath80 gev produce photons with virtuality @xmath2 which then collide with a proton of @xmath81 gev . the momentum transfer to the photon is @xmath82 with @xmath83 . the energy spectrum of the virtual photons is approximated by @xmath84 with @xmath85 being the fraction of the electron energy transferred to the photon , when the virtuality @xmath86 . the momentum of the incoming proton is @xmath87 . the approximation for the virtual photon spectrum is used for the calculation of the resolved and the direct cross section . the expression ( [ ww ] ) factorizes in the cross section for @xmath88 with arbitrary final state @xmath89 if one neglects the longitudinal virtual photon terms . after integration over @xmath2 between @xmath90 with @xmath91 , where @xmath92 is the electron mass , one obtains the weizscker - williams formula as used for calculations with untagged electrons as in @xcite , where @xmath4 dominates . the cross section for the process @xmath88 is then given by the convolution @xmath93 where @xmath94 denotes the cross section for @xmath95 with transversely polarized photons of energy @xmath96 , if the transverse component of @xmath97 is neglected . to have the equivalent conditions as in the zeus analysis we choose @xmath98 and @xmath99 . in the computation of the resolved cross section the approximation @xmath100 is inserted for the energy of the virtual photon , whereas for the direct photon cross section the exact relation for @xmath101 is taken into account through the kinematic relations . all further calculations proceed in the following way . for both , the direct and the resolved cross section , we have a set of two - body contributions and a set of three - body contributions . each set is completely finite , as all singularities have been canceled or absorbed into pdf s . each part depends separately on the phase space slicing parameter @xmath37 . the analytic calculations are valid only for very small @xmath37 , since terms @xmath102 have been neglected in the analytic integrations . for very small @xmath37 , the two separate pieces have no physical meaning . in this case the @xmath103-terms force the two - body contributions to become negative , whereas the three - body cross sections are large and positive . in @xcite we have plotted such a cross sections for direct real photoproduction at @xmath104 . when the two - body and three - body contributions are superimposed to yield a suitable inclusive cross section , as for example the inclusive one- or two - jet cross section , the dependence on the cut - off @xmath37 will cancel . the separation of the two contributions with the slicing parameter @xmath37 is a purely technical device in order to distinguish the phase space regions , where the integrations are done analytically , from those , where they are done numerically . we have checked , by varying @xmath37 between @xmath105 and @xmath106 , that the superimposed two- and three - body contributions are independent of @xmath37 for the inclusive single- and dijet cross sections . first , we study the inclusive single - jet cross section . to calculate it we must choose a jet definition , which combines two nearly collinear partons . we adopt the jet definition of the snowmass meeting @xcite . according to this definition , two partons @xmath107 and @xmath108 are recombined , if @xmath109 , where @xmath110 and @xmath111 are the rapidity and the azimuthal angle of the combined jet respectively , defined as @xmath112 and @xmath113 is chosen as in the experimental analysis . so , two partons are considered as two separate jets or as a single jet depending whether they lie outside or inside the cone with radius @xmath113 around the jet momentum . in nlo , the final state may consist of two or three jets . the three - jet sample contains all three - body contributions , which do not fulfill the cone condition . the cone constraint is evaluated in the hera laboratory system as for real photoproduction ( @xmath5 ) and in the experimental analysis . the calculation of the resolved cross section proceeds as for real photoproduction , i.e. the transverse momentum ( @xmath114 ) of the virtual photon is neglected so that the virtual photon momentum is in the direction of the incoming electron and @xmath100 . the cross section for direct photons , in which the initial state singularity is subtracted , as specified in the previous section , is given in the center - of - mass system @xmath115 and transformed into the hera laboratory system , again neglecting @xmath114 and other small terms proportional to @xmath2 . then the relation between the rapidity @xmath116 of the jet in the c.m .- system and laboratory system is as for real photoproduction , which is @xmath117 this @xmath22 and the corresponding azimuthal angle of the partons in the final state is also used for evaluating the jet definition given above .    in fig . 1a , b , c , the results for @xmath118 are shown integrated over @xmath22 in the interval @xmath119 ( these boundaries are employed in the zeus analysis @xcite ) as a function of @xmath3 for the three values of @xmath120 and @xmath41 gev@xmath8 . the smallest value of @xmath2 has been chosen in such a way that it reproduces the @xmath121 results . inserting @xmath122 gev@xmath8 into the unintegrated weizscker - williams formula corresponds to the one integrated in the region @xmath123 gev@xmath8 . for all three @xmath2 the cross section is dominated by the resolved component at small @xmath3 . near @xmath124 gev the direct component , which is the direct cross section without the subtraction term , denoted by dir@xmath125 , is of the same magnitude as the resolved cross section . as a function of @xmath2 the cross section @xmath118 decreases more or less uniformly in the considered @xmath3 range with increasing @xmath2 . one can see that this decrease is stronger for small @xmath3 as compared to large @xmath3 . in these and the following results the cone radius is @xmath126 . the corresponding rapidity distribution for a fixed @xmath3 will be shown later in fig . 5a , b , c , d for the sas1 m photon pdf .    in fig . 2a , b , c we studied the @xmath22 distribution of the dir@xmath125 contribution at fixed @xmath127 gev and the same @xmath2 values as in fig .  1 in comparison with two approximations , namely the lo cross section and the nlo cross section with @xmath14 only in the photon flux equation ( 17 ) and the rest of the cross section evaluated at @xmath42 as in @xcite . as to be expected this approximation is very good for @xmath122 gev@xmath8 . at the larger @xmath2 however it overestimates the cross section and should not be used . this means that the @xmath2 dependence of the direct part , although the strongest logarithmic @xmath2 dependence has been subtracted , should be taken into account . of course , in the sum of the resolved and the direct cross section the difference is small as long as the resolved part dominates which is true for the smaller @xmath3 s . the lo prediction , which is evaluated with the same structure functions and @xmath61 value as the nlo result , only the hard scattering cross section is calculated in lo , is smaller than the nlo result , the difference decreasing with increasing @xmath2 . of course , this finding depends on the chosen value of @xmath113 . the nlo cross section depends on @xmath113 , whereas the lo curve does not . therefore estimates of the inclusive cross section with lo calculations can be trusted only for large cone radii . it is clear that the resolved and the direct cross sections decrease with increasing @xmath2 for fixed @xmath22 and @xmath3 . it is of interest to know how the ratio of resolved and the direct cross section behaves as a function of @xmath2 . this is a well defined quantity in lo . the variation of this ratio with @xmath2 up to @xmath128 gev@xmath8 for @xmath127 gev and @xmath129 and @xmath130 is plotted in fig . 3a , b , c , d. as we expected this ratio decreases most strongly for @xmath131 , since in the @xmath132 region the resolved component dominates whereas the direct cross section peaks for negative @xmath22 s ( see fig . 2a , b , c ) . if we decrease @xmath22 towards the backward direction this ratio diminishes more or less uniformly for all @xmath2 . in nlo this ratio depends on the scheme chosen for the photon pdf . in the dis@xmath133 subtraction scheme terms in the photon pdf are shifted to the direct cross section as follows from ( 14)(16 ) in section 2 . this necessarily changes the ratio of the resolved to the direct cross section , not only in the absolute value but also in the dependence on @xmath2 . for all @xmath22 s this ratio is quite different from its lo value . the difference between the @xmath27 and dis@xmath50 scheme is small for @xmath134 moderate for @xmath135 and of the order of @xmath136 at @xmath137 . except at @xmath137 the ratio is always larger in the dis@xmath50 scheme than in the @xmath27 scheme . apart from the fact that the ratios plotted in fig . 3a , b , c , d can not be measured directly , they are scheme dependent in nlo and have very large corrections when going from lo to nlo . in fig . 3a , b , c we have also plotted the ratio for the case that the @xmath2 dependence is neglected in the cross section @xmath138 and taken into account only in the photon propagator ( denoted nlo(@xmath42 ) ) . this approximation gives a result very similar to the lo curve showing that the nlo corrections are more important for @xmath14 . of course at @xmath26 this approximation is equal to the nlo result in the @xmath27 scheme . the results shown so far are for a model with three flavors only and therefore should not be compared to the experimental data except when the contribution from the charm quark is added at least in lo . a more realistic approach is to use the photon pdf s sas1 m of ref . @xcite which are constructed for four flavors . the results for @xmath118 integrated over @xmath139 $ ] as a function of @xmath3 for @xmath140 and @xmath141 gev@xmath8 are presented in fig . 4a , b , c. these curves can be compared with the results in fig . 1a , b , c obtained with the grs photon distribution with @xmath71 . the results for the sum of the resolved and direct contributions change between 10% and 30% in the small @xmath3 region and approximately 50% in the large @xmath3 region . as one can see the larger cross section for @xmath78 results primarily from the direct component . since the direct component is more important for large @xmath3 than for small @xmath3 the increase is stronger in the former region .    of interest are also the rapidity distributions for fixed @xmath3 . these are shown for @xmath127 gev as a function of @xmath22 between @xmath142 choosing @xmath143 and @xmath144 gev@xmath8 in fig . 5a , b , c , d. we show the subtracted direct cross section and the resolved cross section and their sum for the photon pdf as in fig . 4a , b , c. expectedly the resolved component has its maximum shifted to positive @xmath22 s as compared to the direct component . the direct component falls off with increasing @xmath22 quite strongly . this comes from the subtraction of the @xmath145 terms as is apparent when we compare with the unsubtracted direct cross section , labeled `` dir '' in fig.5a , b , c , d. the sum of the resolved and the direct ( subtracted ) cross section is more or less constant for the smaller @xmath2 values and decreases with increasing @xmath22 for @xmath146 and @xmath144 gev@xmath8 .    for large enough @xmath2 we expect the unsubtracted direct cross section to be the correct one . in this region the subtraction term ( 12 ) must approximate the pdf of the photon rather well . we have checked this by a direct comparison . with increasing @xmath2 the full direct cross section ( dir ) approaches the sum = res+dir@xmath125 better and better . but even at @xmath147 gev@xmath8 the two cross sections differ at @xmath135 still by approximately 30 % . however , in the backward direction at @xmath137 we also see a difference . in this region , which corresponds in the photon pdf to the region @xmath148 , where the perturbative component dominates , we do not expect a deviation . this may be caused by our approximation of neglecting the transverse momentum @xmath114 of the virtual photon in the relation ( 22 ) and in the calculation of the resolved cross section .    in order to compare with preliminary data for the dijet cross section presented by the zeus collaboration @xcite we calculated the inclusive dijet cross section @xmath149 as a function of @xmath2 . here @xmath3 is the transverse energy of the measured or trigger jet with rapidity @xmath150 . @xmath151 denotes the rapidity of another jet such that in the three - jet sample these two jets have the largest @xmath3 , i.e.  @xmath152 . the calculational procedure is analogous to the inclusive single - jet cross section and follows closely the work presented in @xcite for @xmath42 . since inclusive two - jet cross sections depend on one more variable they constitute a much more stringent test of qcd predictions than inclusive one - jet cross sections . we could predict distributions in @xmath150 and @xmath151 for fixed @xmath3 or distributions in @xmath3 for various values of the two rapidities @xmath150 and @xmath151 in the same way as in @xcite for @xmath42 . since no such information is expected from the experiment in the near future we calculated only the @xmath3 distribution with the two rapidities integrated over the interval @xmath153 following the constraints in the zeus analysis @xcite . the results for @xmath140 and @xmath141 gev@xmath8 are shown in fig . 6a , b , c , where the full curve is @xmath149 as a function of @xmath3 integrated over @xmath150 and @xmath151 in the specified interval and for @xmath154 . the functional dependence on @xmath3 does not change as a function of @xmath2 , only the absolute value of the cross section decreases with increasing @xmath2 .    in fig . 6a , b , c we show also the cross sections for enriched direct and resolved @xmath155 samples . these cross sections are labeled `` dir '' ( dashed curve ) and `` res '' ( dotted curve ) , respectively . the two cross sections are defined with a cut on the variable @xmath156 where @xmath156 is defined by @xmath157 this variable measures the fraction of the proton energy that goes into the production of the two hardest jets . the `` dir '' curve gives the cross section for @xmath158 . this cut on @xmath156 leads to an enrichment of the direct component of the cross section , since exclusive two - jet events from the direct process have @xmath159 . the curve labeled `` res '' gives the cross section in the complementary region @xmath160 which characterizes the enriched resolved @xmath155 sample . the sum of the dir and res curves is equal to the full cross section @xmath149 with no cut on @xmath156 . it must be emphasized that both cross sections , whether @xmath161 or @xmath162 , contain contributions from the direct and resolved part . actually the resolved contribution in the enriched direct sample ( @xmath161 ) is non - negligible . the results in fig . 6a , b , c are for the grs parton distributions in the @xmath27 scheme . as to be expected , with increasing @xmath2 the full cross section is more and more dominated by the dir component , in particular at the larger @xmath3 . this means that the cross section in @xmath160 decreases stronger with @xmath2 than in the @xmath158 region . this could be studied experimentally by measuring the ratio of the two cross sections as a function of @xmath2 for fixed @xmath3 . this has not been done yet . instead , the zeus collaboration @xcite measured the ratio @xmath163res@xmath164dir , where res and dir are the cross sections as defined above , but integrated over @xmath165 gev . with this integration cut on the transverse momenta of the two hardest jets , the transverse momentum of the third jet may vanish which affords particular constraints on the remnant jets . this is treated as in @xcite . furthermore we replaced the grs photon pdf by the sas1 m photon pdf which is for @xmath78 flavors . with these assumptions we calculated the ratio @xmath166 as a function of @xmath2 up to @xmath167 gev@xmath8 and compared it with the zeus @xcite data in fig .  7 in lo ( dotted curve ) and nlo ( full curve ) . the theoretical nlo curve agrees quite well with the data at @xmath168 gev@xmath8 but not with the measured point at @xmath169 gev@xmath8 and at @xmath170 gev@xmath8 corresponding to the photoproduction case . of course the ratio @xmath166 for @xmath4 is much more precise and lies 30 % higher than the predicted cross section . this disagreement at @xmath4 is to be expected since at this value of @xmath2 the measured inclusive dijet cross section for the enriched resolved @xmath155 sample is larger than the predicted cross section for a small @xmath3 cut @xcite . as in @xcite we attribute this difference between theory and experimental data to additional contributions due to multiple interactions with the proton remnant jet in the resolved cross section not accounted for by our nlo predictions . this underlying event contribution is reduced with increasing @xmath171 and/or smaller cone radii @xmath172 . as it seems , for larger @xmath2 , the underlying event contribution is also reduced . of course this could be studied more directly by measuring rapidity distributions for the enriched resolved @xmath173 sample as was done at @xmath4 in @xcite . we have calculated inclusive single jet and dijet cross sections for photoproduction with virtual photons . the direct and resolved contributions were calculated in next - to - leading order qcd using the phase - space slicing method . they were folded with the unintegrated weizscker - williams approximation and existing lo parametrizations for the virtual photon parton densities . the collinear singularity in the direct photon initial state was integrated analytically up to an invariant mass cut - off @xmath37 . contrary to real photons , this specific singularity is not regulated in the dimensional regularization scheme but by the mass of the virtual photon @xmath2 leading to a logarithmic dependence on @xmath2 . this logarithmic term is absorbed into the virtual photon structure function rendering the latter scheme and scale dependent . the remaining finite contribution is constructed in such a way that the corresponding real photon term is obtained in the limit @xmath42 in the @xmath1 scheme . similarly to the construction of virtual photon structure functions by grs and sas , our calculation of the hard scattering cross section then provides a consistent extension from @xmath42 to small , but non - zero @xmath2 . we presented @xmath37-cut independent results using the snowmass jet definition and hera conditions for distributions in the transverse energy and the rapidity of the observed jet and in the photon virtuality . for the case of very small @xmath2 , we found good numerical agreement with the predictions for real photoproduction . at @xmath147 gev@xmath8 , the unsubtracted direct contribution corresponding to the case of deep inelastic scattering approximates the sum of the subtracted direct and resolved contributions quite well . a small discrepancy remains in the forward region , where the resolved contribution is dominant . as in the case of @xmath42 , the resolved component dominates at low @xmath3 and in the forward region . the nlo effects were demonstrated to be important in the ratio of resolved and direct cross sections as a function of @xmath2 . since the theoretical separation between direct and resolved is artificial , some scheme dependence remains here unless both contributions are defined by an experimental cut on @xmath174 in dijet cross sections . then the corresponding ratio shows significant nlo effects and good agreement with the available zeus data for @xmath175 gev@xmath8 . the disagreement below this value can be attributed to additional contributions coming from , e.g. , multiple scattering between the photon and proton remnants . future investigations on virtual photoproduction will require more data on single inclusive jet production as they exist for @xmath42 and at larger transverse energies . with luminosity permitting , a detailed dijet analysis of an infrared safe cross section such as @xmath176 , where the transverse energies of the two jets are not cut at exactly the same value , will provide a much improved insight into the structure of the virtual photon . furthermore , choosing a @xmath177-cluster - like jet definition with smaller cone radii will reduce both the uncertainties in the jet algorithm and in the underlying event . on the theoretical side , one possible improvement is the correct treatment of the transverse momentum of the incoming photon for larger @xmath2 including a correct transformation from the photonic c.m . frame to the hera laboratory system . for a consistent nlo treatment , the inclusion of nlo parton densities for the photon is necessary . these are , however , needed in a parametrized form and should also be studied in correlation with deep inelastic @xmath178 scattering data .                  for reviews see : a. vogt , proceedings of the workshop on two - photon physics at lep and hera , lund , may 1994 , eds . g. jarlskog and l. jnssen ( lund univ . , 1994 ) p. 141 ; t. sjstrand , j.k . storrow , a. vogt , j. phys . g22 ( 1996 ) 893 . utley on behalf of the zeus collaboration , univ . glasgow report glas - ppe/95 - 03 , talk given at eps conference , brussels , 1995 , hep - ex/9508016 ; c. foudas , zeus collaboration , talk given at dis96 workshop , rome , 1996 .
Abstract:

Answer: we present a next - to - leading order calculation for the virtual photoproduction of one and two jets in @xmath0 collisions . soft and collinear singularities are extracted using the phase space slicing method . the collinear photon initial state singularity depends logarithmically on the mass of the virtual photon and is absorbed into the virtual photon structure function . an @xmath1 factorization scheme is defined similarly to the real photon case . numerical results are presented for hera conditions using the snowmass jet definition for inclusive single jet and dijet cross sections . we study the dependence of these cross sections on the transverse energies and rapidities of the jets . finally , we compare the ratio of the experimentally defined resolved and direct cross sections with recent zeus data as a function of the photon virtuality @xmath2 .