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: with the increasing demand of the astronomical community for statistical work , all the major observatories equipped their telescopes with spectrographs having multiple object capabilities . we can distinguish between small field instruments , less than 10 arcmin , and large field ones , more than one degree instruments , on which we will focus here . the large field ones use essentially fibre fed spectrographs . this is true for small telescopes ( like the uk schmidt ) as well as for the 4 m class ones . the first generation instruments used focal plates with manually inserted optical fibres in plugboards with drilled holes made at the object coordinates . this system is still foreseen for dedicated instruments , as the sloan dss project(limmongkol _ et al . _ 1992 ) , but for general use facilities , automatic fibre handling machines are preferred . two main devices are in operation , robot and arm positioners . the robot positioner plugs the different fibres in sequence whereas the arm positioner moves all the fibres simultaneously . the robot is able to take care of a large number of fibres , but the configuration time is long and the positioning is done in a blind manner although the robot is equipped with a vision system . this does not guarantee the dropping process . arms are fast , but more than 30 are hardly feasible . moreover , they are expensive , except if built in - house .    when eso decided to replace its optopus ( lund , 1986 ) starplate system with an automatic one , we opted for the arm approach . in that we were following j. hill and his mx ( hill _ et al . _ 1986 ) and t. ingerson and his argus ( ingerson , 1988 ) , our neighbour instrument at ctio . mx uses 32 arms and argus 24 arms . our fibre positioner mefos , the meudon eso fibre optical system ( gurin _ et al . _ 1993 ) , was built by the observatoire de paris , france , under eso contract . mefos is mounted at the prime focus , providing a large field ( 1 degree ) for the multifibre spectroscopy . the 3.6 m eso telescope has a prime focus triplet corrector delivering a field of one degree , the biggest for a 4 m class telescope until the 2df project will be in full operation at the aat . the focal ratio is f/3.14 , well suited for fibre light input , leading to negligible focal ratio degradation ( gurin _ et al . mefos ( gurin _ et al . _ 1993 ) is put on the red triplet corrector with a trait - point - plan interface , allowing tilt correction , with respect to the optical axis . the telescope z movement of the prime focus is used for focusing . it has 30 arms that can point within the 20 cm , one - degree field . one arm is used for guiding , while the other 29 arms are dedicated to the astronomical objects . figure 1 shows the general arms display . they are distributed around the field as `` fishermen - around - the - pond '' . the arms move radially and in rotation : each arm can cover a 15 degree triangle with its summit at the arm rotation axis and its base in the centre of the field . hence , each arm can access objects at the centre of the field , whereas only one specific arm can reach an object at the field periphery . this situation changes gradually from the centre to the edge of the field . each arm has its individual electronic slave board . the whole instrument is controlled by a pc , independent of the telescope control system ( tcs ) . the arm tips  shown on figure 2  are carrying two fibres separated by 1 arcmin . one is used for the object , the other for the sky recording . both are fed to the spectrograph . the object and sky fibre , each 2.5 arcsec in diameter , can be exchanged and the fibre transmission cancelled . +   + coupled firmly to the arm tip is an imaging fibre bundle that covers a 35 x 35 arcsec@xmath0 area of the sky . the image bundles are projected on a single thomson 1024 x 1024 thick peltier cooled ccd , and are connected to the pc driving the arms . the image bundle of the arms are moved to the object coordinates and an image of the selected objects can be obtained . figure 3 shows the display of the galaxy field as seen through the 29 windows of the ccd . +   +    with this procedure , the objects on which the spectral fibre will be positioned can be seen in advance - contrarily to blind positioning . it is to our knowledge the only such instrument . by analysing the position of the object on the image , in combination with the fixed displacement between the image bundle fibre and the spectral fibre within an arm , a precise offset is determined and the arm sent to the optimal object position . this offset takes care of all imprecision due to the telescope , the instrument and inaccurate coordinates . the poor pointing of the telescope , and the fact that the corrector and the instrument are frequently dismounted , would make blind positioning extremely dangerous . the final positioning accuracy , as measured on stellar sources , is 0.2 arcsec rms .    to avoid any motion between the telescope and the instrument , the guiding is done in the arm environment . a special arm is dedicated for this task . this arm holds a large image conducting optical fibre with a field of view of 1.8 x 2.5 arcmin . to enlarge the field within which a guide star can be selected , this window can be positioned in 3 different ( adjacent ) positions . the guiding camera is sensitive to v@xmath1 . because of this faint magnitude limit , no guide star has to be defined prior to the observations . the guiding fibre can then remain in the outer position at the edge of the field , leaving more space to the positioning of the arms .    in the present stage , the spectral fibres , @xmath2 in diameter and 21 m long , are going from the prime focus down to the cassegrain , where the b&c eso spectrograph is located . this spectrograph is fitted with a f/3 collimator to match the fibre output beam aperture . it has a set of reflection gratings with dispersion range from 2.3 nm / mm to 23.1 nm / mm , and a tek 512 x 512 thin ccd . the first scientific observation was performed in february 1991 with a mefos prototype consisting of nine arms only ( bellenger _ et al . _ 1991 ) . during this period a galaxy of magnitude b=@xmath3 was observed with a resolution of r=370 . the 1@xmath4 exposure resulted in spectrum with a s / n=50 in the red part of the continuum ( 550 nm ) . the first technical run in full configuration at the 3.6 m telescope took place in october 1992 . as only 2 1/2 of the 5 allocated nights were clear , a second run was necessary for the testing and the commissioning . this was carried out in february 1993 during 5 clear nights . some problems with the centering of the image fibres on the object coordinates were noted . it appeared that the optical aberrations correction curve was not accurate enough . we had adopted this curve as it was derived from photography at the prime focus . the mounting of this instrument entails different conditions . by moving an image fibre across the mefos field in real configuration , the correction curve could be adjusted . the objects are now falling within a few arcseconds of the centre of the image fibre helping for their identification . in the 7 nights following this run , mefos was used on two scientific programmes . the instrument worked well and little time was lost due to technical reasons . cuby & mignoli ( 1994 ) studied sky subtraction procedures using mefos with and without beam switching . they were able to measure redshifts for galaxies with magnitudes as faint as b=22@xmath5 from a 90 min exposure ( 2 x 45 min ) . a b=@xmath6 galaxy was observed as well , but the resulting signal to noise ratio of s / n = 7 was too low for a reliable redshift measurement .    since then , mefos was used in 84 nights , without a single failure . since april 1994 , mefos has been used without any further support from the building team . in order to compensate the building team of the daec for their efforts , 32 nights of guaranteed time were allocated to scientists of this laboratory . from the submitted scientific proposals , four programmes were selected . each was granted 8 observing nights with mefos . based on the experience gained in one of these programmes , we describe in the following the preparations before an observing run as well as the typical procedures while observing with mefos at the 3.6 m telescope . in the next section we will then discuss a number of specific problems that we encountered while reducing the data , and , at the end , a number of examples and some statistics will allow some insight into results that can be expected from multifibre spectroscopy with mefos . the described procedures as well as the obtained data are , of course , strongly related to the selected project . it is therefore important to know some details about the project when discussing the performance of mefos . the scientific project `` search in the galactic plane toward the great attractor '' aims at mapping the galaxy distribution in velocity space in the zone of avoidance ( zoa ) in the general direction of the cosmic microwave background dipole ( cmb ) and the great attractor . galaxies hidden behind the zoa are important in derivations of the peculiar motion of the local group relative to the cmb , and in explaining the observed streaming ( bulk flow ) motions , as e.g. in the great attractor region . preliminary scientific results have already been presented ( cayatte _ et al . _ 1994 , kraan - korteweg _ et al . _ 1994 , 1996a , b ) . this programme is a continuation of a long - term programme which made use of optopus , the previous multi - fibre spectrograph of eso . a large part of the extragalactic sky is obscured by the milky way . to bridge this gap , we have embarked on a search for galaxies behind the milky way to the low diameter limit of d=@xmath7 and faint magnitude limit using existing sky surveys ( kraan - korteweg & woudt 1994 ) . to date over 10000 previously unknown galaxies were uncovered in the southern milky way ( @xmath8 , with a strong peak @xmath9 ( cf . fig . 1 in kraan - korteweg & woudt , 1994 ) . this range is well suited to the dynamic range of mefos . galaxies in this magnitude interval have an average number density of the order of 100 - 200 galaxies@xmath10 , suggesting a more efficient use of optopus with 50 fibres in a half - degree field . however , the galaxies of this project lie behind the milky way . they are heavily obscured : around one magnitude at the border @xmath11 of our survey , increasing to 5 - 6 magnitudes towards lower latitudes ( close to the galactic dust equator the galaxies are fully obscured ) . the average number density for @xmath12 magnitude galaxies is therefore considerably lower . in fact , the multifibre capacity of mefos with 29 objects per one - degree field in combination with its capability of obtaining good s / n spectra for these faint obscured galaxies of our survey , make it the ideal instrument to trace large - scale structures close to the plane of the milky way . the galaxies of the zone of avoidance survey not only are very obscured and often of low surface brightness , but since they are located in the milky way , the fields are crowded with foreground stars , many of which are often superimposed on the faint galaxy images . here , the imaging facility of mefos for object identification and optimal centering of the spectral fibre is again a most advantegous property ( as for spectroscopy of lsb objects and/or in crowded fields in general ) .        to prepare an observation , it is necessary to have a list of objects with coordinates precessed to the date of observation and known with a precision of the order of one arcsec . arm - object assignation software is provided which can be obtained by ftp from eso . this work is optimally done at the home institute before the observing run and does not cost telescope time . the software is also implemented at la silla , allowing new or optimization of field assignations during the observing run . the arm assignation software displays the objects of an earlier prepared list on the screen . a circle of one degree ( the size of the mefos field ) can be moved around the screen until optimal centering of the field is attained . the computer then suggests a possible configuration for an arm - object assignation . the software does take the constraints of size and moving capacity of the arms into account , to avoid collisions . the smallest allowed distance between objects depends on the relative arm positions . it lies between 2 and 5.8 arcmin . if the proposed assignation is not optimal from a scientific point of view , it can be modified interactively ( within the above named limits ) . figure 4 shows an example of an assigned field . the arms are displayed with their true extent . +   +      mefos is installed with the corrector at the prime focus . the north - south direction has to be known with a high precision . a set of star fields , equally spaced in right ascension , are used during the set - up for calibration . the pointing of the telescope is essential when observing with mefos . during the set - up , the pointing will be checked in detail . this will lead to a pointing model which is automatically applied during the whole observing run . it is therefore not necessary to recheck the pointing . any change in the pointing should be avoided , as this will destroy the pointing model . the pointing only has to be rechecked if the system had , for some reason , to be rebooted . the instrument is focussed by pointing the guiding arm to a star . this fibre has a larger field of view compared to the object image fibres . the telescope is pointed towards the selected field . as soon as the the telescope is locked on the centre of the mefos field , the guiding has to be switched on . there are 3 adjacent guide fields of 1.8x2.3 arcmin each , in which a guide star can be found . the magnitude of the star can be as faint as 19th magnitude . in the exceptional case that no guide star can be found in any of the 3 possible guide field positions , it is better to switch to another target field in order not to lose valuable observing time . an arm - object reassignment to enhance the likelyhood of finding a star in the guiding window can be done at a later point of time . the arms are then sent to their previously assigned object positions , with the image fibre bundles centered on the objects . during the arm positioning , the pc controlling the mefos operations , continously displays the movements of the arms . it typically takes about 4 minutes for the arms to reach their final positions . as soon as the arms are in position , a ccd image is taken ( see fig . the images are crucial , because they allow object identification and accurate centering of the spectral fibre . these images are raw and can not be used for other scientific purposes like photometry . the readout time for the ccd image is 17 sec . the noise is quickly sky dominated . an exposure time for a v=17@xmath132 star is of the order of 20 sec with 1.7 arcsec seeing . for a b=20@xmath5 galaxy , depending on the shape , about 3 min are needed . we generally exposed as long as 5 or 7 minutes because the objects of our programme are obscured galaxies , often of low surface brightness in fields crowded by foreground stars . the identification of our objects in the 35``x 35 '' field was therefore often difficult . for brighter objects in less crowded fields the advocated 2 to 3 minutes is adequate . after the images are taken and displayed on the screen , the observer has to identify the object on the respective images for final optimal positioning . this can be done automatically or manually . an internal routine searches for the brightest objects within the images and puts numbered boxes around them . the correct object can then be identified by entering the number of the box into the computer . the optimal position can also be indicated manually by clicking with the cursor on the required position . the latter is a more efficient procedure for centering of very low surface brightness objects . depending on the number of objects in the field of view and on the surface brightness of the objects , this step takes between 5 and 10 minutes . if the objects are very faint or multiple , it is advisable to make polaroid copies and mark the position where the spectral fibre will be put . it takes about one minute to send the spectroscopic fibres to the recentered objects .      after the identification step , the observing procedure depends on the scientific programme . a flat field with an internal white lamp allows to identify and to check the position of the fibres on the ccd . the exposure time is about 10 sec . the he and ne lamps are used for the wavelength calibration ( note that the neon lamp has to be requested before the setting ) . it is recommended to have a filter ( bg 28 ) put in front of the neon to reduce the intensity of the red neon lines and avoid saturation . the two lamps can be exposed simultaneously . based on our experience , an exposure time of about 1 min is reasonable . then the scientific exposures can begin . depending on the programme , magnitude and surface brightness of the objects , a sequence of 30 min exposures follows . with more than two consecutive exposures of one field , another calibration was taken to test the stability of the wavelength calibration . this whole procedure has to be repeated for every field . the overhead between one scientific exposure of a field to the next field typically takes 20 to 30 min , calibration included . with two observers at the telescope some of these operations can be done simultaneously , saving therewith on real observing time . for an exposure time of 2x30 min , 5 to 6 fields per night can be observed depending on the length of the nights . we observed 9 fields in 2 1/2 nights in february 1993 , 19 fields in 4 nights in february 1994 and 19 fields in 4 nights in may 1995 . we used the ccd tek # 32 in our observing runs with the grating # 15 . this resulted in a dispersion of 170 / mm and a resolution of about 11 . the wavelength range was chosen to be 3850- 6150 . the data reduction has been carried out on sun stations with the iraf package . the reduction process consists of subtraction of a mean bias , extraction of the spectra , wavelength calibration , fibre transmission correction , sky subtraction , removal of cosmic - ray events and redshift determination . for the last step we used the rvsao package obtained from the smithsonian astrophysical observatory developed by doug mink . with this package , radial velocities can be obtained with the cross - correlation method proposed by tonry and davis ( 1979 ) , and redshifts can also be determined from emission lines if present . in the course of our observing runs , we regularly observed a number of standard velocity stars ( 1 - 2 per night ) which were subsequently used in the cross - correlation of each object spectrum . we will not give a detailed description of the various steps of the the reduction process  this will be given in the forthcoming paper on the scientific programme  but rather highlight some specific problems . efficient extraction of spectra and accurate wavelength calibration are critical in measuring accurate redshifts . this does require deep understanding of the performance of the system such as the illumination of calibration and flat - field sources . in the following we point out some problems that we encountered while reducing the multifibre data .      for each fibre the position of the spectrum on the ccd has to be defined . this is done by measuring the centroid of the profile perpendicular to the dispersion . the curvature in the dispersion direction is corrected by fitting the centroid position for each object as a function of the coordinate along the wavelength axis . the positions of the spectra are determined from a flat field exposure in the same telescope configuration as the objects of the field exposure . no `` cross - talk '' between two adjacent spectra has been evidenced , so an extraction window as large as the spectra separation can be chosen . in our case we use an extraction window of 5 pixels for each aperture ( adding faint signal in the profile wings reduces the s / n ratio for the object spectrum ) . spectra for each aperture in each exposure were then extracted with `` variance '' weighting . the dispersion was determined for each aperture ( e.g. each spectrum ) from the corresponding aperture on the he - ne lamp calibration exposures . for the dispersion solution we used 13 emission lines from the he - ne lamp . we fitted a cubic spline of order 3 to the non - linear component . the residuals were not larger than 0.3  . by measuring the scatter in the skyline position , we always obtained an uncertainty in the wavelength calibration below 0.5  ( e.g. 30@xmath14 ) . the total flexure between the ccd mounted on the b&c spectrograph ( mounted at the cassegrain focus ) and the fiber focus has been checked in both directions . to look for positioning repeatability of the spectra , we have compared two flat field exposures of different telescope pointing . after both the telescope and all arms were moved to other locations , differences between the positions of the spectra for two flat - field exposures were found , though never more than 0.25 pixels ( e.g. 7 ) . flexure with telescope tracking was checked by comparing two flat field exposures taken with the same telescope pointing . the maximum noted flexure in an hour never exceeded 0.25 pixels , i.e. 7 .    in the dispersion direction we compared two lamp calibration exposures taken before and after successive field exposures . the 1993 mefos run revealed a shift of about 1  per hour depending on telescope position . this corresponds to a quarter of a pixel , i.e. is of similar order as in the spatial direction . this effect was seen in almost all the fields . it could be due to a continuous drift with time . the february 1994 mefos run did not exhibit any such effect . while reducing the fields of the february 1993 and 1994 runs , we noticed inconsistencies between the positions of the spectra obtained from the flat field as compared to the objects : the difference in position between the flat and the object is roughly constant up to the middle of the ccd frame and then decreases as a function of pixel position in the spatial direction up to -0.5 pixels ( e.g. 13 ) . this trend is displayed in figure 5 . although the mean curves of the different fields are slightly displaced from each other , they all clearly exhibit the same trend . moreover the mean shift changes as a function of wavelength . it is stronger in the blue part of the image . the problem therefore does not seem related to the mechanical stability and reliability , but rather to the difference in the light path of the lamp used for the flat field exposure and the object exposure . cuby and mignoli ( 1994 ) noticed displacement problems as well . however , they find a difference between flat - field and sky spectra positions which varies linearly with the centroid of the spectra .    because of these difficulties , the mefos - user is adviced to check for such systematic effects and not to blindly use the flat fields for defining the spectra positions . to circumvent these problems , we extracted the spectra positions directly from the object spectra . the two exposures per field were summed up if no shift between the positions on the ccd frames could be found . this addition increases the signal to noise ratio in the less sensitive , blue part of the spectra . the errors in the position of the spectra in the faint blue part are of the order of 0.3 pixels ( as compared to the fwhm of the spectra profile of about 3 pixels ) . a large extracting window , on the same order of the spectra separation , will decrease the loss of signal due to the imprecise determination of centroid of the profile . + for certain fibres the spectra were not properly wavelength calibrated . for these spectra the [ oi ] skyline wavelength position was not close to the actual value of 5577.4  . so we also checked the alignment of the fibres along the slit by comparing he - ne lines positions for different apertures ; the result is shown in the bottom panel of fig . 6 : the mean spectral - shift of all the he - ne lines in the dispersion direction has been plotted versus the aperture position coordinate on the spatial axis . the points seem to be located around a parabola which represents the image of the slit through the spectrograph . in the top panel of fig . 6 , we have plotted the differences between the measured position of the [ oi ] skyline and 5577.4  versus the aperture position coordinate on the spatial axis . the largest displacement of the skyline position occurs for the fibres which are the most misaligned on the slit . as suggested in the previous paragraph , different optical paths for sky and lamp exposure are suspected . here we observe this effect in the dispersion direction . we can not wavelength - calibrate the spectra correctly for objects observed with fibres that are located far from the slit center in the spectral direction . +    the optical paths difference could be explained if misaligned fibres have undergone some bending effect or that the plane of the slit ( in front of the fibers ) is not perfectly perpendicular . this would change the position of the image if the ccd is illuminated by the sky or by a calibration lamp .    to correct the bad wavelength calibration we estimated a mean shift by averaging the [ oi ] skyline offsets . these were measured from each exposure in a field . we then applied an inverse shift of the corresponding amount to each spectrum taken in the `` misaligned '' fibres . the largest discrepancies between measured and actual value was consistently seen for the same fibers in all the 18 observed fields of the 1994 run . the r.m.s error of the mean skyline offset was never larger than 0.17  . generally , it was of the order of 0.12  , corresponding for a wavelength of 5577.4   to an r.m.s . error of 6@xmath14 in velocity . we estimated the relative transmission of each fibre from the [ oi ] 5577.4  skyline flux to correct the fiber efficiencies . after this correction , we averaged all the sky spectra in a given field . this minimizes the noise of the sky spectrum to be subtracted from the object spectra . with mefos , each arm carries its own sky fibre , allowing therewith other sky subtraction methods . detailed investigations on optimal sky subtraction routines have been performed by cuby and mignoli ( 1994 ) . they find averaging over all sky spectra , after fibre efficiency correction as determined from the skyline flux the best solution . this does of course depend on the programme . for spectrophotometry of galactic or extragalactic objects , for instance , precise flux calibration will be required . the fibre relative transmission coefficient deduced from the [ oi ] 5577.4   skyline flux takes into account the sky variations over the field together with the effective fibre efficiency determination , without the possibility of disentangling the two effects . generally , sky variation is rather weak below 6000 even on a one degree field ( see results of cuby and mignoli , 1994 ) . the transmission variations observed for some fibres during the night as well as from night to night , could therefore be largely due to a change in fibre efficiency , which can , for instance occur when light injection conditions vary . nevertheless , the transmission coefficient values obtained from two exposures on the same field show very little scatter . from the statistic on the fibre transmission difference between the two exposures we derived an r.m.s . error on the estimation of this coefficient below 2% . in 1993 and 1994 , 506 objects have been observed on 27 fields . this translates to an average of 19 objects per mefos field . this number is considerably lower than the maximum of 29 objects within a one degree circle . on the one hand , this reflects the choice of a field , on the other hand the constraints of the instrument itself . besides scientific priorities , constraints are given by the arms ; i.e. the width of the arm and the minimum allowed distance . the average number of objects that we could observe in crowded fields was of the order of 22 to 24 galaxies . numbers of 27 to 29 were indeed achieved if the galaxy density reached a hundred or more , as in fields of the cluster abell 3627 ( kraan - korteweg _ et al . _ , 1996a ) .      of 473 spectra in 24 reduced fields , redshifts have been obtained for 93% of the objects , among which 84% are of galaxies and 9% are of stars . + the perturbations of the galaxy spectrum by foreground stars ( of the 47 spectra contaminated by a superimposed star , only 6 redshifts of galaxies could be obtained ) is , of course , a problem inherent to our programme and should not affect observing programmes away from the galactic plane . among the redshifts obtained for galaxies , 77% are reliable ( i.e. with a s / n ratio larger than 20 at 5000  ) with errors in velocities which in the mean are about @xmath15 . for 22% of the objects , the redshift could be attained from emission lines . the distribution of the good s / n redshifts are displayed in fig . 7 for the velocity range of @xmath16 .    the distribution in fig . 7 depends strongly on the extragalactic large - scale structures that are being intercepted in the milky way . for instance , the strong peak at about @xmath17 mainly stems from 4 mefos fields centered on the cluster a3627 at @xmath18 . the velocities from these 4 fields are marked by the dashed histogram . in fact , the velocities from these observations , in combination with the result from our complementary observations at the saao and at parkes , have proven this cluster to be very massive and therefore to be at the bottom of the potential well of the great attractor ( kraan - korteweg _ et al . _ 1996a ) . +    the mean velocity of the observed sample is @xmath19 . the peaks at 12000 , 18000 and 23000@xmath14 are part of filaments that diagonally cross the zone of avoidance . they can be traced to higher latitudes on both sides of the milky way and are indicative of extragalactic large - scale structures of considerable size ( a few hundred mpc ) . for further details see kraan - korteweg _ et al . _ ( 1994 & 1996b ) . +    here , it suffices to maintain that mefos produces good quality spectra for galaxies as distant as @xmath20 . the above results do depend on the properties of the observed galaxies like the magnitude and the mean surface brightness , although the best parameter would be the central surface brightness corresponding to the flux collected at the position of the fiber . one can imagine a very low mean surface brightness galaxy with a relatively bright nucleus , giving a good s / n spectrum .    in fig . 8 the distribution of the apparent magnitudes b@xmath21 of 361 galaxies with redshifts is displayed . the magnitudes are estimates from the iiiaj film copies of the src . although eye estimates , it has been shown that these magnitudes are accurate to @xmath22 with no deviations from linearity up to the faintest galaxies . kraan - korteweg and woudt , 1994 , kraan - korteweg 1996 ) . the mean magnitude of the observed galaxy sample is @xmath23 . if this distribution is compared to the one of the galaxies for which no redshift could be extracted ( dashed histogram in fig . 8) , no clear bias is visible toward the faintest objects . therefore the derivation of a redshift in this magnitude range does not critically depend on magnitude . the distribution shows a strong peak between @xmath24 and @xmath25 which does reflect the magnitude distribution of the galactic plane survey ( cf . fig . 1 in kraan - korteweg & woudt , 1994 ) . + a slightly stronger dependence has been noted with mean surface brightness @xmath26 . the mean surface brightness ranges between @xmath27 with a narrow gaussian distribution around @xmath28 ( @xmath29 ) , whereas the galaxies for which no redshift could be determined have a mean of @xmath30 ( @xmath31 ) . but a redshift could be derived for the faintest few galaxies as well as for the one with the lowest surface brightness . a weak correlation between velocity and surface brightness was noted : while the galaxies with @xmath32 show a stable distribution about the mean , the galaxies with @xmath33 reveal a decrease in mean surface brightness with increasing redshift . + in figure 9 , a few characteristic examples of galaxies from our survey and their resulting spectrum from a 2x30 min exposure with mefos are shown . they illustrate the quality of the obtained spectra as a function of the properties of the observed galaxies and visualize the above discussed tendencies . three of the four objects presented in fig . 9 have a @xmath34 magnitude nearly equal to the mean one of the whole sample , but the mean surface brightness spans a large range . the object shown in the lower - left panel is one of the faintest galaxy of our sample , a 20@xmath35-magnitude spiral galaxy with a very low mean surface brightness ( @xmath36 ) . the ccd images of the galaxies clearly show that the central surface brightness is the most important factor in achieving a high s / n ratio . the velocities of the two objects in the upper panels are of the order of the mean velocity obtained for the sample , but the galaxies in the lower panels are more distant . +    > from the overall analysis of the observed sample , it can be maintained that good signal to noise spectra can be obtained with mefos from 2x30 min exposures for galaxies up to a magnitude of @xmath37 , the approximate magnitude limit of our survey . for these obscured galaxies behind the milky way , we have measured redshifts of galaxies up to @xmath38 . reliable redshifts can be expected for a percentage of over @xmath39 . this efficiency does depend on surface brightness and , especially in the case of this programme , is limited by interfering foreground stars , an effect which would be much stronger were it not for the imaging capability of mefos . in this paper we present the features of mefos , a multi - object fibre positioner coupled to the b&c spectrograph and implemented at the prime focus of the eso 3.6 m telescope at la silla . this instrument is open to the eso community and was used on the sky for 84 nights , until now . mefos is an arm type instrument which is able to position 29 fibres on designed targets and has the same number of sky fibres . very precise fibres positioning with this instrument is only possible on this generation of telescopes with a vision system that is associated to the arms . mefos aim is the observation of moderately crowded fields that need a one degree field for efficient observations . several scientific programmes have been conducted with mefos . the results obtained for one of them , aiming at measuring redshifts of galaxies found in the galactic plane , is given as an example of what can be observed with mefos , and what were the encountered problems with the data reduction . bellenger r. , dreux m. , felenbok p. _ et al . _ : 1991 , _ the messenger _ , * 65 * , 54 + cayatte , v. , balkowski , c. , kraan - korteweg , r.c . : 1994 , in _ unveiling large - scale structures behind the milky way , 4@xmath35 daec meeting _ , eds . c. balkowski & r.c . kraan - korteweg , a.s.p . , * 67 * , 155 + cuby j. g. , mignoli m. : 1994 , _ proc . spie _ , * 2198 * , 98 + gurin j. , bellenger r. , dreux m. _ et al . _ 1993 , in _ fiber optics in astronomy ii _ , ed . p.m. gray , 145 + gurin , j. , felenbok , p. : 1988 , in _ fiber optics in astronomy_,ed . barden , a.s.p . , * 3 * , 52 + hill j. m. : 1986 , _ proc . spie _ , * 627 * , 303 + ingerson t. e. : 1988 , in _ instrumentation for ground based optical astronomy _ l.b robinson ( new york : springer ) , 222 + kraan - korteweg , r.c . : 1996 , in prep . + kraan - korteweg , r.c . , cayatte , v. , balkowski , c. _ et al . _ : 1994 , in _ unveiling large - scale structures behind the milky way , 4@xmath35 daec meeting _ , eds . c. balkowski & r.c . kraan - korteweg , a.s.p . , * 67 * , 99 + kraan - korteweg , r.c . , fairall , a.p . , balkowski , c. : 1995 , _ a&a _ , * 297 * , 617 + kraan - korteweg , r.c . , woudt , p.a . : 1994 , in _ unveiling large - scale structures behind the milky way , 4@xmath35 daec meeting _ , eds . c. balkowski & r.c . kraan - korteweg , a.s.p . , * 67 * , 89 + kraan - korteweg , r.c . , woudt , p.a . , cayatte , v. _ et al . _ : 1996a , _ nature _ , vol . * 379 * , 519 + kraan - korteweg , r.c . , woudt , p.a . , fairall , a.p . _ et al . _ : 1996b , in _ xvth moriond astrophysics meeting on `` clustering in the universe '' _ , eds . s. maurogordato _ _ , 71 + limmongkol , s. , owen , w.a . , siegmund , w.a . , hull , c.l : 1993 , in _ fiber optics in astronomy ii _ , ed . p.m. gray , 127 + lund g. : 1986 , _ eso operating manual _ no . * 6 * +
Abstract:

Answer: we are describing a new multi - fibre positioner , mefos , that was in general use at the la silla observatory , and implemented at the prime focus of the eso 3.6 m telescope . it is an arm positioner using 29 arms in a one degree field . each arm is equipped with an individual viewing system for accurate setting and carries two spectroscopic fibres , one for the astronomical object and the other one for the sky recording needed for sky subtraction . the spectral fibres intercept 2.5 arcsec on the sky and run from the prime focus to the cassegrain , where the b&c spectrograph is located . after describing the observational procedure , we present the first scientific results .


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: dimerized and frustrated spin chains have attracted considerable interest in recent years . this is due to the recent emergence of a variety of experimental quasi - one - dimensional systems containing localized electrons that can be described by spin chains . among these systems there is a sizable number which are gapful due to dimerization in their low temperature phase , such as cugeo@xmath4 @xcite , @xmath5-nav@xmath6o@xmath7 @xcite , ( vo)@xmath6p@xmath6o@xmath8 @xcite , cu(no@xmath4)@xmath9h@xmath6o @xcite , cuwo@xmath10 @xcite , and cu@xmath6(c@xmath6h@xmath11n@xmath6)@xmath6cl@xmath10 @xcite . the first two dimerize as a consequence of the coupling to the lattice degrees of freedom , whereas the other compounds are intrinsically dimerized , i.e. the dimerization does not depend on temperature . the effect of doping on these substances represents an interesting issue , in particular since it pertains also to defects such as missing spins , broken chains and chains of finite length .    as a general rule , defects in low - dimensional antiferromagnetic spin systems which can be described by a rvb - type ground state @xcite involve @xmath0 states in their vicinity @xcite . the appearance of a certain impurity spin at the edges of spin chains can also be discussed in the framework of a nonlinear @xmath12 model for general spin @xmath13 @xcite . we restrict ourselves to @xmath0 , weakly dimerized spin chains . the basic idea there is that without the defect each spin has a partner with which a singlet is formed . if the defect removes one of the partners , the other becomes a free @xmath0 spin which we will henceforth call spinon to distinguish it from the other singlet forming spins . we understand that this spinon comprises also some dressing of the bare @xmath0 spin , i.e. it is not localized at just one site . since the singlet pairing is not static in the rvb - picture the spinon is able to move some distance away from the defect . in a critical , gapless system , it will be delocalized . this delocalization , however , disappears as soon as the couplings are modulated . it is the purpose of this work to elucidate in which way an explicit dimerization acts as a confining potential for the spinon motion leading to its localization . the generic hamiltonian reads @xmath14 where @xmath15 parametrizes the dimerization and @xmath16 the relative frustration by next - nearest neighbor coupling . previous works , e.g. @xcite , viewed dimerization already as confining potential . we understand that there is no confinement without dimerization ( @xmath17 ) . here we will develop a quantum mechanical picture that treats the cases of small and large frustration on equal footing . although this picture is not exact in all details , nevertheless it is able to reproduce the main features of the problem on a semi - quantitative level . the physical quantities considered here are binding energies and local magnetizations . in addition , it will be shown that the confinement depends on the degree of frustration and it will be sublinear in the region of low frustration . for the sake of concreteness , the case of a chain end is illustrated in fig.[fig - idea ] . ( 8.2,5.2 ) ( 0,0.5 )    dimerization localizes spin singlets mostly at the strong bonds . if in a dimerized chain one spin is missing its singlet partner is freed . if the missing spin had a weak bond to the left ( assuming without loss of generality horizontal chains ) and a strong bond to the right , the free spin situation corresponds to the one in fig.[fig - idea](a ) . the configuration in fig.[fig - idea](b ) is then the reflected configuration found on the left of the missing spin . if the missing spin had a strong bond to the left and a weak bond to the right , the free spin situation is found to the left of the missing spin . the configuration in fig.[fig - idea](a ) is then a reflected image whereas the configuration in fig.[fig - idea](b ) is found to the right of the missing spin . fig.[fig - idea](c ) illustrates the situation after two hops of the free spin away from its origin in ( a ) . the crucial point is that the singlets to the left of the spin in fig.[fig - idea](c ) are no longer at the strong bonds . this implies an energy loss which can be viewed as an attractive potential which ties the free spin to its origin . besides chain ends also solitonic modulations , such as kink defects , will be considered in this paper . they can also be viewed as spinon traps with confining potentials . the article is organized in the following way . in the next section , the short - range rvb spinon states that determine the main part of the correlations will be introduced . subsequently , the forms of the kinetic and the potential energies which govern the spinon dynamics will be discussed , i.e. the schrdinger equation of the problem will be setup . in the fourth and the fifth sections the quantitative results will be discussed and compared against computational calculations for chain ends and for kink - modulated chains , respectively . a summary will conclude the article . let us denote by @xmath18 an up - spinon at site @xmath19 ; the other spins are all paired to nearest - neigbor singlets . in this convention , the state in fig.[fig - idea](a ) is denoted by @xmath20 and the state in fig.[fig - idea](c ) is denoted by @xmath21 . these states are not orthogonal but their overlap @xcite is given by @xmath22 this overlap arises from a nel type sequence of up - spins and down - spins between the two spinons . let us consider a state @xmath23 for which we aim to obtain a continuum description . this means that we assume that for the low energy behavior it is sufficient to treat @xmath24 as a slowly varying function of @xmath25 . note that we introduced the factor @xmath26 to focus on the energetically low - lying states . this can be seen most easily for @xmath27 @xcite . in order to define a normalized wave function @xmath28 let us calculate the norm of @xmath29 @xmath30 the approximate step is valid if @xmath24 varies in fact slowly with @xmath25 . eq.([norm ] ) tells us that passing to a normalized continuous wave function @xmath28 means @xmath31 and @xmath32 . the factor 2 introduced here enables the replacement of @xmath33 by the actual site numbers at the end of the calculation . now we calculate the local magnetizations @xmath34 . @xmath35 a contribution is found only if @xmath36 is situated _ between _ the two spinons which equals one half of the overlap ( [ overlap ] ) . for @xmath37 odd one finds @xmath38 where @xmath39 and @xmath40 stand for half the distance to site @xmath37 on the left and on the right . we do not take the chain end into account , i.e. we ignore the fact that summation for small @xmath37 is truncated on the left . for even sites @xmath37 one obtains @xmath41 with the above calculation the local magnetizations have been linked to the probability of finding the spinon at a given site @xmath37 . the sum @xmath42 equals to one half as it has to be for a global @xmath0 state . note that to obtain this result one has to sum even and odd sites separately , and that @xmath37 changes , thus , by two from site to site . summing the moduli @xmath43 one gets immediately @xmath44 . this means that the antiferromagnetic correlations induce total local moments that correspond to the moments of @xmath45 independent spins @xmath0 . this enhancement factor of @xmath46 illustrates why even a low concentration of dopants may induce considerable antiferromagnetism in a compound @xcite . moreover , this enhancement is not due to criticality but it is already built - in in the short - range rvb state . the possible coexistence of dimerization and alternating local magnetization was nicely demonstrated by fukuyama _ by analysis of the corresponding phase hamiltonian @xcite . impurity - induced antiferromagnetism is already extensively investigated for spin ladders @xcite .    in the next paragraphs we will refine the calculation of the norm and the local magnetization by considering indeed a discrete wave function . we would like to remove the non - orthogonality of our basis states in order to be on safer ground for the subsequent reasoning . moreover , certain boundary effects can be captured by this procedure .    by gram - schmidt orthogonalization we find the orthonormal basis    @xmath47 we have    [ lodd ] @xmath48    and the analogous formulae hold for @xmath49 . all other expectation values are zero . from the relations ( [ lodd ] ) one finds for a normalized state @xmath50    [ resodd2 ] @xmath51    repeating the same for even @xmath52 we arrive at    [ leven ] @xmath53    from eqs . ( [ leven ] ) we obtain @xmath54    note that for @xmath55 and slowly varying @xmath56 one re - obtains , of course , the results eqs.([resodd1 ] ) and ( [ reseven1 ] ) . this concludes the calculation of the local magnetizations . in this section we address the important issue of how the simple spinons introduced in the previous paragraph move and how they are attracted by the chain end . we will distinguish two regimes of frustration ( i ) @xmath57 @xcite and ( ii ) @xmath58 .    in regime ( i ) we know that the spinons have a finite mass and that they have a quadratic minimum in their dispersion @xcite . moreover , the strongly frustrated chains display spontaneous symmetry breaking of the translational symmetry , i.e. spontaneous dimerization occurs . the actual value of the gap ( and the dimerization ) , however , is very small up to @xmath59 @xcite ( see also fig.[fig - gap ] ) .    in regime ( ii ) the spinon dispersion is linear in the wave vector @xmath60 @xcite . hence we assume @xmath61 where @xmath62 is the spin wave velocity which is roughly given by @xmath63 according to ref.@xcite . the multiplication by @xmath64 in reciprocal space can be visualized as the consecutive multiplications @xmath65 . to multiply by @xmath60 corresponds to @xmath66 in real space , while a multiplication by @xmath67 corresponds to the convolution with the principal value of @xmath68 . then , @xmath64 in reciprocal space corresponds to the convolution    @xmath69    [ kinetic ]    the symbol @xmath70 stands for the principal value which has to be taken . if one has to discretize the singular operator @xmath71 , the most reliable approach to follow is to go back a few steps using @xmath72 in @xmath60-space and then transform this to real space yielding @xmath73 where @xmath74 is the lattice constant . if we want to treat a chain end , a modification turns out to be necessary . calculating with eq.([kinetic ] ) the spinon wave function for an undimerized odd size chain one finds that the spinon is strongly repelled from the borders . this is similar to the problem of having a quantum mechanical particle in a finite box where the wave functions are sine - like with nodes at the wall positions . it is due to the fact that the kinetic energy can not be fully satisfied close to a wall since hopping processes through the wall are not possible . however , analyzing the numerical results ( see for instance the uniform susceptibility @xmath75 in fig.9(b ) of ref.@xcite ) one realizes that the probability of finding the spinon at a given site is approximately constant . we interprete here the uniform susceptibility @xmath75 as calculated by laukamp _ et al._@xcite as a measure for the probability to find the spinon . from this observation we conclude that the walls not only truncate the hopping processes but that they induce also another effect . to account for the fact that an approximately constant spinon wave function is the ground state wave function for the problem without dimerization we modify eq.([kinetic ] ) to @xmath76 the effect of this modification is most easily understood by going back to eq.([kinetica ] ) and exchanging the sequence of differentiation and convolution . a constant wave function yields a @xmath15-function due to the initial jump . the resulting term is compensated by the @xmath77 term in eq.([bordercorr ] ) . we like to stress that the modification is motivated only phenomenologically . it would be desirable to have a more microscopic justification as well . equivalently , the discrete version of eq.([bordercorr ] ) reads @xmath78 if the site counting starts with @xmath79 at the border . thus , for chain ends with small frustration we will use eq.([bordercorr ] ) as kinetic energy . for bulk problems we will use eq.([kinetic ] ) . the corresponding discrete versions are eq.([bordercorr_discrete ] ) and eq.([kinetic_discrete ] ) , respectively . thus far , we have discussed only the kinetic energy . considering now the potential energy we first consider regime ( i ) and in particular the majumdar - ghosh point @xmath80 for which the ground state is identical to the short - range rvb state . as noted previously , e.g. @xcite , the confining potential is a linear one in this regime . the potential due to the dimerization @xmath81 with @xmath82 is given by @xmath83 by inspecting fig.[fig - idea](a ) and ( c ) we find that @xmath84 increases by @xmath85 if the spinon hops once , since one singlet on the strong bonds is lost on the right side whereas one singlet is inserted on the weak bonds on the left side . thus we have @xmath86 this allows us to propose the following schrdinger equation for the motion of the spinon @xmath87 with the restriction @xmath88 . the value of @xmath40 can be found from a variational investigation of a single spinon and is found to be approximately @xmath89 @xcite . rescaling eq.([schrod1 ] ) by @xmath90 with @xmath91 ( the site spacing is set to unity ) yields @xmath92 the linear differential equation ( [ schrod1 ] ) with the boundary condition @xmath93 is solved by shifted airy functions @xcite    @xmath94    where the @xmath95 are the zeros of @xmath96 which define by @xmath97 the eigenenergies of the rescaled problem . the normalization is given by @xmath98 . summarizing , using eq.([xi1 ] ) we have obtained    [ result1 ] @xmath99    the form of the wave functions are given in fig.[fig - ai ] where @xmath100 is made equal to unity . a linear potential together with a quadratic kinetic energy leads to energies proportional to @xmath101 if @xmath15 measures the potential strength . such a behavior is actually observed in the dependence of the triplet gap @xmath102 on the dimerization @xmath103 , see fig.[fig - gap ] and @xcite . the corresponding calculation is almost equivalent to the above one since the triplet state can be viewed as two parallel spinons bound together by the same potential . then , the coordinate @xmath33 corresponds to the relative coordinate and the mass @xmath40 has to be replaced by the relative mass @xmath104 since one has two kinetic energy contributions , see e.g. @xcite .    in the light of the success of a simple quantum mechanical picture for the case of large frustration we come back to case ( ii ) of subcritical frustration . before entering the discussion we emphasize that without a gap there is _ a priori _ only a less good justification for a one - particle description since the length scale of the object carrying the @xmath0 is the same as the spatial extent of the bound wave functions . yet it is interesting to see that the simple model works well qualitatively and to a certain extent even quantitatively .    for subcritical frustration , we have two important pieces of information . the kinetic energy is linear in @xmath60 , not quadratic . but the exponent of the ( triplet ) gap growth with dimerization is also 2/3 or close to it @xcite ( @xmath105 ) . to illustrate this point we show in fig.[fig - gap ] gap data for four different values of the frustration . power law fits with exponents close to 2/3 describe fairly well the gap growth in the subcritical and the supercritical frustration regime . the appropriate fit parameters , however , depend on the fit interval chosen . it is known that logarithmic corrections make it in general very difficult to observe the asymptotic behavior . our estimations for the binding length ( see below ) provide additional information on how small the dimerization and how large the system have to be in order for the asymptotic behavior to be reached . thus far , we keep the information that the confining potential can not be linear in the subcritical frustration regime . revisiting the scaling that mapped eq.([schrod1 ] ) to the universal form eq.([schrod11 ] ) shows that the potential must grow as a square root @xmath106 in order to retain the exponent 2/3 . this seems a very plausible conjecture which will be corroborated by a direct calculation of the potential .    in order to evaluate the potential energy between the chain end and the spinon we use ideas of talstra et al . these authors have shown that a spinon state can be generated to 98% overlap by inserting a single spin in a spin chain which is otherwise in its ground state . if we take this idea over to our problem of a spinon and a chain end we may assume that the systems between the chain end and the spinon is in its undimerized ground state which is the ground state of a finite piece of chain with open boundary condition . we have to calculate the expectation value of the dimerization operator @xmath107 ( [ hdimer ] ) in this undimerized ground state in order to obtain an estimate for the potential . a small refinement is actually necessary since we are only interested in the expectation value of @xmath107 with respect to the bulk limit . so we set @xmath108 where the subscript @xmath109 refers to the length of the finite piece of chain and @xmath110 from this equation one sees that the reference to the bulk limit yields only a finite offset .    note that calculating @xmath111 with respect to the undimerized ground state yielding a potential @xmath112 is in the spirit of degenerate first order perturbation theory . the effect of the perturbation , here dimerization , in the subspace of the unperturbed elementary excitations is considered . even though these excitations are not really degenerate they are arbitrarily close in energy so that none of them can be discarded . the calculation of @xmath111 is a perfect task for the dmrg approach @xcite , especially since we deal with open boundary conditions . the results are shown in fig.[fig - pot ] . we use the infinite system algorithm @xcite , keeping 100 states in each step . the error in the energies due to basis truncation is found to be smaller than @xmath113 . ( 8.2,9 ) ( 0,0.1 ) as explained in the main text for three values of @xmath16 . the solid lines are square root fits @xmath114 with the prefactors @xmath115 , @xmath116 , and @xmath117 for @xmath118 , @xmath119 , and @xmath120 . the offset @xmath121 is chosen for the improvement of the fits . the dashed lines are power law fits @xmath122 with the exponents and prefactors @xmath123 , @xmath124 for @xmath118 , @xmath125 , @xmath126 for @xmath127 , and @xmath128 , @xmath129 for @xmath59 . inset : square of the same data for longer chains . [ fig - pot],title="fig:",width=302 ]    one clearly observes a sublinear increase of @xmath111 with the chain length . without frustration ( @xmath118 ) a square root power law fits the data perfectly if a small offset on the @xmath33-axis is taken into account . for larger frustration larger exponents yield better fits . the inset , however , shows that for longer chain lengths ( @xmath130 ) , the square root behavior is recovered for @xmath131 which can be seen from the linear behavior of the squared values . the frustration value @xmath59 is included to show that for this value relatively close to the critical one a linear behavior @xmath132 can not yet be seen . another evidence that for @xmath133 the expectation value @xmath111 rises proportional to @xmath134 stems from the dimension 1/2 of the dimerization operator @xcite . by integration @xmath135 follows the conjectured behaviour . so our conjecture of a square root confining potential for lower frustrations is corroborated by direct calculations . we are now in the position to perform an analysis as in eqs.([schrod1 ] -[result1 ] ) . starting from the generalized schrdinger equation with the border effect corrected kinetic part @xmath136 eq.([bordercorr ] ) @xmath137 + \delta a \sqrt{x } \psi\ ] ] we rescale by @xmath138 with @xmath139 where @xmath140 is the spin wave velocity and @xmath141 some constant which can be deduced from fits to data as in fig.[fig - pot ] . eq.([schrod2 ] ) becomes @xmath142 + \sqrt{y } \psi_i\ ] ] from which the energies are found using @xmath143 the resulting two first wave functions are shown in fig.[fig - wfos ] . the dimensionless energies are @xmath144 and @xmath145 . equation ( [ schrod21 ] ) is solved numerically by equidistant discretization and fast fourier transform with up to @xmath146 points . due to the singular behavior of eq.([schrod21 ] ) the precision is only of about @xmath147 . but this is sufficient for the intuitive description that is being proposed here . note that due to the border correction eq.([bordercorr ] ) the wave functions do not vanish at the chain end , i.e. at @xmath148 . furthermore , the wave functions decay rapidly as @xmath149 , but not following an exponential form , e.g. @xcite , but rather a power law @xmath150 . this is due to the highly non - local character of the kinetic energy eq.([kinetic ] ) . thereby , we do not want to claim that the decay for large @xmath33 follows indeed a power law . the picture we are presenting does not contain the generation of further spinon - antispinon pairs since it is a quantum mechanical , _ single _ particle scenario . it is the pair creation possible for energies larger than the gap @xmath151 which eventually leads to the exponential decay of the spinon probability . our picture works best at intermediate distances as we will show in comparison to the dmrg data in the next section . now we are in the position to compare results of bound spinon calculations directly against numerical results . ( 8.2,12.2 ) ( -2.5,0 ) and @xmath15 . @xmath152 is set to @xmath153 . squares : odd sites ; circles : even sites ; filled symbol : dmrg ; open symbols : discrete bound spinon calculation ; dashed lines continuous bound spinon calculation . [ fig - vgl],title="fig:",width=529 ]    the latter are obtained using the dmrg technique @xcite for open boundary conditions with between @xmath154 to 32 states kept in the iterations . the results are shown in fig.[fig - vgl ] . quantum mechanical bound spinon calculations are carried out in two ways . one is the continuum one relying on eqs.([resodd1],[reseven1],[bordercorr ] ) , the input for the spin wave velocity eq.([spinwave ] ) and the potential parameters for a square root potential ( see caption of fig.[fig - pot ] ) . this means that one takes the result for @xmath155 in fig.[fig - wfos ] and rescales @xmath156 with @xmath100 from eq.([xi2 ] ) . the relevant values are @xmath157 , @xmath158 , and @xmath159 for the uppermost to the lowermost panel in fig.[fig - vgl ] . the results are shown as dashed curves . the other one is a discrete treatment relying on eqs . ( [ resodd2],[reseven2],[bordercorr_discrete ] ) . the results are depicted with open symbols . the values of the potential are those read off fig.[fig - pot ] with @xmath160 . for simplicity the very good power law fits ( dashed lines in fig.[fig - pot ] ) are used for the numerics . the agreement between the dmrg results and the bound spinon model is very good , especially in view of the simplicity of the model . the shape of the curves are very well described by the bound spinon results . in particular , the correct binding length is predicted . this is in particular interesting since we know that the quantum mechanical model does not display an exponential decay . this means that the exponential tail matters only at larger lengths . for intermediate lengths , where the main weight is found , the bound spinon calculations work very well . the main feature that is included by the disrete calculation is the decrease of the amplitude for the first four to five sites . it results from the properties of eqs.([resodd2],[reseven2 ] ) and can thus be attributed to border effects which were neglected in eqs.([resodd1],[reseven1 ] ) but properly taken into account in eqs.([resodd2],[reseven2 ] ) . it is _ not _ due to the dynamics of the problem . the very good agreement of the open - symbol curves and the long - dashed curves for larger distances underlines the applicability of the approximation ( [ norm ] ) .    the overall amplitude is larger than the one predicted by the bound spinon model . since the basis we used included only the shortest range singlets it is not surprising that the true antiferromagnetic correlations are in fact larger . thus we view the bound spinon results in this respect as a lower bound for the antiferromagnetic correlations . but it is interesting that even for the amplitudes the agreement is reasonable for correlation lengths @xmath100 around 15 . for larger values of @xmath100 ( middle panel in fig.[fig - vgl ] ) the agreement becomes not as good . another point that can be addressed easily within the bound spinon model are excitation energies , namely the energy difference @xmath161 ( [ schrod22 ] ) between the ground state @xmath162 and the first excitetd state @xmath163 , see fig.[fig - wfos ] . the continuum ( discrete ) calculation yields @xmath164 ( @xmath165 ) , @xmath166 ( @xmath167 ) , and @xmath168 ( @xmath169 ) for the three cases in fig.[fig - vgl ] in descending order . it is not surprising that the continuum calculation provides lower values since it assumes a potential with smaller exponent , see fig.[fig - pot ] ) , so that the distances between the energy levels become smaller . a smaller exponent implies that the potential increase for smaller arguments is larger than for larger arguments . thus the smaller of two consecutive eigenenergies is lifted with respect to the larger one since the wave function belonging to the smaller eigenenergy is more localized . this effect explains also why the deviation is the largest for @xmath59 where the deviation of the optimum fit exponent to @xmath153 is largest . it was mentioned before that the occurrence of the triplet gap in dimerized systems can be viewed also as a binding phenomenon , see e.g. @xcite . in this case , two spinons bind . however , since only the relative coordinate matters , this problem is almost equivalent to the binding of a single spinon to a chain end . the main difference is that the kinetic energy is doubled since both interaction partners move . so it is understandable that the ratio @xmath170 is fairly constant . the triplet gaps @xmath151 as read off fig.[fig - gap ] are @xmath171 , @xmath172 , and @xmath173 . the ratios of the continuum ( discrete ) values are @xmath174 ( @xmath175),@xmath176 ( @xmath177 ) , and @xmath178 ( @xmath179 ) . similar to the result found previously ( @xmath180 ) @xcite the transition from the ground state to the first excited one is at about half the triplet gap . this implies that it should be observable in scattering experiments as a rather sharp feature within the gap . indeed , raman scattering results have revealed such a feature in zn - doped cugeo@xmath4 @xcite even though the energy is higher than the one - dimensional model predicts . it is clear that at energies @xmath181 pair production becomes possible as it is known in qcd . the gedanken experiment of separating quarks in space in spite of the confinement results in the production of new quark pairs . in our situation any excited bound state can decay by producing a pair of spinons at low values of momentum if its energy fulfills @xmath182 . this mechanism gives rise to a continuum which sets in at @xmath181 . thus higher energy bound states are no longer true eigenstates of the problem but have a finite lifetime since they might decay . the discrete calculation actually shows that @xmath183 and higher are no longer stable . this means that besides @xmath184 , also @xmath185 and @xmath186 should exist as distinct sharp modes . experimentally there is so far no indication for these higher modes . but it should be borne in mind that any deviation from pure @xmath187 behavior tends to decrease the binding energies . bound states are shifted closer to the continuum or even disappear in the continuum due to higher dimensional effects , see also @xcite .    in fig.[fig - mg ] we show the same results as in fig.[fig - vgl ] but for the majumdar - ghosh point @xmath80 . they are based on eqs.([schrod1],[resodd1],[reseven1 ] ) in the continuum case . in the discrete case we used eqs.([potmg],[resodd2],[reseven2 ] ) . the discrete kinetic energy used is nearest neighbor hopping such that the quadratic minimum is the same as for eq.([schrod1 ] ) . one should expect an improved agreement since the kind of basis we are working with is particularly suited for the majumdar - ghosh point @xcite . the agreement , however , is not very good . it becomes decisively better if the bound spinon results are shifted by about three lattice spacing . probably dynamic border effects not taken into consideration by the basic assumptions concerning the kinetic and the potential energy ( focused on the low energy behavior ) are important . the binding length @xmath100 for the majumdar - ghosh model according to eq.([xi1 ] ) is only @xmath188 . thus a low - energy approach might be insufficient for such a small binding length . interestingly , eggert and affleck found in a comparison of continuum field theoretical results with numerical data at @xmath17 also that the agreement is much better if the continuum results are shifted by two sites towards the chain end @xcite . an understanding of either of these observations would help to understand the other . spin - peierls systems are characterized by a coupling of magnetic ( quasi ) one - dimensional spin degrees of freedom to the lattice degrees of freedom , i.e.  vibrating distortions ( for reviews see @xcite ) . at low temperature and zero magnetic field the system is in the gapped dimerized ( d ) phase where the coupling strength @xmath189 of the spin chain alternates from site to site . as the magnetic field is increased the gap becomes smaller and at a critical field @xmath190 the phase changes to an incommensurably modulated ( i ) phase . this i phase can be viewed as a soliton lattice where equally spaced @xmath191-like zeros of the distortions @xmath192 occur . in the vicinity of each zero of the distortion a spinon @xmath0 is localized , see @xcite and references therein .    in fig.[fig - kink ] the soliton lattice is illustrated . these results are found from an adiabatic calculation where the distortions are treated statically . most of the investigations of the i phase assume static distortions . the hamiltonian @xmath193 from which fig.[fig - kink ] is computed by minimization of the ground state energy with respect to @xmath194 reads @xcite    [ hamilton ] @xmath195    where @xmath16 denotes the relative frustration and @xmath196 is the @xmath197 component of the total spin of the @xmath198-site chain . the last two terms in eq.([hamilton ] ) are the zeeman energy and the elastic energy associated to the lattice distortion . a site independent spring constant @xmath199 is used . ( 8.2,6.8 ) ( 0,0 ) , with periodic boundary conditions are studied ( parameters @xmath59 , @xmath200 ) . in the upper panel , results from the self - consistent determined distortions @xmath192 ( circles ) which minimize @xmath201 in ( [ hamilton ] ) are shown . they are depicted between the sites since they belong to bonds . the straight lines are fits to the linear rising distortions with slope @xmath202 . the lower panel displays the resulting local magnetizations ( circles , dmrg calculation ) and the bound spinon result ( solid lines ) . for details see main text . [ fig - kink],title="fig:",width=302 ]    there are regions around site , e.g. , 50 where the distortion alternates as in the d phase . the corresponding local magnetization is essentially zero . close to the linear zeros of the distortions at about sites 25 and 75 one observes a strongly alternating local magnetization which adds up to 1/2 : @xmath203 . the distribution of local magnetizations is observable by nmr @xcite and gives interesting information on the i phase and on the validity of the adiabatic treatment @xcite . an interesting observation is that the ratio between the parallel magnetizations to the antiparallel magnetization at ( roughly ) the same sites is 7:4 as implied by eqs.([resodd1],[reseven1 ] ) .    in this paper we are not interested in the self - consistent minimization of ( [ hamilton ] ) but we take the @xmath204 as given and use the quantum mechanical picture developed above to show that the localized alternating magnetization in the vicinity of a zero of the distortion results from a bound spinon . the zero of the distortion is viewed as a kink defect . since the @xmath192 on the odd ( even ) bonds change sign at the zero there is no simple short - range singlet pattern which ensures singlets at the stronger bonds without any free spin . insofar the situation is comparable to the one at chain ends with weak bonds at the ends . the spinon bound to a kink defect is even simpler than the spinon bound to a chain end since no particular boundary effects on the kinetic energy need to be taken into account . thus we use the continuum description based on eqs.([resodd1],[reseven1],[kinetic ] ) . the derivation of the attractive potential requires an additional assumption in order to treat also site dependent distortions . we find it natural to assume that the chain end potential @xmath205 as determined from fig.[fig - pot ] results from expectation values @xmath206 which decay like a power law with increasing distance from the ends @xmath207 where we took into account that the chain end _ and _ the position of spinon constitute borders . the effect of these borders is treated approximatively as additive , i.e. as independent . a general power law with exponent @xmath208 and prefactor @xmath74 is dealt with . the advantage of eq.([general - pot1 ] ) is that the distortion can be made site dependent @xmath209 . in the vicinity of the kink defect , @xmath210 is a good description where @xmath211 is an appropriately determined slope , see upper panel in fig.[fig - kink ] . a short calculation then leads to @xmath212    thus we obtain again a power law confining potential , the exponent of which is increased by one compared with the site - independent distortion . moreover , the potential does not depend on whether the spinon moves to the right or to the left of the zero .    in fig.[fig - wfds ] a generic result for the localized state at a kink defect is shown . all parameters are set to unity . but rescaling permits , as before , to obtain the general relations    @xmath213    the solid lines in the lower panel in fig.[fig - kink ] depict the convincingly good agreement between the bound spinon picture and the dmrg calculation . the value @xmath202 and the parameters of the dashed line fit to the @xmath59 data in fig.[fig - pot ] is used for the binding length @xmath100 in ( [ xi3 ] ) . in this work we have developed a quantum mechanical picture in the spirit of first order degenerate perturbation theory which describes the confinement of spinons due to dimerization . without dimerization there is no confinement independently of the value of the frustration @xmath16 . the picture works below and above the critical frustration of the undimerized chain @xmath214 , even though the justification below @xmath215 is less good . the main ingredients of the calculation are the known facts concerning the kinetic energy of the spinons which is linear below @xmath215 and quadratic ( plus a mass term ) above @xmath215 . from the dimer expectation values close to the borders of finite chains we deduced a potential @xmath205 which is proportional to the dimerization @xmath15 . it was shown that this potential can be nicely fitted for intermediate distances ( @xmath216 ) by sublinear power laws below and just above @xmath214 . for larger distances ( @xmath217 ) our results indicate that @xmath205 increases like a square root for @xmath58 and linearly for @xmath218 . to our knowledge , the fact that the confining potential below @xmath214 has to be sublinear has not been reported before in the literature . the quantum mechanical picture describes very well the local magnetizations close to chain ends with weak bonds for intermediate @xmath100 ( @xmath219 ) . a very good agreement is also found for the spinon bound to a kink defect . the sum of the moduli of the local magnetizations is ( at least ) enhanced by a factor 11/3 over the value for a single spin . for low @xmath100 values the continuum description becomes less reliable . for large values of @xmath100 the form of the site dependence of the local magnetization is still captured by our quantum mechanical picture . in particular , binding lengths can be easily estimated using eqs.([xi1],[xi2],[xi3 ] ) . the overall amplitude of the staggered component , however , is in fact larger for @xmath220 at @xmath221 @xcite .    besides local magnetizations the eigenenergies of the bound states can be estimated easily within our approach . we find that in truly one - dimensional systems transitions to the first three excited states at chain ends should be possible . any degree of higher dimensionality , however , will reduce this number . the formation of gapful triplets can equally be viewed as binding of two spinons @xcite . the essential differences to the case of chain ends are that the site variable denotes the relative coordinate between the spinons and that the kinetic energy is doubled since both interaction partners move . the quantum mechanical picture is scale invariant . this means that for power law potentials the wave functions and energies can all be scaled to one universal case , see figs.[fig - ai ] , [ fig - wfos ] , and [ fig - wfds ] . thus the quantum mechanical description provides a pedestrian approach to perturbations with anomalous dimensions and exponents . the anomalous exponents result naturally from the scaling behavior of the potential and the kinetic energy under changes of the length scale . a dimensional analysis implies already that for @xmath222 and @xmath223 the characteristic lengths scale like @xmath224 and the characteristic excitation energies like @xmath225 . we like to point out that the basis states of our approach can be refined to take longer range antiferromagnetic correlations better into account . to this end longer range singlet pairs must be included and superposed . dispersions and potentials can then be found by concepts similar to those proposed by sutherland @xcite for short range singlets in @xmath226 . ( 8.2,1.5 ) ( 0,0.3 )    as an outlook to higher dimensions we present in fig.[fig - ladder ] the leading sketch of the singlet distribution on a ladder with strong bonds on the rungs . ladders are objects of ample investigations . the alternating magnetizations @xcite and enhanced antiferromagnetic correlations @xcite indicate the similarity to open chains .    in fig.[fig - ladder ] the situation is depicted of the spinon ( arrow ) generated by the insertion of the defect ( circle ) after three hops on the upper leg . one clearly sees that the motion of such a spinon is also confined by a monotonic increasing potential . an increasing distance of the spinon from its original rung induces an increasing number of misplaced singlets . this leads to an increase in energy . from fig.[fig - ladder ] it is clear that the basic concepts of spinon motion and confining potential can be extended also to other rvb - type spin systems . the results discussed in this paper are not restricted to chains only . we thank g. els , p.h.m . van loosdrecht , g. martins , and e. mller - hartmann for helpful discussions . one of us ( gsu ) acknowledges the kind hospitality of the nhmfl , tallahassee , where a large part of this work was carried out . we acknowledge financial support by the dfg through the sfb 341 and by the nsf grant dmr-9520776 .
Abstract:

Answer: a quantum mechanical picture is presented to describe the behavior of confined spinons in a variety of @xmath0 chains . the confinement is due to dimerization and frustration and it manifests itself as a nonlinear potential @xmath1 , centered at chain ends ( @xmath2 ) or produced by modulation kinks ( @xmath3 ) . the calculation extends to weak or zero frustration some previous ideas valid for spinons in strongly frustrated spin chains . the local magnetization patterns of the confined spinons are calculated . a ( minimum ) enhancement of the local moments of about 11/3 over a single @xmath0 is found . estimates for excitation energies and binding lengths are obtained .