Preprint typeset in JHEP style - HYPER VERSION Light-cone sum rules for B π form factors revisited arxiv:8.796v hep-ph] 3 Apr 8 G. Duplančić Max-Planck-Institut fur Physik Werner-Heisenerg-Institut, Föhringer Ring 6, D-885 München, Germany Rudjer Boskovic Institute, Theoretical Physics Division, HR- Zagre, Croatia E-mail: gorand@thphys.ir.hr A. Khodjamirian Theoretische Physik, Fachereich Physik, Universität Siegen, D-5768 Siegen, Germany E-mail: khodjam@hep.physik.uni-siegen.de Th. Mannel Theoretische Physik, Fachereich Physik, Universität Siegen, D-5768 Siegen, Germany E-mail: mannel@hep.physik.uni-siegen.de B. Melić Rudjer Boskovic Institute, Theoretical Physics Division, HR- Zagre, Croatia E-mail: melic@thphys.ir.hr N. Offen Laoratoire de Physique Théorique CNRS/Univ. Paris-Sud, F-945 Orsay, France E-mail: nils.offen@th.u-psud.fr Astract: We reconsider and update the QCD light-cone sum rules for B π form factors. The gluon radiative corrections to the twist- and twist-3 terms in the correlation functions are calculated. The M S -quark mass is employed, instead of the one-loop pole mass used in the previous analyses. The light-cone sum rule for f Bπ q is fitted to the measured q -distriution in B πlν l, fixing the input parameters with the largest uncertainty: the Gegenauer moments of the pion distriution amplitude. For the B π vector form factor at zero momentum transfer we predict f Bπ =.6.4.3. Comining it with the value of the product V u f Bπ extracted from experiment, we otain V u = 3.5 ±.4 ±. ±. 3. In addition, the scalar and penguin B π form factors fbπ q and fbπ T q are calculated. Keywords: B-decays, QCD, Sum rules. Alexander von Humoldt Fellow
Contents. Introduction. Correlation function 3 3. Gluon radiative corrections 5 4. LCSR for B π form factors 7 5. Numerical results 6. Discussion 6 A. Pion distriution amplitudes 9 B. Formulae for gluon radiative corrections B. Amplitudes for f Bπ LCSR B. Amplitudes for f Bπ f Bπ LCSR 6 B.3 Amplitudes for fbπ T LCSR 9 C. Two-point sum rule for f B 3. Introduction The form factors of heavy-to-light transitions at large energies of the final hadrons are among the most important applications of QCD light-cone sum rules LCSR ]. In this paper we concentrate on the B π transition form factors f Bπ, f Bπ and ft Bπ of the electroweak vector u and penguin d currents, respectively. Previously, these form factors have een calculated from LCSR in, 3, 4, 5, 6, 7, 8, 9, ], gradually improving the accuracy. The main advantage of LCSR is the possiility to perform calculations in full QCD, with a finite -quark mass. In the sum rule approach, the B π matrix element is otained from the correlation function of quark currents, rather than estimated directly from a certain factorization ansatz. This correlation function is conveniently designed, so that, at large spacelike external momenta, the operator-product expansion OPE near the light-cone is applicale. Within OPE, the correlation function is factorized in a series of hard-scattering amplitudes convoluted with the pion light-cone distriution amplitudes DA s of growing twist. To otain the B π form factors from the correlation function, one makes use of the hadronic dispersion relation and quark-hadron duality in the B-meson
channel, following the general strategy of QCD sum rules ]. More details can e found in the reviews on LCSR, e.g., in, 3, 4]. A modification of the method, involving B-meson distriution amplitudes and dispersion relation in the pion channel was recently suggested in 5]; the analogous sum rules for B π form factors in soft-collinear effective theory SCET were derived in 6]. LCSR provide analytic expressions for the form factors, including oth hard-scattering and soft end-point contriutions. Because the method is ased on a calculation in full QCD, comined with a rigorous hadronic dispersion relation, the uncertainties in the resulting LCSR are identifiale and assessale. These uncertainties are caused y the truncation of the light-cone OPE, and y the limited accuracy of the universal input, such as the quark masses and parameters of the pion DA s. In addition, a sort of systematic uncertainty is rought y the quark-hadron duality approximation adopted for the contriution of excited hadronic states in the dispersion relation. Importantly, B π form factors are calculale from LCSR in the region of small momentum transfer q large energy of the pion, not yet directly accessile to lattice QCD. The B πlν l decays, with continuously improving experimental data, provide nowadays the most reliale exclusive V u determination. Along with the lattice QCD results, the form factor f Bπ q otained ] from LCSR is used for the V u extraction. Furthermore, the LCSR form factors f, Bπ q can provide inputs for various factorization approaches to exclusive B decays, such as QCD factorization 7], whereas the penguin form factor fbπ T is necessary for the analysis of the rare B πl l decay. Having in mind the importance of B π form factors for the V u determination and for the phenomenological analysis of various exclusive B decays, we decided to reanalyze and update the LCSR for these form factors. One of our motivations was to recalculate the Oα s gluon radiative correction to the twist-3 part of the correlation function. Only a single calculation of this term exists 9, ], whereas the Oα s corrections to the twist- part have een independently otained in 4] and 5]. In what follows, we derive and present the explicit expressions for all Oα s hard-scattering amplitudes and their imaginary parts for the twist- and twist-3 parts of the correlation function and some of these expressions are new. In the OPE of the correlation function the MS mass m µ is used, a natural choice for a virtual -quark propagating in the hard-scattering amplitudes, calculated at large spacelike momentum scales m. Importantly, in the resulting sum rules we keep using the MS mass. Note that the value of m m is rather accurately determined from the ottomonium sum rules. In previous analyses, the one-loop pole mass of the -quark was employed in LCSR. The main motivation was that the pole mass was used also in the two-point sum rule for the B-meson decay constant f B, needed to extract the form factor from LCSR. In the meantime, the f B sum rule is availale also in MS-scheme 8], and we apply this new version here. Furthermore, we fix the most uncertain input parameters, the effective threshold and simultaneously, the Gegenauer moments of the pion twist- DA, y calculating the B- meson mass and the shape of f Bπ q from LCSR and fitting these quantities to their measured values. In addition, the nonperturative parameters of the twist-3,4 pion DA s entering LCSR are updated, using the results of the recent analysis 9].
The paper is organized as follows. In sect. the correlation function is introduced and the leading-order LO terms of OPE are presented, including the contriutions of the pion twist-,3,4 two-particle DA s and twist-3,4 three-particle DA s. In sect. 3 the calculation of the Oα s twist- and twist-3 parts of the correlation function is discussed. In sect. 4 we present LCSR for all three B π form factors. Sect. 5 contains the discussion of the numerical input and results, as well as the estimation of theoretical uncertainties, and finally, the determination of V u. Sect. 6 is devoted to the concluding discussion. App. A contains the necessary formulae and input for the pion DA s. The ulky expressions for the Oα s hard-scattering amplitudes and their imaginary parts are collected in App. B, and the sum rule for f B is given in App. C.. Correlation function The vacuum-to-pion correlation function used to otain the LCSR for the form factors of B π transitions is defined as: F µ p,q = i d 4 x e iq x π p T { ūxγ µ x,m iγ5 d } = { Fq,p q p µ Fq,p q q µ, Γ µ = γ µ F T q,p q p µ q q µ qp ], Γ µ = iσ µν q ν. for the two different u transition currents, For definiteness, we consider the B d π flavour configuration and, for simplicity we use u instead of d in the penguin current, which does not make difference in the adopted isospin symmetry limit. Working in the chiral limit, we neglect the pion mass p = m π = and the u-, d-quark masses, whereas the ratio µ π = m π /m u m d remains finite. At q m and p q m, that is, far from the -flavour thresholds, the quark propagating in the correlation function is highly virtual and the distances near the lightcone x = dominate. It is possile to prove the light-cone dominance, following the same line of arguments as in 5]. Contracting the -quark fields, one expands the vacuum-to pion matrix element in terms of the pion light-cone DA s of growing twist. The light-cone expansion ] of the -quark propagator is used see also 3]: i αx j β = i g s d 4 k π 4 eik x λ dvg µνa a ij vx δ ij /k m m k /k m m k σ µν ] m k vx µγ ν αβ,. where only the free propagator and the one-gluon term are retained. The latter term gives rise to the three-particle DA s in the OPE. Diagrammatically, the contriutions of twoand three-particle DA s to the correlation function are depicted in Fig.. In terms of perturative QCD, these are LO zeroth order in α s contriutions. The Fock components of the pion with multiplicities larger than three, are neglected, as well as the twists higher 3
πp πp d u d u p q q p q q Figure : Diagrams representing the leading-order terms in the correlation function involving the two-particle left and three-particle right pion DA s shown y ovals. Solid, curly and wave lines represent quarks, gluons, and external currents, respectively. than 4. This truncation is justified y the fact that the twist-4 and three-particle corrections to LCSR otained elow turn out to e very small. In addition we include the Oα s gluon radiative corrections to the dominant twist- and twist-3 parts of the correlation function. The OPE result for the invariant amplitude F is then represented as a sum of LO and NLO parts: Fq,p q = F q,p q α sc F 4π F q,p q,.3 and the same for F and F T. The leading-order LO invariant amplitudes F, F, and F T including twist,3,4 contriutions have een otained earlier in 3, 6, 7, ]. We present them here switching to the new notations 9] of the twist-3,4 DA s: { F q,p q = m f du π m ϕ π u µ π uφ p 3π q up m u ] µ π m q 6m m φ σ m q 3πu φ 4πu up m q up u u } dvψ 4π v q up m m f π dv { Dα m q α α 3 vp ] 4m f 3π vq pφ 3π α i Ψ 4π α i Φ 4π α i Ψ 4π α i Φ 4π α i },.4 4
F q,p q = m f π µ π 6 m { du µ π φ p 3π q u up m ] m q φ σ 3π u q up u m m q up u dvψ 4π v },.5 F T q,p q = m f π m f π dv { du m µ π m ϕ π u q up 3m q up φσ 3π u } m q up m φ 4π u q up Dα m q α α 3 vp ] {Ψ 4π α i vφ 4π α i m v Ψ 4π α i Φ 4π α i },.6 where Dα = dα dα dα 3 δ α α α 3, and the definitions of the twist- ϕ π, twist-3 φ p 3π, φσ 3π, Φ 3π and twist-4 φ 4π, ψ 4π, Φ 4π, Ψ 4π, Φ 4π, Ψ 4π pion DA s and their parameters are presented in App. A. Note that all twist-4 terms are suppressed with respect to leading twist- terms, with an additional power of the denominator /m qup compensated y the normalization parameter δπ Λ QCD of the twist-4 DA s. The calculation of the NLO amplitudes F, F,F T will e discussed in the next section. 3. Gluon radiative corrections In the light-cone OPE of the correlation function. each twist component receives gluon radiative corrections. To otain the desired NLO terms, one has to calculate the Oα s one-loop diagrams shown in Fig., convoluting them with the twist- and two-particle twist-3 DA s, respectively. The diagrams are computed using the standard dimensional regularization and M S scheme. In addition, in our calculation the reduction method from ] is employed. The invariant amplitude F in.3 is otained in a factorized form of the convolutions: { F q,p q = f π du T q,p q,uϕ π u µ π m T p q,p q,uφ p 3π u T σ q,p q,uφ σ 3πu] },3. 5
Figure : Diagrams corresponding to the Oα s gluon radiative corrections to the correlation function. where the hard-scattering amplitudes T, T p,σ result from the calculation of the diagrams in Fig.. The two other NLO amplitudes F and m F T have the same expressions with T T, T p,σ p,σ T, and T T T, T T Tp,σ, respectively. The resulting expressions for all hard-scattering amplitudes are presented in App. B. Note that the LO expressions for the correlation functions in.4-.6 also have a factorized, aleit a much simpler form, with the zeroth-order in α s hard-scattering amplitudes stemming from the free propagator of the virtual -quark. In particular, the twist- component in F is a convolution of T = m /m q up ] with ϕ π u. Let us mention some important features of the Oα s terms of OPE. The currents ūγ µ and m iγ5 d in the correlation function are physical and not renormalizale. Hence, the ultraviolet singularities appearing in T and T are canceled y the renormalization of the heavy quark mass. For T T an additional renormalization of the composite qσ µν operator has to e taken into account. Furthermore, in the twist- term in 3. the convolution integral is convergent due to collinear factorization. As explicitly shown in 4, 5], the infrared-collinear divergences of the Oα s diagrams are asored y the well known oneloop evolution 3] of the twist- pion DA. As a result of factorization, a residual dependence on the factorization scale µ f enters the amplitude T and the twist- DA ϕ π. This scale effectively separates the long- and short near the light-cone distances in the correlation function. In the twist-3 part of F, the complete evolution kernel has to include the mixing of two- and three-particle DA s. To avoid these complications, and following 9], the twist-3 pion DA s in 3. are taken in their asymptotic form: φ p u = and φ σ u = 6u u, whereas the nonasymptotic effects in these DA s are only included in the LO part F. We checked that the infrared divergences appearing in the amplitudes T p and T σ cancel in the sum of the φ p and φ σ contriutions with the one-loop renormalization of the parameter µ π i.e., of the quark condensate density. Finally, in accordance with 9, ], all renormalized hard-scattering amplitudes are well ehaved at the end-points u =,, regardless of the 6
form of the DA s. After completing the calculation of OPE terms with the LO NLO accuracy up to twist-4 twist-3, we turn now to the derivation of the sum rules. 4. LCSR for B π form factors In the LCSR approach the B π matrix elements are related to the correlation function. via hadronic dispersion relation in the channel of the γ 5 d current with the fourmomentum squared p q. Inserting hadronic states etween the currents in. one isolates the ground-state B-meson contriutions in the dispersion relations for all three invariant amplitudes: Fq,p q = m B f Bf Bπ q m... B p q Fq,p q = m B f Bf Bπ q f Bπ q ] m B p... q F T q,p q = m B f BfBπ T q m B m π m B p... 4. q where the ellipses indicate the contriutions of heavier states starting from B π. The three B π form factors entering the residues of the B pole in 4. are defined as: π p ūγ µ B d p q = f Bπ q p µ f Bπ q f Bπ q q µ, 4. π p ūσ µν q ν B ] if d p q = q p µ q µ m T B m π q Bπ q µ, m B m π 4.3 and f B = B d m iγ5 d /m B is the B-meson decay constant. Sustituting the OPE results for F, F and F T in l.h.s. of 4., one approximates the contriutions of the heavier states in r.h.s. with the help of quark-hadron duality, introducing the effective threshold parameter s B. After the Borel transformation in the variale p q M, the sum rules for all three B π form factors are otained. The LCSR for the vector form factor reads: ] f Bπ q = em B /M m B f B F q,m,s B α sc F 4π F q,m,s B, 4.4 where F q,m,s B originates from the OPE result for the LO NLO invariant amplitude F q,p q. The LO part of the LCSR has the following expression: { F q,m,s B = m f π due m q ū ϕ π u um u u µ π φ p 3π m u φ σ 3π u m q dφ σ 3π u ] f3π I3π u 6 u m q du m f π u m m u d φ 4π u u } q 4m q du uψ 4π u dvψ 4π v I 4π u, 4.5 7
where ū = u, u = m q /s B q and the short-hand notations introduced for the integrals over three-particle DA s are: u I 3π u = d dα du u I 4π u = d dα du uα /α uα /α Ψ 4π α i Φ 4π α i ] dv Φ 3π α i, α = α α 3, α 3 = u α /v dv Ψ 4π α i Φ 4π α i v α = α α 3, α 3 = u α /v The NLO term in 4.4 is cast in the form of the dispersion relation: F q,m,s B = π s B dse s/m Im s F q,s. 4.6 = f π π s B m m { dse s/m du Im s T q,s,uϕ π u µ π m Im s T p q,s,uφ p 3π u Im st σ q,s,uφ σ 3πu] }, 4.7 where the ulky expressions for the imaginary parts of the amplitudes T,T p,t σ are presented in App. B. The LCSR following from the dispersion relation for the invariant amplitude F in 4. reads: ] f Bπ q f Bπ q = em B /M m B f F q,m,s B α sc F B 4π F q,m,s B, 4.8 where { F q,m,s B = m f π due m q ū µ π φ p um 3π u dφ σ 3π u m u 6u du u } m qψ 4πu. 4.9 Here the contriutions of twist- and of three-particle DA s vanish altogether. Comining 4.4 and 4.8 one is ale to calculate the scalar B π form factor: q fbπq = f Bπ q m B f q. 4. m π 8
Finally, the LCSR for the penguin form factor otained from the third dispersion relation in 4. has the following expression: ] fbπ T q = m B m π e m B /M m B f F T q,m,s B α sc F B 4π F T q,m,s B, 4. where and F T q,m,s B = m f π m q 4 dφ 4π u du u I4π T u = d dα du u due { m q ū ϕ π u um u m u d φ 4π u m q du I4π T u uα /α v Ψ 4π α i Φ 4π α i ] m µ π 3m q } dφ σ 3π u du dv Ψ 4π α i vφ 4π α i v α = α α 3, α 3 = u α /v, 4.. 4.3 The NLO parts F and m F T in LCSR 4.8 and 4., respectively, are represented in the form similar to 4.7, and the corresponding imaginary parts are collected in App. B. For f B entering LCSR we use the well known two-point sum rule 4] otained from the correlator of two m qiγ 5 currents. The latest analyses of this sum rule can e found in 8, 5]; here we employ the M S version 8]. For consistency with LCSR, the sum rule for f B is taken with Oα s accuracy. For convenience, this expression is written down in App. C. Note that the expressions for LCSR in LO are slightly modified as compared to the ones presented in the previous papers. We prefer not to use the so-called surface terms, which originate from the powers of /m q upn with n > in the correlation functions. Instead, we use a completely equivalent ut more compact form, with derivatives of DA s. The twist- NLO part of LCSR for f Bπ, hence, the expressions for T and ImT in App. B, after transition to the pole scheme the additional expressions necessary for this transition are also presented in App. B coincide with the ones otained in 4]. We have also checked an exact numerical coincidence with the twist- NLO part of the sum rule in 5], written in a different analytical form. The explicit expressions for the amplitudes T p,σ, T p,σ, T, and T T Tp,σ,T and their imaginary parts presented in App. B are new. The Oα s spectral density entering the LCSR for f Bπ is given in ] in a different form, that is, with the u-integration performed, making an analytical comparison of our result with this expression very complicated. The numerical comparison is discussed elow, in sect. 6. Furthermore, in 6] the LCSR for the form factor fbπ was otained, and the imaginary 9
s s q m s B q m s B u u Figure 3: Replacing the integration intervals y the contours in the complex planes of u and s variales in the alternative procedure of the numerical integration of NLO amplitudes. part of T was presented. A comparison with our expression for Im T reveals, however, some differences. Since the imaginary parts of the hard-scattering amplitudes have a very cumersome analytical structure, we carried out a special check of these expressions. Each hardscattering amplitude T,... taken as a function of u,q,p q was numerically compared with its dispersion relation in the variale p q = s, where the expression for Im s T,... was sustituted. Note that one has to perform one sutraction in order to render the dispersion integral convergent. In addition, we applied a new method which completely avoids the use of explicit imaginary parts of hard-scattering amplitudes, allowing one to numerically calculate the NLO parts of LCSR, e.g., F q,m,s B in 4.7, analytically continuing integrals to the complex plane. We make use of the fact that the hard-scattering amplitudes T,T p,t σ are analytical functions of the variale s = p q in the upper half of the complex plane, ecause of iǫ s in Feynman propagators. Consider, as an example the twist- part of F given y the integral over s in the second line of 4.7. Since the integration is performed along the real axis, the operation of taking the imaginary part can e moved outside the integral. To proceed, one has to shift the lower limit of the s-integration to any point at q < s < m. This is legitimate ecause all T s are real at s < m. Then one deforms the path of the s-integration, replacing it y a contour in the upper half of the complex plane, as shown schematically in Fig. 3, so that all poles and cuts are away from the integration region. Only when s is approaching the upper limit s B, one nears the pole at u = m q /s B q while performing the integration over u. Because this pole does not touch the limits u =,, it is possile to avoid it y moving the contour of the u-integration into the upper half of the complex u-plane see Fig. 3. After that, oth numerical integrations ecome completely stale. Note, that in oth s- and u-integrations, we integrate over the semi-circle, ut the contour of the integration can e deformed in an aritrary way in the upper half of the complex plane. The numerical integrations of T over these contours yield an imaginary part which represents the desired answer for F. We have checked that the numerical results otained y this alternative method coincide with
twist Parameter Value at µ = GeV Source a π.5 ±.5 average from 9] a π 4 a π. ±. πγγ form factor 3] a π >4 µ π.74.67.38 GeV GMOR relation; m u,d from 8] 3 f 3π.45 ±.5 GeV -point QCD SR 9] ω 3π.5 ±.7 -point QCD SR9] 4 δπ.8 ±.6 GeV -point QCD SR 9] ǫ π 8. ±. -point QCD SR 9] Tale : Input parameters for the pion DA s. the ones otained y the direct integration over the imaginary parts, therey providing an independent check. 5. Numerical results Let us specify the input parameters entering the LCSR 4.4, 4.8 and 4. for B π form factors and the two-point sum rule C. for f B. The value of the -quark mass is taken from one of the most recent determinations 7]: m m = 4.64 ±.5 GeV, 5. ased on the ottomonium sum rules in the four-loop approximation. Note that 5. has a smaller uncertainty than the average over the non-lattice determinations given in 8]: m m = 4. ±.7 GeV. However, as we shall see elow, the uncertainty of m m does not significantly influence the error udget of the final prediction. Furthermore, in our calculation, the scale-dependence m µ m is taken into account in the one-loop approximation which is sufficient for the Oα s -accuracy of the correlation function. Note that using the M S mass inevitaly introduces some scale-dependence of the lower threshold m in the dispersion integrals in oth LCSR and f B sum rule. However, this does not create a prolem, ecause the imaginary part of the OPE correlation function otained from a fixed-order perturative QCD calculation is not an oservale, ut only serves as an approximation for the hadronic spectral density. The QCD coupling α s µ r is otained from α s m Z =.76 ±. 8], with the NLO evolution to the renormalization scale µ r. In addition to µ m and µ r, one encounters the factorization scale µ f in the correlation function, at which the pion DA s are taken. In what follows, we adopt a single scale µ = µ m = µ r = µ f in oth LCSR and two-point SR for f B. The numerical value of µ will e specified elow. The input parameters of the twist- pion DA include f π = 3.7 MeV 8] and the two first Gegenauer moments a π and aπ 4 normalized at a low scale GeV. For the latter
we adopt the intervals presented in Tale. The range for a π GeV is an average 9] over various recent determinations, including, e.g., a π GeV =.6..9 calculated from the two-point sum rule in 9]. For a π 4 we use, following ], the constraint aπ GeV a π 4 GeV =. ±., otained 3] from the analysis of πγγ form factor. Having in mind, that at large scales the renormalization suppresses all higher Gegenauer moments, we set a π >4 = in our ansatz for ϕ πu specified in App. A. The uncertainties of a π,4 GeV remain large, hence we neglect very small effects of their NLO evolution taken into account in 4]. The normalization parameter µ π GeV of the twist-3 two-particle DA s presented in Tale is otained adopting the non-lattice intervals 8] for the light quark masses: m u GeV = 3. ±. MeV, m d GeV = 6. ±.5 MeV. Correspondingly, the quarkcondensate density given y GMOR relation is: qq GeV = f πµ π GeV = 46 8 9 MeV3, 5. where very small Om u,d corrections are neglected. We prefer to use the aove range, rather than a narrower standard interval qq GeV = 4 ± MeV 3 employed in the previous analyses. In fact, 5. is consistent with qq GeV = 54 ± 8 MeV 3 quoted in the review 3], as well as with the recent determination of the light-quark masses from QCD sum rules with Oα 4 s accuracy 3]: m u GeV =.7 ±.4 MeV, m d GeV = 4.8 ±.5 MeV. The remaining parameters of the twist-3 DA s f 3π, ω 3π and twist-4 DA s δπ, ǫ π presented in Tale are taken from 9], where they are calculated from auxiliary twopoint sum rules. The latter are otained from the vacuum correlation functions containing the local quark-gluon operators that enter the matrix elements A.6, A.7 and A., A.3. The one-loop running for all parameters of DA s is taken into account using the scale-dependence relations presented in App. A. Note that the small value of f 3π effectively suppresses all nonasymptotic and three-particle contriutions of the twist-3 DA s. Furthermore, the overall size of the twist-4 contriutions to LCSR is very small. Hence, although the parameters of the twist-3,4 DA s have large uncertainties, only the accuracy of µ π plays a role in LCSR. Finally, in the sum rule C. for f B the gluon condensate density α s /πgg =..6. GeV4 and the ratio of the quark-gluon and quark-condensate densities m =.8 ±. GeV 3] are used, the accuracy of these parameters playing a minor role. The universal parameters listed aove determine the external input for sum rules. The next step is to specify appropriate intervals for the internal parameters: the scale µ, the Borel parameters M and M and the effective thresholds s B and sb. In doing that, we take all external input parameters at their central values, allowing only a π and aπ 4 to vary within the intervals given in Tale. From previous studies 4, 5, 8, ] it is known that an optimal renormalization scale is µ m B m m Λ where Λ does not scale with the heavy quark mass, We also expect that the use of the recently developed renormalon model 33] for the twist-4 DA s, instead of the conventional twist-4 DA s 34] used here, will not noticealy change the numerical results.
and simultaneously, µ has the order of magnitude of the Borel scales defining the average virtuality in the correlation functions. In practice, M and M are varied within the working windows of the respective sum rules, hence one expects that also µ has to e taken in a certain interval. Calculating the total Borel-transformed correlation function that is, the s B limit of LCSR we demand that the contriution of suleading twist-4 terms remains very small, < 3% of the LO twist- term, therey diminishing the contriutions of the higher twists, that are not taken into account in the OPE. This condition puts a lower ound M Mmin = 5 GeV. In addition, in order to keep the α s -expansion in the Boreltransformed correlation function under control, oth NLO twist- and twist-3 terms are kept 3% of their LO counterparts, yielding a lower limit µ.5 GeV. Hereafter a default value µ = 3 GeV is used. Furthermore, we determine the effective threshold parameter s B in LCSR for each M Mmin. We refrain from using equal threshold parameters in LCSR and two-point sum rule for f B, as it was done earlier, e.g. in 4, 8]. Instead, we control the duality approximation y calculating certain oservales directly from LCSR and fitting them to their measured values. Importantly, we include in the fitting procedure not only s B, ut also the two least restricted external parameters a π and aπ 4, under the condition that oth Gegenauer moments remain within the intervals of their direct determination given in Tale. The first oservale used in this analysis is the B-meson mass. In a similar way, as e.g., in, 8], m B is calculated taking the derivative of LCSR over /M and dividing it y the original sum rule. The B-meson mass extracted from LCSR has to deviate from its experimental value m B = 5.79 GeV y less than %. Secondly, we make use of the recent rather accurate measurement of the q -distriution in B πlν y BABAR collaoration 35]. We remind that LCSR for B π form factors are valid up to momentum transfers q m m Λ, typically at < q < 4 5 GeV. To e on the safe side, we take the maximal allowed q slightly lower than in the previous analyses and calculate the slope f Bπ q /f Bπ from LCSR at < q < GeV. The otained ratio is then fitted to the slope of the form factor inferred from the data. We employ the result of 36], where various parameterizations of the form factor f Bπ q are fitted to the measured q -distriution. Since all fits turn out to e almost equally good, we adopt the simplest BK-parameterization 37]: f BK Bπ q f BK Bπ = q /m B α BK q /m B 5.3 with the slope parameter α BK =.53 ±.6 from 36] close to α BK fitted in 35]. After fixing s B for each accessile M, we demand that heavier hadronic states contriute less than 3% of the ground-state B meson contriution to LCSR. This condition yields an upper limit M < Mmax = GeV. The resulting spread of the threshold parameter and Gegenauer moments when M varies etween Mmin and M max is very small: s B = 36 35.5 GeV, a π GeV =.5.7, aπ 4 GeV =.5.3. The 3
.4. f.8 BΠ q f BΠ.6.4. 4 6 8 q GeV Figure 4: The LCSR prediction for the form factor shape f Bπ q /f Bπ fitted to the BK parameterization of the measured q -distriution. The two almost indistinguishale curves are the fit and the parameterization 5.3 at α BK =.53. quality of the fit is illustrated in Fig. 4 where the two curves: the calculated q -shape of the form factor and the BK-parameterization 5.3 are almost indistinguishale. Thus, in our numerical analysis we trade the q -dependence predicted from LCSR for a smaller uncertainty of the Gegenauer moments and for a etter control over the quark-hadron duality approximation. In the final stage of the numerical analysis we turn to the two-point QCD sum rule for f B presented in App. C and find that at the adopted value of the renormalization scale µ = 3 GeV the interval M = 5.±. GeV satisfies the same criteria as the ones imposed in the numerical analysis of LCSR: the smallness of higher power terms in OPE, and suppression of the heavier hadronic contriutions. The threshold parameter s B = 35.6.9. GeV is fixed y calculating m B from this sum rule. This time the deviation from the experimental value is even less than.5 %. For completeness, we quote the resulting interval f B = 4 5 7 MeV. Note that the Oα s correction taken into account in 8] is not included here. As usual, employing the sum rule for f B in order to extract the form factor from the LCSR for the product f B f Bπ turns out to e extremely useful. One oserves a partial cancellation of the α s -corrections in oth LCSR and two-point sum rule and a etter staility with respect to the variation of scales. To demonstrate some important numerical features of the LCSR prediction, in Fig. 5 left we plot the M -dependence of the form factor f Bπ with all other inputs fixed at their central values. The oserved staility, far eyond the adopted working interval in M, serves as a usual criterion of reliaility in QCD sum rule approach. The µ-dependence plotted in Fig. 5 right is very mild from µ =.5 GeV up to µ = 6 GeV. The numerical size of the gluon radiative corrections in LCSR is illustrated in Fig. 6. The numerical analysis yields the following prediction for the vector B π form factor 4
.75.75.7 f BΠ.65.7 f BΠ.65.6.6 4 6 8 4 M GeV.5 3 3.5 4 4.5 5 5.5 6 Μ GeV Figure 5: Dependence of f Bπ on the Borel parameter left and renormalization scale right..8.6.4. NLO f BΠ q -. -.4 -.6 4 6 8 q GeV Figure 6: Gluon radiative corrections to the twist- dotted line and twist-3 solid line parts of LCSR for f Bπ q, as a function of q. The part proportional to φ σ φ p is shown separately y dashed dash-dotted line. at zero momentum transfer:.4 =.63 f Bπ.9.5.4 M,M ±. µ.3. shape ±., 5.4 µπ m where the central value is calculated at µ = 3. GeV, M = 8. GeV, s B = 35.75 GeV, a π GeV =.6, aπ 4 GeV =.4, M = 5. GeV and s B = 35.6 GeV. The percentages of different contriutions to the central value in 5.4 are presented in Tale. In 5.4 the first second uncertainties are due to the variation of the Borel parameters M and M scale µ within the intervals specified aove. The third uncertainty reflects the error of the experimental slope parameter. In addition, we quote the uncertainties due to limited knowledge of the external input parameters. We have estimated them y simply varying these parameters one y one within their intervals and fixing the central values for all internal input parameters. Interestingly, the largest uncertainty of order of % is due to the error in the determination of light-quark masses transformed into the uncertainty of µ π, the coefficient of the large twist-3 LO contriution. The spread 5
caused y the current uncertainty of m m is much smaller, hence, does not influence the resulting total uncertainty, even if one increases the error of m -determination y a factor of two. Remaining theoretical errors caused y the current uncertainties of α s, twist-3,4 DA s parameters and higher-dimensional condensates are very small, and for revity they are not shown in 5.4. Finally, we add all uncertainties in -quark mass MS pole quadrature and otain the interval: input central set II from ] f Bπ.63.58 tw LO 5.5% 39.7% tw NLO 7.4% 7. % tw3 LO 46.7% 4.5 % tw3 NLO -4.4%.4 % tw4 LO -.% -.9% Tale : The form factor f Bπ at zero momentum transfer calculated from LCSR in two different quark-mass schemes and separate contriutions to the sum rule in %. f Bπ =.6.4.3, 5.5 which is our main numerical result. It can e used to normalize the experimentally measured shape, e.g., the one in 5.3, yielding the form factor f Bπ q in the whole q -range of B πlν l. With this prediction at hand, we are in a position to extract V u. For that we use the interval V u f Bπ = 9. ±.6 shape ±.3 BR 4, 5.6 inferred 36] from the measured q -shape 35] and average ranching fraction of B πlν l 38]. We otain: V u = 3.5 ±.4 th ±. shape ±. BR 3, 5.7 where the first error is due to the estimated uncertainty of f Bπ in 5.5, and the two remaining errors originate from the experimental errors in 5.6. A possile small correlation etween the shape uncertainty of our prediction for the form factor and the experimental shape uncertainty is not taken into account. The remaining two B π form factors can now e predicted without any additional input. In particular, we adopt the same Borel parameter M and effective threshold s B, assuming that they only depend on the quantum numers of the interpolating current for B meson. The scalar form factor fbπ q, otained y comining the LCSR for f Bπ and f Bπ f Bπ, and the penguin form factor ft Bπ q are presented in Fig. 7, in comparison with f Bπ q. The predicted interval for the penguin form factor at zero momentum transfer is : fbπ T =.55 ±.35, 5.8 adopting µ = 3 GeV as the renormalization scale of the penguin current. 6. Discussion In this paper, we returned to the LCSR for the B π form factors. We recalculated the Oα s gluon radiative corrections to the twist- and twist-3 hard-scattering amplitudes and 6
.6 f BΠ,T, q.5.4.3 4 6 8 q GeV Figure 7: The LCSR prediction for form factors f Bπ q solid line, f Bπ q dashed line and f T Bπ q dash-dotted line at < q < GeV and for the central values of all input parameters. presented the first complete set of expressions for these amplitudes and their imaginary parts. For the radiative corrections to the twist- part of the LCSR for f Bπ we reproduced the results of 4, 5]. For the radiative corrections to the twist-3 part we confirmed the cancellation of infrared divergences oserved in 9, ] in the case of asymptotic DA s. Including the nonasymptotic effects in these radiative corrections demands taking into account the mixing etween two- and three-particle DA s. In fact, the parameter f 3π determining the size of nonasymptotic twist-3 corrections is numerically small, hence these corrections are not expected to influence the numerical results. Throughout our calculation and in the final sum rule relations we used the MS-mass of the quark, which is the most suitale mass definition for short-distance hard-scattering amplitudes. Indeed, as follows from our numerical analysis, the Oα s corrections to the sum rules turn out to e comparaly small. To demonstrate that, we returned to the pole-mass scheme in LCSR and used exactly the same input as in ] the preferred set with m pole = 4.8 GeV. We calculated the total form factor and separate contriutions to the sum rule in oth quark-mass schemes and compared them in Tale. Note that the twist- NLO correction is distinctively smaller in the MS-mass scheme. In the twist-3 part of the sum rule the α s -correction is small in oth schemes. In the MS scheme, as seen from Fig. 6, this correction is dominated y the contriution of the DA φ p 3π. In the pole scheme there is a partial cancellation etween the contriutions of the two twist-3 DA s. Our numerical results for the form factors f,,t Bπ q in the pole scheme are very close to the ones otained in ]. It is however difficult to compare separate contriutions, ecause they are not presented in ]. We also cannot confirm the numerical values of the twist-3 NLO corrections to the form factor f Bπ q plotted in the figure presented in the earlier pulication 9]. Further improvements of LCSR are possile ut demand sustantial calculational efforts. For example, otaining radiative corrections to the three-particle twist-3,4 contriutions is technically very difficult. Again, we expect no visile change of the predicted form 7
factors ecause three-particle terms are already very small in LO. A more feasile task is to go eyond twist-4 in OPE and estimate the twist-5,6 effects, related to the four-particle pion DA s, at least in the factorization approximation, where one light quark-antiquark pair is replaced y the quark condensate In LCSR for the pion electromagnetic form factor these estimates have een done in 39]. The numerical analysis of LCSR was improved due to the use of the q -shape measurement in B πlν. A smaller theoretical uncertainty of LCSR predictions can e anticipated with additional data on this shape, as well as with more accurate determinations of - and, especially, u, d-quark masses. In this paper, all calculations have een done in full QCD with a finite -quark mass. At the same time, the whole approach naturally relies on the fact that m is a very large scale as compared with Λ QCD and related nonperturative parameters. Our results demonstrate that the twist-hierarchy as well as the perturative expansion of the correlation function work reasonaly well. An interesting prolem is the investigation of the m limit of LCSR and various aspects of this limiting transition, e.g., the hierarchy of radiative and nonasymptotic corrections. This prolem remaining out of our scope was already discussed in several papers: earlier, in 7] at the LO level, in 5],4] at NLO level and more recently, in 6]. ref.] f Bπ q calculation f Bπ q input V u 3 4] lattice n f = 3-3.78±.5±.5 4] lattice n f = 3-3.55±.5±.5 43] - lattice SCET B ππ 3.54 ±.7 ±.44 44] - lattice 3.7 ±. ±. 45] - lattice LCSR 3.47 ±.9 ±.3, 36] LCSR - 3.5 ±.4 ±. this work LCSR - 3.5 ±.4 ±. ±. Tale 3: Recent V u determinations from B πlν l Finally, in Tale 3 we compare our result for V u with the one of the previous LCSR analysis and with the recent lattice QCD determinations otained at large q and extrapolated to small q with the help of various parameterizations. The oserved mutual agreement ensures confidence in the continuously improving V u determination from exclusive B decays. Acknowledgments We are grateful to Th. Feldmann, M. Jamin and R. Zwicky for useful discussions. This work was supported y the Deutsche Forschungsgemeinschaft Project KH 5/-. The work of G.D and B.M was supported y the Ministry of Science, Education and Sport of the Repulic of Croatia, under contract 98-9893-864. The partial support of A. von Humoldt Foundation under the Program of Institute Partnership is acknowledged. The work of N.O. was supported y FLAVIAnet Contract No. MRTN-CT-6-3548. 8
A. Pion distriution amplitudes For convenience, we specify the set of the pion DA s and their parameters used in this paper. The notations and parameters for twist-3 and 4 DA s are taken from 9], where the earlier studies 34, 46] are updated. The two-particle DA s of the pion enter the following decomposition of the ilocal vacuum-pion matrix element for definiteness, π in the final state: π p ū i ωx d j ξ x x = iδij f π du e iup x iūp x /pγ 5 ] ξω ϕ π u γ 5 ] ξω µ π φ p 3π u 6 σ βτγ 5 ] ξω p β x x τ µ π φ σ 3πu 6 /pγ 5] ξω x x φ 4π u i /x /x γ 5 ] ξω u ψ 4π vdv, A. In the aove, the product of the quark fields is expanded near the light-cone, that is, x i = ξ i x, where ξ i are aritrary numers, and x = ; ū = u. The path-ordered gauge-factor Wilson line is omitted assuming the fixed-point gauge for the gluons. The light-cone expansion includes the twist- DA ϕ π, two twist-3 DA s φ p 3π, φσ 3π and two twist- 4 DA s φ 4π and ψ 4π. The usual definitions of DA s are easily otained, multiplying oth parts of A. y the corresponding cominations of γ matrices and taking Dirac and color traces. The decomposition of the three-particle quark-antiquark-gluon matrix element is: π p ū i ω x g s G a µν x 3d j ξ x x = λa ji Dα i e ipα x α x α 3 x 3 3 if 3π σ λρ γ 5 ξω p µ p λ g νρ p ν p λ g µρ Φ 3π α i { f π γ λ γ 5 ξω p ν g µλ p µ g νλ Ψ 4π α i p λp µ x ν p ν x µ } Φ 4π α i Ψ 4π α i p x if { π ǫ µνδργ λ ξω p ρ g δλ p δ g ρλ Ψ 4π α i p λp δ x ρ p ρ x δ p x Φ4π α i Ψ } ] 4π α i. A. including one twist-3 DA Φ 3π and four twist-4 DA s : Φ 4π, Ψ 4π, Φ 4π and Ψ 4π. Here the convention ǫ 3 = is used, which corresponds to Tr{γ 5 γ µ γ ν γ α γ β } = 4iǫ µναβ. The following expressions for the DA s entering the decompositions A. and A. are used: twist- DA: ϕ π u = 6uū a C 3/ u ū a 4 C 3/ 4 u ū, A.3 where, according to our choice, the first two Gegenauer polynomials are included in the nonasymptotic part, with the coefficients having the following LO scale depen- 9
dence: a µ = Lµ,µ ] 5C F 6β a µ, a 4 µ = Lµ,µ ] 9C F 5β a 4 µ A.4 with Lµ,µ = α s µ /α s µ, β = n f /3. twist-3 DA s : Φ 3π α i = 36α α α 3 ω ] 3π 7α 3 3 A.5 with the nonperturative parameters f 3π and ω 3π defined via matrix elements of the following local operators: ] π p ūσ µν γ 5 G αβ d = if 3π p α p µ g βν p β p µ g αν p α p ν g βµ p β p ν g αµ,a.6 π p ūσ µλ γ 5 D β,g αλ ]d 3 7 βūσ µλ γ 5 G αλ d = 3 4 f 3πω 3π p α p β p µ. A.7 The scale dependence of the twist-3 parameters is given y: µ π µ = Lµ,µ ] 4 β µ π µ, f 3π µ = Lµ,µ ] β 7CF f 3π ω 3π µ = Lµ,µ ] β 7CF 6 3 3 f 3π µ,a.8 f 3π ω 3π µ. A.9 The corresponding expressions for the twist-3 quark-antiquark DA are: twist-4 DA s: φ p 3π u = 3 f 3π µ π f π C / u ū 3 f 3πω 3π φ σ 3πu = 6u u 5 f 3π µ π f π Φ 4π α i = δ π ε πα α α α α 3, C / 4 u ū, µ π f π ω 3π C 3/ u ū. A. Ψ 4π α i = 3δ π µα α α 3 3 ε π α 3 ], Φ 4π α i = δ πα α α 3 3 ε π 3α 3 ], Ψ 4π α i = 3δ πα 3 α 3 3 ε π α 3 ], A. are the four three-particle DA s, where the nonperturative parameters δ π and ǫ π are defined as π p ū G αµ γ α d = iδ π f πp µ, A. and up to twist 5 corrections: π p ūd µ, G νξ ]γ ξ d 4 9 µū G νξ γ ξ d = 8 f πδ π ǫ πp µ p ν, A.3 with the scale-dependence: δ π µ = Lµ,µ ] 8C F 3β δ π µ, δ π ǫ πµ = Lµ,µ ] β δ π ǫ πµ. A.4
Note that the twist-4 parameter ω 4π introduced in 9] is replaced y ǫ π = /8ω 4π. Correspondingly, the two-particle DA s of twist 4 are: φ 4π u = { 3 δ π u ū 8δπ ǫ π uū 3uū u 3 5u 6u ln u } ū 3 5ū 6ū ln ū, A.5 ψ 4π u = 3 δ πc / u. A.6 These DA s are related to the original definitions 34] as φ 4π u = 6 g u B. Formulae for gluon radiative corrections u g vdv, ψ 4π u = dg u du. A.7 Here we collect the expressions for the hard-scattering amplitudes entering the factorization formulae 3. and the resulting imaginary parts of these amplitudes determining the radiative correction 4.7 to LCSR 4.4 for f Bπ, as well as the analogous expressions for LCSR 4.8 and 4. for the other two form factors. To compactify the formulae, we use the dimensionless variales r = q m, r = in the imaginary parts r = s/m and the integration variale : ρ = r ur r p q m, B. du = and introduce the cominations of logarithmic functions where Li x = x r Gx = Li x ln x ln x dt t r dρ r r, ln t is the Spence function, and x m x L x = ln x µ, L x = ln x B. ln m µ, B.3 m µ x. B.4 The imaginary parts of the hard-scattering amplitudes are taken at fixed q < m r <, analytically continuing these amplitudes in the variale s = p q or r. The result contains cominations of θ ρ, θρ and δρ and its derivatives. To isolate the spurious infrared divergences which one encounters y taking the imaginary part, we follow 4] and introduce the usual plus-prescription r r dρ { θ ρ θρ } gρ ρ φρ = r r dρ { θ ρ θρ } gρ φρ φ, B.5 ρ
for generic functions φρ, gρ. Furthermore, to make the formulae for imaginary parts more explicit, we partially integrate the derivatives of δρ using, e.g.: r dρδ ρ φρ = r dρδρ d dρ δr φρ, B.6 omitting the terms with δr in all cases where φ =. B. Amplitudes for f Bπ LCSR T = Gr ρ r r r u ρ r r r u r r r ρ r r r u r r r ρ ρ r ρ ρ 3ln m µ ρ r r r u Gρ r r u ln r Gr ln ρ r r u 3ρ, B.7 ρ ] π Im r st = θ ρ r r r ρ L L r r ρ r ρ r r θρ r r r ρ L r ρ r r r ρr ρ L L r L ρ ρ ρ ρ ] ρ r δρ ln r r ln r ln r r Li r 3Li r π r 4 3ln m µ 4 3ln m d dρ µ ], B.8
r r T p = ρ 4r r Gr r r u r ρ r r r r u ρ r r r r u 4r Gρ r r u r ρ r ln r r r u r r r r r r r r r r ρ r r u 33 r 6 r m ρ ρ ln µ r ρ r r r r r r r r u r r 3r r r ρ r ρ ln r r r r u ln ρ r 3 ρ r Gr ρ 4 r ρ, B.9 r = θ ρ L r r L r r r π Im st p r r ρ r r r ρ r r 4r θρ L r r r L ρ r ρ ρ r ρ r r L r r L ρ r ρ r r 3r r r ρ δρ ln r r r ln L r r r lnr r ρ r ρ r r ρ r 3 r r ln r ln r lnr r r π 6 4r m r 3 3ln µ 4 Li r r Li r m ] r 4 3ln d µ dρ δr, B. ] ] 3
3T σ = ρ r r ρ 4r r r u r Gr r r r u r ρ r r ρ r r r r u r r r r u r ρ r r 3r r r r ρ r r r r u 4r r r u r r r r u r Gρ r r r u r ρ r r r r r r r u ln r r ρ ln r 4 ρ r r r r r r r r u 4 r r u 3 r ρ 3r 8r r r 3r r r r r r r r r r ρ r r r u r r r u 5r r r r r r r r r u 3r 4r r r r r u r ρ 3 r r 3r r r ρ r r r r r r r ln ρ r r ρ 6 r ρ 3 5r 3 ρ r 3r 4r r r r r ρ r r r r r u 4r r r r u ln m µ Gr 4r r 3r r r r r r r r ρ r r r r u r r r r u 8 r ρ 3 r ρ r ρ r r r, B. r ρ 4
3 π Im st σ d = θ ρ r L r ρ dρ r L r r ρ r r r ρ r L r r r ρ r r ρ r ] ρ r L ρ d θρ r r L r ρ ρ dρ 3 r L ρ r ρ r ρ r L r r r ρ r 8r r r r r r ρ r L r r ρ r 4r r ρ r r r ρ r r ρ r r r r ρ L ρ r ρ 3 r r ρ 3r r r r r r ρ r r r r r r ρ r r 3r r ρ r r r r r r r ] 4r r r ρ r r r ρ π δρ 3 4r d dρ 3 r 4 r r ln r r ln r Li r r Li r r r r lnr r ln r ln r r r r r ln r lnr r r 5 3r r ln r m ] ln d r µ dρ ln r ln r r r r Li r 3 r Li r r r r r r r 3 r r r r ln r r r lnr ln r 3 r r r r r ln r lnr 4 3 r r r r r ln r r r r m lnr ln r µ m r 4 3ln d µ dρ δr d ]. B. dρ 5
As already mentioned, the aove formulae are otained in the MS scheme. To switch to the one-loop pole mass of quark the following expressions m ρ 3ln 4 µ T = ρ, B.3 m r ρρ 3ln 4 T p = µ r r ρ, B.4 m 3 ρρ r ρ 3 3 3ln T σ = µ 4 6ρ 3. B.5 have to e added to the hard-scattering amplitudes T, T p and T σ. respectively. The corresponding additions to the imaginary parts are m π Im s T = δρ 3ln µ 4 d, B.6 dρ r r π Im s T p m = δρ 3ln 3 r d µ 4 r dρ δr, B.7 3 m π Im st σ = δρ 3ln µ 4 r d dρ d r dρ δr d. B.8 dρ B. Amplitudes for f Bπ f Bπ LCSR T = r r r r r r ln r r ρ r r r ln r r r r u r ln r r ur r ρ r ρln ρ r r uuρ r r, B.9 r ρ π Im s T r = θ ρ r r ρ θρ r ρ δρ ] r r ρ r r r r ρ r r r r ln r r r r r uρ T p = 4Gr r r u r Gr r r r r u r r r ρ ρ m ρ ln µ ln r r ] r ], B. 4 r r u ρ r r u r r ρ r r r r u Gρ ln r r r r r u r r r r u ρ r r r r r r ρ 4 ρ ρ ln ρ ρ 3r r ρ 8 ρ, B. 6
π Im s T p = θ ρ r ] L r ρ r ρ L ρ r r ρ ρ r ρρ r L ρ r r 3r r r r ρ r r ρ r ] ρ δρ lnr 3 lnr r r r r r r ln r m 4 3ln µ L r 4Li r r 3 ln m ] d dρ δr, θρ r L r ρ r ρ µ r r ρ r r ln r 4 3 π B. T σ = r Gr 3u ρ r r r r 3r ρ 3u r r r r 3ρ 3r r u Gr 3 r u r r 3ρ 3 u r r Gρ r 3r r r u 3r u r r 3r ρ r3 r r r r 6r ln r ρ r r r r r r 3r ρ 3r r r u ln r r r r r r r 3r r r u r 3r r u r r 3r r ρ 3r r u 3r r r u 3r r u r r ρ r r 3r 6r ρ r r 3ρ 3 3r 3r u r 3u r r r 3ρ r r r 3 3r 5r r ln ρ 3r ρ r r r ρ 3 3r u 3 r u r r r 3r r 3r 3r r ρ r r 3ρ m ln µ 3r r u 3r u r r r r r r r r r r 3r r ρ r r r r r 3r ρ 7r 7r r r r 6r ρ r r 3ρ 8r r 3 ρ 3, B.3 7
3 r r π Im s T σ L r = θ ρ r ρ L r L ρ L r θρ ρ ρ ] d dρ d dρ L r r ρ L ρ r ρ r ρ r ρ r 3 r r r r ρ r L ρ ρ r r r ρ r r r r r r r ρ r r 3 r r ρ r 3r r ρ r r ρ r r r ρ r ] ρ 3 π δρ d 3 dρ lnr lnr r m L r ln lnr Li r µ ln r r L r ln r r lnr r r r r ln r lnr r r r d 3 r dρ r r 3 ln r lnr r r Li r Li r m 4 3ln d µ dρ δr d dρ r r m ln µ ], B.4 π Im s T p 3 r r π Im s T σ T =, T p = ρ T r r ρ 3 σ = m = δρ 3ln B.5 m 3ln µ 4 ρ, B.6 µ m = δρ 3ln µ m 3ln 4 µ 6ρ 3, B.7 d 4 4 dρ δr, B.8. B.9 d dρ d dρ δr d dρ 8
B.3 Amplitudes for f T Bπ LCSR T T r = r r ρ r r r u r r r r u Gr ρ r r r r u r ln r r ρ r r r r u r ln r r ρ r r r r u ρ r ρ 3 ρ π Im st T = θ ρ ln m µ Gr r r ρ r r r u r r r u ρ ln ρ Gρ r r r u r r r r r r ρ ρ ρ 4 ρ, B.3 r r r ρ r L r L r r ] ρ r ρ r r θρ r r ρ r L r r r ρ r ρr ρ L ρ L ρ L r r ρ ρ r ρ r r r r r r r ρ δρ ln r lnr ln r r r r lnr lnr m r ln Li r 3Li r µ ] ρ π m ] 4 3ln d µ, B.3 dρ T T p 3 = r r u Gr r r ρ r r ρ 3 r r Gr u r r r r ρ r r ln r u r r r ρ r r u ln ρ, r r ρ r r u u ln r 3Gρ r r u u B.3 9