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: precise determination of @xmath2 , a parameter of the cabibbo - kobayashi - maskawa ( ckm ) matrix , is an important goal of many experimental studies . the current experimental limit @xcite @xmath3 is based on measurements of the beauty hadron decays produced at the @xmath4 resonance ( by argus and by cleo ii ) and in @xmath5-boson decays ( by the four experiments at lep ) . in the future large samples of the @xmath1-hadrons collected at @xmath6-factories ( at slac and kek ) and at the hadron colliders will increase the statistical accuracy to a few percent level . to fully exploit the anticipated experimental improvement , the theoretical description of the @xmath1 decay must be known with comparable precision . there are two methods of extracting the value of @xmath2 , based on measurements of the exclusive decay @xmath7 and of the inclusive semileptonic decay width of @xmath1-hadrons @xmath8 . these two methods rely on very different theoretical considerations and experimental procedures and complement each other . their merits and theoretical uncertainties are summarized e.g.  in ref . one of the major sources of the theoretical error are the perturbative qcd corrections at the two loop level . for the exclusive decays at the zero recoil point these corrections have recently been calculated @xcite . this has significantly improved the accuracy of the theoretical prediction for the exclusive method .    regarding the inclusive method , recently the @xmath0 corrections to differential semileptonic decay width of the @xmath1quark at maximal recoil have been calculated @xcite . combined with the previously obtained value at zero ( minimal ) recoil , these results permitted an estimate of the complete @xmath0 correction to the total semileptonic decay width of the @xmath1quark .    in the paper @xcite we presented the results of that calculation and discussed its phenomenological relevance . the purpose of the present paper is a detailed description of the methods employed in that calculation . we would like to note that a complete calculation of two - loop corrections to a fermion decay width has never been performed before , neither in qcd nor in qed ( a longstanding example are the two - loop qed corrections to the muon life time ) . therefore , in the calculation we describe here we had to go beyond the traditional methods used in higher  orders calculation . we hope that a description of some technical aspects and methods will be of interest for the community . let us first mention the difficulties one encounters when trying to compute the fourth order corrections to the fermion decays . one of the problems is an appropriate treatment of the real radiation of one or two gluons . for the virtual radiative corrections there exists a number of algorithms , permitting an efficient , analytical treatment of a large number of complicated diagrams , which typically appear in such calculations . on the contrary , no similar algorithms were available so far for the treatment of the real radiation . the reason why the real radiation at order @xmath9 is difficult to evaluate is that the particle in the initial state ( the decaying @xmath1 quark ) carries a color charge and therefore can radiate . it is the presence of the massive propagator of the initial quark which makes the integrations over the phase space very tough . the kinematical configuration , where the invariant mass of the leptons is equal to zero and the quark in the final state is massive is the first case where a complete analytical evaluation of the real radiation of two gluons in a decay of a fermion turns out possible . another potential problem is the treatment of diagrams which represent one  loop virtual corrections to a single gluon emission in the @xmath1quark decay . the virtual corrections in such situation are one  loop . therefore one might naively expect that this case does not require any sophisticated investigation . unfortunately , the integration of the one  loop formulas ( especially of the boxes ) over the three body phase space is difficult . however , it turns out possible to express the result for the loop as an expansion which can be easily integrated over the phase space . with a systematic algorithm for the expansion , this approach shifts the burden of the calculation to the computer . the idea which permitted us to calculate the contribution of the real radiation of one or two gluons is ( qualitatively speaking ) the expansion in the velocity of the final quark . in the limit @xmath10 the charm quark in the final state is a slowly moving particle , with spatial components of its momentum of the order of @xmath11 , much smaller than its mass . the momenta of gluons and of leptons ( in the case when the invariant mass of the leptons is zero ) are also of the order of @xmath12 . it turns out that by a proper choice of the phase space variables one can systematically expand the amplitudes and the phase space in @xmath13 ( sometimes we also use an equivalent expansion parameter @xmath14 ; throughout this paper we use @xmath15 as the unit of mass , putting @xmath16 ) . the details of the phase space parameterization will be explained in detail . a word of caution is in order here . our technique proved to be very useful in problems where the mass of the quark in the final state does not differ too much from the mass of the quark in the initial state . in this respect an ideal application are semileptonic @xmath17 transitions . on the other hand , for such problems as muon or top quark decay , the expansion parameter may be close to unity and it is not clear at present if our procedure is of any use there . we will , however , show an example where our procedure remains meaningful and delivers reliable predictions even in the case when the expansion parameter equals unity .    the paper is organized as follows . in the next section we discuss how the real radiation of two gluons can be computed . section [ sec : tensors ] is devoted to the treatment of tensor integrals . then we present a detailed study of one  loop virtual corrections to a single gluon emission in the @xmath1decay . we show that in such diagrams a new type of feynman integrals appears and discuss their evaluation . section [ sec : applic ] is devoted to the applications of our techniques . we present the result for the @xmath18 correction to the differential semileptonic decay width @xmath19 up to @xmath20 . we also demonstrate how the known light quark corrections to the top quark decay width can be obtained with good accuracy if the expansion up to high powers of @xmath21 is performed , and discuss an application of our technique to corrections to the zero recoil sum rules . in the last section we present conclusions . we first consider the radiation of two real gluons in a semileptonic @xmath1decay @xmath23 and discuss how its contribution to the width can be computed .    throughout this paper we denote the particles and their four - momenta by the same letters , i.e.  the momentum of the @xmath1quark is @xmath1 etc . moreover , we denote the momentum of the lepton pair @xmath24 as @xmath25 and often speak about the decay @xmath26 . all formulas in this paper apply to the case when the invariant mass of the leptons is zero , i.e. @xmath27 , if not stated otherwise . our aim is to construct an expansion in @xmath28 . in this case the spatial momentum of the @xmath29 quark is of the order of @xmath12 . the momenta of the lepton pair and of the gluons are also of the order of @xmath12 because @xmath30 . the smallness of these momenta permits an expansion in terms of @xmath31 .      to show that such an expansion is indeed possible , we list here all propagators of the virtual particles which appear in a calculation of the two gluon emission in the semileptonic @xmath1decay : @xmath32 all these propagators are shown in the examples in fig . [ fig : propagators ] . if we eliminate the momentum of the @xmath29 quark using momentum conservation , @xmath33 , we can expand the propagators @xmath34 in the taylor series with respect to the `` small '' momenta @xmath35 , @xmath36 , @xmath25 . explicitly , we have @xmath37^j}{[2b(g_1+g_2)-2g_1g_2]^j}.\end{aligned}\ ] ] from these expressions we see that the expansion does not generate any new infrared divergences and , therefore , appears to be permissible . after this expansion only four different types of denominators remain . we list them for clarity : @xmath38 as will become clear from the discussion below , it is more convenient to use @xmath39 as one of the basic propagators than to expand it any further . let us emphasize at this point the advantages of the above expansion . as can be seen from the above expressions , the four  momentum of lepton pair @xmath25 and the momentum of the @xmath29quark do not appear in the denominator . this immediately implies that the integral over the phase space of @xmath25 and @xmath29 factorizes . therefore , the expansion suggests a simple way how the rather nontrivial integration over the phase space of four particles in the final state can be reduced to a more familiar case of the three body phase space . this reduction is described in the next section .      after we have checked that the propagators of all virtual particles can be expanded , the same must be demonstrated for the phase space element this is done in the present section . considering the semileptonic decay of the @xmath1quark and integrating first over the lepton phase space , one obtains the following phase space integration element ( we consider only the point where the invariant mass of the lepton pair @xmath40 is zero ) : @xmath41[dw][dg_1][dg_2 ] \delta^d ( b - c - g_1-g_2-w),\ ] ] where a shorthand notation is @xmath42 = \frac { { \rm d}^{d-1}p}{2p^0},\qquad d=4 - 2\epsilon.\ ] ] it is convenient to introduce an auxiliary vector @xmath43 which equals to the sum of the momenta of the @xmath29quark and @xmath25 boson , @xmath44 . with this notation the phase space integral becomes @xmath45 [ dg_1][dg_2 ]   \delta^d ( b - h - g_1-g_2 ) \int [ dw][dc ] \delta^d ( h - w -c).\ ] ] the integration over the @xmath46 phase space gives ( for now we assume the integrand contains no @xmath25 in the numerator ; the tensor integrals with @xmath25 dependence will be analyzed in section [ sec : tensors ] ) @xmath47[dc ] \delta^d ( h - w -c)=   \frac { \omega _ { d-1}}{2^{d-1 } } \frac { ( h^2-m_c^2)^{d-3}}{(h^2)^{d/2 - 1}}\ ] ] where @xmath48 is the volume of a @xmath49 dimensional sphere of a unit radius . having performed the integration over the @xmath46 phase space we are left with the three particle phase space of one massive ( @xmath43 ) and two massless ( @xmath35 , @xmath36 ) particles and the integration over the square of the momentum @xmath43 . we use the variables ( with @xmath50 ) @xmath51 and get the following expression for the four particle phase space @xcite @xmath52 with @xmath53 and the following limits of integrations @xmath54 in order to expand the phase space element in powers of @xmath55 it is useful to change the variables in the above expression @xmath56 the integration limits for all @xmath57 are the same : @xmath58 expressing the phase  space element through these variables and restoring proper dimensionality we find : @xmath59,\ ] ] where @xmath60    our final aim is to expand @xmath61 in the limit of @xmath62 . the above expression is well suited for this purpose . expanding it in powers of @xmath63 gives rise to simple integrals which can be expressed by beta functions : @xmath64      the basic propagators @xmath65 , @xmath66 , @xmath39 , and @xmath67 , expressed in terms of the variables @xmath68 yield simple expressions which can be integrated over the phase space . we find @xmath69 } , \nonumber \\ p_3 ^ 0 & = & \frac { 1}{2b(g_1+g_2)-2g_1g_2 } = \frac { 1}{\beta x_1},\qquad   p_7 = \frac { 1}{2g_1g_2 } =   \frac { 1-\beta x_1 x_2}{\beta^2 x_1 ^ 2 x_2 ( 1-x_2 ) x_3}.\end{aligned}\ ] ] clearly , the expansion of these quantities in terms of small @xmath70 is readily performed . the resulting integrals are very similar to the integrals in eq . ( [ eq : master ] ) and can be expressed in terms of the beta functions . the same is true for all scalar products of the four momenta which enter the calculation . therefore , the above discussion demonstrates that the part of the semileptonic decay width of the @xmath1quark containing radiation of two gluons can be calculated by expanding the matrix element and the four particle phase space in powers of @xmath70 . what is still missing in our discussion is the study of the integrals over the @xmath46 phase space when the numerator depends on @xmath25 . we discuss this issue in the next section . the expansion of the propagators in terms of the momenta @xmath25 , @xmath35 , @xmath36 , may result in high powers of the scalar products @xmath71 , which have to be integrated over the @xmath46 phase space . in this section we present efficient methods for such integrations . the results of this section are not restricted to the case @xmath27 and are likely to be useful for solving other problems .    the main object we are going to discuss is the following integral over the @xmath46 phase space : @xmath72[dc ] \delta^d ( h - w -c ) ( wg_1)^a ( wg_2)^b   \nonumber \\ & \equiv & \langle   ( wg_1)^a ( wg_2)^b \rangle _ { h = w+c}\end{aligned}\ ] ] for arbitrary @xmath73 . this integral can be rewritten as @xmath74 the important property of the tensor @xmath75 is that it depends on a single external vector only ( @xmath43 ) . this property simplifies the calculation of the necessary integrals . we discuss below two approaches to such calculation ; their relative merits depend on the capabilities of the symbolic manipulation languages , if the algorithm is to be implemented using a computer . the first approach is more general and potentially more efficient . to fully exploit the dependence on the single external momentum we can extract the component of the vector @xmath25 which is parallel to it . for this purpose we write : @xmath76 the value of @xmath77 is then fixed , because in the general case @xmath78 we have @xmath79 therefore , if we use the above substitution for @xmath25 in the tensor integral @xmath80 , we obtain a sum of tensor integrals over transverse components of the @xmath25 ( with respect to @xmath43 ) : @xmath81 evidently , the tensor structure of this integral is trivial : since @xmath82 , @xmath83 must be proportional to the _ absolutely symmetric tensor _ of the @xmath84 dimensional space . the `` building block '' is the metric tensor of @xmath84 dimensional space : @xmath85 . we note that only tensors of _ even _ rank contribute . the generation of the absolutely symmetric tensor of @xmath49 dimensional space can be easily encoded in a symbolic manipulation program . hence , the algorithm described above permits an efficient treatment of the tensor integrals which appear in the problem at hand .      using the same idea as in the previous section we can advance slightly further with the analytical calculation of the integral @xmath86 .    after writing @xmath87 and noticing that @xmath88 and @xmath89 are constants in the @xmath46 phase space , we conclude that the non - trivial integrals to be computed are @xmath90 with @xmath91 ( i=1,2 ) defined in analogy to @xmath92 . the integral @xmath93 is a lorentz  scalar and we can choose an arbitrary frame for its calculation . for instance , we consider @xmath43 to have only the zeroth component and put the @xmath94axis of the @xmath49 dimensional space along the vector @xmath95 . then , rewriting @xmath96 as @xmath97 we find that @xmath98 is expressed as a sum of integrals of the from : @xmath99 this integral is non - zero only if both @xmath100 and @xmath1 are even integers .    since @xmath101 , we can put the @xmath102-axis of the @xmath49 dimensional space along the vector @xmath103 . then we have : @xmath104 the integration element is @xmath105 if we choose @xmath94 and @xmath102 axes as described above , we get @xmath106 the minus sign in the above equations arises because the transverse vectors here are spacelike . in calculating @xmath107 the following formula is useful : @xmath108 this formula is valid for even and positive @xmath70 . using it we get the result for the integral @xmath107 : @xmath109 in order to use these results for the further integration over the phase space @xmath110 it is necessary to express the modulus of the vectors @xmath92 , @xmath95 , @xmath103 in terms of the variables @xmath111 introduced in the previous section . this is a cumbersome but straightforward task and we do not discuss it here in any detail . the steps described above allow an easy implementation in any symbolic manipulation program . they provide an efficient and uniform treatment of the tensor integrals which appear in this calculation . another source of the @xmath112 corrections are the one  loop virtual corrections accompanied by a radiation of an additional real gluon . again , we would like to expand those diagrams in powers of @xmath113 . the purpose of this section is to provide a suitable algorithm . there is a principal difference from the radiation of two real gluons in this part of the corrections . namely , it is not sufficient here to perform a taylor expansion of the relevant diagrams .    to see this , let us analyze a simple example . we consider a diagram where the @xmath1quark radiates a gluon and then a one  loop self  energy correction is inserted on the @xmath1 quark line ( see fig . [ fig : simplest ] ) . taking the limit @xmath114 in such situation implies that the momentum of the real gluon goes to zero . therefore the self  energy diagram should be calculated close to the mass shell of the @xmath1quark . it is well  known that in this limit the self  energy diagram has an on - shell logarithmic singularity of the form @xmath115    since this diagram contains @xmath116 it is clear that it can not be calculated as a taylor expansion with respect to @xmath31 . therefore , a more sophisticated procedure is needed and the method of _ eikonal expansions _ @xcite is used in this case . we illustrate this approach in the following example .      in this subsection we consider a simple example of a one  loop integral , where a result obtained from the eikonal expansion can be compared with an exact formula . we calculate the scalar self  energy diagram : @xmath117}. \label { example}\ ] ] performing feynman parameterization , integrating over the loop momenta , and expanding in @xmath118 we get ( using @xmath119 here ) : @xmath120    the appearance of the @xmath121 term in this result signals that the taylor expansion of the _ integrand _ will be insufficient . in order to get the correct result we add the eikonal expansion . for this purpose we write @xmath122 the taylor expansion term @xmath123 is obtained by expanding the integrand in ( [ example ] ) in a taylor series in @xmath124 . this yields a set of one  loop on  shell integrals which are well known . the validity of such an expansion is determined by the condition @xmath125 . if this is not the case , the taylor expansion breaks down . this breakdown gives rise to an infrared singularity which occurs in the region where @xmath126 and @xmath127 . this spurious divergence can be cancelled by adding an expansion of the integrand in @xmath128 ( for a more detailed discussion see @xcite ) . one should note in addition that if any power of @xmath128 appears in the numerator we get an integral of the form : @xmath129 an important point is that such integrals are zero within the dimensional regularization framework . one simple reason for this is that the integral over transverse components of @xmath130 is scaleless . hence , the only term in the eikonal expansion we should consider is the one with @xmath128 neglected in the second propagator . the eikonal integral therefore reads : @xmath131 it is evident that this integral scales with @xmath124 as @xmath132 . it is this dependence on @xmath133 which gives @xmath134 in the final result . consider @xmath135 in the rest frame of @xmath136 . performing the integration over the time component of the loop momenta first , one gets a simple integral representation over the transverse components of the loop momenta which can be easily calculated . more sophisticated eikonal integrals are considered in the next section , where all necessary formulas can be found .      here we provide a set of formulas for dealing with integrals appearing in the eikonal expansions . as is clear from the example presented above , the region of integration in which we are interested here is characterized by the values of the virtual momentum of the order of the momentum of the real gluon emitted in the decay . we denote the loop momentum by @xmath130 and the momentum of the real gluon by @xmath137 ; explicitly , the eikonal expansion of the propagator @xmath138 is given by @xmath139    the basic eikonal integral in the one  loop correction then reads . ] @xmath140    in the actual calculation @xmath141 is always a _ time  like _ vector which for different diagrams can be @xmath142 .    as we pointed out above , the integrals with any power of @xmath128 in the numerator vanish . there are several ways to see this . first , the integrals over transverse ( with respect to @xmath141 ) components of the loop momenta are scaleless ; also if there is no @xmath128 in the denominator , the poles of the integrand are located only on one side of the @xmath143 integration axis . also , if @xmath144 , the integral in eq . ( [ eik ] ) is zero because it is scaleless . let us describe the most economic way to calculate @xmath145 . we choose the lorentz frame in which @xmath146 . the integration over the time - like component of the vector @xmath130 can be performed using the residues . if the integration contour is closed in the upper half  plane , only that pole in @xmath147 contributes for which @xmath148 . after @xmath149 , the integration over the @xmath49 dimensional space should be performed . to discuss the integration over @xmath150 we first note that after integrating over @xmath149 the denominator of the integral @xmath151 depends on @xmath152 only . therefore , all tensor integrals which appear due to the numerator structure are simply related to the @xmath49 absolutely symmetric tensor of the corresponding rank . hence , the tensor integrals here are similar to the ones discussed in section [ sec : tensors ] . finally , we are left with the radial integration over the modulus of @xmath153 . the following formulas are useful : @xmath154,\end{aligned}\ ] ] with @xmath155 in the above formulas @xmath156 is the volume of the @xmath49 dimensional sphere of unit radius . there is one subtle point concerning the expressions presented above . namely , it is easy to see that the integrals @xmath157 have an additional singularity after the integration over @xmath149 has been performed . this singularity corresponds to the appearance of the imaginary part in some of feynman diagrams which contribute to the result . below we explain how this singularity was treated .    in the singular integrals the integration over the modulus of the transverse momenta is @xmath158 the singularity is located at the point @xmath159 . to make this integral meaningful we rewrite it as @xmath160 rescaling @xmath161 we get : @xmath162 for the derivative we use @xmath163 finally , after performing the remaining integration over @xmath150 using @xmath164 we arrive at the formula for @xmath165 quoted above .      the graphs with a triple gluon coupling lead to the most difficult eikonal integrals . the difficulty originates from the fact that in such graphs one has two massless ( gluon ) propagators . a typical integral reads @xmath166    the taylor expansion of such integral is easy . however , in the eikonal expansion one can not expand massless propagators . therefore , the master integral in this case is @xmath167 where @xmath141 is again a `` large '' time - like vector , as described in the previous section . note also , that in this case the terms with @xmath128 or @xmath168 in the numerator do contribute to the result . let us describe in detail how such integrals can be calculated . for this purpose we first combine the massless denominators : @xmath169 next , we decompose the product of two eikonal propagators @xmath170 using partial fractions . the integral @xmath171 becomes now a sum of the integrals of the following type : @xmath172 ^ 2 ( 2kq-2gq)^n } , \nonumber \\ i_2(n)&=&\int _ { 0}^{1}~{\rm d } x~\int \frac { { \rm d}^dk}{(2\pi)^d }   \frac { 1}{[(k - gx)^2]^2 ( 2kq)^n}.\end{aligned}\ ] ]    in both integrals we shift the integration momentum @xmath173 and get @xmath174^n } , \nonumber \\ i_2(n)&=&\int _ { 0}^{1}~{\rm d } x~\int \frac { { \rm d}^dk}{(2\pi)^d }   \frac { 1}{(k^2)^2 ( 2kq+2gqx)^n}.\end{aligned}\ ] ] now , in the first integral one makes the change of variables @xmath175 and then rescales the integration momenta @xmath176 in both @xmath177 and @xmath178 . the integration over @xmath102 factorizes and can be done using @xmath179    the integrals over @xmath130 differ from the integral in the previous section by the presence of two powers of @xmath128 in the denominator . in this case it is useful to perform a transformation ( similar to the technique of integration by parts @xcite ) , which reduces the integrals with the second power of @xmath128 in the denominator to the integrals with the first power only .    to write down this transformation we introduce the notation @xmath180 then @xmath181.\ ] ] using this formula we end up with integrals with the first power of @xmath128 in the denominator , for which the formulas of the previous section are applicable .      in the treatment of the radiation of one gluon with @xmath182 accuracy , the last step one has to do is to perform the integration over the phase space of the real gluon . the eikonal expansion , described above in detail , gives the virtual corrections in terms of powers and logarithms of the small parameters . this form of the intermediate result simplifies the final integration over the phase space of the decay products . this would not be the case if we had an _ exact _ result for the loop : that would certainly contain dilogarithms of complicated arguments .    the phase space integration with a single gluon is in its general structure similar to the case of two real gluons , discussed in section [ sec : pstwo ] . for completeness , we give a short account of necessary formulas .    in the present case the phase space element is : @xmath183[dw][dg]\delta^d ( b - c - g - w ) , \nonumber \\ & = & \left(\frac { \omega _ { d-1}}{2^{d-1 } } \right)^2 m_b^{2 - 4\epsilon } \beta ^{3 - 4\epsilon } \int _ { 0}^{1 } { \rm d } x j(x ) \frac { 1}{(1-\beta x)^{1-\epsilon } } \nonumber \\ j(x ) & = & x^{1 - 2\epsilon}(1-x)^{1 - 2\epsilon}\end{aligned}\ ] ] with @xmath184    the master integral in this case reads : @xmath185 the results of the eikonal integrals require a slight modification of the master integral . as we saw in section [ eikogen ] , those results contain @xmath186 which results in an extra factor @xmath187 . in this case the master integral is @xmath188 in this section we discuss some applications of the techniques described above . first , we give a complete result for the @xmath189 correction to the differential semileptonic decay width of the @xmath1quark at the maximal recoil point . we next re - derive the blm @xmath0 corrections to the quark top decay width into massless @xmath25boson and @xmath1 quark with the help of the expansions presented above . as will be explained below , the expansion parameter in this case is unity . it not obvious that the techniques described above can be of any use there . therefore it is interesting to check how the procedure works in such extreme limit . finally , we comment on the applications to zero recoil sum rules . the measurement of the inclusive semileptonic decay width permits a determination of the ckm matrix parameter @xmath2 with a small theoretical uncertainty @xcite . the magnitude of the perturbative corrections @xmath0 to this quantity has been subject of discussions in the recent literature . we addressed this problem in ref . @xcite where the exact calculation of the @xmath0 corrections at maximal recoil was used to estimate the total 4th order correction to the @xmath1 quark semileptonic width . here we present more complete results of that calculation . we consider semileptonic decay of the @xmath1quark @xmath190 . the momentum carried away by leptons is denoted by @xmath191 . we write differential semileptonic decay width of the decay @xmath192 at @xmath193 as @xmath194_{q^2=0}\;= \;\gamma _ 0 \left[\delta _ { \rm born } +   \frac { \alpha _ s}{\pi}c_f \delta _ 1 + \left(\frac { \alpha _ s}{\pi}\right)^2 c_f \delta _ 2 \right ] \label{eq : param}\ ] ] where @xmath195 . @xmath196 and @xmath197 is also known in a closed analytical form @xcite . @xmath198 is the new result which gives @xmath0 correction to the differential semileptonic decay width of the @xmath1-quark at the point of zero invariant mass of leptons . for the purpose of the presentation we write the 4th order correction as @xmath199 . \label{eq : deltas}\ ] ]    in the above equation @xmath200 is the expansion parameter . @xmath201 describes the contribution of the massive @xmath1 and @xmath29 quark loops .    for the @xmath202 group the color factors are @xmath203 , @xmath204 , @xmath205 . @xmath206 is the number of light quark flavors whose masses were neglected .    in ref . @xcite we have presented analytical results for the functions @xmath207 up to @xmath208 , though actual calculations have been done up to @xmath209 . our complete formulas for these functions are presented in appendix [ app:8power ] .      as another application of the above techniques we consider the decay of the top quark @xmath211 . it is well known that this decay width ( at least at the born and the one  loop level ) can be well approximated by neglecting the masses of the @xmath25boson and the @xmath1quark . in such case , since @xmath212 , the techniques presented in the main part of this paper can be readily applied . in particular , the formula ( [ eq : param ] ) can be rewritten in such form that it gives the correction to the two - body decay width @xmath213 ( where @xmath141 is a heavy quark ) : @xmath214.\ ] ] the @xmath215 functions here are the same as in eq . ( [ eq : param ] ) , with @xmath15 and @xmath216 replaced by @xmath217 and @xmath218 . the problem , however , is that the procedures described above were based on the expansion of the rate in the mass difference of the final and initial quarks . in the case of the top quark decay with a @xmath1 quark in the final state this means that the expansion parameter @xmath219 is close to 1 for realistic values of @xmath217 and @xmath15 . little is known about convergence properties of the series described in this paper . fortunately , part of the @xmath0 correction to the top quark width is known exactly ; it is the contribution of the massless quarks , calculated for @xmath220 in @xcite ( the relevant diagrams are shown if fig . [ fig : blm ] ) . we can use that limiting case to check if we can reproduce it with our techniques . the exact formula for the massless quark correction reads ( in the @xmath221scheme and for one generation of massless fermions ) @xmath222 , \nonumber \\ \widetilde\gamma _ 0 & = & \frac { g_f~m_t^3}{8\ ; \sqrt { 2}\;\pi}.\end{aligned}\ ] ]    it should be mentioned that the diagrams with real or virtual massless fermions represent the simplest case for the described algorithms . the reason is their simple planar topology , which allows the computer programs to work very fast . for this particular type of diagrams we can expand the width up to a very high power in @xmath31 . we write the result of this expansion as @xmath223 where the function @xmath224 is given as a series of powers and logs of @xmath31 ; we have computed this expansion up to terms of the order @xmath225 . the first few terms of that expansion are given by the function @xmath226 in the appendix . the numerical value of the exact result is @xmath227 the values of the approximate result @xmath224 for @xmath228 are , for three numbers @xmath229 of the summed terms , @xmath230    comparing these numbers with the exact result we see that already the first 11 terms of the expansion give the accuracy of about @xmath231 . unfortunately , the accuracy does not improve significantly with the growing number of terms in the expansion . this slow convergence of the series is caused by the logarithms of the small ( @xmath1 quark ) mass . although these terms are suppressed by powers of that mass ( and should give no contribution to the final result if @xmath232 ) , they spoil the convergence of the expansion . it would be very useful ( and we think it is possible ) to find a systematic way of eliminating those parts of the integrands which give rise to these logs .    on the other hand , from the perspective of practical applications ( like top physics at the next linear collider ) it would be sufficient to know the two  loop qcd correction to the top quark decay width even with @xmath233 percent accuracy . it is tempting to apply our algorithms to this problem and extrapolate the result to the point @xmath234 ; on the other hand , at present we do not know any reliable method for estimating the uncertainty of such result . another useful application of the techniques described above is connected with the so - called zero recoil ( zr ) sum rules @xcite . in this subsection we would like to discuss this point . a more detailed discussion of zr sum rules with the @xmath0 corrections will be given in a future publication @xcite . the zr sum rules are important for the estimates of the zero recoil transition form factors such as @xmath235 for @xmath236 transition or @xmath237 for @xmath238 transition . in turn , the formfactor @xmath235 is a crucial ingredient for @xmath2 determination from the exclusive semileptonic @xmath239 decays .    the sum rules for exclusive heavy  to  heavy flavor transitions are based on the operator product expansion ( ope ) of the hadronic amplitudes in terms of the inverse quark masses . the zeroth order term in this expansion is the parton model where free quarks are substituted for real hadrons both in initial and final state . here we disregard the non - perturbative corrections and discuss how the perturbative corrections @xmath0 to the sum rules can be evaluated . we consider a transition of a @xmath1 quark at rest to a @xmath29 quark and massless partons which occurs under the influence of the external current @xmath240 . the momentum carried away by the external current is @xmath241 , where @xmath242    the quantity which is of primary importance for the perturbative corrections to zr sum rules can be schematically written as : @xmath243 both on the external momenta transfer @xmath244 and on the lorentz structure of the current .    if @xmath245 , the transition is elastic ; the final state is a single @xmath29quark at rest . the @xmath0 corrections in this case reduce to the renormalization of the external current @xmath240 . for vector and axial currents they were calculated in @xcite .    on the other hand , if @xmath246 the @xmath29 quark starts moving and can radiate . the second order corrections @xmath247 were calculated in @xcite . aiming at the fourth order , i.e. at @xmath248 accuracy , one has to consider the final state with the @xmath29quark and two real gluons or light quarks and also the @xmath247 correction to the single gluon emission in @xmath17 transition .    due to the hierarchy of scales , @xmath249 , the final @xmath29quark is moving slowly by definition ; this is therefore a nice place where the techniques described in this paper can be applied . in particular , an algorithm for performing eikonal expansions with the subsequent integration over the phase space appears to be very useful here . this topic will be discussed in detail in ref . in this paper we have described techniques used in the calculation of the @xmath0 corrections to the semileptonic decay width @xmath190 at maximal recoil @xcite . the technical tool we used for that calculation is an expansion of the decay rate in powers and logarithms of the mass difference between the initial and final quarks . we presented a detailed discussion of the algorithm , which enables one to construct such expansion . we treated virtual corrections , emission of one gluon , and emission of two gluons separately . therefore , these algorithms can be used for the analyses of less inclusive quantities than the total decay rate , at least in principle .    in case of two  loop virtual corrections and emission of two real gluons , the expansion in @xmath113 is a taylor expansion . in the case of one  loop corrections to the amplitude of single gluon emission in @xmath1-decay , the taylor expansion is insufficient . an appropriate method is provided by the eikonal expansions , recently introduced in refs . @xcite . when this procedure is used , a new type of feynman integrals appears . these integrals and methods which were used for their evaluation were described here in some detail . the whole construction works well if the mass difference between initial and final state quarks is not too large . this is the case for the semileptonic @xmath17 transitions , where the expansion parameter @xmath250 . in this case we calculated the expansion up to the eleventh power of @xmath31 ; the estimated accuracy of the final result is better than @xmath251 . there is , however , a number of other applications where the initial quark is significantly heavier than the final one . it is not obvious to what extend the present method can be useful in such situation . there is an indication , however , that our procedures give meaningful results even in that limit . as an example , we analyzed the light quark corrections to the width of the top quark decay into massless @xmath25 boson and a @xmath1-quark . we have shown that the first several terms of the expansion in @xmath252 approximate the known exact result with a @xmath253 accuracy . finally , we have argued that the same techniques can be applied to the corrections to the zero recoil sum rules for the heavy flavor transitions . we are grateful to k.  g.  chetyrkin and n.  g.  uraltsev for many helpful discussions . we would like to thank prof . j.  h.  khn for his interest in this work and support . this research has been supported by bmbf 057ka92p and by graduiertenkolleg `` teilchenphysik '' at the university of karlsruhe . 10    , phys . rev . * d54 * , 1 ( 1996 ) . m. shifman , n.  g. uraltsev , and a. vainshtein , phys . rev . * d51 * , 2217 ( 1995 ) , erratum : ibid . d52 , 3149 ( 1995 ) . m. neubert , preprint hep - ph/9702375 , to appear in a. buras and m. lindner ( eds . ) , _ heavy flavours _ , 2nd edition , ( world scientific , singapore ) . a. czarnecki , phys . lett . * 76 * , 4124 ( 1996 ) . a. czarnecki and k. melnikov , hep - ph/9703277 , submitted to nucl . b ( unpublished ) . a. czarnecki and k. melnikov , phys . lett . * 78 * , 3630 ( 1997 ) . a. czarnecki , in _ new computing techniques in physics research iv _ , edited by b. denby and d. perret - gallix ( world scientific , singapore , 1995 ) , pp . 319323 , proceedings of the international workshop ainhep 1995 , pisa , italy . v.  a. smirnov , phys . b394 * , 205 ( 1997 ) .    a. czarnecki and v.  a. smirnov , phys . lett . * b394 * , 211 ( 1997 ) . k.  g. chetyrkin and f. tkachov , nucl . phys . * b192 * , 159 ( 1981 ) . m. jeabek and j.  h. khn , nucl . phys . * b314 * , 1 ( 1989 ) . a. czarnecki and s. davidson , phys . d * 48 * , 4183 ( 1993 ) . b.  h. smith and m.  b. voloshin , phys . b340 * , 176 ( 1994 ) . a. czarnecki , acta phys . pol . * b26 * , 845 ( 1995 ) . j.  g. krner , k. melnikov , and o. yakovlev , z. phys . * c69 * , 437 ( 1996 ) . a. kapustin , z. ligeti , m.  b. wise , and b. grinstein , phys . lett . * b375 * , 327 ( 1996 ) . a.  czarnecki , k.  melnikov ,  n.  g.  uraltsev , in preparation . in this appendix we present the first nine terms of the expansion of the coefficient functions for the @xmath1 quark decay rate at maximal recoil , as defined in ( [ eq : deltas ] ) : @xmath254 with @xmath255 and @xmath256 .
Abstract:

Answer: we have recently presented the complete @xmath0 corrections to the semileptonic decay width of the @xmath1 quark at maximal recoil . here we discuss various technical aspects of that calculation and further applications of similar methods . in particular , we describe an expansion which facilitates the phase space integrations and the treatment of the mixed real  virtual corrections , for which taylor expansion does not work and the so - called eikonal expansion must be employed . several terms of the expansion are given for the @xmath0 qcd corrections to the differential semileptonic decay width of the @xmath1quark at maximal recoil . we also demonstrate how the light quark loop corrections to the top quark decay rate can be obtained using the same methods . we briefly discuss the application of these techniques to the calculation of the @xmath0 correction to zero recoil sum rules for heavy flavor transitions .


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: poly - diacetylene ( pda ) chains dispersed with low concentration in their monomer single - crystal matrix are prototypical quasi one - dimensional materials . @xcite the structural disorder in the chains and their surrounding matrix is tiny , the electronic excitation energies of the diacetylene monomers are much higher than those of the polymer , and the chains electronic excitations in the energy range of visible light can be measured with a very high accuracy.@xcite    exciton - polaritons have been generated that have been shown to be coherent over tens of micrometers , i.e. , several ten thousand monomer units . @xcite this observation was confirmed by weight measurements after dissolving the chains and their monomer matrix.@xcite consequently , the opto - electronic properties of the pdas result from the electrons mutual interaction and their interaction with the lattice potential , while the influence of disorder is negligible . this makes these materials the perfect testing - ground for theoretical model studies which describe interacting electrons on perfectly ordered chains . the typical single - particle gap in pdas is @xmath6 , see sect . [ sec : exp ] . density - functional theory band structure calculations in the local - density approximation for generic pda geometries estimate the bare band - gap to be @xmath7 , or less . @xcite results from various methods are compiled in table  1 of ref . [ ] ; recent calculations using the perdew - burke - ernzerhof global hybrid density functional and the 6 - 311g(2@xmath8,2@xmath9 ) basis set of atom - centered gaussian functions ( geometries from the tpss density functional and the 6 - 31g(@xmath8 ) basis set )  @xcite give @xmath10 . @xcite the comparison shows that electronic exchange and correlations account for a substantial fraction of the single - particle gap . in contrast to inorganic semiconductors , the exciton binding energy in pdas amounts to about 20% of the single - particle gap . such large binding energies suggest that the electron - electron interaction must be treated accurately for the calculation of the optical properties of the pdas .    in order to describe the optical excitations in pdas , two approaches have been taken . the first approach starts from an ab - initio density - functional theory calculation of the bare band structure in local - density approximation ( lda ) , which is then supplemented by an approximate treatment of the residual electron - electron interaction , e.g. , the @xmath11 approximation for the single - particle bands and the bethe - salpeter equation ( bse ) for the excitons ( lda+@xmath11+bse ) . @xcite actual calculations for the pdas often omit the @xmath11  step ( wannier theory ) . @xcite within this approach , a number of experimental data can be reproduced , e.g. , the @xmath12 exciton binding energy and its polarizability . this approach is less successful for the triplet sector . @xcite typically , the energy of the triplet ground state is too high . recall that , if the electron - electron interaction is absent , the energy of the triplet ground state is identical to the single - particle gap . starting from a weak - coupling description of the electron - electron interaction on the chains , it is difficult to obtain the experimentally observed energy renormalization by a factor of almost three . @xcite moreover , polaronic effects are not considered in the lda+bse approach . the second approach to a theoretical description of the primary excitations in polymers starts from a many - particle model hamiltonian that describes only the @xmath13-electrons and their mutual interaction . typically , empirical parameters are used for the tight - binding band structure and for the pariser - parr  pople ( ppp ) potential . @xcite with the help of the density - matrix renormalization - group ( dmrg ) method,@xcite the ground state and elementary excitations for such models can be calculated for large chains with very high accuracy . in this way , the electron - electron interaction is treated without resorting to any approximations .    in a recent study,@xcite we used the hubbard - ohno potential and the tight - binding parameters of ref . [ ] to calculate the binding energy , polarizability and wave function of the singlet exciton , in good agreement with experiment . however , in our previous study we could not reproduce satisfactorily the energy of the triplet ground state . moreover , we did not take polaronic effects into account . here , we shall overcome these shortcomings .    in this work , we perform an extensive dmrg study of the peierls  hubbard - ohno hamiltonian for the @xmath13-electrons on a chain . we start from the tight - binding peierls description of race et al . @xcite but replace their ohno potential  @xcite by the hubbard - ohno potential.@xcite the essential difference between the two parameterizations of the pariser - parr  pople interaction  @xcite lies in the treatment of the coulomb interaction for @xmath13-electrons on the carbon atom . the local hubbard repulsion potential is substantially larger than the corresponding ohno interaction . @xcite    since we are mostly interested in the polaronic effects , we ignore the energetic effects introduced by the ligands  @xmath14 . some preliminary studies show that the energy difference between the singlet and triplet ground states is not influenced by the introduction of a local potential for the carbon atoms which are linked to the side groups . we presume that the dominant influence of the side groups comes from the presence ( or absence ) of strain in the chains . as can be seen from the data for 3bcmu and 4bcmu , different ligands and the resulting strain results in small energetic differences , of about @xmath15 for the nbcmu family . the outline of our work is as follows . in sect . [ sec : exp ] , we summarize the experimental observations on the singlet and triplet in - gap states in pdas . in sect . [ sec : model ] , we define the peierls  hubbard - ohno hamiltonian and the model parameters which provide the basis of our numerical dmrg study . in sect .  [ sec : method ] , we briefly discuss our numerical approach . in sect . [ sec : guessparameters ] , we motivate the parameter regime that we choose for our study . in sect . [ sec : results ] , we present our results . in sect . [ sec : conclusions ] , we summarize and conclude . technical details are deferred to the appendices . we start with an overview of the experimental observations relevant for our study .      first , we collect relevant experimental data on the ground - state properties . the diacetylene monomer building unit is comprised of four carbon atoms . the four outer electrons of each carbon atom are hybridized . three of them form localized bonds . there are @xmath16-bonds between neighboring carbon atoms on the chain . two carbon atoms on the chain share a local @xmath13-bond made by two @xmath17  electrons . the other two carbon atoms share @xmath16-bonds to covalent ligands @xmath18 and @xmath19 , which are several ngstrm long and differ for various members of the pda family . in this work , we focus on poly-(butoxy - carbonyl - nmethylene - urethane ) ( poly - nbcmu ) chains with @xmath20 where the side groups are given by @xmath21 . the fourth carbon electron is delocalized over the carbon backbone in a molecular @xmath13-orbital . due to the peierls effect , in the ground state the @xmath13  electrons dimerize the chain into an alternating sequence of short and long bonds . after dimerization , the four carbon atoms in the unit cell are linked by a triple bond , a single bond , a double bond , and a single bond . the corresponding lewis structure of the ground state is shown in fig .  [ fig : structure ] .    , height=90 ]      for a high - quality single crystal of poly-[1,2-bis-(p - tolylsulfonyloxymethylen)-1-buten-3-inylen ] ( pts ) , the atomic distances at room temperature have been measured  @xcite as @xmath22 , @xmath23 , and @xmath24 [ uncertainties in the last digit in brackets ] for the triple  ( t ) , double  ( d ) , and single  ( s ) bonds , respectively . typical atomic distances for other pda polymer single crystals are  @xcite @xmath25 , @xmath26 , and @xmath27 . the same set of data applies for 3bcmu - pda at low temperatures.@xcite the chain of atoms is not straight ; the single and double bonds alternately form angles of @xmath28 and @xmath29 degrees , with a temperature variation of a few degrees . @xcite in the comparison with our calculations we shall assume that the bond lengths and angles given above for the pda single crystals are representative for chains in their diacetylene monomer matrix .    we shall only deal with planar ( ` blue ' ) pda chains . @xcite the individual polymer chains of 4bcmu are strained in their monomer single crystals but are essentially unstrained in 3bcmu . @xcite strain should be the primary source for differences in the spectra of these two pdas . next , we summarize experimental results on the low - lying electronic excitations . pdas are center - symmetric insulators . their ground state  g is a spin singlet with symmetry @xmath30 under inversion . the charge gap for single - particle excitations , as determined from franz  keldysh oscillations in electro - absorption experiments,@xcite is @xmath31 and @xmath32 in 3bcmu and 4bcmu chains , respectively . the excitation energy of the primary singlet exciton  s ( symmetry @xmath33 ) defines the optical gap , @xmath34 , which amounts to @xmath35 and @xmath36 in 3bcmu and 4bcmu , respectively . therefore , the singlet exciton binding energy , defined by @xmath37 , becomes @xmath38 in 3bcmu , and @xmath39 in 4bcmu , about 24% of the band - gap . the energy levels are sketched in fig . [ fig : energylevels ] .    ) ; s : singlet exciton ( @xmath40 ) ; x@xmath41 , x@xmath42 , x@xmath43 : singlet dark states ( @xmath44 ) ; t : triplet ground state ( @xmath45 ) ; t@xmath46 : optical excitation of the triplet ground state ( @xmath47 ) ; y : dark triplet state ( @xmath48).[fig : energylevels ] ]    the singlet exciton is the energetically lowest state in the spin - singlet sector which can be generated by a single - photon absorption process . in addition , there are further , optically dark states in the gap . the existence of optically dark states  x@xmath41 and x@xmath42 below the optical gap can be inferred from non - radiative decay processes which are monitored via pump - probe spectroscopy , see ref . [ ] for a review . the exciton rapidly populates the states x@xmath49 so that they should have the same spin quantum number .    in principle , the energy of singlet states in the gap with @xmath50-symmetry can be determined via two - photon absorption . two - photon absorption for a single - crystal of the poly - diacetylene paratoluene - sulfonate reveals three gap states with energies @xmath51 , @xmath52 , and @xmath53 . these states should exist in all pdas . note that in the experimental literature ,  @xcite the numbering of the states x@xmath41 and x@xmath42 is reversed . raman scattering reveals vibrational energies which are assigned to the oscillations of the double ( d ) and triple ( t ) bonds . for 3bcmu chains in their monomer matrix they are @xmath54 and @xmath55 . @xcite in accordance with the raman data , the optical spectra of 3bcmu chains show strong exciton replicas at the energies @xmath56 and at @xmath57 , respectively , when the exciton is accompanied by single optical phonons corresponding to the vibrations with frequencies @xmath58 and @xmath59 . electro - absorption measurements  @xcite show that there are more significant single - phonon replicas of the singlet exciton at the energies @xmath60 with @xmath61 and @xmath62 . due to the fact that @xmath63 , the phonon replica at @xmath64 appears in the low - energy flank of the replica at energy @xmath65 . the electro - absorption measurements also permit the identification of multi - phonon replicas , e.g. , at replica energies @xmath66 , @xmath67 , @xmath68 , @xmath69 . @xcite      the triplet sector is more difficult to access experimentally because a transition between the spin - zero ground state  g and the lowest spin - one state  t is optically forbidden . optical pump - probe spectroscopy  @xcite reveals that a small fraction of singlet excitons decays into a long - lived state which can be optically excited by the probe pulse . its long life - time indicates that this in - gap state is the spin - triplet ground state  t. the probe pulse generates transitions from  t to t@xmath46 in the triplet sector . the optical gap in the triplet sector is defined as @xmath70 , and amounts to @xmath71 and @xmath72 , respectively . optical pumping above a threshold @xmath73 very efficiently generates states which show a strong optical absorption with energy  @xmath74 . this can be readily understood if the spin - singlet excitations above @xmath75 fission into triplet pairs . @xcite in turn , these triplet pairs can recombine into singlets and decay optically . if we ignore lattice effects ( bi - polaron formation ) , we obtain a reasonable estimate for the energy of the triplet state  t , @xmath76 . for 3bcmu and 4bcmu this estimate gives @xmath77 and @xmath78 , respectively . the binding energy of the triplet ground state is defined by @xmath79 . it amounts to @xmath80 and @xmath81 in 3bcmu and 4bcmu , respectively , more than 60% of the single - particle gap . the energy of the optically excited triplet ground state  t@xmath46 is found to be @xmath82 and @xmath83 above the ground state , about @xmath15 below the threshold @xmath84 for single - particle excitations . the optically dark singlet states x@xmath41 and x@xmath42 decay non - radiatively into a state  y. its weak population and long life - time indicate that it is reached via intersystem crossing so that it ought to be a spin - triplet state which lies energetically above the triplet ground state  t. it should be of symmetry @xmath85 so that it can not be reached via an optical excitation of the triplet ground state . [ fig : energylevels ] shows the experimentally observed level spectrum for pda chains in the singlet and triplet sectors . we summarize the corresponding values for the in - gap states in table  [ tab : bcmu - energies ] , and compare them to our theoretical results for our best parameter set , @xmath4 , see sect . [ sec : results ] . .first and second column : excitation energies in 3bcmu and 4bcmu at low temperatures . all energies are measured in ev relative to the energy of the ground state , @xmath86 . bold number : directly measured ; italic number : estimate . third column : our results , see sect . [ sec : results ] ; the numbers in square brackets give the excitation energy for the rigid - lattice transition from  g ( @xmath84 , @xmath87 , @xmath88 , @xmath89 ) and from  t ( @xmath74).[tab : bcmu - energies ] [ cols= " < , < , < , < " , ]     the two tables  [ tab:3bcmu - optimization ] and  [ tab:4bcmu - optimization ] indicate that the best set for both 3bcmu and 4bcmu is @xmath90 . then , a look at fig . [ fig : gapsextrapolated - examples ] shows that @xmath91 is the best value for which we have data available . therefore , we shall use the set @xmath92 as our optimal parameter set . the trend shows that @xmath1 might even be a bit smaller , @xmath93 , and @xmath0 a bit larger , @xmath94 , for the optimal case . we plan to perform a more systematic estimate for the optimal parameter set in the near future .    for our optimal parameter set @xmath92 , the corresponding theoretical values for the excitation energies are @xmath95 , @xmath96 , @xmath97 , and @xmath98 . they are also given in table  [ tab : bcmu - energies ] .      in the triplet sector , we find a series of optically dark states  y@xmath99 ( @xmath100 ) at energies just above the triplet ground state . as indicated in fig . [ fig : energylevels ] , the optically dark triplet states  y@xmath99 open decay channels for the states x@xmath41 and x@xmath42 . the structure of the optically dark in - gap triplet states is also interesting from a theoretical point of view . in a band picture , the triplet ground state with spin component @xmath101 consists of a hole in the @xmath102 valence band and an electron in the @xmath103 conduction band , both at momentum @xmath104 . close in energy are the corresponding excitations at finite but small momentum  @xmath105 . the dispersion of these excitations is quadratic as a function of  @xmath105 . , as a function of inverse system size @xmath106 for @xmath107 and @xmath108 . the crossings of the ordinate are equidistant , @xmath109 , @xmath110 , @xmath111 , and @xmath112 , indicating a linear dispersion relation for the excitations . [ fig : conformaltower ] ]    in the interacting case , these excitations may form a spin - flip density wave with momentum  @xmath113 whose dispersion relation at low - energy is given by @xmath114 where @xmath115 is the ( sound ) velocity . for finite chains , we have quantized quasi - momenta , @xmath116 $ ] ( @xmath117 ) , where @xmath8 is the length of the unit cell . therefore , the levels y@xmath99 should obey @xmath118 we confirm this hypothesis in fig . [ fig : conformaltower ] where we show @xmath119 as a function of @xmath106 . as seen from the figure , the energy differences scale to a finite value linearly in @xmath106 . moreover , the extrapolated values are equidistant from which we can read off @xmath120 . using @xmath121 we thus estimate the velocity for the spin - flip density excitations as @xmath122 . for localized spin models , we have @xmath123 so that , approximating @xmath124 , the effective magnetic interaction in our system is of the order of one mev , @xmath125 .      in fig .  [ fig : darksinglets]a we show the ( relaxed ) excitation energies for the four energetically lowest singlet states as a function of @xmath106 for @xmath92 . for our optimal parameter set , there are at least three optically dark singlet states below the singlet exciton . we extrapolate @xmath126 which is about @xmath127 higher than estimated from experiment , see table  [ tab : bcmu - energies ] . this indicates that the local correlations could still be larger than @xmath128 . note , however , that the absolute positions of x@xmath41 and x@xmath42 have not been determined experimentally for nbcmu but it is known from pump - probe spectroscopy  @xcite that there are ( at least ) two optically dark singlet states below the singlet exciton . the next two dark in - gap singlets are almost degenerate in energy , @xmath129 and @xmath130 . they lie below the singlet exciton , as seen in experiment , but about @xmath131 higher than estimated in table  [ tab : bcmu - energies ] . [ fig : darksinglets]b shows the dipole overlap  @xcite for the four singlet states x@xmath41 , x@xmath132 , x@xmath133 , and  s at chain length  @xmath134 . only the exciton  s has a finite dipole overlap . note the energy shift of all in - gap states due to the intrinsic peierls effect . the peierls shift amounts to several tenths of an ev .     for @xmath107 for @xmath108 . after the extrapolation , the states x@xmath132 and x@xmath133 are essentially degenerate , @xmath135 , @xmath136 , @xmath137 , and @xmath138 . the horizontal dashed line gives the value for the charge gap in 3bcmu . ( b )  overlap intensity of the dipole operator for spin - singlet in - gap states after the extrinsic and the intrinsic relaxation for @xmath134 sites . [ fig : darksinglets ] ]    as in the triplet sector , we employ the band picture at fixed particle number to view the elementary excitations of the ground state as a hole in the valence band and an electron in the conduction band . in a wannier picture , the electron - electron interaction forms bound states from these pairs , such as x@xmath41 , the lowest - lying @xmath30 singlet , and  s , the singlet exciton with @xmath139 symmetry . apart from these bound states , there should be a continuum of scattering states . the degeneracy of the states x@xmath132 and x@xmath133 indicates that they are near the threshold to the x - continuum . indeed , two - photon absorption above @xmath73 excites states that can fission into two triplets.@xcite unfortunately , it takes a significant amount of cpu time in our dmrg approach to target more than four in - gap singlet states simultaneously so that a more detailed investigation of the in - gap spectrum remains an open problem .      finally , we show results for the lattice constants in fig .  [ fig : latticeconstants ] . for our optimal set of parameters , the extrapolated values are @xmath140 , @xmath141 , and @xmath142 , for the single , double , and triple bond , respectively . when compared to the experimental values for pda single crystals given in sect . [ subsec : structure ] , the values for the single and double bonds are rather good but the triple bond is too large . this can also be seen from the size of the unit cell . given the lewis structure of fig . [ fig : structure ] , the unit cell in chain direction has the length  @xmath8 with @xmath143 . using @xmath28 , this results in @xmath144 which is slightly larger than the experimental values for 3bcmu chains , @xmath145.@xcite    , @xmath146 and @xmath147 for the single , double , and triple bonds in the chain center as function of inverse system size @xmath106 for @xmath148 for @xmath149 , and @xmath150 . the extrapolated values for the best parameter set ( @xmath91 ) are @xmath140 , @xmath141 , and @xmath142 . [ fig : latticeconstants ] ]    the comparison shows that the spring constants  @xmath151 and  @xmath152 are too large and/or the electron - phonon coupling constant @xmath1 is too small . consequently , the spring constants @xmath151 , @xmath152 should be included as parameters in the optimization procedure . in our study we use quite sizable values for the coulomb parameters . substantial values for the hubbard interaction , @xmath153 , were used by chandross , mazumdar et al .  in their studies of the optical properties of poly - phenylene - vinylene ( ppv ) thin films.@xcite for pda chains immersed in their monomer matrix , we find @xmath128 and confirm their previous result for the ratio between @xmath2 and  @xmath3 , @xmath154.@xcite note that we additionally screen the long - range part of the coulomb interaction by the dielectric constant of the monomer matrix , @xmath155 . a further increase of the hubbard interaction beyond @xmath128 would be problematic , as can be seen from fig .  [ fig : c1]a . for fixed @xmath0 and @xmath1 , the spread of the gaps increases with increasing  @xmath2 . in particular , the calculated singlet - exciton energy would deviate significantly from its experimental value . moreover , the number of dark singlet states below the exciton would become larger than expected from experiment . for the four gaps @xmath84 , @xmath156 , @xmath157 , and @xmath74 ( a )  as a function of  @xmath2 for @xmath158 and @xmath159 , ( b )  as a function of @xmath1 for @xmath158 and @xmath128 , and ( c )  as a function of @xmath0 for @xmath128 and @xmath159.[fig : c1 ] ]    for fixed @xmath158 and @xmath128 , a change from @xmath159 would increase the spread in the gaps , as can be seen from fig . [ fig : c1]b . therefore , we advocate an electron - phonon coupling strength which is some 15% smaller than the value proposed by ehrenfreund et al . , @xmath160.@xcite    in the literature , typical values for the electron transfer integral are @xmath161  @xcite and @xmath158  @xcite . our work indicates that a readjustment of the bare electron transfer integral is not necessary for nbcmu . this can be seen from fig . [ fig : c1]c , where the spread of @xmath162 increases for both smaller and larger  @xmath0 . our study can be extended in several directions . first , the strength of the spring constants  @xmath151 and  @xmath152 must be determined self - consistently . in this way , a better agreement between theory and experiment for the lattice parameters can be obtained . second , the parameter search in the four - dimensional space @xmath163 can be optimized by performing a smooth interpolation between the 81  data sets which we have investigated numerically thus far . moreover , the present analysis does not distinguish between various ligand types ( 3bcmu , 4bcmu ) . the differences are partly due to strain , and a distinction between pdas must therefore be incorporated in the lattice parameters . in addition , different side - groups introduce ( small ) electrostatic potential at the carbon atoms which are linked to the side groups . our numerical analysis shows that the influence of an electrostatic potential of the order of @xmath164 at the side - group sites changes the spectra by @xmath165 so that this contribution can be ignored to first approximation . third , for an optimal parameter set it is interesting to study the properties of the ground and excited states in more detail . for example , the lattice structure of the exciton ( polaron - exciton ) and its polarizability can be calculated  @xcite and be compared with experiment . moreover , further in - gap states should be address , e.g. , the spin-2 ground state @xmath166 , or a second singlet exciton @xmath167 . these tasks are left for a future study . we thank gerhard weiser , michel schott , and benjamin janesko for useful discussions . this work was supported in part by the deutsche forschungsgemeinschaft through grk  790 , and by the hungarian research fund ( otka ) grants nos . k  100908 and k  73455 . acknowledges support from the alexander - von - humboldt foundation and from eth zurich during his time as a visiting professor . following ref .  [ ] , we diagonalize the two - site peierls  hubbard - ohno model in the spin singlet sector to derive the ground - state energy . to this end we consider the two states in site representation @xmath168 in the subspace @xmath169 , the hamilton matrix of the electronic problem has the entries @xmath170 ( @xmath171 ) . a short calculation gives @xmath172 -2t_1(\delta^{\rm e } ) & u/2 -v(\delta^{\rm e } ) \end{array } \label{appeq : matrix } \right ) \ ; , \ ] ] where @xmath173 @xmath174 is the length of the @xmath16-@xmath17 double bond , @xmath175 is the corresponding electron transfer amplitude , @xmath2 is the strength of the electrons local coulomb repulsion , and @xmath176 is their interaction on neighboring sites . note that in ref . [ ] , the unshifted energies were used , i.e. , the previous expressions follow from ours after an energy shift by @xmath177 . the ground - state energy @xmath178 of the two - electron system follows from the diagonalization of the matrix  ( [ appeq : matrix ] ) as @xmath179 ^ 2+[4 t_1(\delta^{\rm e})]^2 } \right ) \;. \label{eq : for - gergo}\ ] ] now that we know the ground - state energy of the electronic hamiltonian explicitly , we do not have to invoke the hellmann  feynman theorem in order to determine the optimal values for the electron transfer amplitudes  @xmath180 . instead , we directly minimize the total ground - state energy @xmath181 with respect to @xmath180 ( @xmath182 ) . we set @xmath183^{3/2 } } \ ; .\ ] ] therefore , the optimization of the electron - lattice problem for two carbon atoms with a double bond leads to the implicit equation @xmath184v'(\delta^{\rm e})-8t_1(\delta^{\rm e } ) } { 2\epsilon_0(\delta^{\rm e})+v(\delta^{\rm e})}\ ; ,   \label{appeq : implicit}\ ] ] which is solved iteratively . assuming the lewis structure of fig . [ fig : structure ] , the position of the atoms in the unit cell are denoted by the two - dimensional vectors @xmath185 , @xmath186 , @xmath187 , and @xmath188 for @xmath189 . their equilibrium positions are denoted as @xmath190 , @xmath191 , @xmath192 , and @xmath193 . the pda structure implies @xmath194 , @xmath195 , and @xmath196 for the singlet , triplet , and doublet bonds , where @xmath197 . we shall treat the atomic motions classically . for convenience , we use periodic boundary conditions , @xmath198 . the kinetic energy of the atoms is given by @xmath199 \ ; .\ ] ] in the spring - constant model , the atoms potential energy is approximated by ( @xmath200 ) @xmath201 \;.\nonumber   \\\end{aligned}\ ] ] to second order in the displacement @xmath202 we can write @xmath203 we define the ( small ) displacements @xmath204 , @xmath205 , @xmath206 , and @xmath207 . then , the lagrange function in the harmonic approximation becomes @xmath208 with @xmath209 \ ; , \nonumber \\ v&= &   \frac{k}{2 } \sum_{l=1}^l \biggl [ \left ( b_l^x - a_l^x\right)^2 + g \left ( c_l^x - b_l^x\right)^2 +   \left ( d_l^x - c_l^x\right)^2 \nonumber   \\ & & + \left [   \cos(\phi ) \left ( a_{l+1}^x - d_l^x\right ) - \sin(\phi ) \left(a_{l+1}^y - d_l^y\right ) \right]^2   \biggr ] \;.   \nonumber\\\end{aligned}\ ] ] the variables @xmath210 and @xmath211 are cyclic and drop out of the problem .      the euler  lagrange equations can be solved using the fourier ansatz @xmath212 with @xmath213 , @xmath214 as the crystal momentum . here , @xmath215 corresponds to @xmath216 . the resulting set of six algebraic equations can be cast into a matrix equation , @xmath217 for the vector @xmath218 . with the abbreviations @xmath219 , @xmath220 , and @xmath221 , the matrix reads @xmath222the equation for the vibrational frequencies @xmath223 as a function of the crystal momentum @xmath105 results from the characteristic equation , @xmath224 . we are interested in the optical phonon modes , @xmath225 . the characteristic equation reduces to @xmath226 for all @xmath227 and @xmath228 . the finite - frequency solutions are denoted as @xmath229 ( @xmath230 ) . we set @xmath231 , and @xmath232 result from the three real roots of @xmath233 in  ( [ thirdorderpolynomial ] ) as @xmath234 . the operator for the kinetic energy  @xmath235 in eq . ( [ eqn : hatt ] ) is readily diagonalized for periodic boundary conditions , @xmath236 . for the ground state of non - interacting electrons , the unit cell consists of four sites , @xmath237 , and the electron transfer amplitudes follow the periodic pattern @xmath238 . thus , we may write @xmath239 we introduce the four operators @xmath240 for electrons with quasi - momentum @xmath241 , @xmath242 via @xmath243 the inverse transformation reads @xmath244 where we used that the sum over the quasi - momenta  @xmath105 generates the orthogonality relation @xmath245 for two lattice indices @xmath246 . in turn , the sum over lattice indices  @xmath247 leads to the orthogonality relation @xmath248 for two quasi - momenta @xmath249 . when we apply the transformation  ( [ trafo ] ) to the kinetic energy  ( [ eqn : tperiodic ] ) , we obtain @xmath250 the remaining task is to band - diagonalize the kinetic energy . to this end we diagonalize the @xmath251-matrix @xmath252 with @xmath253 its eigenvalues in ascending order are @xmath254 , @xmath255 , @xmath256 , @xmath257 with @xmath258 . we find  @xcite @xmath259 ^ 2 & = &   t_{\rm s}^2+t_{\rm d}^2/2+t_{\rm t}^2/2 \nonumber \\ & & \pm   \sqrt{(t_{\rm d}^2-t_{\rm t}^2)^2/4+t_{\rm s}^2 \left[t_{\rm d}^2+t_{\rm t}^2 + 2t_{\rm d}t_{\rm t}\cos(k)\right ] } \nonumber \ ; . \\ \label{appeq : eps}\end{aligned}\]]we write @xmath260 where @xmath261 is a diagonal matrix which contains the eigenvalues of @xmath252 . we define the band - diagonal operators @xmath262 and find @xmath263 \ ; .\end{aligned}\ ] ] the @xmath1- and @xmath264-bands are the valence bands which are filled in the ground state at half band - filling . the @xmath265- and @xmath266-bands are the conduction bands which are empty in the half - filled ground state .    the bare gap at half band - filling obeys @xmath267 with @xmath268 ^ 2 =    t_{\rm s}^2+t_{\rm d}^2/2+t_{\rm t}^2/2 - \sqrt{4t_{\rm s}^2+(t_{\rm t}-t_{\rm d})^2}(t_{\rm d}+t_{\rm t})/2 \ ; .\ ] ] for small deviations , @xmath269 , this simplifies to  @xcite @xmath270                          j .- w . van der horst , p.a . bobbert , m.a.j . michels , g.  brocks , and p.j . kelly , phys . lett .  * 83 * , 4413 ( 1999 ) ; j .- w .  van der horst , p.a . bobbert , p.h.l . de jong , m.a.j . michels , g.  brocks , and p.j . kelly , phys . rev .  b  * 61 * , 15817 ( 2000 ) ; j .- w .  van der horst , p.a . bobbert , and m.a.j . michels , phys . b  * 66 * , 035206 ( 2002 ) .
Abstract:

Answer: we study theoretically poly - diacetylene chains diluted in their monomer matrix . we employ the density - matrix renormalization group method ( dmrg ) on finite chains to calculate the ground state and low - lying excitations of the corresponding peierls hubbard - ohno hamiltonian which is characterized by the electron transfer amplitude  @xmath0 between nearest neighbors , by the electron - phonon coupling constant  @xmath1 , by the hubbard interaction  @xmath2 , and by the long - range interaction  @xmath3 . we treat the lattice relaxation in the adiabatic limit , i.e. , we calculate the polaronic lattice distortions for each excited state . using chains with up to 102  lattice sites , we can safely perform the extrapolation to the thermodynamic limit for the ground - state energy and conformation , the single - particle gap , and the energies of the singlet exciton , the triplet ground state , and the optical excitation of the triplet ground state . the corresponding gaps are known with high precision from experiment . we determine a coherent parameter set @xmath4 from a fit of the experimental gap energies to the theoretical values which we obtain for 81  parameter points in the four dimensional search space @xmath5 . we identify dark in - gap states in the singlet and triplet sectors as seen in experiment . using a fairly stiff spring constant , the length of our unit cell is about one percent larger than its experimental value .