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: ngc  3109 ( ddo  236 ) , a late - type dwarf spiral galaxy , has been classified as sm  iv ( sandage & tammann 1981 ) and is seen almost edge - on with an inclination close to 80@xmath7  ( carignan 1985 ) . since this galaxy has no nucleus it has been assigned the morphological type ir ( van den bergh 1999 ) . it is a southern magellanic dwarf galaxy well resolved into stars and is one of the largest magellanic dwarfs close to the local group at a distance of 1.36@xmath80.10  mpc ( musella et al . 1997 ) with an apparent major - axis diameter of 30  ( 12  kpc ) and minor - axis diameter of 6  ( demers et al . 1985 ) . ngc  3109 has a similar dimension as the lmc . with an absolute magnitude comparable to the smc of @xmath9  mag it is rather underluminous . estimates of the total mass run from 0.6 to @xmath10 . the galaxy is surrounded by a huge hi envelope , quite larger than its optical size ( materne 1980 ; huchtmeier et al . 1980 ; jobin & carignan 1990 ) . jobin & carignan ( 1990 ) derive an hi mass of @xmath11 . the rotation curve requires a large dark matter halo .    on deep exposures ngc  3109 exhibits spiral structure . associations of knots of stars are visible along spiral arms . a globular cluster search revealed ten candidates . bright blue stars are found evenly distributed over the face of the galaxy but particularly along the spiral arms , indicating that star formation is taking place on a galaxy - wide scale ( demers et al . the metallicity of ngc  3109 has been found to be low , similar to the smc ( richer & mccall 1995 ) . minniti et al . ( 1999 ) establish the existence of an extended halo of old and metal - poor stars . ngc  3109 belongs as most luminous member to a subgroup of the local group dwarfs which is relatively isolated ( mateo 1998 ) . another member of this group is the dwarf spheroidal galaxy antila at a distance of only 1.@xmath72 from ngc  3109 . from this apparent separation on the sky the lower limit to the distance between both galaxies is @xmath326  kpc ( whiting et al . 1997 ) .    in this article we derive the x - ray population in the field of ngc  3109 from _ rosat pspc _ observations . we classify a few of the x - ray sources which coincide with the hi extent of ngc  3109 from their x - ray spectral properties . we find two candidate agn which we use to probe the gas content of ngc  3109 . the observations used were carried out with the _ pspc _ detector of the _ rosat _ observatory during two pointed observations in may - june and november  1992 . the satellite , x - ray telescope ( xrt ) and the focal plane detector ( _ pspc _ ) are described in detail in trmper ( 1983 ) and pfeffermann et al . we note that no observations have been performed in the direction of ngc  3109 with the _ hri _ detector of _ rosat_.    lcccc sequence & ra & dec & time interval & expos . + number & & & ( ksec ) + 600174p & 10@xmath123@xmath137.@xmath142&-26@xmath79@xmath1536@xmath16&26-may-92 10:14  10-jun-92 17:27 & 18.5 + 600385p & 10@xmath123@xmath137.@xmath142&-26@xmath79@xmath1536@xmath16&18-nov-92 02:07  18-nov-92 21:27 & 13.4 +    [ tab : pointing ]    the data of the two pointings given in table  [ tab : pointing ] have been retrieved from the public _ rosat _ archive . the data sets from the two observations have been merged to one data set using standard _ exsas _ procedures ( zimmermann et al . 1994 ) . three detection procedures ( local , map , and maximum likelihood ) were applied to the merged pointings using _ exsas _ commands ( zimmermann et al . 1994 ) . in the local and map detection a square shaped detection cell is slid over the image . source counts are determined within the cell , background counts either from an area surrounding the cell ( local background ) or from a ( smoothed ) background image within the same cell . in the maximum likelihood detection the distribution of the detected photons above the background is compared in a maximum likelihood ratio test with the analytical point - spread function . the analysis was performed in the five energy channel ranges soft = ( channel 11 - 41 , 0.1 - 0.4  kev ) , hard = ( channel 52 - 201 , 0.5 - 2.1  kev ) , hard1 = ( channel 52 - 90 , 0.5 - 0.9  kev ) and hard2 = ( channel 91 - 201 , 0.9 - 2.0  kev ) and broad ( 0.1 - 2.4  kev ) . the five source lists were merged to one final source list taking detections at off - axis angles @xmath1750  into account . the maximum likelihood algorithm was used to determine the final source position , the counts in five energy bands and the source extent . a one - dimensional energy and position dependent gaussian distribution was applied in order to obtain the source extent . the source extent ( @xmath18 ) is given as the gaussian @xmath19    @xmath20    hardness ratios @xmath21 and @xmath22 were calculated from the counts in the bands as @xmath23 and @xmath24 . the existence likelihood ratio and the extent likelihood ratio was calculated according to cash ( 1979 ) and cruddace et al . we selected for our final source catalog only detections with an existence likelihood ratio , which is equal to a probability of existence . we give the value for the extent only in case the extent likelihood ratio is . a @xmath25 source error radius was calculated , adding quadratically a systematic error of @xmath26  ( cf . krster 1993 ) . @xmath27    the positional error derived for large off - axis angles may be somewhat underestimated due to the asymmetry of the point - spread - function . but the positional error should not be larger than @xmath31 . we finally screened the catalog of _ rosat _ sources by displaying the positions of these sources on the hard , soft , and broad band _ rosat _ _ pspc _ image . we could confirm 91 of the detected sources . this screened catalog of pointlike and moderately extended sources is given in table  6 . we give in the first column of the catalog for each confirmed source the sequence number and in the second column the source number ( from the catalog of unscreened detections ) . we always refer to the source number in the text .    in fig . [ ps : figros ] we show the _ rosat _ _ pspc _  image of the ngc  3109 field in the hard band ( 0.5  2.0 kev ) . we note that the dwarf spheroidal galaxy antila which is @xmath31.@xmath72  to the south of ngc  3109 is outside of the observed field of view . in fig . [ ps : rossou ] we mark on the _ rosat pspc _ image of the central 20  of ngc  3109 the sources for with accurate hardness ratios @xmath00.2 have been determined . the x - ray colors ( hardness ratios @xmath21 and @xmath22 ) as given in table  6 can be used for a source classification . kahabka et al . ( 1999 , hereafter kpfh99 ) have made a classification of the _ rosat pspc _ x - ray sources in the field of the small magellanic cloud ( smc ) . x - ray binaries have on average harder x - ray spectra than supernova remnants and background agn . the range of hardness ratios @xmath21 and @xmath22 is then different . in addition absorption due to intervening cold gas affects the x - ray colors . kpfh99 have derived the dependence of @xmath21 and @xmath22 on the absorbing column density for a low metallicity galaxy like the smc . the hydrogen column density in the direction of the x - ray source has been derived from the high - resolution 21-cm map of stanimirovic et al . then from the measured x - ray color @xmath22 a classification could be obtained . the strength of this method is that x - ray binaries and agn cover different parts of the diagram for intervening hydrogen columns @xmath28 . agn are located in a band in this diagram with radio loud agn populating the upper regime of the band while radio quiet agn occupy the lower part of the band ( cf . laor et al . 1997 ; brinkmann et al . 1997 ) .    for radio loud agn spectra ( powerlaw photon index @xmath29 ) the dependence of the x - ray colors @xmath21 and @xmath22 on the hydrogen column density @xmath30 ( @xmath31 ) for abundances @xmath30.2 solar ( smc abundances ) has been derived from simulations as    @xmath32    @xmath33    and for radio quiet agn spectra ( powerlaw photon index @xmath34 )    @xmath35    @xmath36    for x - ray binaries with powerlaw photon index @xmath37 spectra the dependence of the x - ray colors @xmath21 and @xmath22 on the hydrogen column density is    @xmath38    @xmath39    equations  3 to 8 can be solved for the hydrogen column density @xmath30 . in combining the equations for @xmath21 and @xmath22 one can derive analytical solutions of x - ray binary and agn tracks in the @xmath21  @xmath22 plane .    for radio loud agn ( powerlaw photon index @xmath29 ) one derives the track    @xmath40    and for radio quiet agn ( powerlaw photon index @xmath34 ) the track    @xmath41    for x - ray binaries ( powerlaw photon index @xmath37 ) we derive the track    @xmath42    these tracks allow a classification of a _ rosat _ source as an agn ( or an x - ray binary ) without the requirement that the intervening hydrogen column density is known . it even is possible to constrain in the @xmath21  @xmath22 plane the intervening hydrogen column density for the source . we note that x - ray binaries are not necessarily seen through the total gas column of the galaxy disk . with a galactic foreground column of @xmath43 towards ngc  3109 we derive with equation  11 for x - ray binaries a lower bound on @xmath22 of @xmath44 . the tracks given above have been derived for a metallicity @xmath30.2 solar which is consistent with the metallicity derived for hii regions in ngc  3109 ( cf . minniti et al . 1999 ) . higher metallicity tracks deviate somewhat from these tracks especially if one approaches the high column density regime @xmath45 . huchtmeier et al . ( 1980 , hereafter hsm80 ) measured the hi distribution across ngc  3109 with the effelsberg radio telescope and a beam size of 9 . we apply the procedure used for the x - ray detections in the field of the smc to the x - ray detections in the field of ngc  3109 . first we make use of the hi map of hsm80 which has an extent of 60  @xmath1 30 . we convert the hi intensity into hydrogen column densities with the equation given in dickey & lockman ( 1990 ) , cf . kahabka ( 1999 ) and we derive a peak of @xmath46 . hsm80 have found that ngc  3109 has a large extent in hi , with an extension ( distortion ) of the hi in the sw . with our x - ray catalog we cover the whole extent of the hi of ngc  3109 . if we restrict the analysis to the inner 20  of this field and to sources with well constrained hardness ratios @xmath00.2 then we can classify 7 sources as agn . a much higher resolution hi map of ngc  3109 with a beam size of 40  has been derived by jobin & carignan ( 1990 , hereafter jc90 ) with the vla . we now make use of this high - resolution hi image which we take from plate  67 of jc90 . the hi distribution of ngc  3109 has an extent of 40  @xmath1 12 . the peak hydrogen column density is @xmath47 , and the lowest column density is @xmath48 .    in fig . [ ps : hijc ] we show the positions of the cataloged _ rosat _ _ pspc _ sources overlaid on the gray scale hi image of ngc  3109 and taken from jc90 . we find 26 _ rosat _ _ pspc _ sources within the hi contours of jc90 . if we restrict the analysis to the inner 20  of this field and to sources with well constrained hardness ratios @xmath00.2 , then we can classify 3 sources ( number 36 , 41 and 63 ) as agn and 2 sources ( number 53 and 73 ) as x - ray binaries . source 53 may be seen through higher gas columns while source 73 is seen through lower gas columns ( see fig . [ ps : hihr2jc ] and equation  11 ) . one source ( number 82 ) can be either class . the two sources 52 and 92 can not be classified as agn or x - ray binaries . source 92 may be a foreground object ( cf . [ ps : hihr2jc ] , fig . [ ps : hr1hr2 ] and table  [ tab : classjc ] ) . source 52 has similar hardness ratios as the lmc snr 0548 - 70.4 ( haberl & pietsch 1999 ) . it could be a ( young ) snr in ngc  3109 . ll source no . & class + 36 , 41 , 63 & a + 53 , 73 & b + 82 & ab + 28 , 68 , 92 & f +    [ tab : classjc ]    ngc  3109 has a mass smaller than , or comparable to , the lmc and 2 x - ray binaries with luminosities above a few times @xmath49 would be in agreement with extrapolations from _ rosat _ findings for the lmc ( cf . haberl & pietsch 1999 ) . there are 14 _ rosat pspc _ sources which are projected onto ngc  3109 intrinsic hydrogen columns of @xmath50 ( cf . [ ps : hijc ] ) . the largest column is derived for source 86 ( @xmath51 ) . this source could be associated with a spiral arm or an hii region of ngc  3109 . two further sources , 53 and 59 , are close to another region of large column density and source 59 is also close to the optical center of ngc  3109 . we further test the correctness of the classification by constructing the distribution of the number of detected x - ray sources @xmath52 with fluxes in excess of a given flux @xmath53 , the @xmath54 ( cf . hasinger et al . we restrict the analysis to sources within the central 20   of ngc  3109 which have been classified as agn or x - ray binaries . we correct the x - ray fluxes for the intervening hydrogen columns by using for the galactic contribution a value of @xmath2 ( dickey & lockman 1990 ) and for the contribution due to ngc  3109 the hi model of jc90 . we find that the observed @xmath54 is in excess of the @xmath54 of the soft extragalactic x - ray background ( hasinger et al . this can be accounted for if an additional less steep component ( e.g. due to x - ray binaries ) is added . if we only construct the @xmath54 of the candidate agn , then the @xmath54 of the soft extragalactic x - ray background is reproduced . this would mean that no significant additional hydrogen is needed to explain the observed sources as background agn and x - ray binaries although only a few candidate agn are within the hi extent of ngc  3109 . the two brightest candidate agn ( source 36 and 63 ) could give an excess in the @xmath54 . we can not exclude that they are at least in part related to ngc  3109 . we generated finding charts in the b - band _ cosmos _ blue plates and in the i - band ( _ ntt emmi _ , cf . [ ps : charts ] ) to search for optical counterparts of the _ rosat _ sources . the b - band plates were used to classify the source using the flux ratio @xmath55 while the i - band image was used to determine the source morphologies from an analysis of their point spread function ( see below ) . we produced finding charts of the _ rosat _ sources given in table  6 using the _ cosmos _ blue plates . we especially investigated the classified _ rosat _ sources given in table  [ tab : classjc ] . we searched for the optical counterparts in the 90% confidence circle of the _ rosat _ source . we list the parameters of the matches in table  [ tab : matches ] . for the two sources 28 and 92 we find bright ( b=11.7  mag and b=12.1  mag ) optical counterparts in the 90% confidence circle of these _ rosat _ sources . this indicates a galactic foreground nature of these sources . for the sources 36 , 41 , 52 , 68 and 73 we find weak ( b=20.4  mag to b=22.3  mag ) unresolved optical counterparts . for source 73 there is also a galaxy ( b=22.7  mag ) in the error circle . for source 53 we can not find an optical counterpart as it is in the crowded field of ngc  3109 and for 63 no optical counterpart exists within the _ rosat _ 90% confidence circle . rccc source & b - mag & @xmath56 & comment + no . & & & + 28 & 11.7 & 3.8 & foreground star + 36 & 20.4 & 0.33 & unresolved , agn + 41 & 21.1 & 0.35 & unresolved , agn + 52 & 20.4 & 0.01 & unresolved + 63 & + 68 & 22.3 & 0.47 & unresolved + & 24.4 & 1.3 & faint + 73 & 21.7 & 0.56 & unresolved + & 22.7 & 0.96 & galaxy + & 25.9 & 2.2 & faint + 92 & 21.2 & 0.48 & unresolved + & 12.1 & 3.2 & foreground star@xmath57 +    note ( a ) : 17  distant from the _ rosat _ source [ tab : matches ]    in addition we calculated the flux ratio @xmath58 with the equation given in haberl & pietsch ( 1999 ) . we make use of the b - magnitude given in table  [ tab : matches ] and the _ pspc _ countrate given in table  6 . we use the abbreviation @xmath59 and obtain the flux ratio from the equation :    @xmath60    using the scheme given in table  3 of haberl & pietsch ( 1999 ) we classify the sources 28 and 92 as candidate foreground stars and the sources 36 and 41 as candidate agn ( cf . table  [ tab : matches ] ) . for the first run of optical matching of x - ray sources in the central part of ngc  3109 we obtained three optical i - band images ( 60.e-0818 ) from the eso archive . all images were taken on 1998 february 2 with the red arm of the eso multi - mode instrument ( _ emmi _ ) at the new technology telescope ( _ ntt _ ) with a total exposure time of 2700 s. the images were taken with the tektronix 2048@xmath12048 pix@xmath61 chip with a scale of 0.268/pix . the resulting field of view covered by the stacked images is 8.8  @xmath1 8.8  and the final image is centered on the galaxy nucleus . the seeing was measured to be @xmath62  throughout all three exposures . the data are used to identify optical couterparts within the region of highest hi column density , which resides at the center of ngc  3109 ( jobin & carignan 1990 ) and is best covered by the field of view of this particular dataset , and to support or reject the classification provided by the hardness - ratio estimate ( see sect .  5 ) . further photometric studies in optical passbands of _ rosat _ sources for most objects in the catalog in table  6 with accurate positions will be given in a subsequent paper .    in order to obtain bona - fide optical counterparts we transformed the pixel coordinates of the i - band to new equatorial coordinates . we used positions of 43 stars out of the usno astrometric catalog ( monet et al . 1998 ) which was obtained from centre de donnes de astronomiques de strasbourg ( cds ) . using the task _ ccmap _ within the iraf environment ( tody et al . 1993 ) we find the mean rms uncertainty for all given optical coordinates being @xmath63ra@xmath64  and @xmath63dec@xmath65 . applying the plate solution to coordinates provided by the _ rosat pspc _ catalog ( see table  6 ) we find 7 optical matches within the field of view . the finding charts for all these regions are given in fig . [ ps : charts ] with the appropriate scaling of the 90%-confidence circle as given in table  6 . the five brightest objects are encircled and labeled according to their luminosity ( label  1 marks the brightest object ) . subsequently , we perform an analysis of the shape parameters of the optical point spread function ( psf ) for the 5 brightest objects within each of the 7 object s finding circles . we used the source - extraction software sextractor v2.1.4 ( bertin & arnouts 1996 ) for determination of ellipticity @xmath66 , full width at half maximum ( fwhm ) , and the star - galaxy - classification parameter cl ( using a neural - network algorithm which was extensively trained , see bertin & arnouts for details ) . we compare our findings to all remaining optical objects most of which are assumed to be point - sources . the code yields @xmath37500 detections for which we plot the three psf - shape parameters ( @xmath66 , fwhm , and cl ) as a function of instrumental i - band magnitude . [ ps : shapepars ] shows the parameters of all detections together with the parameters for the brightest source within the finding circle ( see table  4 ) . note that the brightest object is _ not _ necessarily the optical counterpart of the _ rosat _ source . crrrrr label@xmath67 & i@xmath68 & @xmath63i@xmath68 & fwhm@xmath69 & ell@xmath70 & class@xmath71 +   + 1 & @xmath72 & 0.01 & 0.76 & 0.023 & 0.98 + 2 & 0.01 & 0.03 & 0.90 & 0.076 & 0.96 + 3 & 0.80 & 0.05 & 0.81 & 0.024 & 0.84 + 4 & 0.83 & 0.06 & 0.93 & 0.089 & 0.77 + 5 & 0.93 & 0.06 & 0.79 & 0.065 & 0.66 +   + 1 & @xmath73 & 0.01 & 4.52 & 0.323 & 0.15 + 2 & @xmath74 & 0.01 & 0.87 & 0.337 & 0.96 + 3 & @xmath75 & 0.01 & 1.36 & 0.229 & 0.97 + 4 & @xmath76 & 0.02 & 7.17 & 0.373 & 0.02 + 5 & @xmath76 & 0.02 & 5.69 & 0.622 & 0.01 +   + 1 & @xmath77 & 0.01 & 1.87 & 0.244 & 0.02 + 2 & 0.05 & 0.03 & 2.68 & 0.448 & 0.17 + 3 & 0.39 & 0.04 & 0.75 & 0.065 & 0.88 + 4 & 0.60 & 0.05 & 1.42 & 0.301 & 0.81 + 5 & 0.79 & 0.05 & 0.87 & 0.231 & 0.73 +   +   +   + 1 & @xmath78 & 0.01 & 0.75 & 0.037 & 0.98 + 2 & @xmath79 & 0.02 & 0.95 & 0.066 & 0.81 + 3 & 0.01 & 0.03 & 1.41 & 0.288 & 0.81 + 4 & 0.20 & 0.03 & 1.16 & 0.102 & 0.78 + 5 & 0.65 & 0.05 & 1.00 & 0.015 & 0.84 +   + 1 & @xmath80 & 0.02 & 0.79 & 0.049 & 0.96 + 2 & @xmath81 & 0.02 & 0.73 & 0.059 & 0.96 + 3 & 0.21 & 0.03 & 0.81 & 0.088 & 0.89 + 4 & 0.47 & 0.04 & 0.72 & 0.076 & 0.86 + 5 & 0.97 & 0.06 & 0.77 & 0.075 & 0.67 +   + 1 & @xmath82 & 0.01 & 0.78 & 0.007 & 0.98 + 2 & @xmath83 & 0.02 & 0.80 & 0.033 & 0.96 + 3 & 0.80 & 0.05 & 0.82 & 0.011 & 0.62 + 4 & 1.03 & 0.07 & 1.85 & 0.414 & 0.67 + 5 & 1.77 & 0.13 & 0.72 & 0.245 & 0.49 +    \(a ) numbering in agreement with labels in fig . [ ps : charts ] ( b ) instrumental i magnitude and the photometric error in mag . ( c ) full width at half maximum in arc seconds . ( d ) ellipticity , i.e. 1a / b . ( e ) star - galaxy - classification parameter of sextractor ; objects with cl@xmath84 are point sources while cl@xmath85 classifies objects as extended . [ tab : shapepars ]    the brightest sources 52.1 , 62.1 and 65.1 ( the first number being the source number and the decimal the source label ) in fig . [ ps : charts ] are at least one magnitude brighter than the next fainter object within the finding circle ( see table  [ tab : shapepars ] ) . object 63.1 has about the same magnitude as object 63.2 . therefore no clear preference can be given just from considering the i - band flux . the psf - shape parameters of these 4 objects agree with those found for stars ( see fig . [ ps : shapepars ] ) . the cl classification of sextractor puts all 4 sources in the point - source regime . yet , the optical findings are prone to mis - classification since not every brightest object is most central within the finding circle . although they are not the brightest , the sources 62.5 and 63.2 are the most likely candidate optical counterparts ( infering from their most central position only ) . these were classified as slightly extended and slightly elliptical , respectively .    as can be seen from fig . [ ps : charts ] there are 2 finding charts ( i.e. charts 62 and 65 ) from the edge of the ccd image . thus , we can not exclude brighter objects within the area not covered by the exposure . chart 59 is too crowded to give a reliable identification only on the basis of any photometric psf - shape parameter . however , we note that the center of the _ rosat _ error circle coincides with the galaxy s center . thus , it is very likely that the _ rosat _ source is connected with a x - ray source which is located near the center of the galaxy . in sect.5 we have used the hardness ratios @xmath21 and @xmath22 to classify 3 _ rosat _ sources ( with numbers 36 , 41 and 63 ) and projected onto the hi extent of ngc  3109 as candidate agn . we have searched in the optical b - band _ cosmos _ finding charts of these sources for optical counterparts . in addition we have generated deep i - band finding charts for source 63 and other _ rosat _ sources and we have looked for faint unresolved and extended objects within the positional error circles of these _ rosat _ sources ( cf . [ ps : charts ] ) . we have found several optical candidates in the positional error circle of each _ rosat _ source . therefore optical spectroscopy of all candidates is required to firmly identify the optical counterpart of these _ rosat _ sources . if the background nature of candidate agn `` close to '' ngc  3109 is established then these agn can be used to probe the gas column density of ngc  3109 in the direction of these agn ( taking the additional galactic gas column into account ) . for two sufficiently x - ray bright candidate agn rx  j1003.2 - 2607 ( source 63 ) and rx  j1002.5 - 2602 ( source 36 ) we perform x - ray spectral fitting with the _ exsas _ spectral analysis task ( zimmermann et al . we assume a galactic foreground column density of @xmath86 and determine the additional hydrogen column density ( due to ngc  3109 ) assuming reduced metallicities ( @xmath30.2 solar ) . for both sources we find that the hydrogen column deduced from the x - ray spectral fit is consistent with the hydrogen column deduced from the 21-cm survey of hsm80 and jc90 ( cf . table  [ tab : specpar ] and fig . [ ps : xrayfit ] ) . cccccc @xmath87&flux ( @xmath88 & + & @xmath89 ) & + & & & + & & ( a ) & ( b ) & ( c ) & ( d ) +   + -2.2@xmath90 & 7.9@xmath91 & 11@xmath92 & 11@xmath93 & 10 & 8@xmath82 +   + -2.4@xmath94 & 3.1@xmath95 & 4.1@xmath96 & 4.1@xmath97 & 4 & 1 +    \(a ) no constraint on @xmath87 ; ( b ) constraint on powerlaw photon index : @xmath98 ; ( c ) hi model of hsm80 ; ( d ) hi model of jc90 . [ tab : specpar ]    this result is also consistent with the location of these two sources in the @xmath21  @xmath22 plane ( cf . [ ps : hr1hr2 ] ) . we do not include in the spectral analysis source 53 which has been classified as a candidate x - ray binary in sect.5 . assuming that both agn have canonical powerlaw photon indices @xmath87 = @xmath99(1.75 to 2.25 ) , cf . laor et al . ( 1997 ) , we further tighten the constraints on the hydrogen column density . for rx  j1003.2 - 2607 we find from the x - ray spectral fit @xmath100 and from the 21-map of jc90 @xmath101 . for rx  j1002.5 - 2602 we constrain the hydrogen column density from the x - ray spectral fit to @xmath102 and from the 21-map of jc90 to @xmath103 . both values are in agreement within the uncertainties . the column densities of absorbing material from the x - ray data are slightly larger than those from the 21-cm line , but the uncertainty of the x - ray spectral fit leaves room for some additional hydrogen ( e.g. molecular hydrogen ) in the line of sight of the agn . we assume that the photoionisation cross section is @xmath104 times larger for molecular hydrogen than for atomic hydrogen ( cf . cruddace et al . 1974 ; yan et al . 1998 ) . then we determine the column density due to molecular hydrogen @xmath105 from the total hydrogen column density @xmath106 derived from the x - ray spectral fit and the atomic hydrogen column density @xmath107 derived from the 21-cm observations as    @xmath108    we so constrain the gas column due to molecular hydrogen in ngc  3109 from the x - ray spectral fit of rx  j1003.2 - 2607 to @xmath109 .    towards rx j1003.2 - 2607 the amount of @xmath110 is @xmath3@xmath111 . along this line of sight the molecular mass fraction , @xmath112 . this means @xmath113% of the mass of the total gas is in molecular form . this result can be compared with co observations of ngc  3109 , from which the mass of molecular hydrogen has been determined to @xmath114@xmath115 ( rowan - robinson , philips & white 1980 ) . with a hi mass of @xmath116 a molecular mass fraction of @xmath11410% is obtained . extraplanar absorbing clouds have been found in a high - resolution survey of a sample of 12 edge - on galaxies by howk & savage ( 1999 ) at distances of @xmath117 to @xmath118  kpc from the galaxy plane . for an inclination of 80@xmath7 a radius of the galaxy disk of 12  kpc and a distance to the galaxy of 1.4  mpc similar clouds would be projected in ngc  3109 1  to 4   from the galaxy plane . from these clouds dust absorption has been observed but the dust should coexist with gas in the molecular and atomic phase . similar clouds may exist in ngc  3109 and the candidate agn rx  j1003.2 - 2607 would be seen through gas which is at a height @xmath1191  kpc above the galaxy plane while the candidate agn rx  j1002.5 - 2602 is seen through gas with @xmath1192.4  kpc . from _ rosat _ _ pspc _ observations of the dwarf galaxy ngc  3109 we derive 10 x - ray sources which are contained within the hi extent of this galaxy and which have accurate hardness ratios . we classify 2 of these sources as foreground stars , 2 as candidate x - ray binaries and 3 as candidate background agn . from x - ray spectral fitting we derive for 2 of the agn total hydrogen column densities which we compare with the hi column densities inferred from 21-cm line measurements . we estimate that the molecular mass fraction of the gas is not larger than @xmath360% . upcoming spectroscopy of candidate optical counterparts will help to constrain the nature of these objects and help to understand the hydrogen content of ngc  3109 . furthermore , spectroscopy of globular clusters around ngc  3109 will help to trace the enrichment history ( i.e. the major star formation episodes ) and the kinematics ( i.e. total mass ) of ngc  3109 . the  project is supported by the max - planck - gesellschaft and the bundesministerium fr forschung und technologie ( bmft ) . this research made use of the _ cosmos _ digitized optical survey of the southern sky , operated by the royal observatory edinburgh and the naval research laboratory , with support from nasa . this research has made use of the simbad data base operated at cds , strasbourg , france . we thank k.s . de boer and u. klein for critically reading the manuscript . we thank an anonymous referee for useful comments . bertin e , arnouts , s. , 1996 , a&as 117 , 393 brinkmann w. , yuan w. , & siebert j. , 1997 , a&a 319 , 413 carignan c. , 1985 , apj 299 , 59 cash w. , 1979 , apj 228 , 939 cruddace r. , paresce f. , bowyer s. , & lampton m. , 1974 , apj 187 , 497 cruddace r.g . , hasinger g. , & schmitt j.h.m.m . , 1988 , in : astronomy from large databases , eds . murtagh f. , heck a. , p. 177 demers s. , irwin m.j . , & kunkel w.e . , 1985 , aj 90 , 1967 dickey j.m . , & lockman f.j . , 1990 , ara&a 28 , 215 haberl f. , & pietsch w. , 1999 , a&as 139 , 277 hasinger g. , burg r. , giacconi r. , et al . , 1993 , a&a 275 , 1 howk j.c . , & savage b.d . , 1999 , aj 117 , 2077 huchtmeier w.k . , seiradakis j.h . , & materne j. , 1980 , a&a 91 , 341 [ hsm80 ] jobin m. , & carignan c. , 1990 , aj 100 , 648 [ jc90 ] kahabka p. , 1999 , a&a 344 , 459 kahabka p. , pietsch w. , filipovic m.d . , & haberl f. , 1999 , a&as 136 , 81 [ kpfh99 ] krster m. , 1993 , rosat status report no . 67 laor a. , fiore f. , elvis m. , et al . , 1997 , apj 477 , 93 mateo m. , 1998 , ara&a 36 , 1998 materne j. , 1990 , _ eso / esa workshop on dwarf galaxies _ , ed . m. taringhi & k. kjar , eso garching , p.67 minniti d. , zulstra a.a . , alonso m.v . , 1999 , aj 117 , 881 monet d. , bird a. , canzian b. , dahn c. , guetter h. , harris h. , henden a. , levine s. , luginbuhl c. , monet a.k.b . , rhodes a. , riepe b. , sell s. , stone r. , vrba f. , walker r. , 1998 , u.s . naval observatory flagstaff station ( usnofs ) and universities space research association ( usra ) stationed at usnofs . musella i. , piotto g. , & capaccioli m. , 1997 , aj  114 , 976 pfeffermann e. , briel u.g . , hippmann h. , et al . , 1986 , proc . spie 733 , 519 richer m.g . , & mccall m.l . , 1995 , apj  445 , 642 rowan - robinson m. , philips t.g . , & white g. , 1980 , a&a 82 , 381 sandage a. , & tammann g.a . , 1981 , _ a revised shapley - ames catalogue of bright galaxies _ , carnegie institurion of washington , washington d.c . stanimirovic s. , staveley - smith l. , dickey j.m . , et al . , 1999 mnras 302 , 417 tody , d. 1993 , `` iraf in the nineties '' in astronomical data analysis software and systems ii , a.s.p . conference ser . , vol 52 , eds . hanisch , r.j.v . brissenden , & j. barnes , 173 . trmper j. , 1983 , adv . space res . 2 , no . 4 , 241 van den bergh s. , 1999 , apj 517 , l97 vennes s. , 1999 , apj 525 , 995 voges w. , aschenbach b. , boller t. , et al . , 1999 , a&a 349 , 389 whiting a.b . , irwin m.j . , & hau g.k.t . , 1997 , aj 114 , 996 yan m. , sadbeghpour h.r . , & dalgarno a. , 1998 , apj 496 , 1044 zimmermann h.u . , becker w. , belloni t. , et al . , 1994 , mpe report 257
Abstract:

Answer: we established the catalog of x - ray point sources in the field of the magellanic - type spiral galaxy ngc  3109 ( ddo  236 ) from two _ rosat pspc _ observations . of the 91 x - ray sources 26 are contained within the hi extent of ngc  3109 as derived by jobin & carignan ( 1990 ) with the vla . for 10 of these we can determine accurate hardness ratios @xmath00.2 . we find 3 candidate agn , 2 candidate x - ray binaries and one source which may belong to either class . we also find 2 candidate foreground stars . in a field of 8  @xmath1  8 observed with the ntt in the i - band and centered on the nucleus of ngc  3109 we determine candidates for optical counterparts in the x - ray error circle of 7 _ rosat pspc _ sources . we apply a spectral fit to the _ rosat _ spectrum of the x - ray brightest absorbed candidate agn behind ngc  3109 , rx  j1003.2 - 2607 . assuming a galactic foreground hydrogen column density of @xmath2 we derive from the x - ray spectral fit , assuming reduced metallicities ( @xmath30.2 solar ) , a hydrogen column density due to ngc  3109 of @xmath4 . this value is slightly larger than the hydrogen column density derived from the 21-cm observations of @xmath5 . we estimate that the molecular mass fraction of the gas is not larger than @xmath6% .


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 quantum interface between light and an atomic ensemble has attracted much attention during the past decade @xcite . the reason is rooted in the wide range of possible applications that rely on the coherent coupling between the two systems : implementation of quantum information processing protocols @xcite , slow - light and electromagnetically induced transparency @xcite as well as the creation of deterministic and manipulable photon sources @xcite . in order to create these photon sources , one must be able to create entangled atomic many - body states and map them efficiently into the desired photonic states @xcite . ultracold atoms represent the ideal tool for creation these atomic states under very clean and well - defined conditions . this is due to the advanced techniques that are nowadays available for the trapping of ultracold atoms and for tailoring the interactions between them @xcite . highly excited ( so - called rydberg ) states @xcite have proven to be especially useful in this context due to the strong state - dependent interaction between them , allowing the entanglement of atoms separated by several micrometers @xcite .    in a recent work a laser - driven gas of rydberg atoms trapped on a one - dimensional ring lattice has been proposed as a resource to create collective excitations @xcite . it was shown that these excitations could be converted into few photons following a mapping scheme proposed in refs . @xcite . as a result non - classical states of light such as single - photon sources and entangled pairs of photons could be created @xcite whose properties where determined by the interplay of the ring geometry and the nature of the delocalized atomic excitations . in this work , we extend this study to a gas of atoms trapped in two - dimensional medium - scale square and triangular lattices . we study thoroughly the properties of collective many - body states in these two lattice types and outline methods of selective excitation . we find that one can create single - photon sources in which the photon is emitted in a superposition of a few well - focused beams . the direction of these beams is mainly determined by the geometry and spacing of the lattice . the paper is structured as follows . in section [ sec : system ] we introduce the system we have in mind and present a derivation of the hamiltonian that governs its dynamics . in section [ sec : diagonal ] we discuss the diagonalization of this hamiltonian and explain the approximations considered to do so . the symmetry properties of the many - particle eigenstates are analyzed in section [ sec : excitation ] , together with a proposal of how to excite experimentally these states . we explain how to map these atomic collective excitations into non - classical states of light and discuss their properties in section [ sec : photon ] . finally , in section [ sec : perturbations ] we discuss the validity of the approximations made and take into account disorder caused by quantum uncertainty and/or temperature . we conclude with a summary and outlook in section [ sec : conclusion ] . in the setup we have in mind a gas of atoms is trapped in a regular two - dimensional lattice with @xmath0 sites and a lattice spacing @xmath1 of the order of a few micrometers @xcite ( see fig . [ fig : lattice_and_levels]a ) . we consider the limit of a deep lattice where the vibrational states within each site are well - approximated by the eigenstates of a harmonic oscillator and tunneling between the sites is absent . the lattice shall be uniformly occupied with one atom per site ( mott - insulator state ) and we assume zero temperature such that each atom populates the ground state of its potential well whose spatial width @xmath2 is much smaller than the interparticle separation @xmath1 . the internal ( electronic ) degree of freedom of each atom is approximated by a two - level system . the ground state @xmath3 is coupled to a rydberg ns - state @xmath4 by means of a laser with rabi frequency @xmath5 and detuning @xmath6 as shown in fig . [ fig : lattice_and_levels]b . for the moment , we do not focus on the two auxiliary levels @xmath7 and @xmath8 also shown in this figure . they will be used for the excitation of the many - body states and the mapping of these states into light , that will be explained thoroughly in sections [ sec : excitation ] and [ sec : photon ] , respectively . within the two - level approximation we can identify each atom as a spin-@xmath9 degree of freedom where @xmath10 and @xmath11 @xcite . we can thus employ the pauli matrices ( @xmath12 , @xmath13 and @xmath14 ) to formulate the hamiltonian governing the excitation dynamics of our system . using this the hamiltonian describing the laser excitation process reads in the rotating wave approximation @xmath15,\ ] ] where @xmath16 is the rydberg number operator of the @xmath17-th site . when the atoms are in the rydberg state , they interact strongly via van - der - waals interaction @xmath18 , where @xmath19 is the separation between the atoms and @xmath20 is the van - der - waals coefficient @xcite . the interaction hamiltonian thus becomes @xmath21 with the coefficients @xmath22 where @xmath23 denotes the position of the atom trapped in the @xmath17-th site . adding the interaction with the laser , the complete hamiltonian of the system reads @xmath24.\ ] ]    we assume the external dynamics of the atoms to be frozen . this assumption is well justified because the internal dynamics take place on a much shorter time scale ( typically of the order of a hundred nanoseconds versus milliseconds ) . we will also assume for the moment that the atoms are infinitely localized at the center of the traps , i.e. @xmath25 . this very idealized situation can only be achieved approximately in experiment . we will investigate the effect of a finite width of the wave packets in section [ sec : perturbations ] .     which is considered to be much smaller than the lattice spacing @xmath1 . the largest interaction energy in the system is the one between nearest neighbors @xmath26 . * b * : internal level structure of each of the @xmath0 atoms on the lattice . the two ground states @xmath3 and @xmath7 are coupled to the rydberg state @xmath4 by means of two classical lasers with rabi frequencies @xmath5 and @xmath27 , respectively . a third laser @xmath28 drives the transition @xmath29 . finally , @xmath3 , @xmath7 and an auxiliary state @xmath8 form a lambda scheme where @xmath7 is coupled off - resonantly to @xmath8 and photons are emitted in the transition @xmath30 . ] in this section , we will show how to obtain the eigenvalues and eigenstates of the hamiltonian ( [ eqn : working_hamiltonian ] ) and discuss the approximations used to obtain them . throughout this work , we will focus on the regime where the laser driving @xmath5 is the largest energy scale of the system . in addition , the detuning of the laser @xmath31 is taken to be much smaller than both the rabi frequency ( laser driving ) and the interaction strength such that @xmath32 . the hamiltonian ( [ eqn : working_hamiltonian ] ) contains spin operators that obey anti - commutation and commutation relations when they belong to the same and different sites , respectively . therefore , the underlying algebra is neither bosonic nor fermionic . in order to be able to solve the system at least approximately the spin operators will be expressed in terms of operators that obey a purely bosonic algebra . this is achieved by the holstein - primakoff transformation @xcite . we will make use of a slightly unconventional formulation of this transformation . as @xmath32 , the first term of hamiltonian ( [ eqn : working_hamiltonian ] ) is the dominant one and it is thus convenient to make it diagonal . for this reason , instead of the usual form of the holstein - primakoff transformation , we will use a rotated version of it such that @xmath33 . this is achieved by means of the unitary operator @xmath34 $ ] , after which the holstein - primakoff transformation becomes @xmath35 where @xmath36 ( @xmath37 ) are bosonic operators that create ( annihilate ) non - interacting bosonic excitations within the system at the @xmath17-th site . considering that the eigenstates of the single - atom operator @xmath12 are @xmath38,\ ] ] the equivalent bosonic states are @xmath39 and @xmath40 , with zero and one bosonic excitation on site @xmath17 , respectively .    by making the substitution ( [ eqn : holstein - primakoff ] ) , the hamiltonian ( [ eqn : working_hamiltonian ] ) takes on a rather complicated form that contains square roots of operators , which are difficult to treat in practical calculations . in order to overcome this difficulty , we make an approximation : our study will be focused on many - body states that carry only a few excitations , i.e. whose total number of bosonic excitations @xmath41 is much smaller than the number of sites @xmath0 . in the subspace spanned by these states - denoted by @xmath42 - one has @xmath43 for all @xmath17 . dropping the notation @xmath42 in the following , one can make a taylor expansion of the square root operator such that , after the normal ordering of the operators , one obtains @xmath44 we now substitute the operators of eqs . ( [ eqn : holstein - primakoff ] ) into hamiltonian ( [ eqn : working_hamiltonian ] ) and consistently with the condition @xmath43 keep only the first two terms of the expansion of the square root . furthermore , we neglect all terms of the hamiltonian that are of higher order than quadratic in @xmath36 . the hamiltonian then assumes a quadratic form @xmath45 where the energy offset @xmath46 is given by @xmath47      a closer inspection of the hamiltonian ( [ eqn : hamiltonian_hp ] ) reveals that it is composed of three distinct parts ( @xmath48 ) which effectuate qualitatively different couplings between states with a fixed number of bosonic excitations . the first part is formed by those terms which , applied to a given state , do not change its number of excitations , i.e. @xmath49 as they contain the same number of creation and annihilation operators . second , there are a number of terms that create or annihilate a single bosonic excitation , @xmath50 finally , there is a third class of terms , which change the total number of bosonic excitations in the system by two : @xmath51     of bosonic excitations separated by an energy gap of approximately @xmath52 . the arrows represent schematically the transitions driven by the different parts of the hamiltonian ( [ eqn : hamiltonian_hp]).,width=226 ]    due to the strong laser driving ( @xmath32 ) , the first term of the hamiltonian @xmath53 which is proportional to @xmath5 is clearly dominant and thus determines the coarse structure of the spectrum of the system . this coarse structure is formed by manifolds of quasi - degenerate states separated by energy gaps of approximately @xmath52 . as it is shown schematically in fig . [ fig : spectrum ] , each manifold is formed by states that contain the same total number @xmath54 of bosonic excitations . taking this into account , let us now analyze the rest of the terms of the hamiltonian one by one . the second term in @xmath53 does not change the total number of excitations of the state . its action produces couplings between states inside the same energy manifold , whose strength is proportional to @xmath55 . the hamiltonians @xmath56 and @xmath57 create or destroy one and two excitations every time , respectively . since the states involved in these transitions are separated energetically by approximately @xmath52 and @xmath58 , respectively , they are strongly suppressed . in particular , the corresponding transition rates can be estimated by second order perturbation theory to be proportional to @xmath59 and @xmath60 , which in our parameter regime ( @xmath61 ) are very small quantities . as a consequence it can be affirmed that , to a very good degree of approximation , the hamiltonian ( [ eqn : ham_zero ] ) that drives the intra - manifold dynamics also governs the dynamics of the _ entire system_. the validity of this approximation will be assessed in section [ sec : perturbations ] . we will now diagonalize the hamiltonian @xmath53 [ eq . ( [ eqn : ham_zero ] ) ] . since it is quadratic this is achieved by means of the diagonalization of the matrix @xmath62 , that essentially contains the interaction energies between each pair of atoms . this is achieved by a unitary transformation which reads @xmath63 where @xmath64 represents the @xmath17-th entry of the @xmath65-th eigenvector of @xmath62 and @xmath66 is its corresponding eigenvalue . substituting the last expression into ( [ eqn : ham_zero ] ) , the diagonalized hamiltonian becomes @xmath67 with the operators @xmath68 and the eigenvalues @xmath69 for @xmath70 . hence , the _ collective excitations _ of the system are formed by a superposition of all the bosonic excitations @xmath71 , with the coefficients of the superposition given by the eigenvectors of @xmath62 . we can now write down the eigenstates and corresponding energies of the system . the ground state is @xmath72 and its energy is @xmath46 , given by eq . ( [ eqn : ground ] ) . the first excited manifold is spanned by the @xmath0 states that possess one collective bosonic excitation , i.e. @xmath73 the manifold with @xmath74 is formed by all combinations of two of the @xmath0 eigenexcitations @xmath75 with @xmath76 , so that there exist a total number @xmath77 of these doubly excited states . following the same procedure , one could think of constructing states with higher numbers of excitations by applying more creation operators @xmath78 to the vacuum . however , one has to keep in mind that the hamiltonian ( [ eqn : ham_zero ] ) gives a valid description of the system only if @xmath43 , i.e. in the low - occupation sector of the spectrum where only a few bosonic excitations compared to the total number of atoms are present . our aim is to selectively excite the many - particle states described in the previous section starting from the experimentally relevant initial state in which all atoms are prepared in the electronic ground state @xmath79 to this purpose , we use the excitation scheme proposed in refs . @xcite , that consists of two steps :    first , we want to achieve the ground state of the hamiltonian @xmath53 , @xmath80 , which contains zero bosonic excitations . to do so , we express the single - atom bosonic states in terms of the atomic ones , so that the ground state yields @xmath81 . we also introduce an additional internal single atom stable state @xmath82 ( see fig . [ fig : lattice_and_levels]b ) . in practice @xmath82 could correspond to one of the states of the hyperfine ground state manifold which is different from @xmath83 . using this state we can prepare the desired state @xmath84 by means of the two subsequent laser pulses which are depicted in fig . [ fig : lattice_and_levels]b . we choose the first pulse , with rabi frequency @xmath28 , to be resonant with the single - atom transition @xmath85 . after applying a @xmath86-pulse during a time @xmath87 , this laser is turned off and another one with rabi frequency @xmath27 and resonant with @xmath88 is shone on the system during a time @xmath89 . this amounts to the sequence @xmath90\\\nonumber   & \stackrel{\tau_2}{\rightarrow } & \prod_k \frac{1}{\sqrt{2}}\left[\left|g\right>_k-\left|r\right>_k\right]=\prod_k \left|-\right>_k=\left|0\right>\end{aligned}\ ] ] which results in the desired preparation of @xmath84 . here we would like to remark that the strength of the second laser , given by its rabi frequency @xmath27 , has to be large enough to overcome the interaction energy between atoms in the rydberg state . the required strength will depend strongly on the particular geometry of the system . for example , in a square lattice each atom has four nearest neighbors , so that in this case the condition to be accomplished is roughly @xmath91 , while in a triangular lattice it would be @xmath92 as there are six nearest neighbors for each site . once the ground state of the hamiltonian ( [ eqn : ham_zero ] ) is prepared , the aim is to excite the previously described many - particle states ( [ eqn : first_excited ] ) and ( [ eqn : second_excited ] ) that contain one and two bosonic excitations , respectively , thus belonging to higher energy manifolds . the idea proposed in refs . @xcite is to induce transitions between states with @xmath93 making use of the term of the hamiltonian proportional to the detuning , i.e. the first term of ( [ eqn : ham_one ] ) , @xmath94 which creates or annihilates one bosonic excitation at a time . the problem is that , as we discussed before , these transitions are highly suppressed due to the large value of the energy gap - approximately @xmath52 , see fig . [ fig : spectrum ] - between the manifolds . to overcome this difficulty one can introduce a time - dependent detuning of the form @xmath95 which will therefore act as a coherent oscillating source that couples the two many - body states considered when the oscillation frequency @xmath96 matches the corresponding energy gap between them . the intensity of the transition between the two states considered is calculated in each case by working out the square of the transition matrix elements @xmath97 , where @xmath98 and @xmath99 represent the initial and final state , respectively . let us first consider the excitation from the ground state @xmath84 to the first excited manifold of states ( [ eqn : first_excited ] ) that contain a single collective bosonic excitation . the intensity yields in this case @xmath100 where the rotating wave approximation ( rwa ) has been made and @xmath69 is the corresponding energy gap between the state @xmath101 and the ground state . note that the rwa implies that the value of the energy gap - which is roughly proportional to @xmath5 - is much larger than the detuning , i.e. @xmath102 . following the same procedure , one can calculate the intensity of the transition between the one- and two - excitation states , that yield @xmath103 note that this expression already tells us that the intensity of the transition between the two states will be zero unless one of the excitations in the two - boson state is the same as the single one originally present , i.e. unless @xmath104 or @xmath105 . other selection rules that stem from the particular geometry of the system will be analyzed in the following section .      when calculating the intensity of the transitions between the different manifolds following eqs . ( [ eqn : i1 ] ) and ( [ eqn : i2 ] ) , we observe that only for a few states they have a non - zero value ( see figs . [ fig : i1 ] and [ fig : i2 ] ) . the reason for this can be found considering the particular symmetries of the system . let us consider here two examples of geometries : a square and an equilateral triangular lattice .     that contains the symmetries of a square is formed by 8 elements : rotations of the lattice by @xmath106 and @xmath107 ; flip about the @xmath108 and @xmath109 axes ; flip about the diagonals @xmath110 and @xmath111 . * b * : the group @xmath112 contains the symmetries of an equilateral triangle , formed by 6 elements : rotations of @xmath113 and @xmath114 and flip about the axes @xmath115 and @xmath116 that go through the vertices of the triangle perpendicular to the opposite side . ] the symmetry groups corresponding to the square and the equilateral triangle are @xmath117 and @xmath112 , respectively . these groups are formed by the following transformations ( see fig . [ fig : symmetries]a and b ) :    [ [ square ] ] square + + + + + +    * @xmath118 , rotations of @xmath119 , with @xmath120 * @xmath121 , flip about the vertical and horizontal axes * @xmath122 , flip about the two main diagonals    [ [ equilateral - triangle ] ] equilateral triangle + + + + + + + + + + + + + + + + + + + +    * @xmath123 , rotations of @xmath124 , with @xmath125 * @xmath126 , flip about the three axes trough the vertices    one can observe that both the complete hamiltonian @xmath127 ( [ eqn : working_hamiltonian ] ) and @xmath128 ( [ eqn : ham_detuning ] ) conserve all the symmetries corresponding to the system s geometry . as a consequence , if the initial state is an eigenstate with respect to the previous symmetry operators , the time evolution will be restricted to the subspace spanned by the states with the same quantum number with respect to the cited operators . actually , both the experimentally relevant initial state @xmath129 and the ground state @xmath84 where no bosonic excitation is present , are eigenstates of the operators that represent _ all _ the rotations and flip operations that form the corresponding groups with eigenvalue @xmath130 . this is true for both the square and the triangular lattice . the states that have eigenvalue @xmath130 with respect to all the operators representing the group elements belong to the subspace @xmath131 , called _ totally symmetric_. thus , only the states that belong to this totally symmetric subspace can be accessed in the course of the time evolution under hamiltonians @xmath127 and @xmath128 .    . only the states that belong to the totally symmetric subspace @xmath131 have a non - zero probability of being excited . the three states associated to the highest intensities are shown * a * : in a @xmath132 square lattice and * b * : in a equilateral triangular lattice with @xmath133 sites per side . the energies are given in units of @xmath26 . ]    keeping these selection rules in mind , let us now consider the excitation of the states with a single collective bosonic excitation from the ground state @xmath84 and calculate the intensity of these transitions . an example is given in fig . [ fig : i1]a and b for a square lattice of @xmath134 sites and a triangular one with @xmath135 , respectively . here we observe that , as a consequence of the selection rules , the intensity is always zero except for a very restricted set of final states . these states belong to the subspace @xmath131 , i.e. they are symmetric ( eigenvalue @xmath130 ) with respect to all the symmetry operations corresponding to the particular geometry of the lattice . finally , we also observe that the state with highest eigenenergy is the one with the largest intensity . this is observed for all sizes of the system ( determined by @xmath0 ) that have been studied and for both geometries considered . this state reads @xmath136 , with @xmath137 , and is the closest state to a uniformly shared collective excitation , i.e. a spin wave without spatially varying phase factors .    . only the states that belong to the totally symmetric subspace @xmath131 have a non - zero probability of being excited . the density of bosons of the two states associated to the highest intensities are shown * a * : in a @xmath132 square lattice and * b * : in an equilateral triangular lattice with @xmath133 sites per side . the energies are given in units of @xmath26 . ]    let us now discuss briefly the excitation of the states that carry two bosonic excitations . in particular , let us assume that we have already excited the state @xmath138 which corresponds to the largest intensity @xmath139 and calculate the transitions from that state to the doubly - excited manifold . as we already obtained in eq . ( [ eqn : i2 ] ) , the intensity is zero unless one of the two excitations corresponds to the initial one , e.g. in this case @xmath140 . this corresponds to only @xmath0 possible states out of the total @xmath77 available . moreover , the selection rules will reduce the number of accessible states even more due to the fact that only the states that belong to @xmath131 ( totally symmetric states ) can be accessed from @xmath138 , which also belongs to @xmath131 . in particular , a state that carries two excitations only belongs to the symmetry subset @xmath131 when the two excitations that form it belong to the same symmetry subspace . as a consequence , the number of doubly - excited states with non - zero transition intensity is the same as the number of singly excited ones with @xmath141 , which can also be directly inferred from eq . ( [ eqn : i2 ] ) . to understand this better , let us study an example . for a @xmath142 ( @xmath143 ) square lattice there are only three eigenexcitations that belong to @xmath131 , which are @xmath144 and @xmath145 . let us consider that we are interested in knowing which states can be then accessed from the state @xmath146 . in principle any combination of the eigenexcitations @xmath147 with @xmath148 should be accessible . however , we know from the result ( [ eqn : i2 ] ) that only if @xmath149 or @xmath150 the intensity of the transition is non - zero . finally , the selection rules impose that only the combinations of the three eigenexcitations that belong to @xmath131 , i.e. @xmath151 , @xmath152 and @xmath153 can be accessed from @xmath154 with @xmath155 . in fig . [ fig : i2]a and b we show the corresponding intensities for the same two systems considered before , i.e. @xmath132 square lattice and equilateral triangle with @xmath133 sites on each side , respectively . in this case , in the insets of the figures show the bosonic density in each site @xmath156 , i.e. @xmath157 , for the two states corresponding to the highest intensity . once the many - particle atomic states have been excited , we need to map them into a stable configuration . this is necessary due to the limited lifetime of the rydberg states and if one wants to manipulate the corresponding state further . to this end we make use of the stable hyperfine states introduced in the previous section , @xmath82 ( see fig . [ fig : lattice_and_levels]b ) . let us remind that each of the single - atom bosonic states @xmath158 is equivalent to the atomic single - atom one @xmath159 and @xmath160 , respectively ( see section [ sec : hols - prim ] ) . as a consequence , by inverting the sequence of laser pulses given in ( [ eqn : mapping ] ) , we achieve the mapping @xmath161 so that the excitations are stored now in the two hyperfine ground states . expressed in terms of the atomic states , the many - body states studied and given by ( [ eqn : first_excited ] ) and ( [ eqn : second_excited ] ) become @xmath162 and @xmath163 respectively , with @xmath164 . so far , we found the low - lying eigenexcitations in the lattice and described the properties of the many - body states that can be accessed experimentally by controlling the laser parameters . in this section , we will show how to map the collective atomic excitations into quantum states of light @xcite . the mapping from the atomic excitations encoded in the two hyperfine ground states @xmath83 and @xmath82 to non - classical states of light has been thoroughly investigated in @xcite . in the following we will briefly outline the basic idea . we consider the level scheme shown in the lower part of fig . [ fig : lattice_and_levels]b where the two ground states @xmath83 and @xmath82 together with an auxiliary state @xmath165 form a lambda scheme . the storage state @xmath82 is coupled off - resonantly to the state @xmath165 via a classical laser field with rabi frequency and detuning @xmath166 and @xmath167 respectively . eventually , a photon is emitted to the radiation field by means of the transition @xmath168 . note that we do not consider the decay from the state @xmath165 back to @xmath82 , a condition that can be accomplished by an appropriate choice of the atomic levels @xcite . this setup allows us to establish a direct mapping between any given atomic state - encoded in @xmath83 and @xmath82 - and its corresponding photonic one after a time much longer than the lifetime @xmath169 of the intermediate state @xmath165 with @xmath170 being the corresponding decay rate . in particular , one can show that a single atomic excitation is mapped to a single photon state as @xmath171 where @xmath172 creates a photon with momentum @xmath173 and polarization @xmath174 . the details of the coupling are contained in the function @xmath175 , that - in the limit of @xmath176 - is given by @xmath177 let us clarify the meaning of each term appearing in the previous expression . first , @xmath178 is the momentum of the laser and @xmath179 is the energy of the emitted photon , which also appears in the coefficient @xmath180 here , @xmath181 is the quantization volume , @xmath182 the vacuum permitivity , @xmath183 the dipole operator of the @xmath184 transition and @xmath185 the polarization vector ( @xmath186 ) . finally , in eq . ( [ eqn : g ] ) @xmath187 and @xmath188 represent respectively the eigenvectors and eigenvalues of the operator that governs the atomic dynamics in the process of light emission from an ensemble of atoms . this operator accounts for light scattering at multiple atoms as well as interatomic dipole - dipole interactions and depends mainly on the relative orientation of the atomic transition dipole moments and the ratio between the interparticle separation and the wavelength of the laser , @xmath189 @xcite . we will now use this mapping to obtain the single - photon states that are created when using the atomic states ( [ eqn : first_atomic ] ) that contain one collective excitation as a resource . the produced single - photon states are @xmath190 where @xmath191 represents the photon vacuum . in order to characterize the properties of the photonic states , we will study the angular photonic distribution , i.e. @xmath192 where @xmath193 is the number of photons with momentum @xmath173 and polarization @xmath174 . this quantity provides us with the average photon number in the direction of emission parameterized by @xmath194 and @xmath195 . let us now specify the parameters of our setup . we consider the lattice to be lying in the @xmath196-plane . the dipoles of the transitions @xmath184 are all aligned and oriented perpendicular to the lattice , i.e. @xmath197 . this can be ensured by an appropriate choice of the atomic levels . finally , we consider the situation where the incident direction of the laser momentum @xmath178 is also perpendicular to the lattice plane , i.e. @xmath198 . the only degree of freedom left is - apart from the specific geometry of the lattice - the ratio between the interparticle separation given by the lattice spacing @xmath1 and the wavelength of the laser @xmath199 . we will see in the following how the change of this ratio modifies dramatically the properties of the emitted photon .     in the case of a @xmath132 lattice is mapped into a single photon state by means of a laser whose direction is given by the red arrow . * b * : when @xmath200 , the photon is emitted in a superposition of four beams that are perpendicular to the sides of the square lattice . * c * : the angular distribution is strikingly different in the case @xmath201 , where the emission occurs into two pyramidal - shaped beams with square base , as it can be observed in * d*. ]     in the case of a triangular lattice of side @xmath133 is mapped into a single photon state by means of a laser whose direction is given by the red arrow . * b * : when @xmath200 , the photon is emitted in a superposition of three beams that are perpendicular to the sides of the triangular lattice . * c * : the angular distribution is strikingly different in the case @xmath202 , where the emission occurs into two pyramidal - shaped beams with triangular base , as it can be observed in * d*. ]    in figs . [ fig : light_square ] and [ fig : light_triangle ] we show two examples of these single photon sources . in both cases we have chosen to map to a photonic state the atomic many - particle state with a single collective excitation with the highest energy @xmath203 , which also correspond to the largest transition intensities for each geometry ( see figs . [ fig : i1]a and b ) . the corresponding states are represented in figs . [ fig : light_square]a and [ fig : light_triangle]a for a square and a triangular lattice , respectively . the red arrow represents the direction of the incident mapping laser . for each geometry we show the angular distribution of the emitted photon for two different values of the ratio @xmath189 .    in the case of a square lattice and for @xmath200 ( see fig . [ fig : light_square]b ) we see that the photon is emitted in a superposition of four directions , each one perpendicular to the sides of the square lattice . for the same ratio and a triangular lattice ( fig . [ fig : light_triangle]b ) , the photon is emitted in a superposition of three directions , again perpendicular to the sides of the lattice . with this example we see clearly how the geometry of the lattice plays an important role in the properties of the photonic state . these states are particularly interesting since the photon in this case is not only emitted in a highly directional form but also perpendicularly with respect to the incident laser , which would avoid noise problems in the photon detection process .    in order to illustrate how the ratio @xmath189 also affects the properties of the state , we show in fig . [ fig : light_square]c the mapping of the same atomic state in the square lattice but this time with @xmath201 . we observe that in this case the photon is emitted in a superposition of directions that form two cone - shaped beams that still conserve the information of the geometry of the lattice as can be seen in fig . [ fig : light_square]d . the same numerical experiment has been done in the case of the triangular lattice in figs . [ fig : light_triangle]c and d for @xmath202 with a similar result . so far , we have considered a very idealized situation . in particular , we have focused on a regime where the number of excitations is much smaller than the total number of atoms , the laser driving is much stronger than the interaction between the rydberg states and the atoms are distributed in a regular lattice . in this section we will test the validity of the first two approximations made in order to diagonalize the hamiltonian , i.e. the holstein - primakoff approximation and the assumption that the dynamics is constrained to subspaces with equal excitation number . in addition , we will check the robustness of our results with respect to uncertainties in the atomic positions . this will bring us closer to the experimental realistic situation with quantum uncertainty and/or finite temperature .      throughout section [ sec : diagonal ] where we performed the diagonalization of the hamiltonian ( [ eqn : working_hamiltonian ] ) two approximations were made . first , the holstein - primakoff approximation @xmath204 allowed us to neglect all terms of the hamiltonian of higher than quadratic order . second , we considered that in the regime of strong laser coupling , i.e. @xmath205 , the dynamics of the entire system is approximately driven by the hamiltonian ( [ eqn : ham_zero ] ) , which conserves the number of excitations . in order to test the validity of these approximations , we consider now the complete hamiltonian ( [ eqn : working_hamiltonian ] ) and rewrite it as @xmath206 where @xmath46 is the unperturbed energy of the ground state ( [ eqn : ground ] ) and @xmath53 is the unperturbed hamiltonian ( [ eqn : ham_zero ] ) . the perturbation to this hamiltonian can be divided into three terms @xmath207 where @xmath56 and @xmath57 are given in ( [ eqn : ham_one ] ) and ( [ eqn : ham_two ] ) , respectively and the third contribution @xmath208 contains the terms of third and fourth order in @xmath36 that were neglected by the holstein - primakoff approximation @xmath209\right\}\\    & -&\frac{1}{2}\sum_{k\neq m}v_{km}\left(a^\dag_ka^\dag_ka^\dag_ma_k+a^\dag_ka_ka_ka_m\right.\\    & & + \left.a^\dag_ka^\dag_ka_ka_m+a^\dag_ka^\dag_ma_ka_k\right).\end{aligned}\ ] ] let us now calculate the second order energy shifts caused by the three terms of @xmath210 to the unperturbed energies of the ground and first excited states of the hamiltonian @xmath53 . these contributions vanish only in the limits @xmath211 , @xmath212 and @xmath213 . here we consider that @xmath214 , and we calculate the contributions for a finite value of the other two ratios . we focus first on the corrections to the energy of the ground state . @xmath210 only contains terms that couple states that differ by @xmath215 . as a consequence , only the states of the first and second excited manifold ( [ eqn : first_excited ] ) and ( [ eqn : second_excited ] ) contribute to the second order shift to the energy of the ground state , given by @xmath216 note that , since @xmath217 only @xmath56 and @xmath57 contribute to the correction of the energy @xmath218 . this correction yields @xmath219 and depends only on the values of @xmath5 , @xmath0 and the eigenvectors and eigenvalues of the interaction matrix @xmath62 given by @xmath64 and @xmath66 , respectively . the relative error @xmath220 introduced by these shifts for a square lattice of side @xmath221 and several values of the ratio @xmath222 are given in table [ tab : ground ] . one can observe that the dependence with @xmath222 can be very well approximated to be proportional to @xmath223 , as it was already considered in section [ sec : diagonal ] . note that for @xmath224 the shift of the energy of the ground state is below @xmath225 for all sizes of the lattice considered , so that we can consider the unperturbed energy a very good approximation to the exact solution . the dependence with the size of the lattice has a more complicated form , although one can observe clearly that the perturbation becomes stronger for larger sizes @xmath0 . this is expected since the number of states that are coupled to the ground state also increases , and so do the number of terms of the sum ( [ eqn : correction_ground ] ) . .error ( in @xmath226 ) introduced by the second order energy shifts to the ground state energy @xmath46 given in ( [ eqn : ground ] ) for @xmath214 , different values of @xmath222 and sides of the square lattice given by @xmath221 . [ cols="^,^,^,^,^,^,^,^",options="header " , ]      we have considered so far that the atoms populate the harmonic oscillator ground state and that the atoms were infinitely localized , so that the width of the external wave function was negligible compared to the interparticle separation , i.e. @xmath25 . let us now consider a small but finite value of this ratio given by the inevitable uncertainty of the atomic position due to the finite strength of the confinement and finite temperature of the system . we are particularly interested in how the intensities of the transitions from the ground state to the single excitation manifold @xmath228 given by eq . ( [ eqn : i1 ] ) change when considering this uncertainty . to do so , we distribute randomly the position of the atoms around the center of each site of the lattice with a finite standard deviation @xmath229 . as a consequence , the symmetry of the system is lifted , the selection rules discussed in section [ sec : selection ] are not applicable anymore and thus more states of the single excitation manifold become accessible from the ground state . in fig . [ fig : sigma ] we show the resulting absorption profile after averaging over 10000 realizations up to a maximum value of @xmath230 for a square lattice of @xmath231 sites . we observe that as we increase the uncertainty in the position of the atoms the initial sharp lines ( indicated by the red vertical lines ) become broader and blue - shifted . the asymmetry of the shift can be explained as a consequence of the geometry of the lattice . let us first consider the configuration with @xmath232 and that only one atom in the lattice can move from the center of the site . this atom can move in 6 possible directions , i.e. @xmath233 , @xmath234 or @xmath235 . if it moves in the @xmath196-plane , it will end up closer to at least one of its neighboring atoms . since the interaction grows as @xmath236 , the energy in this case will be increased . only the configurations in which the atom moves in the @xmath235 directions will give rise to a lower interaction energy . in total 4 out of the 6 possibilities give rise to a larger energy , which explains the overall shift to the blue in the profile . when the atom is located in one of the boundary sites , the number of configurations that reduce the energy is 3 out of 6 . as the number of sites on the lattice @xmath0 increases , the relative number of atoms on the boundaries decreases and , as a consequence , the shift to the blue is more pronounced . this phenomenon is consistent with the observations obtained from the numerical simulations . ( in units of @xmath1 ) . the red vertical lines show the transition lines for @xmath232 . the horizontal black line indicates the cut shown in * b * : for @xmath237 . * c * : the same cut for @xmath238 . the results are averaged over 10000 realizations . ] in this work , we have studied the excitation properties of a laser - driven gas of highly excited atoms confined to a two - dimensional lattice . in order to solve the system s hamiltonian we have made two approximations . first , we have restricted our study to the states with a small number of excitations compared with the total number of atoms in the system . moreover , we have focused on the regime where the strength of the laser is much larger than the interaction between the atoms in the rydberg state . after these approximations the resulting hamiltonian could be solved by diagonalizing a small matrix which allowed us to obtain the many - body eigenstates of the system and their corresponding energies . moreover , we have described a protocol of how to experimentally access the many - body states in the low - occupation sector of the spectrum by the use of an oscillating detuning . here we have found that only the states with the same symmetry properties as the ground state of the system can actually be excited . we have shown how to map these collective many - body atomic states into non - classical single photon sources and observed that the interplay between the interatomic distance and the wavelength of the emitted photon changes dramatically the angular distribution of the photon state . in particular , we have identified a regime in which the photon is emitted in a superposition of well - separated beams , whose number and relative direction depends strongly on the geometry ( square or triangular ) of the underlying lattice . this feature has potential applications in quantum information processing , especially to form quantum networks . finally , we have considered a small uncertainty in the positions of the atoms on each site of the lattice due to the finite temperature and the quantum nature of the external atomic wave function , which brings us closer to the actual experimental situation . the intensity of the transition from the ground state to the first excited manifold has been calculated for varying values of the uncertainty in the position and we have found a systematic broadening and blue - shift in the profile . the authors acknowledge funding by epsrc . b. olmos also acknowledges funding by fundacin ramn areces .
Abstract:

Answer: the creation of collective many - body quantum states from a two - dimensional lattice gas of atoms is studied . our approach relies on the van - der - waals interaction that is present between alkali metal atoms when laser excited to high - lying rydberg s - states . we focus on a regime in which the laser driving is strong compared to the interaction between rydberg atoms . here energetically low - lying many - particle states can be calculated approximately from a quadratic hamiltonian . the potential usefulness of these states as a resource for the creation of deterministic single - photon sources is illustrated . the properties of these photon states are determined from the interplay between the particular geometry of the lattice and the interatomic spacing .