Revisiting radiative leptonic B decay

In this paper, we summarize the existing methods of solving the evolution equation of the leading-twist B-meson LCDA. Then, in the Mellin space, we derive a factorization formula with next-to-leading-logarithmic (NLL) resummation for the form factors in the decay at leading power in . Furthermore, we investigate the power suppressed local contributions, factorizable non-local contributions (which are suppressed by and ), and soft contributions to the form factors. In the numerical analysis, which employs the two-loop-level hard function and the jet function, we find that both the resummation effect and the power corrections can sizably decrease the form factors. Finally, the integrated branching ratios are also calculated for comparison with future experimental data.


Introduction
The radiative leptonic decay of the meson is of interest as it is the most important channel to extract the parameters of the -meson light-cone distribution amplitudes (LCDAs) and to test the factorization theorem when the emitted photon is energetic. A precision study of this mode can also be helpful in decreasing the background to the purely leptonic decay process , which is important for determining the CKM matrix element . The radiative leptonic decay amplitude is defined by the QCD matrix element In the rest frame of the meson with momentum , it is convenient to introduce two light-cone vectors and with the definitions At leading order in QED and considering the constraints from the Ward identity, the amplitude can be parameterized as [1,2] A(B → γ ℓ ν) → where the contribution from final-state radiation is accounted for by the redefinition of the axial form factor .
At leading power in , the QCD factorization formula has been derived for the form factors [3,4] and was confirmed under the framework of soft-collin- ear effective theory (SCET) [5,6]. The form factors can be factorized into a convolution of the hard function, jet function, and -meson LCDA. The hard function arises from the matching between heavy-to-light current in the QCD and SCET I operators, and it has been calculated up to the two-loop level [1]. The jet function can be obtained from the matching between SCET I and SCET II [5], and the next-to-next-to-leading-order (NNLO) correction has been recently obtained [7]. The matrix elements of the SCET II operators are actually the definition of -meson LCDA. All the ingredients in the factorization formula depend on the factorization scale, and the radiative corrections lead to large logarithmic terms, which need to be resummed. For the hard function, the three-loop anomalous dimension is known [8][9][10][11], and the two-loop-level anomalous dimensions both for themeson LCDA and the jet function have recently been calculated [7,12]. Therefore, the sufficient condition for a complete NLL resummation is readily available. It was first derived in [13] by performing a Laplace transformation of the -meson LCDA.
Although the factorization formula of decay is well established at leading power, the power corrections are important for finite bottom-quark mass. The power-suppressed corrections of were considered at tree level [1], where a symmetry-preserving form factor was introduced to parameterize the non-local power corrections. The soft contribution from the endpoint region of the momentum of the light quark inside the meson was first studied using dispersion relation and quark-hadron duality in [14]. The QCD correction to the soft contribution at one loop and the contribution from three-particle LCDAs were computed in [2]. In a comprehensive study on the local and non-local power suppressed contributions, the soft contribution and the higher-twist contribution to the decay were presented [15]. The contribution from the hadronic structure of the photon, which can be defined by the matrix elements of power-suppressed SCET operators, was studied in [16,17]. Moreover, based on transverse-momentum-dependent factorization, the power corrections to decay were investigated in [18]. All these studies indicate that the power-suppressed contribution is sizable and should not be neglected in the determination of the first inverse moment of the -meson LCDA.
In this paper, we make improvements from two aspects. The first is to derive the scale-independent factorization formula at the NLL level in the Mellin space, and the second is to perform the phenomenological analysis after combining the NLL leading-power result with the power corrections. This paper is organized as follows. In the next section, we review the evolution of the leadingtwist -meson LCDA. In the third section, we derive the scale-independent factorization formula of the form factors and discuss the power-suppressed contributions; this is followed by the phenomenological analysis. Concluding discussions are presented in the final section. The -meson LCDA is one of the most important ingredients of the QCD factorization formula for exclusive decays. The two-particle LCDAs of the meson in the heavy-quark effective theory (HQET) can be obtained from the coordinate-space matrix elements [19] The LCDAs (in curly brackets) can be expanded around . In the limit, , , themeson LCDAs in the momentum space are defined through the Fourier transformation e iωτΦ± (τ) .
At leading power, only is relevant in the decay, and the evolution equation of is the wellknown Lange-Neubert equation [20]: µ where is the renormalization scale. At the one-loop level, the anomalous dimensions are with the "plus" function defined as ] .

(8) B
In the position space, the evolution equation of themeson LCDA takes the form [21] d d ln µΦ where at the one-loop level, Chinese Physics C Vol. 44, No. 12 (2020) 123106 Whether in the momentum space or in the position space, the evolution equation of the -meson LCDA is the integro-differential equation. It is difficult to solve directly: it must be simplified by an integral transformation. To date, there exist the following treatments: • Performing the Fourier transformation with respect to (or the Mellin transform for ): Then, the evolution kernel of is obtained as [22] φ Through the inverse Fourier transformation, we arrive at the solution to the evolution equation in the momentum space: where , and the functions and take the form • Performing the Mellin transformation to the evolution equation in the position space [21,23]: In the Mellin space, the evolution equation takes a simple form: witĥ where at the one-loop level. The solution in the Mellin space can be obtained directly:φ • It was found that, if the -meson LCDA is transformed into the so-called "dual" space, the evolution kernel is diagonalized [24]. The LCDA in the dual space can be obtained by which satisfies an ordinary differential equation: It is then simple to write the solution: The method mentioned above is equivalent, and the LC-DAs are different expressions of an identical objective. Because the momentum space and the position space are related through a standard Fourier transformation, we are able to derivẽ is related to by definition: Then, we havẽ At the one-loop level, the most convenient method is to work in the dual space since the Bessel function is the eigenfunction of the Lange-Neubert kernel, which is confirmed in [25,26]. The Lange-Neubert kernel can be expressed as a logarithm of the generator of special conformal transformations along the light cone. When the eigenfunction of the generator is transformed to the momentum space, it is simply the Bessel function in Eq. .

B
The two-loop-level anomalous dimension of themeson LCDA was first calculated in the coordinate space in [12]; it is more simply expressed in the dual space: Chinese Physics C Vol. 44, No. 12 (2020) 123106 where and . This equation is also transformed into the momentum space in [7], resulting in the two-loop-level Lange-Neubert equation. The advantage of solving the evolution equation at the twoloop level in the dual space does not hold as the two-loop evolution kernel is not diagonal in this space. On the contrary, the elegant form of the evolution equation (Eq. (16)) in the Mellin space is maintained. Thus, Eq. (18) is still the solution to the evolution equation up to the twoloop level with [23]

B
In a recent paper [13], an alternative approach to solving the evolution equation at the two-loop level was proposed. The essential idea of this approach is to perform a Laplace transformation on the -meson LCDA, where is a fixed reference scale, which can be used to eliminate the logarithmic moment in the factorization formula of . Then, one could derive with the definition where starts from the two-loop level, and the specific expression can be seen in [13]. After the Laplace transformation, the solution to the evolution equation readsφ The normalization depends on the factorization scale through where the quantities , and are given as .
We note that the LCDA is related to through At leading power in , the QCD factorization formula for the form factors can be written as At the one-loop level, the hard function and jet function are given as [5,6,27] The results of the two-loop level hard function and jet function can be found in [1,7]. As the hard function and jet function contain large logarithmic terms, it is important to perform resummation to improve the convergence of the perturbative series. The first complete NLL resummation is given in [13]: where the jet function with . We now derive the scale independent factorization formula at the NLL level in the Mellin space. The evolution of the -meson LCDA is known in the Mellin space. We thus need to perform the Mellin transformation to the jet function, although it is not well-defined.Alternatively, we follow the method in [13] to replace the first argument of by a derivative operator, i.e., Taking advantage of Eq. (22), we have (39)

B
Employing the evolution function of the hard function and -meson LCDA, we obtain where and are the evolution factor of the hard function and -meson decay constant in the HQET, respectively; the specific expression can be found in the appendices of [1,2]. The parameter depends on the factorization scale through µ d dµ The resummed factorization formula in the Mellin space seems more compact than Eq. (37). Now, we turn to the power suppressed contributions. At leading power, and are equal due to the lefthandedness of the weak interaction current and the helicity-conservation of the quark-gluon interaction in the high-energy limit, although this relation might be broken by the power corrections. In [15], the power suppressed contribution is separated into the symmetry-preserving part and symmetry-breaking part , i.e., When the power suppressed contribution is from the region where , in which denotes the separation between the quark-photon vertex and the weak current, it is called the soft contribution. Soft contributions with higher twist -meson LCDA are considered in [15]. Because they are highly suppressed and numerically small, we neglect them in this study. Then, the symmetry breaking part only contains the local contribution and can be written as The symmetry preserving part can be divided into three parts, i.e., , and the explicit expressions for the first two parts are Chinese Physics C Vol. 44, No. 12 (2020) 123106 where and are the three particle twist-3 and twist-4 -meson LCDAs, respectively. is the "genuine" twist-three contribution to the LCDA [15]. The soft contribution with QCD corrections is obtained as where is the factor relating the QCD decay constant of the -meson to the HQET decay constant, and and are the Borel mass and threshold parameter, respectively. The evolution kernel . The effective LCDA takes the form To obtain this function, one must generalize the photon momentum from to , calculate the generalized hard-collinear function in this Euclidean region, perform the dispersion treatment to the convolution of the hard-collinear function with the -meson LCDA, and finally take the limit . The soft contribution actually includes the hadronic effect of the photon; it must overlap with the contribution of the photon LCDA, which will be investigated in future work.

Phenomenological analysis
The fundamental nonperturbative inputs entering the factorization formula of decay include the twoparticle and three-particle -meson distribution amplitudes up to the twist-four accuracy. The decay constant B B of the -meson and the parameters appear in the dispersion approach. In the numerical analysis, we employ the following three-parameter model for the leading twistmeson LCDA [15] ϕ where is the confluent hypergeometric function of the second kind. In dual space, this model has a simpler expression: where is the confluent hypergeometric function of the first kind. In the leading power factorization formula, only the first inverse moment and the logarithmic moments enter the factorization formula; they are defined by For the three-parameter model, the first inverse moment and the first two logarithmic moments are obtained as If the parameter , the three-parameter model is simply the familiar exponential model [19] ϕ which is set as our default model. To estimate the error from the models, we let vary in the region . We then employ two models with and . For the default model, , whose determination has been discussed extensively in the context of exclusive -meson decays (see [28][29][30][31] for further discussion). Here, we employ , which is consistent with the calculations of the semileptonic form factors with -meson LCDAs in the framework of light-cone sum rules [28]. The leading twist -meson LCDA with the three-parameter model is plotted in Fig. 1. In the factorization formula (31), a new parameter is introduced to eliminate : it can be determined once the parameters , and are given. In addition, by utilizing the three-parameter -meson LCDA model, the logarithmic moments are also determined. All the parameters are listed in Table 1.
The higher-twist LCDAs must incorporate the correct low-momentum behaviour and satisfy the equations of motion. All the suggested models can be obtained as particular cases of a more general ansatz: Here, the function obeys the following normalization condition: The following results can then be derived: Taking advantage of the above results, the NLP contribution with and corrections can be obtained as To highlight the influence of the power suppressed contributions, we present the numerical results of the form factors with different contributions: where the photon energy is fixed at GeV, and and are set to be free parameters. The leading power result in the first parentheses is borrowed from [13], and the power suppressed contributions include the symmetry breaking local term, the symmetry preserving term, and the soft contribution. It is obvious that power corrections are sizeable and more important than the terms. Therefore, the power suppressed contribution must play an important role in the determination of . We leave a more detailed study of the subleading power corrections for a future work.
To test the effect of the large logarithm resummation, we plot the dependence of the leading power form factor in Fig. 2. For the leading logarithmic resummation, we employ the tree-level hard function and jet function as well as the one-loop level anomalous dimension and two-loop cusp anomalous dimension. For the NLL resummation case, we follow the convention in [13]. In [13], the contributions from the NNLO hard function and jet function are also considered. Strictly speaking, to resum the logarithmic terms in the NNLO jet function, we need a three-loop anomalous dimension of themeson LCDA, which has not yet been obtained. While it is not phenomenologically as important as the hard-collinear scale, is actually close to the soft scale . From  Fig. 2, we can see that the NLL resummation effect significantly decreases the LL result, and that the NNLO result is approximately 5% smaller than the NLL result.    The form factors including both the LP contribution and NLP corrections are plotted in Fig. 3, where for LP contribution, we adopt the same result as in [13]. Compared with the LP result in Fig. 2, the power corrections significantly decrease the form factors, and the symmetry breaking effect from the NLP local contribution is also sizable. The uncertainties are denoted by the band in Fig.  3. To obtain the uncertainties, we consider various sources, including the following: decay constant ; the first inverse moment ; the hard scale and hard-collinear scale (the same as in [13]); the models of the leading twist -meson LCDA in Table 1 To guarantee the reliability of our calculation, we cut the photon energy at . We integrate over the differential decay width in the interval and then multiply it by the lifetime of the meson to obtain the branching ratio . If we fix , the branching ratio reads where the uncertainties come from the same source as that for the form factors (except for ). The dependence of the branching ratio on is presented in Fig. 4, where the parameter varies in the interval . We can see that the large uncertainty prevents us from precisely determining the parameter . It is thus important to reduce the uncertainty in the parameter estimates, especially the uncertainty in the meson LCDA, and to obtain a more precise prediction of the power suppressed contribution.

Summary
The radiative leptonic decay mode is interesting both theoretically and experimentally. It plays an irreplaceable role in the determination of parameters of the -meson LCDA. The factorization-scale dependence of the -meson LCDA is governed by the LN evolution equation, which is an integro-differential equation and is not easily solved. We summarized the existing method of   solving the LN evolution equation, for both the one-loop and two-loop anomalous dimensions. We then derived a factorization formula with NLL resummation for the form factors appearing at leading power in the Mellin space, which is equivalent to the one obtained in [13] but written in a more compact form. The power corrections to the are sizeable, and much effort has been put into investigating the NLP contributions. In this paper, we included the power suppressed local contributions, factorizable non-local contributions (which are suppressed by and ), and soft contributions. In the numerical analysis, we found that the NLL-re- summation effect significantly decreases the leadingpower form factors, and that the NNLO correction brings approximately 5% additional reduction. The NLP contributions are combined with the leading-power NNLO contributions and also manifestly decrease the form factors. We also calculated the integrated branching fractions of the decay. The large uncertainty from various sources makes it difficult to determine the parameter and other logarithmic moments accurately. In future work, we will consider the NLP corrections more systematically in the framework of the SCET and in the hope of reducing the theoretical uncertainty.