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: distance determinations to galaxies are key to many areas of study in extragalactic astronomy , but remain problematic . the potential for systematic errors lingers above much of the work , particularly at the accuracy level of a percent or less currently desired ( see review by * ? ? ? these concerns can be addressed through cross calibration among techniques that have independent challenges . for example , one estimator might have a hidden dependence on metallicity , another might be more susceptible to internal extinction . one state - of - the - art recent study compares distances obtained using sn ia and cepheid variable stars in the same galaxies @xcite to claim a distance calibration precision of a few percent . such cases , where one has distance measurements from a variety of methods for a specific galaxy , are rare for two reasons . first , because distance estimators are often related to distinct stellar populations , they are not often found within the same stellar system . the classic example is the dichotomy between the galactic environments of rr lyrae and cepheids @xcite . second , and currently most relevant , is that several factors conspire to limit estimators to distinct distance regimes . cepheid observations , like those of rr lyrae , eclipsing binaries and the tip of the red giant branch ( trgb ) , require observations that resolve the stellar populations of the system and are therefore limited to relatively nearby galaxies , even when using the high angular resolution provided by the _ hubble space telescope_. on the other hand , sn ia , due to their rarity , tend to be found in more distant galaxies simply because of the gain in search volume with increasing distance . these effects make it quite difficult to compare the results among certain classes of distance estimators . the @xcite study cited above , which is among the most extensive and complete , is based on only eight galaxies with both cepheid and sn ia measurements . one solution to this limitation relies on identifying a distance estimator that can be applied to _ all _ galaxies , near and far , early or late type , giant or dwarf . even if such an estimator has a large intrinsic distance uncertainty for any individual galaxy , and so is not ultimately superior to the other estimators , the gains reaped by being able to compare the relative precision and accuracy of the preferred estimators over many tens , if not hundreds of galaxies , is critical . distance determinations from scaling relations , such as tully - fisher @xcite or fundamental plane @xcite , span the range of distances probed by many methods and so provide a fiducial against which different distance estimators can be compared . their shortcoming is that they are each applicable to different , narrow classes of galaxies , limited in type and luminosity , and therefore do not resolve the entire problem .    in a series of papers , we have presented a new scaling relation that is applicable to _ all _ galaxies , regardless of luminosity or morphological class @xcite . this relation , termed the fundamental manifold ( fm ) in reference to its most direct antecedent , the fundamental plane , presents an opportunity to uncover hidden systematic differences among the current set of distance estimators . here we use existing data , distances from the ned 1-d database and photometry and kinematics compiled within the hyperleda database @xcite , to test this basic idea . neither the data , due to their inhomogeneity , nor the sample , due to the intractable selection criteria , are optimal for precise cross - calibration of distance estimators and eventual application to related problems , such as the measurement of h@xmath2 . therefore , the results presented here should be viewed as only illustrative of what an optimized treatment using the fm could yield . in  2 we discuss our selection of the sample , in  3 we present the results of this comparison among distance estimators , illustrate how the fm can be used to trace potential internal systematic dependencies , and reverse the calibration and use the fm to obtain distances and estimate @xmath3 . in  4 we summarize our results . we utilize the ned 1-d database , which is a compilation of published distances from a wide ranging set of methods . without prejudice for or against any of the entries in the database , we calculate the weighted average ( weighting inversely by the quoted uncertainty and neglecting the measurement if no uncertainty is given ) for each of the distance estimators , for each galaxy . we require at least two measurements using a particular method to calculate a mean distance measurement for a given galaxy . there are 1153 galaxies for which we retrieve a distance using at least one of the methods we are comparing . in figure [ fig : distances ] we show how distance estimators are often effectively constrained to distinct distance ranges . more critically , the distribution of objects with sn ia distance measurements overlaps very little with the other methods , limiting the number of systems that can be used to directly check for potential systematic errors in the cosmologically - critical ia distances .      to obtain the structural parameters necessary to place a galaxy on the fm , we search the hyperleda database @xcite , which provides morphologies , kinematic measurements , magnitudes and structural parameters , the latter two based on the work of @xcite , for those galaxies for which we obtained distances from ned 1-d . for the rotation measurements , we adopt the homogenized values provided in hyperleda , which include an inclination correction . we calculate weighted averages of both the stellar rotation and velocity dispersion measurements where available , and only accept values only if there are at least two measurements of either quantity , and if either the rotation or dispersion is larger than 10 km sec@xmath4 , to avoid the systems with the largest fractional uncertainties . for the morphologies , we average t - types using uniform weights and require a minimum of only one measurement . the photometric parameters are taken from the archival b - band data , which is the most prevalent . there are various reasons ( star formation , extinction ) to use redder - band photometry , and we will discuss results using the available i - band data . in the future , we will use the s@xmath5 g database @xcite to obtain photometric parameters at 3.6@xmath6 m , which appear to significantly reduce the scatter in scaling relations @xcite . the distribution of morphologies , for the subsample of galaxies for which morphologies are catalogued , is shown in figure [ fig : morphs ] . again , we find that certain distance estimators are limited in their coverage , making comparison on an object - by - object basis difficult . there are 343 galaxies with both distance measurements and the necessary data to place them on the fm .    for a subset of the data , we also obtain an optical color , @xmath7 , from the rc3 catalog @xcite to investigate the role of color , and , by inference , of differences in the stellar mass - to - light ratio , on a galaxy s deviation from the fm . we use observed magnitudes rather than k - corrected rest - frame magnitudes , because , despite the large range in distance of the ned 1-d sample , those galaxies with _ all _ of the necessary data lie at @xmath8 6000 km sec@xmath4 . one evident avenue for progress is obtaining these data for a larger fraction of the galaxies with sn ia distances , and thereby also extending the sample to greater distances . our compilation of the available data is presented in table [ tab : data ] ( complete table available electronically ) . lrrrrrrrrrrrrrrrr ngc7814&11.38&0.83&22.23&118.1&230.7&0.0&172.9& ... & ... & ... &12.70& ... & ... & ... & ... &1047 + ngc0055&8.68&0.54&22.18&400.1&58.9&0.0&0.0&1.92&2.10& ... & ... & ... & ... & ... & ... &142 + ngc0147&10.65&0.78&24.02&376.8&0.0&0.0&26.4& ... &0.64& ... & ... & ... &0.63& ... & ... &@xmath9176 + ngc0150&11.90&0.50&21.40&63.4&164.9&0.0&0.0& ... & ... & ... & ... & ... & ... & ... & ... &1566 + ngc0185&10.38&0.73&22.11&177.1&0.0&0.0&23.3& ... &0.61& ... & ... & ... & ... & ... & ... &@xmath9201 + ngc0205&9.00&0.82&21.82&292.5&0.0&22.0&33.5& ... &0.80& ... & ... & ... & ... & ... & ... &@xmath9258 + ngc0221&9.19&0.88&18.75&65.2&0.0&37.8&80.9& ... &0.69& ... &0.77& ... &0.82& ... & ... &@xmath9203 + ngc0224&4.88&0.68&19.89&802.0&0.0&75.1&183.2&0.80&0.79& ... &0.75&0.72&0.77& ... &0.75&@xmath9294 + ngc0247&9.42&0.54&23.37&492.2&249.8&0.0&0.0&3.76&3.60& ... & ... & ... & ... & ... & ... &142 + ngc0300&8.51&0.58&22.87&594.5&85.1&0.0&0.0&1.99&1.96& ... & ... & ... & ... & ... & ... &150 + [ tab : data ] the fm has been described in detail elsewhere @xcite , but , in summary , it is derived from the virial theorem with the additional information that the mass - to - light ratio within @xmath10 , @xmath11 , can be either 1 ) measured independently ( e.g. , via strong gravitational lensing ) or 2 ) estimated using an empirical function of the internal kinematics , @xmath12 , and the surface brightness within @xmath10 , @xmath13 . in the latter case , we have found that the structural differences ( spatial and kinematic ) among galaxies , other than those captured in the variations in @xmath11 , are insignificant given the observational errors ( @xmath14 dex ) . the fm is a conceptually simpler scaling relation than either tully - fisher ( tf ) or fundamental plane ( fp ) in that it can be expressed independent of wavelength and has only a zero point term to calibrate ( rather than slope in tf or  tilt " in fp ) : @xmath15 where @xmath12 is expressed as the combination of pressure ( @xmath16 ) and rotational support ( @xmath17 ) , @xmath18 ( @xmath19 is expected to lie between 2 and 3 ; * ? ? ? * ; * ? ? ? * ) ) , and @xmath20 is a constant that is effectively the zero point of the fm distance estimator . with measurements ( or estimates ) of all terms on the right hand side of the equation , one solves for @xmath10 , which , in concert with the angular measurement of @xmath10 , provides the distance . because we generally lack independent estimates of @xmath11 , we use the empirical relationship for @xmath21 derived from existing data @xcite . use of the empirical relationship reintroduces wavelength dependences and additional complexity ( analogous to fp  tilt " ) , but the fm still retains the advantage over previous scaling relations in that it is applicable to all galaxies .    as with any empirical relationship , the relationship changes slightly depending on the data used to define it , and uncertainties in the fitting function ultimately translate to systematic uncertainties in the distances estimated using the fm . consequently , comparing fm distances to those obtained using other distance estimators provides a test of both the fm and the independent distance estimators . comparisons of the fm and multiple independent methods will allow us to distinguish between systematic problems in the fm versus problems with the inditvidual independent estimators . here we choose to define the fitting function using only the spheroidal galaxies in the original data @xcite , which are expected to have minimal stellar mass - to - light ratio variations . we exclude globular clusters and ultracompact dwarfs , whose nature is still controversial ( see * ? ? ? * ) and that extend the fm beyond the @xmath10 range needed for this work . although changes in the fitting function for @xmath11 will affect the results , the fitting function is most constrained in exactly the parameter range most relevant here  that of disk and spheroid galaxies with luminosities around @xmath22 . the functional form we use is @xmath23 this equation provides an estimate of @xmath11 in the i - band . there are two related , but distinct , questions that will be addressed in the following section : how much does the use of b - band magnitudes increase the scatter ( there would be no increase if all galaxies had the same @xmath24 color ) and how much scatter is introduced in the @xmath25-band @xmath11 by differences in stellar populations ? to test each distance method , we use the distances provided by ned 1-d to calculate the physical value of @xmath10 and then place the galaxy on the fm . in our original work , from which we take our fitting function for @xmath11 , we used simple hubble flow distances with an adopted hubble parameter . any difference between the true hubble parameter and that adopted previously will affect the calibration of @xmath11 and @xmath20 , manifesting itself as a zero point shift ( _ i.e. _ , as a change in @xmath20 ) . we therefore have freedom to adjust the zero point of the fm as needed ( see below for our treatment of this ) and so our discussion here focuses on relative differences between distance estimators rather than on an absolute calibration . one consideration before proceeding with our comparison of distance estimators and the fm is that of the appropriate choice of @xmath19 in the definition of @xmath12 . in previous work we adopted @xmath26 , as also done by others ( e.g. , * ? ? ? * ; * ? ? ? however , reasonable models suggest that values between 2 and 3 are plausible ( see * ? ? ? * ) . because @xmath27 maps onto morphological type , and morphological type maps onto color , which is presumably related to the stellar mass - to - light ratio , there is the potential to misidentify a relationship between the fm residual ( the deviation of a galaxy from the fm ) and @xmath19 , with one that is truly due to color ( or vice - versa ) . without the proper choice of @xmath19 and some means of avoiding differences in stellar mass - to - light ratios beyond those included in the empirical calibration of @xmath11 , there will be systematic differences when comparing distances to early and late type galaxies . we will mitigate this difficulty by 1 ) proceeding to optimize our choice of @xmath19 , 2 ) removing any residual correlation between fm residuals and color , and 3 ) providing comparisons among galaxies of similar morphological type . there is indeed a strong correlation between @xmath28 and @xmath27 ( a spearman rank correlation test for our sample puts the probability of it arising at random at @xmath29 ) . as such , for the purpose of constraining @xmath19 we work with a subsample for which @xmath27 is independent of color ( figure [ fig : colsel ] ) ; @xmath30 and @xmath31 ( the probability of the correlation arising randomly for this sample is 0.79 ) . this particular part of the analysis , the derivation of @xmath19 , is restricted to early type galaxies because @xmath16 is generally not available for disk galaxies . as a result , the galaxies in figure [ fig : colsel ] are naturally restricted to fairly red colors . if a hidden correlation remains , we expect some of that to influence the best fit value of @xmath19 , again because @xmath32 is expected to correlate with color , and eventually to lead potentially to differences between fm residuals and morphological type .    for this subsample , we calculate the root mean square ( rms ) residual about the fm ( allowing the zero point to float ) and identify the value of @xmath19 that minimizes the rms ( @xmath33 ; see figure [ fig : minalpha ] ) . reassuringly , this value lies within the theoretically plausible range and we adopt this value of @xmath19 for the remainder of this study . however , this calculation is unfortunately quite sensitive to the selected subsample . for example , if we select galaxies that have @xmath34 ( instead of 0.85 ) , then the minimum rms residual occurs at @xmath35 ( although a correlation between color and @xmath27 still exists at @xmath36 confidence for this sample and so this sample is not preferred ) . the robustness of the @xmath19 determination needs to be improved both by obtaining a larger sample of galaxies and by performing this analysis with photometry in a filter band that is less sensitive to morphological type ( and , by relation , to color ) . the effect of varying @xmath19 is principally to alter the relative offset between early and late type galaxies relative to the fm ( see  4 ) , but does not affect the conclusions presented in  5 . any remaining systematic behavior of the fm residuals will point to new phenomenon , including systematic errors in distance estimators . we begin by examining the dependence of fm residual to color ( figure [ fig : color ] ) . the data show a trend of residual with color . it is not surprising that the early - types have near zero mean residual , because the fm relationship we use was defined using spheroidal galaxies ( from a different sample based on a range of sources , see @xcite , but early - types nonetheless ) . neither the sense nor the magnitude of the trend is necessarily straightforward to explain . there is an expectation of a trend with color because the stellar mass - to - light ratio will change with color . however , other factors change as well . for example , the gas mass fraction of a galaxy correlates with color and that mass is unaccounted for in the empirical derivation of @xmath11 because we used only spheroidal systems . we use the mean relationship , shown as the line , to correct all of our data ( and so subsequently restrict ourselves to systems with available color measurements ) . if a particular estimator ( including the fm ) has a systematic error that depends on galaxy type ( color ) , it will contaminate this plot . as with many of the issues raised , we expect this problem to be mitigated using redder photometry , but proceed nonetheless with the data currently available . recall that our expression for @xmath11 provides an estimate of the total @xmath37band mass - to - light ratio within @xmath10 , even though we are using @xmath38band magnitudes . this means that there are two potential sources of scatter introduced by stellar population differences . first , there is the possibility of large variations in @xmath24 within the sample that then introduce scatter because we implicitly assume a constant @xmath24 color ( solar ) when converting @xmath39-band luminosity into the corresponding mass in solar masses . second , scatter can arise from stellar population differences and the associated variations in the stellar mass - to - light ratio . we show , in figure [ fig : bmi ] , the distribution of @xmath24 colors for our galaxy sample , and note that a large fraction of the galaxies are fairly well constrained to colors that vary by about 0.3 mag ( a constant difference from solar color is not a cause for concern because we eventually renormalize the entire relation using the sbf distances ) . more importantly , this range of colors produces small scatter ( 0.046 dex ) in @xmath11 between the values calculated using the @xmath39-band magnitudes ( and the color correction described above ) and the @xmath25-band magnitudes ( with a color correction derived for the @xmath25-band ) . we conclude that variations in @xmath24 colors , beyond those that we remove using a mean trend of fm residual with @xmath28 , are a minor contributor to the scatter we observe in the fm ( @xmath40 dex for these data ) . our @xmath28 color correction is intended to empirically remove such a correlation , although as hinted at previously , this correction may be removing other factors as well that work in parallel to stellar population differences . in figure [ fig : color ] we plot the expected variation in stellar mass - to - light ratio vs. @xmath41 from @xcite , selecting one particular model ( closed box , scaled salpeter imf ) , demonstrating that the effect we are removing is indeed of the same scale as that expected from stellar population variations . finally , we estimate the degree to which going to even redder passbands might reduce the scatter , when one does not have a color measurement . for the same stellar population model as in figure [ fig : color ] , the slope in the logarithmic relation between @xmath41 and the stellar mass - to - light ratio at @xmath42 is about half that in @xmath25 . therefore , if the scatter in the fm is dominated by stellar population differences , and we do not have colors , we still expect to reduce the scatter by two simply by observing at @xmath42 . further gains are possible if extinction is important @xcite . if other causes , not related to the photometric passband , contribute comparably to the scatter , or if we can partially account for stellar population differences , then we expect less than a factor of two improvement by going into the infrared . for example , the scatter in our @xmath39-band fm is the same as in our @xmath25-band fm ( @xmath14 dex ) , because we have accounted for the stellar population differences at @xmath39 . the key to interpreting our results is to determine whether multiple estimators show the same systematic behavior in fm residuals , in which case the problem is most likely to lie with the fm or the various corrections we have implemented , or if the fm residuals are peculiar for specific estimators , in which case the fault lies with that particular estimator . on the other hand , if none of the estimators exhibit systematic problems relative to the fm , then we will conclude that the estimators and the fm are accurate distance estimators to at least within the quoted uncertainties . our aim is to test a variety of popular distance estimators to determine which , if any , may lead to distance estimates with systematic errors . using the fm as a reference , we now compare estimators that we could not previously test against each other . to do this , we plot the galaxies on the fm using the mean distances drawn from the ned 1-d data ( figure [ fig : comp1 ] ) . normalization or slope errors for any particular estimator will suggest that there is a problem with that estimator . on the other hand , if the normalization or slope errors appear endemic , then the problem will lie with the fm .    for each distance estimator with more than 10 galaxies with all the necessary data , we present the results , showing the pressure supported systems , which have velocity dispersions and sometimes also have stellar rotation values , @xmath17(upper panels in the figure ) , and the purely rotationally supported galaxies , which have no quoted stellar velocity dispersion and for which we adopt the width as a measure of @xmath17 ( lower panels ) . in general , if there is a stellar velocity dispersion measurement , the galaxy is a pressure supported systems ( @xmath43 ) . within each panel we also color - code based on morphology , dividing the early and late type galaxies at a t - type of 1 . we have set the zero point of the fm , @xmath20 , using the results from the pressure - supported , sbf sample , which is our largest subsample of galaxies and also predominantly consists of early - type galaxies , which will be less susceptible to extinction and stellar population variations . a cursory examination of the panels shows that there are no disastrous problems with any of the distance estimators as applied to any of the galaxy subsamples , even though some perform better than others ( either in terms of zero point or scatter ) . there are a few individual galaxy outliers , although it is often the same galaxy that is an outlier in multiple panels because distances are available from multiple estimators , suggesting that fault lies not with the distance estimate but rather with one of the other parameters that enters the fm . ngc 4704 , which has an unusually large residual in the sn ia panel , is removed from subsequent analysis of the sn ia distance estimator . in figure [ fig : comp2 ] we show the results for galaxies that have distances measured using eclipsing binaries , rr lyrae , and masers .    in figure [ fig : comp3 ] we compare the mean residual and the uncertainty in the mean between the various distance methods , for pressure and rotationally supported galaxies ( the values are given in table [ tab : comp ] ) . the agreement is strikingly good , with the possible exception of results obtained using the  other " methods , although there the uncertainty is large because of the small number of galaxies in that category . none of the differences are statistically significant , but the magnitude of differences is in the range that is potentially important for distance estimators ( @xmath1 few to tens of % ) . for example , the difference in zero points obtained using the trgb and sn ia methods and pressure supported galaxies is 0.062 , which corresponds to an inferred distance difference of 15% . unfortunately , for the purpose of identifying if there are any real discrepancies at the target level of a few percent , the error bars are still somewhat too large due to the large amount of scatter in the fm ( arising from the heterogeneous b - band data ) and small sample size ( as many galaxies lack the full set of measurements needed to calculate the fm ) . lrrl    cepheids & @xmath44 & @xmath45 + sbf & @xmath46 & @xmath47 + sn ia & @xmath48 & @xmath49 + trgb & @xmath50 & @xmath51 + pne & @xmath52 & @xmath53 + [ tab : comp ]    the uncertainties in our comparison of distance estimators are currently dominated by the number of galaxies in each class and the intrinsic scatter in the fm . the small sample sizes can be addressed relatively easily . many galaxies with extant distance measurements , particularly those with sn ia distances , lack the additional measurements ( typically @xmath12 ) necessary to place them on the fm . these can be obtained with ground - based observing . the fm scatter can be improved by going to redder pass - bands . if all of the data had uncertainties similar to the pressure - supported , sbf measurements , we would easily confirm or refute the level of discrepancy among distance estimators found here . the fm is not only of use in uncovering systematic differences among distance estimators , but is also potentially useful in uncovering systematic errors in any particular distance estimators within a set of galaxies . one long - running concern is the role chemical abundance might play in affecting distance estimators ( see * ? ? ? * ; * ? ? ? * for discussions of the effects on cepheids and sne ia , respectively ) . given the correlation between host galaxy gas - phase metallicity and luminosity @xcite , we use our sample to examine whether there are potential dependencies on metallicity by measuring correlations between fm residuals and luminosity ( a correlation could also implicate another parameter that correlates with luminosity , e.g. , mean stellar age ) . in figure [ fig : abundances ] we show the fm residuals as a function of galaxy absolute magnitude for our five primary estimators . visually , two of the panels appear to show a trend of increasing fm residual with decreasing luminosity ( sbf and sn ia ) . however , this similarity among panels overstates the concordance because some of the galaxies , particularly at the bright end , appear in multiple panels . in table [ tab : metallicity ] , we present measurements of the relationship between the fm residuals and absolute magnitude for each distance estimator . we list both the fit obtained using all of the galaxies and after trimming the poorly populated faint end of the galaxy distribution ( either at @xmath54 or @xmath55 ) . the only two significant correlations are provided by the sbf and sn ia data ( their rank correlation probabilities of arising at random are both @xmath56 , if we confine the sne sample to galaxies with @xmath57 ) . lrrl    cepheids & @xmath58 & @xmath59 + sbf & @xmath60 & @xmath61 + sn ia & @xmath62 & @xmath63 + trgb & @xmath64 & @xmath65 + pne & @xmath66 & @xmath67 + [ tab : metallicity ]    if we attribute the slopes to a metallicity dependence in the estimator , we can quantify the effect and compare to previous determinations . our best fit for the sn ia in galaxies with @xmath57 implies a distance offset of @xmath68 dex/[o / h ] ( for the 4.5 mag/[o / h ] relation of @xcite ) , which agrees with the sense found by @xcite and @xcite in that one would underestimate the distance to sne in more metal rich galaxies , because sne are intrinsically more luminous after light curve correction . it is difficult to quantitatively compare our results to those of @xcite for various reasons , including their use of stellar masses vs. our use of luminosities and our calibration to gas - phase abundances . however , they find distance modulus variations that correspond to about 0.15 over the 1 order of magnitude in mass . if we assume that the order of magnitude in mass correspond to the same magnitude range in luminosity , then we would predict a distance effect of 0.23 dex over the corresponding range , which in turn results in a change of @xmath69 in distance modulus . this is a much larger value than that found by @xcite , but the discrepancy is not significantly dramatic given our large uncertainties . likewise , our value expressed in dex/[o / h ] is larger than that found by @xcite ( @xmath70 ) , but again the discrepancy is @xmath71 . for cepheids , we find no significant correlation , but our measurement ( @xmath72 dex / mag or @xmath73 dex/[o / h ] if the correlation arises from metallicity ) can be compared to previous determinations . given the large uncertainties , this result agrees with previous values , which generally range between @xmath74 and @xmath75 dex/[o / h ] @xcite .      in figure [ fig : fm ] we show the fm for galaxies that have their distances estimated using multiple techniques ( we use the average for the adopted distance ) . the relation is satisfied by galaxies that range widely over luminosity , morphology , and kinematics . it is tempting to invert the arguments above and use the fm not as a fiducial against which to test various distance estimators , but as a distance estimator itself . as stressed in the introduction , these particular data are suboptimal for such a experiment , but we proceed nonetheless .    using the fm now for all of the galaxies in the sample , we estimate a distance to each , and show the relationship between recessional velocity and distance in figure [ fig : hubble ] . we define the fm zero point by setting the mean fm residual among estimators to zero , but define that mean in several different ways . first , we use the sbf calibration . second , we calculate the weighted mean of the residuals in table 2 , excluding our calibration sbf sample . third , we use the zero point from the sn ia for pressure supported galaxies . the results differ by at most 4 km sec@xmath4 mpc@xmath4 in the inferred @xmath3 . we therefore cite a systematic uncertainty of 4 km sec@xmath4 mpc@xmath4 in our @xmath3 measurement below . we choose to plot the results using the intermediate value of the fm zero point derived from the three methods ( the mean of the residuals of the estimators ) . the inset shows the result of fitting the relationship for galaxies with @xmath76 km sec@xmath4 , @xmath77 outliers excluded , binned by 10 , and forcing the fit through the origin . we exclude the low @xmath78 region to minimize the effect of local flows . the resulting slope implies @xmath79 ( random ) @xmath80 4 ( systematic ) km sec@xmath4 mpc@xmath4 . the estimate of the systematic uncertainty does not include a variety of potential problems that are ignored here ( modeling of bulk flows , internal extinction corrections , adjustment for potential biases in the galaxy sample , etc . ) . therefore , while we treat this result cautiously , we present it to show 1 ) the level of precision possible with plausible sample sizes and 2 ) to illustrate how one could determine @xmath3 in a way that utilizes as many ( or as few ) different distance estimators as desired .    one of those sources of uncertainty , internal extinction , can be mitigated by going to the infrared , as shown most recently by @xcite for _ spitzer _ wavelengths where the scatter about the tully - fisher relation is reduced from 0.43 in the extinction - corrected @xmath38band data to 0.31 in the non - extinction - corrected 3.6@xmath6 m data . we expect analogous gains in the fm for galaxies in which internal extinction is important . depending on the fraction of galaxies in a sample for which extinction contributes significantly to the observed fm scatter , proportional gains in the determination of @xmath3 will follow . we demonstrate that the fundamental manifold ( fm ; * ? ? ? * ) provides a fiducial against which any set of distance estimators can be compared . many , if not most , galaxies will only ever have distances measured with one method , making it impossible to use them to test distance estimators directly . the fm fiducial , a benchmark for all distance estimators , constrains whether systematic errors exist among estimators and whether they are creeping into even a single estimator as a function of external parameters , such as metallicity . large sample size is critical for examining selection effects where , for example , one might want to split the sample in various ways ( _ e.g. , _ with environment , luminosity , morphological type , star formation rate ) .    among the distance estimators we examine ( sn ia , cepheids , surface brightness fluctuations , the luminosity of the tip of the red giant branch , circumnuclear masers , eclipsing binaries , rr lyrae stars , and the planetary nebulae luminosity functions ) , we find no statistically significant differences ( @xmath81 ) , but caution that differences that are unacceptably large for state - of - the - art distance uncertainties ( @xmath1 few % ) are still allowed by our data . the lack of precision with which we could discriminate among methods is rooted in sample size , which in certain cases is still too small , and in the less - than - optimal data available ( @xmath39 band and heterogeneous ) . both of these problems will be resolved with a feasible dedication of resources . we illustrate the use of the fm distances to 1 ) test the metallicity sensitivity of different estimators , confirming the sn ia distance dependence on host galaxy metallicity and 2 ) provide an alternative calibration of @xmath3 that replaces the classical ladder approach in the use of extragalactic distance estimators with one that utilizes data over a wide range of distances simultaneously . we thank barry madore and ian steer for compiling and generously providing the ned 1-d distance database without which this project could not have been completed . dz acknowledges financial support for this work from a nasa ltsa award nng05ge82 g and nsf grant ast-0307482 . aiz acknowledges financial support from nasa ltsa award nag5 - 11108 and from nsf grant ast-0206084 . dz and aiz thank the institute of astronomy at cambridge university and the center for cosmology and particle physics at new york university for their hospitality during the completion of this study . we acknowledge the usage of the hyperleda database ( http://leda.univ-lyon1.fr ) and thank all those involved for producing this highly useful resource . this research has made use of the nasa / ipac extragalactic database ( ned ) which is operated by the jet propulsion laboratory , california institute of technology , under contract with the national aeronautics and space administration . we also extend our thanks to those involved in creating and maintaining this excellent resource .
Abstract:

Answer: we demonstrate how the fundamental manifold ( fm ) can be used to cross - calibrate distance estimators even when those  standard candles " are not found in the same galaxy . such an approach greatly increases the number of distance measurements that can be utilized to check for systematic distance errors and the types of estimators that can be compared . here we compare distances obtained using sn ia , cepheids , surface brightness fluctuations , the luminosity of the tip of the red giant branch , circumnuclear masers , eclipsing binaries , rr lyrae stars , and the planetary nebulae luminosity functions . we find no significant discrepancies ( differences are @xmath0 ) between distance methods , although differences at the @xmath1 10% level can not yet be ruled out . the potential exists for significant refinement because the data used here are heterogeneous b - band magnitudes that will soon be supplanted by homogeneous , near - ir magnitudes . we illustrate the use of fm distances to 1 ) revisit the question of the metallicity sensitivity of various estimators , confirming the dependence of sn ia distances on host galaxy metallicity , and 2 ) provide an alternative calibration of h@xmath2 that replaces the classical ladder approach in the use of extragalactic distance estimators with one that utilizes data over a wide range of distances simultaneously .


Instruction: Write an abstract for the following acadamic paper.

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

Paper: it is now possible to realize ideal quantum measurements on individual quantum objects , for instance atoms , ions and photons . beyond being a tool to monitor the system s evolution , projective quantum measurements can be used to prepare it in specific quantum states . for instance , an ideal quantum non - demolition ( qnd ) measurement of a field s photon number projects it onto a photon - number ( fock ) state @xcite . however , due to the basic quantum indetermination of the measurement outcome , measurement - induced state generation is not deterministic . quantum feedback control techniques @xcite make it possible to overcome this limitation and to produce quantum states on demand . these techniques generally combine weak quantum measurements with a real - time correction of the system s state depending on the classical information extracted from the measurements . beyond preparation of specific states , these feedback schemes can also protect them from decoherence , resulting from the coupling of the system with its environment .    in this paper , we propose a quantum feedback scheme for the on - demand preparation of fock states stored in a high - quality superconducting microwave cavity and for their protection against decoherence . this scheme is designed to operate with an existing cavity - qed set - up @xcite . the feedback loop uses three steps . we first extract information on the photon number distribution with a single circular rydberg atom by a qnd process . this atom interacts dispersively with the cavity field and does not exchange energy with it , but experiences a light shift proportional to the photon number . the measurement of this shift , with the help of a ramsey interferometer , provides information on the field intensity . it modifies the field state accordingly through the quantum projection . this information is , in the second step , used to estimate the new cavity field state through a quantum filtering process @xcite . in the third step , we correct the field state , using a coherent field pulse injected in the cavity . the feedback law used to calculate the amplitude of this pulse is chosen by adapting to this discrete situation the lyapunov - based techniques proposed in ref . these three steps represent the three basic components of any feedback loop : a sensor , a controller and an actuator . in contrast to classical closed - loop systems , we use here a quantum sensor . however , both controller and actuator remain classical . by iterating the feedback loop , we steer the cavity towards any target fock state . we also efficiently stabilize it against cavity decay . recently , a similar feedback scheme to generate fock states of light in an optical cavity has been proposed @xcite . it is based on a continuous monitoring of the mean number of photons . instead , we propose here to use a discrete qnd measurement followed by a quantum filtering process providing complete information on the field s density matrix . this precise knowledge of the field s state gives us a much better insight into the feedback action . it could also be an asset in other feedback schemes requiring , for instance , information on the field s phase in addition to its intensity . the paper is organized as follows . in sec . [ sec : setup ] , we describe the elements of the set - up and give realistic experimental parameters . we present the quantum - mechanical operators describing the evolution of the cavity state under measurement , decoherence and pulse injection . section  [ sec : ideal ] is devoted to a detailed analysis of the quantum filter and of the feedback law in an ideal situation . we first present the main elements of the feedback loop and describe the tuning of the controller gain . we give qualitative arguments showing that the proposed strategy is stable . the detailed mathematical proof of convergence and stability , relying on stochastic lyapunov techniques will be published elsewhere @xcite . we finally present quantum monte carlo simulations of closed - loop trajectories of the cavity field exhibiting the feedback performances . in section  [ sec : real ] , we take into account the known experimental imperfections of the existing set - up . we modify accordingly the feedback algorithm and present extensive simulations of its operation . we conclude in sec .  [ sec : conclusion ] . :  high-@xmath0 microwave cavity , @xmath1 :  box producing rydberg atoms , @xmath2 and @xmath3 : low-@xmath0 ramsey cavities , @xmath4 : atomic field - ionization detector , @xmath5 and @xmath6 : microwave sources coupled to @xmath7 and @xmath8 s , respectively . in a quantum filtering process , a real - time control system analyzes the results of qnd measurements of the cavity field and computes the amplitude of a control injection pulse.,scaledwidth=48.0% ]    our quantum feedback algorithm is designed for the ens microwave cavity qed set - up  @xcite . its components are depicted in fig . [ fig : setup ] . the microwave field to be controlled is confined in an ultra - high @xmath0 superconducting cavity @xmath7 ( damping time @xmath9 rydberg atoms , flying one by one at 250 m / s across the cavity mode and dispersively interacting with its field , are used as qnd probes of light  @xcite . they are prepared in a superposition of two circular rydberg states @xmath10 and @xmath11 ( principal quantum numbers 51 and 50 , respectively ) in the low-@xmath0 cavity @xmath2 , driven by the source @xmath6 , in which they experience a resonant @xmath12-pulse . in the bloch sphere representation , the spin corresponding to the two - level atomic system then points along the @xmath13 direction in the equatorial plane ( states @xmath10 and @xmath11 correspond to north and south poles of the sphere , respectively ; we use a frame rotating at the atomic frequency ) .    the light shifts experienced by the non - resonant @xmath14 atomic transition in the cavity field result in a phase shift for the atomic state superposition . at the cavity exit , the atomic spin points along a direction in the equatorial plane at an angle @xmath15 with the @xmath13 axis , correlated to the photon number @xmath16 . the dephasing angle @xmath15 can be controlled by adjusting the atom - cavity detuning @xmath17 and the interaction time . in the large atom - cavity detuning regime , it is a linear function of @xmath16 : @xmath18 , where @xmath19  khz is the vacuum rabi frequency . for the intermediate detuning values used in the experiments , it is a more complex , albeit perfectly known growing function of @xmath16 . a second @xmath12 ramsey pulse in @xmath3 ( phase @xmath20 with respect to that of the pulse in @xmath2 ) followed by the atomic detection ( in the @xmath21 basis ) by the detector @xmath4 amounts to a detection of the atomic spin along an axis at an angle @xmath20 with @xmath13 . it provides information on the cavity field intensity . atoms are sent in the set - up at typical 250 @xmath22s time interval , much shorter than the cavity damping time . note that such a macroscopic time interval is well adapted to elaborate feedback strategies since we have ample time to compute the state estimator and the feedback law between two atomic detections . when no feedback action is performed , the information provided by a few tens of atoms results in a measurement of the dephasing angle @xmath15 and , hence , in a projective qnd measurement of the photon number @xmath16 @xcite . we are interested here instead in the ambiguous information provided by a single atomic detection . detection of the atomic state in @xmath4 projects the field , described initially by the density matrix @xmath23 , onto a new state @xmath24 . depending on the detected atomic state , @xmath25 or @xmath11 , the back - action of the quantum measurement on the field is described by @xmath26 where the operators @xmath27 are given by    [ eq : mgme ] @xmath28 @xmath29    here , @xmath30 is the photon number operator with @xmath31 and @xmath32 the photon annihilation and creation operators , respectively . the measurement operators @xmath27 are diagonal in a photon - number basis ( as is @xmath33 ) . they thus preserve fock states , illustrating the qnd nature of the measurement . the projected state is normalized by the probability @xmath34 of detecting the atom in state @xmath35 .    the cavity field can also be manipulated by injecting into the mode a coherent field pulse generated by the resonant microwave source @xmath5 . its action is described by the displacement operator @xmath36 , where @xmath37 is the complex amplitude of the injected field . the cavity field after displacement thus reads @xmath38    a proper analysis of the feedback scheme requires to take into account all known imperfections of the experimental set - up . besides projective measurement and coherent evolution , the cavity state also evolves due to decoherence . the field is coupled to a reservoir at non - zero temperature ( @xmath39  k ) , and its dynamics is described by the master equation @xcite @xmath40 where @xmath41 is the cavity decay rate and @xmath42 the equilibrium thermal photon number .    moreover , the circular rydberg state preparation is a non - deterministic poisson process . we perform a pulsed excitation from a continuous thermal beam of rubidium atoms , preparing atomic samples at a @xmath43s time interval . the mean number of circular atoms per sample is @xmath44 ( leading to an average time between atoms of 250 @xmath22s ) . the atomic time of flight between @xmath7 and @xmath4 is 350  @xmath22s , corresponding to a delay of @xmath45 atomic samples . the probability for detecting an atom present in a sample is @xmath46 due to finite detection efficiency in @xmath4 . finally , non - ideal atomic state resolution of the field - ionization detector and finite contrast of the ramsey interferometer introduce a probability of erroneous state assignation of @xmath47 . in order to describe the essential elements of the feedback scheme in a simple context , we first consider in this section an idealized experiment with no cavity decay ( @xmath48 ) , a deterministic atomic preparation ( @xmath49 ) , no detection delay ( @xmath50 ) , and a perfect atomic detection ( @xmath51 and @xmath52 ) .      the quantum filtering procedure @xcite provides us with an estimate of the field state by using all available information . right after detection of atom number @xmath53 and before injection of the corresponding control , it includes the initial state of the field @xmath54 , information obtained from all @xmath53 atoms detected so far ( labelled @xmath55 ) and all coherent field injections performed in the former feedback loops ( 1 to @xmath56 ) . each of the elementary processes ( atomic detections and displacements ) can be represented by super - operators acting on the field s density matrix . let us note @xmath57 that associated to the detection of atom @xmath58 and @xmath59 that corresponding to the displacement performed in the @xmath60 iteration of the feedback loop . therefore , after detecting atom @xmath53 , the state of the field  is @xmath61 here and in the following , all super - operator products are ordered by decreasing indices . note that the expression can be computed iteratively from the recurrence @xmath62 setting @xmath63 . by applying this filter recursively at each cycle , we get a reliable real - time estimate of the field state .    based on eqs .   and , the super - operators are defined as @xmath64 here , @xmath65 or @xmath66 depending on the outcome of the @xmath60 atomic measurement . the control amplitude @xmath67 is adjusted after each atom detection according to the feedback law described in the next subsection . we choose to prepare initially the cavity in a coherent state with a mean photon number equal to the target value @xmath68 . the phase of this field is used as a reference . as a consequence , its amplitude @xmath69 is real . the phase @xmath20 and the operator @xmath70 , defining the measurement super - operators @xmath71 through , should be adjusted to optimize the quantum filter  . on the one hand , the mean dephasing angle per photon @xmath72 ( estimated in the useful range of @xmath16 values ) should be large in order to resolve adjacent photon numbers . on the other hand , the measurement should lift the ambiguity between a large range of possible photon numbers from @xmath73 to some @xmath74 , determined by the spread of the initial coherent field . due to the oscillating nature of the measurement operators @xmath27 , this sets a maximum value of @xmath72 of the order of @xmath75 . finally , we need to optimally distinguish @xmath68 from @xmath76 . this is achieved by setting the ramsey phase @xmath20 so that @xmath77 . this setting corresponds to a maximal variation of the atomic state detection probabilities around @xmath68 . in other words , we set the atomic ramsey interferometer at mid - fringe when the cavity contains @xmath68 photons .      the distance between @xmath78 and the target fock state described by a density operator @xmath79 can be conveniently defined  as @xmath80 where @xmath81 is the fidelity of a state @xmath82 with respect to the target state @xmath83 . the function @xmath84 represents a natural choice for the lyapunov function @xcite used to study the stability properties of the quantum controller . this study is presented in more detail in ref . @xcite . let us denote the cavity state after the control displacement in the @xmath85 iteration of the feedback loop as @xmath86 . the optimal convergence of the feedback procedure is obtained by finding at each cycle the injection amplitude @xmath87 which maximizes the fidelity @xmath88 . this search can be performed in a straightforward way by numerical optimization . however , this approach is time - consuming and could not be realistically implemented with existing real - time data analysis systems in the few tens of microseconds time interval between two atomic detections . we overcome this difficulty by developing a faster , analytical calculation of a locally optimal displacement @xmath87 . since the measurement operators , the initial density matrix and the projector on the target fock state are all real operators , we can reasonably choose @xmath87 to be real . the field density matrix then remains always real .    in the limit of small @xmath37 , the baker - campbell - hausdorff formula @xcite yields the following approximation for a displaced state : @xmath89\notag\\ & +          \frac{\alpha^2}{2}[[\rho , a^\dag - a],a^\dag - a ] . \end{aligned}\ ] ] using  , we therefore get for small @xmath87 @xmath90{\rho_{\text{\tiny tag}}}\right ) } \\ + \frac{\alpha_k^2}{2}{\text{tr}\left([[{\rho_{k}},a^\dag\!-\!a],a^\dag\!-\!a]{\rho_{\text{\tiny tag}}}\right)}\ , . \end{gathered}\ ] ] choosing the feedback amplitude as @xmath91{\rho_{\text{\tiny tag}}}\right)},\ ] ] with a gain @xmath92 small enough ensures that @xmath93 as soon as the trace in is strictly positive . furthermore , since @xmath94=[{\rho_{\text{\tiny tag } } } , m_e]=0 $ ] and @xmath95 , the conditional expectation of @xmath96 knowing @xmath97 is given by @xmath98 therefore @xmath99 and , as a result , the expectation value of @xmath100 is a non - decreasing function of @xmath53 : @xmath101    we have only shown so far that @xmath100 is increasing on the average . it does not mean that , for _ all _ individual realizations of the loop sequence , @xmath100 converges to 1 , its maximum value reached when @xmath82 is equal to @xmath83 . indeed , the feedback law does not prevent the convergence towards other fock states , since @xmath102 whenever @xmath78 is the projector on any photon number state . a careful analysis shows that , in each realization , @xmath100 converges either to 1 or , with some small probability , to zero . to overcome this spurious attraction towards wrong photon number states , we propose , as in @xcite , to modify the feedback law   by applying a constant injection  kick " as soon as the cavity state significantly deviates from the target : @xmath103{\rho_{k}}\right ) }      & \mbox { if } f(\rho_{k } ) \geq \varepsilon,\\          \\ c_2\ , \text{sign}\left ( { n_{\text{\tiny tag}}}\!- n_{k}\right ) & \mbox { if } f(\rho_{k } ) <     \varepsilon . \\          \end{array }          \right.\ ] ] here , @xmath104 is the mean photon number in the current state @xmath78 , @xmath105 is a constant kick amplitude . the kick condition is set by @xmath106 with @xmath107 . in this way , by breaking attraction towards other fock states , we make sure that the cavity state converges towards the target one . the controller gain @xmath108 in must be tuned to maximize fidelity @xmath88 at each sampling time @xmath53 . up to third - order terms in @xmath109 , eq . yields @xmath110{\rho_{k}}\right)}\!\big)^2 \left ( c_1\ ! -\ ! \frac{c_1 ^ 2}{2 }              { \text{tr}\left([{\rho_{\text{\tiny tag}}},a^\dag\!-\!a]^2\right)}\!\right)\!.      \end{gathered}\ ] ] by maximizing this expression , we find that the gain @xmath111 ^ 2\right)}^{-1}\ ! = ( 4{n_{\text{\tiny tag}}}\!+2)^{-1}\ ] ] provides the fastest convergence speed in the vicinity of @xmath83 . a detailed mathematical analysis and a convergence proof of this feedback - scheme will be given in @xcite .      . ( a ) atomic states detected at each feedback cycle . ( b ) control injections @xmath87 . ( c ) evolution of photon number probabilities @xmath112 , @xmath113 and @xmath114 shown as dashed - dot blue , solid green and dashed red curve , respectively . ( d ) average over @xmath115 closed - loop quantum trajectories.,scaledwidth=47.0% ]    [ fig : meanideal ]    to assess the performance of the proposed feedback scheme , we have performed quantum monte - carlo numerical simulations . figure  [ fig : trajideal](a)-(c ) presents a single quantum closed - loop trajectory in the ideal measurement setting , with a target fock state @xmath116 . the hilbert space is limited to @xmath117 photons . the atom - cavity detuning is set to @xmath118  khz . the mean dephasing angle per photon around @xmath68 is then about @xmath119 . intuitively , the ramsey phase @xmath20 should be set at the mid - fringe setting for @xmath68 photons , i.e. , @xmath120 . we have numerically observed that the convergence of the process is faster when the ramsey phase alternates between four values in successive feedback cycles : @xmath121 , @xmath122 , @xmath121 and @xmath123 . the phase excursion @xmath124 is set at @xmath125  rad . the kick amplitude @xmath126 is 0.1 with a kick zone defined by @xmath127 . the initial cavity field @xmath128 is the coherent state with 3 photons on the average , @xmath129 . the upper trace in fig .  [ fig : trajideal ] shows a sequence of detected atomic state . the corresponding control inputs @xmath87 are plotted in the next panel . large at the beginning , @xmath87 rapidly converges towards zero . figure  [ fig : trajideal](c ) shows the probability for the number of photons in the cavity field to be smaller , equal or larger than @xmath68 ( dashed - dot blue , solid green and dashed red curve , respectively ) . for this particular quantum trajectory , we observe that @xmath82 converges to @xmath83 in less than @xmath130 feedback cycles . figure  [ fig : meanideal](d ) presents an ensemble average over @xmath115 realizations of the same numerical feedback experiment . the fidelity @xmath100 reaches values above @xmath131 after @xmath132 cycles . after 20 cycles only , about @xmath133 of the trajectories have converged to the target state , illustrating the efficiency of the feedback method . so far , we have neglected experimental imperfections . in this section , we take into account all known imperfections of the present experimental set - up : finite cavity lifetime , poisson distribution of the atom number in atomic samples , non - ideal efficiency and state - selectivity of the detector and , finally , the finite delay between atom - cavity interaction and atomic detection . below , we introduce modifications of the quantum filter  , which make the process efficient in spite of these imperfections . the quantum filtering process must provide us with a density matrix describing all our knowledge of the cavity field at a given time . we chose here this time to be right before the injection of the actuator pulse in the @xmath85 iteration of the loop , after detection of the @xmath85 atom . at this time , our knowledge includes all atomic measurements so far ( from 1 to @xmath53 ) , all displacements performed and the known relaxation of the cavity field in between these events . it must also include the influence of the @xmath134 atomic samples which have interacted with the cavity , but are still flying towards the detector when we record atomic sample @xmath53 . all this information can , as above , be represented by super - operators acting on the field s density matrix . let us note @xmath135 the super - operator describing the detection of an atomic sample . each detection has three possible outcomes , labelled @xmath136 , @xmath137 or @xmath138 ( atom detected in @xmath66 , @xmath137 or no detected atom at all ) . the super - operator @xmath139 describes the action of one of the @xmath134 atomic samples flying from @xmath7 to @xmath4 . the displacement super - operator @xmath140 is given by and the field relaxation during the time interval between atomic samples is described by @xmath141 . the exact expressions of all these super - operators will be given later . therefore , the state of the cavity field before injection  is @xmath142 where @xmath143 includes the information gathered by the first @xmath53 sample detections and @xmath144 includes the influence of the @xmath134 in - flight samples . note that due to the atomic propagation delay between @xmath7 and @xmath4 , the displacement amplitude injected in @xmath7 right after atomic sample @xmath58 has interacted with it is @xmath145 , computed in the @xmath146-th iteration of the loop . hence , the first @xmath134 atomic samples cross the cavity before any displacement . this is formally taken into account by setting @xmath147 to unity for non - positive indices . if @xmath50 , i.e.  for no delay in the detection process , the empty product in equals by convention the identity operator . consequently , in the case of the ideal experimental parameters , the expression naturally reduces to . the density operator @xmath148 is used in the feedback law to calculate the control injection @xmath149 . the product @xmath150 can be easily computed recursively , as in the ideal case . however , we must recalculate @xmath151 at each feedback cycle . we now give the explicit expression of all super - operators considered so far . the measurement process @xmath135 takes into account detector s imperfections and non - ideal atomic source . if the detector clicks and indicates state @xmath152 or @xmath153 , the new quantum state is a statistical weighted mixture of two projected states . the strongest weight corresponds to the projection according to the recorded outcome . with a smaller probability , the atomic detection has been erroneous . the atom was , at detection time , in the state @xmath154 , opposite to @xmath35 , and the field has been projected accordingly . thus , due to the imperfections of state discrimination alone , we get @xmath155 the weights of the two projected states are given by the conditional probability of a wrong detection knowing that the detector has clicked in @xmath156 : @xmath157 with @xmath158 . if the detector does not click ( outcome @xmath159 ) , either there was no atom in the sample or an atom was present but has not been detected . therefore , the estimated state is a mixture of the unperturbed cavity field and of the two projected states corresponding to the two possible states of the undetected atom ( we assume here that these two states are equiprobable ) . we thus get @xmath160 the weights in this mixture are set by the conditional probability to have an undetected atom in a sample : @xmath161 . the interaction with a not - yet - detected sample is described by the super - operator @xmath139 . it is given by with the conditional probability @xmath162 replaced by the atomic sample occupation probability @xmath163 , since the measurement has not been performed yet , resulting in @xmath164    the super - operator @xmath141 describes the evolution due to the cavity field relaxation during the time interval @xmath165 between two atomic samples . using in the approximation of small time interval , @xmath166 , we get @xmath167      . ( a ) sequence of detected atomic states . ( b ) control injections @xmath87 . ( c ) evolution of the photon number probabilities @xmath112 , @xmath113 and @xmath114 ( dashed - dot blue , solid green and dashed red curve , respectively ) of the field state @xmath168 deduced from the monte - carlo simulation . the thin black line gives the estimated state fidelity @xmath169 . ( d ) photon number probabilities for @xmath170 to @xmath171 in @xmath168 . for better clarity , results of only one out of each 10 feedback cycles are shown . ( e ) average over @xmath115 closed - loop quantum trajectories . if the number of quantum trajectories is large enough ( as is the case here ) , the averaged photon number probabilities calculated from the real field states @xmath168 coincide with those averaged over the estimated states . , scaledwidth=45.0% ]    [ fig : meanreal ]    we have performed extensive simulations , similar to those of the idealized measurement in fig . [ fig : trajideal ] , with the feedback algorithm adapted to the realistic experimental parameters , given in sec . we use monte - carlo simulations to follow the state of the cavity along individual realizations of the experiment . the field relaxation is taken into account by simulating random quantum jumps , whose average effect is described by eq . @xcite . for each atomic sample , we choose randomly the number of atoms ( @xmath73 or @xmath172 ) and the real atomic state ( @xmath10 or @xmath11 ) which would be detected by the ideal detector . these results determine uniquely the subsequent field state @xmath173 . we then take into account the detection imperfections to simulate the outcomes of the realistic measurement , which are used in the feedback loop for computing @xmath148 . note that @xmath168 is not accessible in the experiment . here , we use it to simulate the evolution of the field and , by comparing it with the state estimated in the quantum filtering process , to characterize the feedback performance . in the ideal case , considered in sec . [ sec : ideal ] , the estimated state @xmath148 always coincide with the real state of the field @xmath173 . figure  [ fig : trajreal](a)-(d ) shows a typical quantum trajectory , illustrating the robustness of our feedback scheme at the present level of imperfections . the first difference with the ideal case is the presence of many feedback cycles with no detected atoms . nevertheless , after about @xmath174 feedback cycles ( each of @xmath175s duration ) , corresponding to approximatively @xmath176 atoms and @xmath177 detector clicks , the cavity state successfully converges to the target fock state  @xmath178 . the other striking feature is the presence of sudden photon jumps due to the limited cavity lifetime . in this particular trajectory , a photon loss occurred at feedback cycle number @xmath179 ( at about 37  ms ) . one sees that about 60 cycles ( about 5  ms ) are required to detect the jump before applying a large control displacement and about 60 more cycles to restore the target photon number in the cavity field . ( blue solid curve ) and of the state estimated by the quantum filtering process @xmath148 ( red dashed curve ) are defined with respect to the target fock state and are obtained by averaging over @xmath115 quantum trajectories . in each trajectory , the time origin is shifted to the time of a quantum jump.,scaledwidth=48.0% ]    figure  [ fig : meanreal](e ) shows the same type of ensemble average as presented in fig . [ fig : meanideal](d ) . we observe a convergence towards a fidelity of @xmath180 after several hundred cycles resulting in about 100 detector clicks . this convergence , given as the number of actually detected atoms , is similar to that observed in the ideal case . the main difference between the two traces is the reduced asymptotic value of the fidelity of the prepared state : @xmath181 instead of @xmath182 in the ideal case . this reduction is mainly due to the cavity decay , which results in uncontrolled jumps of the photon number . the observed fidelity indicates that for a typical closed - loop trajectory the field stays in the state @xmath183 about @xmath184 of time . since the lifetime of @xmath185 is given by @xmath186  ms @xcite , we estimate that it takes on the average 26  ms to restore the initial fock state after a sudden jump .    in order to observe in more detail the recovery of the system after a sudden photon jump , we recorded @xmath115 quantum trajectories showing sudden photon losses at different times . before averaging the trajectories , we shift their individual time origins to the time of a jump . figure  [ fig : fidelityafterjump ] shows the evolution of the fidelity of the actual ( solid curve ) and estimated ( dashed curve ) cavity states with respect to @xmath83 . as calculated before , it takes about @xmath187  ms to restore the initial photon number state after a sudden photon loss . the maximum fidelity of the quantum feedback , seen in figs . [ fig : meanreal]-[fig : fidelityafterjump ] and corresponding to the ensemble of many realizations , shows how well a chosen photon number state can be preserved on average from its unavoidable decay . however , since in each feedback cycle we acquire actual information on the field , the fidelity of the deterministic generation of fock states in a single realization can be much higher . to reliably produce the target state @xmath83 we keep the feedback process running until we detect its successful convergence to @xmath83 with an estimated fidelity better than @xmath188 . figure [ fig : success](a ) shows a histogram of convergence times for @xmath189 obtained from @xmath115 quantum trajectories . the cumulative distribution , presented by a solid line , shows that after 20  ms the probability of a successful feedback outcome reaches @xmath190 and that it exceeds @xmath191 after 85  ms of the feedback operation . the inset in fig . [ fig : success](a ) shows the average density matrix of converged states with the population in @xmath185 of about @xmath192 .    . the solid line gives the probability of the feedback convergence after a given time . the inset shows the average density matrix of the converged states . ( b ) convergence for several different fock states . these curves result from an average over @xmath115 quantum trajectories.,scaledwidth=48.0% ]    our feedback scheme can be tuned to generate any small photon number states . figure [ fig : success](b ) illustrates the convergence for several fock states from latexmath:[$|{n_{\text{\tiny tag}}}\rangle =    @xmath68 , the initial photon number distribution gets broader as the photon number lifetime gets shorter . therefore , the feedback performance slowly degrades , resulting in the increased convergence time for the same required fidelity  @xmath188 . implementation of this quantum feedback scheme into a real experiment requires fast real - time analysis of the measured data , faster than duration @xmath165 of one feedback cycle . we have tested a real - time processing system @xmath194 ( jger messtechnik ) allowing for fast data acquisition and analysis in combination with complex experiment control . after careful optimization of all feedback computations , the execution of one typical feedback cycle of the proposed scheme requires about 4000 floating - point operations . our preliminary tests have revealed that the adwin system needs for this task about @xmath195s . we are now working on the integration of this control system into our cavity qed set - up . we have presented a quantum feedback protocol designed to deterministically generate small photon number states of a trapped microwave field . the reliable estimation of the cavity state at each feedback cycle allows us to follow in real time the quantum jumps of the photon number and then to efficiently compensate them , thus protecting fock states against decoherence . the two main components of the feedback algorithm are the quantum filter , which estimates the actual state of the field based on the outcomes of qnd measurements , and the feedback law computing the amplitude of a correction microwave pulse whose injection into the cavity mode maximizes the field fidelity with respect to a desired fock state . quantum monte - carlo simulations of the qnd measurements and the quantum feedback response demonstrate the high reliability of our closed - loop scheme even in the presence of realistic experimental imperfections . this scheme is also robust with respect to an imprecise knowledge of the initial state , since the repeated measurements of the state corrects initial wrong or incomplete information . convergence and stability proofs of these feedback schemes will be given elsewhere . for a first set of mathematical results , see ref . @xcite . the presented scheme can be extended to include the ramsey phase @xmath20 and the spin dephasing @xmath70 as additional control parameters of the feedback loop . this should help to optimize the qnd measurement according to our estimate of the actual field state and thus to acquire more information on its dynamics . these studies are under the way . + * acknowledgements * work supported partially by agence nationale de la recherche ( projet anr-05-blan-0200 - 01 and projet blanc cquid 06 - 3 - 13957 ) , japan science and technology agency ( jst ) and eu ( ip project scala ) .
Abstract:

Answer: we propose a quantum feedback scheme for the preparation and protection of photon number states of light trapped in a high-@xmath0 microwave cavity . a quantum non - demolition measurement of the cavity field provides information on the photon number distribution . the feedback loop is closed by injecting into the cavity a coherent pulse adjusted to increase the probability of the target photon number . the efficiency and reliability of the closed - loop state stabilization is assessed by quantum monte - carlo simulations . we show that , in realistic experimental conditions , fock states are efficiently produced and protected against decoherence .