the rare decays @xmath1 are among the most important probes of flavour physics .
they are potentially sensitive to dynamics beyond the standard model ( sm ) and have been intensely studied in the literature @xcite .
measurements have been performed at the @xmath0-meson factories and at the fermilab tevatron @xcite .
excellent future prospects for detailed measurements are provided by the lhc experiments atlas , cms , and lhcb at cern @xcite , and , in the longer run , by super flavour factories based on @xmath7 colliders .
the calculability of @xmath1 decay rates and distributions benefits from the fact that these processes are , to first approximation , semileptonic modes .
correspondingly , the hadronic physics is described by @xmath8 form factors , which multiply a perturbatively calculable amplitude .
this simple picture is not exact because also the nonleptonic weak hamiltonian at scale @xmath9 has @xmath1 matrix elements .
the prominent example is given by hadronic interactions of the form @xmath10 , where the charm quarks annihilate into @xmath11 through a virtual photon .
such charm - loop contributions are more complicated theoretically than the form - factor terms .
even though the charm loops are subdominant numerically in the kinematical regions of interest , they can not be completely neglected . in particular
, the related uncertainty needs to be properly estimated in order to obtain accurate predictions .
we need to distinguish three regions in the dilepton invariant mass @xmath2 , for which the properties of charm loops are markedly different . for @xmath12
the presence of very narrow @xmath13 resonances leads to huge violations of quark hadron duality @xcite and the hadronic backgrounds from @xmath14 , followed by @xmath15 , dominate the short distance rate by two orders of magnitude .
this region in @xmath2 can be removed by experimental cuts .
for @xmath16 the kaon is very energetic and the charm loops can be computed systematically in the heavy - quark limit using qcd factorization for @xmath0 decays into light - like mesons @xcite .
this approach was first employed for @xmath1 in @xcite .
the results have many applications .
a summary with detailed references can be found in @xcite ( see also @xcite for a recent analysis ) .
the high-@xmath2 region , @xmath17 has received comparatively little attention . in this case
the kaon energy is around a @xmath18 or below , and ( soft - collinear ) qcd factorization is less justified , becoming invalid close to the endpoint of the spectrum at @xmath19 . on the other hand , the large value of @xmath2 defines a hard scale for the hadronic contribution to @xmath1 .
consequently an operator product expansion ( ope ) can be constructed , which generates an expansion of the amplitude in powers of @xmath20 ( or @xmath21 ) .
charm loops , and other hadronic contributions , are thus approximated as effective interactions that are local on the soft scales set by @xmath22 and @xmath23 .
this simplifies the computation substantially .
in fact , to leading order in the ope the hadronic contribution reduces to a standard form - factor term .
this picture has been first discussed at lowest order in the ope in @xcite , where it was applied to the endpoint region of @xmath1 and @xmath24 . in @xcite
the ope was considered in some detail , including a discussion of power corrections .    in the present paper
we formulate the ope for the high-@xmath2 region of @xmath1 from the outset .
although our approach is similar in spirit to the analysis of @xcite , the concrete implementation is different .
we will also go beyond the estimates presented in @xcite in several ways .
an important difference is that @xcite combines the ope with heavy - quark effective theory ( hqet ) , whereas we prefer to work with @xmath25-quark fields in full qcd .
the latter formulation has the advantage of a simplified operator basis , which makes the structure of power corrections and their evaluation considerably more transparent .
we also retain the kinematical dependence on @xmath2 in the coefficient functions , rather than expanding it around @xmath26 .
we further discuss the issue of quark - hadron duality , which appears relevant because of the existence of @xmath13 resonance structure in the @xmath2 region of interest .
violations of duality are effects beyond any finite order in the ope . using a resonance model based on a proposal by shifman
, we quantify for the first time the size of duality violations in the high-@xmath2 region of @xmath1 .
further aspects and new results of our analysis will be summarized in sections 8 and 9 .
the main conclusion is that @xmath1 is under very good theoretical control also for @xmath17 .
precise predictions can be obtained in terms of the standard form factors , with essentially negligible effects from the additional hadronic parameters related to power corrections and duality violation .
the paper is organized as follows .
section 2 collects basic expressions for later reference . in section 3 our formulation of the ope for @xmath1 at high @xmath2
is described and the power expansion is constructed explicitly , complete to second order in @xmath27 and with a discussion of weak annihilation as an example of a ( small ) third - order correction . in section 4 we present an estimate of the matrix elements of higher - dimensional operators and quantify their impact on the decay amplitudes for both @xmath28 and @xmath29 transitions .
section 5 discusses the connection between the ope for large @xmath2 and qcd factorization for energetic kaons , which are shown to give consistent results at intermediate @xmath30 . in section 6
we address the subject of duality violation in the context of a toy model analysis .
the estimate is then adapted to the case of @xmath31 in section 7 . in this section
we also address conceptual aspects relevant for the existence of the ope and the notion of quark - hadron duality .
a comparison of our approach with the literature is given in section 8 before we conclude in section 9 .
details on the basis of operators in the ope are described in appendix a and some numerical input is collected in appendix b.
the effective hamiltonian for @xmath32 transitions reads @xcite @xmath33\ ] ] where @xmath34 the operators are given by @xmath35 note that the numbering of @xmath36 is reversed with respect to the convention of @xcite .
our coefficients @xmath37 correspond to @xmath38 in @xcite and we include the factor of @xmath39 in the definition of @xmath40 .
the sign conventions for the electromagnetic and strong couplings correspond to the covariant derivative @xmath41 . with these definitions
the coefficients @xmath42 are negative in the standard model .
we define the kinematic quantities @xmath43 ( where @xmath2 is the dilepton invariant mass squared ) , @xmath44 , and @xmath45 the differential branching fractions for @xmath46 can then be written as @xcite @xmath47 the coefficient @xmath48 contains the wilson coefficient @xmath49 combined with the short - distance parts of the @xmath46 matrix elements of operators @xmath50 . the coefficient @xmath48 multiplies the matrix element of the local operator @xmath51 in the decay amplitude .
the coefficient @xmath52 of the operator @xmath53 is determined by very short distances @xmath54 and is precisely known .
the corresponding formulas for @xmath55 can for instance be found in @xcite .
the amplitudes for the exclusive decays @xmath56 , where @xmath57 , @xmath58 , or a similar meson , are given by the matrix element of the effective hamiltonian in ( [ heff ] ) between the initial @xmath0 meson and the @xmath59 final state .
the dominant contribution comes from the semileptonic operators @xmath40 .
their matrix elements are simple in the sense that all hadronic physics is described by a set of @xmath60 transition form factors .
this is also true for the electromagnetic operator @xmath61 .
the matrix elements of the hadronic operators @xmath62 are more complicated .
they are induced by photon exchange and can be expressed through the matrix element of a correlator between the hadronic part of the effective hamiltonian @xmath63 and the electromagnetic current of the quarks @xmath64 where @xmath65 is the electric charge quantum number of quark flavour @xmath66 and a summation over @xmath66 is understood .
the decay amplitude may thus be written as @xmath67\ ] ] where @xmath68 and @xmath69 are the lepton spinors and @xmath70 for @xmath71 transitions the contribution from @xmath72 is suppressed by the prefactor @xmath73 and can be neglected .    exploiting the presence of the large scale @xmath74 , an operator product expansion ( ope )
can be performed for the non - local term @xmath75 which describes the contribution of 4-quark operators to the @xmath32 amplitude .
such an ope corresponds to integrating out the hard quark loop , leading to a series of local effective interactions for the high-@xmath2 region . to leading order in the large-@xmath2 expansion
this has been presented in @xcite .
a discussion of the ope including higher - order contributions has been given in @xcite .    before going into more detail we discuss the basic structure of the ope for @xmath76 .
the expansion may be written as @xmath77 the operators @xmath78 are composed of quark and gluon fields and have the flavour quantum numbers of @xmath79 .
they are ordered according to their dimension @xmath80 and carry an index @xmath81 labeling different operators with the same dimension .
the @xmath82 are the corresponding wilson coefficients , which can be computed in perturbation theory .
the large scales justifying the expansion are @xmath83 and @xmath2 .
they are counted as quantities of the same order .
the coefficients then scale as @xmath84 in the heavy - quark limit . since the matrix elements @xmath85 scale as @xmath86 , the matrix element of each term in ( [ khco ] ) behaves as @xmath87 .
current conservation implies that all operators are transverse in @xmath66 , @xmath88    it is convenient to work with the @xmath25-quark field in full qcd .
this field could be further expanded within heavy - quark effective theory ( hqet ) , in order to make the @xmath89-dependence fully explicit .
in such an approach many additional operators would arise whose hadronic matrix elements are not readily known .
in contrast , the advantage of using the @xmath25-field in full qcd is that fewer operators appear and that the matrix elements of the leading ones are given by common form factors . in this method
the ope becomes particularly transparent and we will adopt it here .    at leading order in the ope ( @xmath90 ) , illustrated in fig .
[ fig : kh3 ] ,    and in the chiral limit ( @xmath91 ) there are two operators @xmath92 using the equations of motion for the external quarks it can be shown that all possible bilinears @xmath93 and @xmath94 arising from the correlator @xmath76 can be expressed in terms of ( [ o31 ] ) and ( [ o32 ] ) .
consequently , no independent dimension-4 operators of the form @xmath94 can appear in the ope .
the complete proof is given in appendix a. as an example , the operator @xmath95 satisfies the equations - of - motion identity ( for @xmath91 ) @xmath96 for any @xmath97 matrix element with momentum transfer @xmath66 this is equivalent to @xmath98 because of current conservation only the transverse part @xmath99 of such an operator @xmath100 can appear in the ope . from ( [ sdbq ] ) we see that this part can be reduced to a linear combination of ( [ o31 ] ) and ( [ o32 ] ) .
if we keep @xmath101 , two additional operators have to be considered @xmath102 since @xmath103 is small , and numerically similar to @xmath104 , we may formally count these as operators of dimension 4 . because they are absent at order @xmath105
, their impact will be suppressed to the level of @xmath106 , which is negligible .
note that these operators do in any case not introduce new hadronic form factors .    at @xmath107 ( fig .
[ fig : kh5 ] )    we encounter operators with a factor of the gluon field strength @xmath108 , which have the form @xmath109 where the @xmath110 denote dirac and lorentz structures .
we will treat the ope explicitly to the level of @xmath107 , that is including power corrections up to second order in @xmath104 .
although we will not give a full treatment of dimension-6 corrections , we consider as an example the effect of weak annihilation ( fig .
[ fig : kh6 ] ) .
this contribution is characterized by the annihilation of the two valence quarks in the @xmath111 meson in the @xmath112 transition through the weak hamiltonian .
it is described by 4-quark operators , which read schematically @xmath113 with lorentz and dirac structures indicated by @xmath114 , and the light quark field @xmath115 , @xmath80 in the case of non - strange @xmath111 mesons .
weak annihilation provides a mechanism to break isospin symmetry , directly at the level of the transition operator . in @xmath76 weak annihilation , in addition to being a third order power correction , comes only from qcd penguin operators , which have small coefficients . the contribution to isospin breaking from this source
will therefore be strongly suppressed .      in this section
we give explicitly the first few terms in the ope to leading order in renormalization - group improved perturbation theory , that is neglecting relative corrections of @xmath117 .
this order for @xmath118 is relevant in the next - to - leading logarithmic approximation to the @xmath119 amplitude .
we may then write @xmath120 the lower indices of the terms on the r.h.s .
denote the dimension @xmath80 of the corresponding local operators , which come with a coefficient of order @xmath121 .
the first term reads @xmath122\end{aligned}\ ] ] the coefficient in ( [ kh3 ] ) requires a uv renormalization , which has to be consistent with the definition of @xmath123 .
the expression given here corresponds to the ndr scheme used in @xcite .
the function @xmath124 is ( @xmath125 , @xmath126 , @xmath127 ) @xmath128 next we have @xmath129\ , \bar s\gamma_\lambda(1-\gamma_5)gg_{\alpha\rho}b\ , \frac{c_1 q_c}{q^2 } f(x ) \nonumber \\ & & -\frac{q_\lambda}{m_b } \bar s gg_{\alpha\beta } ( g^{\alpha\lambda}\sigma^{\beta\mu } - g^{\alpha\mu}\sigma^{\beta\lambda } ) ( 1+\gamma_5)b\ , \frac{4 c_8 q_b}{q^2}\end{aligned}\ ] ] here @xmath130 with @xmath131 . the charm - loop contribution in ( [ kh5 ] ) , proportional to @xmath132 , can be inferred from @xcite .
note that here we use the convention @xmath133 .
in writing ( [ kh5 ] ) we have neglected terms with the small qcd penguin coefficients @xmath134 .    finally , weak - annihilation diagrams give the dimension-6 term @xmath135 \ , ( \delta_{ij}\delta_{kl } c_4 + \delta_{il}\delta_{kj }
c_3 ) \nonumber\\ & + & \frac{16\pi^2i}{q^4}q_\lambda\sum_{r = u , d }
\bigg [ q_r ( \bar r_i(1-\gamma_5)b_j\ , \bar s_k\sigma^{\mu\lambda}(1+\gamma_5)r_l + \bar r_i\sigma^{\mu\lambda}(1-\gamma_5)b_j\ , \bar s_k(1+\gamma_5)r_l )   \nonumber\\ & & -\frac{1}{3 } ( \bar r_i(1-\gamma_5)b_j\ , \bar s_k\sigma^{\mu\lambda}(1+\gamma_5)r_l -\bar r_i\sigma^{\mu\lambda}(1-\gamma_5)b_j\ , \bar s_k(1+\gamma_5)r_l)\bigg ] \nonumber\\ & & ( \delta_{ij}\delta_{kl } c_6 + \delta_{il}\delta_{kj } c_5)\end{aligned}\ ] ] the terms in ( [ kh6 ] ) only arise from qcd penguin operators , which have small coefficients .
we remark that all operators in ( [ kh3 ] ) , ( [ kh5 ] ) and ( [ kh6 ] ) vanish identically when contracted with @xmath136 , as required by gauge invariance .
the non - factorizable @xmath117 corrections to the charm loop arise from diagrams like the one shown on the r.h.s . of fig .
[ fig : kh3 ] .
the @xmath2-dependence has been recently calculated in analytic form as a taylor expansion in the small parameter @xmath137 @xcite .
analytic results for @xmath138 had been presented in @xcite . in the kinematical range relevant to our considerations
, it has been shown that the convergence of the series is very good .
we therefore use the mathematica input files provided by the authors of @xcite in the online preprint publication for a numerical estimate .
we find that the non - factorizable @xmath117 corrections to the charm loop lead to a 10 - 15% reduction of the real part of @xmath139 and contribute a negative imaginary part of again 10 - 15% relative to the short - distance contribution from @xmath123 ( the precise value is scheme - dependent ) .
this is in agreement with the effect found for the _ inclusive _ @xmath140 rate in the high-@xmath2 region , as discussed in @xcite , and is similar to the effect observed for the low-@xmath2 region in the exclusive decay modes , see table  5 in @xcite .
it is to be stressed that these corrections almost compensate the factorizable charm - loop contribution ( diagram on the l.h.s . in fig .
[ fig : kh3 ] ) .
the reason why the @xmath117 corrections are not suppressed stems from the different colour structure of the diagrams .
whereas the factorizable charm loop comes with a colour - suppressed combination of wilson coefficients , the additional gluon exchange allows the @xmath13-pair to be in a colour - octet state with no such suppression . at
even higher orders in perturbation theory , @xmath141 with @xmath142 , on the other hand , the numerical effect on @xmath139 should really be small , as no new additionally enhanced colour structures will arise .
the computation of the amplitude from the ope requires the evaluation of the matrix elements of the local operators .
we estimate in particular the matrix elements of the higher - dimensional contributions .
this will allow us to quantify power corrections to the @xmath1 amplitude at high @xmath2 .
the cases of @xmath28 and @xmath29 transitions will be considered in turn .
the matrix element of the leading dimension-3 operator is given in terms of the familiar form factors @xmath143 , defined by ( @xmath144 ) @xmath145\ , q^\mu\ ] ]    at the level of the dimension-5 correction in ( [ kh5 ] ) one encounters operators of the form @xmath146 .
their matrix elements introduce , in general , new nonperturbative form factors . using lorentz invariance and the antisymmetry of @xmath147 and @xmath148
one can show that @xmath149 in order to estimate the remaining term we assume @xmath150 for the kaon energy @xmath22 . in this limit
the matrix element can be computed in qcd factorization . to leading order
we then find @xmath151\ ] ] where @xmath152 , @xmath153 is the number of colours , and @xmath154 is the first inverse moment of the @xmath0-meson light - cone distribution amplitude .
this matrix element scales as @xmath155 relative to ( [ bkffpm ] ) in the heavy quark limit , and as @xmath156 for large @xmath153 .    in a similar way we can estimate the weak annihilation term @xmath157\ ] ]
this contribution is power - suppressed as @xmath158 relative to ( [ bkffpm ] ) . the suppression by the small wilson coefficients @xmath159 is partly compensated by a large numerical factor of @xmath160 .
to relative order @xmath158 there is no contribution from the term with @xmath161 and @xmath162 .
note that the result in ( [ bkkh6 ] ) also corresponds to the matrix element obtained when naively factorizing the four - quark operators .
normalized to the amplitude coefficient @xmath163 , the power corrections from @xmath164 and @xmath165 read @xmath166 @xmath167 where @xmath115 , @xmath80 refers to the spectator quark in the @xmath0 meson .
numerically , we find @xmath168 at @xmath169 for central values of the parameters .
this number comes with a substantial uncertainty , in particular from @xmath170 and @xmath171 .
nevertheless , the correction to @xmath172 is very small , most likely below @xmath173 in magnitude .
the correction in ( [ dela9h5 ] ) diminishes further for larger @xmath2 , reaching @xmath174 at the endpoint . towards the endpoint
the kaon becomes soft and the result in ( [ dela9h5 ] ) , based on @xmath175 , can only be viewed as a rough model calculation .
the conclusion that @xmath176 remains negligibly small should however still hold .
the weak annihilation correction is @xmath177 at @xmath169 for @xmath115 and therefore entirely negligible , a consequence also of the small wilson coefficients .      in the case of the decay into a vector meson
the relevant form factors are defined as ( @xmath178 )        a_1(q^2 ) + & - a_2(q^2 )    it is convenient to treat the decay into longitudinally and transversely polarized vector mesons separately . omitting terms proportional to @xmath179 @xmath180\end{aligned}\ ] ]
@xmath181    in the large energy limit @xmath182 , which we may use in the normalization of the power corrections , ( [ bkparff ] ) and ( [ bkperpff ] ) simplify to @xcite @xmath183    the case of a longitudinally polarized @xmath58 is very similar to the case of a pseudoscalar @xmath184 , discussed in section [ subsec : btok ] , and we find @xmath185 where @xmath186 is the decay constant of @xmath187 .    for a @xmath58 with transverse polarization
we obtain @xmath188 where @xmath189 is the decay constant of @xmath190 .
numerically we have @xmath191 and @xmath192 at @xmath193 for our standard set of parameters .
these corrections are of similar size as for the pseudoscalar kaon and they are likewise negligible .
finally , we quote the corrections from weak annihilation    a_9,h6a(k^ * _ ) & = -(c_4 + ) + a_9,h6a(k^ * _ ) & = -(c_6 + ) ( q_r- )    where @xmath115 , @xmath80 refers to the spectator quark in the b meson .
the corrections are tiny , at @xmath193 we have @xmath194 for @xmath115 , and @xmath195 , @xmath196 for @xmath197 .
this is again a consequence of the @xmath198 suppression at large values of @xmath2 , see also @xcite .
the ope for the correlator ( [ khdef ] ) , applied to @xmath1 , is valid as long as the energy of the kaon in the @xmath0-meson rest frame @xmath199 is small compared to @xmath200 .
this condition is certainly fulfilled in the vicinity of the endpoint , but even for @xmath2 as low as @xmath201 , just above the narrow - resonance region , @xmath202 is still fairly small in comparison to the hard scale . on the other hand ,
such a value of @xmath22 is already larger than the qcd scale @xmath23 and one could consider using the factorization methods applicable to the case of energetic kaons . for @xmath2 at @xmath201 , or somewhat above , we have a transition region for the applicability of factorization methods for large @xmath22 and the ope for small @xmath22 .
it is this transition we wish to explore in the present section .
qcd factorization ( qcdf ) at large kaon recoil requires @xmath175 for arbitrary @xmath2 .
the ope method requires @xmath203 , @xmath23 .
both scenarios are consistent with the case @xmath204 which can be realized , at least approximately , as we have seen above . in both , the qcdf and the ope scenario , the leading contributions to the @xmath1 amplitudes are given by short - distance quantities multiplying the standard @xmath8 form factors . in qcdf ,
the first corrections to these form - factor terms come from hard - spectator interactions , which have been computed in @xcite .
normalized to @xmath139 ( called @xmath205 in @xcite ) these corrections may be written as @xmath206 , where the indices refer to the notation of @xcite .
for definiteness we will consider the case of a pseudoscalar kaon in the following .
neglecting the small penguin coefficients @xmath134 and adapting the expressions in @xcite to our notation , the second contribution reads ( @xmath207 ) @xmath208\end{aligned}\ ] ] in the ope limit , defined by treating @xmath209 and taking ( [ q2eklambda ] ) , we have @xmath210 .
this term is therefore subleading with respect to ( [ dela9h5 ] ) , consistent with the absence of a @xmath211 term in @xmath176 .
the remaining term is given by @xmath212 where ( @xmath213 ) @xmath214 @xmath215 @xmath216 @xmath217 @xmath218 the formulas from @xcite have been slightly rewritten here to make them more convenient for application in the high-@xmath2 region where @xmath219 . to make contact with the ope result we eliminate @xmath2 from ( [ tpardef ] ) using @xmath220 and expand in @xmath221 .
the terms proportional to @xmath222 and @xmath223 cancel and we find @xmath224 with the function @xmath225 from ( [ fxdef ] ) . in the limit @xmath226 the qcdf formula ( [ da9p ] )
then reduces to the ope result ( [ dela9h5 ] ) .
we note that the hard - spectator correction ( [ da9p ] ) , which is a leading - power contribution in qcdf , becomes a power - correction @xmath227 in the ope regime .
this behaviour is related to the form factor @xmath228 , which scales as @xmath229 at small @xmath2 and as @xmath230 at large @xmath2 .
we may use the preceding comparison to check the validity of the ope at the relatively low values of @xmath2 around @xmath201 . assuming that the large - energy limit for the kaon is a reasonable approximation , the qcdf expression ( [ da9p ] ) and the ope result ( [ dela9h5 ] ) differ only by the replacement @xmath231 where the numerical values are obtained at @xmath169 for our standard set of parameters , neglecting the kaon mass for consistency .
the real part of the ope approximation is @xmath232 larger in magnitude than the qcdf result , for the imaginary part the discrepancy is about @xmath233 .
such differences are to be expected since the expansion leading from qcdf to the ope is governed by @xmath234 . within this accuracy
, the ope formula still gives a very good estimate of the more complete qcdf result at @xmath169 . for larger values of @xmath2
the numerical difference between the two sides in ( [ da9ph5 ] ) becomes smaller and towards the endpoint the ope is the more appropriate description .
the above exercise suggests that even at @xmath2 as low as @xmath201 the ope is a valid method to obtain the second - order correction ( [ dela9h5 ] ) .
the difference with the qcdf estimate is immaterial in view of the very small overall size of the effect .
the transition from large to small recoil energy has already been considered in @xcite for the @xmath28 form factor ratio @xmath235 .
also in this and in analogous cases , hard - spectator corrections , which are leading - power effects at large recoil , become power suppressed for soft kaons .
similar observations hold for weak annihilation @xcite .
the ope defines a systematic framework to compute the correlator ( [ khdef ] ) at high @xmath2 in qcd . in the minkowski region @xmath236
, there will be uncertainties in the ope - based predictions that go beyond those due to neglected orders in @xmath116 or in @xmath104 .
such effects are refered to as violations of quark - hadron duality .
we investigate their importance first within a toy model for the charm loops in @xmath31 . in section 7
the model is generalized to be closer to the realistic case .
it is illuminating to consider the systematics of duality for the charm - loop contribution to @xmath31 in the simplified context of the toy model introduced in @xcite .
this model assumes the existence of two leptons , @xmath237 with a large mass @xmath238 and @xmath239 with mass @xmath240 , and the effective weak hamiltonian @xmath241\ ] ] all particles have standard strong and electromagnetic interactions .
then @xmath242 gives rise to a loop - induced process @xmath243 via charm- and top - quark penguin diagrams with a gim - like cancellation between them .
the hadronic physics is fully contained in current correlators of the form @xmath244 with @xmath245 .
the decay amplitude can then be written as @xmath246 where @xmath247 is the difference of the charm and top contributions .
we take @xmath248 , and thus @xmath249 comes only from the charm sector .
the correlator @xmath250 fulfills the dispersion relation @xmath251 where the subtraction constant @xmath252 is fixed within the model and can be computed in perturbation theory . to leading order it reads @xmath253 the form of ( [ al12ee ] ) for the amplitude holds to lowest order in @xmath254 and @xmath255 , but to all orders in the strong coupling .
the decay @xmath256 in this model shares important similarities with @xmath257 , but the hadronic dynamics is simplified to the physics of quark - current correlators .
the role of resonances and quark - hadron duality can thus be illustrated in a transparent way .    from ( [ al12ee ] )
we obtain the differential decay rate ( with @xmath258 ) @xmath259 where we defined @xmath260 and @xmath261 to lowest order we have @xmath262 and the partonic expression for @xmath263 in one - loop approximation is @xmath264                             \ln\frac{1+\sqrt{1-x}}{1-\sqrt{1-x } } - i\pi \,;~ & x<1                       \end{array }                                                                     \right.\ ] ] where @xmath131 . for typical values of the parameters ( e.g.  @xmath265gev , @xmath266gev , @xmath267gev ) we have @xmath268 , whereas @xmath269 first rises from 0 at @xmath270 to 4/3 at @xmath271 and then drops again monotonically to the small negative value @xmath272 at @xmath26 .
the coefficient @xmath273 represents the short - distance contribution of the amplitude .
it is real and larger ( parametrically as well as numerically ) than the quark - level charm contribution @xmath274 .
the decay rate ( [ gaml12ee ] ) is proportional to @xmath275 is small compared to @xmath273 , there is a clear hierarchy among the three terms on the r.h.s . of ( [ ccdel ] ) : the first , short - distance term @xmath276 dominates and the next term gives the correction to first - order in @xmath277 , whereas the final term enters at second order .
these features are qualitatively similar in the case of @xmath31 .      in order to investigate the systematics of duality violation in @xmath256
, we find it convenient to consider first a simple model for the quark - current correlator , which has been proposed in @xcite . in this model
the correlator is represented as an infinite sum over resonances , which include finite width effects . in its original form
it applies to massless quarks and we will correspondingly neglect the charm - quark mass in the present section
. a detailed discussion of the model and its use in illustrating duality violation in the @xmath278 ratio and similar quantities has been given in @xcite .
the model has also been used to study duality violation in @xmath279 decays in @xcite .
in shifman s model the correlator @xmath277 in ( [ deldef ] ) reads @xmath280\ ] ] where @xmath281 is the digamma function and @xmath282 @xmath283 is the number of colours , @xmath284 is a scale corresponding to the string tension in qcd and @xmath285 is a ( small ) parameter related to the width - to - mass ratio of the resonances .
the model expression for @xmath277 in ( [ delshm ] ) has the correct analytic behaviour ( a cut for positive @xmath2 but no other singularities on the physical sheet ) and it reproduces the asymptotic result of qcd in the limit of large @xmath2 , @xmath286    using the identity @xmath287_1 + \left[-\pi \cot \pi z + i\pi \right]_2\ ] ] the function in ( [ delshm ] ) can be decomposed into two parts , @xmath288 , corresponding to the two brackets in ( [ psiid ] ) .
@xmath289 is an oscillating function of @xmath2 , exponentially suppressed for large @xmath2 .
it represents the duality violating component of @xmath277 .
the function @xmath290 is a monotonous function , which gives the ope approximation to @xmath277 . the real part of these two functions is shown in fig .
[ fig : rddvope ] .    a plot of the imaginary part of @xmath277 can be found in @xcite .
the duality violating part of @xmath291 can be approximated as @xmath292\approx         -\frac{n\pi}{3 } \exp(-2\pi b r ) \sin(2\pi r)\ ] ] if @xmath293 and @xmath294 .    in the decay rate integrated over the high-@xmath2 part of the spectrum , the duality violating contribution enters proportional to @xmath295 where @xmath296 in ( [ dvrate ] ) we have used the approximation from ( [ deldv ] ) .
typical values of the parameters are @xmath297 the value of @xmath298 corresponds for instance to @xmath299 and @xmath300 . the quantity in ( [ dvrate ] ) is shown in fig .
[ fig : rddvint ] as a function of @xmath301 .
we comment on several important aspects of these results .    *
the duality violating component of @xmath291 in ( [ deldv ] ) exhibits the characteristic oscillating behaviour in @xmath302 with an exponential suppression governed by @xmath303 .
the analogous expression for the duality violating term in @xmath304 , which has a cosine instead of the sine , has been given in @xcite .
( [ dvrate ] ) displays the duality violating contribution from @xmath291 to the partially integrated decay rate .
the integration over @xmath305 extends from a suitably chosen lower limit @xmath301 up to the end of the spectrum .
the parameter @xmath301 should be large enough such that the ope still remains reasonable at the corresponding value of @xmath2 .
using the approximation in ( [ deldv ] ) and expanding in the small quantities @xmath25 and @xmath306 , we find the explicit result written as the last expression in ( [ dvrate ] ) .
this expression is an oscillating function of @xmath301 , with frequency @xmath307 , multiplied by an exponential suppression factor .
the latter is active when the exponent is large , at least @xmath308 .
in addition , the entire term is further suppressed by one power of @xmath309 .
this effect remains even if the exponential suppression is not fully developped .
the power suppression arises because the oscillating contributions average out in the integral , except for a remainder @xmath310 near the lower end of integration @xmath301 . *
the duality violating term ( [ dvrate ] ) is plotted in fig .
[ fig : rddvint ] for a semi - realistic choice of parameters , where @xmath311 and @xmath299 .
the smallest values of @xmath301 shown , @xmath312 , correspond to @xmath313 and @xmath314 , which is already on the low side of the allowed range . from fig .
[ fig : rddvint ] we observe , first , that in the scenario considered here the simple approximation to ( [ dvrate ] ) agrees very well with the full result .
second , the numerical size of the duality violating term is about @xmath315 for the relatively low values of @xmath301 around @xmath316 .
the effect diminishes quickly for larger @xmath301 due to the exponential suppression .
the variation @xmath315 from ( [ dvrate ] ) amounts to @xmath317 when compared with the corresponding ope expression @xmath318    we next turn to a discussion of the @xmath319 term in ( [ ccdel ] ) , where the situation is systematically different from the case of @xmath291 .
the @xmath319 term receives a contribution from the duality violating component given by @xmath320 where in the second step an approximation similar to the one in ( [ deldv ] ) has been used .
there is still an exponential suppression , which makes the entire term negligible for sufficiently large @xmath321 . on the other hand ,
the power suppression with @xmath306 observed in ( [ dvrate ] ) is softened into a behaviour as @xmath322 . for small @xmath25 ( comparable to @xmath306 or smaller )
the violation of duality may become large .
this is in qualitative agreement with the discussion of duality violation for the squared correlator @xmath323 given in @xcite .
the enhanced impact of duality violation is related to the absence of oscillations with alternating sign in the integrand of ( [ delsqdv ] ) .
we conclude that the @xmath319 term in the integrated rate is particularly susceptible to violations of quark - hadron duality , which may lead to substantial ( positive ) deviations from the ope result , unless @xmath2 is large enough for a strong exponential suppression of the effect .
however , the uncertainties in the @xmath319 term may be immaterial if this contribution is only a small correction to the dominant @xmath276 part in ( [ ccdel ] ) .
let us finally illustrate the relative importance of the various contributions to the rate of @xmath256 using the shifman model for the charm loop with the numerical input defined above . in the limit
@xmath324 considered here the function @xmath273 from ( [ cclnmtmc ] ) is modified to @xmath325 the relative contributions to the decay rate ( [ gaml12ee ] ) , ( [ ccdel ] ) then read @xmath326\label{cdnum}\\   \int^1_{0.1 } ds\ , g(s)\ , |\delta_1|^2        & = & 1.25 \qquad [ + 0.10]\label{ddnum}\end{aligned}\ ] ] with @xmath327 .
the central values in ( [ cdnum ] ) and ( [ ddnum ] ) are based on the ope result for @xmath277 , the square brackets indicate the impact of duality violation . the uncertainty in ( [ cdnum ] )
gives the variation due to @xmath328 .
it has the relative size of @xmath317 , which is reduced to @xmath329 in the sum of all contributions .
the shift from duality violation in ( [ ddnum ] ) is positive and essentially negligible in the present example .      in order to obtain a more realistic picture of the @xmath13 spectrum
, we are going to fit the bes data @xcite for the @xmath278ratio in the @xmath13region to a modified version of shifman s model .
the spectra of @xmath13 mesons can be accounted for by linear relations for the squared masses , @xmath330 , @xmath331 , similarly to the case of the light mesons @xcite .
the trajectory of the @xmath332 charmonia , the @xmath333 states @xmath334 , @xmath335 , @xmath336 , @xmath337 , @xmath338 , for instance , follows this pattern . starting form the third resonance ( @xmath339 )
, these states can decay into open charm and have widths of the order of @xmath23 .
the first two are extremely narrow and may be described separately , but their properties are unimportant for duality violation , which is related to the infinite tower of high-@xmath81 resonances .
we therefore choose an ansatz where the sum over resonances begins at @xmath339 rather than @xmath340 .
a finite width is included in analogy to ( [ zdef ] ) and the variable @xmath2 is shifted by a constant into @xmath341 .
this leads to the following expression for the imaginary part of the correlator or , equivalently , the @xmath278 ratio @xmath342 the individual resonances are located at @xmath343 ( @xmath344 ) in the limit @xmath345 .
we observe that a rough description of the bes data @xcite can already be obtained with this formula , where we find @xmath346 from the measured r - ratio below charm threshold , @xmath347 , @xmath348 and @xmath349 .
this yields the result shown in fig .
[ fig : fit0 ] , corresponding to a @xmath350 .
we remark that our values for the fit parameters @xmath351 and @xmath352 are in agreement with the results of @xcite .
there is a second trajectory of @xmath353 charmonia , the @xmath354 states . of these the first two resonances @xmath355 and @xmath356 are known .
the first one is barely above threshold and still rather narrow .
it may be considered separately , similar to @xmath334 and @xmath335 .
note also that @xmath355 is still below our default choice for the lower limit of the high-@xmath2 region , @xmath357 .
the remaining resonances @xmath354 are rather close to the resonances @xmath358 for @xmath142 . for an approximate treatment it appears justified to subsume such a pair of close resonances under a single peak and keep the ansatz given in ( [ step1 ] ) .
the accuracy of this description can be gauged by inspecting fig .
[ fig : fit0 ] . in any case
, the normalization of the second term of @xmath278 in ( [ step1 ] ) is fixed in the large-@xmath2 limit by the free - quark result .
simple fit to bes data for @xmath278 vs. @xmath359 .
right : detailed view.,title="fig:",scaledwidth=45.0% ]   simple fit to bes data for @xmath278 vs. @xmath359 .
right : detailed view.,title="fig:",scaledwidth=45.0% ]    the fit could be refined by including the low - lying , narrow resonances in the fit ansatz .
any finite number of resonances does not change the asymptotic behaviour responsible for duality violation .
also the fit parameters are not much affected by such modifications .
for example , including the @xmath355 reduces the @xmath360 , while the parameters of the continuum ansatz remain almost unchanged , @xmath361 , @xmath347 , @xmath362 and @xmath349
. we will therefore be content with the simple representation given in ( [ step1 ] ) above .
the analytic structure of the matrix element of the operator product in ( [ khdef ] ) can be infered from @xmath363    in order to discuss the properties of the matrix element in ( [ khas ] ) we make the following simplifications .
first , we neglect the small penguin contributions in the weak hamiltonian , that is we take @xmath364 second , in the electromagnetic current we retain only the charm - quark component , @xmath365 , and consider only contributions where the charm fields in @xmath366 are contracted with those in @xmath367 .
this neglects terms where the charm - anticharm pair from @xmath366 annihilates into gluons before connecting to the electromagnetic current .
such contributions are of higher order in @xmath116 and not essential for the problem we want to address .
we consider the matrix element in ( [ khas ] ) as a function of @xmath368 , keeping the kaon energy @xmath369 fixed at a value of order @xmath370 .
this can be achieved by injecting a spurion 4-momentum @xcite @xmath371 into the @xmath366 vertex .
then @xmath372 and for @xmath373 and @xmath369 fixed at their physical values the variable @xmath374 grows with @xmath368 .
the physical kinematics is recovered for @xmath375 . under the simplifying assumptions specified above , the intermediate state @xmath376 in the first sum in ( [ khas ] )
always contains a @xmath13 pair and a strange quark , in general together with other hadronic states , and thus @xmath377 .
the state @xmath376 in the second sum contains a @xmath13 pair and a @xmath25 quark , and @xmath378 .
the matrix element in ( [ khas ] ) is then seen to be an analytic function of @xmath368 in the entire @xmath368-plane , except for two branch cuts at @xmath379 and at @xmath380 .
if we would relax the simplifications above and allowed for intermediate states @xmath376 without @xmath13 pairs , the cuts would extend down to lower values of @xmath381 on the real axis .
the ope of the matrix element in ( [ khas ] ) can be justified for @xmath368 on the imaginary axis , sufficiently far from the origin , that is at @xmath382 , for @xmath383 .
the ope defined in this way in the euclidean can then be analytically continued , term by term , from imaginary @xmath368 onto the positive real axis , corresponding to the minkowskian domain .
terms that are exponentially suppressed in @xmath384 for large positive @xmath385 become oscillating functions of @xmath368 in the minkowskian case , that is for large positive @xmath368 @xcite .
these oscillating terms are invisible at any finite order in the ope and represent the duality violating contribution .      for a quantitative estimate of duality violation
we have to resort to a model of the hadronic correlator in ( [ khas ] ) . to this end
we write the hamiltonian in ( [ hcnopeng ] ) as @xmath386 and assume a factorization of the currents , that is we neglect interactions between @xmath387 and the @xmath388 system .
the coefficient @xmath389 is then treated as a phenomenological parameter . with these simplifications
the correlator in ( [ khas ] ) reduces to @xmath390 where we omitted the longitudinal component @xmath391 .
@xmath392 is the current correlator defined in ( [ pimunu ] ) .
the charm loop in ( [ khfact ] ) contributes to the coefficient @xmath139 in the amplitude of @xmath46 a term @xmath393 in the model of section [ sec : charmbes ] we have @xmath394\ ] ] where @xmath395 and @xmath396 . for the parameters we use the following values @xmath397 we will not employ the model ( [ defd ] ) to describe the entire charm - loop contribution , but only to estimate its duality violating component .
the remainder is more reliably obtained by the ope itself . in order to extract the term that represents duality violation within the model , we decompose @xmath398_1 + \left[-\pi\cot \pi z + i\pi\right]_2\ ] ] this decomposition is useful for @xmath17 , when the second term starts being exponentially suppressed and gives the duality violating contribution . in correspondence with ( [ psiz312 ] ) we have @xmath399 and the duality violating term is @xmath400_2 \approx-\frac{8\pi}{3}\exp(-2\pi br)\left(\sin 2\pi
r - i\cos 2\pi r\right)\ ] ] when integrating the @xmath401 part of the @xmath46 rate over the high-@xmath2 region , from a lower limit @xmath402 to the end of the spectrum , the relative size of the duality violating effect is then given by @xmath403 to first order in the charm - loop contribution .
here we have assumed @xmath389 to be real . using the approximation in ( [ defd2 ] ) and proceeding as in ( [ dvrate ] )
we find @xmath404 this formula gives an excellent approximation of the full result based on ( [ rdv1def ] ) .
we note the ( mild ) exponential suppression and the power suppression by @xmath405 . for @xmath389
we take the value @xmath406 , which is large enough to reproduce the measured @xmath407 branching fraction within the factorization ansatz .
we recall that @xmath172 .
then , for @xmath408 , that is for @xmath409 in the vicinity of @xmath201 , the relative correction @xmath410 is below @xmath411 . in the rate for @xmath46 , @xmath412 is added to @xmath401 , which roughly doubles the result .
the net effect of the uncertainty from ( [ rdv1approx ] ) for the rate integrated over the high-@xmath2 region is therefore only about @xmath413 .    as discussed in section [ sec : dualitytoy ] , the second order effect in @xmath414 is qualitatively different .
it has no cancellations due to oscillating terms , giving a positive correction , and its impact increases with decreasing @xmath25 .
if local duality is at least roughly fulfilled , as is the case for high enough @xmath17 , the duality violation from @xmath415 is still suppressed , being of second order in the small quantity @xmath416 .
the relative size of this component is @xmath417 where the second term uses the approximation in ( [ defd2 ] ) .
the approximate form and the full result in ( [ rdv2 ] ) agree reasonably well . with our set of parameters we find @xmath418 at @xmath419 , about half the size of @xmath410 .
the phenomenological value @xmath406 used above is larger than the perturbative result for the coefficient @xmath389 .
thus it effectively absorbs the ( partly unknown ) effects from factorizable and non - factorizable corrections .
we should emphasize here that the analytic structure of the latter is , in general , more complicated than it is implied by the approximation to @xmath420 in ( [ khfact ] ) .
globally using the larger value of @xmath389 in our numerical estimate of duality violation thus corresponds to the pessimistic scenario where the oscillating terms from the non - factorizable corrections to the charm - loop are added coherently . in reality
, we expect that at least some destructive interference between the various contributions appears , and the amount of duality violation should even be smaller than our estimate .
our model for the charm - loop in ( [ khfact ] ) is similar to the ansatz originally proposed in @xcite and used since then in many phenomenological studies .
however , our motivation for considering this model is essentially different .
in contrast to @xcite it is not our goal to model the hadronic effects on the spectrum point by point in the @xmath2 distribution .
more relevant than the detailed shape of the spectrum is the rate integrated over the entire high-@xmath2 region , which is best described in a model - independent way by the ope , as mentioned above .
we rather employ the model to get an indication of the duality violating effects , which are not captured in an ope calculation .
for this purpose the shifman model for the current correlator is adapted to the charm - quark case @xmath392 , with a choice of parameters consistent with the basic features of the most recent experimental data from bes .
note that the krger - sehgal ansatz for @xmath392 @xcite , consisting of a spectral function with a finite number of resonances and a flat continuum for large @xmath2 , contains no information on duality violation .
the model we use to estimate duality violation is very simple and involves many assumptions .
still we expect it to indicate the systematics and the typical size of the effect , which presumably is closely connected with the @xmath13 resonance structure in the charm - loop contribution .
the high-@xmath2 region of @xmath421 has also been analyzed in the framework of an ope in @xcite and a recent application was presented in @xcite .
our approach differs from the analysis of @xcite in several respects .
we go beyond the work of @xcite by addressing in detail the issue of duality violation and the basis of the ope formalism .
in addition , we extend the ope to include second - order power corrections , which we estimate quantitatively . a list of new results is given in the conclusions .
here we would like to comment further on two conceptual differences between @xcite and our formulation .
this concerns , first , the construction of operators in the ope and , second , the treatment of charm quarks .
the operators in our approach are built from @xmath25-quark fields in full qcd rather than using hqet .
this is convenient because the operator basis is simpler and the ope becomes particularly transparent .
another advantage is that the matrix elements of the leading operators are given by the usual form factors in full qcd . unlike in a hqet framework ,
not all dependence on @xmath89 is made explicit , but this is not essential from a practical point of view and still allows the consistent inclusion of power corrections to a given order in the expansion .
this approach has been used for instance in the ope for inclusive @xmath0 decays applied to the computation of the lifetime difference of @xmath422 mesons @xcite .
these investigations included power corrections @xcite as well as corrections of @xmath117 @xcite .
by contrast , in @xcite a matching onto hqet operators is performed from the start .
this leads to a proliferation of operators , whose matrix elements are not given by the usual form factors .
in fact , to reach a simplification of the resulting hqet expressions , the authors of @xcite partly undo the matching , from hqet back to full qcd expressions .
we prefer to employ the full - qcd formulation throughout , in view of the advantages mentioned above .
we next turn to the second point , the treatment of charm .
the authors of @xcite perform an explicit expansion in @xmath423 , which corresponds to assuming the hierarchy @xmath424 .
this implies that the charm quark is not integrated out at the scale @xmath425 and continues to be an active field below this scale .
operators with charm - quark fields are then present in the ope
.    we will argue that it is conceptually simpler to integrate out charm immediately at the scale @xmath89 , assuming a hierarchy @xmath426 , and that this can be done without any loss in accuracy .
charm - quark effects are then entirely contained in the coefficients of the operators , as is apparent from the formulation of the ope given in section 3 .    to illustrate the point we consider the following example .
let us take the scenario where @xmath424 . in this case
the operator coefficients should be evaluated with @xmath138 .
on the other hand , additional operators involving charm would have to be included in the ope . with explicit @xmath13 fields ,
such operators arise at dimension 6 or higher .
for instance , radiating the virtual photon from the two charm lines in the effective hamiltonian ( [ hcnopeng ] ) and leaving all four quark lines open , gives the dimension-6 operator @xmath427     \ , ( \delta_{ij}\delta_{kl } c_1 + \delta_{il}\delta_{kj } c_2)\ ] ] in analogy to ( [ kh6 ] ) .
this operator has @xmath428 matrix elements contributing to the decay amplitude . assuming next that @xmath429 , the contribution ( [ kh6c ] ) can be simplified by integrating out the charm fields in a further step .
contracting the charm lines into a loop , to which a gluon field is attached , induces below @xmath351 the dimension-6 operator @xmath430 the coefficient of this operator reproduces the logarithmic @xmath351-dependence at this order in the ope .
a similar discussion has been presented in @xcite , illustrating how the logarithmic term @xmath431 in the coefficient of the leading dimension-3 operator ( [ kh3 ] ) is recovered in an effective theory with four - quark operators of the type @xmath10 . whereas the latter effect vanishes for @xmath324
, we note that the coefficient in ( [ kh6d ] ) has a logarithmic divergence in this limit .
this is of no consequence if @xmath429 , where @xmath351 still represents a hard scale .
in fact , for higher - dimensional operators generated by ( [ kh6c ] ) , that is operators with more factors of the gluon field and its derivatives , the coefficients will scale as inverse powers of @xmath351 .
a similar situation exists for inclusive semi - leptonic @xmath432 decays , where four - quark operators with charm , analogous to ( [ kh6c ] ) , also appear at third order in the ope for @xmath424 .
the effect of the corresponding charm loops has been refered to as intrinsic charm in @xcite , where the issue was discussed in great detail .
the above consideration indicates how the nonperturbative @xmath433 matrix element of operators such as ( [ kh6c ] ) can be treated in an expansion in @xmath434 .
however , as demonstrated in @xcite for semi - leptonic @xmath432 decays , the same effects are also described in a framework where charm is integrated out at the @xmath89 scale .
no active charm fields need to be considered in this case and ( [ kh6c ] ) is absent from the ope .
a difference between the ope with or without active charm is that the framework with charm fields and with a strong hierachy @xmath424 assumed , would offer the possibility of efficiently resumming logarithmic terms @xmath435 . since such logarithms are
not very large and in view of the additional power suppression of such terms , such resummations appear not to be necessary in practice .
we also stress that effects such as ( [ kh6c ] ) , irrespective of their detailed treatment , are suppressed at least as @xmath436 .
small differences in the method of their calculation are therefore hardly relevant . we thus conclude that integrating out charm at the scale @xmath89 is entirely justified .
the ope is then constructed in a single step and with a simpler operator basis . for these reasons
the approach appears preferable and we have adopted it here .
the amplitude for @xmath437 , @xmath438 , has the general form given in ( [ abmll ] ) .
it contains the component @xmath439 @xmath440 which receives a nonlocal , hadronic contribution @xmath441 from the matrix element of the nonleptonic weak hamiltonian @xmath366 in addition to the semileptonic form factor terms ( @xmath442 ) .
although the hadronic part @xmath443 is relatively small numerically outside the narrow - resonance region ( @xmath444 of @xmath445 ) , it needs to be reliably computed in order to achieve very accurate predictions .
we have presented a detailed study of the hadronic contribution in the region of large dilepton invariant mass @xmath17 , based on an operator product expansion in inverse powers of the hard scale @xmath200 . working with @xmath25-quark fields in full qcd and factorizing the dependence on @xmath89 , @xmath200 and @xmath351 into the coefficient functions , we obtain the following results :    * to leading order in the ope and to all orders in @xmath116 , @xmath446 is expressed in terms of the standard form factors parametrizing the matrix elements @xmath447 and @xmath448 of dimension-3 operators , up to coefficients calculable in perturbation theory .
to lowest order in @xmath116 only @xmath447 is present . * in the chiral limit ( @xmath91 ) , the first power corrections appear only at second order ( @xmath449 ) and are governed by dimension-5 operators with a gluon field of the form @xmath450 .
the corrections are computed explicitly , using the limit @xmath175 for the hadronic matrix elements , and are shown to be smaller than @xmath173 for @xmath445 . * for @xmath101 the dimension-4 operators @xmath451 and @xmath452 , with right - handed strange quarks , can arise in the ope .
because they are absent at order @xmath105 , their contribution will be suppressed to the negligible level of @xmath106 . besides , no new form factors will be introduced by these operators . * within the ope framework the effect of weak annihilation is a natural ingredient , which we have briefly discussed , mainly for illustration . since in addition to being a third - order power correction it comes with small wilson coefficients , its numerical impact of a few permille is entirely negligible .
* we have clarified the relationship between the ope at high @xmath2 and qcd factorization at low @xmath2 by showing that both descriptions yield consistent results at intermediate values of @xmath30 . *
a relevant topic at high @xmath2 is the issue of quark - hadron duality , which is closely related to the existence of an ope .
we have defined the ope with the help of a spurion momentum that allows for an independent scaling of @xmath2 at fixed @xmath453 and kaon energy .
this allows one to clarify the analytic structure of the matrix elements of the correlator ( [ khconcl ] ) .
we then employed a model based on an infinite series of charm resonances to estimate quantitatively the amount of duality violation , resulting in about @xmath329 for the rate integrated over the high-@xmath2 region .
an important aspect is the different sensitivity to duality violation of the contributions to the rate linear or quadratic in @xmath443 .
the quadratic term is more vulnerable to duality violation , but numerically suppressed as a term of second order in the small charm - loop contribution .
the systematics of duality have further been studied in a toy model for the rare decays , in which the factorization of the charm loop is exact .
the main conclusion is that the high-@xmath2 region of @xmath1 is theoretically under excellent control .
decay rates and distributions are perturbatively calculable up to the nonperturbative effects accounted for by the standard form factors in full qcd .
further nonperturbative corrections are strongly suppressed and negligible within the accuracy of a few percent .
the existence of the ope implies that at high @xmath2 the theory of @xmath1 has an even more solid basis than at low @xmath2 .
an example for applications is the combined analysis of @xmath31 and @xmath454 @xcite . here
the form factor dependence can be essentially eliminated , which leads to precision observables sensitive to new physics effects @xcite . beyond @xmath1 ,
the results of our analysis apply to the high-@xmath2 region of exclusive rare decays with similar final states such as @xmath24 , @xmath455 and , with appropriate modifications , to @xmath456 or @xmath457 .
we show that the ope of @xmath76 in ( [ khdef ] ) can be expressed in terms of @xmath458 at the level of dimension-3 and dimension-4 operators in full qcd , in the chiral limit @xmath91 . for @xmath101 , @xmath459
, right - handed strange quarks can also appear , together with a suppression factor of @xmath103 .
this leads to the two additional operators @xmath460 and @xmath461 in ( [ o41 ] ) and ( [ o42 ] ) , which can be counted as operators of dimension 4 .
the matrix elements of @xmath462 , @xmath463 , @xmath460 and @xmath461 between @xmath111 and @xmath464 are all given by the standard @xmath428 form factors .
together this implies that at leading ( dimension 3 ) and next - to - leading order ( dimension 4 ) in the ope , and to all orders in @xmath116 , standard form factors are the only hadronic matrix elements required .
corrections to these terms arise only at second order ( @xmath465 ) in the ope through operators of dimension 5 . in other words
, there are no first - order power corrections ( @xmath466 ) in the ope of @xmath76 , except for the purely kinematical dependence on @xmath2 and @xmath83 in the coefficient functions .
these functions contain first - order terms such as @xmath467 , which however are calculable and do not introduce unknown hadronic matrix elements .
we demonstrate the completeness of the basis @xmath462 , @xmath463 , @xmath460 and @xmath461 for operators of dimension 3 and 4 by enumerating all possibilities consistent with the relevant symmetries , and using the equations of motion .
we first assume @xmath91 .
operators of dimension 3 built from @xmath468 and @xmath25 have the form @xmath469 where lorentz indices can be contracted with the metric @xmath470 , the @xmath471-tensor and factors of @xmath66 to yield a 4-vector with index @xmath472 .
useful relations are @xmath473 which holds for @xmath474 matrix elements , and @xmath475 or similar identities .
exactly three structures can be formed : @xmath476 this can be seen as follows . for @xmath477 in ( [ opdim3 ] )
the only possibility is the third operator in ( [ opdim3b ] ) .
for @xmath478 we obtain the first operator in ( [ opdim3b ] ) and @xmath479 which is equivalent to the third term in ( [ opdim3b ] ) .
for @xmath480 we can write the second operator in ( [ opdim3b ] ) .
in addition , we could form @xmath481 which leads back to the same structure .
this exhausts the possibilities and proves the completeness of the basis in ( [ opdim3b ] ) at the dimension-3 level .
we next consider operators of dimension 4 .
these have , in addition to the fields in ( [ opdim3 ] ) , one covariant derivative acting on the strange quark .
their general form is @xmath482 because of @xmath483 terms with extra covariant derivatives acting on @xmath25 do not lead to independent dimension-4 operators and may be disregarded .    for @xmath477 in ( [ opdim4 ] ) we can write the operator @xmath484 . using the identity in ( [ sdbq ] )
this operator can be expressed in terms of the basis in ( [ opdim3b ] ) : @xmath485 the other possible structure can be reduced as @xmath486    for @xmath478 the dimension-4 operator has the form @xmath487 we list the possible contractions of ( [ sdgamb ] ) into a 4-vector together with the reduction to the basis in ( [ opdim3b ] ) , neglecting operators of dimension 5 . @xmath488 the r.h.s . of ( [ gammac ] )
is reduced to the basic operators via ( [ sdbop3 ] ) .
the terms on the r.h.s . of ( [ gammad ] )
lead back to the expressions ( [ gammac ] ) and ( [ gammab ] ) .
for @xmath480 the dimension-4 operator becomes @xmath489 we list again the possible contractions and their reduction to the operators in ( [ opdim3b ] ) .
@xmath490 the r.h.s . of ( [ sigmaa ] )
is given by ( [ sdbop3 ] ) .
the r.h.s . of ( [ sigmab ] ) is equivalent to ( [ gammab ] ) and ( [ oneqq ] ) , and ( [ sigmac ] ) reduces to ( [ sigmaa ] ) .
this completes the proof that all operators of dimension 3 and 4 , ( [ opdim3 ] ) and ( [ opdim4 ] ) , can be expressed in terms of ( [ opdim3b ] ) for @xmath91 .
taking current conservation into account leaves us with @xmath462 and @xmath463 in ( [ o31o32 ] ) .
similar arguments hold if @xmath468 is replaced by @xmath491 . since right - handed strange quarks come with a factor @xmath492 , the two further operators @xmath460 and @xmath461 in ( [ o41 ] ) and ( [ o42 ] ) are obtained for @xmath101 .    similarly to the derivations above it
can be shown that dimension-5 operators of the form @xmath493 , contracted to a lorentz vector with index @xmath472 , can always be reduced to linear combinations of ( [ opdim3b ] ) and operators containing a factor of the gluon field strength @xmath147 . the former terms correspond to dimension-3 operators with a purely kinematic power suppression .
genuine operators of dimension 5 therefore have the form @xmath494 , as quoted in ( [ o5n ] ) .
in this appendix we collect input we have used in numerical calculations .
the numbers quoted are our central values .
the form factors have an uncertainty of roughly @xmath495 .
the @xmath29 form factors are parametrized as @xcite @xmath496 @xmath497    the @xmath28 form factor is parametrized as @xcite @xmath498 with the default choice @xmath499 , @xmath500=1.6 and @xmath501 .
further parameters are summarized in table  [ tab : parinput ] .
.[tab : parinput ] input parameters in @xmath18 .
@xmath502 . [ cols="^,^,^,^",options="header " , ]
we thank the galileo galilei institute for theoretical physics for hospitality and the infn for partial support during the initial phase of this project .
g.b . thanks the cern theory division for support and hospitality while this paper was being completed .
we thank alexander khodjamirian and thomas mannel for discussions and critical remarks .
this work was supported in part by the dfg cluster of excellence ` origin and structure of the universe ' and by the dfg graduiertenkolleg gk 1054 .
m.  artuso _ et al .
_ , eur .
j.   c * 57 * ( 2008 ) 309 [ arxiv:0801.1833 [ hep - ph ] ] .
a.  ishikawa _ et al .
_ [ belle collaboration ] , phys .  rev .
* 91 * ( 2003 ) 261601 [ arxiv : hep - ex/0308044 ] .
b.  aubert _ et al . _ [ babar collaboration ] , phys .
d * 73 * ( 2006 ) 092001 [ arxiv : hep - ex/0604007 ] .
b.  aubert _ et al . _
[ babar collaboration ] , phys .  rev .
d * 79 * ( 2009 ) 031102 [ arxiv:0804.4412 [ hep - ex ] ] .
b.  aubert _ et al . _
[ babar collaboration ] , phys .
* 102 * ( 2009 ) 091803 [ arxiv:0807.4119 [ hep - ex ] ] . j.  t.  wei _ et al . _
[ belle collaboration ] , phys .
* 103 * ( 2009 ) 171801 [ arxiv:0904.0770 [ hep - ex ] ] .
e.  pueschel [ cdf collaboration ] , arxiv:1005.5572 [ hep - ex ] ; t.  aaltonen _ et al .
_ [ cdf collaboration ] , arxiv:1101.1028 [ hep - ex ] .    m.  bona _ et al .
_ , arxiv:0709.0451 [ hep - ex ] .
t.  kageyama , aip conf .
* 842 * ( 2006 ) 1064 .
t.  browder _ et al .
_ , jhep * 0802 * ( 2008 ) 110 [ arxiv:0710.3799 [ hep - ph ] ] .
t.  e.  browder , t.  gershon , d.  pirjol , a.  soni and j.  zupan , arxiv:0802.3201 [ hep - ph ] .
m.  beneke , g.  buchalla , m.  neubert and c.  t.  sachrajda , eur .
j.   c * 61 * ( 2009 ) 439 [ arxiv:0902.4446 [ hep - ph ] ] .
m.  beneke , g.  buchalla , m.  neubert and c.  t.  sachrajda , phys .
* 83 * ( 1999 ) 1914 [ arxiv : hep - ph/9905312 ] . m.  beneke , g.  buchalla , m.  neubert and c.  t.  sachrajda , nucl
.  phys .
b * 591 * ( 2000 ) 313 [ arxiv : hep - ph/0006124 ] .
m.  beneke , t.  feldmann and d.  seidel , nucl .
b * 612 * ( 2001 ) 25 [ arxiv : hep - ph/0106067 ] .
a.  khodjamirian , t.  mannel , a.  a.  pivovarov and y.  m.  wang , jhep * 1009 * ( 2010 ) 089 [ arxiv:1006.4945 [ hep - ph ] ] .
g.  buchalla and g.  isidori , nucl .
b * 525 * ( 1998 ) 333 [ arxiv : hep - ph/9801456 ] .
b.  grinstein and d.  pirjol , phys .
d * 70 * ( 2004 ) 114005 [ arxiv : hep - ph/0404250 ] .
g.  buchalla , a.  j.  buras and m.  e.  lautenbacher , rev .
* 68 * ( 1996 ) 1125 [ arxiv : hep - ph/9512380 ] .
a.  j.  buras and m.  mnz , phys .
d * 52 * ( 1995 ) 186 [ arxiv : hep - ph/9501281 ] .
m.  misiak , nucl .
b * 393 * ( 1993 ) 23 [ erratum - ibid .
b * 439 * ( 1995 ) 461 ] .
m.  bartsch , m.  beylich , g.  buchalla and d.  n.  gao , jhep * 0911 * ( 2009 ) 011 [ arxiv:0909.1512 [ hep - ph ] ] ; g.  buchalla , arxiv:1010.2674 [ hep - ph ] . g.  buchalla , g.  hiller and g.  isidori , phys .
d * 63 * ( 2001 ) 014015 [ arxiv : hep - ph/0006136 ] .
g.  buchalla , g.  isidori and s.  j.  rey , nucl .
b * 511 * ( 1998 ) 594 [ arxiv : hep - ph/9705253 ] . c.  greub , v.  pilipp and c.  schpbach , jhep * 0812 * ( 2008 )
040 [ arxiv:0810.4077 [ hep - ph ] ] .
d.  seidel , phys .
d * 70 * ( 2004 ) 094038 [ arxiv : hep - ph/0403185 ] .
m.  beneke , t.  feldmann and d.  seidel , eur .
j.   c * 41 * ( 2005 ) 173 [ arxiv : hep - ph/0412400 ] .
j.  charles _ et al .
_ , phys .
d * 60 * ( 1999 ) 014001 [ arxiv : hep - ph/9812358 ] .
m.  beneke and t.  feldmann , nucl .
b * 592 * ( 2001 ) 3 [ arxiv : hep - ph/0008255 ] .
t.  feldmann and j.  matias , jhep * 0301 * ( 2003 ) 074 [ arxiv : hep - ph/0212158 ] .
b.  blok , m.  a.  shifman and d.  x.  zhang , phys .
d * 57 * ( 1998 ) 2691 [ erratum - ibid .
d * 59 * ( 1999 ) 019901 ] [ arxiv : hep - ph/9709333 ] .
m.  a.  shifman , arxiv : hep - ph/0009131 .
m.  shifman , _ in the proceedings of qcd @ work 2003 : 2nd international workshop on quantum chromodynamics : theory and experiment , conversano , italy , 14 - 18 june 2003 , pp 001_. o.  cata , m.  golterman and s.  peris , phys .  rev .
d * 77 * ( 2008 ) 093006 [ arxiv:0803.0246 [ hep - ph ] ] .
j.  z.  bai _ et al . _
[ bes collaboration ] , phys .
* 88 * ( 2002 )
101802 [ arxiv : hep - ex/0102003 ] . m.  ablikim _ et al .
_ , phys .
* 97 * ( 2006 )
262001 [ arxiv : hep - ex/0612054 ] . j.  h.  kuhn , m.  steinhauser and c.  sturm , nucl .
b * 778 * ( 2007 ) 192 [ arxiv : hep - ph/0702103 ] .
s.  s.  gershtein , a.  k.  likhoded and a.  v.  luchinsky , phys .
d * 74 * ( 2006 ) 016002 [ arxiv : hep - ph/0602048 ] .
b.  chibisov , r.  d.  dikeman , m.  a.  shifman and n.  uraltsev , int .
.  phys .
a * 12 * ( 1997 ) 2075 [ arxiv : hep - ph/9605465 ] . f.  krger and l.  m.  sehgal , phys .
b * 380 * ( 1996 ) 199 [ arxiv : hep - ph/9603237 ] .    c.  bobeth , g.  hiller and d.  van dyk , jhep * 1007 * ( 2010 ) 098 [ arxiv:1006.5013 [ hep - ph ] ] .
m.  beneke , g.  buchalla and i.  dunietz , phys .
d * 54 * ( 1996 ) 4419 [ arxiv : hep - ph/9605259 ] .
m.  beneke , g.  buchalla , c.  greub , a.  lenz and u.  nierste , phys .  lett .
b * 459 * ( 1999 ) 631 [ arxiv : hep - ph/9808385 ] .
i.  i.  bigi , n.  uraltsev and r.  zwicky , eur .
j.   c * 50 * ( 2007 ) 539 [ arxiv : hep - ph/0511158 ] . c.  breidenbach , t.  feldmann , t.  mannel and s.  turczyk , phys .
d * 78 * ( 2008 ) 014022 [ arxiv:0805.0971 [ hep - ph ] ] .
i.  bigi , t.  mannel , s.  turczyk and n.  uraltsev , jhep * 1004 * ( 2010 ) 073 [ arxiv:0911.3322 [ hep - ph ] ] .
t.  hurth and d.  wyler _ in _ j.  l.  .
hewett _ et al .
_ , arxiv : hep - ph/0503261 .
p.  ball and r.  zwicky , phys .
d * 71 * ( 2005 ) 014029 [ arxiv : hep - ph/0412079 ] .
m.  beneke and m.  neubert , nucl .
b * 675 * ( 2003 ) 333 [ arxiv : hep - ph/0308039 ] .
p.  ball , g.  w.  jones and r.  zwicky , phys .
d * 75 * ( 2007 ) 054004 [ arxiv : hep - ph/0612081 ] .