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: over the last two decades , cosmology has entered the stage of ` precision ' science , which involve significant improvements in observational techniques , implementation of more powerful statistical and mathematical tools , and , of course , greatly advanced computational facilities . while these advancements have yielded much new insight into the subject , many important questions still remain unanswered . the nature of dark energy has been the subject of much debate over the past decade @xcite . supernovae data , which gave the first indication of the accelerated expansion of the universe , are expected to elucidate this interesting question further , as the quality of the data steadily improves @xcite . much attention in recent years has focused on determining the properties of dark energy in a model independent manner . this can be done using either parametric @xcite or non - parametric methods @xcite . a comprehensive recent review has been given by @xcite . an earlier work by @xcite suggested a non - parametric method based on smoothing the supernova data over redshift in order to reconstruct cosmological quantities , including the expansion rate , @xmath0 , and the equation of state of dark energy , @xmath5 , in a _ model - independent _ manner . in this approach , the data are dealt with directly , and one does not rely on a parametric functional form for fitting any of the quantities @xmath7 , @xmath0 or @xmath5 . the result obtained by using this approach is , therefore , expected to be model - independent . this method was shown to be successful in discriminating between different models of dark energy if the quality of data is commensurate with that expected from the future supernova acceleration probe ( snap ) . in this paper we improve the smoothing method and apply it to two recent sets of supernovae data : gold @xcite and snls @xcite . we then compare the derived expansion history of universe with the results of baryon acoustic peak observations @xcite . specifically , we use the improved smoothing method to reconstruct the hubble parameter , @xmath0 , and then derive the distance factor , @xmath1 , up to a redshift of @xmath2 independently of the assumption of any cosmological model . this derived value , based on supernovae data , is then compared with the distance factor @xmath1 ( which is also claimed to be relatively independent of dark energy model ) being determined by the detection of the baryon acoustic oscillation peak . one of the main results of this paper is that there is a good agreement between supernovae data ( both gold and snls ) and baryon acoustic peak observations for the values @xmath3 . the derived value of @xmath4 is then used to reconstruct the equation of state of dark energy for both supernovae datasets . we should emphasize here that all the results in this paper are only based on observational data and no theoretical model has been assumed . this is an advantage of this method over the functional fitting methods in which the reconstructed results are biased by an assumed functional form or a theoretical model . the paper is organised as follows . in section ii we briefly explain the smoothing method and we estimate the accuracy of the method based on the quality and the quantity of current datasets . in section iii we apply the smoothing method on the gold supernovae dataset and by using the results of detection of baryon acoustic oscillation peak , we reconstruct @xmath5 . in section iv , we carry out a similar treatment on the snls dataset . finally , in section v we discuss our results along with some concluding remarks . the method of smoothing belongs to the category of non - parametric methods of reconstruction which is complementary to the approach of fitting a parametric ansatz to the dark energy density or the equation of state . most papers using the non - parametric approach exploit a kind of top - hat smoothing in redshift space . instead , we follow a procedure which is well known and frequently used in the analysis of large - scale structure @xcite ; namely , we attempt to smooth noisy data directly using a gaussian smoothing function . in this method we apply gaussian smoothing to supernova data ( which are of the form @xmath8 ) in order to extract information about important cosmological parameters such as @xmath9 and @xmath5 . the smoothing algorithm calculates the luminosity distance at any arbitrary redshift @xmath10 to be    [ eq : bg ] d_l(z,)^s=  d_l(z)^g+n(z ) _ i & & +  exp , & & + n(z)^-1=_i exp  . & &    here @xmath11 is the smoothed luminosity distance at any redshift @xmath10 which depends on luminosity distances of each sne event with the redshift @xmath12 , and @xmath13 is a normalisation parameter . the quantity @xmath14 represents a guess background model which we subtract from the data before smoothing . this approach allows us to smooth noise only , and not the luminosity distance . after noise smoothing , we add back the guess model to recover the luminosity distance . this procedure is helpful in reducing noise in the results . since we do not know which background model to subtract , we may take as a reasonable guess that the data should be close to @xmath6cdm and use @xmath15 as a first approximation and then use a boot - strapping method to find successively better guess models . having obtained the smoothed luminosity distance , we differentiate it once to obtain the hubble parameter , [ eq : hubb ] h(z)= ^-1 , and once again to obtain the equation of state of dark energy @xmath5 , [ eq : w ] w(z)=.    in any kind of smoothing scheme for the luminosity distance , some bias is introduced both in @xmath16 and in derived quantities like @xmath9 and @xmath5 ( see appendix a1 in @xcite to find detailed calculations of the bias ) . it is important to choose a value of @xmath17 which gives a small value of the bias and also reasonably small errors on derived cosmological parameters . to estimate the value of @xmath17 in ( [ eq : bg ] ) , we consider the following relation between the reconstructed results , quality and quantity of the data and the smoothing parameters . one can show that the relative error bars on @xmath9 scale as @xcite [ eq : tegmark ] , where n is the total number of supernovae ( for approximately uniform distribution of supernovae over the redshift range ) and @xmath18 is the noise of the data . from the above equation we see that a larger number of supernovae or larger width of smoothing , @xmath17 , will decrease the error bars on the reconstructed @xmath19 , but as it has been reported earlier @xcite , the bias of the method is approximately related to @xmath20 . this implies that by increasing @xmath17 we will also increase the bias of the results . if we attempt to estimate @xmath17 such that @xmath21 , then for @xmath22 data points ( which is the number of data points in the gold sample ) , we get @xmath23 for a single iteration of our method . however , with each iteration , the errors on the parameters will increase . therefore , using this value of @xmath17 when we use an iterative process to find the guess model will result in such large errors on the cosmological parameters as to render the reconstruction exercise meaningless . it has been shown in @xcite that at the m - th iteration , the error on @xmath24 will be approximately @xmath25 , and the error on @xmath24 scales as @xmath26 . we would like the errors after @xmath27 iterations to be commensurate with the optimum errors obtained for a single iteration , @xmath28 , so we require @xmath29 . therefore , if we wish to stop the boot - strapping after @xmath30 iterations , then @xmath31 . however , after this rough estimation of the values of @xmath17 and @xmath27 , we can still play around these values to find the best combination , by minimising the likelihood of the reconstructed results to the data . in the following , we use @xmath32 and we calculate the @xmath33 of the reconstructed distance moduli to the data after each iteration , and we stop the boot - strapping process after reaching the minimum value of @xmath33 . this effect , that @xmath33 of the reconstructed results goes to a minimum value and increases again with iteration is a reflection of the problem of some iterative reconstruction algorithms which are not error - sensitive . in these cases the noise will be added to the reconstructed results after certain number of iterations and the iterative process should be stopped after reaching the minimum value of @xmath33 to get the best result . similar effect has been reported and studied in the richardson - lucy deconvolution algorithm to reconstruct the form of primordial power spectrum from cmb data @xcite .    in appendix a we show that the results are not sensitive to the chosen value of @xmath17 and also to the assumed initial guess model . the recently released gold sample @xcite , consist of 182 supernovae type ia which have been gathered from five different subsets of data , observed during the last 16 years . the range of redshift for these supernovae are between @xmath34 and @xmath35 . in this section we use this dataset to reconstruct @xmath0 , estimate the value of @xmath4 , and then reconstruct @xmath5 . we choose a flat @xmath6cdm model with @xmath36 as the initial guess model in our calculation and we fix the value of @xmath17 ( width of smoothing ) to be @xmath37 . after each iteration , we compute the @xmath33 and we stop the boot - strapping process once @xmath33 reaches its minimum value .    the @xmath33 at any iteration is calculated from the formula , [ eq : chis ] ^2_rec , j(h_0 ) = _ i and is followed by marginalising over @xmath38 . we have marginalised over @xmath38 by integrating over the probability density @xmath39 for all values of @xmath38 .    in eq.[eq : chis ] , @xmath40 is the reconstructed result at the @xmath41th iteration for the distance moduli at redshift @xmath12 , assuming the value of @xmath38 , and @xmath42 is the gold sample data given by @xcite . in figure [ chi_gold ] we show the @xmath33 of the reconstructed results at different iterations , after marginalising over @xmath38 . as we see , the @xmath33 has a minimum around @xmath43 and after this , @xmath33 is slowly increasing . so we stop the boot - strapping process at this iteration and determine @xmath0 . we can also see that for the initial guess @xmath6cdm model , the @xmath44 of the best recovered result is less than @xmath45 which means that the flat @xmath6cdm model is in agreement with the gold sample within @xmath46 .    by marginalizing over the hubble parameter , we carry out a similar treatment on the data as it has been done by @xcite to calculate the @xmath33 for different cosmological models . as the gold data are based on a hubble parameter of 65 km / sec / mpc , the reconstruction method should be able to recover this value for the hubble constant . in fact the peak of the probability density of the reconstructed result for different values of the hubble parameter should be close to @xmath47 65 km / sec / mpc . in figure [ h_margin_gold ] we show the probability density of the best reconstructed result from gold data for different values of the hubble parameter . we see that the probability density has a sharp peak around @xmath4765 km / sec / mpc . we should also note here , that the reduced @xmath33 of the reconstructed results seems to be consistently below @xmath48 ( however it is not trivial to define the degree of freedom in our smoothing method and hence the reduced @xmath33 , but we can see that the resultant @xmath33 of the reconstructed results is around @xmath49 less than the number of data points ) . we can also see in figure [ chi_gold ] that the reduced @xmath33 of the first initial guess model , which is a @xmath6cdm model , is also below @xmath48 . it shows that the error - bars of the supernovae data points are quite large and many different reconstructed results may have a reduced @xmath33 of less than @xmath48 . in this paper we only calculate the @xmath33 of the reconstructed results and we compare different results by calculating the @xmath44 to the best result with a minimum @xmath33 . @xmath50    \includegraphics[scale=0.38 , angle=-90]{h_gold_riess.ps }    \includegraphics[scale=0.38 , angle=-90]{q_gold.ps } \end{array}$ ]    @xmath50    \includegraphics[scale=0.38 , angle=-90]{bao_gold_err.ps }    \includegraphics[scale=0.38 , angle=-90]{w_gold.ps } \end{array}$ ]    in figure [ h_gold ] , left panel , we show the reconstructed @xmath0 for the gold data set . the red solid line has the highest likelihood and is our best reconstruction . all the other lines are within @xmath51 away from the best recovered result . these lines are recovered results from our smoothing method by using different numbers of iterations in the boot - strapping process . the @xmath44 for all of these lines is less than @xmath48 , and so we can consider them to lie within @xmath51 of the best result . we should note that these green dashed lines in figure [ h_gold ] are in fact a non - exhaustive sample of results which are within @xmath51 away from the best recovered result . as we see in the figure [ h_gold ] , the reconstructed @xmath0 at high redshift has a very big degeneracy . this is expected since there is only a single supernova beyond redshift @xmath52 !    in this figure we can also see 3 uncorrelated and independent measurements of @xmath0 from the gold sample ( blue dotted crosses from @xcite ) for comparison with our results . we can see that these two results are consistent with each other within their @xmath51 limits . however we should mention here that in the @xcite method used by @xcite for uncorrelated estimates of the expansion history , there is a slight bias in the reconstruction of @xmath0 and the higher derivatives of the data . it is mainly because of using the average of the measured quantity @xmath53 , which typically is not a straight line . so as it has been mentioned in @xcite , the measured average of @xmath53 ( and hence @xmath0 ) over a redshift bin will generally lie either slightly above or below the actual curve at the bin center . this can be the reason that why the centres of the crosses in figure [ h_gold ] , left panel , for uncorrelated estimates of the expansion history , are slightly above or below our reconstructed curve for the @xmath0 .    to reconstruct the hubble parameter , @xmath0 , we do not need to know the value of @xmath4 . another important cosmological quantity which we can derive from the reconstructed @xmath0 ( independent of the value of @xmath4 ) , is the deceleration parameter , @xmath54 ,    [ eq : q ] q(z)=(1+z ) -1 .    in figure [ h_gold ] , right panel , we show reconstructed @xmath54 . for the gold data our method shows that the transition between deceleration and acceleration occurred at @xmath55 ( at @xmath51 ) . the best reconstruction shows the redshift of transition to be @xmath56 . this is in agreement with results obtained using parametric methods @xcite .    to derive the equation of state of dark energy @xmath5 , one needs to know the value of @xmath4 , as we see in eq.[eq : w ] . to estimate the value of @xmath4 , without using any parameterisation and in a model - independent way , we can use the results of the detection of the baryon acoustic oscillation peak @xcite . the distance factor @xmath1 up to redshift @xmath2 , measured by observation of luminous red galaxies in detection of baryon acoustic oscillation peak ( which have been claimed to be relatively independent of the model of dark energy ) , can be derived directly for different values of @xmath4 by using the reconstructed @xmath0 , a =  _ 0^z_1  ^2/3  , where the measured value of @xmath1 is @xmath57 at @xmath58 . the 3-year wmap results , when combined with the results of baryon acoustic oscillations , yield @xmath59 for the spectral index of the primordial power spectrum @xcite . by using the best reconstructed results for @xmath0 , we get @xmath60 . in figure [ w_gold ] , left panel , we see the derived value of @xmath61 from supernovae data in comparison with its measured value from observation of lrgs for different values of @xmath4 . it is clear that these two independent observations which are completely different by nature , are very much in agreement if @xmath62 . this derived value of @xmath4 is completely independent of any dark energy model assumption ( within the framework of standard general relativity ) and is in very close agreement with the results from large scale structure measurements from 2df @xcite and sdss @xcite . this derived value of @xmath4 is also in good agreement with the results from @xcite , where a different model - independent method of reconstruction has been used .    now by marginalising over @xmath63 , which is the range of agreement between the two observations , we can reconstruct @xmath5 from our previously reconstructed @xmath0 . in figure [ w_gold ] , right panel , we show the reconstructed @xmath5 , marginalised over @xmath4 for the gold dataset . we see that the data prefer evolving dark energy to the cosmological constant . the degeneracy for the equation of state of dark energy at high redshifts is very large and it is almost impossible to say much about @xmath5 at high redshifts . in this section we use the same procedure as we used in the previous section to deal with snls supernovae data . the snls dataset contains 115 data points in the range of @xmath64 . we use this dataset , first to reconstruct the hubble parameter , @xmath0 , and the deceleration factor , @xmath54 , up to redshift @xmath48 . then by using the results of detection of baryon acoustic oscillation peak we derive the value of @xmath4 , following which we recover the form of @xmath5 . we use the distance modules of the supernovae available in tables 8 and 9 in @xcite as our dataset in this section .    in figure [ chi_snls ] we see the computed @xmath33 for the reconstructed results using smoothing method at each iteration . as we see , the @xmath33 diverges to its minimum value very fast at just the 5th iteration . in figure [ h_snls ] we show the reconstructed @xmath0 ( left panel ) and @xmath54 ( right panel ) for the snls dataset . the red solid line has the best likelihood , which is our best reconstructed result . all the other lines are within @xmath51 away from the best recovered result . we should like to emphasise here that these results ( green dashed lines ) are not representative of all the possibilities which give the likelihood within @xmath51 of the best recovered result . however they can show the overall behaviour of the quantities which we have studied . our results for snls data show that the transition from deceleration to acceleration phase of the universe occurs at redshifts higher than @xmath65 . the fact that we can not put an upper limit to the redshift of the commence of acceleration is due to the absence of supernovae data at @xmath66 in snls dataset . as we have discussed earlier in the previous section , we use the results of detection of baryon acoustic oscillation peak to determine the value of @xmath4 . then by marginalising over the recovered value of @xmath4 , we derive the dynamics of @xmath5 . in figure [ w_snls ] we see the derived value of @xmath4 and the reconstructed form of @xmath5 . we see that the @xmath6cdm model is in much better agreement with snls data than with gold data .    by comparing the recovered results from the snls and gold datasets , we can clearly see an inconsistency between these two supernovae datasets . this inconsistency is obvious by looking at the reconstructed @xmath54 and @xmath5 in the middle and high redshift ranges . gold data suggest the redshift of the commencement of acceleration at @xmath56 while snls data suggest @xmath67 . the reconstructed @xmath5 from these two datasets also shows a very different behavior in the middle and high redshift ranges . the discord between gold and snls supernovae datasets has been reported and studied earlier by @xcite , and a similar discord between gold supernovae data and other cosmological observations like cosmic microwave background observations from wmap and observations of cluster abundance , has also been reported earlier by @xcite . however snls supernovae data seem to be in good agreement with the other cosmological observations . based on all these results and analyses , we may conclude that some significant systematics in the gold data ( or in a part of the data ) might be the reasos for these inconsistencies . interestingly , the recovered values of @xmath4 from gold and snls data ( by using the results of detection of baryon acoustic oscillation peak ) , are in very close agreement . in both cases the derived value of @xmath4 is around @xmath68 . we should note here that the two data sets rely on pretty much the same nearby supernovae samples and that is why the results are similar in this range . it is something which we logically expect to get . but in fact it shows one of the advantages of this method over the functional fitting methods . by using a functional fitting method , the recovered results at any redshift would be equally dependent on the data set in the whole redshift range . but here , by using our smoothing method , we can clearly see that despite the significant differences between the reconstructed results from gold and snls datasets in the middle and high redshift regions , the reconstructed results for the expansion history at low redshifts ( which we use to estimate the value of matter density ) , are not affected by the big differences between the two datasets at the higher redshifts . @xmath50    \includegraphics[scale=0.38 , angle=-90]{h1_snls.ps }    \includegraphics[scale=0.38 , angle=-90]{q1_snls.ps } \end{array}$ ]    @xmath50    \includegraphics[scale=0.38 , angle=-90]{bao_snls_err.ps }    \includegraphics[scale=0.38 , angle=-90]{w_snls.ps } \end{array}$ ] in this paper we have shown that by improving the efficacy of the smoothing method @xcite , we can reconstruct the expansion history of the universe in a model - independent way , using current supernovae data . we have used the smoothing method to reconstruct the expansion history of the universe , @xmath0 , the deceleration parameter , @xmath54 , the value of @xmath4 and the equation of state of dark energy , @xmath5 , independently of any assumption of the theoretical model of the universe , within the framework of standard general relativity . this is an advantage of this method over the functional fitting methods where the results are usually biased by the form of the functional fitting or the assumed theoretical model . we dealt with two recent datasets , gold and snls in our analysis . in determining the value of @xmath4 , we found excellent agreement between gold and snls datasets . this determination is directly related to the supernovae data points at redshifts lower than @xmath69 . we have got @xmath70 for both gold and snls datasets , which is in good agreement with results of sdss and 2df large scale structure observations , and also with results of recent chandra x - ray observations of the relaxed galaxy clusters @xcite . this derived value of @xmath4 also agree with the recent wmap 3 years cmb data , if we assume the broken scale invariant spectrum for the form of the primordial spectrum @xcite . in the derivation of @xmath54 and the stage of transition from deceleration to acceleration in the dynamics of the universe , we found disagreement between gold and snls datasets . gold data suggest the redshift of the commence of acceleration at @xmath56 while snls data suggest @xmath67 . after marginalizing over the derived value of @xmath4 , we have reconstructed @xmath5 . the inconsistency between gold and snls supernovae datasets is also obvious by looking at the reconstructed @xmath5 from these two datasets . the derived form of @xmath5 from snls dataset , is in good concordance with @xmath6cdm model , while gold dataset prefers an evolving form of dark energy ( however @xmath6cdm is still in agreement with the gold dataset to within @xmath46 ) . this discrepancy between gold and snls datasets has been reported earlier by other groups ( @xcite ) . as the gold sample is also relatively in disagreement with the other cosmological observations like cmb and observations of cluster abundance ( @xcite ) , we may conclude that the effect of systematics in the gold dataset ( or at least in a part of the data ) is significant . the large error - bars at the high redshifts for the reconstructed results , reflect the significant lack of data points . this effect may not be seen if we use some of the parametric methods of analysis , but as we deal with the data directly here , we notice that the lack of data points at high redshifts limits our ability to say much about the behavior of the universe at the early stages at high redshifts . this is another important feature of our smoothing method in which the reconstructed results at any redshift rely mostly on the supernovae data points at the same redshift range . i would like to thank varun sahni , alexei starobinsky , tarun souradeep , eric linder , tarun deep saini , stephane fay , arnab kumar ray and ujjaini alam for useful discussions . u. alam , v. sahni , t. d. saini , a. a. starobinsky , 2003 , mon . not . . soc . * 344 * , 1057 . u. alam , v. sahni , a. a. starobinsky , 2007 , j. cosmo . . phys . * 0702 * , 002 . s. w. allen et al . , 2007 , ( astro - ph/0706.0033 ) .    p. astier et al . , 2005 , . astroph . * 447 * , 31 . v. barger , y. gao , d. marfatia , 2007 , phys . b * 648 * 127 .    c. bonvin , r. durrer , m. kunz , 2006 , phys . lett . * 96 * , 191302 . s. m. carrol , 2001 , living rev.rel . m. chevallier and d. polarski , 2001 , int . * d 10 * , 213 . t. chiba and t. nakamura , 2000 , phys . d , * 62 * , 121301(r ) . s. cole , et al . , 2005 , mon . not . . soc . * 362 * , 505 . p. coles and f. lucchin , 1995 , `` cosmology , the origin and evolution of large scale structure '' , john wiley and sons . e. j. copeland , m. sami and s. tsujikawa , 2006 , int . * d 15 * , 1753 . p. s. corasaniti and e. j. copeland , 2003 , phys . d , * 67 * , 063521 .    r. a. daly and s. g. djorgovsky , 2003 , astroph . j. * 597 * , 9 . d. j. eisenstein , et al . , 2005 , astroph . j. * 633 * , 560 . s. fay and r. tavakol , 2006 , phys . d , * 74 * , 083513 . b. gerke and g. efstathiou , 2002 , mon . not . soc.*335 * , 33 .    y. gong , 2005 , int d * 14 * , 599 . y. gong and a. wang , 2007 , phys . d , * 75 * , 043520 . z. k. guo , n. ohta , y. z. zhang , 2005 , phys . d , * 72 * , 023504 . z. k. guo , n. ohta , y. z. zhang , 2007 , mod . a * 22 * , 883 . d. huterer and m. s. turner , 1999 , phys . d , * 60 * , 081301 . d. huterer and g. starkman , 2003 , phys . lett . * 90 * , 031301 . d. huterer and a. cooray , 2005 , phys . d * 71 * , 023506 . d. huterer and h. v. peiris , 2007 , phys d * 75 * , 083503 . h. k. jassal , j.s.bagla , t.padmanabhan , 2005 , phys . d * 72 * , 103503    r. a. knop et al . , 2003 , astroph . j. * 598 * , 102(k ) . legacy archive for microwave background data , web site : _ http://lambda.gsfc.nasa.gov/ _ e. v. linder , 2003 , phys . lett . * 90 * , 091301 . e. v. linder , 2003 , phys . d , * 68 * , 083504 . e. v. linder , d. huterer , 2005 , phys . d , * 72 * , 043509 . i. maor , r. brustein , j. mcmahon and p. j. steinhardt , 2002 , phys . d , * 65 * , 123003 . v. j. martinez and e. saar , 2002 , `` statistics of galaxy distribution '' , chapman and hall . s. nesseris and l. perivolaroupolos , 2004 , phys . d , * 70 * , 043531 . s. nesseris , l. perivolaropoulos , 2007 , j. cosmo . . phys . * 0702 * , 025 . t. padmanabhan , 2003 , phys . * 380 * , 235 . p. j. peebles and b. ratra , 2003 , rev.mod.phys . 75 , 559 . s. j. perlmutter et al . , 1999 , j. * 517 * , 565 . a. g. riess et al . , 1998 , astron . j. * 116 * , 1009 . a. g. riess et al . , 2005 , j. * 607 * , 665(r ) . a. g. riess et al . , 2006 , ( astro - ph/0611572 ) . t. roy choudhury and t. padmanabhan , 2005 , astron . astrophys . * 429 * , 807 . m. sahln , a. r. liddle , d. parkinson , 2007 , phys . d , * 75 * , 023502 . v. sahni and a. a. starobinsky , 2000 , int . * d 9 * , 373 . v. sahni , 2004 , lect . notes phys . 653 , 141 . v. sahni and a. a. starobinsky , 2006 , int . * d 15 * , 2105 . t. d. saini , s. raychaudhury , v. sahni , a. a. starobinsky , 2000 , phys . * 85 * , 1162 . t. d. saini , 2003 , mon . not . . soc . * 344 * , 129 . t. d. saini , j. weller and s. l. bridle , 2004 , mon . not . * 348 * , 603 . a. shafieloo and t. souradeep , 2004 , phys . d , * 70 * , 043523 . a. shafieloo , u. alam , v. sahni , a. a. starobinsky , 2006 , mon . not . . soc . * 366 * , 1081 . a. shafieloo and t. souradeep , 2007 , ( arxiv:0709.1944 ) . d. spergel et al . , 2007 , j. supp * 170 * , 377 .    a. a. starobinsky , 1998 , jetp lett . * 68 * , 757 . m. tegmark , 2002 , phys . d , * 66 * , 103507 . m. tegmark , et al . , 2004 , astroph . j. * 606 * , 702 . j. l. tonry et al . , 2003 , astroph . j. * 594 * , 1 . y. wang and g. lovelace , 2001 , astroph . j. * 562 * , l115 . y. wang and p. mukherjee , 2004 , astroph . j. * 606 * , 654 . y. wang and m. tegmark , 2004 , phys . * 92 * , 241302 . y. wang and m. tegmark , 2005 , phys . d , * 71 * , 103513 . j. weller and a. albresht , 2002 , phys . d , * 65 * , 103512 . in this section we show that the results of the smoothing method are robust against the choice of the initial guess model and also to the chosen value of @xmath17 . we assumed three different cosmological models as our initial guess model and we applied our smoothing method on the gold dataset . the final results by using these three different initial guess models are almost identical with @xmath71 . we have got @xmath72 by using a flat @xmath6cdm model with @xmath36 as the initial guess model after 89 iteration , while we have got @xmath72 for a flat @xmath6cdm model with @xmath73 after 91 iteration , and @xmath74 for a flat quiessence model with @xmath36 and @xmath75 after 104 iteration . in figure [ gt_fig ] we can see the reconstructed @xmath0 and @xmath5 for the gold data set by assuming these three different initial guess models . as we see , the robustness of the method for the choice of the initial guess model is obvious . we have also used different values of @xmath17 ( width of smoothing in eq.[eq : bg ] ) , in our reconstruction process to check the reliability and stability of our results against the changes in the value of @xmath17 . we have used three values of @xmath17 equal to @xmath76 and @xmath77 in our smoothing method and we have applied it on the gold dataset . by using @xmath78 we have got @xmath79 after 9 iteration , while we have got @xmath72 by using @xmath80 after 89 iteration , and @xmath81 by using @xmath82 after 407 iteration . in figure [ wt_fig ] we can see the reconstructed @xmath0 and @xmath5 for the gold dataset by using these three values of @xmath17 . we can clearly see that the results are not sensitive to the given value of @xmath17 . these two examinations confirm the overall robustness of the method for different initial assumptions .
Abstract:

Answer: we have improved upon the method of smoothing supernovae data to reconstruct the expansion history of the universe , @xmath0 , using two latest datasets , gold and snls . the reconstruction process does not employ any parameterisation and is independent of any dark energy model . the reconstructed @xmath0 is used to derive the distance factor @xmath1 up to redshift @xmath2 and the results are compared with the given value of @xmath1 from detection of baryon acoustic oscillation peak ( bao ) . we find very good agreement between supernovae observations and the results from bao for @xmath3 . the estimated values of @xmath4 are completely model - independent and are only based on observational data . the derived values of @xmath4 are then used to reconstruct the equation of state of dark energy , @xmath5 . using our smoothing method we can demonstrate that while snls data are in very good agreement with @xmath6cdm , the gold sample slightly prefers evolving dark energy . we also show that proper estimation of the equation of state of dark energy at the high redshifts would be impossible at the current status of observations . cosmology : theory  cosmological parameters


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: investigations on organic molecular crystals pioneered by pope and cowokers @xcite more than half a century ago made an immense impact on various fields of relevance in physics , chemistry , and materials science . the interplay between geometric and electronic structures has opened up novel materials with unlimited possibilities . the discovery of conducting properties of doped polyacetylene @xcite in 1977 led to a nobel prize in chemistry establishing a whole new field of organic electronics . after decades of research , many families of the organic semiconducting materials have been investigated . among them , oligoacenes @xcite and oligothiophenes @xcite have been examined intensively due to the possibility of making single crystals with few defects via repeated sublimation processes . single crystals of these materials can be used to develop transistors that allow for fundamental research on intrinsic properties as well as performance of organic electronic devices . the holstein molecular crystal model was introduced in the 1950s @xcite to account for the effect of electron - phonon interactions on transport . since then , a lot of work has been carried out to deal with the complexities of the transport phenomenon in organic molecular crystals . as reviewed in ref . @xcite , charge transport in organic semiconductors is governed by electronic coupling and electron - phonon interactions , and the main transport mechanism can be described by polaron and disorder models . many theories on polaron transport are based upon a phenomenological model @xcite or utilize the polaron effective mass approach @xcite . however , validity of these models is restricted to specific parameter ranges . more general microscopic models in understanding polaron transport was developed @xcite from a density matrix approach , which is capable of describing electronic coupling , diagonal and off - diagonal electron - phonon interaction of arbitrary strength over a wide range of temperatures . the effect of off - diagonal coupling poses much more difficulties to tackle than its diagonal counterpart , and as a result , off - diagonal coupling has not been extensively treated before @xcite . more recently , microscopic models from generalized master equation approach @xcite , and dynamical mean - field theory @xcite have been developed . polaron transport theory developed by holstein in his seminal work@xcite was based upon perturbation expansions with its validity sometimes limited to cases of small transfer integrals @xcite . despite its limitation , the theory has been widely used for qualitative interpretations of experimental data , including the temperature dependence of band - narrowing effect as well as the crossover from band - like behavior to hopping transport with increasing temperature . a microscopic transport theory that accounts for simultaneous diagonal and off - diagonal electron - phonon coupling was developed by munn and silbey@xcite based on a canonical transformation of the hamiltonian with an element of optimization . off - diagonal coupling is found to increase the polaron binding energy and introduce new minima and broadening to the polaron band . in contrast , the effect of diagonal coupling is restricted to band narrowing . the yarkony - silbey variational approach@xcite provides an attractive direction to solve the transport problems in organic molecular crystals . this approach has been borrowed to treat two phonon bands in a three - dimensional lattice by parris and kenkre@xcite . however , as the yarkony - silbey approach lacks the flexibility in variational parameters , a more sophisticated approach is needed to capture the correct temperature dependence of diffusion coefficient and mobility . to obtain the optimal basis for a polaron system , finite - temperature variational method by cheng and silbey @xcite combining merrifield s transformation with bogoliubov s theorem has been recently devised bringing substantial improvements to results obtained previously@xcite . the theory is able to capture the universal band - like - to - hopping transition as temperature increases , and more importantly , a temperature - independent mobility at extremely low temperature has been found in agreement with a nonperturbative approach recently developed by ortmann _ however , merrifield s transformation itself is applicable only in the small polaron regime , and the off - diagonal exciton - phonon coupling was not accounted for in ref . @xcite , making the approach insufficient to describe the phonon assisted transport / tunneling@xcite in molecular crystals . thus , further improvements are necessary for cases with simultaneous diagonal and off - diagonal coupling . following up on the munn - silbey transformation method , zhao _ _ et al.__@xcite devised a self - consistent routine to determine the transformation coefficients , demonstrating substantial corrections of the polaron band structure thanks to a built - in momentum dependence of the transformation coefficients . this serves as the starting point towards the evaluation of transport properties in this paper . an attempt to extend the holstein model to higher dimensions with a microscopic model has been made by kenkre _ _ et al.,__@xcite utilizing the generalized master equation approach . this model was able to adequately explain the temperature dependence and anisotropy of measured mobilities in naphthalene . using master equations , wang _ et al . , _ @xcite developed a nonperturbative method to handle the electron - phonon coupling fully quantum - mechanically and quantify charge - carrier transport properties in organic molecular crystals . a more sophisticated microscopic model based on a hamiltonian of the holstein - peierls type has been presented by bobbert and colleagues @xcite . very recently , a theory based on non - perturbative evaluation of kubo formula for the carrier mobility@xcite has been proposed , showing several improvements including the elimination of low - temperature singularity that often appears in theories based on narrow - band approximations . based on a non - perturbative evaluation of the kubo formula , calculations have been carried out for durene crystals@xcite and guanine based materials@xcite , and transport channels of charge carriers are revealed by making use of ab - initio tools . coupling to the intermolecular acoustic modes @xcite was found to play a significant role in charge transport , thus the model presented by bobbert _ _ may be improved further . experimental data of naphthalene can be well reproduced by this model with microscopic parameters obtained from ab - initio calculations . however , this model neglects the intra - molecular modes as well as the coupling of excitons to the intermolecular acoustic phonons , and only the coupling to the optical modes is accounted for . moreover , the models developed by bobbert _ et al . , _ and by munn and silbey , are based on nonlocal canonical transformations with additional approximations . thus the range of validity of these models still requires further exploration , despite some qualitative agreements with experiments . hence , in the light of recent theories by hultell _ _ et al.__@xcite and troisi _ _ et al.__@xcite , it is clear that both local and nonlocal exciton - phonon interactions should be taken into account . a comprehensive theory of the charge transport in organic crystals requires treatments of the hamiltonian having the lattice dealt with quantum mechanically in order to make itself valid for any temperature of interest , and the munn - silbey approach@xcite combined with the self - consistent routine devised by zhao _ _ et al.__@xcite provides a good opportunity to explore the simultaneous effect of diagonal and off - diagonal exciton - phonon coupling on the transport properties of polaron . the rest of the paper is organized as follows . the theoretical models are firstly introduced in section ii , where the essential improvements over the original munn - silbey model have been proposed . in section iii , the numerical results obtained with our approach have been discussed in detail and comparisons with previous results are made to illustrate the success of our approach . finally , the conclusions are drawn in section iv . the generalized holstein hamiltonian and its transformation have been described in detail previously by munn and silbey@xcite . here we will first discuss briefly the canonical transformation of the hamiltonian which yields a weak residual excitation - phonon coupling . evaluation of correlation functions , essential to the final calculation of diffusion coefficient and mobility as a function of temperature , follows next . our major focus in this paper will be the explicit form of the correlation functions , evaluation of critical exponential matrices and especially the transformation coefficients without resorting to the approximations in ref . @xcite . in our study , the transformation coefficients are determined through solving a set of self - consistency equations , rather than an analytical form approximated in ref . @xcite , therefore providing substantial corrections in the transport calculations . in sec . ii a and b , expressions for the scattering and the hopping rates as well as the diffusion coefficients are introduced . correlation functions are presented in sec . ii c.      the generalized holstein hamiltonian incorporating simultaneous diagonal and off - diagonal coupling @xcite has the form : @xmath0 here @xmath1 is the creation ( annihilation ) operator of an excitation ( i.e. , exciton or charge carrier ) with energy @xmath2 at site n , and @xmath3 creates ( destroys ) a phonon with frequency @xmath4 and wave vector _ q_. @xmath5 stands for the electronic transfer integral between site @xmath6 and @xmath7 . in this study , the discrepancies between molecules have been neglected , thus the on - site energies are set to @xmath2 as shown in eq.[hamiltonian0 ] . finally , the excitation - phonon coupling is described by the last term of eq.[hamiltonian0 ] , where the coupling strength @xmath8 must satisfy the relationship : @xmath9 , to make @xmath10 hermitian and translationally invariant . a unitary transformation can be applied to the hamiltonian , as shown below : @xmath11 with @xmath12    the transformation requires @xmath13 in order to keep the translational symmetry of @xmath14 , similar to @xmath8 . details of the transformation have been explicitly introduced in ref . @xcite .    moreover , the hamiltonian in momentum representation is more convenient than the site representation due to the translational symmetry . .  [ hamiltonian0 ] can be written in the momentum representation as : @xmath15 with @xmath16 , and @xmath17 where @xmath18 , and @xmath19 is the phonon frequencies . the excitation - phonon coupling coefficients can also be written as : @xmath20 with following coupling geometry : @xmath21,\ ] ] which indicates that the antisymmetric type of coupling is used throughout the paper . now the transformation takes the following form in the momentum representation : @xmath22 and @xmath23    the transformation can then also be written as : @xmath24 with @xmath25 where @xmath26 is an operator creating a net phonon with momentum @xmath27 . this operator is essential towards the evaluation of a series of important matrices , where @xmath28 with @xmath29_{kk ' } \nonumber \\ \theta_{kk'}^{\dagger}=[{\rm exp}(\textbf{s})]_{kk'}\end{aligned}\ ] ] where @xmath30 stands for matrix @xmath31 . the transformed hamiltonian can be partitioned into @xmath32 terms ( zeroth order part ) and the perturbative @xmath33 terms . thermal average of the @xmath34 operator has the property : @xmath35 which ensures the correct thermal equilibrium behavior of @xmath36 by keeping it diagonal in excitation wave vector . @xmath36 here has following expression : @xmath37 here @xmath38 represents the renormalized transfer integral which can be expressed by : @xmath39    the perturbation parts of the transformed hamiltonian can be written as : @xmath40\}\nonumber\\ & + & n^{-1/2}\sum_{qk}\omega_{q}[a_{-k}^{q}-\sum_{k'}f_{-k'}^{q}\langle\theta_{k'+q , k}^{\dagger}\theta_{k',k - q}\rangle]\nonumber\\ & \times&(b_{q}+b_{-q}^{\dagger})a_{k+q}^{\dagger}a_k,\end{aligned}\ ] ] where @xmath41    however , the last term of @xmath33 has the potential to grow with increasing temperature , motivating the choice of @xmath42 : @xmath43 to curb this possible uncontrolled growth of @xmath34 .    by applying the thermal average routine shown in ref.@xcite , the transformed hamiltonian can be written in a simpler form , which facilitates the numerical calculations . the zeroth order hamiltonian @xmath32 now can be written as : @xmath44 and @xmath45 the renormalized transfer integral in eq . [ epsilontilde ] is given by : @xmath46 where the triadic @xmath47 is a quantity introduced to yield a simpler form of the transformed hamiltonian , with following relationship : @xmath48 here @xmath49 is the bose - einstein distribution , with @xmath50 .    with the form of @xmath51 shown in eq . [ secular ] , the perturbation part @xmath33 of @xmath52 can be written as : @xmath53    in this study , the transformation coefficients @xmath51 ( a - matrix ) will be obtained variationally prior to the calculations of the transport properties . the detailed numerical variational procedures can be found in ref . @xcite . in performing the variational calculations , the transformation coefficients must satisfy following self - consistency equations together with eq . [ selfconsistent1 ] : @xmath54,\ ] ] @xmath55_{kk'}.\ ] ]    the self - consistency equations above follow the munn - silbey secular elimination scheme . in ref . @xcite , this set of equations was solved upon further approximates on @xmath42 , whereas in this paper , it is solved numerically with an accuracy unavailable previously .    in the munn - silbey approach@xcite , the transformation coefficients @xmath56 was at first approximated by the scaling parameters @xmath57 and @xmath58 as : @xmath59.\ ] ]    numerical evaluation of the eqs . [ selfconsistent1]-[selfconsistent3 ] is facilitated by rewriting the eq . [ variational1 ] with real matrices @xmath60 and @xmath61 in the following form : @xmath62.\ ] ]    in the equation above , the sine function leads to zero imaginary part along the line with @xmath63 . in the numerical calculation , the lines with @xmath64 and @xmath65 are treated as removable singularities . moreover , the values of @xmath61 along these lines are chosen to be analytically connect with neighboring values@xcite . the variational matrices @xmath60 and @xmath61 chosen can preserve the symmetry properties of @xmath51 in a way such that @xmath66 and @xmath67 .    the self - consistency approach shown in this paper is largely limited by the secular elimination scheme that is independent on the rigid - lattice tunneling matrix elements @xmath68 . however , since @xmath68 determines the rigid - lattice band structure and related quantities such as effective mass and density of states , this approach only applies for the narrow ( rigid ) band regime where the traditional small polaron theory is still available . the wide ( rigid ) band systems are beyond the scope of the approach discussed in this work . based on the simple introduction and discussion of the hamiltonian and transformation schemes , the parameters used for our numerical calculations are limited to the regime which holds the traditional small polaron scenario . the diffusion coefficient is closely related to the mean square displacement of an excitation in a given period of time . for the transformation scheme shown in this paper , the mean square displacement can be calculated from the excitation density matrix . approximations have been made to reduce the complexity of the exact formulism of the mean diffusion coefficient by neglecting small terms of the perturbation @xmath33 as shown in the section above . the diffusion coefficient can then be expressed as@xcite : @xmath69 where @xmath70 is the nearest neighbor distance , @xmath71 is the wave vector , and @xmath72 . the double bracket denotes the thermal average of the polaron states @xmath73 , where @xmath74 and @xmath75 are the scattering and hopping rates expressed as @xmath76 where @xmath77 , crucial to the calculation of diffusion coefficient , are given by @xmath78 \nonumber \\ & + & \langle{v_{k'+q , k+q}(t)v_{kk'}}\rangle{\rm exp}[-i(e_{k}-e_{k'})t]\},\end{aligned}\ ] ] here the single angular bracket denotes the average of the residual interaction @xmath79 over phonon states . using the formulism above , the evaluation of @xmath80 will be done numerically here without introducing further approximations as in ref . @xcite . with the diffusion coefficient d , the mobility can be easily obtained with the einstein relation @xmath81 .      as inferred from the expression of the diffusion coefficient , the central problem is the evaluation of the correlation functions @xmath82 since both the scattering and the hopping rates are determined by the @xmath83 quantities . it is indicated from the previous formulism that the residual interaction @xmath84 is closely related to the transformation coefficients @xmath51 . hence , the problem of evaluating @xmath80 ultimately hinges on the solution of the self - consistency equations , i.e. , eqs . [ selfconsistent1 ] , [ selfconsistent2 ] and [ selfconsistent3 ] . in this paper , the transformation coefficients that determine the structure of polaron states and energy bands are solved numerically , which provide new information about the transport properties of polarons in molecular crystals .    without the analytical approximation form of @xmath58 carried out in ref . @xcite , the important quantity @xmath85 can be derived as@xcite : @xmath86 and @xmath87 & = & \frac{1}{n}p_{k - k^{\prime } } ( t)(a_{k - q}^{k - k^{\prime } } ) ^{\ast } a_{k}^{k - k^{\prime } } , \\ ( d^{q})_{k , k^{\prime } } [ t ] & = & \frac{1}{n}p_{k - k^{\prime } } ( t)(a_{-k^{\prime } } ^{k - k^{\prime } } ) ^{\ast } a_{-k^{\prime } + q}^{k - k^{\prime } } , \end{aligned}\ ] ] where the expression of @xmath88 is given by eq . [ selfconsistent2 ] . the factors @xmath89 are the characteristics of correlation functions for linear exciton - phonon coupling @xcite , which may take the form @xmath90/n,\end{aligned}\ ] ] where @xmath91 is also the bose - einstein distribution , with @xmath92 and @xmath93 is the mean phonon frequency .    in order for the correlation to decay to zero at long times , @xmath89 should be multiplied by a decay factor @xmath94 , defined as the fourier transform of the phonon density of states , which takes a gaussian form @xcite as given by : @xmath95 where @xmath96 is the phonon bandwidth taken much less than @xmath93 . for the convenience of numerical calculations , the mean phonon frequency @xmath93 is set to 1 . for all calculations performed in this paper , @xmath96 is set to be 0.1 , indicating a rather narrow phonon band@xcite . the inclusion of phonon bandwidth is a natural result of exciton - phonon coupling . for the holstein polaron model , it has been inferred @xcite that phonon dispersion is directly linked to lattice relaxation time , and it is also crucial for the calculation of probabilities of hopping events accompanied simultaneous absorption and emission of phonon . it is also pointed out by holstein@xcite that replacing the phonon dispersion by a single frequency ( einstein model ) yields meaningless results for probabilities of hopping events . a number of factors can result in a finite bandwidth in the phonon density of states , and as pointed out by other authors @xcite , the inclusion of the phonon bandwidth ensures the long - time decay of correlation functions .    in the evaluation of @xmath82 , an important quantity is the four @xmath97 correlation function which is given by : @xmath98_{k_3+r_3-q_1-k_4-r_1,-q_2-r_2 } \nonumber \\ & \times&[\exp ( e^{k_1+r_2-q_3}(t))]_{-q_4-r_3,r_2-q_3 } \nonumber \\ & \times&[\exp ( e^{k_3-q_1-r_1}(t))]_{-q_1,k_3+r_3-q_1-k_4-r_1 } \nonumber \\ & \times&[\exp ( -e^{k_4-q_4-r_3}(t))]_{-q_4,-q_4-r_3 } \nonumber \\ & \times&\delta _ { k_1+q_1+k_4+q_4,k_2+q_2+k_3+q_3}.\end{aligned}\ ] ]    now we can derive the expression for @xmath82 , which is the core item of w. here we first express @xmath99 as follows : @xmath100\}\end{aligned}\ ] ] with @xmath101 and @xmath102 denotes the polaron band with its width determined by @xmath103 . as described in the previous section , the currently applied approach is only applicable for narrow band rigid lattice system . consequently , the value of @xmath103 should not be too large in order to guarantee the accuracy of the transport results obtained . the correlation function is now given by : @xmath104\ } \nonumber\\ & \times&\underset{k_4}{\sum } \{j_{k_4}t_{k_4,q_1;k_4,q_2}-\frac{1}{\sqrt{n}}\underset{q_2}{\sum } \omega _ { q_2}f_{-k_4}^{q_2 } \nonumber\\ & \times&[\frac{1}{\sqrt{n}}(a_{-q_1}^{-q_2}t_{k_4+q_2,q_1;k_4,q_2-q_2}\nonumber\\ & + & a_{-q_1}^{q_2\ast}t_{k_4+q_2,q_1+q_2;k_4,q_2})\nonumber\\ & -&\frac{1}{2}(t_{k_4+q_2,q_1;k_4,q_2}\psi _ { q_2}+\psi_{q_2}t_{k_4+q_2,q_1;k_4,q_2})]\}\rangle\nonumber . \\\end{aligned}\ ] ] the explicit expression of the correlation function has been given in the appendix a.    as read from the expression of the correlation function in eq . [ vvt1 ] , matrices @xmath105 play an essential role towards an accurate evaluation of the correlation function . in ref . @xcite , these matrices were obtained only approximately , while in this paper , they are obtained numerically with unprecedented accuracy . one will find that the accurate evaluation of the @xmath105 has dramatic influence on the overall transport calculations in the following section , especially for the cases with off - diagonal excitation - phonon coupling . ( a)-(c ) and @xmath61 ( d)-(f ) as functions of increasing diagonal coupling strength @xmath106 . @xmath107 for ( a ) and ( d ) , @xmath108 for ( b ) and ( e ) , and @xmath109 for ( c ) and ( f ) . the other parameters are : j=0.1 , @xmath110=0.3 and t=1.0 . ] the numerical evaluation of transformation coefficients @xmath105 involves the solution of the self - consistency equations eqs . [ selfconsistent1 ] , [ selfconsistent2 ] and [ selfconsistent3 ] . details of the evaluation procedures have been discussed in detail in ref . @xcite , here we only revisit some important issues related to the transport properties , especially the role of off - diagonal coupling .    for moderate temperature with @xmath111 , the self - consistency method applied in this paper tends to fail in strong coupling regimes , this is attributed to the failure to minimize errors to desired precision or converging to multiple solutions that are sensitive to procedure initialization . it is also possible for the convergence to fail in strong coupling regimes for solely numerical reasons . this problem increases with increasing excitation - phonon coupling , which lead to large positive and negative eigenvalues for the exponential matrices , thus to the possible failure of obtaining reliable diffusion coefficient . due to the problems of convergence and possible multiple solutions of the method applied , we restricted our studies in the weak and intermediate coupling regimes to ensure the robustness of the results obtained . as demonstrated in ref . @xcite , a - matrix parameters @xmath60 and @xmath61 in eq . [ variational2 ] obtained with the self - consistency procedure have a significant wave vector dependence , which is lacking in the original treatment of munn and silbey @xcite as indicated in eq . [ variational1 ] . in addition , the munn - silbey parametrization of the a matrix , @xmath58 and @xmath57 , can be regarded as the lower bounds of @xmath60 and @xmath61 rather than the average values . more importantly , the wave vector dependence of the variational parameters lead to significant changes in the structure of polaron energy band @xmath73 , where the bimodal variation of @xmath73 can be always found . such changes in the structure of polaron states and the energy band lead to critical modifications of the transport properties , which will be explicitly examined in the next section . in this section , we present detailed structures of a - matrix parameters as well as the corresponding polaron bands obtained with the self - consistency routines , in order to facilitate understanding of calculated transport properties , which are deferred to the next subsection . = 0.1 and ( c ) @xmath110=0.3 at t=0 , where @xmath106 and @xmath103 are both equal to 0.1 for all calculations . the inset of ( a ) is the band structure comparison between the result from ref . @xcite by munn and silbey and that from this work at t=1 with @xmath106 = @xmath110 = @xmath103 = 0.1 . ] ( a)-(c ) and @xmath61 ( d)-(f ) as functions of increasing temperature . @xmath112 for ( a ) and ( d ) , @xmath113 for ( b ) and ( e ) , and @xmath114 for ( c ) and ( f ) . other parameters are j=0.1 , @xmath106=0.1 and @xmath110=0.3 . ]    as shown in fig . [ figure1 ] , the two sets of a - matrix parameters in eq .  [ variational2 ] , @xmath60 and @xmath61 , differ in their wave vector dependence . in the vicinity of @xmath64 , @xmath60 is weakly structured and relatively flat with respect to @xmath115 , while the wave vector dependence of @xmath61 is much more pronounced than that of @xmath60 . regardless of coupling strength , the wave vector dependence of @xmath60 and @xmath61 has a characteristic `` manta ray '' shape , a fact that may be understood through an examination of the real - space structure of the a matrix @xcite . as revealed in fig . [ figure1 ] , the a - matrix parameters tend to have less wave vector dependence for larger @xmath116 . if diagonal coupling strength exceeds that of off - diagonal coupling , the parameters lose most of the @xmath117 dependence . thus , the transformation coefficients in eq .  [ variational1 ] , as approximated by munn and silbey @xcite , represent the a - matrix in eq . [ variational2 ] in the limit of strong diagonal coupling . it is also interesting to look into the change of the polaron band structure ( especially , band narrowing ) as a function of the temperature , which may provide an intuitive understanding of the difference between transport properties obtained here and those by munn and silbey using eq . [ variational1 ] . variance of the polaron bandwidth as well as the entire band with respect to temperature has been obtained with our numerically determined a - matrix , and is illustrated in fig . [ figure2 ] . it is also indicated in the inset of fig . [ figure2](a ) , the polaron band obtained with the transformation coefficients of munn and silbey @xcite lacks a bimodal structure , and has a much larger bandwidth than that from the current approach . it is indeed counterintuitive to see that the greater momentum - space variability we find in our transformation coefficients should translate into a narrower polaron band . the explanation lies in that the binding energy is proportional to the average of @xmath118 . compared with results from ref . @xcite , the average binding energy obtained in the current approach is significantly larger and the average debye - waller factor significantly smaller , which eventually leads to the narrower band .    to lend further support to the polaron band structure calculated here , it is helpful to introduce the toyozawa ansatz @xcite for zero - temperature band comparisons . as shown in figs . [ figure2 ] ( b ) and ( c ) , good agreement between the current method and the toyozawa ansatz has been found . for weak diagonal and off - diagonal coupling as is the case in fig . [ figure2 ] ( b ) , the toyozawa method yields a larger bandwidth , and with increasing @xmath119 , the bimodal structure appears and the agreement between the two approaches improves . the concurrence provides legitimacy to the current approach in finding reliable polaron bands , a step that is essential for the transport calculations . moreover , the band narrowing effect is in good qualitative agreement with those obtained with the holstein - peierls model and ab - initio calculations in ref .  @xcite . as can be readily seen from fig . [ figure2 ] ( a ) , the off - diagonal coupling strength helps boost the bandwidth for temperatures lower than @xmath93 , @xmath107 , and @xmath120 .    to better understand the bandwidth narrowing effect , the parameter matrices @xmath60 and @xmath61 have been plotted in fig . [ figure3 ] for three temperatures , @xmath121 , and @xmath114 in unit of @xmath93 . other control parameters in fig . [ figure3 ] are set to be the same as those for fig . [ figure2](c ) . as temperature increases , the wave vector variations grow for both @xmath60 and @xmath61 , while in fig . [ figure1 ] , the variations are reduced as @xmath116 increases . as explained in ref . @xcite , a larger wave vector variation in @xmath122 and @xmath61 is transformed into a weaker distortion of the energy band , which explains the temperature dependence in fig . [ figure2 ] ( a ) . in contrast to an overall enhancement of polaron band - narrowing by including off - diagonal coupling in ref . @xcite , our model demonstrate the possibility that at low temperatures off - diagonal coupling may increase the bandwidth while at high temperature the effect is reversed . such an effect can be understood from eqs . [ epsilontilde]-[selfconsistent3 ] , where the bandwidth is shown to be closely related to the renormalized transfer integral @xmath38 and the polaron binding energy , which is the third term on the right hand side of eq . [ epsilontilde ] . at lower temperatures , the binding energy dominates in eq . [ epsilontilde ] in the presence of off - diagonal coupling strength @xmath119 , resulting in a bandwidth that increases with @xmath119 . at high temperatures , however , the renormalized transfer integral @xmath38 dominates at weak off - diagonal coupling , and the bandwidth derived from the @xmath38 term shrinks with increasing @xmath119 . stronger off - diagonal coupling leads to stronger correlation which helps to form more spatial extended polaron structure@xcite , which naturally leads to higher coherent phonon - assisted exicton transfer , thus the polaron bandwidth at low temperature for stronger @xmath119 case is lager than the smaller one . however as temperature increases , the inclusion of stronger off - diagonal coupling strength helps to increase the thermal fluctuation of exciton ( or polaron ) through the renormalized j term in eq . [ epsilontilde ] , thus in this case phonon has more significant role in reducing the bandwidth . the influence of the diagonal coupling strength @xmath116 and off - diagonal coupling strength @xmath119 ( especially the latter ) is substantial in determining momentum space modulations of transformation coefficients @xmath105 , which in turn determine the structure of polaron states and the energy bands . for all calculations here , the number of site is equal to 6 and the hopping integral j is 0.1 in unit of @xmath93 unless specified otherwise . comparisons between the results obtained here and those from ref . @xcite have been made to illustrate the importance of the a - matrix momentum dependence as well as the straightforward numerical evaluation of transport properties sans additional approximations . for a 1 dimensional chain with 4 , 6 and 8 sites , where the diagonal coupling strength @xmath106 = 0.1 for ( a)-(c ) , and @xmath106 = 0.5 for ( d)-(f ) . the off - diagonal coupling strength @xmath110 is equal to @xmath123 for ( a ) and ( d ) , @xmath108 for ( b ) and ( e ) , and @xmath124 for ( c ) and ( f ) . finally , the transfer integral j=0.1 for all cases . ]    in the calculations of the transport properties , computation time increases dramatically with the number of sites . our existing computational capabilities allow only calculations for less than 10 sites without overshooting memory requirements . to ensure the robustness of the results obtained with current method , we computed the diffusion coefficient d for 1d lattices of 4 , 6 and 8 sites . the results with various lattice sizes and control parameters are displayed in fig . [ figure4 ] , indicating that almost identical results are obtained for lattices of 6 and 8 sites in the absence of off - diagonal coupling . hence , in such a scenario , calculations with 6 sites are expected to give sufficiently reliable values of d at affordable computational cost . as @xmath110 increases , however , the number of sites becomes an issue in determining the value of d due to excitation - phonon correlation spans over longer distances . the effect of the lattice size is more prominent in the low temperature regime as shown in figs . [ figure4](b ) and ( c ) , where the transport properties are dominated by coherent or band - like contributions . at higher temperatures , the hopping contribution eventually dominates the transport and the results become less sensitive to the system size . this can also be understood from polaron band narrowing at high temperatures and a reduced effective transfer integral with increasing temperature . as the diagonal coupling strength @xmath116 increases , the difference between the results obtained with 4 sites and 6 sites become less obvious than the weak coupling cases . even better agreement is reached for @xmath110 = 0 , and the results are found to be almost identical from 4 to 8 sites as shown in fig . [ figure4](d ) . it is also interesting to spot nonmonotonic changes of the diffusion coefficient ( i.e. , a `` hump '' ) as a function of the temperature , in the absence of any size dependence , in the low temperature regime for@xmath125 and nonzero @xmath110 , as shown in figs . [ figure4 ] ( e ) and ( f ) . such a hump also appears in fig . [ figure4 ] ( c ) for @xmath107 and @xmath126 despite a slightly visible size dependence of the diffusion coefficient . for a fixed value of @xmath106 , the `` hump '' becomes more pronounced with increasing @xmath110 , underlying the partial role of off - diagonal coupling .     on ( a ) total diffusion coefficient , ( b ) band - like contribution and ( c ) hopping contribution . the other control parameters are @xmath106=0.5 , @xmath110=0.0 and j=0.1 . ]    as shown in fig . [ figure5 ] , increasing the phonon bandwidth @xmath96 leads to an increase in the band - like contribution to the diffusion coefficient , and a decrease of the relative importance of the hopping contribution . these findings from our numerical calculations for a case of diagonal exciton - phonon coupling are in good agreement with analytical results from ref . qualitatively similar results can be obtained with the inclusion of off - diagonal coupling .        in the original work of munn and silbey , the diffusion coefficient as a function of temperature have been studied with structureless a - matrix parameters @xmath57 and @xmath58 . the detailed structure for the a - matrix parameters is found to have a substantial impact on the transport properties , as revealed in fig . [ figure6 ] by comparing the diffusion coefficients from ref . @xcite with those obtained here . moreover , we would like to note that in the absence of off - diagonal coupling our results are in good agreement with those from refs . @xcite and @xcite . the results in this paper can be compared with those obtained using other methods such as a green s function approach@xcite and variational exact diagonalization with a better construction of phonon states@xcite . for example , with an improved numerical technique in ref . @xcite , the 1d polaron mass and radius are studied for an expanded parameter regime , where comparisons with our approach are especially convenient . despite the discrepancies between our results and those of munn and silbey in high temperature regime where hopping transport dominates , agreement is found in low temperature regime where band - like transport prevails . this can be understood as follows . as discussed in ref . @xcite , the transformation coefficients obtained self - consistently are shown to have stronger wave vector dependencies as @xmath119 increases . for weak off - diagonal coupling , however , the average values of @xmath60 and @xmath61 generated by the self - consistency procedure are in reasonable agreement with those by munn and silbey @xcite despite some systematic deviations , which leads to the aforementioned low - temperature agreement . as shown in figs . [ figure2 ] and [ figure3 ] , a higher temperature induces a higher amplitude of wave vector variations of the a - matrix parameters together with a much flatter polaron band as compared with its munn - silbey counterpart , resulting in a diminished diffusion coefficient . finally , the diffusion coefficients obtained here in the absence of the off - diagonal coupling are shown to have qualitative agreements with those from ref . the differences between the two can be ascribed to a series approximations made in ref . @xcite leading to the calculation of the diffusion coefficient , which in this work are substituted with numerically exact evaluations . within the framework of the munn - silbey approach to polaronic transport in organic molecular crystals , diffusion coefficients have been computed with numerical means while taking into account simultaneous diagonal and off - diagonal exciton - phonon coupling for a wide range of temperatures . with the transformation coefficients determined self - consistently without the additional approximations found in refs . @xcite , calculated transport properties exhibit substantial corrections on those obtained previously thanks to the added momentum - dependence of the a - matrix parameters and numerically exact evaluation of many essential quantities . the general transport properties obtained are also in good qualitative agreement with those obtained with the holstein - peierls model and ab - initio calculations@xcite demonstrating robustness of the current approach . comparisons between this work and the one in ref . @xcite reveal the importance of the electron - phonon correlations and the polaron band structure to transport properties of the extended holstein molecular crystal model . by treating the diagonal and off - diagonal exciton - phonon coupling on an equal footing in our self - consistency procedure , the current approach is one of the few capable to model realistically the transport phenomenon in organic molecular crystals , for which the effect of off - diagonal coupling often can not be neglected . off - diagonal coupling , in itself , is a transport mechanism as it is an agent for polaronic localization . inclusion of off - diagonal coupling in the current approach allows us to see the critical role it plays in determining both polaron structures and transport properties over a wide range of temperatures . of course , just as any other treatment of transport , the current approach can be improved , especially in the low temperature regime , where the temperature - independent mobility has recently been found . work in this direction is in progress . support from the singapore ministry of education through the academic research fund ( tier 2 ) under project no . t207b1214 and the singapore national research foundation through the competitive research programme ( crp ) under project no . nrf - crp5 - 2009 - 04 is gratefully acknowledged . the expression of the correlation functions @xmath82 is mainly composed of three parts as : @xmath127 first of all , the expression of @xmath128 is given by : @xmath129,\end{aligned}\ ] ] where @xmath130 and the expression for @xmath131 and @xmath132 has been given by eq . [ fourtheta ] and eq . [ twotheta ] , respectively . here and onwards , @xmath133 represents the total number of sites used in the transport calculation . + now we turn to the expression of @xmath134 , which is written as : @xmath135 \nonumber\\ & + & \frac{1}{\sqrt{n}}[\underset{q_1}{\sum } j_{k_4}\omega _ { q_1}f_{-k_3}^{q_1}\cdot\langle \psi _ { q_1}z_{q_1,k_4;k_4,q_2}(t)\rangle \langle t_{k_3+q_1,k_1;k_3,k_2}t_{k_4,q_1;k_4,q_2}(t)\rangle -\frac{1}{n}\underset{q_1,q_2}{\sum } \omega_{q_1}\omega_{q_2}f_{-k_3}^{q_1}f_{-k_4}^{q_2}\nonumber\\ & \times&(a_{-q_1}^{-q_2}\langle \psi _ { q_1}z_{q_1,k_4+q_2;k_4,q_2-q_2}(t)\rangle\langle t_{k_3+q_1,k_1;k_3,k_2}t_{k_4+q_2,q_1;k_4,q_2-q_2}(t)\rangle\nonumber\\ & + & a_{-q_1}^{q_2\ast } \langle \psi _ { q_1}z_{q_1+q_2,k_4+q_2;k_4,q_2}(t)\rangle\langle t_{k_3+q_1,k_1;k_3,k_2}t_{k_4+q_2,q_1+q_2;k_4,q_2}(t)\rangle)]\},\end{aligned}\ ] ] where @xmath136\delta _ { -k_1+k_2,q}-a_{-k_4}^{k_3-k_4}[ne^{-i\omega t}-(n+1)e^{i\omega t}]\delta _ { -k_3+k_4,q}\},\ ] ] and @xmath137\delta_{q ,- k_1+k_2}-a_{-k_4}^{k_3-k_4}[(n+1)e^{i\omega t}-ne^{-i\omega t}]\delta _ { q ,- k_3+k_4}\},\\\ ] ] here @xmath6 is the bose - einstein distribution , with @xmath138 . + finally we give the expression of @xmath139 :
Abstract:

Answer: improved results using a method similar to the munn - silbey approach have been obtained on the temperature dependence of transport properties of an extended holstein model incorporating simultaneous diagonal and off - diagonal exciton - phonon coupling . the hamiltonian is partially diagonalized by a canonical transformation , and optimal transformation coefficients are determined in a self - consistent manner . calculated transport properties exhibit substantial corrections on those obtained previously by munn and silbey for a wide range of temperatures thanks to a numerically exact evaluation and an added momentum - dependence of the transformation matrix . results on the diffusion coefficient in the moderate and weak coupling regime show distinct band - like and hopping - like transport features as a function of temperature .